An In-Depth Investigation of Micropillar Array Columns for the Separation of Finite-Sized Particles by Hydrodynamic Chromatography.

The exact moment approach (EMA) is adopted to predict, without any fitting parameters, the plate height curves for polystyrene microparticles of different sizes in micropillar array columns performed by hydrodynamic chromatography. The EMA allows us to decouple the contribution of horizontal and vertical dispersion terms and thus investigate the influence of pillar height and interpillar distance on separation performance. In the convection-controlled regime, we found that axial dispersion is mainly controlled by the vertical dispersion term, the latter being due to the flow-arresting top and bottom walls. This vertical contribution can be estimated from the axial dispersion in rectangular, open tubular channels formed between the pillars. Henceforth, plate height curves can be accurately predicted by simply adding the estimated vertical term to the horizontal dispersion term evaluated from 2D simulations. This finding allowed us to understand that, to improve separation performance, it is advisible to decrease the interpillar distance (expected result) and decrease the pillar height (counterintuitive result).

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.

Theoretical computation of the band broadening in micro-pillar array columns.

We have investigated the possibility to establish a theoretical plate height expression for the band broadening in the most widely used micro-pillar array column format, i.e., a cylindrical pillar array wherein the pillar walls and the channel bottom are coated with a thin layer of meso­porous material. Assuming isotropic diffusion in the shell-layer, it was found that the vertical diffusive transport along the porous shell-layer covering the pillar walls significantly suppresses the band broadening originating from the vertical migration velocity gradients. As the vertical transport in the shell-layer increases linearly with the retention equilibrium constant K, this leads to an anomalous dependency on the retention factor. Indeed, instead of increasing with k'' and following the classic (1+ak''+bk''2)/(1 + k'')2-dependency governing a classic Taylor-Aris system, the variation of the mobile zone mass transfer resistance term hCm in a 3D pillar array with bottom-wall retention goes through a maximum (resp. factor 1.5 (k''=4) and 2 (k''=16) difference between observed and classic Taylor-Aris behaviour). This effect increases with increasing pillar heights and increasing reduced velocities. Because of this complex k''-dependency, it proves very cumbersome to establish a general plate height equation covering all conditions. Instead, a plate height expression was established that is limited up to k''=4, but remains accurate for higher k''-values for cases where the ratio of pillar height over inter-pillar distance remains below 5. It can however be anticipated the proposed analytical model is only valid in a rather limited range around the presently considered external porosity of ε=0.5.

Nonadditivity and Nonlinearity of Mobile and Stationary Zone Mass Transfer Resistances in Chromatography.

Using a two-zone moment analysis (TZMA) method based on Brenner's generalized dispersion theory for two-dimensional (2D) and three-dimensional (3D) periodic media, we investigated the mechanisms for dispersion in particulate media for liquid chromatography. This was done using a set of plate height data covering an unprecedented wide range of retention factors, diffusion coefficients, and velocities, all computed with unequaled accuracy. Applying Giddings' additivity test, based on alternatingly making the diffusion coefficient in the mobile and stationary zones infinitely large, the dispersion data clearly indicate a lack of additivity. Although this lack could be directly understood by identifying the existence of multiple parallel mass transfer paths, the additivity assumption interestingly overestimates the true C term band broadening (typically by more than 10%, depending on conditions and dimensionality of the system). However, Giddings originally asserted the occurrence of parallel paths would always lead to an underestimation of the dispersion. The origin of the lack of additivity is analyzed in detail and qualitatively explained. Finally, we also established a generic framework for the modeling of the effect of the reduced velocity and the retention coefficient on the C term in ordered chromatographic media. This led to the introduction of a new expression for the mobile zone mass transfer term, which, unlike the currently used literature expression, contains the complete k″ dependency.

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.

Impact of transversal vortices on the performance of open-tubular liquid chromatography.

