A model for the packing of irregularly shaped grains

In this paper we will attempt to address the problem of the packing properties of granular materials composed of irregularly shaped grains (using configurational statistical mechanics). In particular, we will develop a model for a system of irregular grains based upon perturbing a packing of mono- or poly-disperse spheres. In the mono-disperse case we will show that the system packs less densely than a packing of perfect spheres, except when local correlations between configurations of grains are taken into account. The opposite is found to be true for a perturbation expansion based upon poly-disperse spheres. Finally we will show that for a bi-disperse packing of spheres phase segregation occurs for any size ratio and discuss whether this is to be expected.     

Coarsening in granular systems

We review a few representative examples of granular experiments or models where phase separation, accompanied by domain coarsening, is a relevant phenomenon. We first elucidate the intrinsic non-equilibrium, or athermal, nature of granular media. Thereafter, dilute systems, the so-called “granular gases”, are discussed: idealized kinetic models, such as the gas of inelastic hard spheres in the cooling regime, are the optimal playground to study the slow growth of correlated structures, e.g., shear patterns, vortices, and clusters. In fluidized experiments, liquid–gas or solid–gas separations have been observed. In the case of monolayers of particles, phase coexistence and coarsening appear in several different setups, with mechanical or electrostatic energy input. Phenomenological models describe, even quantitatively, several experimental measures, both for the coarsening dynamics and for the dynamic transition between different granular phases. The origin of the underlying bistability is in general related to negative compressibility from granular hydrodynamics computations, even if the understanding of the mechanism is far from complete. A relevant problem, with important industrial applications, is related to the demixing or segregation of mixtures, for instance in rotating tumblers or on horizontally vibrated plates. Finally, the problem of compaction of highly dense granular materials, which is relevant in many practical situations, is usually described in terms of coarsening dynamics: there, bubbles of misaligned grains evaporate, allowing the coalescence of optimally arranged islands and a progressive reduction of the total occupied volume.     

Effective elasticity coefficients of native rocks and consolidated granular matter

The elastic coefficients of binary heterogeneous materials, such as several native rock materials or consolidated granular matter will be determined in terms of a perturbation expansion. Furthermore, in order to check the validity of the obtained results, these are compared with numerical investigations using Boole's model of randomly distributed spheres. Finally, we apply the results on several classes of native rocks and consolidated granular materials.     

Collisional source of the second moment and pressure tensor for inhomogeneous rapid granular flows of highly inelastic spheres

By extending the constitutive theories for homogeneous granular flows of highly inelastic spheres by Richman (J. Rheol 33 (1989) 1293), Chou (J. CSME 16-6 (1995) 577), and Chou and Richman (Physica A 259 (1998) 430), the collisional source of the second moment of fluctuation velocity and pressure tensor were developed in this study for inhomogeneous rapid granular flows of identical smooth highly inelastic spheres. The important mean fields in this flow are the solid fraction, mean velocity, and full second moment of fluctuation velocity. The collisional source of second moment and the collisional flux of momentum are based upon an anisotropic Maxwellian velocity distribution function. The constitutive theory was combined with the experimental results measured by Hsiau and Jang (Exp. Thermal Fluid Sci. 17 (1998) 202) so as to determine the profiles of pressure tensor and collision source of second moment in the inhomogeneous rapid granular shear flows of inelastic spheres. The normal pressure discrepancies were also observed.     

Three-dimensional clogging structures of granular spheres near hopper orifice

The characteristic clogging structures of granular spheres blocking three-dimensional granular flow through hopper outlet are analyzed based on packing structures reconstructed using magnetic resonance imaging techniques.Spheres in clogging structures are arranged in a way with typical features of load-bearing,such as more contacting bonds close to the horizontal plane and more mutually-stabilized contact configurations than packing structures away from the orifice.The requirement of load-bearing inevitably leads to the cooperativity of clogging structures with a correlation length of several particle diameters.This correlation length being comparable with the orifice diameter suggests that a clogging structure is composed of several mutually-stabilized structural motifs to span the orifice perimeter,instead of a collection of independent individual spheres to cover the whole orifice area.Accordingly,we propose a simple geometric model to explain the unexpected linear dependence of the average size of three-dimensional clogging structures on orifice diameter.     

