Diffusiophoresis and electrophoresis of a charged sphere parallel to one or two plane walls

The diffusiophoretic and electrophoretic motions of a dielectric spherical particle in an electrolyte solution located between two infinite parallel plane walls are studied theoretically. The imposed electrolyte concentration gradient or electric field is constant and parallel to the two plates, which may be either impermeable to the ions/charges or prescribed with the far-field concentration/potential distribution. The electrical double layer at the particle surface is assumed to be thin relative to the particle radius and to the particle-wall gap widths, but the polarization effect of the mobile ions in the diffuse layer is incorporated. The presence of the neighboring walls causes two basic effects on the particle velocity: first, the local electrolyte concentration gradient or electric field on the particle surface is enhanced or reduced by the walls, thereby speeding up or slowing down the particle; second, the walls increase the viscous retardation of the moving particle. To solve the conservative equations, the general solution is constructed from the fundamental solutions in both rectangular and spherical coordinates. The boundary conditions are enforced first at the plane walls by the Fourier transforms and then on the particle surface by a collocation technique. Numerical results for the diffusiophoretic and electrophoretic velocities of the particle relative to those of a particle under identical conditions in an unbounded solution are presented for various values of the relevant parameters including the relative separation distances between the particle and the two plates. For the special case of motions of a spherical particle parallel to a single plate and in the central plane of a slit, the collocation results agree well with the approximate analytical solutions obtained by using a method of reflections. The presence of the lateral walls can reduce or enhance the particle velocity, depending on the properties of the particle-solution system, the relative particle-wall separation distances, and the electrochemical boundary condition at the walls. In general, the boundary effects on diffusiophoresis and electrophoresis are quite significant and complicated, and they no longer vary monotonically with the separation distances for some situations.     

Particle Interactions in Diffusiophoresis and Electrophoresis of Colloidal Spheres with Thin but Polarized Double Layers

The diffusiophoretic and electrophoretic motions of two colloidal spheres in the solution of a symmetrically charged electrolyte are analyzed using a method of reflections. The particles are oriented arbitrarily with respect to the electrolyte gradient or the electric field, and they are allowed to differ in radius and in zeta potential. The thickness of the electric double layers surrounding the particles is assumed to be small relative to the radius of each particle and to the gap width between the particles, but the effect of polarization of the mobile ions in the diffuse layer is taken into account. A slip velocity of fluid and normal fluxes of solute ions at the outer edge of the thin double layer are used as the boundary conditions for the fluid phase outside the double layers. The method of reflections is based on an analysis of the electrochemical potential and fluid velocity disturbances produced by a single dielectric sphere placed in an arbitrarily varying electrolyte gradient or electric field. The solution for two-sphere interactions is obtained in expansion form correct to O(r(12)(-7)), where r(12) is the distance between the particle centers. Our analytical results are found to be in good agreement with the available numerical solutions obtained using a boundary collocation method. On the basis of a model of statistical mechanics, the results of two-sphere interactions are used to analytically determine the first-order effect of the volume fraction of particles of each type on the mean diffusiophoretic and eletrophoretic velocities in a bounded suspension. For a suspension of identical spheres, the mean diffusiophoretic velocity can be decreased or increased as the volume fraction of the particles is increased, while the mean electrophoretic velocity is reduced with the increase in the particle concentration. Generally speaking, the particle interaction effects can be quite significant in typical situations. Copyright 2000 Academic Press.     

Diffusiophoresis in a suspension of charge-regulating colloidal spheres

