Density inversion in rapid granular flows: the supported regime

This paper presents numerical findings on rapid 2D and 3D granular flows on a bumpy base. In the supported regime studied here, a strongly sheared, dilute and agitated layer spontaneously appears at the base of the flow and supports a compact packing of grains moving as a whole. In this regime, the flow behaves like a sliding block on the bumpy base. In particular, for flows on a horizontal base, the average velocity decreases linearly in time and the average kinetic energy decreases linearly with the travelled distance, those features being characteristic of solid-like friction. This allows us to define and measure an effective friction coefficient, which is independent of the mass and velocity of the flow. This coefficient only loosely depends on the value of the micromechanical friction coefficient whereas the infuence of the bumpiness of the base is strong. We give evidence that this dilute and agitated layer does not result in significantly less friction. Finally, we show that a steady regime of supported flows can exist on inclines whose angle is carefully chosen.     

Cluster formation, pressure and density measurements in a granular medium fluidized by vibrations

We report experimental results on the behavior of an ensemble of inelastically colliding particles, excited by a vibrated piston in a vertical cylinder. When the particle number is increased, we observe a transition from a regime where the particles have erratic motions (“granular gas”) to a collective behavior where all the particles bounce like a nearly solid body. In the gas-like regime, we measure the density of particles as a function of the altitude and the pressure as a function of the number N of particles. The atmosphere is found to be exponential far enough from the piston, and the “granular temperature”, T, dependence on the piston velocity, V, is of the form , where is a decreasing function of N. This may explain previous conflicting numerical results. Received 1 February 1999     

A quasi-elastic regime for vibrated granular gases

Using simple scaling arguments and two-dimensional numerical simulations of a granular gas excited by vibrating one of the container boundaries, we study a double limit of small 1-r and large L, where r is the restitution coefficient and L the size of the container. We show that if the particle density n₀ and (1-r²)(n₀ Ld) where d is the particle diameter, are kept constant and small enough, the granular temperature, i.e. the mean value of the kinetic energy per particle, 〈E 〉/N, tends to a constant whereas the mean dissipated power per particle, 〈D 〉/N, decreases like when N increases, provided that (1-r²)(n₀ Ld)² < 1. The relative fluctuations of E, D and the power injected by the moving boundary, I, have simple properties in that regime. In addition, the granular temperature can be determined from the fluctuations of the power I(t) injected by the moving boundary.     

Correlations and aggregate statistics in granular packs

We study how the aggregate statistical properties for density fluctuations in granular aggregates scale with the sample size and how such a scaling is associated with the correlations between grains. Correlations are studied both between grain positions and between Vorono? cell volumes, showing distinct behaviors and properties. A non-linear scaling in the aggregate volume fluctuations as function of the sample size is discovered and the connection between such anomalous scaling and correlations is explained. It emerges that volume fluctuations might be described by means of a single universal equation for all samples at all cluster sizes.     

Improving mixture of grains by using bi-dimensional Galton boards

This paper presents an experimental and numerical study that deals with the problem of mixing grains falling down through a bi-dimensional Galton board (BGB). The special issue addressed here is the influence of the presence of lateral walls in the BGB. Disks of equal diameters but different species are launched from the top of the device. During the fall, disks collide with obstacles (arranged to form a triangular lattice) and with the lateral walls. The exit distribution of particles at the bottom of the board is determined and the incidence of the presence of walls in the mixing quality is studied as a function of W the relative separation between lateral walls. Two types of indexes are evaluated to characterize the efficiency in the obtained mixture. The presence of walls has proven to be crucial to enhance the quality of the mixture of particles.     

Strength and failure of cemented granular matter

