Numerical simulations of dense granular suspensions in laminar flow under constant and varying shear rates

Using an immersed boundary-lattice Boltzmann method, we investigated the response of dense granular suspensions to time-varying shear rates and flow reversals. The evolution of the relative apparent viscosity and particle structures was analysed. The concentration of solids (phi(v)) and particle Reynolds numbers (Re-p) were varied over the ranges 6% <= phi(v) <= 47% and 0.105 <= Re-p <= 0.529. The simulations included sub-grid scale corrections for unresolved lubrication forces and torques (normal and tangential). When phi(v) surpasses 30%, the contribution of the tangential lubrication corrections to the shear stress is dominant. While for intermediate solids fractions we find weak shear-thinning, we see weak shear-thickening for phi(v) > 40%. We show how the structure and apparent viscosity of a suspension evolves after a reversal of the shear direction. For 47% solids, simulations with step changes in the shear rate show the effects of the previous shear history on the viscosity of the suspension.

Automated summary

The research utilized an immersed boundary-lattice Boltzmann method to study the response of dense granular suspensions to time-varying shear rates and flow reversals. The results showed that the contribution of tangential lubrication corrections to shear stress is dominant when solid concentration exceeds 30%. Weak shear-thinning was found for intermediate solids fractions, while weak shear-thickening was observed for solid concentrations above 40%.