Transient electroosmotic flow induced by DC or AC electric fields in a curved microtube

This study investigates transient electroosmotic flow in a rectangular curved microtube in which the fluid is driven by the application of an external DC or AC electric field. The resultant flow-field evolutions within the microtube are simulated using the backwards-Euler time-stepping numerical method to clarify the relationship between the changes in the axial-flow velocity and the intensity of the applied electric field. When the electric field is initially applied or varies, the fluid within the double layer responds virtually immediately, and the axial velocity within the double layer tends to follow the varying intensity of the applied electric field. The greatest net charge density exists at the corners of the microtube as a result of the overlapping electrical double layers of the two walls. It results in local maximum or minimum axial velocities in the corners during increasing or decreasing applied electric field intensity in either the positive or negative direction. As the fluid within the double layer starts to move, the bulk fluid is gradually dragged into motion through the diffusion of momentum from the double layer. A finite time is required for the full momentum of the double layer to diffuse to the bulk fluid; hence, a certain phase shift between the applied electric field and the flow response is inevitable. The patterns of the axial velocity contours during the transient evolution are investigated in this study. It is found that these patterns are determined by the efficiency of momentum diffusion from the double layer to the central region of the microtube.     

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.     

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.     

Start‐up of electrophoresis of an arbitrarily oriented dielectric cylinder

An analytical study is presented for the transient electrophoretic response of a circular cylindrical particle to the step application of an electric field. The electric double layer adjacent to the particle surface is thin but finite compared with the radius of the particle. The time‐evolving electroosmotic velocity at the outer boundary of the double layer is utilized as a slip condition so that the transient momentum conservation equation for the bulk fluid flow is solved. Explicit formulas for the unsteady electrophoretic velocity of the particle are obtained for both axially and transversely applied electric fields, and can be linearly superimposed for an arbitrarily‐oriented applied field. If the cylindrical particle is neutrally buoyant in the suspending fluid, the transient electrophoretic velocity is independent of the orientation of the particle relative to the applied electric field and will be in the direction of the applied field. If the particle is different in density from the fluid, then the direction of electrophoresis will not coincide with that of the applied field until the steady state is attained. The growth of the electrophoretic mobility with the elapsed time for a cylindrical particle is substantially slower than for a spherical particle.     

Evaluation of electroosmotic drag coefficient of water in hydrated sodium perfluorosulfonate electrolyte polymer

The electroosmotic drag coefficient of water molecules in hydrated sodium perfluorosulfonate electrolyte polymer is evaluated on the basis of the velocity distribution functions of the sodium cations and water molecules with an electric field applied using molecular dynamics simulations. The simulation results indicate that both velocity distribution functions of water molecules and of sodium cations agree well with the classic Maxwellian velocity distribution functions when there is no electric field applied. If an electric field is applied, the distribution functions of velocity component in directions perpendicular to the applied electric field still agree with the Maxwellian velocity distribution functions but with different temperature parameters. In the direction of the applied electric field, the electric drag causes the velocity distribution function to deviate from the Maxwellian velocity distribution function; however, to obey the peak shifted Maxwellian distribution function. The peak shifting velocities coincide with the average transport velocities induced by the electric field, and could be applied to the evaluation of the electroosmotic drag coefficient of water. By evaluation of the transport velocities of water molecules in the first coordination shells of sodium cations, sulfonate anion groups, and in the bulk, it is clearly shown that the water molecules in the first coordination shell of sodium cations are the major contribution to the electroosmotic drag and momentum transfer from water molecules within the first coordination shell to the other water molecules also contributes to the electroosmotic drag. © 2008 Wiley Periodicals, Inc. J Comput Chem 2009     

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.     

Electrokinetic flow in a capillary with a charge-regulating surface polymer layer

An analytical study of the steady electrokinetic flow in a long uniform capillary tube or slit is presented. The inside wall of the capillary is covered by a layer of adsorbed or covalently bound charge-regulating polymer in equilibrium with the ambient electrolyte solution. In this solvent-permeable and ion-penetrable surface polyelectrolyte layer, ionogenic functional groups and frictional segments are assumed to distribute at uniform densities. The electrical potential and space charge density distributions in the cross section of the capillary are obtained by solving the linearized Poisson-Boltzmann equation. The fluid velocity profile due to the application of an electric field and a pressure gradient through the capillary is obtained from the analytical solution of a modified Navier-Stokes/Brinkman equation. Explicit formulas for the electroosmotic velocity, the average fluid velocity and electric current density on the cross section, and the streaming potential in the capillary are also derived. The results demonstrate that the direction of the electroosmotic flow and the magnitudes of the fluid velocity and electric current density are dominated by the fixed charge density inside the surface polymer layer, which is determined by the regulation characteristics such as the dissociation equilibrium constants of the ionogenic functional groups in the surface layer and the concentration of the potential-determining ions in the bulk solution.     

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.     

Diffusioosmosis and electroosmosis in a capillary slit with surface charge layers

