Role of particle size in the kinematic properties of silo flow

We experimentally analyze the effect that particle size has on the mass flow rate of a quasi two-dimensional silo discharged by gravity. In a previous work, Janda et al. [Phys. Rev. Lett. 108, 248001 (2012)PRLTAO0031-900710.1103/PhysRevLett.108.248001] introduced a new expression for the mass flow rate based on a detailed experimental analysis of the flow for 1-mm diameter beads. Here, we aim to extend these results by using particles of larger sizes and a variable that was not explicitly included in the proposed expression. We show that the velocity and density profiles at the outlet are self-similar and scale with the outlet size with the same functionalities as in the case of 1-mm particles. Nevertheless, some discrepancies are evidenced in the values of the fitting parameters. In particular, we observe that larger particles lead to higher velocities and lower packing fractions at the orifice. Intriguingly, both magnitudes seem to compensate giving rise to very similar flow rates. In order to shed light on the origin of this behavior we have computed fields of a solid fraction, velocity, and a kinetic-stress like variable in the region above the orifice.

Simulating competitive egress of noncircular pedestrians

We present a numerical framework to simulate pedestrian dynamics in highly competitive conditions by means of a force-based model implemented with spherocylindrical particles instead of the traditional, symmetric disks. This modification of the individuals' shape allows one to naturally reproduce recent experimental findings of room evacuations through narrow doors in situations where the contact pressure among the pedestrians was rather large. In particular, we obtain a power-law tail distribution of the time lapses between the passage of consecutive individuals. In addition, we show that this improvement leads to new features where the particles' rotation acquires great significance.

Effect of obstacle position in the flow of sheep through a narrow door

In a recent work [Phys. Rev. E 91, 022808 (2015)PLEEE81539-375510.1103/PhysRevE.91.022808] it was reported that placing an obstacle in front of a gate has a beneficial effect in the flow of sheep through it. Here, we extend such results by implementing three different obstacle positions. We have observed that the flow is improved in two cases, while it worsens in the other one; the last instance happens when the obstacle is too close to the door. In this situation, the outcomes suggest that clogging develops between the doorjamb and the obstacle, contrary to the cases when the obstacle is farther, in which case clogging always occurs at the very door. The effectiveness of the obstacle (a strategy put forward to alleviate clogging in emergency exits) is therefore quite sensitive to its location. In addition, the study of the temporal evolution of the flow rate as the test develops makes evident a steady behavior during the entire duration of the entrance. This result is at odds with recent findings in human evacuation tests where the flow rate varies over time, therefore challenging the fairness of straightforward comparisons between pedestrian behavior and animal experimental observations.

Clogging of granular materials in silos: effect of gravity and outlet size

By means of extensive numerical simulations we disclose the role of the driving force in the clogging of inert particles passing through a constriction. We uncover the effect of gravity and outlet size on the flow rate and kinetic energy within the system, and use these quantities to deepen our understanding of the blocking process. First, we confirm the existence of a finite avalanche size when the driving force tends to zero. The magnitude of this limit avalanche size grows with the outlet size, as expected due to geometrical reasons. In addition, there is an augment of the avalanche size when the driving force is increased, an effect that is enhanced by the outlet size. This phenomenology is explained by assuming that in order to get a stable clog developed, two conditions must be fulfilled: (1) an arch spanning the outlet size should be formed; (2) the arch should resist until the complete dissipation of the kinetic energy within the system. From these assumptions, we are able to obtain the probability that an arch gets destabilized, which is shown to primarily depend on the square root of the kinetic energy. A minor additional dependence of the outlet size is also observed which is explained in the light of recent results of the arch resistance in vibrated silos.

Experimental proof of faster-is-slower in systems of frictional particles flowing through constrictions

