Modeling of dispersion in a polymeric chromatographic monolith

Multiple-open-tubular column enabling transverse diffusion. Part 2: Channel size distribution and structure optimization

Multiple-open-tubular columns enabling transverse diffusion (MOTTD) are made of straight, parallel, and cylindrical flow channels separated by a mesoporous stationary phase. In Part 1, a model of band broadening along MOTTD columns accounting for longitudinal diffusion, the trans-channel velocity bias, and mass transfer resistance in the stationary phase was proposed and validated. In this Part 2, the model is completed by considering the impact of short-range inter-channel velocity biases on the MOTTD plate number. These velocity biases are caused by the wide distribution of the channel diameters. Different ratios, ρ, of the average inner diameter, 2<r_c>, of the flow channels to their closest center-to-center distance d (d= 5 μm, ρ= 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9) with a relative standard deviation (RSD) increasing from 0 to 50% are considered. The zone retention factor k₁ was increased from 1 to 25. The complete model of band broadening is validated after adjustment to dispersion data obtained by 1) the lattice-Boltzmann method for modeling fluid flow, 2) a random-walk particle-tracking (RWPT) technique to address advective-diffusive transport, and 3) by considering two distinct populations of flow channels (inner radii r_c,1=<r_c>(1-RSD) and r_c,2=<r_c>(1+RSD)) arranged at the nodes of a hexagonal compact array. The completed model of band broadening in MOTTD columns reveals that the RSD of the channel diameters has only a moderate impact on the optimum plate number of MOTTD columns: the relative increase of the minimum plate height do not exceed 30% even for the largest RSDs. However, when the mass transfer of the analyte is governed by its slow rate of transverse diffusion across the MOTTD column, the plate height can be increased by up to 100% at high average velocities. Regarding the best trade-off between analysis speed and column performance at a fixed pressure drop of 400 bar, irrespective of the zone retention factor and RSD of the distribution of the channel diameters, the fastest analyses are recommended for MOTTD columns having a small structural parameter ρ. In contrast, for the longest analysis times, the largest values of ρ are required to maximize the performance of MOTTD columns.

Transport properties and size exclusion effects in wide-pore superficially porous particles

The effects of hydrodynamic radius on the transport of solute molecules in packed beds of wide-pore superficially porous particles (SPP) are studied using pore-scale simulation. The free molecular diffusion rate varies with radius through the Stokes-Einstein relation. Lattice Boltzmann and Langevin methods are used to model fluid motion and the transport of an ensemble of solute molecules in the fluid, providing statistics on solute concentration, flux, molecule age and residence time, as a function of depth in the SPP. Intraparticle effective diffusion and bed dispersion coefficients are calculated and correlated with the hydrodynamic radius and accessible porosity. The relative importance of convection and diffusion are found to depend on the molecule (tracer) size through the diffusion rate, and convection effects are more significant for larger, slower-diffusing molecules. When larger molecules are utilized, the intraparticle concentration is reduced in proportion to the local particle porosity, leading to a natural definition of the accessible porosity used in size exclusion chromatography (SEC). Although the pore shape is complex, the SEC constant K can be calculated directly from simulation. Simulation demonstrates that the effective diffusion coefficient is elevated near the particle hull, which is largely open to interstitial flow, and decreases with depth into the particle. All molecules studied here have transport access to the entire particle depth, although the accessible volume at a given depth depends on their size. The first passage time into the particle is well predicted by the diffusion rate, but residence time is influenced by convection, shortening the average visit duration. These results are of interest in "perfusion" chromatography where convection is thought to increase separation efficiency for large biomolecules.

Impact of diffusion on transverse dispersion in two-dimensional ordered and random porous media

Solute dispersion in fluid flow results from the interaction between advection and diffusion. The relative contributions of these two mechanisms to mass transport are characterized by the reduced velocity ν, also referred to as the Péclet number. In the absence of diffusion (i.e., when the solute diffusion coefficient D_{m}=0 and ν→∞), divergence-free laminar flow of an incompressible fluid results in a zero-transverse dispersion coefficient (D_{T}=0), both in ordered and random two-dimensional porous media. We demonstrate by numerical simulations that a more realistic realization of the condition ν→∞ using D_{m}≠0 and letting the fluid flow velocity approach infinity leads to completely different results for ordered and random two-dimensional porous media. With increasing reduced velocity, D_{T} approaches an asymptotic value in ordered two-dimensional porous media but grows linearly in disordered (random) structures depending on the geometrical disorder of a structure: a higher degree of heterogeneity results in a stronger growth of D_{T} with ν. The obtained results reveal that disorder in the geometrical structure of a two-dimensional porous medium leads to a growth of D_{T} with ν even in a uniform pore-scale advection field; however, lateral diffusion is a prerequisite for this growth. By contrast, in ordered two-dimensional porous media the presence of lateral diffusion leads to a plateau for the transverse dispersion coefficient with increasing ν.

