Force-induced diffusion in suspensions of hydrodynamically interacting colloids

2016 ◽  
Vol 795 ◽  
pp. 739-783 ◽  
Author(s):  
N. J. Hoh ◽  
R. N. Zia

We study the influence of hydrodynamic, thermodynamic and interparticle forces on the diffusive motion of a Brownian probe driven by a constant external force through a dilute colloidal dispersion. The influence of these microscopic forces on equilibrium self-diffusivity (passive microrheology) is well known: all three act to hinder the short- and long-time self-diffusion. Here, via pair-Smoluchowski theory, we explore their influence on self-diffusion in a flowing suspension, where particles and fluid have been set into motion by an externally forced probe (active microrheology), giving rise to non-equilibrium flow-induced diffusion. The probe’s motion entrains background particles as it travels through the bath, deforming the equilibrium suspension microstructure. The shape and extent of microstructural distortion is set by the relative strength of the external force $F^{\mathit{ext}}$ to the entropic restoring force $kT/a_{\mathit{th}}$ of the bath particles, defining a Péclet number $\mathit{Pe}\equiv F^{\mathit{ext}}/(2kT/a_{\mathit{th}})$; and also by the strength of hydrodynamic interactions, set by the range of interparticle repulsion ${\it\kappa}=(a_{\mathit{th}}-a)/a$, where $kT$ is the thermal energy and $a_{\mathit{th}}$ and $a$ are the thermodynamic and hydrodynamic sizes of the particles, respectively. We find that in the presence of flow, the same forces that hinder equilibrium diffusion now enhance it, with diffusive anisotropy set by the range of interparticle repulsion ${\it\kappa}$. A transition from hindered to enhanced diffusion occurs when diffusive and advective forces balance, $\mathit{Pe}\sim 1$, where the exact value is a sensitive function of the strength of hydrodynamics, ${\it\kappa}$. We find that the hindered to enhanced transition straddles two transport regimes: in hindered diffusion, stochastic forces in the presence of other bath particles produce deterministic displacements (Brownian drift) at the expense of a maximal random walk. In enhanced diffusion, driving the probe with a deterministic force through an initially random suspension leads to fluctuations in the duration of probe–bath particle entrainment, giving rise to enhanced, flow-induced diffusion. The force-induced diffusion is anisotropic for all $\mathit{Pe}$, scaling as $D\sim \mathit{Pe}^{2}$ in all directions for weak forcing, regardless of the strength of hydrodynamic interactions. When probe forcing is strong, $D\sim \mathit{Pe}$ in all directions in the absence of hydrodynamic interactions, but the picture changes qualitatively as hydrodynamic interactions grow strong. In this nonlinear regime, microstructural asymmetry weakens while the anisotropy of the force-induced diffusion tensor increases dramatically. This behaviour owes its origins to the approach toward Stokes flow reversibility, where diffusion along the direction of probe force scales advectively while transverse diffusion must vanish.

1996 ◽  
Vol 312 ◽  
pp. 223-252 ◽  
Author(s):  
Jeffrey F. Morris ◽  
John F. Brady

Self-diffusion in a suspension of spherical particles in steady linear shear flow is investigated by following the time evolution of the correlation of number density fluctuations. Expressions are presented for the evaluation of the self-diffusivity in a suspension which is either raacroscopically quiescent or in linear flow at arbitrary Peclet number $Pe = \dot{\gamma}a^2/2D$, where $\dot{\gamma}$ is the shear rate, a is the particle radius, and D = kBT/6πa is the diffusion coefficient of an isolated particle. Here, kB is Boltzmann's constant, T is the absolute temperature, and η is the viscosity of the suspending fluid. The short-time self-diffusion tensor is given by kBT times the microstructural average of the hydrodynamic mobility of a particle, and depends on the volume fraction $\phi = \frac{4}{3}\pi a^3n$ and Pe only when hydrodynamic interactions are considered. As a tagged particle moves through the suspension, it perturbs the average microstructure, and the long-time self-diffusion tensor, D∞s, is given by the sum of D0s and the correlation of the flux of a tagged particle with this perturbation. In a flowing suspension both D0s and D∞ are anisotropic, in general, with the anisotropy of D0s due solely to that of the steady microstructure. The influence of flow upon D∞s is more involved, having three parts: the first is due to the non-equilibrium microstructure, the second is due to the perturbation to the microstructure caused by the motion of a tagged particle, and the third is by providing a mechanism for diffusion that is absent in a quiescent suspension through correlation of hydrodynamic velocity fluctuations.The self-diffusivity in a simply sheared suspension of identical hard spheres is determined to O(øPe3/2) for Pe ≤ 1 and ø ≤ 1, both with and without hydro-dynamic interactions between the particles. The leading dependence upon flow of D0s is 0.22DøPeÊ, where Ê is the rate-of-strain tensor made dimensionless with $\dot{\gamma}$. Regardless of whether or not the particles interact hydrodynamically, flow influences D∞s at O(øPe) and O(øPe3/2). In the absence of hydrodynamics, the leading correction is proportional to øPeDÊ. The correction of O(øPe3/2), which results from a singular advection-diffusion problem, is proportional, in the absence of hydrodynamic interactions, to øPe3/2DI; when hydrodynamics are included, the correction is given by two terms, one proportional to Ê, and the second a non-isotropic tensor.At high ø a scaling theory based on the approach of Brady (1994) is used to approximate D∞s. For weak flows the long-time self-diffusivity factors into the product of the long-time self-diffusivity in the absence of flow and a non-dimensional function of $\bar{P}e = \dot{\gamma}a^2/2D^s_0(\phi)$. At small $\bar{P}e$ the dependence on $\bar{P}e$ is the same as at low ø.


