Role of vibrations in the jamming and unjamming of grains discharging from a silo

Flow dynamics of different particle shapes in a rectangular silo

The present work investigates the effect of using six different particle shapes of equal volume on the discharge process of a rectangular silo with adjustable width, equipped with a flat bottom orifice opening of varying size. We find that the discharge rate decreases with the increasing aspect ratio of the particles for both lentil-shaped (oblate) and rice-shaped (prolate ellipsoidal) particles and macaroni-shaped particles show the lowest discharge rate among all the particle shapes. In addition, the silo width influences the discharge in such a way that the rates at which different particle shapes flow out from the system become more distinguishable at smaller silo widths. We observe that the velocity profile near the orifice opening becomes narrower and less sharp with increasing aspect ratio for both lentil- and rice-shaped particles. Moreover, the silo width does not have a significant influence on the velocity profile very near to the orifice, but, its influence becomes more noticeable with increasing height within the silo.

Reduction of clogging of vibrated grains passing through a narrow aperture by the addition of low-friction particles

We study the flow of grains under vibration passing through a small aperture in two dimensions using discrete element method simulations. Such systems are prone to clogging and strategies to ease the flow are desirable in multiple applications. We show that the addition of low-friction particles to the system can reduce clogging and lead to an enhancement of the net flux of the original species. Along with the role of the particle friction, we analyze the influence of both the size of the added particles and the mixing ratio. We consider systems with constant height of the granular column (using particle reinjection) as well as processes of fully emptying the containing hopper.

Discharge of rice-shaped particles from a monolayer flat-bottom silo

In this work, we performed experiments regarding the outflow of spheres and two different types of rice-shaped particles in a quasi-two-dimensional monolayer silo with a flat bottom. We investigate the velocity and solid fraction profiles at the orifice and test whether the profiles for nonspherical particles have similar self-similar properties as in the spherical case. We find that the magnitude and shape of the velocity profiles for all three particle types are in a similar range. In contrast, the solid fraction at the orifice has a dome-shaped profile for both rice particles, whereas the profile for spherical particles is rather flat. The discharge rate determined from the velocity and solid fraction profiles describes the independently measured experimental discharge rate very well for all three investigated particle types.

Nonlinear effect of grain elongation on the flow rate in silo discharge

By means of two-dimensional numerical simulations based on contact dynamics, we present a systematic analysis of the joint effects of grain shape (i.e., grain elongation) and system size on silo discharge for increasing orifice sizes D. Grains are rounded-cap rectangles whose aspect ratio are varied from 1 (disks) to 7. In order to clearly isolate the effect of grain shape, the mass of the grains is keeping constant as well as the condition of the discharge by reintroducing the exiting grains at the top of the silo. In order to quantify the possible size effects, the thickness W of the silos is varied from 7 to 70 grains diameter, while keeping the silos aspect ratio always equal to 2. We find that, as long as size effects are negligible, the flow rate Q increases as a Beverloo-like function with D, also for the most elongated grains. In contrast, the effects of grain elongation on the flow rate depend on orifice size. For small normalized orifice sizes, the flow rate is nearly independent with grain elongation. For intermediate normalized orifice sizes the flow rate first increases with grain elongation up to a maximum value that depends on the normalized size of the orifice and saturates as the grains become more elongated. For larger normalized orifice size, the flow rate is an increasing function of grains' aspect ratio. Velocity profiles and packing fraction profiles close to the orifice turn out to be self-similar for all grain shapes and for the whole range of orifice and system sizes studied. Following the methodology introduced by Janda et al. [Phys. Rev. Lett. 108, 248001 (2012)PRLTAO0031-900710.1103/PhysRevLett.108.248001], we explain the nonlinear variation of Q with grain elongation, and for all orifice sizes, from compensation mechanisms between the velocity and packing fraction measured at the center of the orifice. Finally, an equation to predict the evolution of Q as a function of the aspect ratio of the grains is deduced.

