Shear viscosity scaling of granular suspensions across dilute to dense regimes

Zaohui Zhang,Teng Man,Herbert E. Huppert,Sergio Andres Galindo-Torres
2023-07-17
Abstract:In this letter, following an extensive experimental validation, we perform constant-volume shearing simulations of non-Brownian granular suspensions using the discrete element method coupled with the lattice Boltzmann method. We choose a wide range of solid fractions, shear rates, fluid viscosities, particle sizes, and inter-particle frictional coefficients to obtain a scaling solution for the viscous behavior of suspensions in both dilute and dense regimes. This letter demonstrates that, with a proposed dilute-dense transitional solid fraction, $\phi_d$, there exists a strong correlation between the inverse relative viscosity and the shear stress. This work incorporates both the $\phi$-dependence and the $\dot{\gamma}$-dependence of suspension viscosity in a universal framework, which provides a scaling solution for granular suspensions across dilute and dense regimes and sheds light on the dilute-dense transition mechanisms.
Soft Condensed Matter
What problem does this paper attempt to address?
This paper aims to solve the problem of shear viscosity scaling of non - Brownian particle suspensions in the sparse - to - dense state. Specifically, through experimental verification and numerical simulation, the researchers explored the effects of different solid - phase fractions (\(\phi\)), shear rates (\(\dot{\gamma}\)), fluid viscosities (\(\eta_f\)), particle sizes (\(D_p\)) and inter - particle friction coefficients (\(\mu_p\)) on the viscosity behavior of particle suspensions. The main goal of the paper is to establish a general model that can describe the viscosity change of particle suspensions from the sparse to the dense state and reveal the sparse - dense transition mechanism. ### Main research content 1. **Numerical simulation method**: The method of coupling the discrete element method (DEM) with the lattice Boltzmann method (LBM) was used for constant - volume shear simulation, and the simulation object was non - Brownian particle suspensions. 2. **Parameter range**: A wide range of solid - phase fractions, shear rates, fluid viscosities, particle sizes and inter - particle friction coefficients were selected to obtain the scaling solution of the viscosity behavior of suspensions. 3. **Key findings**: - A critical solid - phase fraction \(\phi_d\) was proposed to distinguish between the sparse and dense states. - A strong correlation was found between the shear stress \(\tau\) and the reciprocal of the relative viscosity \(\eta_r\). - A general framework was established, taking into account the effects of both the solid - phase fraction \(\phi\) and the shear rate \(\dot{\gamma}\) on the viscosity of suspensions. ### Research background Particle suspensions are of great significance in industrial production and environmental protection, but it is still challenging to establish a general model to describe their behavior. During the transition from the sparse to the dense state of particle suspensions, they will experience a change from Newtonian fluids to non - Newtonian fluids, and finally, the jamming phenomenon will occur at the critical solid - phase fraction \(\phi_c\). Traditional viscosity - scaling methods can usually only describe the behavior in the sparse state and cannot fully explain the transition mechanism from the sparse to the dense state. ### Research methods 1. **Numerical simulation**: The DEM - LBM coupling method was used for simulation, and the simulation parameters included the solid - phase fraction \(\phi\), the shear rate \(\dot{\gamma}\), the fluid viscosity \(\eta_f\), the particle size \(D_p\) and the inter - particle friction coefficient \(\mu_p\). 2. **Experimental verification**: The accuracy of the numerical simulation results was verified by experimental data. 3. **Data analysis**: The variation laws of the shear stress \(\tau\) and the relative viscosity \(\eta_r\) with the solid - phase fraction \(\phi\) and the shear rate \(\dot{\gamma}\) were analyzed. ### Main conclusions 1. **Critical solid - phase fraction \(\phi_d\)**: A critical solid - phase fraction \(\phi_d\) was proposed to distinguish between the sparse and dense states. 2. **Viscosity - scaling relationship**: The scaling relationship between the shear stress \(\tau\) and the reciprocal of the relative viscosity \(\eta_r\) was established, that is: \[ \eta_r^{-1} \sim (\phi_d - \phi)^{\beta_\eta}, \quad \phi < \phi_d \] \[ \tau_y \sim (\phi - \phi_d)^{\beta_\tau}, \quad \phi > \phi_d \] where \(\beta_\eta = 1.434\) and \(\beta_\tau = 1.476\). 3. **General framework**: A general framework was proposed, which can describe the viscosity behavior of particle suspensions in both the sparse and dense states, providing a new perspective for understanding the sparse - dense transition mechanism. In conclusion, through numerical simulation and experimental verification, this paper established a general model for describing the viscosity behavior of non - Brownian particle suspensions in the sparse - to - dense state, providing an important theoretical basis for understanding and predicting the behavior of particle suspensions.