Morris, Jeffrey F. and Brady, John F. (1996) Self-diffusion in sheared suspensions. Journal of Fluid Mechanics, 312 . pp. 223-252. ISSN 0022-1120. http://resolver.caltech.edu/CaltechAUTHORS:20120208-141206043
- Published Version
See Usage Policy.
Use this Persistent URL to link to this item: http://resolver.caltech.edu/CaltechAUTHORS:20120208-141206043
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 = ẏa^2/2D, where ẏ is the shear rate, a is the particle radius, and D = k_BT/6πηa is the diffusion coefficient of an isolated particle. Here, k_B 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 k_BT times the microstructural average of the hydrodynamic mobility of a particle, and depends on the volume fraction ø = 4/3π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 D_0^s and the correlation of the flux of a tagged particle with this perturbation. In a flowing suspension both D_0^s and D_∞^s are anisotropic, in general, with the anisotropy of D_0^s 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(φPe^(3/2)) for Pe « 1 and ø « 1, both with and without hydro-dynamic interactions between the particles. The leading dependence upon flow of D_0^s is 0.22DøPeÊ, where Ê is the rate-of-strain tensor made dimensionless with ẏ. Regardless of whether or not the particles interact hydrodynamically, flow influences D_∞^s at O(øPe) and O(øPe^(3/2)). In the absence of hydrodynamics, the leading correction is proportional to øPeDÊ. The correction of O(øPe^(3/2)), which results from a singular advection-diffusion problem, is proportional, in the absence of hydrodynamic interactions, to øPe^(3/2)DI; 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 Pe = ẏa^2/2D^s_0(φ)$. At small Pe the dependence on Pe is the same as at low ø.
|Additional Information:||© 1996 Cambridge University Press. Received 20 June 1995 and in revised form 20 October 1995. This work was supported in part by the National Science Foundation through grant CTS-9420415 and by the Office of Naval Research through grant N00014-95-1-0423. The suggestions of Dr Chris Harris of Shell Research, Amsterdam, for the evaluation of the integrals involved in the diffusivity calculation were very helpful.|
|Official Citation:||Self-diffusion in sheared suspensions preview Jeffrey F. Morris and John F. Brady Journal of Fluid Mechanics / Volume 312 / pp 223 - 252 Published online: 26 April 2006 DOI:10.1017/S002211209600198|
|Usage Policy:||No commercial reproduction, distribution, display or performance rights in this work are provided.|
|Deposited By:||Ruth Sustaita|
|Deposited On:||08 Feb 2012 23:14|
|Last Modified:||26 Dec 2012 14:49|
Repository Staff Only: item control page