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.

It is well known that granular mixtures that differ in size or shape segregate when sheared. In the past, two mechanisms have been proposed to describe this effect, and it is unclear if both exist. To settle this question, we consider a bidisperse mixture of spheroids of equal volume in a rotating drum, where the two mechanisms are predicted to act in opposite directions. We present evidence that there are two distinct segregation mechanisms driven by relative overstress. Additionally, we showed that, for nonspherical particles, these two mechanisms (kinetic and gravity) can act in different directions leading to a competition between the effects of the two. As a result, the segregation intensity varies nonmonotonically as a function of aspect ratio (AR), and, at specific points, the segregation direction changes for both prolate and oblate spheroids, explaining the surprising segregation reversal previously reported. Consistent with previous results, we found that the kinetic mechanism is dominant for (almost) spherical particles. Furthermore, for moderate aspect ratios, the kinetic mechanism is responsible for the spherical particles' segregation to the periphery of the drum, and the gravity mechanism plays only a minor role. Whereas, at the extreme values of AR, the gravity mechanism notably increases and overtakes its kinetic counterpart.

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.