AbstractThis paper presents three-dimensional numerical simulations of non-Brownian concentrated suspensions in a Couette flow at zero Reynolds number using a fictitious domain method. Contacts between particles are modelled using a discrete element method (DEM)-like approach, which allows for a more physical description, including roughness and friction. This work emphasizes the effect of friction between particles and its role on rheological properties, especially on normal stress differences. Friction is shown to notably increase viscosity and second normal stress difference $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}|N_2|$ and decrease $|N_1|$, in better agreement with experiments. The hydrodynamic and contact contributions to the overall particle stress are particularly investigated. This shows that the effect of friction is mostly due to the additional contact stress since the hydrodynamic stress remains unaffected by friction. Simulation results are also compared with experiments, such as normal stresses or effective friction coefficient $\mu (I_v)$, and the agreement is improved when friction is accounted for. This suggests that friction is operative in actual suspensions.

AbstractIn this paper, we conduct a numerical investigation of sheared suspensions of non-colloidal spherical particles on which a torque is applied. Particles are mono-dispersed and neutrally buoyant. Since the torque modifies particle rotation, we show that it can indeed strongly change the effective viscosity of semi-dilute or even more concentrated suspensions. We perform our calculations up to a volume fraction of 28 %. And we compare our results to data obtained at 40 % by Yeo and Maxey (Phys. Rev. E, vol. 81, 2010, p. 62501) with a totally different numerical method. Depending on the torque orientation, one can increase (decrease) the rotation of the particles. This results in a strong enhancement (reduction) of the effective shear viscosity of the suspension. We construct a dimensionless number $\Theta $ which represents the average relative angular velocity of the particles divided by the vorticity of the fluid generated by the shear flow. We show that the contribution of the particles to the effective viscosity can be suppressed for a given and unique value of $\Theta $ independently of the volume fraction. In addition, we obtain a universal behaviour (i.e. independent of the volume fraction) when we plot the relative effective viscosity divided by the relative effective viscosity without torque as a function of $\Theta $. Finally, we show that a modified Faxén law can be equivalently established for large concentrations.