Parametrical studies of electroosmotic transport characteristics in submicrometer channels

Spatially two-dimensional nonequilibrium mathematical model describing electroosmotic flow through a submicrometer channel with an electric charge fixed on the channel walls is presented. This system is governed by the hydrodynamic, electrostatic, and mass transport phenomena. The model is based on the coupled mass balances, Poisson, Navier-Stokes, and Nernst-Planck equations. Nonslip boundary conditions are employed. The effect of an imposed electric field on the system behavior is studied by means of a numerical analysis of the model equations. We have obtained the following findings. If the channel width is comparable to the thickness of the electric double layer, the system behaves as an ion-exchange membrane and the dependence of the electric current passing through the channel on the applied voltage is strongly nonlinear. In the case of negatively (positively) charged walls, a narrow region of very low conductivity (so-called ionic gate) is formed in the free electrolyte near the channel entry facing the anode (cathode) side. For a wide channel, the electric current is proportional to the applied voltage and the velocity of electrokinetic flow is linearly proportional to the electric field strength. Complex hydrodynamics (eddy formation and existence of ionic gates) is the most interesting characteristics of the studied system. Hence, current-voltage and velocity-voltage curves and the corresponding spatial distributions of the model variables at selected points are studied and described in detail.     

Non-orthogonal micro-free flow electrophoresis: From theory to design concept

Micro-free flow electrophoresis (μFFE) is a technique that facilitates continuous separation of molecules in a shallow channel with a hydrodynamic flow and an electric field at an angle to the flow. We recently developed a general theory of μFFE that suggested that an electric field non-orthogonal to the flow could improve resolution. Here, we used computer modeling to study resolution as a function of the electric field strength and the angle between the electric field and the hydrodynamic flow. In addition we used our general theory of μFFE to investigate other important influences on resolution, which include the velocity of the hydrodynamic flow, the height of the separation channel, and the magnitude and direction of the electroosmotic flow. Finally, we propose four designs that could be used to generate non-orthogonal electric fields and discuss their relative merits.     

The bis(p-sulfonatophenyl)phenylphosphine-assisted synthesis and phase transfer of ultrafine gold nanoclusters

The behavior of Couette flow of nanofluids composed of negatively-charged nanoparticles dispersed in aqueous NaCl solutions is studied theoretically. The equation for calculating the Couette flow velocity profiles is derived. The induced electric fields and velocity profiles are calculated as a function of key parameters including nanoparticle size and volume fraction. We have found for the first time that the velocity profile of nanofluids containing charged nanoparticles deviates significantly from the classical linear velocity profile for Couette flow. This previously unseen flow phenomenon is attributed to the dominance of the electric field strength induced by the flow of charged nanoparticles. This new mechanism of nanoparticle-induced microfluidic transport could lead to novel microfluidic and tribological applications.     

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.     

Transient electroosmotic flow induced by AC electric field in micro-channel with patchwise surface heterogeneities

This paper investigates two-dimensional, time-dependent electroosmotic flow driven by an AC electric field via patchwise surface heterogeneities distributed along the micro-channel walls. The time-dependent flow fields through the micro-channel are simulated for various patchwise heterogeneous surface patterns using the backwards-Euler time stepping numerical method. Different heterogeneous surface patterns are found to create significantly different electrokinetic transport phenomena. The transient behavior characteristics of the generated electroosmotic flow are then discussed in terms of the influence of the patchwise surface heterogeneities, the direction of the applied AC electric field, and the velocity of the bulk flow. It is shown that the presence of oppositely charged surface heterogeneities on the micro-channel walls results in the formation of localized flow circulations within the bulk flow. These circulation regions grow and decay periodically in phase with the applied periodic AC electric field intensity. The location and rotational direction of the induced circulations are determined by the directions of the bulk flow velocity and the applied electric field.     

Diffusioosmosis of electrolyte solutions in a fine capillary slit

