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   
                    This content will become publicly available on January 1, 2026
                            
                            Stress and flow inhomogeneity in shear-thickening suspensions
                        
                    
    
            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   
        
    
                            - Award ID(s):
- 2226485
- PAR ID:
- 10563490
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- Journal of Colloid and Interface Science
- Volume:
- 678
- Issue:
- PA
- ISSN:
- 0021-9797
- Page Range / eLocation ID:
- 218 to 225
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            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
- 
            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
- 
            Nearly, all dense suspensions undergo dramatic and abrupt thickening transitions in their flow behavior when sheared at high stresses. Such transitions occur when the dominant interactions between the suspended particles shift from hydrodynamic to frictional. Here, we interpret abrupt shear thickening as a precursor to a rigidity transition and give a complete theory of the viscosity in terms of a universal crossover scaling function from the frictionless jamming point to a rigidity transition associated with friction, anisotropy, and shear. Strikingly, we find experimentally that for two different systems—cornstarch in glycerol and silica spheres in glycerol—the viscosity can be collapsed onto a single universal curve over a wide range of stresses and volume fractions. The collapse reveals two separate scaling regimes due to a crossover between frictionless isotropic jamming and frictional shear jamming, with different critical exponents. The material-specific behavior due to the microscale particle interactions is incorporated into a scaling variable governing the proximity to shear jamming, that depends on both stress and volume fraction. This reformulation opens the door to importing the vast theoretical machinery developed to understand equilibrium critical phenomena to elucidate fundamental physical aspects of the shear thickening transition.more » « less
- 
            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
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
