Granular flow down an inclined plane: Bagnold scaling and rheology

Unsteady granular flows down an inclined plane

The continuum description of granular flows is still a challenge despite their importance in many geophysical and industrial applications. We extend previous works, which have explored steady flow properties, by focusing on unsteady flows accelerating or decelerating down an inclined plane in the simple shear configuration. We solve the flow kinematics analytically, including predictions of evolving velocity and stress profiles and the duration of the transient stage. The solution shows why and how granular materials reach steady flow on slopes steeper than the angle of repose and how they decelerate on shallower slopes. The model might facilitate development of natural hazard assessment and may be modified in the future to explore unsteady granular flows in different configurations.

A dilation-driven vortex flow in sheared granular materials explains a rheometric anomaly

Granular flows occur widely in nature and industry, yet a continuum description that captures their important features is yet not at hand. Recent experiments on granular materials sheared in a cylindrical Couette device revealed a puzzling anomaly, wherein all components of the stress rise nearly exponentially with depth. Here we show, using particle dynamics simulations and imaging experiments, that the stress anomaly arises from a remarkable vortex flow. For the entire range of fill heights explored, we observe a single toroidal vortex that spans the entire Couette cell and whose sense is opposite to the uppermost Taylor vortex in a fluid. We show that the vortex is driven by a combination of shear-induced dilation, a phenomenon that has no analogue in fluids, and gravity flow. Dilatancy is an important feature of granular mechanics, but not adequately incorporated in existing models.

Flowing granular materials exhibit fluid-like features, but not all of them can be described by extending fluid mechanics. Here, the authors show vortex flow in a granular layer sheared between coaxial cylinders, and attribute it to the effect that the material moves away from the sheared region.

Depth-averaged analytic solutions for free-surface granular flows impacting rigid walls down inclines

In the present paper, flows of granular materials impacting wall-like obstacles down inclines are described by depth-averaged analytic solutions. Particular attention is paid to extending the existing depth-averaged equations initially developed for frictionless and incompressible fluids down a horizontal plane. The effects of the gravitational acceleration along the slope, and of the retarding acceleration caused by friction as well, are systematically taken into account. The analytic solutions are then used to revisit existing data on rigid walls impacted by granular flows. This approach allows establishing a complete phase diagram for granular flow-wall interaction.

Divergence of Viscosity in Jammed Granular Materials: A Theoretical Approach

A theory for jammed granular materials is developed with the aid of a nonequilibrium steady-state distribution function. The approximate nonequilibrium steady-state distribution function is explicitly given in the weak dissipation regime by means of the relaxation time. The theory quantitatively agrees with the results of the molecular dynamics simulation on the critical behavior of the viscosity below the jamming point without introducing any fitting parameter.

Cooperativity and segregation in confined flows of soft binary glasses

When a suspension containing particles of different sizes flows through a confined geometry a size gradient can be established, with large particles accumulating in the channel center. Such size separation driven by hydrodynamic interactions is expected to facilitate membrane filtration and may lead to the design of novel and innovative separation techniques. For this, a wide range of particle concentrations has to be investigated, in order to clarify whether shear-induced migration can be utilized at concentrations close to or above the colloidal glass transition, where particle motion is severely hindered and hydrodynamic interactions are screened. We explore this scenario by studying the flow of binary mixtures of soft colloidal microgels, well above their liquid-solid transition, through narrow microchannels. We find that, even though the flow becomes strongly heterogeneous, in both space and time, characterized by a large cooperativity length, size segregation still occurs. This suggests that even above the glass transition shear-induced diffusion could still be used as a fractionation mechanism, which is of great relevance for process intensification purposes.

Granular flows on a dissipative base

We study inclined channel flows of sand over a sensor-enabled composite geotextile fabric base that dissipates granular fluctuation energy. We record strain of the fabric along the flow direction with imbedded fiber-optic Bragg gratings, flow velocity on the surface by correlating grain position in successive images, flow thickness with the streamwise shift of an oblique laser light sheet, velocity depth profile through a transparent side wall using a high-speed camera, and overall discharge rate. These independent measurements at inclinations between 33∘ and 37∘ above the angle of repose at 32.1±0.8∘ are consistent with a mass flow rate scaling as the 3/2 power of the flow depth, which is markedly different than flows on a rigid bumpy boundary. However, this power changes to 5/2 when flows are forced on the sand bed below its angle of repose. Strain measurements imply that the mean solid volume fraction in the flowing layer above the angle of repose is 0.268±0.033, independent of discharge rate or inclination.

Helical inner-wall texture prevents jamming in granular pipe flows