The steady diffusioosmotic flows of an electrolyte solution along a charged plane wall and in a capillary channel between two identical parallel charged plates generated by an imposed tangential concentration gradient are theoretically investigated. The plane walls may have either a constant surface potential or a constant surface charge density. The electrical double layers adjacent to the charged walls may have an arbitrary thickness and their electrostatic potential distributions are determined by the Poisson-Boltzmann equation. Solving a modified Navier-Stokes equation with the constraint of no net electric current arising from the cocurrent diffusion, electric migration, and diffusioosmotic convection of the electrolyte ions, the macroscopic electric field and the fluid velocity along the tangential direction induced by the imposed electrolyte concentration gradient are obtained semianalytically as a function of the lateral position in a self-consistent way. The direction of the diffusioosmotic flow relative to the concentration gradient is determined by the combination of the zeta potential (or surface charge density) of the wall, the properties of the electrolyte solution, and other relevant factors. For a given concentration gradient of an electrolyte along a plane wall, the magnitude of fluid velocity at a position in general increases with an increase in its electrokinetic distance from the wall, but there are exceptions. The effect of the lateral distribution of the induced tangential electric field and the relaxation effect in the double layer on the diffusioosmotic flow are found to be very significant.     

Deformation and motion of a charged conducting drop in a dielectric liquid under a nonuniform electric field

As a tool for transporting a drop inside another fluid, a charged conducting drop driven by Coulombic force is considered. Specifically, deformation and motion of a charged conducting drop under nonuniform electric fields are studied using the perturbation method. For simplicity in analysis, the applied electric field is assumed to be expressed as the sum of a uniform field and a linear field and the flow is assumed to be in the Stokes flow range. The deformed drop shape due to electrical stress is computed to the first order of the electrical Weber number (W). Then the electric force and the hydrodynamic drag are computed to derive the formula of the translation velocity, which is valid up to O(W). Several important results have also been obtained for the effect of drop deformation on the electric and hydrodynamic forces exerted on the drop.     

Boundary effects on electrophoresis of a colloidal cylinder with a nonuniform zeta potential distribution

The electrophoretic motion of a long dielectric circular cylinder with a general angular distribution of its surface potential under a transversely imposed electric field in the vicinity of a large plane wall parallel to its axis is analyzed. The thickness of the electric double layers adjacent to the solid surfaces is assumed to be much smaller than the particle radius and the gap width between the surfaces, but the applied electric field can be either perpendicular or parallel to the plane wall. The presence of the confining 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 viscous retardation of the moving particle; (3) an electroosmotic flow of the suspending fluid may exist due to the interaction between the charged wall and the tangentially imposed electric field. Through the use of cylindrical bipolar coordinates, the Laplace and Stokes equations are solved analytically for the two-dimensional electric potential and velocity fields, respectively, in the fluid phase, and explicit formulas for the quasisteady electrophoretic and angular velocities of the cylindrical particle are obtained. To apply these formulas, one has only to calculate the multipole moments of the zeta potential distribution at the particle surface. It is found that the existence of a plane wall near a nonuniformly charged particle can cause its translation or rotation which does not occur in an unbounded fluid with the same applied electric field.     

Analysis of electrokinetic transport of a spherical particle in a microchannel

Electrokinetically driven microfluidic devices that are used for biological cell/particle manipulation (e.g., cell sorting, separation) involve electrokinetic transport of these particles in microchannels whose dimension is comparable with particles' size. This paper presents an analytical study on electrokinetic transport of a charged spherical particle in a charged parallel-plate microchannel. Under the thin electric double-layer assumption, solutions in closed-form solutions for the particle velocity and disturbed electrical and fluid velocity fields are obtained for plane-symmetric (along the channel centerline) and asymmetric (off the channel centerline) motions of a sphere in a parallel-plate microchannel. The effects of relative particle size and eccentricity (i.e., off the centerline distance) on a particle's translational and rotational velocities are analyzed.     

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.     

Flow of an Anisotropic Dielectric Liquid in a Diverging Channel in a Strong Transverse Electric Field