The "faster-is-slower" (FIS) effect was first predicted by computer simulations of the egress of pedestrians through a narrow exit [D. Helbing, I. J. Farkas, and T. Vicsek, Nature (London) 407, 487 (2000)]. FIS refers to the finding that, under certain conditions, an excess of the individuals' vigor in the attempt to exit causes a decrease in the flow rate. In general, this effect is identified by the appearance of a minimum when plotting the total evacuation time of a crowd as a function of the pedestrian desired velocity. Here, we experimentally show that the FIS effect indeed occurs in three different systems of discrete particles flowing through a constriction: (a) humans evacuating a room, (b) a herd of sheep entering a barn, and (c) grains flowing out a 2D hopper over a vibrated incline. This finding suggests that FIS is a universal phenomenon for active matter passing through a narrowing.

Experimental proof of faster-is-slower in systems of frictional particles flowing through constrictions

The "faster-is-slower" (FIS) effect was first predicted by computer simulations of the egress of pedestrians through a narrow exit [D. Helbing, I. J. Farkas, and T. Vicsek, Nature (London) 407, 487 (2000)]. FIS refers to the finding that, under certain conditions, an excess of the individuals' vigor in the attempt to exit causes a decrease in the flow rate. In general, this effect is identified by the appearance of a minimum when plotting the total evacuation time of a crowd as a function of the pedestrian desired velocity. Here, we experimentally show that the FIS effect indeed occurs in three different systems of discrete particles flowing through a constriction: (a) humans evacuating a room, (b) a herd of sheep entering a barn, and (c) grains flowing out a 2D hopper over a vibrated incline. This finding suggests that FIS is a universal phenomenon for active matter passing through a narrowing.

Stability of clogging arches in a silo submitted to vertical vibrations

We present experimental results on the endurance of arches that block the outlet of a two-dimensional silo when subjected to vertical vibration. In a recent paper [C. Lozano et al., Phys. Rev. Lett. 109, 068001 (2012)], it was shown that the arch resistance against vibrations is determined by the maximum angle among those formed between each particle in the bridge and its two neighbors: the larger the maximum angle is, the weaker the bridge. It has also been reported that the breaking time distribution shows a power-law tail with an exponent that depends on the outlet size, the vibration intensity, and the load [I. Zuriguel et al., Sci. Rep. 4, 7324 (2014)]. Here we connect these previous works, demonstrating the importance of the maximum angle in the arch on the exponent of the breaking time distribution. Besides, we find that the acceleration needed to break an arch does not depend on the ramp rate of the applied acceleration, but it does depend on the outlet size above which the arch is formed. We also show that high frequencies of vibration reveal a change in the behavior of the arches that endure very long times. These arches have been identified as a subset with special geometrical features. Therefore, arches that cannot be broken by means of a given external excitation might exist.

Granular segregation driven by particle interactions

We report the results of an experimental study of particle-particle interactions in a horizontally shaken granular layer that undergoes a second order phase transition from a binary gas to a segregation liquid as the packing fraction C is increased. By focusing on the behavior of individual particles, the effect of C is studied on (1) the process of cluster formation, (2) cluster dynamics, and (3) cluster destruction. The outcomes indicate that the segregation is driven by two mechanisms: attraction between particles with the same properties and random motion with a characteristic length that is inversely proportional to C. All clusters investigated are found to be transient and the probability distribution functions of the separation times display a power law tail, indicating that the splitting probability decreases with time.

Flow and clogging of a sheep herd passing through a bottleneck

We present an experimental study of a flock passing through a narrow door. Video monitoring of daily routines in a farm has enabled us to collect a sizable amount of data. By measuring the time lapse between the passage of consecutive animals, some features of the flow regime can be assessed. A quantitative definition of clogging is demonstrated based on the passage time statistics. These display broad tails, which can be fitted by power laws with a relatively large exponent. On the other hand, the distribution of burst sizes robustly evidences exponential behavior. Finally, borrowing concepts from granular physics and statistical mechanics, we evaluate the effect of increasing the door size and the performance of an obstacle placed in front of it. The success of these techniques opens new possibilities regarding their eventual extension to the management of human crowds.

