Influence of moderate Joule heating on electroosmotic flow velocity, retention, and efficiency in capillary electrochromatography

The influence of Joule heating on electroosmotic flow velocity, the retention factor of neutral analytes, and separation efficiency in capillary electrochromatography was investigated theoretically and experimentally. A plot of electrical current against the applied electrical field strength was used to evaluate the Joule heating effect. When the mobile phase concentration of Tris buffer exceeded 5.0 mM in the studied capillary electrochromatography systems using particulate and monolithic columns (with an accompanying power level of heat dissipation higher than 0.35 W/m), the Joule heating effect became clearly noticeable. Theoretical models for describing the variation of electroosmotic flow velocity with increasing applied field strength and the change of retention factors for neutral analytes with electrical field strength at higher Tris buffer concentrations were analyzed to explain consequences of Joule heating in capillary electrochromatography. Qualitative agreement between experimental data and implications of the theoretical model analysis was observed. The decrease of separation efficiency in capillary electrochromatography with macroporous octadecylsilica particles at high buffer concentration can be also attributed to Joule heating mainly via the increased axial diffusion of the analyte molecules and dispersion of solute bands by a nonuniform electroosmotic flow profile over the column cross-section. However, within a moderate temperature range, the contribution of the macroscopic velocity profile in the column arising from radial temperature gradients is insignificant.     

Focused electrophoretic motion and selected electrokinetic dispensing of particles and cells in cross-microchannels

The electrokinetic focusing and the resultant accelerated electrophoretic motion of polystyrene particles and red blood cells were visualized in microfluidic cross-channels. The experimentally measured width of the focused stream and the measured velocity increase of particles and cells at different voltage ratios follow the proposed analytical formula within the experimental error. The attained velocity increase is insensitive to the particle size, particle property (i.e., particle or cell), and particle trajectory. By solving the electrical potential field in the cross-channel at the experimental conditions, we demonstrate that the squeezed electrical field lines in the channel intersection determine the shape of the focused stream, and the nonuniform distribution of axial electrical field strength underlies the variation of particle/cell electrophoretic velocity through the focusing region. However, the dielectrophoretic force resulting from the nonuniform electrical field in the intersection seems to push the acceleration region of particles and cells slightly in the downstream direction. We have also achieved the single particle/cell dispensing by instantly triggering an electrical pulse perpendicular to the focused particulate flow in a double-cross microchannel. The electrokinetic manipulation of particle/cell in microchannels demonstrated in this work can be used for developing integrated lab-on-a-chip devices for studies of cells.     

Assessment of Joule heating and its effects on electroosmotic flow and electrophoretic transport of solutes in microfluidic channels

Joule heating is inevitable when an electric field is applied across a conducting medium. It would impose limitations on the performance of electrokinetic microfluidic devices. This article presents a 3-D mathematical model for Joule heating and its effects on the EOF and electrophoretic transport of solutes in microfluidic channels. The governing equations were numerically solved using the finite-volume method. Experiments were carried out to investigate the Joule heating associated phenomena and to verify the numerical models. A rhodamine B-based thermometry technique was employed to measure the solution temperature distributions in microfluidic channels. The microparticle image velocimetry technique was used to measure the velocity profiles of EOF under the influence of Joule heating. The numerical solutions were compared with experimental results, and reasonable agreement was found. It is found that with the presence of Joule heating, the EOF velocity deviates from its normal "plug-like" profile. The numerical simulations show that Joule heating not only accelerates the sample transport but also distorts the shape of the sample band.     

Streaming potential and electroosmotic flow in heterogeneous circular microchannels with nonuniform zeta potentials: requirements of flow rate and current continuities

Real surfaces are typically heterogeneous, and microchannels with heterogeneous surfaces are commonly found due to fabrication defects, material impurities, and chemical adsorption from solution. Such surface heterogeneity causes a nonuniform surface potential along the microchannel. Other than surface heterogeneity, one could also pattern the various surface potentials along the microchannels. To understand how such variations affect electrokinetic flow, we proposed a model to describe its behavior in circular microchannels with nonuniform surface potentials. Unlike other models, we considered the continuities of flow rate and electric current simultaneously. These requirements cause a nonuniform electric field distribution and pressure gradient along the channel for both pressure-driven flow (streaming potential) and electric-field-driven flow (electroosmosis). The induced nonuniform pressure and electric field influence the electrokinetic flow in terms of the velocity profile, the flow rate, and the streaming potential.     

Joule heating effects on electrokinetic flows with conductivity gradients

