Self-similarity of pressure profiles during forced granular flows.

We present measurements of the vertical stress profile σ on the base of flat-bottomed cylindrical silos discharged through an orifice centered on its base. An overweight forces the material on top of the free surface. The mean bottom pressure σ(z,D,W), with z the height of the granular column, D the silo diameter, and W the mass of the overweight, increases significantly at the end of the discharge. Inspired by early models of stress distribution, we show that σ measured at z=0 can be rescaled to yield a collapse of the data, as a function of z/D, for all D and W explored. We also show that the profile σ(r) is self-similar as a function of the radial coordinate r and can be rescaled to collapse the data for different z,D, and W. Although the model correctly predicts the functional dependences, it fails in quantitative terms. These results challenge our understanding of free and forced granular flows through orifices.

Evolution of force networks during stick-slip motion of an intruder in a granular material: Topological measures extracted from experimental data.

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 do not include precise information about individual contact forces, we use an approach developed in our previous work [Basak et al., J. Eng. Mech. 147, 04021100 (2021)0733-939910.1061/(ASCE)EM.1943-7889.0002003] 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 slip precursors. In particular, the number of generators describing the structure and complexity of force networks increases consistently before slips. 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.

Universal features of the stick-slip dynamics of an intruder moving through a confined granular medium.

Experiments and simulations of an intruder dragged by a spring through a two-dimensional annulus of granular material exhibit robust force fluctuations. At low packing fractions (ϕ<ϕ_{0}), the intruder clears an open channel. Above ϕ_{0}, stick-slip dynamics develop, with an average energy release that is independent of the particle-particle and particle-base friction coefficients but does depend on the width W of the annulus and the diameter D of the intruder. A simple model predicts the dependence of ϕ_{0} on W and D, allowing for a data collapse for the average energy release as a function of ϕ/ϕ_{0}. These results pose challenges for theories of mechanical failure in amorphous materials.

Pedestrian dynamics at the running of the bulls evidence an inaccessible region in the fundamental diagram.

We characterize the dynamics of runners in the famous "Running of the Bulls" Festival by computing the individual and global velocities and densities, as well as the crowd pressure. In contrast with all previously studied pedestrian systems, we unveil a unique regime in which speed increases with density that can be understood in terms of a time-dependent desired velocity of the runners. Also, we discover the existence of an inaccessible region in the speed-density state diagram that is explained by falls of runners. With all these ingredients, we propose a generalization of the pedestrian fundamental diagram for a scenario in which people with different desired speeds coexist.

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.

Differential equation for the flow rate of discharging silos based on energy balance.

Since the early work of Hagen [G. H. L. Hagen, Ber. Verhandl. K. Preuss. Akad. Wiss. Berlin 17, 35 (1852)] and Beverloo et al. [W. Beverloo et al., Chem. Eng. Sci. 15, 260 (1961)CESCAC0009-250910.1016/0009-2509(61)85030-6], 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 to 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 due to a compensation of energy terms. Moreover, this equation can be used to explore different 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.

Intruder in a two-dimensional granular system: Effects of dynamic and static basal friction on stick-slip and clogging dynamics.

We present simulation results for an intruder pulled through a two-dimensional granular system by a spring using a model designed to mimic the experiments described by Kozlowski et al. [Phys. Rev. E 100, 032905 (2019)2470-004510.1103/PhysRevE.100.032905]. 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. We then use simulations to explore a broader range of parameter space, focusing on the friction between the particles and the base. We consider 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 behavior when the coefficient is decreased below 0.1, while static friction plays only a marginal role.

Dynamics of a grain-scale intruder in a two-dimensional granular medium with and without basal friction.

We report on a series of experiments in which a grain-sized intruder is pushed by a spring through a two-dimensional granular material composed 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.

Structure of force networks in tapped particulate systems of disks and pentagons. I. Clusters and loops.

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 and the companion paper [Kondic et al., Phys. Rev. E 93, 062903 (2016)10.1103/PhysRevE.93.062903], 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. Such differences in network properties are expected to lead to different macroscale response of the considered systems, despite the fact that averaged measures (such as force probability density function) do not show any obvious differences. 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.

Structure of force networks in tapped particulate systems of disks and pentagons. II. Persistence analysis.

In the companion paper [Pugnaloni et al., Phys. Rev. E 93, 062902 (2016)10.1103/PhysRevE.93.062902], we use classical measures based on force probability density functions (PDFs), as well as Betti numbers (quantifying the number of components, related to force chains, and loops), to describe the force networks in tapped systems of disks and pentagons. In the present work, we focus on the use of persistence analysis, which allows us to describe these networks in much more detail. This approach allows us not only to describe but also to quantify the differences between the force networks in different realizations of a system, in different parts of the considered domain, or in different systems. We show that persistence analysis clearly distinguishes the systems that are very difficult or impossible to differentiate using other means. One important finding is that the differences in force networks between disks and pentagons are most apparent when loops are considered: the quantities describing properties of the loops may differ significantly even if other measures (properties of components, Betti numbers, force PDFs, or the stress tensor) do not distinguish clearly or at all the investigated systems.