Elongated particles discharged with a conveyor belt in a two-dimensional silo

The flow of elliptical particles out of a two-dimensional silo when extracted with a conveyor belt is analyzed experimentally. The conveyor belt-placed directly below the silo outlet-reduces the flow rate, increases the size of the stagnant zone, and it has a very strong influence on the relative velocity fluctuations as they strongly increase everywhere in the silo with decreasing belt speed. In other words, instead of slower but smooth flow, flow reduction by belt leads to intermittent flow. Interestingly, we show that this intermittency correlates with a strong reduction of the orientational order of the particles at the orifice region. Moreover, we observe that the average orientation of the grains passing through the outlet is modified when they are extracted with the belt, a feature that becomes more evident for large orifices.

Flow rate resonance of actively deforming particles

Lymphoid organs are unusual multicellular tissues: they are densely packed, but the lymphocytes trafficking through them are actively moving. We hypothesize that the intriguing ability of lymphocytes to avoid jamming and clogging is in part attributable to the dynamic shape changes that cells undergo when they move. In this work, we test this hypothesis by investigating an idealized system, namely, the flow of self-propelled, oscillating particles passing through a narrow constriction in two dimensions (2D), using numerical simulations. We found that deformation allows particles with these properties to flow through a narrow constriction in conditions when non-deformable particles would not be able to do so. Such a flowing state requires the amplitude and frequency of oscillations to exceed threshold values. Moreover, a resonance leading to the maximum flow rate was found when the oscillation frequency matched the natural frequency of the particle related to its elastic stiffness. To our knowledge, this phenomenon has not been described previously. Our findings could have important implications for understanding and controlling flow in a variety of systems in addition to lymphoid organs, such as granular flows subjected to vibration.

Discharge of elongated grains in silos under rotational shear

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)2470-004510.1103/PhysRevE.103.062905]. 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.

Silo discharge of mixtures of soft and rigid grains

We study the outflow dynamics and clogging phenomena of mixtures of soft, elastic low-friction spherical grains and hard frictional spheres of similar size in a quasi-two-dimensional (2D) silo with narrow orifice at the bottom. Previous work has demonstrated the crucial influence of elasticity and friction on silo discharge. We show that the addition of small amounts, even as low as 5%, of hard grains to an ensemble of soft, low-friction grains already has significant consequences. The mixtures allow a direct comparison of the probabilities of the different types of particles to clog the orifice. We analyze these probabilities for the hard, frictional and the soft, slippery grains on the basis of their participation in the blocking arches, and compare outflow velocities and durations of non-permanent clogs for different compositions of the mixtures. Experimental results are compared with numerical simulations. The latter strongly suggest a significant influence of the inter-species particle friction.

Enhancement of the flow of vibrated grains through narrow apertures by addition of small particles

We analyze the flow and clogging of circular grains passing through a small aperture under vibration in two dimensions. Via discrete element method simulations, we show that when grains smaller than the original ones are introduced in the system as an additive, the net flow of the original species can be significantly increased. Moreover, there is an optimal radius of the additive particles that maximizes the effect. This finding may constitute the basis for technological applications not only concerning the flow of granular materials but also regarding active matter, including pedestrian evacuation.

Discharge of a granular silo under mechanical vibrations

In this paper, we study the flow rate of model granular material in a silo under the influence of mechanical vibrations. Experimental measurements and discrete element simulations (DEM) are performed in a quasi-2D silo. The influence on the flow rate of the opening size and the vibration applied on the entire silo is studied. Two distinct regimes are evidenced, governed by the Froude number Fr and the relative frequency Ω. In the first regime, a decreased flow rate is observed when increasing the vibration intensity. This behavior is explained by the presence of reorganizations induced by the vibration, leading to a more homogeneous but also slower flow. In the second regime, an increased flow rate is evidenced when increasing the vibration intensity. We find this behavior comes from the intermittent nature of the flow, where the flow rate is directly controlled by the propagation of shock waves all along the silo.

