skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Strain-Controlled Critical Slowing Down in the Rheology of Disordered Networks
Networks and dense suspensions frequently reside near a boundary between soft (or fluidlike) and rigid (or solidlike) regimes. Transitions between these regimes can be driven by changes in structure, density, or applied stress or strain. In general, near the onset or loss of rigidity in these systems, dissipation-limiting heterogeneous nonaffine rearrangements dominate the macroscopic viscoelastic response, giving rise to diverging relaxation times and power-law rheology. Here, we describe a simple quantitative relationship between nonaffinity and the excess viscosity. We test this nonaffinity-viscosity relationship computationally and demonstrate its rheological consequences in simulations of strained filament networks and dense suspensions. We also predict critical signatures in the rheology of semiflexible and stiff biopolymer networks near the strain stiffening transition.  more » « less
Award ID(s):
2019745 2224030
PAR ID:
10512524
Author(s) / Creator(s):
; ;
Publisher / Repository:
American Physical Society
Date Published:
Journal Name:
Physical Review Letters
Volume:
131
Issue:
17
ISSN:
0031-9007
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. As a paradigmatic model of active fluids, bacterial suspensions show intriguing rheological responses drastically different from their counterpart colloidal suspensions. Although the flow of bulk bacterial suspensions has been extensively studied, the rheology of bacterial suspensions under confinement has not been experimentally explored. Here, using a microfluidic viscometer, we systematically measure the rheology of dilute Escherichia coli suspensions under different degrees of confinement. Our study reveals a strong confinement effect: the viscosity of bacterial suspensions decreases substantially when the confinement scale is comparable or smaller than the run length of bacteria. Moreover, we also investigate the microscopic dynamics of bacterial suspensions including velocity profiles, bacterial density distributions, and single bacterial dynamics in shear flows. These measurements allow us to construct a simple heuristic model based on the boundary layer of upstream swimming bacteria near confining walls, which qualitatively explains our experimental observations. Our study sheds light on the influence of the boundary layer of collective bacterial motions on the flow of confined bacterial suspensions. Our results provide a benchmark for testing different rheological models of active fluids and are useful for understanding the transport of microorganisms in confined geometries. 
    more » « less
  2. 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
  3. 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
  4. Debris flows are dense and fast-moving complex suspensions of soil and water that threaten lives and infrastructure. Assessing the hazard potential of debris flows requires predicting yield and flow behavior. Reported measurements of rheology for debris flow slurries are highly variable and sometimes contradictory due to heterogeneity in particle composition and volume fraction ( ฯ• ) and also inconsistent measurement methods. Here we examine the composition and flow behavior of source materials that formed the postwildfire debris flows in Montecito, CA, in 2018, for a wide range of ฯ• that encapsulates debris flow formation by overland flow. We find that shear viscosity and yield stress are controlled by the distance from jamming, ฮ” ฯ• = ฯ• m โˆ’ ฯ• , where the jamming fraction ฯ• m is a material parameter that depends on grain size polydispersity and friction. By rescaling shear and viscous stresses to account for these effects, the data collapse onto a simple nondimensional flow curve indicative of a Bingham plastic (viscoplastic) fluid. Given the highly nonlinear dependence of rheology on ฮ” ฯ• , our findings suggest that determining the jamming fraction for natural materials will significantly improve flow models for geophysical suspensions such as hyperconcentrated flows and debris flows. 
    more » « less
  5. 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