Instability occurs in the electrokinetic flow of fluids with conductivity and/or permittivity gradients if the applied electric field is beyond a critical value. Understanding such an electrokinetic instability is significant for both improved transport (via the suppressed instability) and enhanced mixing (via the promoted instability) of liquid samples in microfluidic applications. This work presents the first study of Joule heating effects on electrokinetic microchannel flows with conductivity gradients using a combined experimental and numerical method. The experimentally observed flow patterns and measured critical electric fields under Joule heating effects to different extents are reasonably predicted by a depth-averaged numerical model. It is found that Joule heating increases the critical electric field for the onset of electrokinetic instability because the induced fluid temperature rise and in turn the fluid property change (primarily the decreased permittivity) lead to a smaller electric Rayleigh number.     

Analytical study of Joule heating effects on electrokinetic transportation in capillary electrophoresis

Electric fields are often used to transport fluids (by electroosmosis) and separate charged samples (by electrophoresis) in microfluidic devices. However, there exists inevitable Joule heating when electric currents are passing through electrolyte solutions. Joule heating not only increases the fluid temperature, but also produces temperature gradients in cross-stream and axial directions. These temperature effects make fluid properties non-uniform, and hence alter the applied electric potential field and the flow field. The mass species transport is also influenced. In this paper we develop an analytical model to study Joule heating effects on the transport of heat, electricity, momentum and mass species in capillary-based electrophoresis. Close-form formulae are derived for the temperature, applied electrical potential, velocity, and pressure fields at steady state, and the transient concentration field as well. Also available are the compact formulae for the electric current and the volume flow rate through the capillary. It is shown that, due to the thermal end effect, sharp temperature drops appear close to capillary ends, where sharp rises of electric field are required to meet the current continuity. In order to satisfy the mass continuity, pressure gradients have to be induced along the capillary. The resultant curved fluid velocity profile and the increase of molecular diffusion both contribute to the dispersion of samples. However, Joule heating effects enhance the sample transport velocity, reducing the analysis time in capillary electrophoretic separations.     

Effect of Joule heating on electrokinetic transport

The Joule heating (JH) is a ubiquitous phenomenon in electrokinetic flow due to the presence of electrical potential gradient and electrical current. JH may become pronounced for applications with high electrical potential gradients or with high ionic concentration buffer solutions. In this review, an in-depth look at the effect of JH on electrokinetic processes is provided. Theoretical modeling of EOF and electrophoresis (EP) with the presence of JH is presented and the important findings from the previous studies are examined. A numerical study of a fused-silica capillary PCR reactor powered by JH is also presented to extend the discussion of favorable usage of JH.     

Numerical analysis of the thermal effect on electroosmotic flow and electrokinetic mass transport in microchannels

Joule heating is present in electrokinetically driven flow and mass transport in microfluidic systems. Nowadays, there is a trend of replacing costly glass-based microfluidic systems by the disposable, cheap polymer-based microfluidic systems. Due to poor thermal conductivity of polymer materials, the thermal management of the polymer-based microfluidic systems may become a problem. In this study, numerical analysis is presented for transient temperature development due to Joule heating and its effect on the electroosmotic flow (EOF) and mass species transport in microchannels. The proposed model includes the coupling Poisson-Boltzmann (P-B) equation, the modified Navier-Stokes (N-S) equations, the conjugate energy equation, and the mass species transport equation. The results show that the time development for both the electroosmotic flow field and the Joule heating induced temperature field are less than 1 s. The Joule heating induced temperature field is strongly dependent on channel size, electrolyte concentration, and applied electric field strength. The simulations reveal that the presence of the Joule heating can result in significantly different characteristics of the electroosmotic flow and electrokinetic mass transport in microchannels.     

A lattice Boltzmann algorithm for electro-osmotic flows in microfluidic devices

In this paper, a finite-difference-based lattice Boltzmann (LB) algorithm is proposed to simulate electro-osmotic flows (EOF) with the effect of Joule heating. This new algorithm enables a nonuniform mesh to be adapted, which is desirable for handling the extremely thin electrical double layer in EOF. The LB algorithm has been validated by simulating a problem with an available analytical solution and it is found that the numerical results predicted by the algorithm are in good agreement with the analytical solution. The LB algorithm is also applied to modeling a mixed electro-osmotic/pressure driven flow in a channel. The numerical results show that Joule heating plays an important role in EOF.     

Rapid concentration of deoxyribonucleic acid via Joule heating induced temperature gradient focusing in poly-dimethylsiloxane microfluidic channel

This paper reports rapid microfluidic electrokinetic concentration of deoxyribonucleic acid (DNA) with the Joule heating induced temperature gradient focusing (TGF) by using our proposed combined AC and DC electric field technique. A peak of 480-fold concentration enhancement of DNA sample is achieved within 40 s in a simple poly-dimethylsiloxane (PDMS) microfluidic channel of a sudden expansion in cross-section. Compared to a sole DC field, the introduction of an AC field can reduce DC field induced back-pressure and produce sufficient Joule heating effects, resulting in higher concentration enhancement. Within such microfluidic channel structure, negative charged DNA analytes can be concentrated at a location where the DNA electrophoretic motion is balanced with the bulk flow driven by DC electroosmosis under an appropriate temperature gradient field. A numerical model accounting for a combined AC and DC field and back-pressure driven flow effects is developed to describe the complex Joule heating induced TGF processes. The experimental observation of DNA concentration phenomena can be explained by the numerical model.     