Particle flow rate in silos under rotational shear

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)10.1103/PhysRevE.100.012906]. 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 nonmonotonic behavior of Q when varying the rotational shear.

Intermittent flow and transient congestions of soft spheres passing narrow orifices

Soft, low-friction particles in silos show peculiar features during their discharge. The outflow velocity and the clogging probability both depend upon the momentary silo fill height, in sharp contrast to silos filled with hard particles. The reason is the fill-height dependence of the pressure at the orifice. We study the statistics of silo discharge of soft hydrogel spheres. The outflow is found to become increasingly fluctuating and even intermittent with decreasing orifice size, and with decreasing fill height. In orifices narrower than two particle diameters, outflow can stop completely, but in contrast to clogs formed by rigid particles, these congestions may dissolve spontaneously. We analyze such non-permanent congestions and attribute them to slow reorganization processes in the container, caused by viscoelasticity of the material.

Power-Law Scaling in Granular Rheology across Flow Geometries

Based on discrete element method simulations, we propose a new form of the constitutive equation for granular flows independent of packing fraction. Rescaling the stress ratio μ by a power of dimensionless temperature Θ makes the data from a wide set of flow geometries collapse to a master curve depending only on the inertial number I. The basic power-law structure appears robust to varying particle properties (e.g., surface friction) in both 2D and 3D systems. We show how this rheology fits and extends frameworks such as kinetic theory and the nonlocal granular fluidity model.

Transition from clogging to continuous flow in constricted particle suspensions

When suspended particles are pushed by liquid flow through a constricted channel, they might either pass the bottleneck without trouble or encounter a permanent clog that will stop them forever. However, they may also flow intermittently with great sensitivity to the neck-to-particle size ratio D/d. In this Rapid Communication, we experimentally explore the limits of the intermittent regime for a dense suspension through a single bottleneck as a function of this parameter. To this end, we make use of high time- and space-resolution experiments to obtain the distributions of arrest times (T) between successive bursts, which display power-law tails (∝T^{-α}) with characteristic exponents. These exponents compare well with the ones found for as disparate situations as the evacuation of pedestrians from a room, the entry of a flock of sheep into a shed, or the discharge of particles from a silo. Nevertheless, the intrinsic properties of our system (i.e., channel geometry, driving and interaction forces, particle size distribution) seem to introduce a sharp transition from a clogged state (α≤2) to a continuous flow, where clogs do not develop at all. This contrasts with the results obtained in other systems where intermittent flow, with power-law exponents above two, were obtained.

Granular flow through an aperture: Influence of obstacles near the outlet

We study how the presence of obstacles in a confined system of monodisperse disks affects their discharge through an aperture. The disks are driven by a horizontal conveyor belt that moves at constant velocity. The mean packing fraction at the outlet decreases as the distance between the obstacles and the aperture decreases. The obstacles organize the dynamics of the stagnant zones in two characteristic behaviors that differ mainly in the magnitude of the fluctuations of the fraction of stagnant disks in the system. It is shown that the effective aperture is reduced by the presence of obstacles.

Nonergodicity in silo unclogging: Broken and unbroken arches

We report an experiment on the unclogging dynamics in a two-dimensional silo submitted to a sustained gentle vibration. We find that arches present a jerking motion where rearrangements in the positions of their beads are interspersed with quiescent periods. This behavior occurs for both arches that break down and those that withstand the external perturbation: Arches evolve until they either collapse or get trapped in a stable configuration. This evolution is described in terms of a scalar variable characterizing the arch shape that can be modeled as a continuous-time random walk. By studying the diffusivity of this variable, we show that the unclogging is a weakly nonergodic process. Remarkably, arches that do not collapse explore different configurations before settling in one of them and break ergodicity much in the same way than arches that break down.

Granular flow from silos with rotating orifice

