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.     

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.     

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.     

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.     

A model for Joule heating-induced dispersion in microchip electrophoresis

This paper presents an analytical and parameterized model for analyzing the effects of Joule heating on analyte dispersion in electrophoretic separation microchannels. We first obtain non-uniform temperature distributions in the channel resulting from Joule heating, and then determine variations in electrophoretic velocity, based on the fact that the analyte's electrophoretic mobility depends on the buffer viscosity and hence temperature. The convection-diffusion equation is then formulated and solved in terms of spatial moments of the analyte concentration. The resulting model is validated by both numerical simulations and experimental data, and holds for all mass transfer regimes, including unsteady dispersion processes that commonly occur in microchip electrophoresis. This model, which is given in terms of analytical expressions and fully parameterized with channel dimensions and material properties, applies to dispersion of analyte bands of general initial shape in straight and constant-radius-turn channels. As such, the model can be used to represent analyte dispersion in microchannels of more general shape, such as serpentine- or spiral-shaped channels.     

Numerical modeling of the Joule heating effect on electrokinetic flow focusing

In electrokinetically driven microfluidic systems, the driving voltage applied during operation tends to induce a Joule heating effect in the buffer solution. This heat source alters the solution's characteristics and changes both the electrical potential field and the velocity field during the transport process. This study performs a series of numerical simulations to investigate the Joule heating effect and analyzes its influence on the electrokinetic focusing performance. The results indicate that the Joule heating effect causes the diffusion coefficient of the sample to increase, the potential distribution to change, and the flow velocity field to adopt a nonuniform profile. These variations are particularly pronounced under tighter focusing conditions and at higher applied electrical intensities. In numerical investigations, it is found that the focused bandwidth broadens because thermal diffusion effect is enhanced by Joule heating. The variation in the potential distribution induces a nonuniform flow field and causes the focused bandwidth to tighten and broaden alternately as a result of the convex and concave velocity flow profiles, respectively. The present results confirm that the Joule heating effect exerts a considerable influence on the electrokinetic focusing ratio.     

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.     

Finite sample effect in temperature gradient focusing

Temperature gradient focusing (TGF) is a new and promising equilibrium gradient focusing method which can provide high concentration factors for improved detection limits in combination with high-resolution separation. In this technique, temperature-dependent buffer chemistry is employed to generate a gradient in the analyte electrophoretic velocity. By the application of a convective counter-flow, a zero-velocity point is created within a microchannel, at which location the ionic analytes accumulate or focus. In general, the analyte concentration is small when compared with buffer ion concentrations, such that the focusing mechanism works in the ideal, linearized regime. However, this presumption may at times be violated due to significant sample concentration growth or the use of a low-concentration buffer. Under these situations the sample concentration becomes non-negligible and can induce strong nonlinear interactions with buffer ions, which eventually lead to peak shifting and distortion, and the loss of detectability and resolution. In this work we combine theory, simulation, and experimental data to present a detailed study on nonlinear sample-buffer interactions in TGF. One of the key results is the derivation of a generalized Kohlrausch regulating function (KRF) that is valid for systems in which the electrophoretic mobilities are not constant but vary spatially. This generalized KRF greatly facilitates analysis, allowing reduction of the problem to a single equation describing sample concentration evolution, and is applicable to other problems with heterogeneous electrophoretic mobilities. Using this sample evolution equation we have derived an understanding of the nonlinear peak deformation phenomenon observed experimentally in TGF. We have used numerical simulations to validate our theory and to quantitatively predict TGF. Our simulation results demonstrate excellent agreement with experimental data, and also indicate that the proper inclusion of Taylor dispersion is important for the accurate modeling of TGF. This work is an important first step towards the understanding and prediction of the more complex, nonlinear, and multi-species interactions which often occur in on-chip electrophoretic assays such as TGF.     

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.     

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.     

Electrochemical response and separation in cyclic electric field-flow fractionation