1999 ◽  
Vol 401 ◽  
pp. 243-274 ◽  
Author(s):  
DAVID R. FOSS ◽  
JOHN F. BRADY

The behaviour of the long-time self-diffusion tensor in concentrated colloidal dispersions is studied using dynamic simulation. The simulations are of a suspension of monodisperse Brownian hard spheres in simple shear flow as a function of the Péclet number, Pe, which measures the relative importance of shear and Brownian forces, and the volume fraction, ϕ. Here, Pe = &γdot;a2/D0, where &γdot; is the shear rate, a the particle size and D0 = kT/6πηa is the Stokes–Einstein diffusivity of an isolated particle of size a with thermal energy kT in a solvent of viscosity η. Two simulations algorithms are used: Stokesian Dynamics for inclusion of the many-body hydrodynamic interactions, and Brownian Dynamics for suspensions without hydrodynamic interactions. A new procedure for obtaining high-quality diffusion data based on averaging the results of many short simulations is presented and utilized. At low shear rates, low Pe, Brownian diffusion due to a random walk process dominates and the characteristic scale for diffusion is the Stokes–Einstein diffusivity, D0. At zero Pe the diffusivity is found to be a decreasing function of ϕ. As Pe is slowly increased, O(Pe) and O(Pe3/2) corrections to the diffusivity due to the flow are clearly seen in the Brownian Dynamics system in agreement with the theoretical results of Morris & Brady (1996). At large shear rates, large Pe, both systems exhibit diffusivities that grow linearly with the shear rate by the non-Brownian mechanism of shear-induced diffusion. In contrast to the behaviour at low Pe, this shear-induced diffusion mode is an increasing function of ϕ. Long-time rotational self-diffusivities are of interest in the Stokesian Dynamics system and show similar behaviour to their translational analogues. An off-diagonal long-time self-diffusivity, Dxy, is reported for both systems. Results for both the translational and rotational Dxy show a sign change from low Pe to high Pe due to different mechanisms in the two regimes. A physical explanation for the off-diagonal diffusivities is proposed.


1994 ◽  
Vol 281 ◽  
pp. 51-80 ◽  
Author(s):  
Chingyi Chang ◽  
Robert L. Powell

We study the average mobilities and long-time self-diffusion coefficients of a suspension of bimodally distributed spherical particles. Stokesian dynamics is used to calculate the particle trajectories for a monolayer of bimodal-sized spheres. Hydrodynamic forces only are considered and they are calculated using the inverse of the grand mobility matrix for far-field many-body interactions and lubrication formulae for near-field effects. We determine both the detailed microstructure (e.g. the pair-connectedness function and cluster formation) and the macroscopic properties (e.g. viscosity and self-diffusion coefficients). The flow of an ‘infinite’ suspension is simulated by considering 25, 49, 64 and 100 particles to be one ‘cell’ of a periodic array. Effects of both the size ratio and the relative fractions of the different-sized particles are examined. For the microstructures, the pair-connectedness function shows that the particles form clusters in simple shearing flow due to lubrication forces. The nearly symmetric angular structures imply the absence of normal stress differences for a suspension with purely hydrodynamic interactions between spheres. For average mobilities at infinite Péclet number, Ds0, our simulation results suggest that the reduction of Ds0 as concentration increases is directly linked to the influence of particle size distribution on the average cluster size. For long-time self-diffusion coefficients, Ds∞, we found good agreement between simulation and experiment (Leighton & Acrovos 1987 a; Phan and Leighton 1993) for monodispersed suspensions. For bimodal suspensions, the magnitude of Ds∞, and the time to reach the asymptotic diffusive behaviour depend on the cluster size formed in the system, or the viscosity of the suspension. We also consider the effect of the initial configuration by letting the spheres be both organized (size segregated) and randomly placed. We find that it takes a longer time for a suspension with an initially organized structure to achieve steady state than one with a random structure.


