Shear thickening and history-dependent rheology of monodisperse suspensions with finite inertia via an immersed boundary lattice Boltzmann method

Three-dimensional direct numerical simulations of dense suspensions of monodisperse spherical particles in simple shear flow have been performed at particle Reynolds numbers between 0.1 and 0.6. The particles translate and rotate under the influence of the applied shear. The lattice Boltzmann method was used to solve the flow of the interstitial Newtonian liquid, and an immersed boundary method was used to enforce the no-slip boundary condition at the surface of each particle. Short range spring forces were applied between colliding particles over sub-grid scale distances to prevent overlap. We computed the relative apparent viscosity for solids volume fractions up to 38% for several shear rates and particle concentrations and discuss the effects of these variables on particle rotation and cluster formations. The apparent viscosities increase with increasing particle Reynolds number (shear thickening) and solids fraction. As long as the particle Reynolds number is low (0.1), the computed viscosities are in good agreement with experimental measurements, as well as theoretical and empirical equations. For higher Reynolds numbers, we find much higher viscosities, which we relate to slower particle rotation and clustering. Simulations with a sudden change in shear rate also reveal a history (or hysteresis) effect due to the formation of clusters. We quantify the changes in particle rotation and clustering as a function of the Reynolds number and volume fraction.