The discharge of elongated particles from a silo with rotating bottom is investigated numerically. The introduction of a slight transverse shear reduces the flow rate Q by up to 70% compared with stationary bottom, but the flow rate shows a modest increase by further increasing the external shear. Focusing on the dependency of flow rate Q on orifice diameter D, the spheres and rods show two distinct trends. For rods, in the small-aperture limit Q seems to follow an exponential trend, deviating from the classical power-law dependence. These macroscopic observations are in good agreement with our earlier experimental findings [Phys. Rev. E 103, 062905 (2021)]. With the help of the coarse-graining methodology we obtain the spatial distribution of the macroscopic density, velocity, kinetic pressure, and orientation fields. This allows us detecting a transition from funnel to mass flow pattern caused by the external shear. Additionally, averaging these fields in the region of the orifice reveals that the strong initial decrease in Q is mostly attributed to changes in the flow velocity, while the weakly increasing trend at higher rotation rates is related to increasing packing fraction. Similar analysis of the grain orientation at the orifice suggests a correlation of the flow rate magnitude with the vertical orientation and the packing fraction at the orifice with the order of the grains. Lastly, the vertical profile of mean acceleration at the center of the silo denotes that the region where the acceleration is not negligible shrinks significantly due to the strong perturbation induced by the moving wall.

In the present work, we investigate experimentally and numerically the motion of solid macroscopic spheres (Brownian and colloidal effects are negligible) when settling from rest in a quiescent fluid toward a solid wall under confined and unconfined configurations. Particle trajectories for spheres of two types of materials are measured using a high-speed digital camera. For unconfined configurations, our experimental findings are in excellent agreement with well-established analytical frameworks, used to describe the forces acting on the sphere. Besides, the experimental values of the terminal velocity obtained for different confinements are also in very good agreement with previous theoretical formulations. Similar conditions are simulated using a resolved CFD-DEM approach. After adjusting the parameters of the numerical model, we analyze the particle dynamics under several confinement conditions. The simulations results are contrasted with the experimental findings, obtaining a good agreement. We analyze several systems varying the radius of the bead and show the excellent agreement of our results with previous analytical approaches. However, the results indicate that confined particles have a distinct dynamics response when approaching the wall. Consequently, their motion cannot be described by the analytical framework introduced for the infinite system. Indeed, the confinement strongly affects the spatial scale where the particle is affected by the bottom wall and, accordingly, the dimensionless results cannot be collapsed in a single master curve, using the particle size as a characteristic length. Alternatively, we rationalize our findings using a kinematic approximation to highlight the relevant scale of the problem. Our outcomes suggest it is possible to determine a new spatial scale to describe the collisional process, depending on the specific confining conditions.

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 emergence of coherent vortices has been observed in a wide variety of many-body systems such as animal flocks, bacteria, colloids, vibrated granular materials or human crowds. Here, we experimentally demonstrate that pedestrians roaming within an enclosure also form vortex-like patterns which, intriguingly, only rotate counterclockwise. By implementing simple numerical simulations, we evidence that the development of swirls in many-particle systems can be described as a phase transition in which both the density of agents and their dissipative interactions with the boundaries play a determinant role. Also, for the specific case of pedestrians, we show that the preference of right-handed people (the majority in our experiments) to turn leftwards when facing a wall is the symmetry breaking mechanism needed to trigger the observed global counterclockwise rotation.

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.

Granular gases are interesting multiparticle systems which, irrespective of the apparent simplicity of particle interactions, exhibit a rich scenario of so far only little understood features. We have numerically investigated a dense granular gas composed of frictional spherocylinders which are excited mechanically by lateral vibrating container walls. This study was stimulated by experiments in microgravity on parabolic flights. The formation of spatial inhomogeneities (clusters) was observed in a region near the corners of the container, about halfway from the excitation plates. The particles in the clusters show a tendency to align parallel to the container walls, seemingly increasing the stabilizing effect of friction. The simulation results provide hints that the phase difference of the vibrations of the two excitation walls might affect the cluster dynamics.

Some years ago, Harth et al. experimentally explored the steady state dynamics of a heated granular gas of rod-like particles in microgravity [K. Harth et al. Phys. Rev. Lett. 110, 144102 (2013)]. Here, we report numerical results that quantitatively reproduce their experimental findings and provide additional insight into the process. A system of sphero-cylinders is heated by the vibration of three flat side walls, resulting in one symmetrically heated direction, one non-symmetrically heated direction, and one non-heated direction.

We present a critical comparative analysis between numerical and experimental results of quasi-two-dimensional silo and hopper flows. In our approach, the Discrete Element Method was employed to describe a single-layer mono-disperse sphere confined by two parallel walls with an orifice at the bottom. As a first step, we examined the discharge process, varying the size of the outlet and the hopper angle. Next, we set the simulation parameters fitting the experimental flow rate values obtained experimentally. Remarkably, the numerical model captured the slight non-monotonic dependence of the flow rate with the hopper angle, which was detected experimentally. Additionally, we analyzed the vertical velocity and solid fractions profiles at the outlet numerically and experimentally. Although numerical results also agreed with the experimental observations, a slight deviation appeared systematically between both approaches. Finally, we explored the impact of the system's confinement on this process, examining the consequences of particle-particle and particle-wall friction on the system macroscopic response. We mainly found that the degree of confinement and particle-wall friction have a relevant impact on the outflow dynamics. Our analysis demonstrated that the naive 2D approximation of this 3D flow process fails to describe it accurately.

