Inclusion of DLVO forces in simulations of non-Brownian solid suspensions: Rheology and structure

The understanding of the rheological behaviour of suspensions in aqueous electrolytes is necessary for the optimal design of hydraulic transport lines. In these applications, particle size is at least 10 micron, and the particle Reynolds number, Rep, is finite: O(10−1). Although there are some experimental and numerical data on the rheology of such suspensions, the number of detailed analyses is limited. Therefore, 3-D direct numerical simulations of dense suspensions in aqueous electrolytes are conducted to assess the dynamics of the relative apparent viscosity and particle structures. The solid–liquid interfaces are resolved, and the flow is simulated, employing an in-house immersed boundary-lattice Boltzmann method code. In addition to the hydrodynamics resolved in the computational grid, our simulations include unresolved sub-grid scale lubrication corrections and non-contact electric double layer (EDL) and Van der Waals forces for a wide range of particle volume fractions, ϕv, at a single Rep=0.1. Under these conditions, the contribution of the Van der Waals force was found to be weak. With an increase in ϕv, the effect of EDL forces decreased the relative apparent viscosity. Particle layering and structural arrangements were analysed for ϕv=43 and 52%. As the Debye length (i.e., the thickness of EDL) decreases, the particle layers near the walls weakened. The analyses reveal how at these high volume fractions, chain-like assemblies are transformed into hexagonal arrangements. •For the parameter values considered, Van der Waals forces have a weak effect on suspension rheology.•Electric double layer (EDL) forces reduce both the relative apparent viscosity and the rate of cluster formation.•EDL forces contribute to particle layering and alter morphology from chain-like structures to regular hexagons.