On the modelling of the effective longitudinal diffusion in bi-continuous chromatographic beds

We report on the possibility to extend to bi-continuous packings the two models for the effective longitudinal diffusion Deff, or B-term band broadening, recently proposed for discontinuous chromatographic beds. In bi-continuous packings, like monolithic columns, solutes experience a connected end-to-end pathway in both the mobile and stationary zones, as opposed to discontinuous packings, wherein the stationary adsorptive zone is distributed over a set of isolated elements. Since it is unclear whether a densely packed bed of spherical particles should be treated as a continuous or a bi-continuous medium, this extension is also crucial to fully understand the behaviour of packed particle beds. The proposed models for the effective longitudinal diffusion Deff originate from the adoption of the Two Zone Moment Analysis (TZMA) method by which Deff can be expressed as a linear combination of two essential quantities γm and γs, referred to as effective zone-diffusion factors. In the present work we propose two analytical models for γm and γs that now cover both the discontinuous and the bi-continuous case. To validate the theory, several bi-continuous packings are investigated, including the tetrahedral skeleton model (TSM), six different Triple Periodic Minimal Surface (TPMS) monoliths and randomly packed beds of spheres. For all of these, the models provide highly accurate results for Deff over a wide range of porosities and zone retention factors k″. The comparison with literature experimental data for both monolithic silica columns and columns packed with fully porous and porous-shell particles is also presented.

On the contribution of the top and bottom walls in micro-pillar array columns and related high-aspect ratio chromatography systems

We report on a steady-state based, and hence highly accurate numerical modelling study of the effect of the top and bottom wall in the current generation of micro-pillar array columns. These have a mesoporous retention layer that not only covers the pillar walls but also the bottom wall. Our results show that the performance of these columns can in general not be improved by also covering the top wall with the same layer, despite the increased column symmetry this approach would offer. The reason for this is that the local species retardation caused by a retentive layer is much stronger than the pure flow arresting effect of an uncovered wall. At least, this has a crucial impact in high aspect-ratio systems such as micro-pillar array columns because these require a small inter-pillar distance to promote mass transfer together with a large channel depth to enable a sufficiently high flow rate. On the other hand, a notable improvement could be made if micro-pillar array would be produced without having a retentive layer at the bottom. At Péclet number Pe = 50 and aspect ratio AR = 5 for flow-channels, this gain amounts up to about 4.5 h-units at a zone retention factor k'' = 2 and 1.75 h-units at k'' = 16 (gain scales almost linearly with Pe). To verify these results, we also considered another high aspect-ratio system with a simplified geometry: the open-tubular channel with a flat-rectangular cross-section. This led to very similar observations, thus confirming the findings for the micro-pillar array. The results produced in the present study also allow us to conclude that the classic modelling paradigm adopted in chromatography, which is based on the independency and hence additivity of the hCm- and hCs-contributions, can lead to large modelling errors in chromatographic systems with a high aspect-ratio, even when their geometry is so simple as that of a straight open-tubular channel with constant cross-section. Indeed, when both zones are treated independently, the analysis misses how the vertical diffusion through the retentive layer helps suppressing the vertical gradients in the mobile zone. The diffusion through this layer occurs in a ratio of k''Ds/Dm (Dm being the diffusion coefficient in mobile phase zone and Ds being the diffusion coefficient in stationary phase zone), such that at high retention factors this diffusion contribution even becomes the dominant one.

An alternative general model for the effective longitudinal diffusion in chromatographic beds filled with ordered porous particles