We report experimental results of the speed-density relation emerging in pedestrian dynamics when individuals keep a prescribed safety distance among them. To this end, we characterize the movement of a group of people roaming inside an enclosure varying different experimental parameters: (i) global density, (ii) prescribed walking speed, and (iii) suggested safety distance. Then, by means of the Voronoi diagram we are able to compute the local density associated to each pedestrian, which is afterward correlated with its corresponding velocity at each time. In this way, we discover a strong dependence of the speed-density relation on the experimental conditions, especially with the (prescribed) free speed. We also observe that when pedestrians walk slowly, the speed-density relation depends on the global macroscopic density of the system, and not only on the local one. Finally, we demonstrate that for the same experiment, each pedestrian follows a distinct behavior, thus giving rise to multiple speed-density curves.

In this work, we numerically examine the steady-state granular flow of 3D non-spherical particles down an inclined plane. We use a hybrid CPU/GPU implementation of the discrete element method of nonspherical elongated particles. Thus, a systematic study of the system response is performed varying the particle aspect ratio and the plane inclination. Similar to the case of spheres, we observe three well-defined regimes: arresting flows, steady uniform flows and accelerating flows. Both steady and dynamic macroscopic macroscopic fields are derived from microscopic data, by time-averaging and spatial smoothing (coarse-graining), including density, velocity, as well as the kinetic and contact stress tensors. The internal morphology of the flow was quantified by exploring the solid fraction profiles and the particle orientation distribution. Furthermore, the system¿s characteristic time and length scales are investigated in detail. Our aim is to achieve a continuum mechanical description of granular flows composed of non-spherical particles based on the micromechanical details. Thus, to evaluate the influence of particle shape on the constitutive response in granular of those systems. However, to meet the proceeding¿s page restrictions here we will only discuss the dependence of some terms of the continuum averaged equations on the coarse-graining scale, specifically the case of the kinetic part of the stress tensor.

The effect of the filling mechanism on the packing of faceted particles with different aspect ratios has been examined. We have experimentally measured the particle angular distribution and the packing fraction of ensembles of faceted particles deposited in a two-dimensional box. The granular system has been numerically simulated using a two-dimensional model of faceted particles. We found that increasing the feeding rate results in an enhancement of the disorder in the final deposit and, consequently, in a reduction of the number of particles oriented in their most stable configuration. In this regime, the final packing fraction monotonically decreases as the feeding rate increases. The correlations between the final packing morphology and the stress transmission were examined by describing the micromechanical properties of the deposits. For the case of elongated particles, increasing the feeding rate leads to an enhancement of the stress transmission towards the sides of the box. On the contrary, for the case of square particles, increasing the feeding rate promotes vertical transmission of the stress.

General-purpose computation on Graphics Processing Units (GPU) on staff computers has recently become an attractive alternative to parallel computing on clusters and supercomputers. We present the GPU-implementation of an accurate molecular dynamics algorithm for a system of spheres. accurate molecular dynamics algorithm for a system of spheres. The new hybrid CPU-GPU implementation takes into account all the degrees of freedom, including the quaternion representation of 3D rotations. For additional versatility, the contact interaction between particles is defined using a force law of enhanced generality, which accounts for the elastic and dissipative interactions, and the hard-sphere interaction parameters are translated to the soft-sphere parameter set. We prove that the algorithm complies with the statistical mechanical laws by examining the homogeneous cooling of a granular gas with rotation. The results are in excellent agreement with well established mean-field theories for low-density hard sphere systems. This GPU technique dramatically reduces user waiting time, compared with a traditional CPU implementation.

We use particle-based simulations to examine the flow of particles through an exit. Simulations involve both gravity-driven particles (representing granular material) and velocity-driven particles (mimicking pedestrian dynamics). Contact forces between particles include elastic, viscous, and frictional forces; and simulations use bunker geometry. Power laws are observed in the relation between flow rate and exit width. Simulations of granular flow showed that the power law has little dependence on the coefficient of friction. Polydisperse granular systems produced higher flow rates than those produced by monodisperse ones. We extend the particle model to include the main features of pedestrian dynamics: thoracic shape, shoulder rotation, and desired velocity oriented towards the exit. Higher desired velocity resulted in higher flow rate. Granular simulations always give higher flow rate than pedestrian simulations, despite the values of aspect ratio of the particles. In terms of force distribution, pedestrians and granulates share similar properties with the non-democratic distribution of particles. with the non-democratic distribution of forces that poses high risks of injuries in a bottleneck situation.

The role of density and velocity profiles in the flow of particles through apertures has been recently put on evidence in a two-dimensional experiment (Phys. Rev. Let. 108, 108, 118001). in a two-dimensional experiment (Phys. Rev. Lett. 108, 248001). For the whole range of apertures studied, both velocity and density profiles are selfsimilar and the obtained scaling functions allow to derive the relevant scales of the problem. Indeed, by means of the functionality obtained for these profiles, an exact expression for the mass flow rate was proposed. Such expression showed a perfect agreement with the experiential data. In this work, we generalize this study to the three dimensional case.We perform numerical simulations of a 3D silo in which the velocity and volume fraction profiles are determined. Both profiles shows that the scaling obtained for 2D can be generalized to the 3D case. Finally, the scaling of the mass flow rate with the outlet radius is discussed. outlet radius is discussed.