Viscoinertial regime of immersed granular flows.

By means of extensive coupled molecular dynamics-lattice Boltzmann simulations, accounting for grain dynamics and subparticle resolution of the fluid phase, we analyze steady inertial granular flows sheared by a viscous fluid. We show that, for a broad range of system parameters (shear rate, confining stress, fluid viscosity, and relative fluid-grain density), the frictional strength and packing fraction can be described by a modified inertial number incorporating the fluid effect. In a dual viscous description, the effective viscosity diverges as the inverse square of the difference between the packing fraction and its jamming value, as observed in experiments. We also find that the fabric and force anisotropies extracted from the contact network are well described by the modified inertial number, thus providing clear evidence for the role of these key structural parameters in dense suspensions.

Collapse dynamics and runout of dense granular materials in a fluid.

We investigate the effect of an ambient fluid on the dynamics of collapse and spread of a granular column simulated by means of the contact dynamics method interfaced with computational fluid dynamics. The runout distance is found to increase as a power law with the aspect ratio of the column, and, surprisingly, for a given aspect ratio and packing fraction, it may be similar in the grain-inertial and fluid-inertial regimes but with considerably longer duration in the latter case. We show that the effect of fluid in viscous and fluid-inertial regimes is to both reduce the kinetic energy during collapse and enhance the flow by lubrication during spread. Hence, the runout distance in a fluid may be below or equal to that in the absence of fluid due to compensation between those effects.

Tensile strength and fracture of cemented granular aggregates.

Cemented granular aggregates include a broad class of geomaterials such as sedimentary rocks and some biomaterials such as the wheat endosperm. We present a 3D lattice element method for the simulation of such materials, modeled as a jammed assembly of particles bound together by a matrix partially filling the interstitial space. From extensive simulation data, we analyze the mechanical properties of aggregates subjected to tensile loading as a function of matrix volume fraction and particle-matrix adhesion. We observe a linear elastic behavior followed by a brutal failure along a fracture surface. The effective stiffness before failure increases almost linearly with the matrix volume fraction. We show that the tensile strength of the aggregates increases with both the increasing tensile strength at the particle-matrix interface and decreasing stress concentration as a function of matrix volume fraction. The proportion of broken bonds in the particle phase reveals a range of values of the particle-matrix adhesion and matrix volume fraction for which the cracks bypass the particles and hence no particle damage occurs. This limit is shown to depend on the relative toughness of the particle-matrix interface with respect to the particles.