Electroosmotic Flow Phenomena in Packed Capillaries: From the Interstitial Velocities to Intraparticle and Boundary Layer Mass Transfer

Inherent peak compression of charged analytes in electrochromatography

This work resolves peak compression of charged analytes in CEC with strong cation-exchange stationary phase particles. By combining electrochromatographic peak shape analysis with the results of numerical simulations and confocal laser scanning microscopy in the packed capillaries, we identify electrical field-induced concentration polarization as the key physical phenomenon responsible for the inherent existence of local electrical field gradients on the scale of an individual support particle. Consequently, positive and negative field gradients exist between and inside the particles along the whole packing. Their intensity depends on the particles cation-selectivity (governed by the particles volume charge density and the mobile phase ionic strength) and the applied field strength. The interplay of these local field gradients with the analytes retention (intraparticle adsorption) determines whether fronting, tailing, or spiked analyte peaks are observed, and it provides a mechanism by which strongly retained analytes can be eluted over long distances with little zone dispersion. Our analysis explains the "anomalous" peak compression effects with strong cation-exchange particles, which have been reported more than a decade ago (Smith, N. W., Evans, M. B., Chromatographia 1995, 41, 197-203) and since then remained largely unresolved.

Modeling wall effects in capillary electrochromatography

A volume averaging technique for modeling electroosmotic and pressure driven flow in microchannels is applied to a packed capillary electrochromatography column. This model can be applied to fluid flow in both porous and open channels and can account for porosity variation in the column due to packing and zeta potential mismatches between the wall and the packing material. Numerical results are presented and compared with experimental results from the literature. Several different porosity models are shown to produce similar concentration profiles and allow the model to reproduce wall effects observed in experimental columns.

Fluid dynamics in capillary and chip electrochromatography

This review is concerned with the phenomenological fluid dynamics in capillary and chip electrochromatography (EC) using high-surface-area random porous media as stationary phases. Specifically, the pore space morphology of packed beds and monoliths is analyzed with respect to the nonuniformity of local and macroscopic EOF, as well as the achievable separation efficiency. It is first pointed out that the pore-level velocity profile of EOF through packed beds and monoliths is generally nonuniform. This contrasts with the plug-like EOF profile in a single homogeneous channel and is caused by a nonuniform distribution of the local electrical field strength in porous media due to the continuously converging and diverging pores. Wall effects of geometrical and electrokinetic nature form another origin for EOF nonuniformities in packed beds which are caused by packing hard particles against a hard wall with different zeta potential. The influence of the resulting, systematic porosity fluctuations close to the confining wall over a distance of a few particle diameters becomes aggravated at low column-to-particle diameter ratio. Due to the hierarchical structure of the pore space in packed beds and silica-based monoliths which are characterized by discrete intraparticle (intraskeleton) mesoporous and interparticle (interskeleton) macroporous spatial domains, charge-selective transport prevails within the porous particles and the monolith skeleton under most general conditions. It forms the basis for electrical field-induced concentration polarization (CP). Simultaneously, a finite and -- depending on morphology -- often significant perfusive EOF is realized in these hierarchically structured materials. The data collected in this review show that the existence of CP and its relative intensity compared to perfusive EOF form fundamental ingredients which tune the fluid dynamics in EC employing monoliths and packed beds as stationary phases. This addresses the (electro)hydrodynamics, associated hydrodynamic dispersion, as well as the migration and retention of charged analytes.

Correlated displacement-T2 MRI by means of a Pulsed Field Gradient-Multi Spin Echo Method

A method for correlated displacement-T2 imaging is presented. A Pulsed Field Gradient-Multi Spin Echo (PFG-MSE) sequence is used to record T2 resolved propagators on a voxel-by-voxel basis, making it possible to perform single voxel correlated displacement-T2 analyses. In spatially heterogeneous media the method thus gives access to sub-voxel information about displacement and T2 relaxation. The sequence is demonstrated using a number of flow conducting model systems: a tube with flowing water of variable intrinsic T2's, mixing fluids of different T2's in an "X"-shaped connector, and an intact living plant. PFG-MSE can be applied to yield information about the relation between flow, pore size and exchange behavior, and can aid volume flow quantification by making it possible to correct for T2 relaxation during the displacement labeling period Delta in PFG displacement imaging methods. Correlated displacement-T2 imaging can be of special interest for a number of research subjects, such as the flow of liquids and mixtures of liquids or liquids and solids moving through microscopic conduits of different sizes (e.g., plants, porous media, bioreactors, biomats).

Pore-Scale Dispersion in Electrokinetic Flow through a Random Sphere Packing