Cemented granular materials (CGMs) consist of densely packed solid particles and a pore-filling solid matrix sticking to the particles. We use a sub-particle lattice discretization method to investigate the particle-scale origins of strength and failure properties of CGMs. We show that jamming of the particles leads to highly inhomogeneous stress fields. The stress probability density functions are increasingly wider for a decreasing matrix volume fraction, the stresses being more and more concentrated in the interparticle contact zones with an exponential distribution as in cohesionless granular media. Under uniaxial loading, pronounced asymmetry can occur between tension and compression both in strength and in the initial stiffness as a result of the presence of bare contacts (with no matrix interposed) between the particles. Damage growth is analyzed by considering the evolution of stiffness degradation and the number of broken bonds in the particle phase. A brutal degradation appears in tension as a consequence of brittle fracture in contrast to the more progressive nature of damage growth in compression. We also carry out a detailed parametric study in order to assess the combined influence of the matrix volume fraction and particle-matrix adherence. Three regimes of crack propagation can be distinguished corresponding to no particle damage, particle abrasion and particle fragmentation, respectively. We find that particle damage scales well with the relative toughness of the particle-matrix interface with respect to the particle toughness. This relative toughness is a function of both matrix volume fraction and particle-matrix adherence and it appears therefore to be the unique control parameter governing transition from soft to hard behavior.     

Ripple formation in weakly turbulent flow

The formation of granular ripples under liquid shear flow in an annular channel is studied experimentally. The erodible granular bed is subject to weakly turbulent flows without a defined sharp boundary layer close to the granular bed. The flow field and the degree of turbulence is characterized quantitatively by using a particle image velocimeter and a laser-Doppler velocimeter, respectively. A new range of particle Reynolds numbers at the lower limit of the Shields diagram were explored. Quantitative measurements of the granular flow on the surface reveal that the threshold for particle motion coincides within the order of one percent with the threshold for ripple formation. In fully developed ripples it was found that on the leeward side of the ripples regions of low-velocity gradients exist where granular motion is scarce, indicating that the coupling between the ripples is mainly caused by the flow field of the liquid.     

Cellular model for the compaction of a vertically tapped granular column

A cellular model for the compaction of granular material is described. It takes into account horizontal redistribution as well as vertical transfer of particles. Parameters are the width of the horizontal redistribution and the settling probability. Numerical simulations of the behaviour of a granular column in a container are shown as an example, and the evolution of some characteristic features over time has been followed for some typical configurations. Experimental results for the time evolution of the density can be reproduced for a settling probability proportional to the unoccupied spaces for particles in the lower cells. Received 3 March 2000 and Received in final form 21 July 2000     

Upward penetration of grains through a granular medium

We study experimentally the creeping penetration of guest (percolating) grains through densely packed granular media in two dimensions. The evolution of the system of the guest grains during the penetration is studied by image analysis. To quantify the changes in the internal structure of the packing, we use Voronoï tessellation and a certain shape factor which is a clear indicator of the presence of different underlying substructures (domains). We first consider the impact of the effective gravitational acceleration on upward penetration of grains. It is found that the higher effective gravity increases the resistance to upward penetration and enhances structural organization in the system of the percolating grains. We also focus our attention on the dependence of the structural rearrangements of percolating grains on some parameters like polydispersity and the initial packing fraction of the host granular system. It is found that the anisotropy of penetration is larger in the monodisperse case than in the bidisperse one, for the same value of the packing fraction of the host medium. Compaction of initial host granular packing also increases anisotropy of penetration of guest grains. When a binary mixture of large and small guest grains is penetrated into the host granular medium, we observe size segregation patterns.     

Granular flows in a rotating drum: the scaling law between velocity and thickness of the flow

The flow of dry granular material in a half-filled rotating drum is studied. The thickness of the flowing zone is measured for several rotation speeds, drum sizes and beads sizes (size ratio between drum and beads ranging from 47 to 7400). Varying the rotation speed, a scaling law linking mean velocity vs. thickness of the flow, v∼h^m, is deduced for each couple (beads, drum). The obtained exponent m is not always equal to 1, the value previously reported for a drum in litterature, but varies with the geometry of the system. For small size ratios, exponents higher than 1 are obtained due to a saturation of the flowing zone thickness. The exponent of the power law decreases with the size ratio, leading to exponents lower than 1 for high size ratios. These exponents imply that the velocity gradient of a dry granular flow in a rotating drum is not constant. More fundamentally, these results show that the flow of a granular material in a rotating drum is very sensible to the geometry, and that the deduction of the “rheology” of a granular medium flowing in such a geometry is not obvious.     