Granular segregation in a thin drum rotating with periodic modulation

We present the results of an experimental investigation into the effects of a sinusoidal modulation of the rotation rate on the segregation patterns formed in thin drum of granular material. The modulation transforms the base pattern formed under steady conditions by splitting or merging the initial streaks. Specifically, the relation between the frequency of modulation and the rotation rate determines the number of streaks which develop from the base state. The results are in accord with those of Fiedor and Ottino [J. Fluid. Mech. 533, 223 (2005)10.1017/S0022112005003952], and we show that their ideas apply over a wide range of parameter space. Furthermore, we provide evidence that the observed relationship is maintained for filling fractions far from 50% and generalize the result in terms of the geometry of the granular deposit.

Topological analysis of tapped granular media using persistent homology

We use the first Betti number of a complex to analyze the morphological structure of granular samples in mechanical equilibrium. We investigate two-dimensional granular packings after a tapping process by means of both simulations and experiments. States with equal packing fraction obtained with different tapping intensities are distinguished after the introduction of a filtration parameter which determines the particles (nodes in the network) that are joined by an edge. This is accomplished by just using the position of the particles obtained experimentally and no other information about the possible contacts, or magnitude of forces.

Role of driving force on the clogging of inert particles in a bottleneck

We present numerical results of the effect that the driving force has on the clogging probability of inert particles passing through a bottleneck. When the driving force is increased by four orders of magnitude, the mean avalanche size remains almost unaltered (increases 1.6 times) while the flow rate and the avalanche duration display strong dependence on this magnitude. This indicates that in order to characterize the ability of a system to clog, the right variable to consider is the number of particles that pass through the outlet. The weak dependence of this magnitude on the driving force is explained in terms of the average kinetic energy of the flowing grains that has to be dissipated in order to get an arch stabilized.

Vibrot, a simple device for the conversion of vibration into rotation mediated by friction: preliminary evaluation

While “vibrational noise” induced by rotating components of machinery is a common problem constantly faced by engineers, the controlled conversion of translational into rotational motion or vice-versa is a desirable goal in many scenarios ranging from internal combustion engines to ultrasonic motors. In this work, we describe the underlying physics after isolating a single degree of freedom, focusing on devices that convert a vibration along the vertical axis into a rotation around this axis. A typical Vibrot (as we label these devices) consists of a rigid body with three or more cantilevered elastic legs attached to its bottom at an angle. We show that these legs are capable of transforming vibration into rotation by a “ratchet effect”, which is caused by the anisotropic stick-slip-flight motion of the leg tips against the ground. Drawing an analogy with the Froude number used to classify the locomotion dynamics of legged animals, we discuss the walking regime of these robots. We are able to control the rotation frequency of the Vibrot by manipulating the shaking amplitude, frequency or waveform. Furthermore, we have been able to excite Vibrots with acoustic waves, which allows speculating about the possibility of reducing the size of the devices so they can perform tasks into the human body, excited by ultrasound waves from the outside.

Contact network topology in tapped granular media

We analyze the contact network of simulated two-dimensional granular packings in different states of mechanical equilibrium obtained by tapping. We show that topological descriptors of the contact network allow one to distinguish steady states of the same mean density obtained with different tap intensities. These equal-density states were recently proven to be distinguishable through the mean force moment tensor. In contrast, geometrical descriptors, such as radial distribution functions, bond order parameters, and Voronoi cell distributions, can hardly discriminate among these states. We find that small-order loops of contacts-the polygons of the network-are especially sensitive probes for the contact structure.

Influence of the feeding mechanism on deposits of square particles

