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

Geometric origin of excess low-frequency vibrational modes in amorphous solids

Glasses have a large excess of low-frequency vibrational modes in comparison with crystalline solids. We show that such a feature is a necessary consequence of the geometry generic to weakly connected solids. In particular, we analyze the density of states of a recently simulated system, comprised of weakly compressed spheres at zero temperature. We account for the observed a) constancy of the density of modes with frequency, b) appearance of a low-frequency cutoff, and c) power-law increase of this cutoff with compression. We predict a length scale below which vibrations are very different from those of a continuous elastic body.Comment: 4 pages, 2 figures. Argument rewritten, identical result

Statistics of the contact network in frictional and frictionless granular packings

Simulated granular packings with different particle friction coefficient mu are examined. The distribution of the particle-particle and particle-wall normal and tangential contact forces P(f) are computed and compared with existing experimental data. Here f equivalent to F/F-bar is the contact force F normalized by the average value F-bar. P(f) exhibits exponential-like decay at large forces, a plateau/peak near f = 1, with additional features at forces smaller than the average that depend on mu. Computations of the force-force spatial distribution function and the contact point radial distribution function indicate that correlations between forces are only weakly dependent on friction and decay rapidly beyond approximately three particle diameters. Distributions of the particle-particle contact angles show that the contact network is not isotropic and only weakly dependent on friction. High force-bearing structures, or force chains, do not play a dominant role in these three dimensional, unloaded packings.Comment: 11 pages, 13 figures, submitted to PR

Gravity-driven Dense Granular Flows

We report and analyze the results of numerical studies of dense granular flows in two and three dimensions, using both linear damped springs and Hertzian force laws between particles. Chute flow generically produces a constant density profile that satisfies scaling relations suggestive of a Bagnold grain inertia regime. The type of force law has little impact on the behavior of the system. Bulk and surface flows differ in their failure criteria and flow rheology, as evidenced by the change in principal stress directions near the surface. Surface-only flows are not observed in this geometry.Comment: 4 pages, RevTeX 3.0, 4 PostScript figures (5 files) embedded with eps

Microscopic elasticity of complex systems

Lecture Notes for the Erice Summer School 2005 Computer Simulations in Condensed Matter: from Materials to Chemical Biology. Perspectives in celebration of the 65th Birthday of Mike Klein organized by Kurt Binder, Giovanni Ciccotti and Mauro Ferrari

Continuum theory of partially fluidized granular flows

A continuum theory of partially fluidized granular flows is developed. The theory is based on a combination of the equations for the flow velocity and shear stresses coupled with the order parameter equation which describes the transition between flowing and static components of the granular system. We apply this theory to several important granular problems: avalanche flow in deep and shallow inclined layers, rotating drums and shear granular flows between two plates. We carry out quantitative comparisons between the theory and experiment.Comment: 28 pages, 23 figures, submitted to Phys. Rev.

Constitutive relations for compressible granular flow in the inertial regime

Granular flows occur in a wide range of situations of practical interest to industry, in our natural environment and in our everyday lives. This paper focuses on granular flow in the so-called inertial regime, when the rheology is independent of the very large particle stiffness. Such flows have been modelled with the μ(I),Φ(I)-rheology, which postulates that the bulk friction coefficient μ (i.e. the ratio of the shear stress to the pressure) and the solids volume fraction ϕ are functions of the inertial number I only. Although the μ(I),Φ(I)-rheology has been validated in steady state against both experiments and discrete particle simulations in several different geometries, it has recently been shown that this theory is mathematically ill-posed in time-dependent problems. As a direct result, computations using this rheology may blow up exponentially, with a growth rate that tends to infinity as the discretization length tends to zero, as explicitly demonstrated in this paper for the first time. Such catastrophic instability due to ill-posedness is a common issue when developing new mathematical models and implies that either some important physics is missing or the model has not been properly formulated. In this paper an alternative to the μ(I),Φ(I)-rheology that does not suffer from such defects is proposed. In the framework of compressible I-dependent rheology (CIDR), new constitutive laws for the inertial regime are introduced; these match the well-established μ(I) and Φ(I) relations in the steady-state limit and at the same time are well-posed for all deformations and all packing densities. Time-dependent numerical solutions of the resultant equations are performed to demonstrate that the new inertial CIDR model leads to numerical convergence towards physically realistic solutions that are supported by discrete element method simulations

