Hostname: page-component-cd9895bd7-gvvz8 Total loading time: 0 Render date: 2024-12-20T07:16:28.510Z Has data issue: false hasContentIssue false

Dynamic simulation of suspensions of non-Brownian hard spheres

Published online by Cambridge University Press:  26 April 2006

D. I. Dratler
Affiliation:
Department of Chemical Engineering, University of Illinois, Urbana, IL 61801, USA Present address: Exxon Production Research Company, P.O. Box 2189, Houston, TX 77252-2189, USA.
W. R. Schowalter
Affiliation:
Department of Chemical Engineering, University of Illinois, Urbana, IL 61801, USA

Abstract

In this work, we investigate the suitability of models based solely on continuum hydrodynamics for Stokesian Dynamics simulations of sheared suspensions of non-Brownian hard spheres. The suspensions of interest consist of monolayers of uniform rigid spheres subjected to a linear shear field. Areal fractions ranged from ϕa = 0.2 to 0.6. For these suspensions, two sets of Stokesian Dynamics simulations were performed. For the first set, particle interactions were assumed to be strictly hydrodynamic in nature. These simulations are analogous to those of Brady & Bossis (1985) and Chang & Powell (1993). For the second set of simulations, particles were subjected to both hydrodynamic and short-range repulsive forces. The repulsion serves as a qualitative model of non-hydrodynamic effects important when particle separation distances are small. Results from both sets of simulations were found to be within the range of established experimental data for viscosities of suspensions. However, simulations employing the pure hydrodynamic model lead to very small separation distances between particles. These small separations give rise to particle overlaps that could not be eliminated by time-step refinement. The instantaneous number of overlaps increased with density and typically exceeded the number of particles at the highest densities considered. More critically, for very dense suspensions the simulations failed to approach a long-time asymptotic state. For simulations employing a short-range repulsive force, these problems were eliminated. The repulsion had the effect of preventing extremely small separations, thereby eliminating particle overlaps. At high concentrations, viscosities computed using the two methods are significantly different. This suggests that the dynamics of particles near contact have a significant impact on bulk properties. Furthermore, the results suggest that a critical aspect of the physics important at small particle separation distances is missing from the pure hydrodynamic model, making it unusable for computing microstructures of dense suspensions. In contrast, simulations employing a short-range repulsive force appear to produce more realistic microstructures, and can be performed even at very high densities.

Type
Research Article
Copyright
© 1996 Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

Allen, M. P. & Tildesley, D. J. 1987 Computer Simulation of Liquids. Clarendon Press.
Batchelor, G. K. 1970 The stress system in a suspension of force-free particles. J. Fluid Mech. 41, 545570.Google Scholar
Bonnecaze, R. T. & Brady, J. F. 1992 Dynamic simulation of an electrorheological fluid. J. Chem. Phys. 96, 21832202.Google Scholar
Bossis, G. & Brady, J. F. 1984 Dynamic simulation of sheared suspensions. I. General method. J. Chem. Phys. 80, 51415154.Google Scholar
Bossis, G. & Brady, J. F. 1989 The rheology of Brownian suspensions. J. Chem Phys. 91, 18661874.Google Scholar
Brady, J. F. & Bossis, G. 1985 The rheology of concentrated suspensions of spheres in simple shear flow by numerical simulation. J. Fluid Mech. 155, 105129.Google Scholar
Brady, J. F. & Bossis, G. 1988 Stokesian Dynamics. Ann. Rev. Fluid Mech. 20, 111157.Google Scholar
Brady, J. F., Phillips, R. J., Lester, J. C. & Bossis, G. 1988 Dynamic simulation of hydrodynamically interacting suspensions. J. Fluid Mech. 195, 257280.Google Scholar
Chang, C. & Powell, R. L. 1993 Dynamic simulation of bimodal suspensions of hydrodynamically interacting spherical particles. J. Fluid Mech. 253, 125.Google Scholar
Chang, C. & Powell, R. L. 1994 The rheology of bimodal hard-sphere dispersions. Phys. Fluids 6, 16281636.Google Scholar
Conte, S. D. & Boor, C. de 1980 Elementary Numerical Analysis: An Algorithmic Approach. McGraw-Hill.
Dratler, D. I., Schowalter, W. R. & Hoffman, R. L. 1996 Dynamic simulation of shear thickening in concentrated colloidal suspensions. J. Fluid Mech. (submitted).Google Scholar
Durlofsky, L. J. & Brady, J. F. 1989 Dynamic simulation of bounded suspensions of hydrodynamically interacting particles. J. Fluid Mech. 200, 3967.Google Scholar
Durlofsky, L., Brady, J. F. & Bossis, G. 1987 Dynamic simulation of hydrodynamically interacting particles. J. Fluid Mech. 180, 2149.Google Scholar
Heyes, D. M. & Melrose, J. R. 1993 Brownian dynamics simulations of model hard-sphere suspensions. J. Non-Newtonian Fluid Mech. 46, 128.Google Scholar
Lewis, T. B. & Nielsen, L. E. 1968 Viscosity of dispersed and aggregated suspensions of spheres. Trans. Soc. Rheol. 12, 421443.Google Scholar
Phung, T. & Brady, J. F. 1991 Microstructured fluids: structure, diffusion and rheology of colloidal dispersions. In Slow Dynamics of Condensed Matter (ed. K. Kawasaki, M. Tokuyama & T. Kawakatsu). AIP.
Russel, W. B., Saville, D. A. & Schowalter, W. R. 1989 Colloidal Dispersions. Cambridge University Press.
Sengen, M. Z. & Probstein, R. F. 1989 High-shear-limit viscosity and the maximum packing fraction in concentrated monomodal suspensions. PhysicoChem. Hydrodyn. 11, 229241.Google Scholar
Thomas, D. G. 1965 Transport characteristics of suspension: VIII. A note on the viscosity of Newtonian suspensions of uniform spherical particles. J. Colloid Sci. 20, 267277.Google Scholar