Brady, John F. and Phillips, Ronald J. and Lester, Julia C. and Bossis, Georges (1988) Dynamic simulation of hydrodynamically interacting suspensions. Journal of Fluid Mechanics, 195 . pp. 257-280. ISSN 0022-1120 http://resolver.caltech.edu/CaltechAUTHORS:20120531-135717062
- Published Version
See Usage Policy.
Use this Persistent URL to link to this item: http://resolver.caltech.edu/CaltechAUTHORS:20120531-135717062
A general method for computing the hydrodynamic interactions among an infinite suspension of particles, under the condition of vanishingly small particle Reynolds number, is presented. The method follows the procedure developed by O'Brien (1979) for constructing absolutely convergent expressions for particle interactions. For use in dynamic simulation, the convergence of these expressions is accelerated by application of the Ewald summation technique. The resulting hydrodynamic mobility and/or resistance matrices correctly include all far-field non-convergent interactions. Near-field lubrication interactions are incorporated into the resistance matrix using the technique developed by Durlofsky, Brady & Bossis (1987). The method is rigorous, accurate and computationally efficient, and forms the basis of the Stokesian-dynamics simulation method. The method is completely general and allows such diverse suspension problems as self-diffusion, sedimentation, rheology and flow in porous media to be treated within the same formulation for any microstructural arrangement of particles. The accuracy of the Stokesian-dynamics method is illustrated by comparing with the known exact results for spatially periodic suspensions.
|Additional Information:||© 1988 Cambridge University Press. Received 8 October 1987. Published Online April 21 2006. We would like to thank our colleague L. Durlofsky for many useful discussions concerning Stokesian dynamics. The work was supported in part by the Centre National pour la Recherche Scientifique and by the National Science Foundation under grants CBT-8696067 and INT-8413695. Computer time was provided by Centre de Calcul Vectoriel pour la Recherche and by the San Diego Supercomputer Center. R. J. P. would also like to acknowledge the support of an NSF fellowship.|
|Usage Policy:||No commercial reproduction, distribution, display or performance rights in this work are provided.|
|Deposited By:||Tony Diaz|
|Deposited On:||01 Jun 2012 23:35|
|Last Modified:||26 Dec 2012 15:17|
Repository Staff Only: item control page