2 results
Dynamic simulation of suspensions of non-Brownian hard spheres
- D. I. Dratler, W. R. Schowalter
-
- Journal:
- Journal of Fluid Mechanics / Volume 325 / 25 October 1996
- Published online by Cambridge University Press:
- 26 April 2006, pp. 53-77
-
- Article
- Export citation
-
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.
Dynamic simulation of shear thickening in concentrated colloidal suspensions
- D. I. DRATLER, W. R. SCHOWALTER, R. L. HOFFMAN
-
- Journal:
- Journal of Fluid Mechanics / Volume 353 / 25 December 1997
- Published online by Cambridge University Press:
- 25 December 1997, pp. 1-30
-
- Article
- Export citation
-
Stokesian Dynamics has been used to investigate the origins of shear thickening in concentrated colloidal suspensions. For this study, we considered a monolayer suspension composed of charge-stabilized non-Brownian monosized rigid spheres dispersed at an areal fraction of ϕa=0.74 in a Newtonian liquid. The suspension was subjected to a linear shear field. In agreement with established experimental data, our results indicate that shear thickening in this system is associated with an order–disorder transition of the suspension microstructure. Below the critical shear rate at which this transition occurs, the suspension microstructure consists of two-dimensional analogues of experimentally observed sliding layer configurations. Above this critical shear rate, suspensions are disordered, contain particle clusters, and exhibit viscosities and microstructures characteristic of suspensions of non-Brownian hard spheres. In addition, suspensions possessing the sliding layer microstructure at the beginning of supercritical shearing tend to retain this microstructure for a period of time before disordering. The onset of this disorder is due to the formation of particle doublets within the suspension. Once formed, these doublets rotate, due to the bulk motion, and disrupt the long-range order of the suspension. The cross-stream component of the centre-to-centre separation vector associated with the two particles forming a doublet, which is zero when the doublet is perfectly aligned with the bulk velocity vector, grows exponentially with time. This strongly suggests that the evolution of these doublets is due to a change in the stability of the sliding layer configurations, with this type of ordered microstructure being linearly unstable above a critical shear rate. This contention is supported by results of a stability analysis. The analysis shows that a single string of particles is subject to a linear instability leading to the formation of particle doublets. Simulations were repeated with different numbers of particles in the computational domain, with the results found to be qualitatively independent of system size.