Magneto-vibratory separation of glass and bronze granular mixtures immersed in a paramagnetic liquid

A fluid-immersed granular mixture may spontaneously separate when subjected to vertical vibration, separation occurring when the ratio of particle inertia to fluid drag is sufficiently different between the component species of the mixture. Here, we describe how fluid-driven separation is influenced by magneto-Archimedes buoyancy, the additional buoyancy force experienced by a body immersed in a paramagnetic fluid when a strong inhomogeneous magnetic field is applied. In our experiments glass and bronze mixtures immersed in paramagnetic aqueous solutions of MnCl₂ have been subjected to sinusoidal vertical vibration. In the absence of a magnetic field the separation is similar to that observed when the interstitial fluid is water. However, at modest applied magnetic fields, magneto-Archimedes buoyancy may balance the inertia/fluid-drag separation mechanism, or it may dominate the separation process. We identify the vibratory and magnetic conditions for four granular configurations, each having distinctive granular convection. Abrupt transitions between these states occur at well-defined values of the magnetic and vibrational parameters. In order to gain insight into the dynamics of the separation process we use computer simulations based on solutions of the Navier-Stokes' equations. The simulations reproduce the experimental results revealing the important role of convection and gap formation in the stability of the different states.     

Bidisperse granular avalanches on inclined planes: A rich variety of behaviors

Experiments were performed to provide insight into the flow behavior and structure of bimodal mixtures of grains in gravity-driven, free-surface flows. Unsteady unconfined flows were produced by releasing instantaneously a dry granular mass, composed of two particle sizes, over a rough inclined plane. As a result of size segregation, the small particles are found at the bottom of the flow and final deposit, the large particles are found at the free surface, but also on the lateral borders and at the front of the flow. The lateral and vertical inhomogeneous repartitions of particles lead to two main effects that are completely absent in monodispersed flows. The outline effect results from the accumulation of large beads on the periphery of the flow depending on the value of the relative friction of each particle species on the plane. This effect in turn causes a narrowing of the flow and/or an increase of length of the final deposit. The interface effect results of the interaction between layers of different size particles and causes the modification of the thickness of the deposit. These effects occur simultaneously and their combination leads to a great variety of behaviors. In this investigation, evidence of the diversity of behaviors is presented as the size ratio, relative friction and concentration of each particle species are varied.     

Monodisperse dry granular flows on inclined planes: Role of roughness

Recent studies have pointed out the importance of the basal friction on the dynamics of granular flows. We present experimental results on the influence of the roughness of the inclined plane on the dynamics of a monodisperse dry granular flow. We found experimentally that there exists a maximum of the friction for a given relative roughness. This maximum is shown to be independent of the slope angle. This behavior is observed for four planes with different bump sizes (given by the size of the beads glued on the plane) from 200 m to 2 mm. The relative roughness corresponding to the maximum of the friction can be predicted with a geometrical model of stability of one single bead on the plane. The main parameters are the size of the bumps and the size of the flowing beads. In order to obtain a higher precision, the model also takes into account the spacing between the bumps of the rough plane. Experimental results and model are in good agreement for all the planes we studied. Other parameters, like the sphericity of the beads, or irregularities in the thickness of the layer of glued particles, are shown to be of influence on the friction.     

Stress fluctuations and macroscopic stick-slip in granular materials

