Flow of wet granular materials: a numerical study

Abstract: We simulate dense assemblies of frictional spherical grains in steady shear flow under controlled normal stress $P$ in the presence of a small amount of an interstitial liquid, which gives rise to capillary menisci, assumed isolated (pendular regime), and to attractive forces. The system behavior depends on two dimensionless control parameters: inertial number $I$ and reduced pressure $P^*=aP/(\pi\Gamma)$, comparing confining forces $\sim a^2P$ to meniscus tensile strength $F_0=\pi\Gamma a$, for grains of diameter $a$ joined by menisci with surface tension $\Gamma$. We pay special attention to the quasi-static limit of slow flow and observe systematic, enduring strain localization in some of the cohesion-dominated ($P^*\sim 0.1$) systems. Homogeneous steady flows are characterized by the dependence of internal friction coefficient $\mu^*$ and solid fraction $\Phi$ on $I$ and $P^*$. We record fairly small but not negligible normal stress differences and the moderate sensitivity of the system to saturation within the pendular regime. Capillary forces have a significant effect on the macroscopic behavior of the system, up to $P^*$ values of several units. The concept of effective pressure may be used to predict an order of magnitude for the strong increase of $\mu^*$ as $P^*$ decreases but such a crude approach is unable to account for the complex structural changes induced by capillary cohesion. Likewise, the Mohr-Coulomb criterion for pressure-dependent critical states is, at best, an approximation valid within a restricted range of pressures, with $P^*\ge 1$. At small enough $P^*$, large clusters of interacting grains form in slow flows, in which liquid bonds survive shear strains of several units. This affects the anisotropies associated to different interactions, and the shape of function $\mu^*(I)$, which departs more slowly from its quasistatic limit than in cohesionless systems.