Search

The flow of dry granular materials from silos is of great practical interest in industry and of theoretical import for understanding multiphase dynamics. Recent studies have demonstrated that one way to control the rate of flow from a silo is to tilt it. However, this may not be practical in many industrial applications. Here, we demonstrate in experiments of quasi-2D silo discharge of monodisperse grains that the flow rate can be modulated by rotating the orifice - through elevating and shifting one side of the base - instead of tilting the entire silo. We use high-speed image analysis to track the average motion of grains in the silo. We first show that the flow rate decreases with orifice angle, but that this decrease is not as strong as when a silo is tilted or when a lateral orifice is used. However, with the addition of a grain-sized ridge on each side of the orifice, the flow rate collapses with prior tilted-silo results. We then characterize the flow velocity of grains exiting the orifice and highlight key features of the stagnant zones and slip zones on each orifice side. Finally, we model our results based on these measurements, demonstrating the importance of horizontal creep along slip zones next to the orifice and the narrowest opening cross-section through which the material flows. These findings reveal a simple method for controlling both flow rate and direction, and highlight the importance of both dynamics within and geometry of the stagnant zones near the orifice., Comment: 25 pages total: 2 pages with graphical abstract and highlights for preprint, 18 pages for main text, and 5 pages for supplemental materials. 16 figures in main text, 4 figures in supplemental materials. (Videos will be provided with Supplemental Materials upon peer-reviewed publication.)

In quasi-two-dimensional experiments with photoelastic particles confined to an annular region, an intruder constrained to move in a circular path halfway between the annular walls experiences stick-slip dynamics. We discuss the response of the granular medium to the driven intruder, focusing on the evolution of the force network during sticking periods. Because the available experimental data does not include precise information about individual contact forces, we use an approach developed in our previous work (Basak et al, J. Eng. Mechanics (2021)) based on networks constructed from measurements of the integrated strain magnitude on each particle. These networks are analyzed using topological measures based on persistence diagrams, revealing that force networks evolve smoothly but in a nontrivial manner throughout each sticking period, even though the intruder and granular particles are stationary. Characteristic features of persistence diagrams show identifiable changes as a slip is approaching, indicating the existence of slip precursors. Key features of the dynamics are similar for granular materials composed of disks or pentagons, but some details are consistently different. In particular, we find significantly larger fluctuations of the measures computed based on persistence diagrams, and therefore of the underlying networks, for systems of pentagonal particles.

The flow of granular materials through constricted openings is important in many natural and industrial processes. These complex flows - featuring dense, dissipative flow in the bulk but low-dissipation, low density outpouring in the vicinity of the orifice - have long been characterized empirically by the Beverloo rule and, recently, modeled successfully using energy balance. The dependence of flow rate on the silo's angle with respect to gravity, however, is not captured by current models. We experimentally investigate the role of tilt angle in this work using a quasi-2D monolayer of grains in a silo. We measure mass flow rate, the average exit velocities of grains, and the packing fraction along the orifice with varying tilt angles. We propose a model that describes our results (and earlier findings with 3D systems [H. G. Sheldon and D. J. Durian, Granul. Matter 12, 579 (2010)]) by considering the dependence of outpouring speed and angle with respect to the orifice angle and, importantly, the angle of stagnant zones adjacent to the orifice. We conclude by posing questions about possible extensions of our model in order to describe spatial variations of exit velocity and density along the orifice cross section., Comment: 10 pages and 11 figures

The flow of grains through narrow apertures posses an extraordinary challenge: clogging. Strategies to alleviate the effect of clogging, such as the use of external vibration, are always part of the design of machinery for the handling of bulk materials. It has recently been shown that one way to reduce clogging is to use a small fraction of small particles as an additive. Besides, several works reported that self-repelling magnetic grains can flow through narrow apertures with little clogging, which suggest these are excellent candidates as "lubricating" additives for other granular materials. In this work, we study the effect of adding self-repelling magnetic particles to a sample of grains in two-dimensions. We find that, in contrast with intuition, the added magnetic grains not necessarily aid the flow of the original species.

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.

