Velocity correlations in dense granular flows

Velocity fluctuations of grains flowing down a rough inclined plane are experimentally studied. The grains at the free surface exhibit fluctuating motions, which are correlated over few grains diameters. The characteristic correlation length is shown to depend on the inclination of the plane and not on the thickness of the flowing layer. This result strongly supports the idea that dense granular flows are controlled by a characteristic length larger than the particle diameter

Rapid granular flows on a rough incline: phase diagram, gas transition, and effects of air drag

We report experiments on the overall phase diagram of granular flows on an incline with emphasis on high inclination angles where the mean layer velocity approaches the terminal velocity of a single particle free falling in air. The granular flow was characterized by measurements of the surface velocity, the average layer height, and the mean density of the layer as functions of the hopper opening, the plane inclination angle and the downstream distance x of the flow. At high inclination angles the flow does not reach an x-invariant steady state over the length of the inclined plane. For low volume flow rates, a transition was detected between dense and very dilute (gas) flow regimes. We show using a vacuum flow channel that air did not qualitatively change the phase diagram and did not quantitatively modify mean flow velocities of the granular layer except for small changes in the very dilute gas-like phase.Comment: 10 pages, 16 figures, accepted to Phys. Rev.

The S shape of a granular pile in a rotating drum

The shape of a granular pile in a rotating drum is investigated. Using Discrete Elements Method (DEM) simulations we show that the "S shape" obtained for high rotation speed can be accounted for by the friction on the end plates. A theoretical model which accounts for the effect of the end plates is presented and the equation of the shape of the free surface is derived. The model reveals a dimensionless number which quantifies the influence of the end plates on the shape of the pile. Finally, the scaling laws of the system are discussed and numerical results support our conclusions

Flow of wet granular materials