Granular pipe flows are characterized by intermittent behavior and large, potentially destructive solid fraction variations in the transport direction. By means of particle-based numerical simulations of gravity-driven flows in vertical pipes, we show that it is possible to obtain steady material transport by adding a helical texture to the inner-wall of the pipe. The helical texture leads to a more homogeneous mass flux along the pipe, prevents the emergence of large density waves and substantially reduces the probability of plug formation thus avoiding jamming of the particulate flow. We show that the granular mass flux Q through a pipe of diameter D with a helical texture of wavelength λ follows the equation Q = Q0·{1 - B sin[arctan(2πD/λ)]}, where Q0 is the flow without helix, predicted from the well-known Beverloo equation. Our new expression yields, thus, a modification of the Beverloo equation with only one additional fit parameter, B, and describes the particle mass flux with the helical texture with excellent quantitative agreement with simulation results. Future application of the method proposed here has the potential to improve granular pipe flows in a broad range of processes without the need for energy input from any external source.

Deformation-driven diffusion and plastic flow in amorphous granular pillars

We report a combined experimental and simulation study of deformation-induced diffusion in compacted quasi-two-dimensional amorphous granular pillars, in which thermal fluctuations play a negligible role. The pillars, consisting of bidisperse cylindrical acetal plastic particles standing upright on a substrate, are deformed uniaxially and quasistatically by a rigid bar moving at a constant speed. The plastic flow and particle rearrangements in the pillars are characterized by computing the best-fit affine transformation strain and nonaffine displacement associated with each particle between two stages of deformation. The nonaffine displacement exhibits exponential crossover from ballistic to diffusive behavior with respect to the cumulative deviatoric strain, indicating that in athermal granular packings, the cumulative deviatoric strain plays the role of time in thermal systems and drives effective particle diffusion. We further study the size-dependent deformation of the granular pillars by simulation, and find that different-sized pillars follow self-similar shape evolution during deformation. In addition, the yield stress of the pillars increases linearly with pillar size. Formation of transient shear lines in the pillars during deformation becomes more evident as pillar size increases. The width of these elementary shear bands is about twice the diameter of a particle, and does not vary with pillar size.

Macroscopic force experienced by extended objects in granular flows over a very broad Froude-number range : Macroscopic granular force on extended object

This paper revisits a great number of data from previous studies about the macroscopic force experienced by either objects moving at constant speed and depth inside static granular materials or motionless objects subject to steady granular flows. It focuses on extended objects whose immersed height is equal or close to the thickness of the surrounding granular medium. A simple scaling argument allows demarcating quasi-static from speed-squared force contributions for all the data from different geometries over a very broad range of Froude number. However, a wide scatter of the data is observed in the quasi-static regime. In the first step, a mean-field model is proposed to describe the average force. Mass and momentum balances are applied to a control volume, namely the expected volume of grains disturbed by the object, which is assumed to extend across the whole width and the entire height of the granular system. This allows defining an equivalent length scale which is computed by fitting the force predicted by the model to the available force data. In the second step, a circular shape is assumed for the effective mobilized domain and the associated diameter can be directly extracted from the computed equivalent length scale. This effective diameter is found to vary linearly with both the object width and the thickness of the granular layer moving around the extended object or the immersed depth of the object. The scaling highlights the key role played by the geometry which may enhance the force in the quasi-static regime.

Average balance equations, scale dependence, and energy cascade for granular materials

A new averaging method linking discrete to continuum variables of granular materials is developed and used to derive average balance equations. Its novelty lies in the choice of the decomposition between mean values and fluctuations of properties which takes into account the effect of gradients. Thanks to a local homogeneity hypothesis, whose validity is discussed, simplified balance equations are obtained. This original approach solves the problem of dependence of some variables on the size of the averaging domain obtained in previous approaches which can lead to huge relative errors (several hundred percentages). It also clearly separates affine and nonaffine fields in the balance equations. The resulting energy cascade picture is discussed, with a particular focus on unidirectional steady and fully developed flows for which it appears that the contact terms are dissipated locally unlike the kinetic terms which contribute to a nonlocal balance. Application of the method is demonstrated in the determination of the macroscopic properties such as volume fraction, velocity, stress, and energy of a simple shear flow, where the discrete results are generated by means of discrete particle simulation.

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.

Origin of a depth-independent drag force induced by stirring in granular media