Theoretical and numerical analysis of temperature gradient focusing via Joule heating

We present a detailed theoretical and numerical analysis of temperature gradient focusing (TGF) via Joule heating-an analytical species concentration and separation technique relying upon the dependence of an analyte's velocity on temperature due to the temperature dependence of a buffer's ionic strength and viscosity. The governing transport equations are presented, analyzed, and implemented into a quasi-1D numerical model to predict the resulting temperature, velocity, and concentration profiles along a microchannel of varying width under an applied electric field. Numerical results show good agreement with experimental trials presented in previous work. The model is used to analyze the effects of varying certain geometrical and experimental parameters on the focusing performance of the device. Simulations also help depict the separation capability of the device, as well as the effectiveness of different buffer systems used in the technique. The analysis provides rule-of-thumb methodology for implementation of TGF into analytical systems, as well as a fundamental model applicable to any lab-on-a-chip system in which Joule heating and temperature-dependent electrokinetic transport are to be analyzed.     

Numerical modeling of Joule heating-induced temperature gradient focusing in microfluidic channels

Temperature gradient focusing (TGF) is a recently developed technique for spatially focusing and separating ionic analytes in microchannels. The temperature gradient required for TGF can be generated either by an imposed temperature gradient or by Joule heating resulting from an applied electric field that also drives the flow. In this study, a comprehensive numerical model describing the Joule heating induced temperature development and TGF is developed. The model consists of a set of governing equations including the Poisson-Boltzmann equation, the Laplace equation, the Navier-Stokes equations, the energy equations and the mass transport equation. As the thermophysical and electrical properties including the liquid dielectric constant, viscosity, and electric conductivity are temperature-dependent, these governing equations are coupled, and therefore the coupled governing equations are solved numerically by using a CFD-based numerical method. The numerical simulations agree well with the experimental results, suggesting the valid mathematical model presented in this study.     

On the surface conductance, flow rate, and current continuities of microfluidics with nonuniform surface potentials

The characteristics of electrokinetic flow in a microchannel depend on both the nature of surface potentials, that is, whether it is uniform or nonuniform, and the electrical potential distribution along the channel. In this paper, the nonlinear Poisson-Boltzmann equation is used to model the electrical double layer and the lattice Boltzmann model coupled with the constraint of current continuity is used to simulate the microfluidic flow field in a rectangular microchannel with a step variation of surface potentials. This current continuity, including surface conduction, convection, and bulk conduction currents, has often been neglected in the literature for electroosmotic flow with nonuniform (heterogeneous) microchannels. Results show that step variation of ion distribution caused by step variation surface potential will influence significantly the electrical potential distribution along the channel and volumetric flow rate. For the system considered, we showed that the volumetric flow rate could have been overestimated by as much as 70% without consideration of the current continuity constraint.     

Towards high concentration enhancement of microfluidic temperature gradient focusing of sample solutes using combined AC and DC field induced Joule heating

It is challenging to continuously concentrate sample solutes in microfluidic channels. We present an improved electrokinetic technique for enhancing microfluidic temperature gradient focusing (TGF) of sample solutes using combined AC and DC field induced Joule heating effects. The introduction of an AC electric field component services dual functions: one is to produce Joule heat for generating temperature gradient; the other is to suppress electroosmotic flow. Consequently the required DC voltages for achieving sample concentration by Joule heating induced TGF are reduced, thereby leading to smaller electroosmotic flow (EOF) and thus backpressure effects. As a demonstration, the proposed technique can lead to concentration enhancement of sample solutes of more than 2500-fold, which is much higher than the existing literature reported microfluidic concentration enhancement by utilizing the Joule heating induced TGF technique.     

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.     

Analytical and numerical study of Joule heating effects on electrokinetically pumped continuous flow PCR chips

Joule heating is an inevitable phenomenon for microfluidic chips involving electrokinetic pumping, and it becomes a more important issue when chips are made of polymeric materials because of their low thermal conductivities. Therefore, it is very important to develop methods for evaluating Joule heating effects in microfluidic chips in a relatively easy manner. To this end, two analytical models have been established and solved using the Green's function for evaluating Joule heating effects on the temperature distribution in a microfluidic-based PCR chip. The first simplified model focuses on the understanding of Joule heating effects by ignoring the influences of the boundary conditions. The second model aims to consider practical experimental conditions. The analytical solutions to the two models are particularly useful in providing guidance for microfluidic chip design and operation prior to expensive chip fabrication and characterization. To validate the analytical solutions, a 3-D numerical model has also been developed and the simultaneous solution to this model allows the temperature distribution in a microfluidic PCR chip to be obtained, which is used to compare with the analytical results. The developed numerical model has been applied for parametric studies of Joule heating effects on the temperature control of microfluidic chips.     

