Elegance of disordered granular packings: a validation of Edward's hypothesis

We have found a way to analyze Edwards' density of states for static granular packings in the special case of round, rigid, frictionless grains assuming a constant coordination number. It obtains the most entropic density of single grain states, which predicts several observables including the distribution of contact forces. We compare these results against empirical data obtained in dynamic simulations of granular packings. The agreement is quite good, helping validate the use of statistical mechanics methods in granular physics. The differences between theory and empirics are mainly related to the coordination number, and when the empirical data are sorted by that number we obtain several insights that suggest an underlying elegance in the density of states.     

Entropy and temperature of a static granular assembly: an ab initio approach

We construct a statistical framework for static assemblies of deformable grains which parallels that of equilibrium statistical mechanics but with a conservation principle based on the mechanical stress tensor. We define a state function that has all the attributes of entropy. In particular, maximizing this function leads to a well-defined granular temperature and the equivalent of the Boltzman distribution for ensembles of grain packings. Predictions of the ensemble are verified against simulated packings of frictionless, deformable disks.     

Simulation and Statistical Analysis of Random Packings of Ellipsoids

The paper describes an algorithm for the generation of random packings of ellipsoids of revolution, which is a natural generalization of an older algorithm for the case of spherical particles and belongs to the class of collective rearrangement algorithms. It yields packings with a broad spectrum of densities and geometrical properties, aspect ratio which depend on and the various parameters controlling the motions of ellipsoids during the simulation. The ellipsoid systems obtained are characterized by second‐order characteristics, namely pair correlation and orientation correlation functions. Furthermore, the constructed packings are described in a phase diagram taken from statistical physics.     

Voronoï tessellation reveals the condensed matter character of folded proteins

The packing geometry of amino acids in folded proteins is analyzed via a modified Vorono? tessellation method which distinguishes bulk and surface. From a statistical analysis of the Vorono? cells over 40 representative proteins, it appears that the packings are in average similar to random packings of hard spheres encountered in condensed matter physics, with a quite strong fivefold local symmetry. Moreover, the statistics permits one to establish a classification of amino acids in terms of increasing propensity to be buried in agreement with what is known from chemical considerations.     

Force network ensemble: a new approach to static granular matter

An ensemble approach for force distributions in static granular packings is developed. This framework is based on the separation of packing and force scales, together with an a priori flat measure in the force phase space under the constraints that the contact forces are repulsive and balance on every particle. We show how the formalism yields realistic results, both for disordered and regular triangular "snooker ball" configurations, and obtain a shear-induced unjamming transition of the type proposed recently for athermal media.     

Dynamical Janssen effect on granular packing with moving walls

Apparent mass (M(app)) measurements at the bottom of granular packings inside a vertical tube in relative motion are reported. They demonstrate that Janssen's model is valid over a broad range of velocities v. The variability of the measurements is lower than for static packings and the theoretical exponential increase of M(app) with the height of the packing is precisely followed (the corresponding characteristic screening length is of the order of the tube diameter). The limiting apparent mass at large heights is independent of v and significantly lower than the static value.     

Granular statistical mechanics – a personal perspective

The science of granular matter has expanded from an activity for specialised engineering applications to a fundamental field in its own right. This has been accompanied by an explosion of research and literature, which cannot be reviewed in one paper. A key to progress in this field is the formulation of a statistical mechanical formalism that could help develop equations of state and constitutive relations. This paper aims at reviewing some milestones in this direction. An essential basic step toward the development of any static and quasi-static theory of granular matter is a systematic and useful method to quantify the grain-scale structure and we start with a review of such a method. We then review and discuss the ongoing attempt to construct a statistical mechanical theory of granular systems. Along the way, we will clarify a number of misconceptions in the field, as well as highlight several outstanding problems.     

Contact forces in a granular packing

We present the results of a systematic numerical investigation of force distributions in granular packings. We find that all the main features of force transmission previously established for two-dimensional systems of hard particles hold in three-dimensional systems and for soft particles, too. In particular, the probability distribution of normal forces falls off exponentially for forces above the mean force. For forces below the mean, this distribution is either a decreasing power law when the system is far from static equilibrium, or nearly uniform at static equilibrium, in agreement with recent experiments. Moreover, we show that the forces below the mean do not contribute to the shear stress. The subnetwork of the contacts carrying a force below the mean thus plays a role similar to a fluid surrounding the solid backbone composed of the contacts carrying a force above the mean. We address the issue of the computation of contact forces in a packing at static equilibrium. We introduce a model with no local simplifying force rules, that allows for an exact computation of contact forces for given granular texture and boundary conditions. (c) 1999 American Institute of Physics.     

