Effects of surface heterogeneity on flow circulation in electroosmotic flow in microchannels

The characteristics of electroosmotic flow in a cylindrical microchannel with non-uniform zeta potential distribution are investigated in this paper. Two-dimensional full Navier–Stokes equation is used to model the flow field and the pressure field. The numerical results show the distorted electroosmotic velocity profiles and various kinds of flow circulation resulting from the axial variation of the zeta potential. The influences of heterogeneous patterns of zeta potential on the velocity profile, the induced pressure distribution and the volumetric flow rate are discussed in this paper. This work shows that using either heterogeneous patterns of zeta potential or a combination of a heterogeneous zeta potential distribution and an applied pressure difference over the channel can generate local flow circulations and hence provide effective means to improve the mixing between different solutions in microchannels.     

Expressions for evaluating the possibility of slip at the liquid-solid interface in open tube capillary electrochromatography

In this work, expressions are constructed and solved that describe the velocity field of electroosmotic flow (EOF) in open tube capillary electrochromatography (CEC) systems when the possibility of having unequal tangential velocities at the liquid-solid interface is considered and a slip condition is employed as a boundary condition for the velocity of the EOF at the capillary wall. The coupled equations of hydrodynamics (momentum balance equation) and electrostatics (Poisson equation) are solved numerically in order to obtain the distribution of the velocity field as well as the value of the volumetric flow rate in the open tube. Also, expressions for the velocity field and the volumetric flow rate of the EOF are presented that are valid for certain electrolytic systems and for certain parameter values for which analytical solutions to the momentum balance and Poisson equations could be obtained. The results presented in this work indicate that having slip in the velocity of the EOF at the wall of the capillary could (i) substantially increase the electroosmotic velocity in the plug-flow region of the radial domain of the open capillary tube and (ii) increase the portion of the radial domain of the open capillary tube where the velocity of the EOF has a plug-flow profile, which in turn could increase the average velocity and volumetric flow rate of the EOF in the open capillary tube. Furthermore, the modeling approach and the results presented in this work indicate a method for experimentally evaluating the possibility of having slip in the velocity of the EOF at the capillary wall.     

Helmholtz-Smoluchowski velocity for viscoelastic electroosmotic flows

Many biofluids such as blood and DNA solutions are viscoelastic and exhibit extraordinary flow behaviors, not existing in Newtonian fluids. Adopting appropriate constitutive equations these exotic flow behaviors can be modeled and predicted reasonably using various numerical methods. However, the governing equations for viscoelastic flows are not easily solvable, especially for electroosmotic flows where the streamwise velocity varies rapidly from zero at the wall to a nearly uniform velocity at the outside of the very thin electric double layer. In the present investigation, we have devised a simple method to find the volumetric flow rate of viscoelastic electroosmotic flows through microchannels. It is based on the concept of the Helmholtz-Smoluchowski velocity which is widely adopted in the electroosmotic flows of Newtonian fluids. It is shown that the Helmholtz-Smoluchowski velocity for viscoelastic fluids can be found by solving a simple cubic algebraic equation. The volumetric flow rate obtained using this Helmholtz-Smoluchowski velocity is found to be almost the same as that obtained by solving the governing partial differential equations for various viscoelastic fluids.     

Frequency-dependent laminar electroosmotic flow in a closed-end rectangular microchannel

This article presents an analysis of the frequency- and time-dependent electroosmotic flow in a closed-end rectangular microchannel. An exact solution to the modified Navier-Stokes equation governing the ac electroosmotic flow field is obtained by using the Green's function formulation in combination with a complex variable approach. An analytical expression for the induced backpressure gradient is derived. With the Debye-Hückel approximation, the electrical double-layer potential distribution in the channel is obtained by analytically solving the linearized two-dimensional Poisson-Boltzmann equation. Since the counterparts of the flow rate and the electrical current are shown to be linearly proportional to the applied electric field and the pressure gradient, Onsager's principle of reciprocity is demonstrated for transient and ac electroosmotic flows. The time evolution of the electroosmotic flow and the effect of a frequency-dependent ac electric field on the oscillating electroosmotic flow in a closed-end rectangular microchannel are examined. Specifically, the induced pressure gradient is analyzed under effects of the channel dimension and the frequency of electric field. In addition, based on the Stokes second problem, the solution of the slip velocity approximation is presented for comparison with the results obtained from the analytical scheme developed in this study.     

