An official website of the United States government
Discrete-particle simulations of bidisperse shear thickening suspensions are reported. The work considers two packing parameters, the large-to-small particle radius ratio ranging from [Formula: see text] (nearly monodisperse) to [Formula: see text], and the large particle fraction of the total solid loading with values [Formula: see text], 0.5, and 0.85. Particle-scale simulations are performed over a broad range of shear stresses using a simulation model for spherical particles accounting for short-range lubrication forces, frictional interaction, and repulsion between particles. The variation of rheological properties and the maximum packing fraction [Formula: see text] with shear stress [Formula: see text] are reported. At a fixed volume fraction [Formula: see text], bidispersity decreases the suspension relative viscosity [Formula: see text], where [Formula: see text] is the suspension viscosity and [Formula: see text] is the suspending fluid viscosity, over the entire range of shear stresses studied. However, under low shear stress conditions, the suspension exhibits an unusual rheological behavior: the minimum viscosity does not occur as expected at [Formula: see text], but instead decreases with further increase of [Formula: see text] to [Formula: see text]. The second normal stress difference [Formula: see text] acts similarly. This behavior is caused by particles ordering into a layered structure, as is also reflected by the zero slope with respect to time of the mean-square displacement in the velocity gradient direction. The relative viscosity [Formula: see text] of bidisperse rate-dependent suspensions can be predicted by a power law linking it to [Formula: see text], [Formula: see text] in both low and high shear stress regimes. The agreement between the power law and experimental data from literature demonstrates that the model captures well the effect of particle size distribution, showing that viscosity roughly collapses onto a single master curve when plotted against the reduced volume fraction [Formula: see text].
more »
« less
Ordered domains in sheared dense suspensions: The link to viscosity and the disruptive effect of friction
Monodisperse suspensions of Brownian colloidal spheres crystallize at high densities, and ordering under shear has been observed at densities below the crystallization threshold. We perform large-scale simulations of a model suspension containing over [Formula: see text] particles to quantitatively study the ordering under shear and to investigate its link to the rheological properties of the suspension. We find that at high rates, for [Formula: see text], the shear flow induces an ordering transition that significantly decreases the measured viscosity. This ordering is analyzed in terms of the development of layering and planar order, and we determine that particles are packed into hexagonal crystal layers (with numerous defects) that slide past each other. By computing local [Formula: see text] and [Formula: see text] order parameters, we determine that the defects correspond to chains of particles in a squarelike lattice. We compute the individual particle contributions to the stress tensor and discover that the largest contributors to the shear stress are primarily located in these lower density, defect regions. The defect structure enables the formation of compressed chains of particles to resist the shear, but these chains are transient and short-lived. The inclusion of a contact friction force allows the stress-bearing structures to grow into a system-spanning network, thereby disrupting the order and drastically increasing the suspension viscosity.
more »
« less
- Award ID(s):
- 2026842
- PAR ID:
- 10376796
- Date Published:
- Journal Name:
- Journal of Rheology
- Volume:
- 66
- Issue:
- 5
- ISSN:
- 0148-6055
- Page Range / eLocation ID:
- 1055 to 1065
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We introduce a novel approach to reveal ordering fluctuations in sheared dense suspensions, using line scanning in a combined rheometer and laser scanning confocal microscope. We validate the technique with a moderately dense suspension, observing modest shear-induced ordering and a nearly linear flow profile. At high concentration ( ϕ = 0.55) and applied stress just below shear thickening, we report ordering fluctuations with high temporal resolution, and directly measure a decrease in order with distance from the suspension’s bottom boundary as well as a direct correlation between order and particle concentration. Higher applied stress produces shear thickening with large fluctuations in boundary stress which we find are accompanied by dramatic fluctuations in suspension flow speeds. The peak flow rates are independent of distance from the suspension boundary, indicating that they likely arise from transient jamming that creates solid-like aggregates of particles moving together, but only briefly because the high speed fluctuations are interspersed with regions flowing much more slowly, suggesting that shear thickening suspensions possess complex internal structural dynamics, even in relatively simple geometries.more » « less
-
We investigate the rheological behavior of athermal particle suspensions using experiments and theory. A generalized version of the homogenization estimates of Ponte Castañeda and Willis [J. Mech. Phys. Solids, 43(12), 1919–1951 (1995)] is presented for the effective viscosity of athermal suspensions accounting for additional microstructural features (e.g., polydispersity) via an empirical parameter, [Formula: see text]. For the case of identically sized spheres dispersed with statistical isotropy in a Newtonian fluid, the parameter [Formula: see text] is estimated from the results of Batchelor and Green [J. Fluid Mech. 56(2), 375–400 (1972)] for the Huggins coefficient. Predictions for the macroscopic viscosity are found to be in good agreement with measurements for monodisperse polymethyl methacrylate (PMMA) spheres in glycerol, as well as for the empirical Krieger–Dougherty equation for the shear viscosity. The proposed estimates have the added benefit that they can also be used to get information on the statistics of the stress and strain-rate fields in the fluid and particle phases. In addition, results for the effective shear viscosity are used in combination with the linear comparison method of Ponte Castañeda [J. Mech. Phys. Solids 39(1), 45–71 (1991)] to generate the corresponding estimates for the effective macroscopic behavior and field statistics of particle suspensions in (viscoplastic) yield stress fluids. Good agreement is also found between the theoretical estimates and experimental results for the effective yield and flow stress of suspensions with monodisperse PMMA spheres in Carbopol. Finally, it is argued that the results for the phase averages and fluctuations of the stress and strain-rate fields can be used to provide a physical interpretation for the parameter [Formula: see text] in terms of the polydispersity of the suspension and its implications for the percolation threshold.more » « less
-
Dense particulate suspensions demonstrate a significant increase in average viscosity beyond a material-specific critical shear stress. Here, we analyze the steady-state structure of a suspension of monodisperse silica microspheres in the shear-thickening regime. Using dynamic measurement of boundary stress, we show that the flow is characterized by a cluster of high-stress fronts that propagate in the flow direction at a speed of 1/2 relative to the top plate of the rheometer. We apply high-speed line scan imaging to reveal dramatic fluctuations in particle speed, ordering, and concentration associated with the fronts and show that the structure is consistent with transiently jammed networks that contain high interparticle stresses that percolate across the rheometer gap, but which are present only briefly during the passage of the high-stress fronts.more » « less
-
Hypothesis: The viscosity of dense suspensions surges when the applied stress surpasses a material-specific critical threshold. There is growing evidence that the thickening transition involves non-uniform flow and stress with considerable spatiotemporal complexity. Nevertheless, it is anticipated that dense suspensions of calcium carbonate particles with purely repulsive interactions may not conform to this scenario, as indicated by local pressure measurements with millimeter spatial resolution. Experiment: Here we utilize Boundary Stress Microscopy (BSM), a technique capable of resolving stresses down to the micron scale, to search for evidence of stress heterogeneity. In addition, we measure the flow field at the lower boundary of the suspension where the boundary stress is measured. Findings: We find localized regions of high-stresses that are extended in the vorticity direction and propagate in the flow direction at a speed approximately half that of the rheometer’s top plate. These high-stress regions proliferate with the applied stress accounting for the increased viscosity. Furthermore, the velocity of particles at the lower boundary of the suspension shows a significant and complex nonaffine flow that accompanies regions of high-stresses. Hence, our findings demonstrate that stress and flow inhomogeneity are intrinsic characteristics of shear-thickening suspensions, regardless of the nature of interparticle interactions.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1122/8.0000453