For dry granular materials falling through a circular exit at the bottom of a silo, no continuous flow can be sustained when the diameter D of the exit is less than five times the characteristic size of the grains. If the bottom of the silo rotates horizontally with respect to the wall of the silo, finite flow rate can be sustained even at small D. We investigate the effect of bottom rotation to the flow rate of monodisperse plastic beads of d=6mm diameter from a cylindrical silo of 19 cm inner diameter. We find that the flow rate W follows Beverloo law down to D=1.3d and that W increases with the rotation speed ω in the small exit regime. If the exit is at an off-center distance R from the axis of the silo, W increases with the rate of area swept by the exit. On the other hand, when the exit diameter is large, W decreases with increasing ω at small ω but increases with ω at large ω. Such nonmonotonic dependence of flow rate on rotation speed may be explained as a gradual change from funnel flow to mass flow due to the shear at the bottom of the silo.

Flow and clogging of particles in shaking random obstacles

Transport of three types of particles (passive particles, active particles, and polar particles) is investigated in a random obstacle array in the presence of a dc drift force. The obstacles are static or synchronously shake along the given direction. When the obstacles are static, the average velocity is a peaked function of the dc drift force (negative differential mobility) for low particle density, while the average velocity monotonically increases with the dc drift force (positive differential mobility) for high particle density. Under the same conditions, passive particles are most likely to pass through the obstacles, while polar particles are easily trapped by the obstacles. The polar alignment can strongly reduce the overall mobility of particles. When the obstacles shake along the given direction, the optimal shaking frequency or amplitude can maximize the average velocity. It is more effective to reduce clogging for the transverse shaking than that for the longitudinal shaking.

Controlled Fluidization, Mobility, and Clogging in Obstacle Arrays Using Periodic Perturbations

We show that the clogging susceptibility and flow of particles moving through a random obstacle array can be controlled with a transverse or longitudinal ac drive. The flow rate can vary over several orders of magnitude, and we find both an optimal frequency and an optimal amplitude of driving that maximizes the flow. For dense arrays, at low ac frequencies, a heterogeneous creeping clogged phase appears in which rearrangements between different clogged configurations occur. At intermediate frequencies, a high-mobility fluidized state forms, and, at high frequencies, the system reenters a heterogeneous frozen clogged state. These results provide a technique for optimizing flow through heterogeneous media that could also serve as the basis for a particle separation method.

Ergodicity breaking dynamics of arch collapse

Flows in hoppers and silos are susceptible to clogging due to the formation of arches at the exit. The failure of these arches is the key to reinitiation of flow, yet the physical mechanism of failure is not well understood. Experiments on vibrated hoppers exhibit a broad distribution of the duration of clogs. Using numerical simulations of a hopper in two dimensions, we show that arches become trapped in locally stable shapes that are explored dynamically under vibrations. The shape dynamics, preceding failure, break ergodicity and can be modeled as a continuous-time random walk with a broad distribution of waiting, or trapping, times. We argue that arch failure occurs as a result of this random walk crossing a stability boundary, which is a first-passage process that naturally gives rise to a broad distribution of unclogging times.

Trap Model for Clogging and Unclogging in Granular Hopper Flows

Granular flows through narrow outlets may be interrupted by the formation of arches or vaults that clog the exit. These clogs may be destroyed by vibrations. A feature which remains elusive is the broad distribution p(τ) of clog lifetimes τ measured under constant vibrations. Here, we propose a simple model for arch breaking, in which the vibrations are formally equivalent to thermal fluctuations in a Langevin equation; the rupture of an arch corresponds to the escape from an energy trap. We infer the distribution of trap depths from experiments made in two-dimensional hoppers. Using this distribution, we show that the model captures the empirically observed heavy tails in p(τ). These heavy tails flatten at large τ, consistently with experimental observations under weak vibrations. But, here, we find that this flattening is systematic, which casts doubt on the ability of gentle vibrations to restore a finite outflow forever. The trap model also replicates recent results on the effect of increasing gravity on the statistics of clog formation in a static silo. Therefore, the proposed framework points to a common physical underpinning to the processes of clogging and unclogging, despite their different statistics.

