Analysis of electrokinetic pumping efficiency in capillary tubes

A mathematical model for predicting the maximum pumping efficiency and pressure difference generation by an electrokinetic-driven fluid pumping system through a capillary tube is presented in this study. Both the maximum pumping efficiency and optimum pressure difference generation are found to depend on a single variable. This single variable is termed as the figure of merit since it determines the performance of electrokinetic pumping. The figure of merit is found to depend on three dimensionless parameters, the normalized Debye length, zeta potential, and Levine number indicating the nominal ratio of convective current to conductive current. All three parameters can be related to the pH value and concentration of aqueous salt solution by the introduction of concentration-dependent electrical conductivity and pH-dependent zeta potential. By presenting the maximum pumping efficiency and optimum pressure difference generation as functions of pH value, salt concentration, and capillary tube radius, it is found that both maximum pumping efficiency and optimum pressure difference generation increase with the decrease in capillary radius and salt concentration. The optimum pH values at which the maximum pumping efficiency and optimum pressure difference generation occur are found to be in the range between 6 and 9. For the salt concentration of 10(-6) M, pH 6.9, and a capillary tube radius value of 0.5 micro m, the predicted maximum pumping efficiency is 5.4% which is close to the experimental measurement reported in the literature.     

Calculation of the dynamic impedance of the double layer on a planar electrode by the theory of electrokinetics

Applications of microelectromechanical systems in the biotechnological arena (bioMEMS) are a subject of great current interest. Accurate calculation of electric field distribution in these devices is essential to the understanding and design of processes such as dielectrophoresis and AC electroosmosis that drive MEMS-based devices. In this paper, we present the calculation of the electrical double-layer impedance (Z(el)) of an ideally polarizable plane electrode using the standard model of colloidal electrokinetics. The frequency variation of the electrical potential drop across the double layer above a planar electrode in a general electrolyte solution is discussed as a function of the electrode zeta potential zeta, the Debye length kappa(-1), the electrolyte composition and the bulk region thickness L.     

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.     

Streaming potential in cylindrical pores of poly(ethylene terephthalate) track-etched membranes: variation of apparent zeta potential with pore radius

Streaming potential variation with pressure measured through poly(ethylene terephthalate) track-etched membranes of different pore sizes led to the determination of an apparent interfacial potential zetaa in the presence of 10-2 M KCl. The variation of zetaa with the pore radius r0 is interpreted by (i) the electric double layer overlap effect and (ii) the presence of a conductive gel layer. We propose a method which integrates both effects by assuming a simple model for the conductive gel at the pore wall. We observed the following three domains of pore size: (i) r0 > 70 nm, where surface effects are negligible; (ii) approximately 17 nm < r0 < 70 nm, where the pore/solution interface could be described as a conductive gel of thickness around 1 nm; (iii) r0 < approximately 17 nm, which corresponds to the region strongly damaged by the ion beam and is not analyzed here. The first one (zeta = -36.2 mV) corresponds to the raw material when etching has completely removed the ion beam predamaged region, which corresponds to the second intermediate domain (zeta = -47.3 mV). There the conductance of the gel layer deduced from the treatment of streaming potential data was found to be compatible with the number of ionic sites independently determined by the electron spin resonance technique.     

Conditions for similitude and the effect of finite Debye length in electroosmotic flows

Under certain conditions, the velocity field is similar to the electric field for electroosmotic flow (EOF) inside a channel. There was a disagreement between investigators on the necessity of the infinitesimal-Reynolds-number condition for the similarity when the Helmholtz-Smoluchowski relation is applied throughout the boundaries. What is puzzling is a recent numerical result that showed, contrary to the conventional belief, an evident Reynolds number dependence of the EOF. We show here that the notion that the infinitesimal-Reynolds-number condition is required originates from the misunderstanding that the EOF is the Stokes flow. We point out that the EOF becomes the potential flow when the Helmholtz-Smoluchowski relation is applied at the boundaries. We carry out a numerical simulation to investigate the effect of finiteness of the Debye length and the vorticity layer inherently existing at the channel wall. We show that the Reynolds number dependence of the previous numerical simulation resulted from the finiteness of the Debye length and subsequent convective transport of vorticity toward the bulk flow. We discuss in detail how the convection of vorticity occurs and what factors are involved in the transport process, after carrying out the simulation for different Reynolds numbers, Debye lengths, corner radii, and geometries.     