Experiments have shown that when a horizontal cylinder rotates around the vertical axis in a granular medium, the drag force in the stationary regime becomes independent of the depth, in contradiction with the frictional picture stipulating that the drag should be proportional to the hydrostatic pressure. The goal of this study is to understand the origin of this depth independence of the granular drag. Intensive numerical simulations using the discrete element method are performed giving access to the stress distribution in the packing during the rotation of the cylinder. It is shown that the rotation induces a strong anisotropy in the stress distribution, leading to the formation of arches that screen the hydrostatic pressure in the vicinity of the cylinder and create a bubble of low pressure.

Precisely cyclic sand: self-organization of periodically sheared frictional grains

The disordered static structure and chaotic dynamics of frictional granular matter has occupied scientists for centuries, yet there are few organizational principles or guiding rules for this highly hysteretic, dissipative material. We show that cyclic shear of a granular material leads to dynamic self-organization into several phases with different spatial and temporal order. Using numerical simulations, we present a phase diagram in strain-friction space that shows chaotic dispersion, crystal formation, vortex patterns, and most unusually a disordered phase in which each particle precisely retraces its unique path. However, the system is not reversible. Rather, the trajectory of each particle, and the entire frictional, many-degrees-of-freedom system, organizes itself into a limit cycle absorbing state. Of particular note is that fact that the cyclic states are spatially disordered, whereas the ordered states are chaotic.

Non-monotonic dependence of the friction coefficient on heterogeneous stiffness

The complexity of the frictional dynamics at the microscopic scale makes difficult to identify all of its controlling parameters. Indeed, experiments on sheared elastic bodies have shown that the static friction coefficient depends on loading conditions, the real area of contact along the interfaces and the confining pressure. Here we show, by means of numerical simulations of a 2D Burridge-Knopoff model with a simple local friction law, that the macroscopic friction coefficient depends non-monotonically on the bulk elasticity of the system. This occurs because elastic constants control the geometrical features of the rupture fronts during the stick-slip dynamics, leading to four different ordering regimes characterized by different orientations of the rupture fronts with respect to the external shear direction. We rationalize these results by means of an energetic balance argument.

Spheronization process particle kinematics determined by discrete element simulations and particle image velocimentry measurements

Spheronization is an important pharmaceutical manufacturing technique to produce spherical agglomerates of 0.5-2mm diameter. These pellets have a narrow size distribution and a spherical shape. During the spheronization process, the extruded cylindrical strands break in short cylinders and evolve from a cylindrical to a spherical state by deformation and attrition/agglomeration mechanisms. Using the discrete element method, an integrated modeling-experimental framework is presented, that captures the particle motion during the spheronization process. Simulations were directly compared and validated against particle image velocimetry (PIV) experiments with monodisperse spherical and dry γ-Al2O3 particles.

RESULT

demonstrate a characteristic torus like flow pattern, with particle velocities about three times slower than the rotation speed of the friction plate. Five characteristic zones controlling the spheronization process are identified: Zone I, where particles undergo shear forces that favors attrition and contributes material to the agglomeration process; Zone II, where the static wall contributes to the mass exchange between particles; Zone III, where gravitational forces combined with particle motion induce particles to collide with the moving plate and re-enter Zone I; Zone IV, where a subpopulation of particles are ejected into the air when in contact with the friction plate structure; and Zone V where the low poloidal velocity favors a stagnant particle population and is entirely controlled by the batch size. These new insights in to the particle motion are leading to deeper process understanding, e.g., the effect of load and rotation speed to the pellet formation kinetics. This could be beneficial for the optimization of a manufacturing process as well as for the development of new formulations.

Attractive particle interaction forces and packing density of fine glass powders

We study the packing of fine glass powders of mean particle diameter in the range (4-52) μm both experimentally and by numerical DEM simulations. We obtain quantitative agreement between the experimental and numerical results, if both types of attractive forces of particle interaction, adhesion and non-bonded van der Waals forces are taken into account. Our results suggest that considering only viscoelastic and adhesive forces in DEM simulations may lead to incorrect numerical predictions of the behavior of fine powders. Based on the results from simulations and experiments, we propose a mathematical expression to estimate the packing fraction of fine polydisperse powders as a function of the average particle size.

Impact of the timestep in some molecular dynamics simulations on compression of granular systems

We conduct two-dimensional molecular dynamics simulations to study the statistical distribution of the force-moment (defined as stress multiplied by volume) of static granular packings under external isotropic compression. To that end, we generate packings by compressing initially ordered lattices using irregular, randomly generated, walls. Velocity-Verlet algorithm and linear spring-dashpot interactions are employed. With this specific method, the obtained statistical distributions of the force-moment are similar for different initial packings. However they depend on the timestep selection within a range of values. This shows that inadequate molecular dynamic simulations may provide different classes of solutions for the same physical process, and this could cause problems to validate theoretical approaches based on statistical mechanics.