Slow relaxation dynamics of clogs in a vibrated granular silo

We experimentally explore the vibration-induced unclogging of arches halting the flow in a two-dimensional silo. The endurance of arches is determined by carrying out a survival analysis of their breaking times. By analyzing the dynamics of two morphological variables, we demonstrate that arches evolve toward less regular structures and it seems that there may exist a certain degree of irregularity that the arch reaches before collapsing. Moreover, we put forward that σ (the standard deviation of all angles between consecutive beads) describes faithfully the morphological evolution of the arch. Focusing on long-lasting arches, we study σ calculating its two-time autocorrelation function and its mean-squared displacement. In particular, the apparent logarithmic increase of the correlation and the decrease of the mean-squared displacement of σ when the waiting time is increased reveal a slowing down of the dynamics. This behavior is a clear hallmark of aging phenomena and confirms the lack of ergodicity in the unclogging dynamics. Our findings provide new insights on how an arch tends to destabilize and how the probability that it breaks with a long sustained vibration decreases with time.

Effect of form of obstacle on speed of crowd evacuation

This paper investigates the effect of the form of an obstacle on the time that a crowd takes to evacuate a room, using a toy model. Pedestrians are modeled as active soft matter moving toward a point with intended velocities. An obstacle is placed in front of the exit, and it has one of four shapes: a cylindrical column, a triangular prism, a quadratic prism, or a diamond prism. Numerical results indicate that the evacuation-completion time depends on the shape of the obstacle. Obstacles with a circular cylinder (C.C.) shape yield the shortest evacuation-completion time in the proposed model.

Clogging Transition of Vibration-Driven Vehicles Passing through Constrictions

We report experimental results on the competitive passage of elongated self-propelled vehicles rushing through a constriction. For the chosen experimental conditions, we observe the emergence of intermittencies similar to those reported previously for active matter passing through narrow doors. Noteworthy, we find that, when the number of individuals crowding in front of the bottleneck increases, there is a transition from an unclogged to a clogged state characterized by a lack of convergence of the mean clog duration as the measuring time increases. It is demonstrated that this transition-which was reported previously only for externally vibrated systems such as colloids or granulars-appears also for self-propelled agents. This suggests that the transition should also occur for the flow through constrictions of living agents (e.g., humans and sheep), an issue that has been elusive so far in experiments due to safety risks.

Flow and clog in a silo with oscillating exit

When grains flow out of a silo, flow rate W increases with exit size D. If D is too small, an arch may form and the flow may be blocked at the exit. To recover from clogging, the arch has to be destroyed. Here we construct a two-dimensional silo with movable exit and study the effects of exit oscillation (with amplitude A and frequency f) on flow rate, clogging, and unclogging of grains through the exit. We find that, if exit oscillates, W remains finite even when D (measured in unit of grain diameter) is only slightly larger than one. Surprisingly, while W increases with oscillation strength Γ≡4π^{2}Af^{2} as expected at small D, W decreases with Γ when D≥5 due to induced random motion of the grains at the exit. When D is small and oscillation speed v≡2πAf is slow, temporary clogging events cause the grains to flow intermittently. In this regime, W depends only on v-a feature consistent to a simple arch breaking mechanism, and the phase boundary of intermittent flow in the D-v plane is consistent to either a power law: D∝v^{-7} or an exponential form: D∝e^{-D/0.55}. Furthermore, the flow time statistic is Poissonian whereas the recovery time statistic follows a power-law distribution.

An experimental investigation on the speed of sand flow through a fixed porous bed

