The rheology of systems containing rigid spheres suspended in both viscous and viscoelastic media, studied by Stokesian dynamics simulations

An extensive Stokesian dynamics study of the rheological behavior of suspensions of rigid spheres subjected to an oscillating shear strain is presented. Two types of suspensions are considered: (1) rigid spheres in a viscous medium, and (2) rigid spheres in a viscoelastic medium. For this latter system we need to extend the Stokesian dynamics method, which was originally developed by Brady and Bossis for particles suspended in a viscous medium. The derivation of the necessary equations for these extended Stokesian dynamics simulations is given. In this derivation we use the well known correspondence principle and apply it to the set of equations for the hydrodynamic forces and stresslets that are valid for spheres suspended in a viscous medium. Then, using Fourier transformation we obtain differential equations for these forces and stresslets in the viscoelastic case. The contribution of the Brownian motion of the spheres to the bulk stress is found to be independent of the viscoelastic properties of the suspending medium. As an example of a viscoelastic medium we have chosen the Maxwell fluid. In our computer simulations we have calculated both elastic and viscous moduli and compared these results with experimental data. For the case of spheres suspended in a viscous medium we find that the elastic modulus reaches a plateau at high frequencies. Finally we present a simple analytical model which accurately reproduces the hydrodynamic contribution to the viscosity of spheres suspended in a Maxwell medium. This model is used to interpret the experimental results of Aral and Kaylon [Aral, B. K. and D. M. Kaylon, J. Rheol. 41, 599–620 (1977)] that were obtained at Peclet numbers currently inaccessible to Stokesian dynamics simulations.

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

[2]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

[3]  J. Mellema,et al.  Stress tensor expressions for dispersions , 1995 .

[4]  W. Gleissle,et al.  The characterization of the steady‐state shear and normal stress functions of highly concentrated suspensions formulated with viscoelastic liquids , 1993 .

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

[6]  C. Kruif,et al.  Linear viscoelasticity in dispersions of adhesive hard spheres , 1994 .

[7]  R. Phillips.  Dynamic simulation of hydrodynamically interacting spheres in a quiescent second-order fluid , 1996, Journal of Fluid Mechanics.

[8]  G. Nägele,et al.  On the dynamics and structure of charge-stabilized suspensions , 1996 .

[9]  C. Gardiner.  Handbook of Stochastic Methods , 1983 .

[10]  James J. Feng,et al.  Dynamic simulation of sedimentation of solid particles in an Oldroyd-B fluid , 1996 .

[11]  H. Hasegawa,et al.  On the usefulness of stochastic calculus applied to adiabatic elimination in Langevin equations , 1980 .

[12]  J. M. Sancho,et al.  Adiabatic elimination for systems of Brownian particles with nonconstant damping coefficients , 1982 .

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

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

[15]  D. Kalyon,et al.  Viscoelastic material functions of noncolloidal suspensions with spherical particles , 1997 .

[16]  Robert Lionberger,et al.  High frequency modulus of hard sphere colloids , 1994 .

[17]  J. Koelman,et al.  Dynamic simulations of hard-sphere suspensions under steady shear , 1993 .

[18]  D. Heyes,et al.  Linear viscoelasticity of concentrated hard-sphere dispersions , 1994 .

[19]  C. Kruif,et al.  Hard sphere colloidal dispersions: Mechanical relaxation pertaining to thermodynamic forces , 1987 .

[20]  R. Powell,et al.  Rheology of suspensions of spherical particles in a Newtonian and a non-Newtonian fluid , 1984 .

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

[22]  R. Verberg,et al.  Viscosity of colloidal suspensions , 1996, chao-dyn/9606008.

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

[24]  J. Brady.  The rheological behavior of concentrated colloidal dispersions , 1993 .

[25]  Evgenii A. Novikov,et al.  Functionals and the random-force method in turbulence theory , 1965 .

[26]  D. Bousfield,et al.  Prediction of suspension viscoelasticity through particle motion modeling , 1995 .

[27]  P. Senjanović,et al.  Toward a systematic description of heavy quark systems in QCD , 1978 .

[28]  A. Wagner,et al.  Simulation of a cusped bubble rising in a viscoelastic fluid with a new numerical method , 1999, cond-mat/9904029.

[29]  J. Brady,et al.  The rheology of Brownian suspensions , 1989 .

[30]  G. Vinogradov,et al.  Theory of dilute polymer solutions in viscoelastic fluid with a single relaxation time , 1984 .

[31]  G. Batchelor,et al.  The stress system in a suspension of force-free particles , 1970, Journal of Fluid Mechanics.

[32]  John F. Brady,et al.  Dynamic simulation of sheared suspensions. I. General method , 1984 .

[33]  Daniel D. Joseph,et al.  Direct simulation of the sedimentation of elliptic particles in Oldroyd-B fluids , 1998, Journal of Fluid Mechanics.

[34]  C. Broeck,et al.  Stochastic descriptions of the dynamics of interacting brownian particles , 1986 .

[35]  F. C. MacKintosh,et al.  Determining Microscopic Viscoelasticity in Flexible and Semiflexible Polymer Networks from Thermal Fluctuations , 1997 .

[36]  Dominique d'Humières,et al.  A lattice Boltzmann model for Jeffreys viscoelastic fluid , 1998 .

[37]  D. Heyes,et al.  Molecular simulations of colloidal liquids , 1994 .

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

[39]  G. Batchelor.  The effect of Brownian motion on the bulk stress in a suspension of spherical particles , 1977, Journal of Fluid Mechanics.

[40]  Peter V. Coveney,et al.  Simulating the rheology of dense colloidal suspensions using dissipative particle dynamics , 1997 .

[41]  R. Christensen.  Theory of viscoelasticity : an introduction , 1971 .

[42]  Howard Brenner,et al.  The Stokes resistance of an arbitrary particle , 1964 .