This work is devoted to study numerically the self-diffusion of spherocylindrical particles flowing down an inclined plane, using the discrete element method (DEM). This system is challenging due to particles being non-spherical and because they are subjected to a non-uniform shear rate. We performed simulations for several aspect ratios and inclination angles, tracking individual particle trajectories. Using the simulation data, we computed the diffusion coefficients D, and a coarse-graining methodology allowed accessing the shear rate spatial profiles (y)over dot (z). This data enabled us to identify the spatial regions where the diffusivity strongly correlates with the local shear rate. Introducing an effective particle size d(perpendicular to), we proposed a well-rationalized scaling law between D and (y)over dot. Our findings also identified specific locations where the diffusivity does not correlate with the shear rate. This observation corresponds to zones where has non-linear spatial variation, and the velocity probability density distributions exhibit asymmetric shapes.

The time evolution of silo discharge is investigated for different granular materials made of spherical or elongated grains in laboratory experiments and with discrete element model (DEM) calculations. For spherical grains, we confirm the widely known typical behavior with constant discharge rate (except for initial and final transients). For elongated particles with aspect ratios between 2 < L/d < 6.1, we find a peculiar flow rate increase for larger orifices before the end of the discharge process. While the flow field is practically homogeneous for spherical grains, it has strong gradients for elongated particles with a fast-flowing region in the middle of the silo surrounded by a stagnant zone. For large enough orifice sizes, the flow rate increase is connected with a suppression of the stagnant zone, resulting in an increase in both the packing fraction and flow velocity near the silo outlet within a certain parameter range.

We report experimental evidence of clogging due to the spontaneous development of hanging arches when a granular sample composed of spherical particles flows down a narrow vertical pipe. These arches, akin to the ones responsible for silo clogging, can only be possible due to the role of frictional forces; otherwise they will be unstable. We find that, contrary to the silo case, the probability of clogging in vertical narrow tubes does not decrease monotonically with the ratio of the pipe-to-particle diameters. This behavior is related to the clogging prevention caused by the spontaneous ordering of particles apparent in certain aspect ratios. More importantly, by means of numerical simulations, we discover that the interparticle normal force distributions broaden in systems with higher probability of clogging. This feature, which has been proposed before as a distinctive feature of jamming in sheared granular samples, suggests that clogging and jamming are connected in pipe flow.

Very recently, To et al. have experimentally explored granular flow in a cylindrical silo, with a bottom wall that rotates horizontally with respect to the lateral wall [Phys. Rev. E 100, 012906 (2019)]. Here we numerically reproduce their experimental findings, in particular, the peculiar behavior of the mass flow rate Q as a function of the frequency of rotation f. Namely, we find that for small outlet diameters D the flow rate increased with f, while for larger D a nonmonotonic behavior is confirmed. Furthermore, using a coarse-graining technique, we compute the macroscopic density, momentum, and the stress tensor fields. These results show conclusively that changes in the discharge process are directly related to changes in the flow pattern from funnel flow to mass flow. Moreover, by decomposing the mass flux (linear momentum field) at the orifice into two main factors, macroscopic velocity and density fields, we obtain that the nonmonotonic behavior of the linear momentum is caused by density changes rather than by changes in the macroscopic velocity. In addition, by analyzing the spatial distribution of the kinetic stress, we find that for small orifices increasing rotational shear enhances the mean kinetic pressure < p(k)> and the system dilatancy. This reduces the stability of the arches, and, consequently, the volumetric flow rate increases monotonically. For large orifices, however, we detected that < p(k)> changes nonmonotonically, which might explain the ...

Projectile impact into a light granular material composed of expanded polypropylene (EPP) particles is investigated systematically with various impact velocities. Experimentally, the trajectory of an intruder moving inside the granular material is monitored with a recently developed non-invasive microwave radar system. Numerically, discrete element simulations together with coarse-graining techniques are employed to address both dynamics of the intruder and response of the granular bed. Our experimental and numerical results of the intruder dynamics agree with each other quantitatively and are in congruent with existing phenomenological model on granular drag. Stepping further, we explore the 'microscopic' origin of granular drag through characterizing the response of granular bed, including density, velocity and kinetic stress fields at the mean-field level. In addition, we find that the dynamics of cavity collapse behind the intruder changes significantly when increasing the initial speed . Moreover, the kinetic pressure ahead of the intruder decays exponentially in the co-moving system of the intruder. Its scaling gives rise to a characteristic length scale, which is in the order of intruder size. This finding is in perfect agreement with the long-scale inertial dissipation type that we find in all cases.