The statistical mechanics of granular systems composed of spheres and elongated grains

In this paper we will apply a statistical mechanical theory of granular materials to a system composed of a mixture of elongated rods and spheres (both mono- and polydisperse). We will show for a system of rods and spheres that the system's constituents do not completely phase segregate even at minimum close packing. The generalisation of this to packings of rods and bi-disperse spheres produces a similar result except that the two species of sphere do phase segregate at some particular packing.     

Dynamically tessellating algorithm for analysis of pore size distribution in particle agglomerates

We describe a novel physical application of the OctTree data structure [P. Meagher, Comput. Graphics Image Process 19(2) (1982) 129–147] in a dynamically tessellating algorithm, in conjunction with an object-oriented, constructive solid geometry library (DOC), to efficiently determine pore size distributions in large multi-particle systems. We apply the DOC library to investigate the evolving dynamics of pore formation in multi-particle systems, such as a mixture of smooth hard cubes and spheres and a collection of frictional soft spheres. We demonstrate that the algorithm is able to provide insight into the effect of structural changes on the porosity network; for example, during the uniaxial compaction of soft spheres, we find the number density of pores increases while the mean volume of the pores decreases. This trend is responsible for a shift in the distribution of the pore volumes to favour smaller volumes. We anticipate that the DOC method will have wider applications in the area of granular materials for studying the changes in pore structure in both experimental and numerical systems as a complement to the analysis of particle packing.     

Vibration-induced granular segregation: a phenomenon driven by three mechanisms

The segregation of large spheres in a granular bed under vertical vibrations is studied. In our experiments, we systematically measure rise times as a function of density, diameter, and depth, for two different sinusoidal excitations. The measurements reveal that, at low frequencies, inertia and convection are the only mechanisms behind segregation. Inertia (convection) dominates when the relative density is greater (less) than one. At high frequencies, where convection is suppressed, fluidization of the granular bed causes either buoyancy or sinkage and segregation occurs.     

Microstructure characterization of granular materials

This paper presents techniques and algorithms to compute microstructure properties of irregular-shaped granulate assemblies utilizing 3D images. The techniques are capable of extracting microstructure properties of particles such as centeroid, particle size distribution, shape indices (i.e., sphericiy and angularity), number of contacts (i.e., distribution of coordination numbers), contact network, packing efficiency, distribution of local void ratio and radial distribution function. Such properties are critical parameters for micromechanical-based numerical models to capture micro- and macromechanical behavior of geomaterials. X-ray microtomography was used to reconstruct high-resolution 3D image of a natural sand system to represent granular materials. Microstructure properties of the sand system were computed and compared with properties of a computer-simulated image of periodic random spheres. Findings indicate that the use of simplified systems of idealized spheres to model micro- and macromechanical behavior of granular systems can lead to inaccurate results due to the differences in the microstructure between both systems. Methods presented in this paper enabled capturing a more realistic microstructure that can be incorporated in micromechanical models to better simulate, understand, or explain macroscale behavior of granular materials based on their actual microstructure.     

Numerical simulation of the sedimentation of a sphere in a sheared granular fluid: a granular Stokes experiment

We study, computationally, the sedimentation of a sphere of higher mass in a steady, gravity-driven granular flow of otherwise identical spheres, on a rough inclined plane. Taking a hydrodynamic approach at the scale of the particle, we find the drag force to be given by a modified Stokes law and the buoyancy force by the Archimedes principle, with excluded volume effects taken into account. We also find significant differences between the hydrodynamic case and the granular case, which are highlighted.     

Probability distribution of inherent states in models of granular media and glasses