Electric field-flow fractionation (EFFF) is a separation technique that couples a lateral electric field with axial Poiseuille flow to separate particles on the basis of size and/or mobility. In unidirectional EFFF, the field rapidly decreases in time due to charging of the double layer. The field strength could be increased by performing EFFF with cyclic electric fields. In cyclic electric field-flow fractionation (CEFFF), a periodic voltage, which can be either sinusoidal or square-wave, is applied in the lateral direction. In this paper, we measure the electrochemical response of CEFFF, i.e., the current-time response for a given time-dependent voltage and then utilize this electrochemical response in a transport model to predict separation. The CEFFF device studied here comprises two gold-coated glass plates separated by a spacer. The transient current profiles are measured for a step change and cyclic square-shaped voltage. The current profile is compared with the equivalent circuit model, and is fitted to a sum of two decaying exponentials. The dependence of the electrochemical response on voltage, frequency, channel thickness, and salt concentration is studied. Next, the electrochemical data are utilized in the convection-diffusion equation to develop a model for separation by CEFFF. The equations are solved by using a combination of analytical and numerical techniques to determine the mean velocity and the dispersion coefficient of molecules, and to determine the effect of various parameters on the separation efficiency of the EFFF device. Also, the model predictions are compared with experimental data available in the literature.     

Estimation of Joule heating effect on temperature and pressure distribution in electrokinetic-driven microchannel flows

In this study we present simple analytical models that predict the temperature and pressure variations in electrokinetic-driven microchannel flow under the Joule heating effect. For temperature prediction, a simple model shows that the temperature is related to the Joule heating parameter, autothermal Joule heating parameter, external cooling parameter, Peclet number, and the channel length to channel hydraulic diameter ratio. The simple model overpredicted the thermally developed temperature compared with the full numerical simulation, but in good agreement with the experimental measurements. The factors that affect the external cooling parameters, such as the heat transfer coefficient, channel configuration, and channel material are also examined based on this simple model. Based on the mass conservation, a simple model is developed that predicts the pressure variations, including the temperature effect. An adverse pressure gradient is required to satisfy the mass conservation requirement. The temperature effect on the pressure gradient is via the temperature-dependent fluid viscosity and electroosmotic velocity.     

Taylor-Aris dispersion in temperature gradient focusing

Microfluidic temperature gradient focusing (TGF) uses an axial temperature gradient to induce a gradient in electrophoretic flux within a microchannel. When balanced with an opposing fluid flow, charged analytes simultaneously focus and separate according to their electrophoretic mobilities. We present a theoretical and experimental study of dispersion in TGF. We model the system using generalized dispersion analysis that yields a 1-D convection-diffusion equation that contains dispersion terms particular to TGF. We consider analytical solutions for the model under uniform temperature gradient conditions. Using a custom TGF experimental setup, we compare focusing measurements with the theoretical predictions. We find that the theory well represents the focusing behavior for flows within the Taylor-Aris dispersion regime.     

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.     

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.     

Quantitative temperature gradient focusing performed using background electrolytes at various pH values

Temperature gradient focusing (TGF) has previously been shown to be a practical technique for simultaneous concentration and separation of a variety of samples. In this paper, we demonstrate that TGF can be conducted at a wide range of pH values. Techniques for first-order prediction of the suitability of a given BGE for focusing are discussed. Buffer suitability for TGF is assessed experimentally by simultaneously concentrating and separating a pair of fluorescent analytes. One analyte is held at constant concentration for use as an internal standard while the concentration of the other dye is varied. Peak area is shown to vary linearly with the input dye concentration. A high degree of resolution (R(s) >3) of fluorescein and carboxyfluorescein, as well as for two LysoSensor-based dyes, is also observed. Foucusing and separation by TGF was successfully conducted quantitatively in BGE solutions of pH from 3.0 to 10.5.     

Joule heating and heat transfer in poly(dimethylsiloxane) microfluidic systems

