Consititutive modeling of time-dependent flows in frictional suspensions

This paper summarizes recent joint work towards a constitutive modelling framework for dense granular suspensions. The aim is to create a time-dependent, tensorial theory that can implement the physics described in steady state by the Wyart-Cates model. This model of shear thickening suspensions supposes that lubrication films break above a characteristic normal force so that frictional contact forces come into play: the resulting non-sliding constraints can be enough to rigidify a system that would flow freely at lower stresses [1]. Implementing this idea for time-dependent flows requires the introduction of new concepts including a configuration-dependent ‘jamming coordinate’, alongside a decomposition of the velocity gradient tensor into compressive and extensional components which then enter the evolution equation for particle contacts in distinct ways. The resulting approach [2, 3] is qualitatively successful in addressing (i) the collapse of stress during flow reversal in shear flow, and (ii) the ability of transverse oscillatory flows to unjam the system. However there is much work required to refine this approach towards quantitative accuracy, by incorporating more of the physics of contact evolution under flow as determined by close interrogation of particle-based simulations.

Transients in pressure-imposed shearing of dense granular suspensions

Granular materials whether dry or immersed in fluid show dilation or compaction depending upon the initial conditions, solid fraction and normal stress. Here we probe the transient response of a dense granular suspension subjected to change of applied normal stress under simple shear. In this aim, normal-stress-imposed discrete element particle simulations are developed considering the contributions arising from the drag induced on the particles by fluid phase. These pressure-imposed simulations show transient behaviors of dense granular suspensions such as dilation or compaction before reaching a steady state following the µ(J) rheology. Less expectedly, the transient behavior, in particular the height of the system as a function of applied strain, can also be described by assuming that the system follows the steady µ(J) rheology at all times.