Stress transmission in cemented bidisperse granular materials.

We analyze stress distributions in a two-dimensional bidisperse cemented granular packing for a broad range of the values of particle-size ratio, the volumes of large and small particles, and the amount of cementing matrix. In such textured porous materials, the stress concentration, which controls the fracture and fragmentation of the material under tensile loading or in grinding processes, reflects not only the porosity but also the contact network of the particle phase and the resulting stress chains. By means of peridynamic simulations under tensile loading, we show how both the texture and stress distribution depend on size ratio, volume ratio, and the amount of the cementing matrix. In particular, the volume fraction of the class of small particles plays a key role in homogenizing stresses across the system by reducing porosity. Interestingly, the texture controls not only the porosity but also the distribution of pores inside the system with its statistical variability, found to be strongly correlated with the homogeneity of stresses inside the large particles. The most homogeneous stress distribution occurs for the largest size ratio and largest volume fraction of small particles, corresponding to the lowest pore size dispersion and the cushioning effect of small particles and its similar role to the binding matrix for stress redistribution across the packing. At higher porosity, the tensile stresses above the mean stress fall off exponentially in all phases with an exponent that strongly depends on the texture. The exponential part broadens with decreasing matrix volume fraction and particle-size ratio. These correlations reveal the strong interplay between size polydispersity and the cohesive action of the binding matrix for stress distribution, which is significant for the behavior of textured materials in grinding operations.

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.

Rheology of granular materials composed of nonconvex particles.

By means of contact dynamics simulations, we investigate the shear strength and internal structure of granular materials composed of two-dimensional nonconvex aggregates. We find that the packing fraction first grows as the nonconvexity is increased but declines at higher nonconvexity. This unmonotonic dependence reflects the competing effects of pore size reduction between convex borders of aggregates and gain in porosity at the nonconvex borders that are captured in a simple model fitting nicely the simulation data both in the isotropic and sheared packings. On the other hand, the internal angle of friction increases linearly with nonconvexity and saturates to a value independent of nonconvexity. We show that fabric anisotropy, force anisotropy, and friction mobilization, all enhanced by multiple contacts between aggregates, govern the observed increase of shear strength and its saturation with increasing nonconvexity. The main effect of interlocking is to dislocate frictional dissipation from the locked double and triple contacts between aggregates to the simple contacts between clusters of aggregates. This self-organization of particle motions allows the packing to keep a constant shear strength at high nonconvexity.

Multiscale force networks in highly polydisperse granular media.

We investigate highly polydisperse packings subjected to simple shear by contact dynamics simulations. A major unsolved issue is how granular texture and force chains depend on the size polydispersity and how far they influence the shear strength. The numerical treatment was made possible by ensuring the statistical representation of particle size classes. An unexpected finding is that the internal friction angle is independent of polydispersity. We show that this behavior is related to two mechanisms underlying the stability of force chains: (i) The class of largest particles captures strong force chains, and (ii) these chains are equilibrated by weak forces carried by increasingly smaller particles as the size span broadens. In the presence of adhesion between particles, the Coulomb cohesion increases with size polydispersity as a result of enhanced force anisotropy.

Short-time dynamics of a packing of polyhedral grains under horizontal vibrations.

We analyze the dynamics of a 3D granular packing composed of particles of irregular polyhedral shape confined inside a rectangular box with a retaining wall subjected to horizontal harmonic forcing. The simulations are performed by means of the contact dynamics method for a broad set of loading parameters. We explore the vibrational dynamics of the packing, the evolution of solid fraction and the scaling of dynamics with the loading parameters. We show that the motion of the retaining wall is strongly anharmonic as a result of jamming and grain rearrangements. It is found that the mean particle displacement scales with inverse square of frequency, the inverse of the force amplitude and the square of gravity. The short-time compaction rate grows in proportion to frequency up to a characteristic frequency, corresponding to collective particle rearrangements between equilibrium states, and then it declines in inverse proportion to frequency.

Space-filling properties of polydisperse granular media.

We present a systematic investigation of the morphology and space-filling properties of polydisperse densely packed granular media in two dimensions. A numerical procedure is introduced to generate collections of circular particles with size distributions of variable shape and span constrained by explicit criteria of statistical representativity. We characterize the domain of statistically accessible distribution parameters for a bounded number of particles. This particle generation procedure is used with two different deposition protocols in order to build large close-packed samples of prescribed polydispersity. We find that the solid fraction is a strongly nonlinear function of the size span, and the highest levels of solid fraction occur for the uniform distribution by volume fractions. As the span is increased, a transition occurs from a regime of topological disorder where the packing properties are governed by particle connectivity to a regime of metric disorder where pore-filling small particles prevail. The polydispersity manifests itself in the first regime through the variability of local coordination numbers. We observe a continuous decrease of the number of particles with four contacts and the growth of two populations of particles with three and five contacts. In the second regime, radial distribution functions show that the material is homogeneous beyond only a few average particle diameters. We also show that the packing orientational order is linked with fabric anisotropy and it declines with size span.

Strength and failure of cemented granular matter.