The three-dimensional velocity field and corresponding hydrodynamic dispersion in electrokinetic flow through a random bulk packing of impermeable, nonconducting spheres are studied by quantitative numerical analysis. First, a fixed bed with interparticle porosity of 0.38 is generated using a parallel collective-rearrangement algorithm. Then, the interparticle velocity field is calculated using the lattice-Boltzmann (LB) method, and a random-walk particle-tracking method is finally employed to model advection-diffusion of an inert tracer in the LB velocity field. We demonstrate that the pore-scale velocity profile for electroosmotic flow (EOF) is nonuniform even under most ideal conditions, including a negligible thickness of the electrical double layer compared to the mean pore size, a uniform distribution of the electrokinetic potential at the solid-liquid interface, and the absence of applied pressure gradients. This EOF dynamics is caused by a nonuniform distribution of the local electrical field strength in the sphere packing and engenders significant hydrodynamic dispersion compared to pluglike EOF through a single straight channel. Both transient and asymptotic dispersion behaviors are analyzed for EOF in the context of packing microstructure and are compared to pressure-driven flow in dependence of the average velocity through the bed. A better hydrodynamic performance of EOF originates in a still much smaller amplitude of velocity fluctuations on a mesoscopic scale (covering several particle diameters), as well as on the microscopic scale of an individual pore.

Modeling electroosmotic and pressure-driven flows in porous microfluidic devices: Zeta potential and porosity changes near the channel walls

This work presents analytical solutions for both pressure-driven and electroosmotic flows in microchannels incorporating porous media. Solutions are based on a volume-averaged flow model using a scaling of the Navier-Stokes equations for fluid flow. The general model allows analysis of fluid flow in channels with porous regions bordering open regions and includes viscous forces, permitting consideration of porosity and zeta potential variations near channel walls. To obtain analytical solutions problems are constrained to the linearized Poisson-Boltzmann equation and a variation of Brinkman's equation [Appl. Sci. Res., Sect. A 1, 27 (1947); 1, 81 (1947)]. Cases include one continuous porous medium, two adjacent regions of different porosities, or one open channel adjacent to a porous region, and the porous material may have a different zeta potential than that of the channel walls. Solutions are described for two geometries, including flow between two parallel plates or in a cylinder. The model illustrates the relative importance of porosity and zeta potential in different regions of each channel.

Diffusion NMR spectroscopy in supramolecular and combinatorial chemistry: an old parameter--new insights

Intermolecular interactions in solution play an important role in molecular recognition, which lies at the heart of supramolecular and combinatorial chemistry. Diffusion NMR spectroscopy gives information over such interactions and has become the method of choice for simultaneously measuring diffusion coefficients of multicomponent systems. The diffusion coefficient reflects the effective size and shape of a molecular species. Applications of this technique include the estimation of association constants and mapping the intermolecular interactions in multicomponent systems as well as investigating aggregation, ion pairing, encapsulation, and the size and structure of labile systems. Diffusion NMR spectroscopy can also be used to virtually separate mixtures and screen for specific ligands of different receptors, and may assist in finding lead compounds.

Electrohydrodynamics in hierarchically structured monolithic and particulate fixed beds

We have investigated the basic dependence of electroosmotic flow (EOF) velocity and hydrodynamic dispersion in capillary electrochromatography (CEC) on the variation of applied field and mobile phase ionic strengths employing silica-based particulate and monolithic fixed beds. These porous media have a hierarchical structure characterized by discrete intraparticle (intraskeleton) mesoporous and interparticle (interskeleton) macroporous spatial domains. While the macroporous domains contain quasi-electroneutral electrolyte solution, the ion-permselectivity (charge-selectivity) of the mesoporous domains determines the co-ion exclusion and counter-ion enrichment at electrochemical equilibrium (without superimposed electrical field) which depends on mesopore-scale electrical double layer (EDL) overlap and surface charge density. This adjustable, locally charge-selective transport realized under most general conditions forms the basis for concentration polarization (CP) induced by electrical fields superimposed in CEC. CP characterizes the formation of convective diffusion boundary layers with reduced (depleted CP zone) and increased (enriched CP zone) electrolyte concentration, respectively, at the anodic and cathodic interfaces in fixed beds containing the cation-selective, silica-based particles (or monolith skeleton). CP originates in the electrical field-induced coupled mass and charge transport normal to the charge-selective interfaces and has consequences for the EOF dynamics, hydrodynamic dispersion, and analyte retention in CEC. A secondary EDL with mobile counter-ionic space charge can be induced in the depleted CP zone leading to induced-charge EOF in the macroporous domains. It is characterized by a nonlinear dependence of the average EOF velocities on applied field strength and strong local velocity components tangential to the surface which enhance lateral pore-scale dispersion, thereby decreasing (axial) zone spreading. Differences in the pore space morphology of random-close sphere packings and monoliths criticially affect the intensity of CP and induced-charge EOF in these materials. CP is identified as a key phenomenon in CEC which also influences effective migration and the retention of charged analytes because the local intensity of CP inherently depends on applied field and mobile phase ionic strengths.