This study analytically examines the steady diffusioosmotic and electroosmotic flows of an electrolyte solution in a fine capillary slit with each of its inside walls covered by a layer of adsorbed polyelectrolytes. In this solvent-permeable and ion-penetrable surface charge layer, idealized polyelectrolyte segments are assumed to distribute at a uniform density. The electric double layer and the surface charge layer may have arbitrary thicknesses relative to the gap width between the slit walls. The electrostatic potential distribution on a cross section of the slit is obtained by solving the linearized Poisson–Boltzmann equation, which applies to the case of low potentials or low fixed-charge densities. Explicit formulas for the fluid velocity profile due to the imposed electrolyte concentration gradient or electric field through the slit are derived as the solution of a modified Navier–Stokes/Brinkman equation. The results demonstrate that the structure of the surface charge layer can lead to an augmented or a diminished electrokinetic flow (even a reversal in direction of the flow) relative to that in a capillary with bare walls, depending on the characteristics of the capillary, of the surface charge layer, and of the electrolyte solution. For the diffusioosmotic flow with an induced electric field, competition between electroosmosis and chemiosmosis can result in more than one reversal in direction of the flow over a range of the Donnan potential of the adsorbed polyelectrolyte in the capillary.     

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.     

Electroosmotic flow in a water column surrounded by an immiscible liquid

In this paper, we conducted numerical simulation of the electroosmotic flow in a column of an aqueous solution surrounded by an immiscible liquid. While governing equations in this case are the same as that in the electroosmotic flow through a microchannel with solid walls, the main difference is the types of interfacial boundary conditions. The effects of electric double layer (EDL) and surface charge (SC) are considered to apply the most realistic model for the velocity boundary condition at the interface of the two fluids. Effects on the flow field of ?-potential and viscosity ratio of the two fluids were investigated. Similar to the electroosmotic flow in microchannels, an approximately flat velocity profile exists in the aqueous solution. In the immiscible fluid phase, the velocity decreases to zero from the interface toward the immiscible fluid phase. The velocity in both phases increases with ?-potential at the interface of the two fluids. The higher values of ?-potential also increase the slip velocity at the interface of the two fluids. For the same applied electric field and the same ?-potential at the interface of the two fluids, the more viscous immiscible fluid, the slower the system moves. The viscosity of the immiscible fluid phase also affects the flatness of the velocity profile in the aqueous solution.     

An AC electroosmotic micropump for circular chromatographic applications

Flow rates of up to 50 microm s(-1) have been successfully achieved in a closed-loop channel using an AC electroosmotic pump. The AC electroosmotic pump is made of an interdigitated array of unequal width electrodes located at the bottom of a channel, with an AC voltage applied between the small and the large electrodes. The flow rate was found to increase linearly with the applied voltage and to decrease linearly with the applied frequency. The pump is expected to be suitable for circular chromatography for the following reasons: the driving forces are distributed over the channel length and the pumping direction is set by the direction of the interdigitated electrodes. Pumping in a closed-loop channel can be achieved by arranging the electrode pattern in a circle. In addition the inherent working principle of AC electroosmotic pumping enables the independent optimisation of the channel height or the flow velocity.     

Tradeoff between mixing and transport for electroosmotic flow in heterogeneous microchannels with nonuniform surface potentials

Electroosmotic flow (EOF) is a phenomenon associated with the movement of an aqueous solution induced by the application of an electric field in microchannels. The characteristics of EOF depend on the nature of the surface potential, i.e., whether it is uniform or nonuniform. In this paper, a lattice Boltzmann model (LBM) combined with the Poisson-Boltzmann equation is used to simulate flow field in a rectangular microchannel with nonuniform (step change) surface potentials. The simulation results indicate that local circulations can occur near a heterogeneous region with nonuniform surface potentials, in agreement with those by other authors. Largest circulations, which imply a highest mixing efficiency due to convection and short-range diffusion, were found when the average surface potential is zero, regardless of whether the distribution of the heterogeneous patches is symmetric or asymmetric. In this work, we have illustrated that there is a trade-off between the mixing and liquid transport in EOF microfluidics. One should not simply focus on mixing and neglect liquid transport, as performed in the literature. Excellent mixing could lead to a poor transport of electroosmotic flow in microchannels.     

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.     

Molecular dynamics simulations of electroosmosis in perfluorosulfonic acid polymer

An atomistic MD simulation method has been developed to study the electroosmotic drag in the hydrated perfluorosulfonic acid polymer. The transport characteristics of the hydroniums and water molecules are evaluated from their velocity distribution functions with an electric field applied. It is shown that the microstructure of the hydrated perfluorosulfonic acid polymer is not perturbed significantly by the electric field up to 2 V/microm, and the velocity distribution functions obey the peak shifted Maxwell velocity distribution functions. The evaluated peak shifting velocities are only about 1% of the average thermal motion. The hydronium flow and water flow are evaluated from the average transport velocities or the peak shifting velocities. The electroosmotic drag coefficients from the MD simulations are in good correspondence with the experimental values. It is also shown that the electroosmotic drag coefficient has no or weak temperature dependence.     