Cemented granular materials (CGMs) consist of densely packed solid particles and a pore-filling solid matrix sticking to the particles. We use a sub-particle lattice discretization method to investigate the particle-scale origins of strength and failure properties of CGMs. We show that jamming of the particles leads to highly inhomogeneous stress fields. The stress probability density functions are increasingly wider for a decreasing matrix volume fraction, the stresses being more and more concentrated in the interparticle contact zones with an exponential distribution as in cohesionless granular media. Under uniaxial loading, pronounced asymmetry can occur between tension and compression both in strength and in the initial stiffness as a result of the presence of bare contacts (with no matrix interposed) between the particles. Damage growth is analyzed by considering the evolution of stiffness degradation and the number of broken bonds in the particle phase. A brutal degradation appears in tension as a consequence of brittle fracture in contrast to the more progressive nature of damage growth in compression. We also carry out a detailed parametric study in order to assess the combined influence of the matrix volume fraction and particle-matrix adherence. Three regimes of crack propagation can be distinguished corresponding to no particle damage, particle abrasion and particle fragmentation, respectively. We find that particle damage scales well with the relative toughness of the particle-matrix interface with respect to the particle toughness. This relative toughness is a function of both matrix volume fraction and particle-matrix adherence and it appears therefore to be the unique control parameter governing transition from soft to hard behavior.

Stress transmission in wet granular materials.

We analyze stress transmission in wet granular media in the pendular state by means of three-dimensional molecular-dynamics simulations. We show that the tensile action of capillary bonds induces a self-stressed particle network organized in two percolating "phases" of positive and negative particle pressures. Various statistical descriptors of the microstructure and bond force network are used to characterize this partition. Two basic properties emerge: 1) the highest particle pressure is located in the bulk of each phase; 2) the lowest pressure level occurs at the interface between the two phases, involving also the largest connectivity of the particles via tensile and compressive bonds. When a confining pressure is applied, the number of tensile bonds falls off and the negative phase breaks into aggregates and isolated sites.

Multi-scale analysis of the stress state in a granular slope in transition to failure.

By means of contact dynamics simulations, we analyze the stress state in a granular bed slowly tilted toward its angle of repose. An increasingly large number of grains are overloaded in the sense that they are found to carry a stress ratio above the Coulomb yield threshold of the whole packing. Using this property, we introduce a coarse-graining length scale at which all stress ratios are below the packing yield threshold. We show that this length increases with the slope angle and jumps to a length comparable to the depth of the granular bed at an angle below the angle of repose. This transition coincides with the onset of dilation in the packing. We map this transition into a percolation transition of the overloaded grains, and discuss it in terms of long-range correlations and granular slope metastability.

Self-stresses and crack formation by particle swelling in cohesive granular media.

We present a molecular-dynamics study of force patterns, tensile strength, and crack formation in a cohesive granular model where the particles are subjected to swelling or shrinkage gradients. Nonuniform particle size change generates self-equilibrated forces that lead to crack initiation as soon as the strongest tensile contacts begin to fail. We find that the tensile strength is well below the theoretical strength as a result of inhomogeneous force transmission in granular media. The cracks propagate either inward from the edge upon shrinkage or outward from the center upon swelling. We show that the coarse-grained stresses are correctly predicted by an elastic model that incorporates particle size change as metric evolution.

Intermittent flow of a collection of rigid frictional disks in a vertical pipe.

We study the quasistatic flow of a collection of rigid frictional disks pushed upward (against the gravity) inside a narrow vertical pipe by a compliant mechanism. The contact dynamics method was used for the numerical simulations in combination with a friction law at disk-disk and wall-disk contacts characterized by discontinuous velocity weakening from a static threshold to a dynamic coefficient of friction. The material is sheared by the rolling of particles at the walls inducing a convective motion in the bulk. We observe a transition from constant flow to an intermittent flow when the driving velocity is reduced below a characteristic velocity that scales as k(-1/2) with the stiffness k of the pushing mechanism. The intermittent flow is composed of alternating phases of creep motion, where the pressure at the bottom of the granular column rises nonlinearly with time, and sudden slip, corresponding to a fast pressure drop. We show that the mean static pressure is correctly predicted by the Janssen model. The interplay between friction mobilization at the walls and structural changes in the bulk gives rise to a broad distribution of slip amplitudes characterized by a power law with an exponent approximately -1.7 that appears to be robust with respect to our system parameters.

Model for granular texture with steric exclusion.

We propose a method to characterize the geometrical texture of a granular packing at the particle scale including the steric hindrance effect. This method is based on the assumption of a maximum directional disorder (statistical entropy) compatible with both the strain-induced anisotropy of the contact network and steric exclusions. We show that the predicted statistics for the local configurations are in fairly good agreement with our numerical data.

Dynamic rearrangements and packing regimes in randomly deposited two-dimensional granular beds.

We study the structural properties of two-dimensional granular packings prepared by random deposition from a source line. We consider a class of random ballistic deposition models based on single-particle relaxation rules controlled by a critical angle, and we show that these local rules can be formulated as rolling friction in the framework of dynamic methods for the simulation of granular materials. We find that a packing prepared by random deposition models is generically unstable, and undergoes dynamic rearrangements. As a result, the dynamic method leads systematically to a higher solid fraction than the geometrical model for the same critical angle. We characterize the structure of the packings generated by both methods in terms of solid fraction, contact connectivity, and anisotropy. Our analysis provides evidence for four packing regimes as a function of solid fraction, the mechanisms of packing growth being different in each regime.