Strain stiffening in random packings of entangled granular chains

Random packings of granular chains are presented as a model system to investigate the contribution of entanglements to strain stiffening. The chain packings are sheared in uniaxial compression experiments. For short chain lengths, these packings yield when the shear stress exceeds the scale of the confining pressure, similar to granular packings of unconnected particles. In contrast, packings of chains which are long enough to form loops exhibit strain stiffening, in which the effective stiffness of the material increases with strain, similar to many polymer materials. The latter packings can sustain stresses orders-of-magnitude greater than the confining pressure, and do not yield until the chain links break. X-ray tomography measurements reveal that the strain-stiffening packings contain system-spanning clusters of entangled chains.     

Stress propagation through frictionless granular material

We examine the network of forces to be expected in a static assembly of hard, frictionless spherical beads of random sizes, such as a colloidal glass. Such an assembly is minimally connected: the ratio of constraint equations to contact forces approaches unity for a large assembly. However, the bead positions in a finite subregion of the assembly are underdetermined. Thus to maintain equilibrium, half of the exterior contact forces are determined by the other half. We argue that the transmission of force may be regarded as unidirectional, in contrast to the transmission of force in an elastic material. Specializing to sequentially deposited beads, we show that forces on a given buried bead can be uniquely specified in terms of forces involving more recently added beads. We derive equations for the transmission of stress averaged over scales much larger than a single bead. This derivation requires the ansatz that statistical fluctuations of the forces are independent of fluctuations of the contact geometry. Under this ansatz, the d(d+1)/2-component stress field can be expressed in terms of a d-component vector field. The procedure may be generalized to nonsequential packings. In two dimensions, the stress propagates according to a wave equation, as postulated in recent work elsewhere. We demonstrate similar wave-like propagation in higher dimensions, assuming that the packing geometry has uniaxial symmetry. In macroscopic granular materials we argue that our approach may be useful even though grains have friction and are not packed sequentially.     

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

The series of equilibrium states reached by disordered packings of rigid, frictionless disks 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 magnitude, the dependence of the resulting strain on stress direction is a Levy flight in the thermodynamic limit.     

Structural transitions in granular packs: statistical mechanics and statistical geometry investigations

We investigate equal spheres packings generated from several experiments and from a large number of different numerical simulations. The structural organization of these disordered packings is studied in terms of the network of common neighbours. This geometrical analysis reveals sharp changes in the network’s clustering occurring at the packing fractions (fraction of volume occupied by the spheres respect to the total volume, ρ) corresponding to the so called Random Loose Packing limit (RLP, ρ ～ 0.555) and Random Close Packing limit (RCP, ρ ～ 0.645). At these packing fractions we also observe abrupt changes in the fluctuations of the portion of free volume around each sphere. We analyze such fluctuations by means of a statistical mechanics approach and we show that these anomalies are associated to sharp variations in a generalized thermodynamical variable which is the analogous for these a-thermal systems to the specific heat in thermal systems.     

Elasticity from the force network ensemble in granular media

Transmission of forces in static granular materials is studied within the framework of the force network ensemble, by numerically evaluating the mechanical response of hexagonal packings of frictionless grains and rectangular packings of frictional grains. In both cases, close to the point of application of the overload, the response is nonlinear and displays two peaks, while at larger length scales, it is linear and elasticlike. The crossover between these two behaviors occurs at a depth that increases with the magnitude of the overload and decreases with increasing friction.     

Discrete Element Method studies of the collision of one rapid sphere on 2D and 3D packings

We performed numerical simulations of one-bead collision on the surface of a static granular medium. The simulations have been done for two- and three-dimensional packings of beads. The effect of the incident bead velocity, the shot angle, the mechanical parameters and the packing structure are analyzed for ordered and disordered 2D packings and only disordered 3D packings. The 2D results are in good agreement with experimental available data. The 3D simulations give good preliminaries results about the shock-wave propagation through the stacking and provides new insights in the ejection process (“splash function”).     

Correlations in the elastic response of dense random packings