The present paper develops a Statistical Mechanics approach to the inherent states of glassy systems and granular materials by following the original ideas proposed by Edwards for granular media. We consider three lattice models (a diluted spin glass, a system of hard spheres under gravity and a hard-spheres binary mixture under gravity) introduced to describe glassy and granular systems. They are evolved using a “tap dynamics” analogous to that of experiments on granular media. We show that the asymptotic states reached in such a dynamics are not dependent on the particular sample history and are characterized by a few thermodynamical parameters. We assume that under stationarity these systems are distributed in their inherent states satisfying the principle of maximum entropy. This leads to a generalized Gibbs distribution characterized by new “thermodynamical” parameters, called “configurational temperatures” (related to Edwards compactivity for granular materials). Finally, we show by Monte Carlo calculations that the average of macroscopic quantities over the tap dynamics and over such distribution indeed coincide. In particular, in the diluted spin glass and in the system of hard spheres under gravity, the asymptotic states reached by the system are found to be described by a single “configurational temperature”. Whereas in the hard-spheres binary mixture under gravity the asymptotic states reached by the system are found to be described by two thermodynamic parameters, coinciding with the two configurational temperatures which characterize the distribution among the inherent states when the principle of maximum entropy is satisfied under the constraint that the energies of the two species are independently fixed. Received 19 March 2002 and Received in final form 14 June 2002     

Subdiffusion and cage effect in a sheared granular material

We investigate experimentally the diffusion properties of a bidimensional bidisperse dry granular material under quasistatic cyclic shear. The comparison of these properties with results obtained both in computer simulations of hard spheres systems and Lennard-Jones liquids and experiments on colloidal systems near the glass transition demonstrates a strong analogy between the statistical behavior of granular matter and these systems, despite their intrinsic microscopic differences (thermal vs athermal). More specifically, we study in detail the cage dynamics responsible for the subdiffusion in the slow relaxation regime, and obtain the values of relevant time and length scales.     

Scaling laws for melting ice avalanches

This Letter describes an investigation of interfacial melting in ice-bearing granular flows. It is proposed that energy associated with granular collisions causes melting at an ice particle's surface, which can thus occur at temperatures well below freezing. A laboratory experiment has been designed that allows quantification of this process and its effect on the dynamics of a granular shear flow of ice spheres. This experiment employs a rotating drum, half filled with ice particles, situated in a temperature controlled laboratory. Capillary forces between the wetted melted particle surfaces lead to the clumping of particles and enhanced flow speeds, in turn leading to further melting. Dimensional analysis defines a parameter space for further experimentation.     

Velocity correlations and the structure of nonequilibrium hard-core fluids

A model for the pair-distribution function of nonequilibrium hard-core fluids is proposed based on a model for the effect of velocity correlations on the structure. Good agreement is found with molecular dynamics simulations of granular fluids and of sheared elastic hard spheres. It is argued that the incorporation of velocity correlations are crucial to correctly modeling atomic scale structure in nonequilibrium fluids.     

Moment Equations for a Granular Material in a Thermal Bath

We compute the moment equations for a granular material under the simplifying assumption of pseudo-Maxwellian particles approximating dissipative hard spheres. We obtain the general moment equations of second and third order and the isotropic moment equations of any order. Our equations describe, in the space homogeneous case, the granular system described by a Boltzmann-like collision term and subject to a Brownian motion due to the interaction with a bath, described by a Fokker–Planck term. The trend to equilibrium is studied in detail.     

On the properties of “short” granular magnets with disordered chains of grains: A field between grains

The values and functional exponential form of the demagnetizing factor of “short” cylindrical samples of the polyspherical medium and samples (elements) of this medium, i.e., chains of spheres, have been established. The tortuosity factor of chains formed by spheres (as effective elementary conductors of the magnetic flux in a magnetized granular medium) has been analyzed. It has been demonstrated that the tortuosity factor for chains in a filling medium is insignificant. The strength of the field between spheres in the chain, like the magnetic permeability of the effective magnetization channel, has an extremum (bell-shaped) profile. It has been revealed (by using magnetophoresis of finely dispersed ferroparticles) that the gradient of this strength exhibits a two-extremum profile.     

