Normal stress and diffusion in a dilute suspension of hard spheres undergoing simple shear

Citation
Ie. Zarraga et Dt. Leighton, Normal stress and diffusion in a dilute suspension of hard spheres undergoing simple shear, PHYS FLUIDS, 13(3), 2001, pp. 565-577
Citations number
17
Categorie Soggetti
Physics
Journal title
PHYSICS OF FLUIDS
ISSN journal
10706631 → ACNP
Volume
13
Issue
3
Year of publication
2001
Pages
565 - 577
Database
ISI
SICI code
1070-6631(200103)13:3<565:NSADIA>2.0.ZU;2-R
Abstract
The complete set of normal stresses in a dilute suspension of hard spheres undergoing simple shear at low Reynolds number is calculated using a path i ntegration approach for the cases where the concentration is uniform and wh ere a small gradient in concentration is present. As expected, the normal s tresses are seen to be a strong function of epsilon (s)=2(b-a)/a, where b i s the hard sphere radius and a is the particle radius. The normal stress di fferences N-1 and N-2, are negative while the osmotic pressure is large and positive, with Pi > parallel toN(2)parallel to and N-1-->0 as epsilon (s)- -> infinity. For epsilon (s)much less than1, the asymmetry in the pair dist ribution function due to a depletion of particles in the extensional side o f a pair interaction leads to \N-1\> \N-2\. On the other hand, for epsilon (s)--> infinity, the additional stresslet induced when hard sphere radii to uch dominates the stress generated in the suspension, and N-2 becomes the p revailing normal stress difference. The self and gradient diffusivities are calculated using da Cunha and Hinch's [J. Fluid Mech. 309, 211 (1996)] tra jectory method. Numerical results for the diffusivities are in agreement wi th those obtained by da Cunha and Hinch for epsilon (s)less than or equal t o0.08 while matching the analytically obtained diffusivities for large epsi lon (s). Finally, we calculate the normal stress in the presence of a small concentration gradient and compare two models of migration for this case, namely the suspension balance model of Nott and Brady [J. Fluid Mech. 275, 157 (1994)] and the diffusive flux model first introduced by Leighton and A crivos [J. Fluid Mech. 181, 415 (1987)]. The results show that although the two models equally describe migration in the presence of a concentration g radient for the case where b much greater thana (or epsilon (s)--> infinity ), the two models are shown to be quantitatively different when near-field hydrodynamics are relevant. (C) 2001 American Institute of Physics.