Results are presented for the autocorrelation function of the vortexlike nonaffine piece of the linear elastic displacement field in dense random bidisperse packings of harmonically repulsive disks in 2D. The autocorrelation function is shown to scale precisely with the length of the simulation cell in systems ranging from 20 to 100 particles across. It is shown that, to first order, the displacement fields can be thought to arise from the action of uncorrelated local random forcing of a homogeneous elastic sheet, and a theory is presented which gives excellent quantitative agreement with the form of the correlation functions. These results suggest measurements to be made in many types of densely packed, random materials where the elastic displacement fields are accessible experimentally such as granular materials, dense emulsions, colloidal suspensions, etc.     

Physics of particulate flows: From sand avalanche to active suspensions in plants

Flows of granular media in air or in a liquid have been a research field for physicists for several decades. Sometimes solid, sometimes liquid, these particulate materials exhibit peculiar behaviors, which have motivated many studies at the frontiers between nonlinear physics, soft matter physics and fluid mechanics. This paper presents a summary of the recent advances in the field, with a focus on the development of continuous approaches, which make it possible to treat granular media as a complex fluid and to develop a granular hydrodynamics. We also discuss how the better understanding of granular flows we have today may help to address more complex materials, such as colloidal suspensions or some biological systems.     

Jamming III: Characterizing randomness via the entropy of jammed matter

The nature of randomness in disordered packings of frictional and frictionless spheres is investigated using theory and simulations of identical spherical grains. The entropy of the packings is defined through the force and volume ensemble of jammed matter and this is shown to be difficult to calculate analytically. A mesoscopic ensemble of isostatic states is then utilized in an effort to predict the entropy through the definition of a volume function that is dependent on the coordination number. Equations of state are obtained relating entropy, volume fraction and compactivity characterizing the different states of jammed matter, and elucidating the phase diagram for jammed granular matter. Analytical calculations are compared to numerical simulations using volume fluctuation analysis and graph theoretical methods, with reasonable agreement. The entropy of the jammed system reveals that random loose packings are more disordered than random close packings, allowing for an unambiguous interpretation of both limits. Ensemble calculations show that the entropy vanishes at random close packing (RCP), while numerical simulations show that a finite entropy remains in the microscopic states at RCP. The notion of a negative compactivity, which explores states with volume fractions below those achievable by existing simulation protocols, is also explored, expanding the equations of state. The mesoscopic theory reproduces the simulations results in shape well, though a difference in magnitude implies that the entire entropy of the packing may not be captured by the methods presented herein. We discuss possible extensions to the present mesoscopic approach describing packings from random loose packing (RLP) to RCP to the ordered branch of the equation of state in an effort to understand the entropy of jammed matter in the full range of densities from RLP to face-centered cubic (FCC) packing.     

Definition and relevance of nonequilibrium intensive thermodynamic parameters

We show that intensive thermodynamic parameters associated to additive conserved quantities can be naturally defined from a statistical approach in far-from-equilibrium steady-state systems, under few assumptions, and without any detailed balance requirement. It may apply, e.g., to dissipative systems such as granular gases where volume or mass is still conserved or to systems with periodic boundary conditions where fluxes of conserved quantities are present. We emphasize the usefulness of this concept to characterize the coexistence of different nonequilibrium phases and discuss the influence of the contact between two different systems, in relation with measurement issues.     

Stresses in isostatic granular systems and emergence of force chains

Progress is reported on several questions that bedevil understanding of granular systems: (i) Are the stress equations elliptic, parabolic, or hyperbolic? (ii) How can the often-observed force chains be predicted from a first-principles continuous theory? (iii) How do we relate insight from isostatic systems to general packings? Explicit equations are derived for the stress components in two dimensions including the dependence on the local structure. The equations are shown to be hyperbolic and their general solutions, as well as the Green function, are found. It is shown that the solutions give rise to force chains, and the explicit dependence of the force chains trajectories and magnitudes on the local geometry is predicted. Direct experimental tests of the predictions are proposed. Finally, a framework is proposed to relate the analysis to nonisostatic and more realistic granular assemblies.     

Apollonian networks: simultaneously scale-free, small world, euclidean, space filling, and with matching graphs

We introduce a new family of networks, the Apollonian networks, that are simultaneously scale-free, small-world, Euclidean, space filling, and with matching graphs. These networks describe force chains in polydisperse granular packings and could also be applied to the geometry of fully fragmented porous media, hierarchical road systems, and area-covering electrical supply networks. Some of the properties of these networks, namely, the connectivity exponent, the clustering coefficient, and the shortest path are calculated and found to be particularly rich. The percolation, the electrical conduction, and the Ising models on such networks are also studied and found to be quite peculiar. Consequences for applications are also discussed.