Joule heating is a significant problem in electrokinetically driven microfluidic chips, particularly polymeric systems where low thermal conductivities amplify the difficulty in rejecting this internally generated heat. In this work, a combined experimental (using a microscale thermometry technique) and numerical (using a 3D "whole-chip" finite element model) approach is used to examine Joule heating and heat transfer at a microchannel intersection in poly(dimethylsiloxane)(PDMS), and hybrid PDMS/Glass microfluidic systems. In general the numerical predictions and the experimental results agree quite well (typically within +/- 3 degree C), both showing dramatic temperature gradients at the intersection. At high potential field strengths a nearly five fold increase in the maximum buffer temperature was observed in the PDMS/PDMS chips over the PDMS/Glass systems. The detailed numerical analysis revealed that the vast majority of steady state heat rejection is through lower substrate of the chip, which was significantly impeded in the former case by the lower thermal conductivity PDMS substrate. The observed higher buffer temperature also lead to a number of significant secondary effects including a near doubling of the volume flow rate. Simple guidelines are proposed for improving polymeric chip design and thereby extend the capabilities of these microfluidic systems.     

Development of a temperature gradient focusing method for in situ extraterrestrial biomarker analysis

Scanning temperature gradient focusing (TGF) is a recently described technique for the simultaneous concentration and separation of charged analytes. It allows for high analyte peak capacities and low LODs in microcolumn electrophoretic separations. In this paper, we present the application of scanning TGF for chiral separations of amino acids. Using a mixture of seven carboxyfluorescein succinimidyl ester-labeled amino acids (including five chiral amino acids) which constitute the Mars7 standard, we show that scanning TGF is a very simple and efficient method for chiral separations. The modulation of TGF separation parameters (temperature window, pressure scan rate, temperature range, and chiral selector concentration) allows optimization of peak efficiencies and analyte resolutions. The use of hydroxypropyl-beta-CD at low concentration (1-5 mmol/L) as a chiral selector, with an appropriate pressure scan rate ( -0.25 Pa/s) and with a low temperature range (3-25 degrees C over 1 cm) provided high resolution between enantiomers (Rs >1.5 for each pair of enantiomers) using a short, 4 cm long capillary. With these new results, the scanning TGF method appears to be a viable method for in situ trace biomarker analysis for future missions to Mars or other solar system bodies.     

Theoretical background for clinical and biomedical applications of electromigration techniques

This review presents the theoretical principles of analytical electrophoresis. The basic rules which control the movement of ionic species in electrostatic fields, together with the phenomenological theory of the resulting mass transport are analysed. The separation principles and capabilities of zone electrophoresis, isotachophoresis, isoelectric focusing, micellar electrokinetic chromatography and gel electrophoresis are evaluated. The most important effects accompanying electrophoresis, such as the production of Joule heat, electroosmosis and diffusion are discussed.     

Temperature gradient focusing of bio-analyte in a microfluidic channel dealing with non-Newtonian electrolyte considering temperature-dependent zeta potential

Temperature gradient focusing (TGF) relies on establishing a precise balance between the electrophoretic motility of a target analyte and the advective flow of the background electrolyte (BGE) to locally concentrate the analyte in a microfluidic configuration. This paper presents a finite-element-based numerical analysis where the coupled electric field and the transport equations are solved to describe the effects of the shear-dependent apparent viscosity of a non-Newtonian BGE on the localized concentration buildup of a charged bio-sample inside a microchannel by TGF via Joule heating. Effects of the temperature-dependent nature of the wall zeta potential and the flow behavior index (n) of BGE on the flow, thermal, and species concentration profiles inside the microchannel have been investigated. Study using a fluorescein-Na analyte sample shows that the maximum normalized analyte concentration (C_max/C₀) reduces as the zeta potential increases linearly with temperature. The maximum concentration enhancement is achieved when the BGE displays the Newtonian rheology. For example, C_max/C₀ increases 134- to 280-fold when n is increased from 0.8 to 1 (pseudoplastic regime) and again reduces to 190-fold when n increases further from 1 to 1.2 (dilatant regime).