Simultaneous estimation of zeta potential and slip coefficient in hydrophobic microchannels

Electroosmotic flows through hydrophobic microchannels experience velocity slip at the channel wall, which increases the volumetric flow rate at a given electric potential gradient. The conventional method of zeta potential estimation using the volumetric flow rate may yield quite inaccurate zeta potential unless the velocity slip is appropriately taken care of. In the present investigation we develop a method for simultaneous estimation of zeta potential and velocity slip coefficient in the electroosmotic flow through a hydrophobic microchannel using velocity measurements. The relevant inverse problem is solved through the minimization of a performance function utilizing a conjugate gradient method. The present method is found to estimate the zeta potential and slip coefficient accurately even with noisy velocity measurements.     

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.     

Electroosmotic flow in channels with step changes in zeta potential and cross section

We consider the effects that step changes in zeta potential and cross section have on electroosmosis in long-and-narrow channels with arbitrary cross-sectional shapes. The Stokes equation of flow is solved analytically utilizing the thin Debye layer approximation to provide effective slip velocities on the channel walls. The effects of channel dimensions, surface potentials, applied pressure drop, and applied voltage are discussed. One anecdotal case, a two-region rectangular channel, is presented to illustrate the solution. The flow in each region is a combination of a uniform electroosmotic flow and a nonuniform pressure-driven flow. The electroosmotic pumping causes the pressure gradient in each region to adjust so that the flow rate is the same in each region and the overall applied pressure drop is met, resulting in convex velocity profiles in some regions and concave velocity profiles in other regions. By appropriate choice of the applied pressure drop, flat velocity profiles may be achieved in one or more regions.     

Analysis of electroosmotic flow of power-law fluids in a slit microchannel

Electroosmotic flow of power-law fluids in a slit channel is analyzed. The governing equations including the linearized Poisson-Boltzmann equation, the Cauchy momentum equation, and the continuity equation are solved to seek analytical expressions for the shear stress, dynamic viscosity, and velocity distribution. Specifically, exact solutions of the velocity distributions are explicitly found for several special values of the flow behavior index. Furthermore, with the implementation of an approximate scheme for the hyperbolic cosine function, approximate solutions of the velocity distributions are obtained. In addition, a generalized Smoluchowski velocity is introduced by taking into account contributions due to the finite thickness of the electric double layer and the flow behavior index of power-law fluids. Calculations are performed to examine the effects of kappaH, flow behavior index, double layer thickness, and applied electric field on the shear stress, dynamic viscosity, velocity distribution, and average velocity/flow rate of the electroosmotic flow of power-law fluids.     

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.     

调控双电层电渗流数值研究

在通道壁面垂直施加一个调控电场可以改变双电层电荷密度和Zeta电位势,实现对电渗流的调控.采用电场Poisson方程、动量守恒的Navier-Stokes方程、电解质离子输运的Nernst-Planck方程和液体混合反应的组分浓度输运方程,本文对微通道壁面离散布置调控电极的情况进行了数值分析.数值算例包括单电极、双电极...     

Electrokinetic-vortex formation near a two-part cylinder with same-sign zeta potentials in a straight microchannel