Internal structure of inertial granular flows

We analyze inertial granular flows and show that, for all values of the inertial number I, the effective friction coefficient μ arises from three different parameters pertaining to the contact network and force transmission: (1) contact anisotropy, (2) force chain anisotropy, and (3) friction mobilization. Our extensive 3D numerical simulations reveal that μ increases with I mainly due to an increasing contact anisotropy and partially by friction mobilization whereas the anisotropy of force chains declines as a result of the destabilizing effect of particle inertia. The contact network undergoes topological transitions, and beyond I≃0.1 the force chains break into clusters immersed in a background "soup" of floating particles. We show that this transition coincides with the divergence of the size of fluidized zones characterized from the local environments of floating particles and a slower increase of μ with I.

Homogeneous steady state of a confined granular gas

The nonequilibrium statistical mechanics and kinetic theory for a model of a confined quasi-two-dimensional gas of inelastic hard spheres is presented. The dynamics of the particles includes an effective mechanism to transfer the energy injected in the vertical direction to the horizontal degrees of freedom. The Enskog approximation is formulated and used as the basis to investigate the temperature and the distribution function of the steady state eventually reached by the system. An exact scaling of the distribution function of the system having implications on the form of its moments is pointed out. The theoretical predictions are compared with numerical results obtained by a particle simulation method, and a good agreement is found.

Evolution of the force distributions in jammed packings of soft particles

The evolution of the force distributions during the isotropic compression of two-dimensional packings of soft frictional particles is investigated numerically. Regardless of the applied deformation, the normal contact force distribution P(f(n)) can be fitted by the product of a power law, and a stretched exponential, while the tangential force distribution P(f(t)) is fitted well by a Gaussian. With increasing strain, the asymptotic behavior at large forces does not change, but both P(f(n)) and P(f(t)) exhibit a broadening, even though, when scaled with the average forces, their widths decrease. Furthermore, the distribution of friction mobilization P(η) is a decreasing function of η=|f(t)|/(μf(n)), except for an increased probability of fully mobilized contacts (η=1). The excess coordination number of the packings increases with the applied strain, indicating that the more a packing is compressed the more stable it becomes.

The influence of rolling friction on the shear behaviour of non-cohesive pharmaceutical granules--an experimental and numerical investigation

Granule shear behaviour was investigated experimentally and numerically to evaluate the reliability of the numerical model. Additionally, parameters affecting the ensuing flow regimes - elastic quasi-static and inertial non-collisional - were highlighted. Furthermore, the influence of using the Lees-Edwards periodic boundary conditions or the standard boundary conditions was studied. Experiments were performed with microcrystalline cellulose granules of three size distributions using the FT4 powder rheometer. The numerical parameters, particle size, effective density, and particle stiffness were selected to match the experimental conditions. Experimentally, an unexpected particle size effect was evident where the resistance to shear increased with particle size. Numerically, combining rolling friction and increased shear rate enabled a transition from the inertial non-collisional to the elastic quasi-static regime at a reduced sliding friction coefficient. Presumably, this is an effect of increased particle overlap creating stronger contacts and facilitating force chain formation. Both boundary conditions provided comparable results provided a correction of system size was made, where larger systems were required for the standard boundary conditions. A satisfactory qualitative agreement between the experimentally and numerically determined yield loci emphasised the predictive capacity of the DEM. Rolling friction was in addition concluded to be an essential model parameter for obtaining an improved quantitative agreement.

Shallow granular flows down flat frictional channels: steady flows and longitudinal vortices

Granular flows down inclined channels with smooth boundaries are common in nature and industry. Nevertheless, flat boundaries have been much less investigated than bumpy ones, which are used by most experimental and numerical studies to avoid sliding effects. Using numerical simulations of each grain and of the side walls we recover quantitatively experimental results. At larger angles we predict a rich behavior, including granular convection and inverted density profiles suggesting a Rayleigh-Bénard type of instability. In many aspects flows on a flat base can be seen as flows over an effective bumpy base made of the basal rolling layer, giving Bagnold-type profiles in the overburden. We have tested a simple viscoplastic rheological model [Nature (London) 441, 727 (2006)] in average form. The transition between the unidirectional and the convective flows is then clearly apparent as a discontinuity in the constitutive relation.

Numerical Simulation of Particle Flow Motion in a Two-Dimensional Modular Pebble-Bed Reactor with Discrete Element Method

