Using particle dynamics simulations, we investigate the strength and microstructure of agglomerates of wet frictional particles subjected to axial compression. The numerical model accounts for the cohesive and viscous effects of the binding liquid up to a debonding distance with the liquid assumed to be distributed homogeneously inside the agglomerate. We show that wet agglomerates undergo plastic deformation due to the rearrangements of primary particles during compression. The compressive strength is thus characterized by the plastic threshold before the onset of failure by the irreversible loss of wet contacts between primary particles. We find that the agglomerate plastic threshold is proportional to the characteristic cohesive stress defined from the liquid-vapor surface tension and the mean diameter of primary particles, with a prefactor that is a nearly linear function of the debonding distance and increases with size span. We analyze the agglomerate microstructure and, considering only the cohesive capillary forces at all bonds between primary particles, we propose an expression of the plastic strength as a function of the texture parameters such as the wet coordination number and packing fraction. This expression is shown to be consistent with our simulations up to a multiplicative factor reflecting the distribution of the capillary bridges. (C) 2018 Published by Elsevier Ltd.

VL - 92 UR - https://www-sciencedirect-com.libproxy.mit.edu/science/article/pii/S0093641318301216 JO - Mechanics Research Communications ER - TY - JOUR T1 - Effect of particle size distribution on 3D packings of spherical particles JF - EPJ Web of Conferences Y1 - 2017 A1 - Taiebat, Mahdi A1 - Patrick Mutabaruka A1 - Roland Jean-Marc Pellenq A1 - Farhang Radjaï ED - Saeid Nezamabadi ED - Luding, S. ED - Jean-Yves Delenne AB -We use molecular dynamics simulations of frictionless spherical particles to investigate a class of polydisperse granular materials in which the particle size distribution is uniform in particle volumes. The particles are assembled in a box by uniaxial compaction under the action of a constant stress. Due to the absence of friction and the nature of size distribution, the generated packings have the highest packing fraction at a given size span, defined as the ratio *α* of the largest size to the smallest size. We find that, up to *α* = 5, the packing fraction is a nearly linear function of *α*. While the coordination number is nearly constant due to the isostatic nature of the packings, we show that the connectivity of the particles evolves with *α*. In particular, the proportion of particles with 4 contacts represents the largest proportion of particles mostly of small size. We argue that this particular class of particles occurs as a result of the high stability of local configurations in which a small particle is stuck by four larger particles.

The shear behavior of granular materials immersed in a viscous fluid depends on fluid properties (viscosity, density), particle properties (size, density) and boundary conditions (shear rate, confining pressure). Using computational fluid dynamics simulations coupled with molecular dynamics for granular flow, and exploring a broad range of the values of parameters, we show that the parameter space can be reduced to a single parameter that controls the packing fraction and effective friction coefficient. This control parameter is a modified inertial number that incorporates viscous effects.

VL - 140 UR - http://www.epj-conferences.org/10.1051/epjconf/201714009044 JO - EPJ Web Conf. ER - TY - JOUR T1 - Strength of wet agglomerates of spherical particles: effects of friction and size distribution JF - EPJ Web of Conferences Y1 - 2017 A1 - Vo, Thanh-Trung A1 - Patrick Mutabaruka A1 - Jean-Yves Delenne A1 - Saeid Nezamabadi A1 - Farhang Radjaï ED - Luding, S. AB -We investigate the mechanical behavior of wet granular agglomerates composed of spherical particles by means of molecular dynamics simulations. The capillary cohesion force is modeled as an attraction force at the contact between two particles and expressed as an explicit function of the gap and volume of the liquid bridge. We are interested in the effect of the friction coefficient between primary particles. The agglomerates are subjected to diametrical compression tests. We find that the deformation is ductile involving particle rearrangements. However, a well-defined stress peak is observed and the peak stress is used as a measure of the compressive strength of the agglomerate. The strength increases with friction coefficient but levels off at friction coefficients above 0.4. Furthermore, the compressive strength is an increasing function of particle size span.

VL - 140 UR - http://www.epj-conferences.org/10.1051/epjconf/201714008021 JO - EPJ Web Conf. ER - TY - JOUR T1 - Viscoinertial regime of immersed granular flows JF - Physical Review E Y1 - 2017 A1 - L. Amarsid A1 - Jean-Yves Delenne A1 - Patrick Mutabaruka A1 - Yann Monerie A1 - Perales, F. A1 - Farhang Radjaï AB -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.