Patterns and Collective Behavior in Granular Media: Theoretical Concepts

Granular materials are ubiquitous in our daily lives. While they have been a subject of intensive engineering research for centuries, in the last decade granular matter attracted significant attention of physicists. Yet despite a major efforts by many groups, the theoretical description of granular systems remains largely a plethora of different, often contradicting concepts and approaches. Authors give an overview of various theoretical models emerged in the physics of granular matter, with the focus on the onset of collective behavior and pattern formation. Their aim is two-fold: to identify general principles common for granular systems and other complex non-equilibrium systems, and to elucidate important distinctions between collective behavior in granular and continuum pattern-forming systems.Comment: Submitted to Reviews of Modern Physics. Full text with figures (2Mb pdf) avaliable at http://mti.msd.anl.gov/AransonTsimringReview/aranson_tsimring.pdf Community responce is appreciated. Comments/suggestions send to [email protected]

Metastable liquid-liquid phase transition in a single-component system with only one crystal phase and no density anomaly

We investigate the phase behavior of a single-component system in 3 dimensions with spherically-symmetric, pairwise-additive, soft-core interactions with an attractive well at a long distance, a repulsive soft-core shoulder at an intermediate distance, and a hard-core repulsion at a short distance, similar to potentials used to describe liquid systems such as colloids, protein solutions, or liquid metals. We showed [Nature {\bf 409}, 692 (2001)] that, even with no evidences of the density anomaly, the phase diagram has two first-order fluid-fluid phase transitions, one ending in a gas--low-density liquid (LDL) critical point, and the other in a gas--high-density liquid (HDL) critical point, with a LDL-HDL phase transition at low temperatures. Here we use integral equation calculations to explore the 3-parameter space of the soft-core potential and we perform molecular dynamics simulations in the interesting region of parameters. For the equilibrium phase diagram we analyze the structure of the crystal phase and find that, within the considered range of densities, the structure is independent of the density. Then, we analyze in detail the fluid metastable phases and, by explicit thermodynamic calculation in the supercooled phase, we show the absence of the density anomaly. We suggest that this absence is related to the presence of only one stable crystal structure.Comment: 15 pages, 21 figure

The Influence of the Degree of Heterogeneity on the Elastic Properties of Random Sphere Packings

The macroscopic mechanical properties of colloidal particle gels strongly depend on the local arrangement of the powder particles. Experiments have shown that more heterogeneous microstructures exhibit up to one order of magnitude higher elastic properties than their more homogeneous counterparts at equal volume fraction. In this paper, packings of spherical particles are used as model structures to computationally investigate the elastic properties of coagulated particle gels as a function of their degree of heterogeneity. The discrete element model comprises a linear elastic contact law, particle bonding and damping. The simulation parameters were calibrated using a homogeneous and a heterogeneous microstructure originating from earlier Brownian dynamics simulations. A systematic study of the elastic properties as a function of the degree of heterogeneity was performed using two sets of microstructures obtained from Brownian dynamics simulation and from the void expansion method. Both sets cover a broad and to a large extent overlapping range of degrees of heterogeneity. The simulations have shown that the elastic properties as a function of the degree of heterogeneity are independent of the structure generation algorithm and that the relation between the shear modulus and the degree of heterogeneity can be well described by a power law. This suggests the presence of a critical degree of heterogeneity and, therefore, a phase transition between a phase with finite and one with zero elastic properties.Comment: 8 pages, 6 figures; Granular Matter (published online: 11. February 2012