Liquid transport in rectangular microchannels by electroosmotic pumping

The characteristics of electroosmotic flow in rectangular microchannels were investigated in this paper. A 2D Poisson–Boltzmann equation and the 2D momentum equation were used to model the electric double layer field and the flow field. The numerical solutions show significant influences of the channel cross-section geometry (i.e. the aspect ratio) on the velocity field and the volumetric flow rate. Also, the numerical simulation of the electroosmotic flow reveals how the velocity field and the volumetric flow rate depend on the ionic concentration, zeta potential, channel size and the applied electrical field strength.     

Electroosmotically enhanced mass transfer through polyacrylamide gels

We present an internal pumping strategy to enhance solute fluxes in polymer gels. The method is based on electroosmotic flow driven by an electric field applied across a gel that has been doped with charged colloidal inclusions. This work is motivated by the need to enhance the transport in gel-based biosensor devices whose response dynamics are often mass transfer limited. In this case, polyacrylamide gel slabs were doped with immobilized, charged silica colloids, and the flux of a fluorescent tracer was measured as a function of applied field strength, the volume fraction and size of the colloidal silica inclusions, and the bulk electrolyte composition. Significant flux enhancements were achieved with applied electric currents on the order of a few mA. Control experiments indicated that the flux enhancement was not due to any distortion of the gel diffusional properties in response to the presence of the inclusions. At a constant inclusion volume fraction, the electroosmotic solute flux enhancement was strongest for the smallest particle sizes that provide the highest total surface area, consistent with the electroosmotic mechanism whereby fluid flow is generated along the solid/liquid interface.     

Electrophoresis of a colloidal sphere in a spherical cavity with arbitrary zeta potential distributions

An analytical study is presented for the quasi-steady electrophoretic motion of a dielectric sphere situated at the center of a spherical cavity when the surface potentials are arbitrarily nonuniform. The applied electric field is constant, and the electric double layers adjacent to the solid surfaces are assumed to be much thinner than the particle radius and the gap width between the surfaces. The presence of the cavity wall causes three basic effects on the particle velocity: (1) the local electric field on the particle surface is enhanced or reduced by the wall; (2) the wall increases the viscous retardation of the moving particle; and (3) a circulating electroosmotic flow of the suspending fluid exists because of the interaction between the electric field and the charged wall. The Laplace and Stokes equations are solved analytically for the electric potential and velocity fields, respectively, in the fluid phase, and explicit formulas for the electrophoretic and angular velocities of the particle are obtained. To apply these formulas, one has to calculate only the monopole, dipole, and quadrupole moments of the zeta-potential distributions at the particle and cavity surfaces. It is found that the contribution from the electroosmotic flow developing from the interaction of the imposed electric field with the thin double layer adjacent to the cavity wall and the contribution from the wall-corrected electrophoretic driving force to the particle velocities can be superimposed as a result of the linearity of the problem.     

Electrokinetic flow and electric conduction of salt-free solutions in a capillary

The electrokinetic flow and accompanied electric conduction of a salt-free solution in the axial direction of a charged circular capillary are analyzed. No assumptions are made about the surface charge density (or surface potential) and electrokinetic radius of the capillary, which are interrelated. The Poisson–Boltzmann equation and modified Navier–Stokes equation are solved for the electrostatic potential distribution and fluid velocity profile, respectively. Closed-form formulas for the electroosmotic mobility and electric conductivity in the capillary are derived in terms of the surface charge density. The relative surface potential, electroosmotic mobility, and electric conductivity are monotonic increasing functions of the surface charge density and electrokinetic radius. However, the rises of the relative surface potential and electroosmotic mobility with an increase in the surface charge density are suppressed substantially when it is high due to the effect of counterion condensation. The analytical prediction that the electroosmotic mobility grows with increases in the surface charge density and electrokinetic radius agrees with the experimental results for salt-free solutions in circular microchannels in the literature.     

Numerical simulation of bubble dynamics in a micro-channel under a nonuniform electric field

A numerical method is used to simulate the motion and coalescence of air bubbles in a micro-channel under a nonuniform electric field. The channel is equipped with arrays of electrodes embedded in its wall and voltages are applied on the electrodes to generate a specified electric field gradient in the longitudinal direction. In the study, the Navier-Stokes equations are solved by using the level set method handling the deformable/moving interfaces between the bubbles and the ambient liquid. Both the polarization Coulomb force and the dielectrophoresis force are considered as the force source of the Navier-Stokes equations by solving the Maxwell's equations. The flow field equations and the electric field equations are coupled and solved by using the finite element method. The electric field characteristics and the dynamic behavior of a bubble are analyzed by studying the distributions of the electric field and the force, the deformation and the moving velocity of the air bubble. The result suggests that the model of dispersed drops suspended in the immiscible dielectric liquid and driven by a nonuniform electric field is an effective method for the transportation and coalescence of micro-drops.