Competing Timescales Lead to Oscillations in Shear-Thickening Suspensions

On pressure-driven Poiseuille flow with non-monotonic rheology

Shear thickening fluids are liquids that stiffen as the applied stress increases. If many of these types of fluids follow a monotonic rheological curve, some experimental and numerical studies suggest that certain fluids, like cornstarch, may exhibit a non-monotonic, S-shaped rheology. Such non-monotonic behavior has however proved very difficult to observe experimentally in classical rheometer. To explain such difficulties, the possible presence of vorticity banding in the rheometer has been considered. To prevent such instabilities, we use a capillary rheometer, which is a cylindrical tube, measuring the flow rate versus the applied pressure drop. With this setup, we indeed observe a non-monotonic behavior: the flow rate increases monotonically at low pressure drops up to a maximum, after which it abruptly decreases to an almost constant flow rate regardless of further increases in pressure drop. This maximum-jump-plateau behavior occurs over a wide range of concentrations and is reproducible without hysteresis, which is in agreement with an S-shaped rheology. However, the obtained flow versus pressure difference function Q ( Δ P ) does not agree with the classical Wyart-Cates rheological model, which predicts an S-shaped non-monotonic function, but with neither a jump nor a plateau. To understand this jump-plateau behavior, we remark that any rheological model would establish a relationship between the flow rate and the local pressure gradient, but not the total pressure drop. We thus discuss and analyze the implications of having an S-shaped non-monotonic flow rate-pressure gradient in Poiseuille flow. In particular, we discuss the possibility of a non-uniform pressure gradient in the direction of the flow, i.e., a kind of streamwise banding. The key issue is then the selection of the gradient pressure distribution along the tube. One solution could arise from an analogy of this problem with the spinodal decomposition. It, however, leads to an increase in flow rate with∂ x P up to a plateau between two values of∂ x P as determined by the Maxwell construction. To account for the bump-jump behavior, we have implemented a simple dynamical stochastic version of the Wyart-Cates model, where the thickening occurs with a characteristic time. As a result, with increasing the total pressure drop, the flow rate increases monotonically up to a maximum value. Beyond this point, the flow rate drops abruptly to a lower value, forming a slowly decreasing plateau. This behavior is likely to account for the maximum-jump-plateau observed in the experiments. We also show that in such a system, the final state is quite sensitive to the initial state of the fluid, especially its homogeneity. Our results then demonstrate that the mere presence of a non-monotonic rheological curve is sufficient to predict the occurrence of stress banding in the streamwise direction and a plateau flow rate, even if the suspension remains homogeneous.

A frictional soliton controls the resistance law of shear-thickening suspensions in pipes

Pipe flows are commonly found in nature and industry as an effective mean of transporting fluids. They are primarily characterized by their resistance law, which relates the mean flow rate to the driving pressure gradient. Since Poiseuille and Hagen, various flow regimes and fluid rheologies have been investigated, but the behavior of shear-thickening suspensions, which jam above a critical shear stress, remains poorly understood despite important applications (e.g., concrete or food processing). In this study, we build on recent advances in the physics of shear-thickening suspensions to address their flow through pipes and establish their resistance law. We find that for discontinuously shear-thickening suspensions (large particule volume fractions), the flow rate saturates at high driving stress. Local pressure and velocity measurements reveal that this saturation stems from the emergence of a frictional soliton: a unique, localized, superdissipative, and backpropagating flow structure coexisting with the laminar frictionless flow phase observed at low driving stress. We characterize the remarkably steep effective rheology of the frictional soliton and show that it sets the resistance law at the whole pipe scale. These findings offer an unusual perspective on low-Reynolds suspension flows through pipes, intriguingly reminiscent of the transition to turbulence for simple fluids. They also provide a predictive law for the transport of such suspensions in pipe systems, with implications for a wide range of applications.

Optimizing non-Newtonian fluids for impact protection of laminates