The interactions between particles in particulate systems are organized in `force networks', mesoscale features that bridge between the particle scale and the scale of the system as a whole. While such networks are known to be crucial in determining the system wide response, extracting their properties, particularly from experimental systems, is difficult due to the need to measure the interparticle forces. In this work, we show by analysis of the data extracted from simulations that such detailed information about interparticle forces may not be necessary, as long as the focus is on extracting the most dominant features of these networks. The main finding is that a reasonable understanding of the time evolution of force networks can be obtained from incomplete information such as total force on the particles. To compare the evolution of the networks based on the completely known particle interactions and the networks based on incomplete information (total force each grain) we use tools of algebraic topology. In particular we will compare simple measures defined on persistence diagrams that provide useful summaries of the force network features.

A granular gripper is a device used to hold objects by taking advantage of the phenomenon of Reynold's dilatancy. A membrane containing a granular sample is allowed to deform around the object to be held and then vacuum is used to jam the granular material inside the membrane. This allows to hold the object against external forces since deformation of the granular material is prevented by not allowing the system to increase its volume. The maximum holding force supported by the gripper depends on a number of variables. In this work, we show that in the regime of frictional holding (where the gripper does not interlock with the object), the maximum holding force does not depend on the granular material used to fill the membrane. Results for a variety of granular materials can be collapsed into a single curve if maximum holding force is plotted against the penetration depth achieved. The results suggest that the most important feature in selecting a particular granular material is its deformability to ensure an easy flow during the initial phase of the gripping process., Comment: 5 pages, 6 figures

We study, via DEM simulations, the apparent mass $m$ and loss factor $\eta$ of particle dampers (PD) attached to a vertically driven, single degree of freedom mechanical system. Although many studies focus on $\eta$, less work has been devoted to $m$. It has been recently demonstrated [M. Masmoudi \textit{et al}. Granular Matter 18 (2016) 71.] that $m$ non-linearly depends on the driving acceleration $\gamma$ according to a power law, $m\propto\gamma^k$. Experiments using 3D packings of particles suggest $k=-2$. However, simulations with 1D columns of particles on a vibrating plate and theoretical predictions based on the inelastic bouncing ball model (IBBM) suggest that $k=-1$. These findings left open questions whether m may depend on the dimensionality of the packing or on lateral interactions between walls and grains. In turn, $\eta$ was shown to follow a universal curve, $\eta\propto\gamma^{-1}$, whatever the dimensionality and the constraints in the motion of the grains. In this work, we consider PD under different confinement conditions in the motion of the particles (1D, quasi-1D, quasi-2D and 3D). We find that the dynamical response of the PD ($m$ and $\eta$) is not sensitive to the lateral confinement or dimensionality. However, we have observed two distinct regimes: (i) In the inertial regime, $\eta$ decays according to the IBBM for all dimensions, $\eta\propto\gamma^{-1}$, while $m$ falls with an apparent power law behaviour that matches Masmoudi's experiments, $m\propto\gamma^{-2}$, for all dimensions but only in the range of moderate acceleration, before becoming negative for very high accelerations. (ii) In the quasi-static regime, both $m$ and $\eta$ display a complex behavior as functions of the excitation amplitude, but tend to the IBBM prediction, $m\propto\gamma^{-1}$ and $\eta\propto\gamma^{-1}$., Comment: 21 pages, 5 figures

Since the early work of Hagen in 1852 and Beverloo et al. in 1961, the flow rate of granular material discharging through a circular orifice from a silo has been described by means of dimensional analysis and experimental fits, and explained through the "free fall arch" model. Here, in contrast with the traditional approach, we derive a differential equation based on the energy balance of the system. This equation is consistent with the well known Beverloo rule thanks to a compensation of energy terms. Moreover, this new equation can be used to explore new conditions for silo discharges. In particular, we show how the effect of friction on the flow rate can be predicted. The theory is validated using discrete element method simulations., Comment: 11 pages, 8 figures

When a granular material is freely discharged from a silo through an orifice at its base, the flow rate remains constant throughout the discharge. However, it has been recently shown that, if the discharge is forced by an overweight, the flow rate increases at the final stages of the discharge, in striking contrast to viscous fluids [Madrid et al. Europhys. Lett. (2018)]. {Although the general mechanism that drives this increase in the flow rate has been discussed, there exist yet a number of open questions regarding this phenomenon. One such questions is to what extent is the internal velocity profile affected, beyond the trivial overall increase consistent with the increasing flow rate.} We study via Discrete Element Method simulations the internal velocity profiles during forced silo discharges and compare them with those of free discharges. The changes in velocity profiles are somewhat subtle. Interestingly, during free discharges, while the velocity profiles are steady at the silo base and above a height equivalent to one silo diameter, there exists a transition region where the profile evolves in time, despite the constant flow rate. In contrast, forced discharges present steady profiles at all heights of the granular column during the initial constant flow phase, followed by an overall increase of the velocities when the acceleration phase develops., Comment: 6 pages, 5 figures

We discuss the results of simulations of an intruder pulled through a two-dimensional granular system by a spring, using a model designed to lend insight into the experimental findings described by Kozlowski et al. [Phys. Rev. E, 100, 032905 (2019)]. In that previous study the presence of basal friction between the grains and the base was observed to change the intruder dynamics from clogging to stick-slip. Here we first show that our simulation results are in excellent agreement with the experimental data for a variety of experimentally accessible friction coefficients governing interactions of particles with each other and with boundaries. Then, we use simulations to explore a broader range of parameter space, focusing on the friction between the particles and the base. We consider a range of both static and dynamic basal friction coefficients, which are difficult to vary smoothly in experiments. The simulations show that dynamic friction strongly affects the stick-slip behaviour when the coefficient is decreased below 0.1, while static friction plays only a marginal role in the intruder dynamics., Comment: 19 pages, 10 figures, 2 tables

We report on a series of experiments in which a grain-sized intruder is pushed by a spring through a 2D granular material comprised of photoelastic disks in a Couette geometry. We study the intruder dynamics as a function of packing fraction for two types of supporting substrates: a frictional glass plate and a layer of water for which basal friction forces are negligible. We observe two dynamical regimes: intermittent flow, in which the intruder moves freely most of the time but occasionally gets stuck, and stick-slip dynamics, in which the intruder advances via a sequence of distinct, rapid events. When basal friction is present, we observe a smooth crossover between the two regimes as a function of packing fraction, and we find that reducing the interparticle friction coefficient causes the stick-slip regime to shift to higher packing fractions. When basal friction is eliminated, we observe intermittent flow at all accessible packing fractions. For all cases, we present results for the statistics of stick events, the intruder velocity, and the force exerted on the intruder by the grains. Our results indicate the qualitative importance of basal friction at high packing fractions and suggest a possible connection between intruder dynamics in a static material and clogging dynamics in granular flows., Comment: 10 pages, 11 figures, 1 table. V2: What is now Fig. 5 was added to clarify the raw data measured. References [27] and [37] were updated since they were published only recently

Static granular packs have been studied in the last three decades in the frame of a modified equilibrium statistical mechanics that assumes ergodicity as a basic postulate. The canonical example on which this framework is tested consists in the series of static configurations visited by a granular column subjected to taps. By analyzing the response of a realistic model of grains, we demonstrate that volume and stress variables visit different regions of the phase space at low tap intensities in different realizations of the experiment. We show that the tap intensity beyond which sampling by tapping becomes ergodic coincides with the forcing necessary to break all particle-particle contacts during each tap. These results imply that the well-known "reversible" branch of tapped granular columns is only valid at relatively high tap intensities., Comment: 7 pages, 5 figures

The Beverloo scaling for the gravity flow of granular materials through orifices has two distinct universal features. On the one hand, the flow rate is independent of the height of the granular column. On the other hand, less well-known yet more striking, the flow rate is fairly insensitive to the material properties of the grains (density, Young's modulus, friction coefficient, etc.). We show that both universal features are lost if work is done on the system at a high rate. In contrast to viscous fluids, the flow rate increases during discharge if a constant pressure is applied to the free surface of a granular column. Moreover, the flow rate becomes sensitive to the material properties. Nevertheless, a new universal feature emerges: the dissipated power scaled by the mean pressure and the flow rate follows a master curve for forced and unforced conditions and for all material properties studied. We show that this feature can be explained if the granular flow in the silo is assumed to be a quasistatic shear flow under the $\mu(I)$-rheology., Comment: 6 pages, 4 figures

We present a differential equation for the flow rate of granular materials during the discharge of a silo. This is based in the energy balance of the variable mass system in contrast with the traditional derivations based on heuristic postulates such as the free fall arch. We show that this new equation is consistent with the well known Beverloo rule, providing an independent estimate for the universal Beverloo prefactor. We also find an analytic expression for the pressure under discharging conditions., Comment: 4 pages, 2 figures

The force network of a granular assembly, defined by the contact network and the corresponding contact forces, carries valuable information about the state of the packing. Simple analysis of these networks based on the distribution of force strengths is rather insensitive to the changes in preparation protocols or to the types of particles. In this paper we consider two dimensional simulations of tapped systems built from frictional disks and pentagons, and study the structure of the force networks of granular packings by considering network's topology as force thresholds are varied. We show that the number of clusters and loops observed in the force networks as a function of the force threshold are markedly different for disks and pentagons if the tangential contact forces are considered, whereas they are surprisingly similar for the network defined by the normal forces. In particular, the results indicate that, overall, the force network is more heterogeneous for disks than for pentagons. Additionally, we show that the states obtained by tapping with different intensities that display similar packing fraction are difficult to distinguish based on simple topological invariants.

We implement the Wang-Landau algorithm to sample with equal probabilities the static configurations of a model granular system. The "non-interacting rigid arch model" used is based on the description of static configurations by means of splitting the assembly of grains into sets of stable arches. This technique allows us to build the entropy as a function of the volume of the packing for large systems. We make a special note of the details that have to be considered when defining the microstates and proposing the moves for the correct sampling in these unusual models. We compare our results with previous exact calculations of the model made at moderate system sizes. The technique opens a new opportunity to calculate the entropy of more complex granular models., Comment: 9 pages, 3 figures

The "faster-is-slower" (FIS) effect was first predicted by computer simulations of the egress of pedestrians through a narrow exit [Helbing D, Farkas I J, Vicsek T, Nature 407:487-490 (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.

We report two-dimensional simulations of circular and polygonal grains passing through an aperture at the bottom of a silo. The mass flow rate for regular polygons is lower than for disks as observed by other authors. We show that both the exit velocity of the grains and the packing fraction are lower for polygons, which leads to the reduced flow rate. We point out the importance of the criteria used to define when two objects of different shape are considered to be of the same size. Depending on this criteria, the mass flow rate may vary significantly for some polygons. Moreover, the particle flow rate is non-trivially related to a combination of mass flow rate, particle shape and particle size. For some polygons, the particle flow rate may be lower or higher than that of the corresponding disks depending on the size comparison criteria., Comment: 9 pages, 8 figures