In this study, we document several experiments that investigate the speed of the flow of fine sand through a fixed porous bed of packed glass beads under various conditions, including the height of the sand column (H) and porous bed (h) and the diameter of the glass beads (D) and sand grains (d). The experiments are conducted with glass beads packed at a constant density and sand at a different dry bulk density. The results show that the height of the sand does not affect the speed of the sand flow. The speed of the sand flow (v) decreases with an increase in h until h approaches a certain value. An equation [Formula: see text] is proposed, where a and k are the experimentally determined coefficients. Moreover, the flow of sand through a fixed porous bed could be regarded as parallel flow through multiple pipes, and therefore, the relationship between D and the number and diameter of pipes, N and D _p, are discussed. Further investigations are needed for the result that the flow of sand through a porous bed or multiple parallel pipes cannot be simplified to flow through one orifice with a certain diameter.

Clogging and jamming transitions in periodic obstacle arrays

We numerically examine clogging transitions for bidisperse disks flowing through a two-dimensional periodic obstacle array. We show that clogging is a probabilistic event that occurs through a transition from a homogeneous flowing state to a heterogeneous or phase-separated jammed state where the disks form dense connected clusters. The probability for clogging to occur during a fixed time increases with increasing particle packing and obstacle number. For driving at different angles with respect to the symmetry direction of the obstacle array, we show that certain directions have a higher clogging susceptibility. It is also possible to have a size-specific clogging transition in which one disk size becomes completely immobile while the other disk size continues to flow.

Outflow and clogging of shape-anisotropic grains in hoppers with small apertures

Outflow of granular material through a small orifice is a fundamental process in many industrial fields, for example in silo discharge, and in everyday's life. Most experimental studies of the dynamics have been performed so far with monodisperse disks in two-dimensional (2D) hoppers or spherical grains in 3D. We investigate this process for shape-anisotropic grains in 3D hoppers and discuss the role of size and shape parameters on avalanche statistics, clogging states, and mean flow velocities. It is shown that an increasing aspect ratio of the grains leads to lower flow rates and higher clogging probabilities compared to spherical grains. On the other hand, the number of grains forming the clog is larger for elongated grains of comparable volumes, and the long axis of these blocking grains is preferentially aligned towards the center of the orifice. We find a qualitative transition in the hopper discharge behavior for aspect ratios larger than ≈6. At still higher aspect ratios >8-12, the outflowing material leaves long vertical holes in the hopper that penetrate the complete granular bed. This changes the discharge characteristics qualitatively.

Intermittency and velocity fluctuations in hopper flows prone to clogging

We study experimentally the dynamics of granular media in a discharging hopper. In such flows, there often appears to be a critical outlet size D_{c} such that the flow never clogs for D>D_{c}. We report on the time-averaged velocity distributions, as well as temporal intermittency in the ensemble-averaged velocity of grains in a viewing window, for both D<D_{c} and D>D_{c}, near and far from the outlet. We characterize the velocity distributions by the standard deviation and the skewness of the distribution of vertical velocities. We propose a measure for intermittency based on the two-sample Kolmogorov-Smirnov D_{KS} statistic for the velocity distributions as a function of time. We find that there is no discontinuity or kink in these various measures as a function of hole size. This result supports the proposition that there is no well-defined D_{c} and that clogging is always possible. Furthermore, the intermittency time scale of the flow is set by the speed of the grains at the hopper exit. This latter finding is consistent with a model of clogging as the independent sampling for stable configurations at the exit with a rate set by the exiting grain speed [C. C. Thomas and D. J. Durian, Phys. Rev. Lett. 114, 178001 (2015)PRLTAO0031-900710.1103/PhysRevLett.114.178001].

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.

Flux of granular particles through a shaken sieve plate