Electrokinetic Effects on the Transport of Charged Analytes in Biporous Media with Discrete Ion-Permselective Regions

The influence of external electrical fields on local concentration distributions and the mass transport of ionic background (buffer) species, as well as eluting co- and counterionic tracer molecules, was investigated in a fixed bed of native glass beads by confocal laser scanning microscopy and numerical simulations. Due to the negative surface charge of the porous glass beads and significant electrical double layer overlap, the intraparticle mesopore space becomes ion-permselective. This cation selectivity and the externally superimposed electrical fields induce concentration polarization in the bulk electrolyte solution adjacent to the particles. At the anodic hemisphere of a bead, the actual interplay of convection, diffusion, and electromigration leads to the formation of a convective-diffusion boundary layer with reduced ion concentrations relative to the bulk solution. At the opposite, cathodic hemisphere where counterions leave a bead in the direction of the applied field, electrolyte concentrations increase generating an enriched concentration polarization zone. Complementary data from quantitative confocal laser scanning microscopy and numerical simulations provide insight into the spatial variations of chemical and electrical potential gradients in the hierarchically structured material, including molar flux densities of the background ionic species, and reveal the elution dynamics of co- and counterionic analytes. These results demonstrate that concentration polarization in the external fluid domain, as well as the magnitude and sign of electrophoretic with respect to electroosmotic mobility in the ion-permselective domain, are major local contributions to coupled mass and charge transport, reflecting analyte retention, migration, and dispersion on a macroscopic scale.

Numerical analysis of electroosmotic flow in dense regular and random arrays of impermeable, nonconducting spheres

We present a numerical scheme for analyzing steady-state isothermal electroosmotic flow (EOF) in three-dimensional random porous media, involving solution of the coupled Poisson, Nernst-Planck, and Navier-Stokes equations. While traditional finite-difference methods were used to resolve the Poisson-Nernst-Planck problem, the (electro)hydrodynamics has been addressed with high efficiency using the lattice-Boltzmann method. The developed model allows simulation of electrokinetic transport under most general conditions, including arbitrary value and distribution of electrokinetic potential at the solid-liquid interface, electrolyte composition, and pore space morphology. The approach provides quantitative information on a spatial distribution of simulated velocities. This feature was utilized to characterize EOF fields in regular and random, confined and bulk packings of hard (i.e., impermeable, nonconducting) spheres. Important aspects of pore space morphology (sphere size distribution), surface heterogeneity (mismatch in electrokinetic potentials at confining wall and sphere surface), and fluid phase properties (electrical double layer thickness) were investigated with respect to their influence on the EOF dynamics over microscopic and macroscopic spatial domains. Most important is the observation of a generally nonuniform pore-level EOF velocity profile in the sphere packings (even in the thin double layer limit) which is caused by pore space morphology and which is in contrast to the pluglike velocity distribution in a single, straight capillary under the same conditions.

Perspective on concentration polarization effects in electrochromatographic separations

This work illustrates the appearance and electrohydrodynamic consequences of concentration polarization in the particulate and monolithic fixed beds used in capillary electrochromatography and related electrical-field assisted processes. Key property of most porous materials is the co-existence of bulk, quasi-electroneutral macroporous regions and mesoporous compartments which are ion-permselective (due to electrical double-layer overlap) causing different transport numbers for co-ionic and counterionic species, e.g., background electrolyte components, or the analytes. For a cathodic electroosmotic flow the (cation) permselectivity, together with diffusive and electrokinetic transport induces depleted and enriched concentration polarization zones at the anodic and cathodic interfaces, respectively, in dependence of the mobile phase ionic strength and applied electrical fields. At high field strength a secondary, nonequilibrium electrical double layer may be created in the depleted concentration polarization zones of a material stimulating electroosmosis of the second kind. The potential of this induced-charge electroosmosis with respect to nonlinear flow velocities and electrokinetic instability mixing (basically destroying the concentration polarization zones) is analyzed in view of the pore space morphology in random-close packings of spherical-shaped, porous particles and hierarchically structured monoliths. Possible applications based on a fine-tuning of the illustrated effects emerge for microfluidic pumping and mixing, or the intensification of sample recovery in adsorption processes. With this perspective we want to focus the attention on concentration polarization in electrochromatographic systems by presenting and discussing original data acquired on relevant microscopic as well as macroscopic scales, and point towards the importance of related effects in colloid and membrane science.*