Particles suspended in liquids, known as suspensions, are prevalent worldwide, from everyday products like creams and pastes to critical industries like pharmaceuticals and vaccines. Suspension rheology (that is, what is the suspension stress in a given flow condition) is key to many industrial processes, and is currently still an open problem. Suspensions may show many complex flow behaviors such as shear thickening (viscosity increasing with flow rate), thinning (viscosity decreasing with flow rate), and jamming (complete loss of flowability). Often, particles are large enough to not be affected by Brownian motion. External forces or driving are necessary for their motion, regardless of whether they are in a jammed or unjammed state. Therefore these concentrated suspensions exhibit mechanical properties characterized by significant dependence on previous deformation history. Suspensions of hard particles show dramatic viscosity reduction upon reversal of the flow direction (known as “shear reversal”), or under high-frequency orthogonal superposition, that is, when applying a secondary oscillating shear orthogonal to the main flow direction. This is due to the presence of a shear-induced anisotropic microstructure, which is particularly resistant to the kind of flow that gave rise to it, but not so much to flows in other directions.
Recent experiments and numerical simulations started to characterize the mechanical response of non-Brownian suspensions of hard particles to a change of flow direction. This project aims at fully characterizing the relation between flow history and mechanical response using numerical simulations and at using it to develop a constitutive model, that is, a set of equations able to predict the stress in arbitrary flow conditions. The main expected result is a quantitative constitutive model for non-Brownian suspensions, which is currently still lacking, despite decades of attempts and progress in this direction. The originality of our approach is to guide, test and validate our modelling with the rich data sets obtained from our simulations, which probe the tensorial response of the suspension in a much more complete manner than what was done previously, with e.g. shear reversal data. Such a model would enable a complete fluid mechanical model for suspensions, allowing to simulate and predict suspension flows in arbitrary geometries and flow conditions. This model would also be a basis and reference point to build richer models for more complex suspensions, e.g. for soft or adhesive particles.