Effects of overlapping electric double layer on mass transport of a macro‐solute across porous wall of a micro/nanochannel for power law fluid

Effects of overlapping electric double layer and high wall potential on transport of a macrosolute for flow of a power law fluid through a microchannel with porous walls are studied in this work. The electric potential distribution is obtained by coupling the Poisson's equation without considering the Debye–Huckel approximation. The numerical solution shows that the center line potential can be 16% of wall potential at pH 8.5, at wall potential −73 mV and scaled Debye length 0.5. Transport phenomena involving mass transport of a neutral macrosolute is formulated by species advective equation. An analytical solution of Sherwood number is obtained for power law fluid. Effects of fluid rheology are studied in detail. Average Sherwood number is more for a pseudoplastic fluid compared to dilatant upto the ratio of Poiseuille to electroosmotic velocity of 5. Beyond that, the Sherwood number is independent of fluid rheology. Effects of fluid rheology and solute size on permeation flux and concentration of neutral solute are also quantified. More solute permeation occurs as the fluid changes from pseudoplastic to dilatant.     

Analysis of the response of suspended colloidal soft particles to a constant electric field

A network model, originally designed for an electrokinetic study of soft particle suspensions, has been used for an in-depth analysis of the physical behavior of these systems under the action of an externally applied DC electric field. The versatility of the network simulation method used makes it possible to obtain information readily not only about the electrophoretic mobility, but also about any physical variable of interest at all points around the suspended particle: electric potential, ion concentrations, fluid velocity. The field-induced polarization of the double layer is described in terms of the dependence of these and other derived variables (volume charge density, electric field components, ion flux components) on the distance to the membrane-solution interface. In contrast to colloidal suspensions of hard particles, which basically depend on just two parameters (the reciprocal Debye length multiplied by the particle radius, kappaa, and the zeta potential, zeta), soft particle suspensions require a wider parameter set. First, there are two characteristic diffusion lengths in the system (one inside the membrane and the other in the solution) and two geometrical lengths (the core radius a and the membrane thickness (b-a)). Furthermore, there is the fixed charge density inside the membrane (and possibly a surface charge density over the core) that cannot be represented by a zeta potential. Finally, the parameter that characterizes the interaction between the fluid and the permeable membrane, gamma, strongly influences the behavior of the system. Dependences on all these parameters (except the geometrical ones) are included in this study.     

Determination of the dynamic electrophoretic mobility of a spherical colloidal particle through a novel numerical solution of the electrokinetic equations

The standard equations developed to describe the electrophoretic motion of a charged particle immersed in an electrolyte subjected to an oscillating electric field are solved numerically with a new technique suitable for stiff systems. The focus of this work is to use this solution to determine the dynamic particle mobility, one of several quantities that can be extracted from these equations. This solution is valid from low frequencies to indefinitely high frequencies and has no restriction on zeta potential, double-layer thickness, or electrolyte composition. The solution has been used to calculate the dynamic electrophoretic mobility of a particle for a wide range of double-layer thicknesses and zeta potentials. The solution agrees with analytic approximations obtained previously by other authors under the conditions of a thin double layer and low zeta potential. The results are also consistent with calculations valid at frequencies where the ion diffusion length extends a significant distance beyond the double layer as obtained by another numerical technique.     

Analysis of electroosmotic flow with step change in zeta potential

The term electroosmotic flow refers to the bulk flow of an aqueous solution induced by the application of the electric field to the zeta potential. The characteristics of EOF in a microchannel depend upon the nature of the zeta potential, i.e., whether it is uniform or nonuniform. In this study, the full Navier-Stokes equation and the Nernst-Planck equation are used to model the change in EOF characteristics that occur when a step change in zeta potential is applied. It is found that the thickness of the electrical double layer gradually increases downstream from the location at which the zeta potential is increased. The results indicate that a step change in zeta potential causes a significant variation in the velocity profile and in the pressure distribution.     

Changes in zeta potential caused by a dc electric current for thin double layers

