An official website of the United States government
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
Flow induced rigidity percolation in shear thickening suspensions
Discontinuous shear thickening (DST) is associated with a sharp rise in a suspension’s viscosity with increasing applied shear rate or stress. Key signatures of DST, highlighted in recent studies, are the very large fluctuations of the measured stress as the suspension thickens with increasing rate. A clear link between microstructural development and the dramatic increase in stress fluctuations has not been established yet. To identify the microstructural underpinnings of this behavior, we perform simulations of sheared dense suspensions. Through an analysis of the particle contact network, we identify a subset of constrained particles that contributes directly to the rapid rise in viscosity and large stress fluctuations. Indeed, both phenomena can be explained by the growth and percolation of constrained particle networks—in direct analogy to rigidity percolation. A finite size scaling analysis confirms this to be a percolation phenomenon and allows us to estimate the critical exponents. Our findings reveal the specific microstructural self-organization transition that underlies DST.
more »
« less
- Award ID(s):
- 2026842
- PAR ID:
- 10573824
- Publisher / Repository:
- AIPP
- Date Published:
- Journal Name:
- Journal of Rheology
- Volume:
- 68
- Issue:
- 2
- ISSN:
- 0148-6055
- Page Range / eLocation ID:
- 219 to 228
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Dense suspensions of particles in viscous liquid often demonstrate the striking phenomenon of abrupt shear thickening, where their viscosity increases strongly with increase of the imposed stress or shear rate. In this work, discrete-particle simulations accounting for short-range hydrodynamic, repulsive, and contact forces are performed to simulate flow of shear thickening bidisperse suspensions, with the packing parameters of large-to-small particle radius ratio δ = 3 and large particle fraction ζ = 0.15, 0.50, and 0.85. The simulations are carried out for volume fractions 0.54 ≤ ϕ ≤ 0.60 and a wide range of shear stresses. The repulsive forces, of magnitude F R , model the effects of surface charge and electric double-layer overlap, and result in shear thinning at small stress, with shear thickening beginning at stresses σ ∼ F R a −2 . A crossover scaling analysis used to describe systems with more than one thermodynamic critical point has recently been shown to successfully describe the experimentally-observed shear thickening behavior in suspensions. The scaling theory is tested here on simulated shear thickening data of the bidisperse mixtures, and also on nearly monodisperse suspensions with δ = 1.4 and ζ = 0.50. Presenting the viscosity in terms of a universal crossover scaling function between the frictionless and frictional maximum packing fractions collapses the viscosity for most of the suspensions studied. Two scaling regimes having different exponents are observed. The scaling analysis shows that the second normal stress difference N 2 and the particle pressure Π also collapse on their respective curves, with the latter featuring a different exponent from the viscosity and normal stress difference. The influence of the fraction of frictional contacts, one of the parameters of the scaling analysis, and its dependence on the packing parameters are also presented.more » « less
-
Lithium-ion battery cathodes are manufactured by coating slurries, liquid suspensions that typically include carbon black (CB), active material, and polymer binder. These slurries have a yield stress and complex rheology due to CB’s microstructural response to flow. While optimizing the formulation and processing of slurries is critical to manufacturing defect-free and high-performance cathodes, engineering the shear rheology of cathode slurries remains challenging. In this study, we conducted simultaneous rheo-electric measurements on 3 wt% CB suspensions in N-methyl-2-pyrrolidone containing various loadings of active material NMC811 and polyvinylidene difluoride. Accounting for the changes in the infinite shear viscosity, the yield stress, and the medium viscosity due to the presence of NMC and polymers, we defined the differential relative viscosity. This differential relative viscosity, Δ𝜂𝑟, is a measure of the distance from the infinite shear rate, where carbon black agglomerates are fully broken down. We find that Δ𝜂𝑟 collapses all flow curves regardless of formulation with an empirical relationship Δ𝜂𝑟=2.18𝑀𝑛𝑓−0.92, indicating a quantitative prediction of the flow curve of cathode slurries across a wide range of formulation space. We then used electrical conductivity to identify and quantify shear-induced structure memory, evidenced in the ratio of the shear conductivity over the post-shear quiescent conductivity. We find that similar to the changes in the yield stress, increasing NMC concentration increases memory retention, and in contrast, the addition of PVDF erases memory effects. Our findings here will provide valuable insight into engineering the formulation and processing conditions of lithium-ion battery cathodes.more » « less
-
Abstract Fragmentation plays a critical role in eruption explosivity by influencing the eruptive jet and plume dynamics that may initiate hazards such as pyroclastic flows. The mechanics and progression of fragmentation during an eruption are challenging to constrain observationally, limiting our understanding of this important process. In this work, we explore seismic radiation associated with unsteady fragmentation. Seismic force and moment tensor fluctuations from unsteady fragmentation arise from fluctuations in fragmentation depth and wall shear stress (e.g., from viscosity variations). We use unsteady conduit flow models to simulate perturbations to a steady‐state eruption from injections of heterogeneous magma (specifically, variable magma viscosity due to crystal volume fraction variations). Changes in wall shear stress and pressure determine the seismic force and moment histories, which are used to calculate synthetic seismograms. We consider three heterogeneity profiles: Gaussian pulse, sinusoidal, and stochastic. Fragmentation of a high‐crystallinity Gaussian pulse produces a distinct very‐long‐period seismic signature and associated reduction in mass eruption rate, suggesting joint use of seismic, infrasound, and plume monitoring data to identify this process. Simulations of sinusoidal injections quantify the relation between the frequency or length scale of heterogeneities passing through fragmentation and spectral peaks in seismograms, with velocity seismogram amplitudes increasing with frequency. Stochastic composition variations produce stochastic seismic signals similar to observed eruption tremor, though computational limitations restrict our study to frequencies less than 0.25 Hz. We suggest that stochastic fragmentation fluctuations could be a plausible eruption tremor source.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1122/8.0000786