Open-Tubular Liquid Chromatography (OTLC) is currently limited by two shortcomings, namely the low ratio of adsorbing area to the channel volume and the large values of the Height Equivalent of the Theoretical Plate (HETP) due to Taylor-Aris dispersion. Previous work focusing on axial dispersion of nonadsorbing solutes showed how it is possible to tame the Taylor-Aris effect by inducing transversal velocity components acting alongside the main pressure-driven axial flow. We here analyze the impact of transversal flow on the separation resolution in OTLC, where simultaneous equilibrium adsorption at the channel walls is superimposed to the analyte transport in the mobile phase. A three-dimensional steady flow generated by the combination of a pressure-driven flow and an electroosmotically-induced transversal flow is used as case study. Flows geometries possessing regular and chaotic streamlines are created by axially-invariant and periodically-alternate arrangements of the electrodes along the channel walls, respectively. By enforcing Brenner's macrotransport approach, we predict the column length achieving a prescribed level of resolution as a function of the Péclet number and of the species affinity towards the stationary adsorbing phase. Results show that the presence of transversal flows can lower sensitively the dependence of the column length on the Péclet number. Flows possessing chaotic streamlines prove the most efficient choice at large eluent velocities and low values of the column adsorption constant.

Taming Taylor-Aris dispersion through chaotic advection.

Taylor-Aris dispersion represents an undesired phenomenon in pressure-driven liquid chromatography, often responsible for the unchecked increase of the Height Equivalent of the Theoretical Plate (HETP) when high throughput operating conditions are sought. Previous work on the subject showed how it is possible to contain the augmented dispersion in empty microchannels by inducing cross-sectional velocity components yielding an overall helical structure of the flow streamlines. Here, we explore the possibility of further reducing axial dispersion by devising flow conditions that give rise to chaotic streamlines. A three-dimensional steady flow generated by the combination of a pressure-driven Poiseuille flow and an electroosmotically-induced spatially periodic flow is used as a case study. Brenner's macrotransport approach is used to predict the axial dispersion coefficient of a diffusing solute in flows possessing regular, partially chaotic and globally chaotic kinematic features. Accurate Lagrangian-stochastic simulations of particle ensembles are used to validate the predictions obtained through Brenner's paradigm. Our findings suggest that the Taylor-Aris phenomenon can be altogether suppressed in the limit of globally chaotic kinematics. A theoretical interpretation of this outcome is developed by combining macrotransport theory with established results of the spectral approach to mixing in advecting-diffusing chaotic flows.

Invariant manifold approach for quantifying the dynamics of highly inertial particles in steady and time-periodic incompressible flows.

The dynamics of finite-sized particles with large inertia are investigated in steady and time-dependent flows through the numerical solution of the invariance equation, describing the spatiotemporal evolution of the slow/inertial manifold representing the effective particle velocity field. This approach allows for an accurate reconstruction of the effective particle divergence field, controlling clustering/dispersion features of particles with large inertia for which a perturbative approach is either inaccurate or not even convergent. The effect of inertia on heavy and light particles is quantified in terms of the rate of contraction/expansion of volume elements along a particle trajectory and of the maximum Lyapunov exponent for systems exhibiting chaotic orbits, such as the time-periodic sine-flow on the 2D torus and the time-dependent 2D cavity flow.

Swelling and Drug Release in Polymers through the Theory of Poisson-Kac Stochastic Processes.

Experiments on swelling and solute transport in polymeric systems clearly indicate that the classical parabolic models fail to predict typical non-Fickian features of sorption kinetics. The formulation of moving-boundary transport models for solvent penetration and drug release in swelling polymeric systems is addressed hereby employing the theory of Poisson-Kac stochastic processes possessing finite propagation velocity. The hyperbolic continuous equations deriving from Poisson-Kac processes are extended to include the description of the temporal evolution of both the Glass-Gel and the Gel-Solvent interfaces. The influence of polymer relaxation time on sorption curves and drug release kinetics is addressed in detail.

A Non-Isothermal Moving-Boundary Model for Continuous and Intermittent Drying of Pears.

A non-isothermal moving-boundary model for food dehydration, accounting for shrinkage and thermal effects, is proposed and applied to the analysis of intermittent dehydration in which air temperature, relative humidity, and velocity vary cyclically in time. The convection-diffusion heat transport equation, accounting for heat transfer, water evaporation, and shrinkage at the sample surface, is coupled to the convection-diffusion water transport equation. Volume shrinkage is not superimposed but predicted by the model through the introduction of a point-wise shrinkage velocity. Experimental dehydration curves, in continuous and intermittent conditions, are accurately predicted by the model with an effective water diffusivity Deff(T) that depends exclusively on the local temperature. The non-isothermal model is successfully applied to the large set of experimental data of continuous and intermittent drying of Rocha pears.