An asymptotic solution was obtained to describe one-dimensional, steady-state transport of a symmetric binary electrolyte normal to two large parallel electrodes, in the limit in which the Debye length is infinitesimal compared to the distance separating the two electrodes. Despite the nonzero ion flux, Boltzmann's equation continues to describe the relationship between either ion concentration and the electrostatic potential inside the diffuse part of the double layer, while local electroneutrality applies outside, even for current densities approaching the limiting value. In the absence of ion adsorption or dissociation reactions at the electrodes, the magnitude of any charge or zeta potential arising on the electrodes at zero current is determined by the equilibrium constant for the redox reactions which would exchange ionic charge carriers for electric charge carriers at the electrode surface. Nonzero current causes the ionic strength of the bulk to vary with position. This perturbs the Debye length of the diffuse cloud on either electrode: it is the local ionic strength just outside the cloud which determines the Debye length for that cloud. Nonzero current also changes the zeta potential. The dimensionless rate of change dζ/dJ was as large as 30.     

Zeta potential of ion-conductive membranes by streaming current measurements

Surface charge properties have a significant influence on membrane retention and fouling performance. As a key parameter describing the surface charge of membranes used in aqueous applications, zeta potential measurements on membranes of various types have attracted great attention. During the zeta potential characterization of a series of ion-conductive sulfonated poly(sulfone) membranes, it was found that the measured streaming current varied with the thickness of the sample, which is not predicted by the classical Smoluchowski equation. Moreover, for higher conductivity membranes with an increased concentration of sulfonate groups, the zeta potential tended toward zero. It was determined that the influence of membrane bulk conductance on the measured streaming current must be taken into account in order to correctly interpret the streaming current data for ion-conductive polymers and understand the relationship between membrane chemical composition and zeta potential. Extrapolating the measured streaming current to a membrane thickness of zero has proven to be a feasible method of eliminating the error associated with measuring the zeta potential on ion conductive polymer membranes. A linear resistance model is proposed to account for the observed streaming currents where the electrolyte channel is in parallel with the ion-conductive membranes.     

Diffusiophoresis of concentrated suspensions of spherical particles with distinct ionic diffusion velocities

The diffusiophoresis of a concentrated spherical dispersion of colloidal particles subject to a small electrolyte gradient is analyzed theoretically for an arbitrary zeta potential and double layer thickness. In particular, the influence of the difference in the diffusivities of cations and anions is discussed. A unit cell model is used to simulate a spherical dispersion, and a pseudospectral method is adopted to solve the equations governing the phenomenon under consideration. We show that, as in the case of an infinitely dilute dispersion, when the diffusivities of cations and anions are different, the diffusiophoretic mobility is no longer an even function of the zeta potential or double layer thickness. In contrast to the case of identical diffusivity of cations and anions, a local electric field is induced in the present case due to an unbalanced charge distribution between higher and lower concentration regions. Depending upon the direction of this induced electric field, the diffusiophoretic mobility can be larger or smaller than that for the case of identical diffusivity. The diffusiophoretic mobility is influenced mainly by the induced electric field arising from the difference in the ionic diffusivities, the concentration gradient, and the effect of double layer polarization.     

Effect of double-layer polarization on the forces that act on a nanosized cylindrical particle in an ac electrical field

The polarization of, the forces acting on, and the electroosmotic flow field around a cylindrical particle of radius a* and uniform zeta potential zeta* submerged in an electrolyte solution and subjected to alternating electric fields are computed by solving the Poisson-Nernst-Planck (PNP) equations (the standard model). The dipole coefficient and the electrostatic and hydrodynamic forces are calculated as functions of the electric field's frequency, the solute concentration, and the particle's surface charge. The calculations are not restricted to small Debye screening lengths (lambdaD*). At relatively low frequencies, the polarization coefficient is nearly frequency-independent. As the frequency increases above D*/a*(2), where D* is the effective diffusion coefficient, the polarization coefficient initially increases, attains a maximum, and then decreases to an asymptotic value (when the frequency exceeds (1+Du)D*/lambdaD(*2), where Du is the Dukhin number). At low frequencies, when (lambdaD*/a*)(2)e(|zeta*F*/(2R*T*)|) < 1, the PNP calculations are in excellent agreement with the predictions of the Dukhin-Shilov (DS) low-frequency theory. At high frequencies, when lambda D*/a* < 1, the PNP calculations are in excellent agreement with the Maxwell-Wagner-O'Konski (MWO) theory.     

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.     

Diffusiophoresis of a cylindrical colloidal particle oriented parallel to an electrolyte concentration gradient field