Morphological Properties of Methacrylate-Based Polymer Monoliths: From Gel Porosity to Macroscopic Inhomogeneities

Shaping chemical interfaces of hard and soft matter materials into physical morphologies that guarantee excellent transport properties is of central importance for technologies relying on adsorption, separation, and reaction at the interface. Polymer monoliths with a hierarchically structured pore space, for example, are widely used in flow-driven processes, whose efficiency depends on the morphology of the support material over several length scales. Compared with alternative support structures, particularly silica monoliths, polymer monoliths yield lower efficiency, which suggests a suboptimal morphology. Based on physical reconstruction by serial block-face scanning electron microscopy we evaluate the structural features of a methacrylate-based polymer monolith from the pore scale to the column scale. The morphological data reveal a homogeneous polymer skeleton with a solute-impenetrable core-porous shell architecture and a heterogeneous macropore space that suffers from inhomogeneities at the short-range and the transcolumn scale. Although the morphology of the polymer phase is favorable to efficient mass transport, the performance of the polymer monolith is limited by severe transcolumn gradients in macroporosity and macropore size. We propose to overcome these morphological limitations by pursuing a preparation strategy that involves active rather than passive shaping of the macropore space, for example, by using silica monoliths as templating structures for polymer monolith preparation.

Topological analysis of non-granular, disordered porous media: determination of pore connectivity, pore coordination, and geometric tortuosity in physically reconstructed silica monoliths

Different approaches are applied and compared, which are universally applicable to quantify pore coordination, pore and pore-throat connectivity, and geometric tortuosity.

Predictive model of solute transport with reversible adsorption in spatially periodic hierarchical porous media

Solute transport in hierarchical porous media with reversible adsorption is predicted by the volume averaging method. A transient macroscopic advection-diffusion equation is derived to describe the multiscale solute transport problem. The theoretical expression of the dispersion tensor is obtained which is the function of pore-scale velocity profile. Steady closure equations are derived to calculate the dispersion tensor. With the aid of pore-scale simulation in unit cells of the hierarchical porous media, the dispersion tensor can be calculated. The model is verified by comparing its predictions and obtaining favorable agreement with results of direct numerical simulations and with experimental data for columns comprised of ordered, porous pillars. It is straightforward to use the model to predict the solute transport behavior in fixed beds packed with particles.

Analytical silica monoliths with submicron macropores: current limitations to a direct morphology-column efficiency scaling

Shrinking the structural elements of a particulate bed or monolith (i.e., the particle or domain size) yields more efficient columns only when the homogeneity of the bed can be conserved in that process. We investigate this complex issue for a set of 2nd generation analytical silica monoliths with macropores reaching submicron dimensions using chromatographic methods, mercury intrusion porosimetry, scanning electron microscopy, and confocal laser scanning microscopy (CLSM), and present eddy dispersion simulations and a chord length distribution analysis for the CLSM-based physical reconstructions at macropore resolution. The combined results allow us to identify relevant morphological advances made from 1st to 2nd generation monoliths and additionally highlight the current limitations to a direct morphology-efficiency scaling with respect to the performance that can be accomplished in HPLC practice with these columns. Whereas the improvement in radial homogeneity from 1st to 2nd generation silica monoliths is represented by a dramatic increase in column efficiency, the further reduction of macropore size in the 2nd generation monoliths does not lead to the expected improvement of plate height data, although these monoliths realize submicron macropores at a simultaneously conserved bulk macropore space homogeneity and negligible radial heterogeneity. Our study implies that limitations to further improved column efficiency arise from the intrinsic border effects of the used 4.6mm i.d. analytical columns. This includes the sample distribution onto the monoliths and asynchronous sample collection through the endfittings at the column inlet and outlet, respectively. Only when these effects are reduced will additionally improved 2nd generation monoliths live up to column efficiencies, which are envisioned for them based on their morphological properties.

Geometrical and topological measures for hydrodynamic dispersion in confined sphere packings at low column-to-particle diameter ratios

At low column-to-particle diameter (or aspect) ratio (d(c)/d(p)) the kinetic column performance is dominated by the transcolumn disorder that arises from the morphological gradient between the more homogeneous, looser packed wall region and the random, dense core. For a systematic analysis of this morphology-dispersion relation we computer-generated a set of confined sphere packings varying three parameters: aspect ratio (d(c)/d(p)=10-30), bed porosity (ɛ=0.40-0.46), and packing homogeneity. Plate height curves were received from simulation of hydrodynamic dispersion in the packings over a wide range of reduced velocities (v=0.5-500). Geometrical measures derived from radial porosity and velocity profiles were insufficient as morphological descriptors of the plate height data. After Voronoi tessellation of the packings, topological information was obtained from the statistical moments of the free Voronoi volume (V(free)) distributions. The radial profile of the standard deviation of the V(free) distributions in the form of an integral measure was identified as a quantitative scalar measure for the transcolumn disorder. The first morphology-dispersion correlation for confined sphere packings deepens our understanding of how the packing microstructure determines the kinetic column performance.