Stress partition and microstructure in size-segregating granular flows.

When a granular mixture involving grains of different sizes is shaken, sheared, mixed, or left to flow, grains tend to separate by sizes in a process known as size segregation. In this study, we explore the size segregation mechanism in granular chute flows in terms of the pressure distribution and granular microstructure. Therefore, two-dimensional discrete numerical simulations of bidisperse granular chute flows are systematically analyzed. Based on the theoretical models of J. M. N. T. Gray and A. R. Thornton [Proc. R. Soc. A 461, 1447] and K. M. Hill and D. S. Tan [J. Fluid Mech. 756, 54 (2014)], we explore the stress partition in the phases of small and large grains, discriminating between contact stresses and kinetic stresses. Our results support both gravity-induced and shear-gradient-induced segregation mechanisms. However, we show that the contact stress partition is extremely sensitive to the definition of the partial stress tensors and, more specifically, to the way mixed contacts (i.e., involving a small grain and a large grain) are handled, making conclusions on gravity-induced segregation uncertain. By contrast, the computation of the partial kinetic stress tensors is robust. The kinetic pressure partition exhibits a deviation from continuum mixture theory of a significantly higher amplitude than the contact pressure and displays a clear dependence on the flow dynamics. Finally, using a simple approximation for the contact partial stress tensors, we investigate how the contact stress partition relates to the flow microstructure and suggest that the latter may provide an interesting proxy for studying gravity-induced segregation.

Continuum simulation of the discharge of the granular silo: a validation test for the µ(I) visco-plastic flow law.

Using a continuum Navier-Stokes solver with the µ(I) flow law implemented to model the viscous behavior, and the discrete Contact Dynamics algorithm, the discharge of granular silos is simulated in two dimensions from the early stages of the discharge until complete release of the material. In both cases, the Beverloo scaling is recovered. We first do not attempt a quantitative comparison, but focus on the qualitative behavior of velocity and pressure at different locations in the flow. A good agreement for the velocity is obtained in the regions of rapid flows, while areas of slow creep are not entirely captured by the continuum model. The pressure field shows a general good agreement, while bulk deformations are found to be similar in both approaches. The influence of the parameters of the µ(I) flow law is systematically investigated, showing the importance of the dependence on the inertial number I to achieve quantitative agreement between continuum and discrete discharge. However, potential problems involving the systems size, the configuration and "non-local" effects, are suggested. Yet the general ability of the continuum model to reproduce qualitatively the granular behavior is found to be very encouraging.

Friction and the oscillatory motion of granular flows.

This contribution reports on numerical simulations of two-dimensional granular flows on erodible beds. The broad aim is to investigate whether simple flows of model granular matter exhibit spontaneous oscillatory motion in generic flow conditions, and in this case, whether the frictional properties of the contacts between grains may affect the existence or the characteristics of this oscillatory motion. The analysis of different series of simulations shows that the flow develops an oscillatory motion with a well-defined frequency which increases like the inverse of the velocity's square root. We show that the oscillation is essentially a surface phenomenon. The amplitude of the oscillation is higher for lower volume fractions and can thus be related to the flow velocity and grains' friction properties. The study of the influence of the periodic geometry of the simulation cell shows no significant effect. These results are discussed in relation to sonic sands.

Memory of the unjamming transition during cyclic tiltings of a granular pile.

Discrete numerical simulations are performed to study the evolution of the microstructure and the response of a granular packing during successive loading-unloading cycles, consisting of quasistatic rotations in the gravity field between opposite inclination angles. We show that internal variables--e.g., stress and fabric of the pile--exhibit hysteresis during these cycles due to the exploration of different metastable configurations. Interestingly, the hysteretic behavior of the pile strongly depends on the maximal inclination of the cycles, giving evidence of the irreversible modifications of the pile state occurring close to the unjamming transition. More specifically, we show that for cycles with maximal inclination larger than the repose angle, the weak-contact network carries the memory of the unjamming transition. These results demonstrate the relevance of a two-phase description--strong- and weak-contact networks--for a granular system, as soon as it has approached the unjamming transition.

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.