2014 ◽  
Vol 140 (24) ◽  
pp. 244116 ◽  
Author(s):  
Claudio Contreras-Aburto ◽  
César A. Báez ◽  
José M. Méndez-Alcaraz ◽  
Ramón Castañeda-Priego

2002 ◽  
Vol 719 ◽  
Author(s):  
Ian D. Sharp ◽  
Hartmut A. Bracht ◽  
Hughes H. Silvestri ◽  
Samuel P. Nicols ◽  
Jeffrey W. Beeman ◽  
...  

AbstractIsotopically controlled silicon multilayer structures were used to measure the enhancement of self- and dopant diffusion in extrinsic boron doped silicon. 30Si was used as a tracer through a multilayer structure of alternating natural Si and enriched 28Si layers. Low energy, high resolution secondary ion mass spectrometry (SIMS) allowed for simultaneous measurement of self- and dopant diffusion profiles of samples annealed at temperatures between 850°C and 1100°C. A specially designed ion-implanted amorphous Si surface layer was used as a dopant source to suppress excess defects in the multilayer structure, thereby eliminating transient enhanced diffusion (TED) behavior. Self- and dopant diffusion coefficients, diffusion mechanisms, and native defect charge states were determined from computer-aided modeling, based on differential equations describing the diffusion processes. We present a quantitative description of B diffusion enhanced self-diffusion in silicon and conclude that the diffusion of both B and Si is mainly mediated by neutral and singly positively charged self-interstitials under p-type doping. No significant contribution of vacancies to either B or Si diffusion is observed.


2018 ◽  
Author(s):  
Kazutaka Takao ◽  
Kazunori Takamiya ◽  
Da-Qiao Ding ◽  
Tokuko Haraguchi ◽  
Yasushi Hiraoka ◽  
...  

AbstractHomologous sets of parental chromosomes must pair during meiosis to produce recombined sets of chromosomes for their progeny. This is accompanied by nuclear oscillatory movements. This study aimed to elucidate the significance of these movements with a model, wherein external force was applied to the oscillating nucleus and via hydrodynamic interactions within the nucleus. Simulations revealed that a major force for aligning homologous chromosomes is length-dependent sorting during chromosomal torsional turning, which occur when the nucleus reverses the direction of its movement.


2021 ◽  
Vol 925 ◽  
Author(s):  
Emma Gonzalez ◽  
Christian Aponte-Rivera ◽  
Roseanna N. Zia

We present a computational study of the equilibrium dynamics of a polydisperse hard-sphere colloidal dispersion confined in a spherical cavity. We account for many-body hydrodynamic and lubrication interactions between particles and with the confining cavity utilizing our confined Stokesian dynamics model, expanded here for size polydispersity. We find that, even though the tendency of polydispersity to homogenize structure in a suspension is still present in confinement, strong correlations induced by the cavity resist homogenization. Although seemingly opposite, these two effects have a common driver, which is to maximize configurational entropy of particles in the cavity interior. These structural effects couple with the hydrodynamics to change the particle dynamics: polydispersity weakens lubrication effects near the cavity wall, allowing small (large) particles to diffuse faster (slower) than in a monodisperse suspension. As a small (large) particle gets farther from the wall, polydispersity weakens many-body hydrodynamic couplings, driving diffusivity up (down). While the local cage dynamics dominates short-time self-diffusion, long-time dynamics is also affected. In the concentrated regime, polydispersity and confinement combine to induce radial de-mixing into size-segregated populations. The cavity becomes the most influential ‘nearest neighbour’, setting the length scale of and dynamics within these radial domains. This intermediate length-scale caging makes the angular dynamics insensitive to polydispersity but leads to radial long-time mean-square displacement that changes qualitatively with volume composition. These results hold promise for explaining colloidal-scale physics implicated in the functioning of biological cells, and the engineering of non-living confined colloids where size de-mixing could be useful in the design of encapsulated micro-reactors and therapeutic vesicles.


2020 ◽  
Vol 101 (4) ◽  
Author(s):  
Nima H. Siboni ◽  
Alice L. Thorneywork ◽  
Alicia Damm ◽  
Roel P. A. Dullens ◽  
Jürgen Horbach

Sign in / Sign up

Export Citation Format

Share Document