We experimentally investigate a discharging flux of granular particles through a sieve plate subject to vertical vibrations. The mean mass flux shows a non-monotonic relation with the vibration strength. High-speed photography reveals that two stages, the free flight of the particles’ bulk over the plate and the adhesion of the particles’ bulk with the plate, alternately appear, where only the adhesion stage contributes to the flow. With two independent methods, we then measure the adhesion time under different vibration conditions, and define an adhesion flux. The adhesion flux monotonically increases with increasing vibration strength. By rescaling the adhesion flux, we find that the adhesion flux is approximately determined by the peak vibration velocity of the shaker. The conclusion is examined with other sieve geometries.

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.

On-and-off dynamics of a creeping frictional system

We report on the dynamics of a model frictional system submitted to minute external perturbations. The system consists of a chain of sliders connected through elastic springs that rest on an incline. By introducing cyclic expansions and contractions of the rest length of the springs, we induce the reptation of the chain. Decreasing the amplitude of the perturbation below a critical value, we observe an intermittent creep regime characterized by alternated periods of reptation (flowing state) and rest (quiescent state). A further decrease of the perturbation leads to the disappearance of the reptation. The width of the transition region between the continuous creep and the full stop (i.e., the range of excitation amplitudes where the intermittent creep is observed) is shown to depend on the difference between the static (μ(s)) and the dynamic (μ(d)) friction coefficients. For μ(s) = μ(d) the intermittent creep is not observed. Studying the statistical features of the intermittent creep regime for any given perturbation amplitude, we find that the time the system resides in each state (flowing or quiescent) suggests that: i) reptation events are uncorrelated, and ii) rest events are history dependent. We show that this latter history dependence is consistent with the aging of the stress state inside the chain of sliders during the quiescent periods.

Granular flow through an aperture: influence of the packing fraction

For the last 50 years, the flow of a granular material through an aperture has been intensely studied in gravity-driven vertical systems (e.g., silos and hoppers). Nevertheless, in many industrial applications, grains are horizontally transported at constant velocity, lying on conveyor belts or floating on the surface of flowing liquids. Unlike fluid flows, that are controlled by the pressure, granular flow is not sensitive to the local pressure but rather to the local velocity of the grains at the outlet. We can also expect the flow rate to depend on the local density of the grains. Indeed, vertical systems are packed in dense configurations by gravity, but, in contrast, in horizontal systems the density can take a large range of values, potentially very small, which may significantly alter the flow rate. In the present article, we study, for different initial packing fractions, the discharge through an orifice of monodisperse grains driven at constant velocity by a horizontal conveyor belt. We report how, during the discharge, the packing fraction is modified by the presence of the outlet, and we analyze how changes in the packing fraction induce variations in the flow rate. We observe that variations of packing fraction do not affect the velocity of the grains at the outlet, and, therefore, we establish that flow-rate variations are directly related to changes in the packing fraction.

Critical bottleneck size for jamless particle flows in two dimensions

We propose a simple microscopic model for arching phenomena at bottlenecks. The dynamics of particles in front of a bottleneck is described by a one-dimensional stochastic cellular automaton on a semicircular geometry. The model reproduces oscillation phenomena due to the formation and collapsing of arches. It predicts the existence of a critical bottleneck size for continuous particle flows. The dependence of the jamming probability on the system size is approximated by the Gompertz function. The analytical results are in good agreement with simulations.

Geometry dependence of the clogging transition in tilted hoppers

We report the effects of system geometry on the clogging of granular material flowing out of flat-bottomed hoppers with variable aperture size D and with variable angle θ of tilt of the hopper away from horizontal. In general, larger tilt angles make the system more susceptible to clogging. To quantify this effect for a given θ, we measure the distribution of mass discharged between clogging events as a function of aperture size and extrapolate to the critical size at which the average mass diverges. By repeating for different angles, we map out a clogging phase diagram as a function of D and θ that demarcates the regimes of free flow (large D, small θ) and clogging (small D, large θ). We do this for both circular holes and long rectangular slits. Additionally, we measure four types of grain: smooth spheres (glass beads), compact angular grains (beach sand), disklike grains (lentils), and rodlike grains (rice). For circular apertures, the clogging phase diagram is found to be the same for all grain types. For narrow slit apertures and compact grains, the shape is also the same as for circular holes when expressed in terms of projected area of the aperture against the average flow direction. For lentils and rice discharged from slits, the behavior differs and may be due to alignment between grain and slit axes.

