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: Resistance and mobility functions for the near-contact motion of permeable particles
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
Award ID(s):
1603806
PAR ID:
10558303
Author(s) / Creator(s):
;
Publisher / Repository:
Cambridge University Press
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
938
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract We study higher uniformity properties of the Möbius function$$\mu $$, the von Mangoldt function$$\Lambda $$, and the divisor functions$$d_k$$on short intervals$$(X,X+H]$$with$$X^{\theta +\varepsilon } \leq H \leq X^{1-\varepsilon }$$for a fixed constant$$0 \leq \theta < 1$$and any$$\varepsilon>0$$. More precisely, letting$$\Lambda ^\sharp $$and$$d_k^\sharp $$be suitable approximants of$$\Lambda $$and$$d_k$$and$$\mu ^\sharp = 0$$, we show for instance that, for any nilsequence$$F(g(n)\Gamma )$$, we have$$\begin{align*}\sum_{X < n \leq X+H} (f(n)-f^\sharp(n)) F(g(n) \Gamma) \ll H \log^{-A} X \end{align*}$$ when$$\theta = 5/8$$and$$f \in \{\Lambda , \mu , d_k\}$$or$$\theta = 1/3$$and$$f = d_2$$. As a consequence, we show that the short interval Gowers norms$$\|f-f^\sharp \|_{U^s(X,X+H]}$$are also asymptotically small for any fixedsfor these choices of$$f,\theta $$. As applications, we prove an asymptotic formula for the number of solutions to linear equations in primes in short intervals and show that multiple ergodic averages along primes in short intervals converge in$$L^2$$. Our innovations include the use of multiparameter nilsequence equidistribution theorems to control type$$II$$sums and an elementary decomposition of the neighborhood of a hyperbola into arithmetic progressions to control type$$I_2$$sums. 
    more » « less
  3. Abstract A result of Gyárfás [12] exactly determines the size of a largest monochromatic component in an arbitrary$$r$$-colouring of the complete$$k$$-uniform hypergraph$$K_n^k$$when$$k\geq 2$$and$$k\in \{r-1,r\}$$. We prove a result which says that if one replaces$$K_n^k$$in Gyárfás’ theorem by any ‘expansive’$$k$$-uniform hypergraph on$$n$$vertices (that is, a$$k$$-uniform hypergraph$$G$$on$$n$$vertices in which$$e(V_1, \ldots, V_k)\gt 0$$for all disjoint sets$$V_1, \ldots, V_k\subseteq V(G)$$with$$|V_i|\gt \alpha$$for all$$i\in [k]$$), then one gets a largest monochromatic component of essentially the same size (within a small error term depending on$$r$$and$$\alpha$$). As corollaries we recover a number of known results about large monochromatic components in random hypergraphs and random Steiner triple systems, often with drastically improved bounds on the error terms. Gyárfás’ result is equivalent to the dual problem of determining the smallest possible maximum degree of an arbitrary$$r$$-partite$$r$$-uniform hypergraph$$H$$with$$n$$edges in which every set of$$k$$edges has a common intersection. In this language, our result says that if one replaces the condition that every set of$$k$$edges has a common intersection with the condition that for every collection of$$k$$disjoint sets$$E_1, \ldots, E_k\subseteq E(H)$$with$$|E_i|\gt \alpha$$, there exists$$(e_1, \ldots, e_k)\in E_1\times \cdots \times E_k$$such that$$e_1\cap \cdots \cap e_k\neq \emptyset$$, then the smallest possible maximum degree of$$H$$is essentially the same (within a small error term depending on$$r$$and$$\alpha$$). We prove our results in this dual setting. 
    more » « less
  4. We explored the settling dynamics of vertically aligned particles in a quiescent, stratified two-layer fluid using particle tracking velocimetry. Glass spheres of$$d=4\,{\rm mm}$$diameter were released at frequencies of 4, 6 and 8 Hz near the free surface, traversing through an upper ethanol layer ($$H_1$$), whereHis height or layer thickess, varying from$$10d$$to$$40d$$and a lower oil layer. Results reveal pronounced lateral particle motion in the ethanol layer, attributed to a higher Galileo number ($$Ga = 976$$, ratio of buoyancy–gravity to viscous effects), compared with the less active behaviour in the oil layer ($$Ga = 16$$). The ensemble vertical velocity of particles exhibited a minimum just past the density interface, becoming more pronounced with increasing$$H_1$$, and suggesting that enhanced entrainment from ethanol to oil resulted in an additional buoyancy force. This produced distinct patterns of particle acceleration near the density interface, which were marked by significant deceleration, indicating substantial resistance to particle motion. An increased drag coefficient occurred for$$H_1/d = 40$$compared with a single particle settling in oil; drag reduced as the particle-release frequency ($$\,f_p$$) increased, likely due to enhanced particle interactions at closer proximity. Particle pair dispersions, lateral ($$R^2_L$$) and vertical ($$R^2_z$$), were modulated by$$H_1$$, initial separation$$r_0$$and$$f_p$$. The$$R^2_L$$dispersion displayed ballistic scaling initially, Taylor scaling for$$r_0 < H_1$$and Richardson scaling for$$r_0 > H_1$$. In contrast,$$R^2_z$$followed a$$R^2_z \sim t^{5.5}$$scaling under$$r_0 < H_1$$. Both$$R^2_L$$and$$R^2_z$$plateaued at a distance from the interface, depending on$$H_1$$and$$f_p$$. 
    more » « less
  5. 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