Focusing particles by induced charge electrokinetic flow in a microchannel

A novel method of sheathless particle focusing by induced charge electrokinetic flow in a microchannel is presented in this paper. By placing a pair of metal plates on the opposite walls of the channel and applying an electrical field, particle focusing is achieved due to the two pairs of vortex that constrain the flow of the particle solution. As an example, the trajectories of particles under different electrical fields with only one metal plate on one side channel wall were numerically simulated and experimentally validated. Other flow focusing effects, such as the focused width ratio (focused width/channel width) and length ratio (focused length/half‐length of metal plate) of the sample solution, were also numerically studied. The results show that the particle firstly passes through the gaps between the upstream vortices and the channel walls. Afterwards, the particle is focused to pass through the gap between the two downstream vortices that determine the focused particle position. Numerical simulations show that the focused particle stream becomes thin with the increases in the applied electrical field and the length of the metal plates. As regards to the focused length ratio of the focused stream, however, it slightly increases with the increase in the applied electrical field and almost keeps constant with the increase in the length of the metal plate. The size of the focused sample solution, therefore, can be easily adjusted by controlling the applied electrical field and the sizes of the metal plates.     

Joule heating induced transient temperature field and its effects on electroosmosis in a microcapillary packed with microspheres

The Joule heating induced transient temperature field and its effect on the electroosmotic flow in a capillary packed with microspheres is analyzed numerically using the control-volume-based finite difference method. The model incorporates the coupled momentum equation for the electroosmotic velocity, the energy equations for the Joule heating induced temperature distributions in both the packed column and the capillary wall, and the mass and electric current continuity equations. The temperature-dependent physical properties of the electrolyte solution are taken into consideration. The characteristics of the Joule heating induced transient development of temperature and electroosmotic flow fields are studied. Specifically, the simulation shows that the presence of Joule heating causes a noticeable axial temperature gradient in the thermal entrance region and elevates a significant temperature increment inside the microcapillary. The temperature changes in turn greatly affect the electroosmotic velocity by means of the temperature-dependent fluid viscosity, dielectric constant, and local electric field strength. Furthermore, the model predicts an induced pressure gradient to counterbalance the axial variation of the electroosmotic velocity so as to maintain the fluid mass continuity. In addition, under specific conditions, the present model is validated by comparing with the existing analytical model and experimental data from the literature.     

Effect of the Joule heating and of the material voids on free-convective transport in porous or fibrous media with applied electrical fields

The effect of the geometry of the soil in electrokinetic application has been studied by using capillary models of annular geometry. The Joule heating generation has been included as a primary effect of temperature development leading to buoyancy flows. The heat transfer model has been formulated for conduction-dominated regime. The results of this model have been coupled with the motion equation to obtain the analytical hydrodynamic velocity profile. Numerical illustrations, demonstrating the effect of the cross-sectional area of the annular region on the velocity field, have been included. It is observed that a substantial effect on the magnitude of such velocity field for different parameters of the system. The results are useful to obtain better understanding of the role of the soil geometry in potential soil cleaning field operations.     

Peak dispersion and contributions to plate height in nonaqueous capillary electrophoresis at high electric field strengths: propanol as background electrolyte solvent

Peak dispersion effects in nonaqueous capillary electrophoretic separations of aromatic anionic analytes were investigated in a propanolic background electrolyte solution. Poly(glycidylmethacrylate-co-N-vinylpyrrolidone) coating was applied to the capillary to suppress the electroosmotic flow and to improve the repeatability of the migration times. Electrical field strengths up to 2000 Vcm(-1) were applied in separations and the separation efficiencies were compared with theoretical values calculated on the basis of plate height theory. The contributions to the total plate height were calculated for injection plug length, diffusion, Joule heating, electromigration dispersion, analyte adsorption to the capillary wall, and detector slit aperture length. Analyte diffusion coefficients were measured by Taylor dispersion method, while distribution constants were measured chromatographically. Agreement between the calculated and empirical results was fairly good even though some approximations were required. In most cases the longitudinal diffusion contribution governed the total plate height, while the contribution of Joule heating was insignificant even at exceptionally high field strengths used. The relatively long detection slit aperture was found to influence the separation efficiency strongly, while the other dispersion sources that were investigated were of minor importance, except for adsorption in the case of one analyte. With all analytes, the dispersive effect of longitudinal diffusion was reduced as the field strength was increased, leading to enhanced migration velocities and faster separations.