An analytical study of diffusiophoresis in a homogeneous suspension of identical spherical charge-regulating particles with an arbitrary thickness of the electric double layers in a solution of a symmetrically charged electrolyte with a uniform prescribed concentration gradient is presented. The charge regulation due to association/dissociation reactions of ionogenic functional groups on the particle surface is approximated by a linearized regulation model, which specifies a linear relationship between the surface charge density and the surface potential. The effects of particle-particle electrohydrodynamic interactions are taken into account by employing a unit cell model, and the overlap of the double layers of adjacent particles is allowed. The electrokinetic equations that govern the electric potential profile, the ionic concentration distributions, and the fluid flow field in the electrolyte solution surrounding the particle in a unit cell are linearized assuming that the system is only slightly distorted from equilibrium. Using a regular perturbation method, these linearized equations are solved with the equilibrium surface charge density (or zeta potential) of the particle as the small perturbation parameter. Closed-form formulas for the diffusiophoretic velocity of the charge-regulating sphere correct to the second order of its surface charge density or zeta potential are derived. Our results indicate that the charge regulation effect on the diffusiophoretic mobility is quite sensitive to the boundary condition for the electric potential specified at the outer surface of the unit cell. For the limiting cases of a very dilute suspension and a very thin or very thick electric double layer, the particle velocity is independent of the charge regulation parameter.     

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.     

Electroviscous sphere-wall interactions

A theoretical analysis is presented to determine the forces of interaction between an electrically charged spherical particle and a charged plane wall when the particle translates parallel to the wall and rotates around its axis in a symmetric electrolyte solution at rest. The electroviscous effects, arising from the coupling between the electrical and hydrodynamic equations, are determined as a solution of three partial differential equations, derived from Cox's general theory [R.G. Cox, J. Fluid Mech. 338 (1997) 1], for electroviscous ion concentration, electroviscous potential and electroviscous flow field. It is a priori assumed that the double layer thickness surrounding each charged surfaces is much smaller than the particle size. Using the matched asymptotic expansion technique, the electroviscous forces experienced by the sphere are explicitly determined analytically for small particle-wall distances, but low and intermediate Peclet numbers.     

Diffusiophoresis of a nonuniformly charged sphere in an electrolyte solution

The diffusiophoresis of a rigid, nonuniformly charged spherical particle in an electrolyte solution is analyzed theoretically focusing on the influences of the thickness of double layer, the surface charge distribution, the effect of electrophoresis, and the effect of double-layer polarization. We show that the nonuniform charge distribution on the particle surface yields complicated effect of double-layer polarization, leading to interesting diffusiophoretic behaviors. For example, if the sign of the middle part of the particle is different from that of its left- and right-hand parts, then depending upon the charge density and the fraction of the middle part, the particle can move either to the high-concentration side or to the low-concentration side. Both the diffusiophoretic velocity and its direction can be manipulated by the distribution of the surface charge density. In particular, if the electrophoresis effect is significant, then those properties are governed by the averaged surface charge density of the particle. A dipolelike particle, where its left- (right-) hand half is negatively (positively) charged, always migrates toward the low-concentration (left-hand) side, that is, it has a negative diffusiophoretic velocity. In addition, that diffusiophoretic velocity has a negative local minimum as the thickness of double layer varies.     

Transient electrophoresis of dielectric spheres

The dynamic electrophoretic response of a spherical dielectric particle suspended in an electrolyte solution to a step change in the applied electrics field is analytically studied. The electrical double layer surrounding the particle may have either a small but finite thickness or a very large thickness relative to the particle radius. For the case of electrophoresis of a particle with a thin double layer, the local electroosmotic velocity at the outer edge of the double layer evolving with time after the external field is imposed is used as an apparent slip boundary condition at the particle surface so that the unsteady equation of motion for the fluid flow outside the double layer is solved. Closed-form formulas for the transient electrophoretic mobility of the particle are derived as functions of relevant parameters. The results demonstrate that, when the double layer surrounding the particle is relatively thin, the normalized electrophoretic mobility at a given dimensionless time decreases monotonically with a decrease in the parameter kappaa, where kappa(-1) is the Debye screening length and a is the particle radius. When the double layer of the particle is relatively thick, the particle mobility can have magnitudes comparable to those for a particle with a thin double layer in the initial stage, but will become much smaller afterward. In general, the effect of the relaxation time for transient electrophoresis is negligible, regardless of the value of kappaa.     

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.     