This paper deals with the quasi-static regime of deformation of granular matter. It investigates the size of the Representative Elementary Volume (REV), which is the minimum packing size above which the macroscopic mechanical behaviour of granular materials can be defined from averaging. The first part uses typical results from recent literature and finds that the minimum REV contains in general 10 grains; this result holds true either for most experiments or for Discrete Element Method (DEM) simulation. This appears to be quite small. However, the second part gives a counterexample, which has been found when investigating uniaxial compression of glass spheres which exhibit stick-slip; we show in this case that the minimum REV becomes 10⁷ grains. This makes the system not computable by DEM. Moreover, similarity between the Richter law of seism and the exponential statistics of stick-slip is stressed. Received 8 March 2002 and Received in final form 12 July 2002     

Effective temperature and jamming transition in dense, gently sheared granular assemblies

We present extensive computational results for the effective temperature, defined by the fluctuation-dissipation relation between the mean square displacement and the average displacement of grains, under the action of a weak, external perturbation, of a sheared, bi-disperse granular packing of compressible spheres. We study the dependence of this parameter on the shear rate and volume fractions, the type of particle and the observable in the fluctuation-dissipation relation. We find the same temperature for different tracer particles in the system. The temperature becomes independent on the shear rate for slow enough shear suggesting that it is the effective temperature of the jammed packing. However, we also show that the agreement of the effective temperature for different observables is only approximate, for very long times, suggesting that this defintion may not capture the full thermodynamics of the system. On the other hand, we find good agreement between the dynamical effective temperature and a compactivity calculated assuming that all jammed states are equiprobable. Therefore, this definition of temperature may capture an instance of the ergodic hypothesis for granular materials as proposed by theoretical formalisms for jamming. Finally, our simulations indicate that the average shear stress and apparent shear viscosity follow the usual relation with the shear rate for complex fluids. Our results show that the application of shear induces jamming in packings whose particles interact by tangential forces.     

Gas-driven subharmonic waves in a vibrated two-phase granular material

Vibrated powders exhibit striking phenomena: subharmonic waves, oscillons, convection, heaping, and even bubbling. We demonstrate novel rectangular profile subharmonic waves for vibrated granular material, that occur uniquely in the two-phase case of grains, and a fluid, such as air. These waves differ substantially from those for the gas-free case, exhibit different dispersion relations, and occur for specific shaking parameters and air pressure, understandable with gas-particle flow models. These waves occur when the gas diffusively penetrates the granular layer in a time comparable to the shaker period. As the pressure is lowered towards P =0, the granular-gas system exhibits a Knudsen regime. This instability provides an opportunity to quantitatively test models of two-phase flow.     

Granular flow from a silo: Discrete-particle simulations in three dimensions

Molecular (or granular) dynamics methods are used to study the gravity-driven flow of granular material through a horizontal aperture in three dimensions. The grains are spherical and modeled using a short-range repulsive interaction, together with normal and tangential frictional damping forces. The material is contained in a rough-walled cylindrical container with a circular hole in its base, and to permit flow measurements under steady-state conditions a continuous feed approach is employed in which exiting grains are replaced at the upper surface of the material. The dependence of flow velocity and discharge rate on aperture diameter is found to agree with experiment; other quantities such as the kinetic energy and pressure distributions are also examined. Received 5 June 2000 and Received in final form 21 September 2000     

Instabilities in vertically vibrated fluid-grain systems

When a bed of fluid-immersed fine grains is exposed to vertical vibration a wealth of phenomena may be observed. At low frequencies a horizontal bed geometry is generally unstable and the bed breaks spatial symmetry, acquiring a tilt. At the same time it undergoes asymmetric granular convection. Fine binary mixtures may separate completely into layers or patterns of stripes. The separated regions may exhibit instabilities in which they undergo wave-like motion or exhibit quasi-periodic oscillations. We briefly review these and a number of related behaviours, identifying the physical mechanisms behind each. Finally, we discuss the magneto-vibratory separation of binary mixtures which results from exposing each component to a different effective gravity and describe the influence of a background fluid on this process.