The two-zone moment-analysis method for the determination of the dispersion tensor in hierarchical retentive porous media has been adopted to compute and model the effective longitudinal diffusion Deff, or equivalently the B-term band broadening, in chromatographic beds filled with ordered porous particles. On the one hand, this approach offers accurate numerical results for Deff while keeping computational expenses low. On the other hand, it also gives direct insight for the analytical modelling, readily revealings the two main essential quantities (resp. referred to as the mobile-zone and stationary-zone effective diffusion factors γm and γs) that contribute to Deff. Modelling these two main parameters provided us with two new analytical models for Deff: a general one, valid for diluted and concentrated packings and accurate in the whole range of relevant intra-particle diffusion coefficient Dpz, and an approximate one, reliable for diluted packings and accurate also for concentrated packings with low to intermediate values of Dpz. The large advantage of both models is that they do not need any fitting parameter because all the required information is incorporated into the experimentally accessible geometric obstruction factor in the mobile phase originating from the tortuosity of the through-pore space (limiting case of fully solid particles without any retention). These models hence serve as an alternative to the Effective Medium Theory (EMT) models used so far in the literature. To validate the theory, five ordered geometries have been investigated. The accuracy of the general model proposed has been quantified and found to be comparable with that of the 3rd order approximate Torquato model for four geometries, even for macro-porosities close to the close-packing limit. The case of a 2-d triangular array of ellipsoidal particles with different elongations is also investigated to show the general validity and applicability of the models.

Moment analysis for predicting effective transport properties in hierarchical retentive porous media

We report on a new homogenization approach to solve, with drastically improved speed and accuracy, the general advection-diffusion equation in hierarchical porous media with localized diffusion and adsorption/desorption processes, thus opening the way to a much deeper understanding of the band broadening process in chromatographic systems. The proposed robust and efficient moment-based approach allows us to compute the exact local and integral concentration moments and hence provides exact solutions for the effective velocity and dispersion coefficients of migrating solute particles. Innovative to the proposed method is also that it not only produces the exact effective transport parameters of the long-time asymptotic solution, but also their entire transient. The analysis of the transient behaviour can be used, for example, to properly identify the time and length scales needed to achieve the macro-transport conditions. If the hierarchical porous media can be represented as the periodic repetition of a unit lattice cell, the method only requires the solution of the time-dependent advection-diffusion equations for the zeroth order and first-order exact local moments, exclusively on the unit cell. This implies an enormous reduction of the computational efforts and a significant improvement of the accuracy of the results when compared to the direct numerical simulation (DNS) approaches which require flow domains that are long enough to achieve steady-state conditions, and hence often cover tens to hundreds of unit cells. The reliability of the proposed method is verified by comparing its predictions with DNS results, in one, two and three dimensions, in both transient and asymptotic conditions. The influence of top and bottom no-slip walls on the separation performance of chromatographic columns with micromachined porous and nonporous pillars is discussed in detail.

A comparison of 2 micron inner diameter open tubular column liquid chromatography with pressure-driven isocratic, slip-flow, and electrochromatographic modes of operation: a theoretical study

Modifications to the flow profile used in open tube capillary liquid chromatography (OT-CLC) include using slip-flow walls and using electroosmosis as a fluid pump as practiced in electrochromatography. These modifications are implemented experimentally by changing the capillary surface and solvent conditions which results in the change of boundary conditions at the capillary wall. In this paper we employ a theory-based study and compare the zone broadening of simple solutes using parabolic flow from a liquid pump, slip-flow from a highly hydrophobic inner surface with water eluent, and electroosmosis for the conditions of pure water and dilute salt utilizing 2 µm inner diameter OT capillaries. In general, the two types of flow other than parabolic exhibit thin zones in the early part of the chromatogram, consistent with previous studies of slip-flow and electroosmotic flow used in electrochromatography. Electrochromatography is shown to yield higher efficiency and less zone broadening than parabolic and slip-flow conditions used in this study. Nonetheless, it is found that the zone standard deviations are shown to be similar for these flow profiles as is the number of plates for these different flow profiles under the conditions utilized here. It is revealed that these modifications do not warrant the effort to maintain the special solvent conditions when compared to gradient elution OT-CLC, which gives a nearly constant peak width throughout the chromatogram, is easiest to implement, and is the method of choice for complex analysis.