Electrophoretic Motion of a Spherical Particle with a Thick Double Layer in Bounded Flows

The electrophoretic mobility of a spherical colloidal particle with low zeta potential near a solid charged boundary is calculated numerically for arbitrary values of the double layer thickness by a generalization of Teubner's method to the case of bounded flow. Three examples are considered: a sphere near a nonconducting planar wall with electric field parallel to the wall, near a perfectly conducting planar wall with electric field perpendicular to the wall, and on the axis of a cylindrical pore with electric field parallel to the axis. The results are compared with recent analytical calculations using the method of reflections. For the case of a charged sphere near a neutral surface, the reflection results are quite good, provided there is no double layer overlap, in which case there can be extra effects for constant potential particles that are entirely missed by the analytical expressions. For a neutral sphere near a charged surface, the reflection results are less successful. The main reason is that the particle feels the profile of the electroosmotic flow, an effect ignored by construction in the method of reflections. The general case is a combination of these, so that the reflections are more reliable when the electrophoretic motion dominates the electroosmotic flow. The effect on particle mobility of particle-wall interactions follows the trend expected on geometric grounds in that sphere-plane interactions are stronger than sphere-sphere interactions and the effect on a sphere in a cylindrical pore is stronger still. In the latter case, particle mobility can fall by more than 50% for thick double layers and a sphere half the diameter of the pore. The agreement between numerical results and analytical results follows the same trend, being worst for the sphere in a pore. Nevertheless, the reflections can be reliable for some geometries if there is no double layer overlap. This is demonstrated for a specific example where reflection results have previously been compared with experiments on protein mobility through a membrane (J. Ennis et al., 1996, J. Membrane Sci. 119, 47). Copyright 1999 Academic Press.     

Diffusioosmosis and electroosmosis of electrolyte solutions in fibrous porous media

The steady diffusioosmotic and electroosmotic flows of an electrolyte solution in the fibrous porous medium constructed by a homogeneous array of parallel charged circular cylinders are analyzed under conditions of small Peclet and Reynolds numbers. The imposed electrolyte concentration gradient or electric field is constant and can be oriented arbitrarily with respect to the axes of the cylinders. The thickness of the electric double layers surrounding the cylinders is assumed to be small relative to the radius of the cylinders and to the gap width between two neighboring cylinders, but the polarization effect of the diffuse ions in the double layers is incorporated. Through the use of a unit cell model, the appropriate equations of conservation of the electrochemical potential energies of ionic species and the fluid momentum are solved for each cell, in which a cylinder is envisaged to be surrounded by a coaxial shell of the fluid. Analytical expressions for the diffusioosmotic and electroosmotic velocities of the bulk electrolyte solution as functions of the porosity of the ordered array of cylinders are obtained in closed form for various cases. Comparisons of the results of the cell model with different conditions at the outer boundary of the cell are made. In the limit of maximum porosity, these results can be interpreted as the diffusiophoretic and electrophoretic velocities of an isolated circular cylinder caused by the imposed electrolyte concentration gradient or electric field.     

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.     

Electrophoretic mobility of a charged spherical colloidal particle covered with an uncharged polymer layer

A general expression is derived for the electrophoretic mobility of a spherical charged colloidal particle covered with an uncharged polymer layer in an electrolyte solution in an applied electric field for the case where the particle zeta potential is low. It is assumed that electrolyte ions as well as water molecules can penetrate the polymer layer. Approximate analytic expressions for the electrophoretic mobility of particles carrying low zeta potentials are derived for the two extreme cases in which the particle radius is very large or very small.     

Transient electrophoresis of spherical particles at low potential and arbitrary double-layer thickness