A steady plane flow of an anisotropically polarizable liquid in a channel with nonparallel walls was considered. One of the walls was grounded, and the other was under a high electric potential. The polarization anisotropy was described in terms of a unit vector whose direction was determined by a relaxation equation. The dependence of the polarization of the liquid on the strength of the electric field and the anisotropy vector was specified using an equilibrium relation. Such a model can describe, for example, a suspension of anisotropically polarizable particles in a highly insulating liquid. The velocity, pressure, polarization, anisotropy vector, and electric field distributions in the liquid were determined and investigated. It was shown that, at some critical Reynolds number, backflows are initiated near the channel walls. The dependence of the critical Reynolds number on the diverging angle of the channel and on the properties of a liquid in a strong electric field was determined. The applied electric field increases the critical Reynolds number, which provides a means of controlling the regime of the considered flow using electrical methods.     

Transport of charged samples in fluidic channels with large zeta potentials

In this article, we present an analysis on the transport of charged samples through micro- and nanofluidic channels with large zeta potentials (|zeta| > (kBT)/e). Using the Method of Moments formulation, the diffusion-convection equation has been solved to evaluate the mean velocity and the dispersion of analyte bands in a parallel-plate device under electrokinetically- and pressure-driven flow conditions. The effect of electromigration induced by the lateral electric field within the Debye layer has been quantified in our work using a Peclet number (Pe t) based on the characteristic electrophoretic velocity of the solute molecules in the transverse direction. It has been shown that while the effects of transverse electromigration on analyte transport only depends on the product Pe t zeta* for |zeta*| = (ezeta)/kBT < 1, both these parameters independently affect the flow of charged species in large zeta potential systems. For a given value of Pe t zeta*, the mean velocity and the slug dispersivity can vary by as much as an order of magnitude in going from a small zeta potential system (|zeta*| < 1) to a channel with |zeta*| = 4.     

Voltage-controlled separation of proteins by electromobility focusing in a dialysis hollow fiber

Electromobility focusing (EMF) is a relatively new protein separation technique that utilizes an electric field gradient and a hydrodynamic flow. Proteins are focused in order of electrophoretic mobility at points where their electrophoretic migration velocities balance the hydrodynamic flow velocity. Steady state bands are formed along the separation channel when equilibrium is reached. Further separation and detection can be easily achieved by changing the electric field profile. In this paper. we describe an EMF system with on-line UV absorption detection in which the electric field gradient was formed using a dialysis hollow fiber. Protein focusing and preconcentration were performed with this system. Voltage-controlled separation was demonstrated using bovine serum albumin and myoglobin as model proteins. The limitations of the current method are discussed, and possible solutions are proposed.     

Prediction of permeate flux during osmotic pressure-controlled electric field-enhanced cross-flow ultrafiltration

Electric field-enhanced cross-flow ultrafiltration has been carried out to separate protein, bovine serum albumin, from aqueous solution using a 30,000 molecular weight cutoff membrane. A theoretical model is developed to predict permeate flux under a laminar flow regime including the effects of external d.c. electric field and suction through the membrane for osmotic pressure-controlled ultrafiltration. The governing equations of the concentration profile in the developing mass transfer boundary layer in a rectangular channel are solved using a similarity solution method. The effect of d.c. electric field on the variation of membrane surface concentration and permeate flux along the length of the channel is quantified using this model. The expression of Sherwood number relation for estimation of mass transfer coefficient is derived. The analysis revealed that there is a significant effect of electric field on the mass transfer coefficient. A detailed parametric study has been carried out to observe the effect of feed concentration, electric field, cross-flow velocity, and pressure on the permeate flux. For 1 kg/m3 BSA solution, by applying a d.c. electric field of 1000 V/m, the permeate flux increases from 42 to 98 L/m2 h compared to that with zero electric field. The experimental results are successfully compared with the model predicted results.     

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.     

Electrohydrodynamics of a drop under nonaxisymmetric electric fields

