We use cookies to ensure that we give you the best experience on our website. By continuing to browse this repository, you give consent for essential cookies to be used. You can read more about our Privacy and Cookie Policy.

Durham Research Online
You are in:

Shear reversal in dense suspensions : the challenge to fabric evolution models from simulation data.

Chacko, Rahul N. and Mari, Romain and Fielding, Suzanne M. and Cates, Michael E. (2018) 'Shear reversal in dense suspensions : the challenge to fabric evolution models from simulation data.', Journal of fluid mechanics., 847 . pp. 700-734.


Dense suspensions of hard particles are important as industrial or environmental materials (e.g. fresh concrete, food, paint or mud). To date, most constitutive models developed to describe them are, explicitly or effectively, ‘fabric evolution models’ based on: (i) a stress rule connecting the macroscopic stress to a second-rank microstructural fabric tensor Q; and (ii) a closed time-evolution equation for Q. In dense suspensions, most of the stress comes from short-ranged pairwise steric or lubrication interactions at near-contacts (suitably defined), so a natural choice for Q is the deviatoric second moment of the distribution P(p) of the near-contact orientations p. Here we test directly whether a closed time-evolution equation for such a Q can exist, for the case of inertialess non-Brownian hard spheres in a Newtonian solvent. We perform extensive numerical simulations accessing high levels of detail for the evolution of P(p) under shear reversal, providing a stringent test for fabric evolution models. We consider a generic class of these models as defined by Hand (J. Fluid Mech., vol. 13, 1962, pp. 33–46) that assumes little as to the micromechanical behaviour of the suspension and is only constrained by frame indifference. Motivated by the smallness of microstructural anisotropies in the dense regime, we start with linear models in this class and successively consider those increasingly nonlinear in Q. Based on these results, we suggest that no closed fabric evolution model properly describes the dynamics of the fabric tensor under reversal. We attribute this to the fact that, while a second-rank tensor captures reasonably well the microstructure in steady flows, it gives a poor description during significant parts of the microstructural evolution following shear reversal. Specifically, the truncation of P(p) at second spherical harmonic (or second-rank tensor) level describes ellipsoidal distributions of near-contact orientations, whereas on reversal we observe distributions that are markedly four-lobed; moreover, P˙(p) has oblique axes, not collinear with those of Q in the shear plane. This structure probably precludes any adequate closure at second-rank level. Instead, our numerical data suggest that closures involving the coupled evolution of both a fabric tensor and a fourth-rank tensor might be reasonably accurate.

Item Type:Article
Full text:(AM) Accepted Manuscript
Download PDF
Publisher Web site:
Publisher statement:This article has been published in a revised form in Journal of fluid mechanics This version is free to view and download for private research and study only. Not for re-distribution, re-sale or use in derivative works. © 2018 Cambridge University Press.
Date accepted:28 March 2018
Date deposited:05 April 2018
Date of first online publication:29 May 2018
Date first made open access:29 November 2018

Save or Share this output

Look up in GoogleScholar