A theoretical study is presented for the dynamic electrophoretic response of a charged spherical particle in an unbounded electrolyte solution to a step change in the applied electric field. The electric double layer surrounding the particle may have an arbitrary thickness relative to the particle radius. The transient Stokes equations modified with the electrostatic effect which govern the fluid velocity field are linearized by assuming that the system is only slightly distorted from equilibrium. Semianalytical results for the transient electrophoretic mobility of the particle are obtained as a function of relevant parameters by using the Debye-Huckel approximation. The results demonstrate that the electrophoretic mobility of a particle with a constant relative mass density at a specified dimensionless time normalized by its steady-state quantity decreases monotonically with a decrease in the parameter kappaa, where kappa(-1) is the Debye screening length and a is the particle radius. For a given value of kappaa, a heavier particle lags behind a lighter one in the development of the electrophoretic mobility. In the limits of kappaa --> infinity and kappaa = 0, our results reduce to the corresponding analytical solutions available in the literature. The electrophoretic acceleration of the particle is a monotonic decreasing function of the time for any fixed value of kappaa. In practical applications, the effect of the relaxation time for the transient electrophoresis is negligible, regardless of the value of kappaa or the relative mass density of the particle.     

Boundary Effects on Diffusiophoresis of Cylindrical Particles in Nonelectrolyte Gradients

The diffusiophoretic motion of a long circular cylinder in a transversely imposed gradient of a nonionic solute near a large plane wall parallel to its axis is analyzed. The range of the interaction between the solute and the solid surfaces is assumed to be small relative to the particle radius and to the gap width between the particle and the wall, but the polarization effect of the mobile solute in the thin diffuse layers adjacent to the solid surfaces caused by the strong adsorption of the solute is incorporated. A normal flux of the solute and a slip velocity of the fluid at the outer edge of the diffuse layers are used as the boundary conditions for the fluid domain outside the diffuse layers. Through the use of cylindrical bipolar coordinates along with these boundary conditions, a set of transport equations governing this problem is solved in the quasisteady situation and the wall effects on the diffusiophoresis of the cylinder are computed for various cases. For the diffusiophoretic motion of a cylinder normal to a plane, the particle mobility decreases monotonically with the decrease of the distance of the particle axis from the wall. The stronger the polarization effect in the diffuse layer, the weaker the wall effect on the diffusiophoresis. The effect of the normal plane on the diffusiophoresis of a cylinder is much more significant than that for a sphere at the same separation. For the diffusiophoresis of a cylinder parallel to a plane, the boundary effect is a complicated function of the relevant parameters (not necessarily varies monotonically with the extent of separation) mainly due to the existence of a diffusio-osmotic flow caused by the tangential fluid velocity at the plane wall. Copyright 2000 Academic Press.     

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.     

Diffusiophoresis of an Elongated Cylindrical Nanoparticle along the Axis of a Nanopore

The translation of a charged, elongated cylindrical nanoparticle along the axis of a nanopore driven by an imposed axial salt concentration gradient is investigated using a continuum theory, which consists of the ionic mass conservation equations for the ionic concentrations, the Poisson equation for the electric potential in the solution, and the modified Stokes equations for the hydrodynamic field. The diffusiophoretic motion is driven by the induced electrophoresis and chemiphoresis. The former is driven by the generated overall electric field arising from the difference in the ionic diffusivities and the double layer polarization, while the latter is generated by the induced osmotic pressure gradient around the charged particle. The induced diffusiophoretic motion is investigated as functions of the imposed salt concentration gradient, the ratio of the particle’s radius to the double layer thickness, the cylinder’s aspect ratio (length/radius), the ratio of the nanopore size to the particle size, the surface charge densities of the nanoparticle and the nanopore, and the type of the salt used. The induced diffusiophoretic motion of a nanorod in an uncharged nanopore is mainly governed by the induced electrophoresis, driven by the induced electric field arising from the double layer polarization. The induced particle motion is driven by the induced electroosmotic flow, if the charges of the nanorod and nanopore wall have the same sign.     

Diffusiophoresis of a moderately charged spherical colloidal particle