We derive the general expression for the diffusiophoretic mobility of a cylindrical particle oriented parallel to an applied electrolyte concentration gradient field in a symmetrical electrolyte solution. From the general mobility expression as combined with an approximate analytic expression with negligible error for the electric potential distribution around a cylinder, an accurate analytic mobility expression is obtained, which is applicable for arbitrary values of the particle zeta potential and the electrical double layer thickness. It is also found that the low zeta potential approximation is an excellent approximation for low-to-moderate values of the particle zeta potential.     

Symmetrical radial distribution functions in the potentially theory of electrolyte solutions

A symmetrical radial distribution function g_ij in electric double layer theory has been proposed by Féat and Levine. We adopt their work to the electrolyte bulk to show how a symmetrical g_ij may be obtained in the potential theory of electrolyte solutions.     

Estimation of zeta potential by electrokinetic analysis of ionic fluid flows through a divergent microchannel

The streaming potential is generated by the electrokinetic flow effect within the electrical double layer of a charged solid surface. Surface charge properties are commonly quantified in terms of the zeta potential obtained by computation with the Helmholtz-Smoluchowski (H-S) equation following experimental measurement of streaming potential. In order to estimate a rigorous zeta potential for cone-shaped microchannel, the correct H-S equation is derived by applying the Debye-Hückel approximation and the fluid velocity of diverging flow on the specified position. The present computation provides a correction ratio relative to the H-S equation for straight cylindrical channel and enables us to interpret the effects of the channel geometry and the electrostatic interaction. The correction ratio decreases with increasing of diverging angle, which implies that smaller zeta potential is generated for larger diverging angle. The increase of Debye length also reduces the correction ratio due to the overlapping of the Debye length inside of the channel. It is evident that as the diverging angle of the channel goes to nearly zero, the correction ratio converges to the previous results for straight cylindrical channel.     

Electrophoresis and electric conduction in a salt-free suspension of charged particles

The electrophoresis and electric conduction of a suspension of charged spherical particles in a salt-free solution are analyzed by using a unit cell model. The linearized Poisson-Boltzmann equation (valid for the cases of relatively low surface charge density or high volume fraction of the particles) and Laplace equation are solved for the equilibrium electric potential profile and its perturbation caused by the imposed electric field, respectively, in the fluid containing the counterions only around the particle, and the ionic continuity equation and modified Stokes equations are solved for the electrochemical potential energy and fluid flow fields, respectively. Explicit analytical formulas for the electrophoretic mobility of the particles and effective electric conductivity of the suspension are obtained, and the particle interaction effects on these transport properties are significant and interesting. The scaled zeta potential, electrophoretic mobility, and effective electric conductivity increase monotonically with an increase in the scaled surface charge density of the particles and in general decrease with an increase in the particle volume fraction, keeping each other parameter unchanged. Under the Debye-Hückel approximation, the dependence of the electrophoretic mobility normalized with the surface charge density on the ratio of the particle radius to the Debye screening length and particle volume fraction in a salt-free suspension is same as that in a salt-containing suspension, but the variation of the effective electric conductivity with the particle volume fraction in a salt-free suspension is found to be quite different from that in a suspension containing added electrolyte.     

Motion of a colloidal sphere with interfacial self-electrochemical reactions induced by a magnetic field

The motion of a spherical colloidal particle with spontaneous electrochemical reactions occurring on its surface in an ionic solution subjected to an applied magnetic field is analyzed for an arbitrary zeta potential distribution. The thickness of the electric double layer adjacent to the particle surface is assumed to be much less than the particle radius. The solutions of the Laplace equations governing the magnetic scalar potential and electric potential, respectively, lead to the magnetic flux and electric current density distributions in the particle and fluid phases of arbitrary magnetic permeabilities and electric conductivities. The Stokes equations modified with the Lorentz force contribution for the fluid motion are dealt by using a generalized reciprocal theorem, and closed-form formulas for the translational and angular velocities of the colloidal sphere induced by the magnetohydrodynamic effect are obtained. The dipole and quadrupole moments of the zeta potential distribution over the particle surface cause the particle translation and rotation, respectively. The induced velocities of the particle are unexpectedly significant, and their dependence on the characteristics of the particle-fluid system is physically different from that for electromagnetophoretic particles or phoretic swimmers.