Evolution of pressure profiles during the discharge of a silo

Granular flow of 3D mixtures of soft and hard spheres

The discharge of granular mixtures composed of hard frictional beads and soft low-friction beads was investigated in a cylindrical silo in experiments and numerical simulations. In the two limits, we find a fill height dependent flow rate for 100% low friction soft grains and a height independent flow rate for 100% hard frictional grains. When mixing the two types of grains, the transition between the two limiting cases occurs rather abruptly. Namely, adding only 20% of hard frictional grains to a sample of low friction soft grains changes the dependence of the flow rate on the discharged mass significantly, i.e. causes the slope of the curve to decrease by 50-70%. Our numerical simulations reveal that the main factor leading to the strong change in the flow rate behavior at low hard grain concentration is the high sensitivity of the stress conditions in the orifice region to the mixture composition. Since frictional dissipation can be an important factor influencing the flow rate, we also analyze the frictional properties of our samples in two additional experiments: (i) quasistatic shear tests in a split-bottom shear cell and (ii) drag force measurements on an object moved in the mixture. The mixtures show increasing dissipation as a function of increasing hard grain concentration in both of these measurements, but the increase is rather modest in the low concentration range, thus it does not explain the abrupt change in the silo discharge rate.

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.

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.

Flow-induced surface crystallization of granular particles in cylindrical confinement

An interesting phenomenon that a layer of crystallized shell formed at the container wall during an orifice flow in a cylinder is observed experimentally and is investigated in DEM simulation. Different from shear or vibration driven granular crystallization, our simulation shows during the flow the shell layer is formed spontaneously from stagnant zone at the base and grows at a constant rate to the top with no external drive. Roughness of the shell surface is defined as a standard deviation of the surface height and its development is found to disobey existed growth models. The growth rate of the shell is found linearly proportional to the flow rate. This shell is static and served as a rough wall in an orifice flow with frictionless sidewall, which changes the flow profiles and its stress properties, and in turn guarantees a constant flow rate.

Discharge of a 2D magnetic silo

We investigate experimentally the discharge of a 2D-silo, a Hele-Shaw cell, filled with a mono layer of ferromagnetic grains submitted to an external magnetic field B perpendicular to the cell plane. In this case the magnetic pair interactions are repulsive. We show that the granular flow rate decreases systematically with the amplitude of the external magnetic field applied. Interestingly, while the output flow rate remains constant during an experiment, we reveal very large spatio-temporal fluctuations of the packing density within the cell, particularly evident for magnetic field of high amplitudes.

Velocity scaling in the region of orifice influence in silo draining under gravity

This study utilizes computations based on a soft-particle discrete element method for investigating the scaling of velocity in the region of orifice influence situated directly above and in proximity to the outlet in a two-dimensional silo. The velocity at the exit scales with the outlet size (D), in striking agreement with the earlier studies. However, the scaling of velocity upstream of the outlet with D as the length scale does not exist. Consequently, we present a scaling with a length parameter h_{e} being the height of an equi-inertial curve, which is defined to be a curve on which the inertial number is constant, thereby consolidating the coexisting different flow regimes in a discharging silo. The velocity corresponding to an equi-inertial curve, when measured relative to the velocity at the outlet, scales very well with h_{e} for low inertial numbers belonging to the dense flow regime. However, such scaling does not hold for high inertial numbers corresponding to the rapid flow regime in the region located closer to the orifice. We tie this scaling breakdown to the velocity fluctuations in light of the similarity between the profiles of scaled relative velocity and the scaled kinetic pressure, suggesting h_{e} to be a promising candidate for unifying the kinematics of granular flow near the outlet in the silo.

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.

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.

Silo collapse under granular discharge

We investigate, at a laboratory scale, the collapse of cylindrical shells of radius R and thickness t induced by a granular discharge. We measure the critical filling height for which the structure fails upon discharge. We observe that the silos sustain filling heights significantly above an estimation obtained by coupling standard shell-buckling and granular stress distribution theories. Two effects contribute to stabilize the structure: (i) below the critical filling height, a dynamical stabilization due to granular wall friction prevents the localized shell-buckling modes to grow irreversibly; (ii) above the critical filling height, collapse occurs before the downward sliding motion of the whole granular column sets in, such that only a partial friction mobilization is at play. However, we notice also that the critical filling height is reduced as the grain size d increases. The importance of grain size contribution is controlled by the ratio d/√[Rt]. We rationalize these antagonist effects with a novel fluid-structure theory both accounting for the actual status of granular friction at the wall and the inherent shell imperfections mediated by the grains. This theory yields new scaling predictions which are compared with the experimental results.

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.

Continuum simulation of the discharge of the granular silo: a validation test for the μ(I) visco-plastic flow law

Using a continuum Navier-Stokes solver with the μ(I) flow law implemented to model the viscous behavior, and the discrete Contact Dynamics algorithm, the discharge of granular silos is simulated in two dimensions from the early stages of the discharge until complete release of the material. In both cases, the Beverloo scaling is recovered. We first do not attempt a quantitative comparison, but focus on the qualitative behavior of velocity and pressure at different locations in the flow. A good agreement for the velocity is obtained in the regions of rapid flows, while areas of slow creep are not entirely captured by the continuum model. The pressure field shows a general good agreement, while bulk deformations are found to be similar in both approaches. The influence of the parameters of the μ(I) flow law is systematically investigated, showing the importance of the dependence on the inertial number I to achieve quantitative agreement between continuum and discrete discharge. However, potential problems involving the systems size, the configuration and "non-local" effects, are suggested. Yet the general ability of the continuum model to reproduce qualitatively the granular behavior is found to be very encouraging.

Local stresses in the Janssen granular column

We study experimentally the distribution of local stresses in a granular material confined inside a vertical cylinder. We use an image correlation technique to measure the displacement field of the container induced by the forces exerted by the grains on the inner wall. We describe an optimization procedure based on the linear theory of elastic shells to deduce the distribution of these forces from the measured displacement field. They correspond to the stress field of the granular material close to the container's inner wall. We first confirm the validity of Janssen's description for various experiments, including the influence of the bead diameter and the effect of an additional mass on top of the granular column. We then apply this method to determine the stress field during the gravity-driven discharge of a silo through an aperture.