The transition from frictional to lubricated flow of a dense suspension of non-Brownian particles is studied. The pertinent parameter characterizing this transition is the Leighton number $Le = \frac{\eta_s \dot{\gamma}}{\sigma}$, which represents the ratio of lubrication to frictional forces. The Leighton number $Le$ defines a critical shear rate below which no steady flow without localization exists. In the frictional regime the shear flow is localized. The lubricated regime is not simply viscous: the ratio of shear to normal stresses remains constant, as in the frictional regime; moreover the velocity profile has a single universal form in both frictional and lubricated regimes. Finally, a discrepancy between local and global measurements of viscosity is identified, which suggests inhomogeneity of the material under flow.Comment: Accepted for publication by Physical Review Letters (december 2004

Scaling of the critical slip distance in granular layers

We investigate the nature of friction in granular layers by means of numerical simulation focusing on the critical slip distance, over which the system relaxes to a new stationary state. Analyzing a transient process in which the sliding velocity is instantaneously changed, we find that the critical slip distance is proportional to the sliding velocity. We thus define the relaxation time, which is independent of the sliding velocity. It is found that the relaxation time is proportional to the layer thickness and inversely proportional to the square root of the pressure. An evolution law for the relaxation process is proposed, which does not contain any length constants describing the surface geometry but the relaxation time of the bulk granular matter. As a result, the critical slip distance is scaled with a typical length scale of a system. It is proportional to the layer thickness in an instantaneous velocity change experiment, whereas it is scaled with the total slip distance in a spring-block system on granular layers.Comment: 4 papge

A constitutive law for dense granular flows

A continuum description of granular flows would be of considerable help in predicting natural geophysical hazards or in designing industrial processes. However, the constitutive equations for dry granular flows, which govern how the material moves under shear, are still a matter of debate. One difficulty is that grains can behave like a solid (in a sand pile), a liquid (when poured from a silo) or a gas (when strongly agitated). For the two extreme regimes, constitutive equations have been proposed based on kinetic theory for collisional rapid flows, and soil mechanics for slow plastic flows. However, the intermediate dense regime, where the granular material flows like a liquid, still lacks a unified view and has motivated many studies over the past decade. The main characteristics of granular liquids are: a yield criterion (a critical shear stress below which flow is not possible) and a complex dependence on shear rate when flowing. In this sense, granular matter shares similarities with classical visco-plastic fluids such as Bingham fluids. Here we propose a new constitutive relation for dense granular flows, inspired by this analogy and recent numerical and experimental work. We then test our three-dimensional (3D) model through experiments on granular flows on a pile between rough sidewalls, in which a complex 3D flow pattern develops. We show that, without any fitting parameter, the model gives quantitative predictions for the flow shape and velocity profiles. Our results support the idea that a simple visco-plastic approach can quantitatively capture granular flow properties, and could serve as a basic tool for modelling more complex flows in geophysical or industrial applications.Comment: http://www.nature.com/nature/journal/v441/n7094/abs/nature04801.htm

Force transmission in a packing of pentagonal particles

We perform a detailed analysis of the contact force network in a dense confined packing of pentagonal particles simulated by means of the contact dynamics method. The effect of particle shape is evidenced by comparing the data from pentagon packing and from a packing with identical characteristics except for the circular shape of the particles. A counterintuitive finding of this work is that, under steady shearing, the pentagon packing develops a lower structural anisotropy than the disk packing. We show that this weakness is compensated by a higher force anisotropy, leading to enhanced shear strength of the pentagon packing. We revisit "strong" and "weak" force networks in the pentagon packing, but our simulation data provide also evidence for a large class of "very weak" forces carried mainly by vertex-to-edge contacts. The strong force chains are mostly composed of edge-to-edge contacts with a marked zig-zag aspect and a decreasing exponential probability distribution as in a disk packing

A two-dimensional depth-averaged μ(I)-rheology for dense granular avalanches

Steady uniform granular chute flows are common in industry and provide an important test case for new theoretical models. This paper introduces depth-integrated viscous terms into the momentum-balance equations by extending the recent depth-averaged μ(I) -rheology for dense granular flows to two spatial dimensions, using the principle of material frame indifference or objectivity. Scaling the cross-slope coordinate on the width of the channel and the velocity on the one-dimensional steady uniform solution, we show that the steady two-dimensional downslope velocity profile is independent of scale. The only controlling parameters are the channel aspect ratio, the slope inclination angle and the frictional properties of the chute and the sidewalls. Solutions are constructed for both no-slip conditions and for a constant Coulomb friction at the walls. For narrow chutes, a pronounced parabolic-like depth-averaged downstream velocity profile develops. However, for very wide channels, the flow is almost uniform with narrow boundary layers close to the sidewalls. Both of these cases are in direct contrast to conventional inviscid avalanche models, which do not develop a cross-slope profile. Steady-state numerical solutions to the full three-dimensional μ(I) -rheology are computed using the finite element method. It is shown that these solutions are also independent of scale. For sufficiently shallow channels, the depth-averaged velocity profile computed from the full solution is in excellent agreement with the results of the depth-averaged theory. The full downstream velocity can be reconstructed from the depth-averaged theory by assuming a Bagnold-like velocity profile with depth. For wide chutes, this is very close to the results of the full three-dimensional calculation. For experimental validation, a laser profilometer and balance are used to determine the relationship between the total mass flux in the chute and the flow thickness for a range of slope angles and channel widths, and particle image velocimetry (PIV) is used to record the corresponding surface velocity profiles. The measured values are in good quantitative agreement with reconstructed solutions to the new depth-averaged theory

Internal states of model isotropic granular packings. I. Assembling process, geometry and contact networks

This is the first paper of a series of three, reporting on numerical simulation studies of geometric and mechanical properties of static assemblies of spherical beads under an isotropic pressure. Frictionless systems assemble in the unique random close packing (RCP) state in the low pressure limit if the compression process is fast enough, slower processes inducing traces of crystallization, and exhibit specific properties directly related to isostaticity of the force-carrying structure. The different structures of frictional packings assembled by various methods cannot be classified by the sole density. While lubricated systems approach RCP densities and coordination number z^*~=6 on the backbone in the rigid limit, an idealized "vibration" procedure results in equally dense configurations with z^*~=4.5. Near neighbor correlations on various scales are computed and compared to available laboratory data, although z^* values remain experimentally inaccessible. Low coordination packings have many rattlers (more than 10% of the grains carry no force), which should be accounted for on studying position correlations, and a small proportion of harmless "floppy modes" associated with divalent grains. Frictional packings, however slowly assembled under low pressure, retain a finite level of force indeterminacy, except in the limit of infinite friction.Comment: 29 pages. Published in Physical Review

Using PIV to measure granular temperature in saturated unsteady polydisperse granular flows

The motion of debris flows, gravity-driven fast moving mixtures of rock, soil and water can be interpreted using the theories developed to describe the shearing motion of highly concentrated granular fluid flows. Frictional, collisional and viscous stress transfer between particles and fluid characterizes the mechanics of debris flows. To quantify the influence of collisional stress transfer, kinetic models have been proposed. Collisions among particles result in random fluctuations in their velocity that can be represented by their granular temperature, T. In this paper particle image velocimetry, PIV, is used to measure the instantaneous velocity field found internally to a physical model of an unsteady debris flow created by using “transparent soil”—i.e. a mixture of graded glass particles and a refractively matched fluid. The ensemble possesses bulk properties similar to that of real soil-pore fluid mixtures, but has the advantage of giving optical access to the interior of the flow by use of plane laser induced fluorescence, PLIF. The relationship between PIV patch size and particle size distribution for the front and tail of the flows is examined in order to assess their influences on the measured granular temperature of the system. We find that while PIV can be used to ascertain values of granular temperature in dense granular flows, due to increasing spatial correlation with widening gradation, a technique proposed to infer the true granular temperature may be limited to flows of relatively uniform particle size or large bulk