skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Nonlinear electrophoretic velocity of a spherical colloidal particle
Electrophoresis is the motion of a charged colloidal particle in an electrolyte under an applied electric field. The electrophoretic velocity of a spherical particle depends on the dimensionless electric field strength$$\beta =a^*e^*E_\infty ^*/k_B^*T^*$$, defined as the ratio of the product of the applied electric field magnitude$$E_\infty ^*$$and particle radius$$a^*$$, to the thermal voltage$$k_B^*T^*/e^*$$, where$$k_B^*$$is Boltzmann's constant,$$T^*$$is the absolute temperature, and$$e^*$$is the charge on a proton. In this paper, we develop a spectral element algorithm to compute the electrophoretic velocity of a spherical, rigid, dielectric particle, of fixed dimensionless surface charge density$$\sigma$$over a wide range of$$\beta$$. Here,$$\sigma =(e^*a^*/\epsilon ^*k_B^*T^*)\sigma ^*$$, where$$\sigma ^*$$is the dimensional surface charge density, and$$\epsilon ^*$$is the permittivity of the electrolyte. For moderately charged particles ($$\sigma ={O}(1)$$), the electrophoretic velocity is linear in$$\beta$$when$$\beta \ll 1$$, and its dependence on the ratio of the Debye length ($$1/\kappa ^*$$) to particle radius (denoted by$$\delta =1/(\kappa ^*a^*)$$) agrees with Henry's formula. As$$\beta$$increases, the nonlinear contribution to the electrophoretic velocity becomes prominent, and the onset of this behaviour is$$\delta$$-dependent. For$$\beta \gg 1$$, the electrophoretic velocity again becomes linear in field strength, approaching the Hückel limit of electrophoresis in a dielectric medium, for all$$\delta$$. For highly charged particles ($$\sigma \gg 1$$) in the thin-Debye-layer limit ($$\delta \ll 1$$), our computations are in good agreement with recent experimental and asymptotic results.  more » « less
Award ID(s):
2002120
PAR ID:
10502186
Author(s) / Creator(s):
;
Publisher / Repository:
Cambridge University Press
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
968
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Given a family$$\mathcal{F}$$of bipartite graphs, theZarankiewicz number$$z(m,n,\mathcal{F})$$is the maximum number of edges in an$$m$$by$$n$$bipartite graph$$G$$that does not contain any member of$$\mathcal{F}$$as a subgraph (such$$G$$is called$$\mathcal{F}$$-free). For$$1\leq \beta \lt \alpha \lt 2$$, a family$$\mathcal{F}$$of bipartite graphs is$$(\alpha,\beta )$$-smoothif for some$$\rho \gt 0$$and every$$m\leq n$$,$$z(m,n,\mathcal{F})=\rho m n^{\alpha -1}+O(n^\beta )$$. Motivated by their work on a conjecture of Erdős and Simonovits on compactness and a classic result of Andrásfai, Erdős and Sós, Allen, Keevash, Sudakov and Verstraëte proved that for any$$(\alpha,\beta )$$-smooth family$$\mathcal{F}$$, there exists$$k_0$$such that for all odd$$k\geq k_0$$and sufficiently large$$n$$, any$$n$$-vertex$$\mathcal{F}\cup \{C_k\}$$-free graph with minimum degree at least$$\rho (\frac{2n}{5}+o(n))^{\alpha -1}$$is bipartite. In this paper, we strengthen their result by showing that for every real$$\delta \gt 0$$, there exists$$k_0$$such that for all odd$$k\geq k_0$$and sufficiently large$$n$$, any$$n$$-vertex$$\mathcal{F}\cup \{C_k\}$$-free graph with minimum degree at least$$\delta n^{\alpha -1}$$is bipartite. Furthermore, our result holds under a more relaxed notion of smoothness, which include the families$$\mathcal{F}$$consisting of the single graph$$K_{s,t}$$when$$t\gg s$$. We also prove an analogous result for$$C_{2\ell }$$-free graphs for every$$\ell \geq 2$$, which complements a result of Keevash, Sudakov and Verstraëte. 
    more » « less
  2. Large-eddy simulation was used to model turbulent atmospheric surface layer (ASL) flow over canopies composed of streamwise-aligned rows of synthetic trees of height,$$h$$, and systematically arranged to quantify the response to variable streamwise spacing,$$\delta _1$$, and spanwise spacing,$$\delta _2$$, between adjacent trees. The response to spanwise and streamwise heterogeneity has, indeed, been the topic of a sustained research effort: the former resulting in formation of Reynolds-averaged counter-rotating secondary cells, the latter associated with the$$k$$- and$$d$$-type response. No study has addressed the confluence of both, and results herein show secondary flow polarity reversal across ‘critical’ values of$$\delta _1$$and$$\delta _2$$. For$$\delta _2/\delta \lesssim 1$$and$$\gtrsim 2$$, where$$\delta$$is the flow depth, the counter-rotating secondary cells are aligned such that upwelling and downwelling, respectively, occurs above the elements. The streamwise spacing$$\delta _1$$regulates this transition, with secondary cell reversal occurring first for the largest$$k$$-type cases, as elevated turbulence production within the canopy necessitates entrainment of fluid from aloft. The results are interpreted through the lens of a benchmark prognostic closure for effective aerodynamic roughness,$$z_{0,{Eff.}} = \alpha \sigma _h$$, where$$\alpha$$is a proportionality constant and$$\sigma _h$$is height root mean square. We report$$\alpha \approx 10^{-1}$$, the value reported over many decades for a broad range of rough surfaces, for$$k$$-type cases at small$$\delta _2$$, whereas the transition to$$d$$-type arrangements necessitates larger$$\delta _2$$. Though preliminary, results highlight the non-trivial response to variation of streamwise and spanwise spacing. 
    more » « less
  3. A lubrication analysis is presented for the resistances between permeable spherical particles in near contact,$$h_0/a\ll 1$$, where$$h_0$$is the minimum separation between the particles, and$$a=a_1 a_2/(a_1+a_2)$$is the reduced radius. Darcy's law is used to describe the flow inside the permeable particles and no-slip boundary conditions are applied at the particle surfaces. The weak permeability regime$$K=k/a^{2} \ll 1$$is considered, where$$k=\frac {1}{2}(k_1+k_2)$$is the mean permeability. Particle permeability enters the lubrication resistances through two functions of$$q=K^{-2/5}h_0/a$$, one describing axisymmetric motions, the other transverse. These functions are obtained by solving an integral equation for the pressure in the near-contact region. The set of resistance functions thus obtained provide the complete set of near-contact resistance functions for permeable spheres and match asymptotically to the standard hard-sphere resistances that describe pairwise hydrodynamic interactions away from the near-contact region. The results show that permeability removes the contact singularity for non-shearing particle motions, allowing rolling without slip and finite separation velocities between touching particles. Axisymmetric and transverse mobility functions are presented that describe relative particle motion under the action of prescribed forces and in linear flows. At contact, the axisymmetric mobility under the action of oppositely directed forces is$$U/U_0=d_0K^{2/5}$$, where$$U$$is the relative velocity,$$U_0$$is the velocity in the absence of hydrodynamic interactions and$$d_0=1.332$$. Under the action of a constant tangential force, a particle in contact with a permeable half-space rolls without slipping with velocity$$U/U_0=d_1(d_2+\log K^{-1})^{-1}$$, where$$d_1=3.125$$and$$d_2=6.666$$; in shear flow, the same expression holds with$$d_1=7.280$$. 
    more » « less
  4. An analysis is presented for the axisymmetric lubrication resistance between permeable spherical particles. Darcy's law is used to describe the flow in the permeable medium and a slip boundary condition is applied at the interface. The pressure in the near-contact region is governed by a non-local integral equation. The asymptotic limit$$K=k/a^{2} \ll 1$$is considered, where$$k$$is the arithmetic mean permeability, and$$a^{-1}=a^{-1}_{1}+a^{-1}_{2}$$is the reduced radius, and$$a_1$$and$$a_2$$are the particle radii. The formulation allows for particles with distinct particle radii, permeabilities and slip coefficients, including permeable and impermeable particles and spherical drops. Non-zero particle permeability qualitatively affects the axisymmetric near-contact motion, removing the classical lubrication singularity for impermeable particles, resulting in finite contact times under the action of a constant force. The lubrication resistance becomes independent of gap and attains a maximum value at contact$$F=6{\rm \pi} \mu a W K^{-2/5}\tilde {f}_c$$, where$$\mu$$is the fluid viscosity,$$W$$is the relative velocity and$$\tilde {f}_c$$depends on slip coefficients and weakly on permeabilities; for two permeable particles with no-slip boundary conditions,$$\tilde {f}_c=0.7507$$; for a permeable particle in near contact with a spherical drop,$$\tilde {f}_c$$is reduced by a factor of$$2^{-6/5}$$. 
    more » « less
  5. Abstract What proportion of integers$$n \leq N$$may be expressed as$$x^2 + dy^2$$for some$$d \leq \Delta $$, with$$x,y$$integers? Writing$$\Delta = (\log N)^{\log 2} 2^{\alpha \sqrt {\log \log N}}$$for some$$\alpha \in (-\infty , \infty )$$, we show that the answer is$$\Phi (\alpha ) + o(1)$$, where$$\Phi $$is the Gaussian distribution function$$\Phi (\alpha ) = \frac {1}{\sqrt {2\pi }} \int ^{\alpha }_{-\infty } e^{-x^2/2} dx$$. A consequence of this is a phase transition: Almost none of the integers$$n \leq N$$can be represented by$$x^2 + dy^2$$with$$d \leq (\log N)^{\log 2 - \varepsilon }$$, but almost all of them can be represented by$$x^2 + dy^2$$with$$d \leq (\log N)^{\log 2 + \varepsilon}\kern-1.5pt$$. 
    more » « less