VL - 96 IS - 1 JO - Phys. Rev. E ER - TY - JOUR T1 - Transient dynamics of a 2D granular pile JF - The European Physical Journal E Y1 - 2015 A1 - Patrick Mutabaruka A1 - Krishna Kumar A1 - Kenichi Soga A1 - Farhang Radjaï A1 - Jean-Yves Delenne AB -We investigate by means of Contact Dynamics simulations the transient dynamics of a 2D granular pile set into motion by applying shear velocity during a short time interval to all particles. The spreading dynamics is directly controlled by the input energy whereas in recent studies of column collapse the dynamics scales with the initial potential energy of the column. As in column collapse, we observe a power-law dependence of the runout distance with respect to the input energy with nontrivial exponents. This suggests that the power-law behavior is a generic feature of granular dynamics, and the values of the exponents reflect the distribution of kinetic energy inside the material. We observe two regimes with different values of the exponents: the low-energy regime reflects the destabilization of the pile by the impact with a runout time independent of the input energy whereas the high-energy regime is governed by the input energy. We show that the evolution of the pile in the high-energy regime can be described by a characteristic decay time and the available energy after the pile is destabilized.

By means of coupled molecular dynamics-computational fluid dynamics simulations, we analyze the initiation of avalanches in a granular bed of spherical particles immersed in a viscous fluid and inclined above its angle of repose. In quantitative agreement with experiments, we find that the bed is unstable for a packing fraction below 0.59 but is stabilized above this packing fraction by negative excess pore pressure induced by the effect of dilatancy. From detailed numerical data, we explore the time evolution of shear strain, packing fraction, excess pore pressures, and granular microstructure in this creeplike pressure redistribution regime, and we show that they scale excellently with a characteristic time extracted from a model based on the balance of granular stresses in the presence of a negative excess pressure and its interplay with dilatancy. The cumulative shear strain at failure is found to be ≃ 0.2, in close agreement with the experiments, irrespective of the initial packing fraction and inclination angle. Remarkably, the avalanche is triggered when dilatancy vanishes instantly as a result of fluctuations while the average dilatancy is still positive (expanding bed) with a packing fraction that declines with the initial packing fraction. Another nontrivial feature of this creeplike regime is that, in contrast to dry granular materials, the internal friction angle of the bed at failure is independent of dilatancy but depends on the inclination angle, leading therefore to a nonlinear dependence of the excess pore pressure on the inclination angle. We show that this behavior may be described in terms of the contact network anisotropy, which increases with a nearly constant connectivity and levels off at a value (critical state) that increases with the inclination angle. These features suggest that the behavior of immersed granular materials is controlled not only directly by hydrodynamic forces acting on the particles but also by the influence of the fluid on the granular microstructure.

VL - 89 IS - 5 JO - Phys. Rev. E ER - TY - Generic T1 - Shear flow of dense granular suspensions by computer simulations T2 - 3rd International Symposium on Geomechanics from Micro to Macro Y1 - 2014 A1 - L. Amarsid A1 - Patrick Mutabaruka A1 - Jean-Yves Delenne ED - Kenichi Soga ED - Krishna Kumar ED - Giovanna Biscontin AB -We analyze the shear flow of dense granular materials composed of circular particles immersed in a viscous fluid by means of Molecular Dynamics simulations interfaced with the Lattice Boltzmann Method. A homogeneous flow of the suspension is obtained through periodic boundary conditions and by directly applying a confining pressure on the granular phase and shearing the fluid phase. The stead-state rheology can be described in terms of effective friction coefficient and packing fraction of the suspension as a function of the ratio of viscous shear stress to confining pressure (frictional description), on one hand, and in terms of normal and shear viscosities of the suspension as a function of the packing fraction (viscous description), on the other hand. We show that the simulation data are consistent with both descriptions and in close agreement with the corresponding scaling laws observed in recent experiments.

JF - 3rd International Symposium on Geomechanics from Micro to Macro PB - CRC Press CY - SEP 01-03-2014 Univ Cambridge, Cambridge, ENGLAND VL - Geomechanics from Micro to Macro UR - http://prodinra.inra.fr/record/369542 ER -