An analytic expression is obtained for the diffusiophoretic mobility of a charged spherical colloidal particle in a symmetrical electrolyte solution. The obtained expression, which is expressed in terms of exponential integrals, is correct to the third order of the particle zeta potential so that it is applicable for colloidal particles with low and moderate zeta potentials at arbitrary values of the electrical double-layer thickness. This is an improvement of the mobility formula derived by Keh and Wei, which is correct to the second order of the particle zeta potential. This correction, which is related to the electrophoresis component of diffusiophoresis, becomes more significant as the difference between the ionic drag coefficients of electrolyte cations and anions becomes larger and vanishes in the limit of thin or thick double layer. A simpler approximate mobility expression is further obtained that does not involve exponential integrals.     

Diffusiophoresis in suspensions of highly charged soft particles

Diffusiophoresis phenomenon of aoft particles suspended in binary electrolyte solutions is explored theoretically in this study based on the spherical cell model, focusing on the chemiphoresis component in absence of diffusion potential. Both the electrostatic and hydrodynamic aspects of the boundary confinement, or steric effect, due to the presence of neighboring particles are examined extensively under various electrokinetic conditions. Significant local extrema are found in mobility profiles expressed as functions of the Debye length in general, synchronized with the strength of the motion-inducing double layer polarization. Moreover, a seemingly peculiar phenomenon is observed that the soft particles may move faster in more concentrated suspensions. The competition between the simultaneous enhancement of the motion-inducing electric driving force and the motion-retarding hydrodynamic drag force from the boundary confinement effect of the neighboring particles is found to be responsible for it. The above findings are also demonstrated experimentally in a very recent study on the diffusiophoretic motion of soft particles through porous collagen hydrogels. The results presented here are useful in various practical applications of soft particles like drug delivery.     

Electroviscous cylinder-wall interactions

A theoretical analysis is presented to determine the forces of interaction between an electrically charged cylindrical particle and a charged plane boundary wall when the particle translates parallel to the wall and rotates around its axis in a symmetric electrolyte solution at rest. The electroviscous effects, arising from the coupling between the electrical and hydrodynamic equations, are determined as a solution of three partial differential equations, derived from R.G. Cox's general theory [J. Fluid Mech. 338 (1997) 1], for electroviscous ion concentration, electroviscous potential, and electroviscous flow field. It is assumed a priori that the double layer thickness surrounding each charged surface is much smaller than the length scale of the problem. Using the matched asymptotic expansion technique, the electroviscous forces experienced by the cylinder are explicitly determined analytically for small particle-wall distances for low and intermediate Peclet numbers. It is found that the tangential force usually increases the drag above the purely hydrodynamic drag, although for certain conditions the drag can be reduced. Similarly the normal force is usually repulsive, i.e., it is an electrokinetic lift force, but under certain conditions the normal force can be attractive.     

Photophoresis of an aerosol sphere normal to a plane wall

A combined analytical-numerical study is presented for the quasisteady photophoretic motion of a spherical aerosol particle of arbitrary thermal conductivity and surface properties exposed to a radiative flux perpendicular to a large plane wall. The Knudsen number is assumed to be so small that the fluid flow is described by a continuum model with a temperature jump, a thermal slip, and a frictional slip at the surface of the radiation-absorbing particle. In the limit of small Peclet and Reynolds numbers, the appropriate equations of conservation of energy and momentum for the system are solved using a boundary collocation method and numerical results for the photophoretic velocity of the particle are obtained for various cases. The presence of the neighboring wall causes two basic effects on the particle velocity: first, the local temperature gradient on the particle surface is enhanced or reduced by the wall, thereby speeding up or slowing down the particle; second, the wall increases viscous retardation of the moving particle. The net effect of the wall can decrease or increase the particle velocity, depending upon the relative conductivity and surface properties of the particle as well as the relative particle-wall separation distance. In general, the boundary effect of a plane wall on the photophoresis of an aerosol particle can be quite significant in some situations. In most aerosol systems, the boundary effect on photophoresis is weaker than that on the motion driven by a gravitational field.