Vortex formation near a two-part cylinder with zeta potentials of different values but the same sign under an external DC electric field is numerically investigated in this paper. The cylinder, inserted in a straight microchannel filled with an aqueous solution, is composed of an upstream part and a downstream part. When a DC electric field is applied in the channel, under certain conditions, the vortex will form near the cylinder due to the different velocities of electroosmotic flow generated on the cylinder surface. The numerical results reveal that the larger the velocity difference of electroosmotic flow generated on the two-part cylinder and the smaller the channel width, the more conducive to vortex formation in the channel. In addition, if the zeta potential ratios of cylinder downstream part to upstream part and channel wall to cylinder upstream part are unchanged, the DC electric field strength and the zeta potential value do not affect the pattern of vortices formed in the channel. This study provides a way for vortex formation in microchannels and has the potential application in microfluidic devices.     

Electroosmotic Flow in Microchannels

The electroosmotic flow induced by an applied electrostatic potential field through microchannels between two parallel plates and a 90 degrees bend is analyzed in this work. A nonlinear, two-dimensional Poisson-Boltzmann equation governing the electrical double-layer field and the Laplace equation governing the electrostatic field distribution in microchannels are numerically solved using a finite-difference method. A body force caused by the interaction between the electrical double-layer field and the applied electrostatic field is included in the full Navier-Stokes equations. The effects of the electrical double-layer field and the applied electrostatic field on the fluid velocity distribution, pressure drop, and skin friction are discussed. A small pressure drop along the parallel plates is detected, although it is always neglected in the literature. Pressure is not a constant across the channel height. The axial velocity profile is no longer flat across the channel height when the Reynolds number is large. A separation bubble is detected near the 90 degrees junction when the Reynolds number is large. Copyright 2001 Academic Press.     

Remotely powered distributed microfluidic pumps and mixers based on miniature diodes

We demonstrate new principles of microfluidic pumping and mixing by electronic components integrated into a microfluidic chip. The miniature diodes embedded into the microchannel walls rectify the voltage induced between their electrodes from an external alternating electric field. The resulting electroosmotic flows, developed in the vicinity of the diode surfaces, were utilized for pumping or mixing of the fluid in the microfluidic channel. The flow velocity of liquid pumped by the diodes facing in the same direction linearly increased with the magnitude of the applied voltage and the pumping direction could be controlled by the pH of the solutions. The transverse flow driven by the localized electroosmotic flux between diodes oriented oppositely on the microchannel was used in microfluidic mixers. The experimental results were interpreted by numerical simulations of the electrohydrodynamic flows. The techniques may be used in novel actively controlled microfluidic-electronic chips.     

Numerical study of electroosmotic slip flow of fractional Oldroyd-B fluids at high zeta potentials

In this paper, an investigation of the electroosmotic flow of fractional Oldroyd-B fluids in a narrow circular tube with high zeta potential is presented. The Navier linear slip law at the walls is considered. The potential field is applied along the walls described by the nonlinear Poisson–Boltzmann equation. It's worth noting here that the linear Debye–Hückel approximation can't be used at the condition of high zeta potential and the exact solution of potential in cylindrical coordinates can't be obtained. Therefore, the Matlab bvp4c solver method and the finite difference method are employed to numerically solve the nonlinear Poisson–Boltzmann equation and the governing equations of the velocity distribution, respectively. To verify the validity of our numerical approach, a comparison has been made with the previous work in the case of low zeta potential and the excellent agreement between the solutions is clear. Then, in view of the obtained numerical solution for the velocity distribution, the numerical solutions of the flow rate and the shear stress are derived. Furthermore, based on numerical analysis, the influence of pertinent parameters on the potential distribution and the generation of flow is presented graphically.     

Electroosmotic flow in microchannels with arbitrary geometry and arbitrary distribution of wall charge