We consider the electrohydrodynamics of a spherical drop in a nonaxisymmetric electric field, which can be approximated by the sum of a uniform field and a linear straining field. We obtain the analytic solution of the three-dimensional flow fields inside and outside a drop for the Stokes flow regime by using Lamb's general solution and the leaky dielectric model. With the analytic solution, the dielectrophoretic migration velocity of a drop is obtained as a function of the type and the frequency of the imposed electric field. The direction of drop motion is found to be parallel to the dielectrophoretic force. The analytic solution is also used to investigate the characteristics of the interfacial flow under various nonaxisymmetric electric fields. While investigating the interfacial flow, we find a surface vortex structure under certain nonaxisymmetric electric fields, which is found to be related to the chaotic mixing inside the drop. Finally, we consider the chaotic features of three-dimensional flows inside the drop under static nonaxisymmetric electric fields.     

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.     

Electro‐osmotic flow in nanochannels with voltage‐controlled polyelectrolyte brushes: Dependence on grafting density and normal electric field

Molecular dynamics simulations were performed for electro‐osmotic flow (EOF) confined in a polyelectrolyte‐grafted nanochannel under variable grafting density and normal electric field. With decreasing the value of the normal electric field, the brush undergoes a collapse transition, and the ion distribution is changed significantly. The brush thickness increases on increasing the grafting density at positive and weak negative electric fields, whereas a reduced brush thickness is observed at strong negative electric field. Our results further reveal that the flow velocity is not only dependent on conformational transition of the brush but also related to the cation and anion distributions. At low grafting density, the EOF is almost completely quenched at high electric field strength due to strong surface friction between ions and walls. For the case of very dense grafting, the flow velocity is influenced weakly within the brush when varying the grafting density. Additionally, a bidirectional flow occurs at an intermediate electric field. The investigation on fluid flux indicates that the fluid flux is insensitive to the grafting density, when the normal electric field is removed. For nonzero normal electric fields, a significant change in the fluid flux is observed at low grafting densities. © 2012 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2012     

Simulation and experiment of asymmetric electrode placement for electrophoretic exclusion in a microdevice

Electrophoretic exclusion (EE) is a counterflow gradient technique that exploits hydrodynamic flow and electrophoretic forces to exclude, enrich, and separate analytes. Resolution for this technique has been theoretically examined and the smallest difference in electrophoretic mobilities that can be completely separated is estimated to be 10^?13 cm²/Vs. Traditional and mesoscale systems have been used, whereas microfluidics offers a greater range of geometries and configurations towards approaching this theoretical limit. To begin to understand the impact of seemingly subtle changes to the entrance flow and the electric field configurations, three closely related microfluidic interfaces were modeled, fabricated, and tested. These interfaces consisted of systematically varying placement of an asymmetric electrode relative to a channel entrance: leading electrode placed outside the channel entrance, leading electrode aligned with the channel, and leading electrode placed within the channel. A charged fluorescent dye is used as a sensitive and accurate probe for the model and to test the concentration variation at these interfaces. Models and experiments focused on visualizing the concentration profile in areas of high temporal dynamics, thus providing a severe test of the models. Experimental data and simulation results showed strong qualitative agreement. The complexity of the electric and flow fields about this interface and the agreement between models and testing suggests the theoretical assessment capabilities can be used to faithfully design novel and more efficient interfaces.     

Influence of an electric field on the non-Newtonian response of a hybrid-aligned nematic cell under shear flow

The authors study shear flow in hybrid-aligned nematic cells under the action of an applied electric field by solving numerically a hydrodynamic model. The authors apply this model to a flow-aligning nematic liquid crystal (4'-n-pentyl-4-cyanobiphenyl) and obtain the director's configuration and the velocity profile at the stationary state. The authors calculate the local and apparent viscosities of the system and found that the competition between the shear flow and the electric field gives rise to an interesting non-Newtonian response with regions of shear thickening and thinning. The results also show an important electrorheological effect ranging from a value a bit larger than the Miesowicz viscosity etab [Nature (London) 17, 261 (1935)] for small electric fields and large shear flows to etac for large electric fields and small shear flows. The analysis of the first normal stress difference shows that for small negative shear rates, the force between the plates of the cell is attractive, while it is repulsive for all other values of shear rates. However, under the application of the electric field, one can modify the extent of the region of attraction. Finally, the authors have calculated the dragging forces on the plates of the cell and found that it is easier to shear in one direction than in the other.