A 2-D asymmetric exclusion model for granular flows

A 2-D version of the asymmetric exclusion model for granular sheared flows is presented. The velocity profile exhibits two qualitatively different behaviors, dependent on control parameters. For low friction, the velocity profile follows an exponential decay while for large friction the profile is more accurately represented by a Gaussian law. The phase transition occurring between these two behavior is identified by the appearance of correlations in the cluster size distribution. Finally, a mean--field theory gives qualitative and quantitative good agreement with the numerical results.Comment: 13 pages, 5 figures; typos added, one definition change

Pre-avalanche instabilities in a granular pile

We investigate numerically the transition between static equilibrium and dynamic surface flow of a 2D cohesionless granular system driven by a continuous gravity loading. This transition is characterized by intermittent local dynamic rearrangements and can be described by an order parameter defined as the density of critical contacts, e.g. contacts where the friction is fully mobilized. Analysis of the spatial correlations of critical contacts shows the occurence of ``fluidized'' clusters which exhibit a power-law divergence in size at the approach of the stability limit. The results are compatible with recent models that describe the granular system during the static/dynamic transition as a multi-phase system.Comment: 9 pages, 6 figures, submitted to Phys. Rev. Let

Statistical Mechanics of Stress Transmission in Disordered Granular Arrays

We give a statistical-mechanical theory of stress transmission in disordered arrays of rigid grains with perfect friction. Starting from the equations of microscopic force and torque balance we derive the fundamental equations of stress equilibrium. We illustrate the validity of our approach by solving the stress distribution of a homogeneous and isotropic array.Comment: 4 pages, to be published in PR

Creep motion in a granular pile exhibiting steady surface flow

We investigate experimentally granular piles exhibiting steady surface flow. Below the surface flow, it has been believed exisitence of a `frozen' bulk region, but our results show absence of such a frozen bulk. We report here that even the particles in deep layers in the bulk exhibit very slow flow and that such motion can be detected at an arbitrary depth. The mean velocity of the creep motion decays exponentially with depth, and the characteristic decay length is approximately equal to the particle-size and independent of the flow rate. It is expected that the creep motion we have seeen is observable in all sheared granular systems.Comment: 3 pages, 4 figure

Strain versus stress in a model granular material: a Devil's staircase

The series of equilibrium states reached by disordered packings of rigid, frictionless discs in two dimensions, under gradually varying stress, are studied by numerical simulations. Statistical properties of trajectories in configuration space are found to be independent of specific assumptions ruling granular dynamics, and determined by geometry only. A monotonic increase in some macroscopic loading parameter causes a discrete sequence of rearrangements. For a biaxial compression, we show that, due to the statistical importance of such events of large magnitudes, the dependence of the resulting strain on stress direction is a Levy flight in the thermodynamic limit.Comment: REVTeX, 4 pages, 5 included PostScript figures. New version altered throughout text, very close to published pape

Cooperativity in sandpiles: statistics of bridge geometries

Bridges form dynamically in granular media as a result of spatiotemporal inhomogeneities. We classify bridges as linear and complex, and analyse their geometrical characteristics. In particular, we find that the length distribution of linear bridges is exponential. We then turn to the analysis of the orientational distribution of linear bridges and find that, in three dimensions, they are {\it vertically diffusive but horizontally superdiffusive}; thus, when they exist, long linear bridges form `domes'. Our results are in good accord with Monte Carlo simulations of bridge structure; we make predictions for quantities that are experimentally accessible, and suggest that bridges are very closely related to force chains.Comment: 15 pages, 10 figures. Minor changes and update

Granular Rheology in Zero Gravity

We present an experimental investigation on the rheological behavior of model granular media made of nearly elastic spherical particles. The experiments are performed in a cylindrical Couette geometry and the experimental device is placed inside an airplane undergoing parabolic flights to cancel the effect of gravity. The corresponding curves, shear stress versus shear rate, are presented and a comparison with existing theories is proposed. The quadratic dependence on the shear rate is clearly shown and the behavior as a function of the solid volume fraction of particles exhibits a power law function. It is shown that theoretical predictions overestimate the experiments. We observe, at intermediate volume fractions, the formation of rings of particles regularly spaced along the height of the cell. The differences observed between experimental results and theoretical predictions are discussed and related to the structures formed in the granular medium submitted to the external shear.Comment: 10 pages, 6 figures to be published in Journal of Physics : Condensed Matte

Random Packings of Frictionless Particles

We study random packings of frictionless particles at T=0. The packing fraction where the pressure becomes nonzero is the same as the jamming threshold, where the static shear modulus becomes nonzero. The distribution of threshold packing fractions narrows and its peak approaches random close-packing as the system size increases. For packing fractions within the peak, there is no self-averaging, leading to exponential decay of the interparticle force distribution.Comment: 4 pages, 3 figure

Force Distribution in a Granular Medium

We report on systematic measurements of the distribution of normal forces exerted by granular material under uniaxial compression onto the interior surfaces of a confining vessel. Our experiments on three-dimensional, random packings of monodisperse glass beads show that this distribution is nearly uniform for forces below the mean force and decays exponentially for forces greater than the mean. The shape of the distribution and the value of the exponential decay constant are unaffected by changes in the system preparation history or in the boundary conditions. An empirical functional form for the distribution is proposed that provides an excellent fit over the whole force range measured and is also consistent with recent computer simulation data.Comment: 6 pages. For more information, see http://mrsec.uchicago.edu/granula

Density of states in random lattices with translational invariance

We propose a random matrix approach to describe vibrational excitations in disordered systems. The dynamical matrix M is taken in the form M=AA^T where A is some real (not generally symmetric) random matrix. It guaranties that M is a positive definite matrix which is necessary for mechanical stability of the system. We built matrix A on a simple cubic lattice with translational invariance and interaction between nearest neighbors. We found that for certain type of disorder phonons cannot propagate through the lattice and the density of states g(w) is a constant at small w. The reason is a breakdown of affine assumptions and inapplicability of the elasticity theory. Young modulus goes to zero in the thermodynamic limit. It strongly reminds of the properties of a granular matter at the jamming transition point. Most of the vibrations are delocalized and similar to diffusons introduced by Allen, Feldman et al., Phil. Mag. B v.79, 1715 (1999).Comment: 4 pages, 5 figure