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
Dynamics and memory of boundary stresses in discontinuous shear thickening suspensions during oscillatory shear
We report direct measurements of spatially resolved surface stresses of a dense suspension during large amplitude oscillatory shear (LAOS) in the discontinuous shear thickening regime using boundary stress microscopy. Consistent with previous studies, bulk rheology shows a dramatic increase in the complex viscosity above a frequency-dependent critical strain. We find that the viscosity increase is coincident with that appearance of large heterogeneous boundary stresses, indicative of the formation of transient solid-like phases (SLPs) on spatial scales large compared to the particle size. The critical strain for the appearance of SLPs is largely determined by the peak oscillatory stress, which depends on the peak shear rate and the frequency-dependent suspension viscosity. The SLPs dissipate and reform on each cycle, with a spatial pattern that is highly variable at low frequencies but remarkably persistent at the highest frequency measured ( ω = 10 rad s −1 ).
more »
« less
- Award ID(s):
- 1809890
- NSF-PAR ID:
- 10281696
- Date Published:
- Journal Name:
- Soft Matter
- Volume:
- 17
- Issue:
- 5
- ISSN:
- 1744-683X
- Page Range / eLocation ID:
- 1337 to 1345
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Understanding the orientation dynamics of anisotropic colloidal particles is important for suspension rheology and particle self-assembly. However, even for the simplest case of dilute suspensions in shear flow, the orientation dynamics of non-spherical Brownian particles are poorly understood. Here we analytically calculate the time-dependent orientation distributions for non-spherical axisymmetric particles confined to rotate in the flow–gradient plane, in the limit of small but non-zero Brownian diffusivity. For continuous shear, despite the complicated dynamics arising from the particle rotations, we find a coordinate change that maps the orientation dynamics to a diffusion equation with a remarkably simple ratio of the enhanced rotary diffusivity to the zero shear diffusion: $D_{eff}^{r}/D_{0}^{r}=(3/8)(p-1/p)^{2}+1$ , where $p$ is the particle aspect ratio. For oscillatory shear, the enhanced diffusion becomes orientation dependent and drastically alters the long-time orientation distributions. We describe a general method for solving the time-dependent oscillatory shear distributions and finding the effective diffusion constant. As an illustration, we use this method to solve for the diffusion and distributions in the case of triangle-wave oscillatory shear and find that they depend strongly on the strain amplitude and particle aspect ratio. These results provide new insight into the time-dependent rheology of suspensions of anisotropic particles. For continuous shear, we find two distinct diffusive time scales in the rheology that scale separately with aspect ratio $p$ , as $1/D_{0}^{r}p^{4}$ and as $1/D_{0}^{r}p^{2}$ for $p\gg 1$ . For oscillatory shear flows, the intrinsic viscosity oscillates with the strain amplitude. Finally, we show the relevance of our results to real suspensions in which particles can rotate freely. Collectively, the interplay between shear-induced rotations and diffusion has rich structure and strong effects: for a particle with aspect ratio 10, the oscillatory shear intrinsic viscosity varies by a factor of ${\approx}2$ and the rotational diffusion by a factor of ${\approx}40$ .more » « less
-
We present a numerical study of non-colloidal spherical and rigid particles suspended in Newtonian, shear thinning and shear thickening fluids employing an immersed boundary method. We consider a linear Couette configuration to explore a wide range of solid volume fractions ( $0.1\leqslant \unicode[STIX]{x1D6F7}\leqslant 0.4$ ) and particle Reynolds numbers ( $0.1\leqslant Re_{p}\leqslant 10$ ). We report the distribution of solid and fluid phase velocity and solid volume fraction and show that close to the boundaries inertial effects result in a significant slip velocity between the solid and fluid phase. The local solid volume fraction profiles indicate particle layering close to the walls, which increases with the nominal $\unicode[STIX]{x1D6F7}$ . This feature is associated with the confinement effects. We calculate the probability density function of local strain rates and compare the latter’s mean value with the values estimated from the homogenisation theory of Chateau et al. ( J. Rheol. , vol. 52, 2008, pp. 489–506), indicating a reasonable agreement in the Stokesian regime. Both the mean value and standard deviation of the local strain rates increase primarily with the solid volume fraction and secondarily with the $Re_{p}$ . The wide spectrum of the local shear rate and its dependency on $\unicode[STIX]{x1D6F7}$ and $Re_{p}$ point to the deficiencies of the mean value of the local shear rates in estimating the rheology of these non-colloidal complex suspensions. Finally, we show that in the presence of inertia, the effective viscosity of these non-colloidal suspensions deviates from that of Stokesian suspensions. We discuss how inertia affects the microstructure and provide a scaling argument to give a closure for the suspension shear stress for both Newtonian and power-law suspending fluids. The stress closure is valid for moderate particle Reynolds numbers, $O(Re_{p})\sim 10$ .more » « less
-
We report direct measurements of spatially resolved stress at the boundary of a shear-thickening cornstarch suspension revealing persistent regions of high local stress propagating in the flow direction at the speed of the top boundary. The persistence of these propagating fronts enables precise measurements of their structure, including the profile of boundary stress measured by boundary stress microscopy (BSM) and the nonaffine velocity of particles at the bottom boundary of the suspension measured by particle image velocimetry (PIV). In addition, we directly measure the relative flow between the particle phase and the suspending fluid (fluid migration) and find the migration is highly localized to the fronts and changes direction across the front, indicating that the fronts are composed of a localized region of high dilatant pressure and low particle concentration. The magnitude of the flow indicates that the pore pressure difference driving the fluid migration is comparable to the critical shear stress for the onset of shear thickening. The propagating fronts fully account for the increase in viscosity with applied stress reported by the rheometer and are consistent with the existence of a stable jammed region in contact with one boundary of the system that generates a propagating network of percolated frictional contacts spanning the gap between the rheometer plates and producing strong localized dilatant pressure.more » « less
-
Ice at depth in ice-stream shear margins is thought to commonly be temperate, with interstitial meltwater that softens ice. Models that include this softening extrapolate results of a single experimental study in which ice effective viscosity decreased by a factor of ∼3 over water contents of ∼0.01–0.8%. Modeling indicates this softening by water localizes strain in shear margins and through shear heating increases meltwater at the bed, enhancing basal slip. To extend data to higher water contents, we shear lab-made ice in confined compression with a large ring-shear device. Ice rings with initial mean grain sizes of 2–4 mm are kept at the pressure-melting temperature and sheared at controlled rates with peak stresses of ∼0.06–0.20 MPa, spanning most of the estimated shear-stress range in West Antarctic shear margins. Final mean grain sizes are 8–13 mm. Water content is measured by inducing a freezing front at the ice-ring edges, tracking its movement inward with thermistors, and fitting the data with solutions of the relevant Stefan problem. Results indicate two creep regimes, below and above a water content of ∼0.6%. Comparison of effective viscosity values in secondary creep with those of tertiary creep from the earlier experimental study indicate that for water contents of 0.2–0.6%, viscosity in secondary creep is about twice as sensitive to water content than for ice sheared to tertiary creep. Above water contents of 0.6%, viscosity values in secondary creep are within 25% of those of tertiary creep, suggesting a stress-limiting mechanism at water contents greater than 0.6% that is insensitive to ice fabric development in tertiary creep. At water contents of ∼0.6–1.7%, effective viscosity is independent of water content, and ice is nearly linear-viscous. Minimization of intercrystalline stress heterogeneity by grain-scale melting and refreezing at rates that approach an upper bound as grain-boundary water films thicken might account for the two regimes.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/D0SM01917H