General solutions are developed for direct current (DC) and alternating current (AC) electroosmotic flows in microfluidic channels with arbitrary cross-sectional geometry and arbitrary distribution of wall charge (zeta potential). The applied AC electric field can also be of arbitrary waveform. By proposing a nondimensional time scale varpi defined as the ratio of the diffusion time of momentum across the electric double-layer thickness to the period of the applied electric field, we demonstrate analytically that the Helmholtz-Smoluchowski electroosmotic velocity is an appropriate slip condition for AC electroosmotic flows in typical microfluidic applications. With this slip condition approach, electroosmotic flows in rectangular and asymmetric trapezoidal microchannels with nonuniform wall charge, as examples, are investigated. The unknown constants in the proposed general solutions are numerically determined with a least-squares method through matching the boundary conditions. We find that the wall charge affects significantly the electroosmotic flow while the channel geometry does not. Moreover, the flow feature is characterized by another nondimensional time scale Omega defined as the ratio of the diffusion time of momentum across the channel hydraulic radius to the period of the applied electric field. The onset of phase shift between AC electroosmotic velocity and applied electric field is also examined analytically.     

Electrokinetic mixing vortices due to electrolyte depletion at microchannel junctions

Due to electric field leakage across sharp corners, the irrotational character of Ohmic electroosmotic flow is violated. Instead, we demonstrate experimentally and theoretically evidence of electrolyte depletion and vortex separation in electroosmotic flow around a junction between wide and narrow channels. When the penetration length of the electric field exceeds the width of the narrow channel and if the electric field is directed from the narrow to the wide channel, the electromigration of ions diminishes significantly at the junction end of the narrow channel due to this leakage. Concentration depletion then develops at that location to maintain current balance but it also increases the corner zeta potential and the local electroosmotic slip velocity. A back pressure gradient hence appears to maintain flow balance and, at a sufficient magnitude, generates a pair of vortices.     

Influence of the three-dimensional heterogeneous roughness on electrokinetic transport in microchannels

Surface roughness has been considered as a passive means of enhancing species mixing in electroosmotic flow through microfluidic systems. It is highly desirable to understand the synergetic effect of three-dimensional (3D) roughness and surface heterogeneity on the electrokinetic flow through microchannels. In this study, we developed a three-dimensional finite-volume-based numerical model to simulate electroosmotic transport in a slit microchannel (formed between two parallel plates) with numerous heterogeneous prismatic roughness elements arranged symmetrically and asymmetrically on the microchannel walls. We consider that all 3D prismatic rough elements have the same surface charge or zeta potential, the substrate (the microchannel wall) surface has a different zeta potential. The results showed that the rough channel's geometry and the electroosmotic mobility ratio of the roughness elements' surface to that of the substrate, epsilon(mu), have a dramatic influence on the induced-pressure field, the electroosmotic flow patterns, and the electroosmotic flow rate in the heterogeneous rough microchannels. The associated sample-species transport presents a tidal-wave-like concentration field at the intersection between four neighboring rough elements under low epsilon(mu) values and has a concentration field similar to that of the smooth channels under high epsilon(mu) values.     

Electroosmotic flow in a capillary annulus with high zeta potentials

The electroosmotic flow through an annulus is analyzed under the situation when the two cylindrical walls carry high zeta potentials. The analytical solutions for the electric potential profile and the electroosmotic flow field in the annulus are obtained by solving the Poisson-Boltzmann equation and the Stokes equation under an analytical scheme for the hyperbolic sine function. A mathematical expression for the average electroosmotic velocity is derived in a fashion similar to the Smoluchowski equation. Hence, a correction formula is introduced to modify the Smoluchowski equation, taking into account contributions due to the finite thickness of the electric double layer (EDL) and the geometry ratio-dependent correction. Specifically, under a circumstance when the two annular walls are oppositely charged, the flow direction can be determined from the sign of such correction formula, and there exists a zero-velocity plane inside the annulus. With the assumption of large electrokinetic diameters, the location of the zero-velocity plane can be estimated from the analytical expression for the velocity distribution. In addition, the characteristics of the electroosmotic flow through the annulus are discussed under the influences of the EDL parameters and geometric ratio of the inner radius to the outer radius of the annulus.     

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.