On the force of electrostatic interaction between two conducting spheres

A compact solution is obtained to the problem on the force of interaction between two conducting spheres with preset charges on their surfaces in zero external field. The derivation is based on exact solution of the problem of the potential distribution in the bispherical coordinate system. The expression for the force was derived by differentiating the potential energy of interaction between the spheres with respect to the distance between their centers. It is shown using numerical calculations that with decreasing distance between the spheres, the ratio of their charges for which the forces of interaction between the charges are zero tends to the ratio of the charges of contacting spheres. It follows hence that for any ratio of charges of the same polarity, which differs from the ratio of charges of the contacting spheres, there always exists a small distance between the spheres, at which they attract each other.     

Jumping Platonic Solids on a Vibrating Plate

The mechanical excitation of solid bodies by a vibrating plate is not only an interesting fundamental problem, but also has relevance in many mechanical systems, for example, in the preparation of granular gases in microgravity. Herein, the energy input by an oscillating plate is numerically investigated as a function of the excitation parameters for selected Platonic solids and their behavior compared to that of a jumping sphere. The most important additional features, not relevant for spheres, are the excitation of body rotations and a permanent energy exchange between the rotational and translational degrees of freedom during the collisions with the plate. The distribution of kinetic energies is analyzed by a numerical simulation of the dynamics, using structures which emulate the mechanical behavior of regular polyhedra.     

Jamming II: Edwards’ statistical mechanics of random packings of hard spheres

The problem of finding the most efficient way to pack spheres has an illustrious history, dating back to the crystalline arrays conjectured by Kepler and the random geometries explored by Bernal in the 1960s. This problem finds applications spanning from the mathematician’s pencil, the processing of granular materials, the jamming and glass transitions, all the way to fruit packing in every grocery. There are presently numerous experiments showing that the loosest way to pack spheres gives a density of ∼55% (named random loose packing, RLP) while filling all the loose voids results in a maximum density of ∼63%-64% (named random close packing, RCP). While those values seem robustly true, to this date there is no well-accepted physical explanation or theoretical prediction for them. Here we develop a common framework for understanding the random packings of monodisperse hard spheres whose limits can be interpreted as the experimentally observed RLP and RCP. The reason for these limits arises from a statistical picture of jammed states in which the RCP can be interpreted as the ground state of the ensemble of jammed matter with zero compactivity, while the RLP arises in the infinite compactivity limit. We combine an extended statistical mechanics approach ‘a la Edwards’ (where the role traditionally played by the energy and temperature in thermal systems is substituted by the volume and compactivity) with a constraint on mechanical stability imposed by the isostatic condition. We show how such approaches can bring results that can be compared to experiments and allow for an exploitation of the statistical mechanics framework. The key result is the use of a relation between the local Voronoi volumes of the constituent grains (denoted the volume function) and the number of neighbors in contact that permits us to simply combine the two approaches to develop a theory of volume fluctuations in jammed matter. Ultimately, our results lead to a phase diagram that provides a unifying view of the disordered hard sphere packing problem and further sheds light on a diverse spectrum of data, including the RLP state. Theoretical results are well reproduced by numerical simulations that confirm the essential role played by friction in determining both the RLP and RCP limits. The RLP values depend on friction, explaining why varied experimental results can be obtained.     

Crystallization of confined non-Brownian spheres by vibrational annealing

We introduce an experimental method to crystallize ensembles of non-Brownian spheres confined in narrow containers. The method is based on programmed vibrations and a cooling procedure (annealing). Starting with a granular gas, the system slowly relaxes into a solid ordered structure: Body-centered-tetragonal and face-centered-cubic single crystals are obtained depending on the dimensions of the capillaries. Dry and lubricated beads behave differently, indicating that a sticking coefficient between the particles is important in the dynamics of the crystallization.