We explore the rheology predicted by a recently proposed constitutive model for jammed suspensions of soft elastic particles derived from microscopic dynamics [Cuny et al., arXiv:2102.05938]. Our model predicts that the orientation of the anisotropy
of the microstructure, governed by an interplay between flow vorticity and contact elasticity, plays a key role at yielding and in flow. It generates normal stress differences contributing significantly to the yield criterion and Trouton ratio. It gives rise to non-trivial transients such as stress overshoots in step increases of shear rates, residual stresses after flow cessation and power law decay of the shear rate in creep. Finally, it explains the collapse of storage modulus as measured in parallel superposition for a yielded suspension.
Dense suspensions of particles are relevant to many applications and are a key platform for developing a fundamental physics of out-of-equilibrium systems. They present challenging flow properties, apparently turning from liquid to solid upon small c
hanges in composition or, intriguingly, in the driving forces applied to them. The emergent physics close to the ubiquitous jamming transition (and to some extent the glass and gelation transitions) provides common principles with which to achieve a consistent interpretation of a vast set of phenomena reported in the literature. In light of this, we review the current state of understanding regarding the relation between the physics at the particle scale and the rheology at the macroscopic scale. We further show how this perspective opens new avenues for the development of continuum models for dense suspensions.
The origin of the abrupt shear thickening observed in some dense suspensions has been recently argued to be a transition from frictionless (lubricated) to frictional interactions between immersed particles. The Wyart-Cates rheological model, built on
this scenario, introduced the concept of fraction of frictional contacts $f$ as the relevant order parameter for the shear thickening transition. Central to the model is the equation-of-state relating $f$ to the applied stress $sigma$, which is directly linked to the distribution of the normal components of non-hydrodynamics interparticle forces. Here, we develop a model for this force distribution, based on the so-called $q$-model that we borrow from granular physics. This model explains the known $f(sigma)$ in the simple case of sphere contacts displaying only sliding friction, but also predicts strong deviation from this usual form when stronger kinds of constraints are applied on relative motion. We verify these predictions in the case of contacts with rolling friction, in particular a broadening of the stress range over which shear thickening occurs. We finally discuss how a similar approach can be followed to predict $f(sigma)$ in systems with other variations from the canonical system of monodisperse spheres with sliding friction, in particular the case of large bidispersity.
Discrete particle simulations are used to study the shear rheology of dense, stabilized, frictional particulate suspensions in a viscous liquid, toward development of a constitutive model for steady shear flows at arbitrary stress. These suspensions
undergo increasingly strong continuous shear thickening (CST) as solid volume fraction $phi$ increases above a critical volume fraction, and discontinuous shear thickening (DST) is observed for a range of $phi$. When studied at controlled stress, the DST behavior is associated with non-monotonic flow curves of the steady-state stress as a function of shear rate. Recent studies have related shear thickening to a transition between mostly lubricated to predominantly frictional contacts with the increase in stress. In this study, the behavior is simulated over a wide range of the dimensionless parameters $(phi,tilde{sigma}$, and $mu)$, with $tilde{sigma} = sigma/sigma_0$ the dimensionless shear stress and $mu$ the coefficient of interparticle friction: the dimensional stress is $sigma$, and $sigma_0 propto F_0/ a^2$, where $F_0$ is the magnitude of repulsive force at contact and $a$ is the particle radius. The data have been used to populate the model of the lubricated-to-frictional rheology of Wyart and Cates [Phys. Rev. Lett.{bf 112}, 098302 (2014)], which is based on the concept of two viscosity divergences or textquotedblleft jammingtextquotedblright points at volume fraction $phi_{rm J}^0 = phi_{rm rcp}$ (random close packing) for the low-stress lubricated state, and at $phi_{rm J} (mu) < phi_{rm J}^0$ for any nonzero $mu$ in the frictional state; a generalization provides the normal stress response as well as the shear stress. A flow state map of this material is developed based on the simulation results.
Dynamic particle-scale numerical simulations are used to show that the shear thickening observed in dense colloidal, or Brownian, suspensions is of a similar nature to that observed in non-colloidal suspensions, i.e., a stress-induced transition from
a flow of lubricated near-contacting particles to a flow of a frictionally contacting network of particles. Abrupt (or discontinuous) shear thickening is found to be a geometric rather than hydrodynamic phenomenon; it stems from the strong sensitivity of the jamming volume fraction to the nature of contact forces between suspended particles. The thickening obtained in a colloidal suspension of purely hard frictional spheres is qualitatively similar to experimental observations. However, the agreement cannot be made quantitative with only hydrodynamics, frictional contacts and Brownian forces. Therefore the role of a short-range repulsive potential mimicking the stabilization of actual suspensions on the thickening is studied. The effects of Brownian and repulsive forces on the onset stress can be combined in an additive manner. The simulations including Brownian and stabilizing forces show excellent agreement with experimental data for the viscosity $eta$ and the second normal stress difference $N_2$.
The discontinuous shear thickening (DST) of dense suspensions is a remarkable phenomenon in which the viscosity can increase by several orders of magnitude at a critical shear rate. It has the appearance of a first order phase transition between two
hypothetical states that we have recently identified as Stokes flows with lubricated or frictional contacts, respectively. Here we extend the analogy further by means of novel stress-controlled simulations and show the existence of a non-monotonic steady-state flow curve analogous to a non-monotonic equation of state. While we associate DST with an S-shaped flow curve, at volume fractions above the shear jamming transition the frictional state loses flowability and the flow curve reduces to an arch, permitting the system to flow only at small stresses. Whereas a thermodynamic transition leads to phase separation in the coexistence region, we observe a uniform shear flow all along the thickening transition. A stability analysis suggests that uniform shear may be mechanically stable for the small Reynolds numbers and system sizes in a rheometer.
Particles suspended in a Newtonian fluid raise the viscosity and also generally give rise to a shear-rate dependent rheology. In particular, pronounced shear thickening may be observed at large solid volume fractions. In a recent article (R. Seto, R.
Mari, J. F. Morris, and M. M. Denn., Phys. Rev. Lett., 111:218301, 2013) we have considered the minimum set of components to reproduce the experimentally observed shear thickening behavior, including Discontinuous Shear Thickening (DST). We have found frictional contact forces to be essential, and were able to reproduce the experimental behavior by a simulation including this physical ingredient along with viscous lubrication. In the present article, we thoroughly investigate the effect of friction and express it in the framework of the jamming transition. The viscosity divergence at the jamming transition has been a well known phenomenon in suspension rheology, as reflected in many empirical laws for the viscosity. Friction can affect this divergence, and in particular the jamming packing fraction is reduced if particles are frictional. Within the physical description proposed here, shear thickening is a direct consequence of this effect: as the shear rate increases, friction is increasingly incorporated as more contacts form, leading to a transition from a mostly frictionless to a mostly frictional rheology. This result is significant because it shifts the emphasis from lubrication hydrodynamics and detailed microscopic interactions to geometry and steric constraints close to the jamming transition.
Discontinuous shear thickening (DST) observed in many dense athermal suspensions has proven difficult to understand and to reproduce by numerical simulation. By introducing a numerical scheme including both relevant hydrodynamic interactions and gran
ularlike contacts, we show that contact friction is essential for having DST. Above a critical volume fraction, we observe the existence of two states: a low viscosity, contactless (hence, frictionless) state, and a high viscosity frictional shear jammed state. These two states are separated by a critical shear stress, associated with a critical shear rate where DST occurs. The shear jammed state is reminiscent of the jamming phase of granular matter. Continuous shear thickening is seen as a lower volume fraction vestige of the jamming transition.