Modular pebble-bed nuclear reactor (MPBNR) technology is promising due to its attractive features such as high fuel performance and inherent safety. Particle motion of fuel and graphite pebbles is highly associated with the performance of pebbled-bed modular nuclear reactor. To understand the mechanism of pebble’s motion in the reactor, we numerically studied the influence of number ratio of fuel and graphite pebbles, funnel angle of the reactor, height of guide ring on the distribution of pebble position, and velocity by means of discrete element method (DEM) in a two-dimensional MPBNR. Velocity distributions at different areas of the reactor as well as mixing characteristics of fuel and graphite pebbles were investigated. Both fuel and graphite pebbles moved downward, and a uniform motion was formed in the column zone, while pebbles motion in the cone zone was accelerated due to the decrease of the cross sectional flow area. The number ratio of fuel and graphite pebbles and the height of guide ring had a minor influence on the velocity distribution of pebbles, while the variation of funnel angle had an obvious impact on the velocity distribution. Simulated results agreed well with the work in the literature.

Frictional dependence of shallow-granular flows from discrete particle simulations

A shallow-layer model for granular flows down inclines is completed with a closure relation for the macroscopic bed friction obtained from micro-scale, discrete particle simulations of steady flows over geometrically rough bases with contact friction. Microscopic friction can be different between bulk particles and with particles at the base, where the latter is systematically varied. When extending the known friction closure relation to be a function of both bulk flow and bed properties, surprisingly, we find that the macroscopic bed friction is only weakly dependent on the contact friction of the bed particles and is predominantly determined by the properties of the flowing particles. Implications for constitutive modelling and possible experiments to better understand the bulk rheology are discussed.

Shear-induced alignment and dynamics of elongated granular particles

The alignment, ordering, and rotation of elongated granular particles was studied in shear flow. The time evolution of the orientation of a large number of particles was monitored in laboratory experiments by particle tracking using optical imaging and x-ray computed tomography. The experiments were complemented by discrete element simulations. The particles develop an orientational order. In the steady state the time- and ensemble-averaged direction of the main axis of the particles encloses a small angle with the streamlines. This shear alignment angle is independent of the applied shear rate, and it decreases with increasing grain aspect ratio. At the grain level the steady state is characterized by a net rotation of the particles, as dictated by the shear flow. The distribution of particle rotational velocities was measured both in the steady state and also during the initial transients. The average rotation speed of particles with their long axis perpendicular to the shear alignment angle is larger, while shear aligned particles rotate slower. The ratio of this fast/slow rotation increases with particle aspect ratio. During the initial transient starting from an unaligned initial condition, particles having an orientation just beyond the shear alignment angle rotate opposite to the direction dictated by the shear flow.

Friction and the oscillatory motion of granular flows

This contribution reports on numerical simulations of two-dimensional granular flows on erodible beds. The broad aim is to investigate whether simple flows of model granular matter exhibit spontaneous oscillatory motion in generic flow conditions, and in this case, whether the frictional properties of the contacts between grains may affect the existence or the characteristics of this oscillatory motion. The analysis of different series of simulations shows that the flow develops an oscillatory motion with a well-defined frequency which increases like the inverse of the velocity's square root. We show that the oscillation is essentially a surface phenomenon. The amplitude of the oscillation is higher for lower volume fractions and can thus be related to the flow velocity and grains' friction properties. The study of the influence of the periodic geometry of the simulation cell shows no significant effect. These results are discussed in relation to sonic sands.

Discrete simulation of dense flows of polyhedral grains down a rough inclined plane

The influence of grain angularity on the properties of dense flows down a rough inclined plane are investigated. Three-dimensional numerical simulations using the nonsmooth contact dynamics method are carried out with both spherical (rounded) and polyhedral (angular) grain assemblies. Both sphere and polyhedra assemblies abide by the flow start and stop laws, although much higher tilt angle values are required to trigger polyhedral grain flow. In the dense permanent flow regime, both systems show similarities in the bulk of the material (away from the top free surface and the substrate), such as uniform values of the solid fraction, inertial number and coordination number, or linear dependency of the solid fraction and effective friction coefficient with the inertial number. However, discrepancies are also observed between spherical and polyhedral particle flows. A dead (or nearly arrested) zone appears in polyhedral grain flows close to the rough bottom surface, reflected by locally concave velocity profiles, locally larger coordination number and solid fraction values, smaller inertial number values. This dead zone disappears for smooth bottom surfaces. In addition, unlike sphere assemblies, polyhedral grain assemblies exhibit significant normal stress differences, which increase close to the substrate.