Accelerated Stokesian Dynamics simulations

A new implementation of the conventional Stokesian Dynamics (SD) algorithm, called accelerated Stokesian Dynamics (ASD), is presented. The equations governing the motion of N particles suspended in a viscous fluid at low particle Reynolds number are solved accurately and efficiently, including all hydrodynamic interactions, but with a significantly lower computational cost of O(N ln N). The main differences from the conventional SD method lie in the calculation of the many-body long-range interactions, where the Ewald-summed wave-space contribution is calculated as a Fourier transform sum and in the iterative inversion of the now sparse resistance matrix. The new method is applied to problems in the rheology of both structured and random suspensions, and accurate results are obtained with much larger numbers of particles. With access to larger N, the high-frequency dynamic viscosities and short-time self-diffusivities of random suspensions for volume fractions above the freezing point are now studied. The ASD method opens up an entire new class of suspension problems that can be investigated, including particles of non-spherical shape and a distribution of sizes, and the method can readily be extended to other low-Reynolds-number-flow problems.

[1]  A. Ladd.  Numerical simulations of particulate suspensions via a discretized Boltzmann equation. Part 2. Numerical results , 1993, Journal of Fluid Mechanics.

[2]  H. Lekkerkerker,et al.  Measurement of the short‐time self‐diffusion coefficient in dilute and concentrated suspensions: Influence of direct particle interactions , 1987 .

[3]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[4]  R. Ottewill,et al.  Study of particle motion in concentrated dispersions by tracer diffusion , 1987, Nature.

[5]  David J. Jeffrey,et al.  Calculation of the resistance and mobility functions for two unequal rigid spheres in low-Reynolds-number flow , 1984, Journal of Fluid Mechanics.

[6]  Wiley,et al.  Numerical simulation of the dense random packing of a binary mixture of hard spheres: Amorphous metals. , 1987, Physical review. B, Condensed matter.

[7]  Thomas E. Cheatham,et al.  Particle-Mesh Based Methods for Fast Ewald Summation in Molecular Dynamics Simulations , 1997, PPSC.

[8]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[9]  J. Brady,et al.  Structure, diffusion and rheology of Brownian suspensions by Stokesian Dynamics simulation , 2000, Journal of Fluid Mechanics.

[10]  Thomas M Truskett,et al.  Is random close packing of spheres well defined? , 2000, Physical review letters.

[11]  J. Brady,et al.  Stokesian Dynamics simulation of Brownian suspensions , 1996, Journal of Fluid Mechanics.

[12]  D. I. Dratler,et al.  Dynamic simulation of suspensions of non-Brownian hard spheres , 1996, Journal of Fluid Mechanics.

[13]  Erik Karl Guckel.  Large scale simulations of particulate systems using the PME method , 1999 .

[14]  D. S. Pearson,et al.  Viscoelastic behavior of concentrated spherical suspensions , 1994 .

[15]  J. Pasciak,et al.  Computer solution of large sparse positive definite systems , 1982 .

[16]  E. Cuthill,et al.  Reducing the bandwidth of sparse symmetric matrices , 1969, ACM '69.

[17]  H. Hasimoto.  On the periodic fundamental solutions of the Stokes equations and their application to viscous flow past a cubic array of spheres , 1959, Journal of Fluid Mechanics.

[18]  Segrè,et al.  Short-time Brownian motion in colloidal suspensions: Experiment and simulation. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[19]  George M. Homsy,et al.  Stokes flow through periodic arrays of spheres , 1982, Journal of Fluid Mechanics.

[20]  Rintoul,et al.  Metastability and Crystallization in Hard-Sphere Systems. , 1996, Physical review letters.

[21]  J. Brady,et al.  Dynamic simulation of hydrodynamically interacting suspensions , 1988, Journal of Fluid Mechanics.

[22]  A. Sangani,et al.  An O(N) algorithm for Stokes and Laplace interactions of particles , 1996 .

[23]  Louis J. Durlofsky,et al.  Dynamic simulation of hydrodynamically interacting particles , 1987, Journal of Fluid Mechanics.

[24]  de Kruif CG,et al.  Linear viscoelastic behavior of dense hard-sphere dispersions. , 1989, Physical review. A, General physics.

[25]  R. J. Speedy.  Pressure of the metastable hard-sphere fluid , 1997 .

[26]  John F. Brady,et al.  Hydrodynamic transport properties of hard-sphere dispersions. II. Porous media , 1988 .

[27]  Salvatore Torquato,et al.  Computer simulations of dense hard‐sphere systems , 1996 .

[28]  G. Batchelor,et al.  The determination of the bulk stress in a suspension of spherical particles to order c2 , 1972, Journal of Fluid Mechanics.

[29]  W. V. Megen,et al.  Tracer diffusion in concentrated colloidal dispersions. III. Mean squared displacements and self‐diffusion coefficients , 1989 .

[30]  John F. Brady,et al.  Many-body effects and matrix inversion in low-Reynolds-number hydrodynamics , 2001 .

[31]  J. Hofman.  Hydrodynamic properties of dense colloidal crystals , 1999 .

[32]  P. Pusey,et al.  Measurement of the wave‐vector dependent diffusion coefficient in concentrated particle dispersions , 1985 .

[33]  A. Ladd.  Numerical simulations of particulate suspensions via a discretized Boltzmann equation. Part 1. Theoretical foundation , 1993, Journal of Fluid Mechanics.

[34]  A. Sangani,et al.  A method for computing Stokes flow interactions among spherical objects and its application to suspensions of drops and porous particles , 1994 .

[35]  H. Brenner,et al.  Spatially periodic suspensions of convex particles in linear shear flows. III. Dilute arrays of spheres suspended in Newtonian fluids , 1983 .

[36]  C. Beenakker.  Ewald sum of the Rotne-Prager tensor , 1986 .

[37]  P. Saffman,et al.  On the Settling Speed of Free and Fixed Suspensions , 1973 .

[38]  J. Brady,et al.  Hydrodynamic transport properties of hard-sphere dispersions. I. Suspensions of freely mobile particles , 1988 .

[39]  H. G. Petersen.  Accuracy and efficiency of the particle mesh Ewald method , 1995 .

[40]  Anthony J. C. Ladd,et al.  Hydrodynamic transport coefficients of random dispersions of hard spheres , 1990 .

[41]  R. J. Speedy,et al.  On the reproducibility of glasses , 1994 .

[42]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[43]  John F. Brady,et al.  Accelerated Stokesian dynamics: Brownian motion , 2003 .

[44]  David J. Jeffrey,et al.  The pressure moments for two rigid spheres in low-Reynolds-number flow , 1993 .

[45]  P. Pusey,et al.  Measurement of the short-time self-mobility of particles in concentrated suspension. Evidence for many-particle hydrodynamic interactions , 1983 .

[46]  Joseph B. Keller,et al.  Effective viscosity of a periodic suspension , 1984, Journal of Fluid Mechanics.

[47]  I. Snook,et al.  Tracer diffusion in concentrated colloidal dispersions , 1986 .