Relevance of system size to the steady-state properties of tapped granular systems.

We investigate the steady-state packing fraction ϕ and force moment tensor Σ of quasi-two-dimensional granular columns subjected to tapping. Systems of different height h and width L are considered. We find that ϕ and Σ, which describe the macroscopic state of the system, are insensitive to L for L>50d (with d the grain diameter). However, results for granular columns of different heights cannot be conciliated. This suggests that comparison between results of different laboratories on this type of experiments can be done only for systems of same height. We show that a parameter ɛ=1+(Aω)2/(2gh), with A and ω the amplitude and frequency of the tap and g the acceleration of gravity, can be defined to characterize the tap intensity. This parameter is based on the effective flight of the granular bed, which takes into account the h dependency. When ϕ is plotted as a function of ɛ, the data collapses for systems of different h. However, this parameter alone is unable to determine the steady state to be reached since different Σ can be observed for a given ɛ if different column heights are considered.

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.

Clogging transition of many-particle systems flowing through bottlenecks.

When a large set of discrete bodies passes through a bottleneck, the flow may become intermittent due to the development of clogs that obstruct the constriction. Clogging is observed, for instance, in colloidal suspensions, granular materials and crowd swarming, where consequences may be dramatic. Despite its ubiquity, a general framework embracing research in such a wide variety of scenarios is still lacking. We show that in systems of very different nature and scale -including sheep herds, pedestrian crowds, assemblies of grains, and colloids- the probability distribution of time lapses between the passages of consecutive bodies exhibits a power-law tail with an exponent that depends on the system condition. Consequently, we identify the transition to clogging in terms of the divergence of the average time lapse. Such a unified description allows us to put forward a qualitative clogging state diagram whose most conspicuous feature is the presence of a length scale qualitatively related to the presence of a finite size orifice. This approach helps to understand paradoxical phenomena, such as the faster-is-slower effect predicted for pedestrians evacuating a room and might become a starting point for researchers working in a wide variety of situations where clogging represents a hindrance.

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.

Arching during the segregation of two-dimensional tapped granular systems: mixtures versus intruders.

We present numerical simulations of binary mixtures of granular disks subjected to tapping. We consider the size segregation process in terms of the arches formed by small and big particles. Although arching has been proposed as one of the chief mechanisms that determines size segregation in non-convecting systems, there is no direct data on arching to support the existing proposals. The pseudo-dynamic approach chosen for this work allows for a straightforward identification of arches in the bulk of the column. We find that, indeed, arch formation and breakage are crucial to the segregation process. Our results show that the presence of large particles induce the formation of more arches than found in mono-sized samples. However, tapping leads to the progressive breakage of big arches where large particles are involved as the segregation process takes place. Interestingly, isolated intruders may or may not rise under tapping depending not only on the size ratio (as it is well known) but also on the degree of ordering of the environment.

First-order phase transition during displacement of amphiphilic biomacromolecules from interfaces by surfactant molecules.

The adsorption of surfactants onto a hydrophobic interface, already laden with a fixed number of amphiphilic macromolecules, is studied using the self consistent field calculation method of Scheutjens and Fleer. For biopolymers having unfavourable interactions with the surfactant molecules, the adsorption isotherms show an abrupt jump at a certain value of surfactant bulk concentration. Alternatively, the same behaviour is exhibited when the number of amphiphilic chains on the interface is decreased. We show that this sudden jump is associated with a first-order phase transition, by calculating the free energy values for the stable and the metastable states at both sides of the transition point. We also observe that the transition can occur for two approaching surfaces, from a high surfactant coverage phase to a low surfactant coverage one, at sufficiently close separation distances. The consequence of this finding for the steric colloidal interactions, induced by the overlap of two biopolymer + surfactant films, is explored. In particular, a significantly different interaction, in terms of its magnitude and range, is predicted for these two phases. We also consider the relevance of the current study to problems involving the competitive displacement of proteins by surfactants in food colloid systems.

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.

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.

Creep motion of a model 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 springs we observe a reptation of the chain. We account for the average reptation velocity theoretically. The velocity of small systems exhibits a series of plateaus as a function of the incline angle. Due to elastic effects, there exists a critical amplitude below which the reptation is expected to cease. However, rather than a full stop of the creep, we observe in numerical simulations a transition between a continuous-creep and an irregular-creep regime when the critical amplitude is approached. The latter transition is reminiscent of the transition between the continuous and the irregular compaction of granular matter submitted to periodic temperature changes.