Rheology of non-Brownian suspensions of rough frictional particles under shear reversal: A numerical study

Abstract

We perform particle scale simulations of suspensions submitted to shear reversal. The simulations are based on the Force Coupling method, adapted to account for short range lubrication interactions together with direct contact forces between particles, including surface roughness, contact elasticity, and solid friction. After shear reversal, three consecutive steps are identified in the viscosity transient: An instantaneous variation, followed by a rapid contact force relaxation, and finally a long time evolution. The separated contributions of hydrodynamics and contact forces to the viscosity are investigated during the transient, allowing a qualitative understanding of each step. In addition, the influence of the contact law parameters (surface roughness height and friction coefficient) on the transient is evaluated. Concerning the long time transient, the difference between the steady viscosity and minimum viscosity is shown to be proportional to the contact contribution to the steady viscosity, allowing in principle easy determination of the latter in experiments. The short time evolution is studied as well. After the shear reversal, the contact forces vanish over a strain that is very short compared to the typical strain of the long time transient, allowing to define an apparent step between the viscosity before shear reversal and after contact force relaxation. This step is shown to be an increasing function of the friction coefficient between particles. Two regimes are identified as a function of the volume fraction. At low volume fraction, the step is small compared to the steady contact viscosity, in agreement with a particle pair model. As the volume fraction increases, the value of the viscosity step increases faster than the steady contact viscosity, and, depending on the friction coefficient, may approach it.