In a previous paper [Hidalgo et al., Phys. Rev. Lett. 103, 118001 (2009)] it was shown that square particles deposited in a silo tend to align with a diagonal parallel to the gravity, giving rise to a deposit with very particular properties. Here we explore, both experimentally and numerically, the effect on these properties of the filling mechanism. In particular, we modify the volume fraction of the initial configuration from which the grains are deposited. Starting from a very dilute case, increasing the volume fraction results in an enhancement of the disorder in the final deposit characterized by a decrease of the final packing fraction and a reduction of the number of particles oriented with their diagonal in the direction of gravity. However, for very high initial volume fractions, the final packing fraction increases again. This result implies that two deposits with the same final packing fraction can be obtained from very different initial conditions. The structural properties of such deposits are analyzed, revealing that, although the final volume fraction is the same, their micromechanical properties notably differ.

Flow and clogging in a silo with an obstacle above the orifice

In a recent paper [Zuriguel et al., Phys. Rev. Lett. 107, 278001 (2011)] it has been shown that the presence of an obstacle above the outlet can significatively reduce the clogging probability of granular matter pouring from a silo. The amount of this reduction strongly depends on the obstacle position. In this work, we present new measurements to analyze different outlet sizes, extending foregoing results and revealing that the effect of the obstacle is enhanced as the outlet size is increased. In addition, the effect of the obstacle position on the flow rate properties and in the geometrical features of arches is studied. These results reinforce previous evidence of the pressure reduction induced by the obstacle. In addition, it is shown how the mean avalanche size and the average flow rate are not necessarily linked. On the other hand, a close relationship is suggested between the mean avalanche size and the flow rate fluctuations.

Breaking arches with vibrations: the role of defects

We present experimental and numerical results regarding the stability of arches against external vibrations. Two-dimensional strings of mutually stabilizing grains are geometrically analyzed and subsequently submitted to a periodic forcing at fixed frequency and increasing amplitude. The main factor that determines the granular arch resistance against vibrations is the maximum angle among those formed between any particle of the arch and its two neighbors: the higher the maximum angle is, the easier it is to break the arch. On the basis of an analysis of the forces, a simple explanation is given for this dependence. From this, interesting information can be extracted about the expected magnitudes of normal forces and friction coefficients of the particles composing the arches.

Flow rate of particles through apertures obtained from self-similar density and velocity profiles

"Beverloo's law" is considered as the standard expression to estimate the flow rate of particles through apertures. This relation was obtained by simple dimensional analysis and includes empirical parameters whose physical meaning is poorly justified. In this Letter, we study the density and velocity profiles in the flow of particles through an aperture. We find that, for the whole range of apertures studied, both profiles are self-similar. Hence, by means of the functionality obtained for them the mass flow rate is calculated. The comparison of this expression with the Beverloo's one reveals some differences which are crucial to understanding the mechanism that governs the flow of particles through orifices.

Silo clogging reduction by the presence of an obstacle

We present experimental results on the effect that inserting an obstacle just above the outlet of a silo has on the clogging process. We find that, if the obstacle position is properly selected, the probability that the granular flow is arrested can be reduced by a factor of 100. This dramatic effect occurs without any remarkable modification of the flow rate or the packing fraction above the outlet, which are discarded as the cause of the change in the clogging probability. Hence, inspired by previous results of pedestrian crowd dynamics, we propose that the physical mechanism behind the clogging reduction is a pressure decrease in the region of arch formation.

Stress distribution of faceted particles in a silo after its partial discharge

We present experimental and numerical results of the effect that a partial discharge has on the morphological and micro-mechanical properties of non-spherical, convex particles in a silo. The comparison of the particle orientation after filling the silo and its subsequent partial discharge reveals important shear-induced orientation, which affects stress propagation. For elongated particles, the flow induces an increase in the packing disorder which leads to a reduction of the vertical stress propagation developed during the deposit generated prior to the partial discharge. For square particles, the flow favors particle alignment with the lateral walls promoting a behavior opposite to the one of the elongated particles: vertical force transmission, parallel to gravity, is induced. Hence, for elongated particles the flow developed during the partial discharge of the silo leads to force saturation with depth whereas for squares the flow induces hindering of the force saturation observed during the silo filling.