Non-Newtonian fluids can be used for the protection of flexible laminates. Understanding the coupling between the flow of the protecting fluid and the deformation of the protected solids is necessary in order to optimize this functionality. We present a scaling analysis of the problem based on a single coupling variable, the effective width of a squeeze flow between flat rigid plates, and predict that impact protection for laminates is optimized by using shear-thinning, and not shear-thickening, fluids. The prediction is verified experimentally by measuring the velocity and pressure in impact experiments. Our scaling analysis should be generically applicable for non-Newtonian fluid-solid interactions in diverse applications.

Shear thickening in presence of adhesive contact forces: The singularity of cornstarch

HYPOTHESIS

A number of dense particle suspensions experience a dramatic increase in viscosity with the shear stress, up to a solid-like response. This shear-thickening process is understood as a transition under flow of the nature of the contacts - from lubricated to frictional - between initially repellent particles. Most systems are now assumed to fit in with this scenario, which is questionable.

EXPERIMENT

Using an in-house pressure sensor array, we provide a spatio-temporal map of the normal stresses in the flows of two shear-thickening fluids: a stabilized calcium carbonate suspension, known to fit in with the standard scenario, and a cornstarch suspension, which spectacular thickening behavior remains poorly understood.

FINDINGS

We evidence in cornstarch a unique, stable heterogeneous structure, which moves in the velocity direction and does not appear in calcium carbonate. Its nature changes from a stress wave to a rolling solid jammed aggregate at high solid fraction and small gap width. The modeling of these heterogenities points to an adhesive force between cornstarch particles at high stress, also evidenced in microscopic measurements. Cornstarch being also attractive at low stress, it stands out of the classical shear-thickening frame, and might be part of a larger family of adhesive and attractive shear-thickening fluids.

Structure of propagating high-stress fronts in a shear-thickening suspension

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.

Origin of Two Distinct Stress Relaxation Regimes in Shear Jammed Dense Suspensions

Many dense particulate suspensions show a stress induced transformation from a liquidlike state to a solidlike shear jammed (SJ) state. However, the underlying particle-scale dynamics leading to such striking, reversible transition of the bulk remains unknown. Here, we study transient stress relaxation behaviour of SJ states formed by a well-characterized dense suspension under a step strain perturbation. We observe a strongly nonexponential relaxation that develops a sharp discontinuous stress drop at short time for high enough peak-stress values. High resolution boundary imaging and normal stress measurements confirm that such stress discontinuity originates from the localized plastic events, whereas system spanning dilation controls the slower relaxation process. We also find an intriguing correlation between the nature of transient relaxation and the steady-state shear jamming phase diagram obtained from the Wyart-Cates model.

Viscoelasticity of non-colloidal hydrogel particle suspensions at the liquid-solid transition

Suspensions of soft particles transition from a viscous fluid to a soft material upon increases in phase volume. The criteria defining the transition to this jammed state are difficult to define due to the porous and deformable nature of soft particles. Here, we characterise the rheology of aqueous suspensions of industrially relevant non-colloidal, polydisperse, frictional agarose microgels and evaluate shear and viscoelastic behaviour across a range of phase volumes from the dilute regime to the highly concentrated regime. In order to model the viscoelastic response of suspensions without free fitting parameters, the random close packing volume fraction (φrcp) and the particle modulus are determined, respectively, from particle size distribution measurements and direct measurements of reduced elastic modulus of individual particles (Erp) using Atomic Force Microscopy. It is found that at φrcp, previously shown to correspond to divergence of the viscosity, also corresponds to the suspension transition from a viscous to viscoelastic fluid. However, the transition to a jammed solid-like state (φj) occurs at phase volumes exceeding this value (i.e. φj > φrcp). The suspension modulus and its sudden growth at φj are well-predicted by the Evans and Lips model that incorporates the Erp of the hydrogel particles. This rheological behaviour showing a dual transition is reminiscent of two families of systems: (i) colloidal suspensions and (ii) frictional-adhesive non-colloidal suspensions. However, it does not strictly follow either case. We propose that the width of the transition region is dictated by frictional contact, particle size distribution and particle modulus, and plan to further probe this in future work.