Jamming of particles in a two-dimensional fluid-driven flow

The jamming of particles under flow is of critical importance in a broad range of natural and industrial settings, such as the jamming of ice in rivers, or the plugging of suspended solids in pipeline transport. Relatively few studies have been carried out on jamming of suspended particles under flow, in comparison to the many studies on jamming in gravity-driven flows that have revealed various features of the jamming process. Fluid-driven particle flows differ in several aspects from gravity-driven flows, particularly in being compatible with a range of particle concentrations and velocities. Additionally, there are fluid-particle interactions and hydrodynamic effects. To investigate particle jamming in fluid-driven flows, we have performed both experiments and computer simulations on the flow of circular particles floating over water in an open channel with a restriction. We determined the flow-rate boundary for a dilute-to-dense flow transition, similar to that seen in gravity-driven flows. The maximum particle throughput increased for larger restriction sizes consistent with a Beverloo equation form over the entire range of particle mixtures and restriction sizes. The exponent of ~3/2 in the Beverloo equation is consistent with approximately constant acceleration of grains due to fluid drag in the immediate region of the opening. We verified that the jamming probability from the dense flow gave a geometric distribution in the number of particles escaping before a jam. The probability of jamming in both experiments and simulations was found to be dependent on the ratio of channel opening to particle size, but only weakly dependent on the fluid flow velocity. Flow entrance effects were measured and observed to affect the jamming probability, and dependence on particle friction coefficient was determined from simulation. A comprehensive model for the jamming probability integrating these observations from the different flow regimes was shown to be in good agreement for experimental data on average times before jamming.

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.

How relative humidity affects random packing experiments

The influence of relative humidity (RH) on the extremely slow compaction dynamics of a granular assembly has been experimentally investigated. Millimeter-sized glass beads are considered. Compaction curves are fitted by stretched exponentials with characteristic time τ and exponent δ, which are seen to be deeply affected by the moisture content. A kinetic model, taking into account both triboelectric and capillary effects, is in excellent agreement with our results. It confirms the existence of an optimal condition at a relative humidity ≈45% for minimizing cohesive interactions between glass beads. The exponent δ is seen to depend strongly on the diffusive character of grains and voids inside the packing: diffusion for cohesiveless particles and subdiffusion when cohesion plays a role. As a consequence, the RH represents a relevant parameter that should be reported for every experimental work on a slowly driven dense random packing.

Evolution of pressure profiles during the discharge of a silo

We report measurements of the pressure profile in the outlet plane of a discharging silo. We observe that, whatever the preparation of the granular system, a dynamic Janssen effect is at play: the apparent mass of the grains (i.e., the part of their mass sustained by the base) is significantly smaller than their actual mass because of the redirection of the weight to the lateral wall of the container. The pressure profiles reveal a significant decrease in pressure in the vicinity of the outlet as the system discharges, whereas the flow rate remains constant. The measurements are thus a direct experimental proof that the flow rates of granular material through an aperture are not controlled by the local stress conditions.

Granular flow through an aperture: pressure and flow rate are independent

We simultaneously measure the flow rate and the normal force on the base, near the outlet, during the discharge through an orifice of a dense packing of monosized disks driven by a conveyor belt. We find that the normal force on the base decreases even when a constant flow rate is measured. In addition, we show, by changing the mass of the disks, that pressure can be changed while the flow rate remains constant. Conversely, we are able, by changing the belt velocity, to set different flow rates for the same pressure. The experiment confirms that, contrary to what has been implicitly assumed in numerous works, the flow rate through an aperture is not controlled by the pressure in the outlet region.