A granular material is a conglomeration of discrete solid, macroscopic particles characterized by a loss of energy whenever the particles interact (the most common example would be friction when grains collide). The constituents that compose granular material are large enough such that they are not subject to thermal motion fluctuations. Thus, the lower size limit for grains in granular material is about 1 μm. On the upper size limit, the physics of granular materials may be applied to ice floes where the individual grains are icebergs and to asteroid belts of the Solar System with individual grains being asteroids. Some examples of granular materials are snow, nuts, coal, sand, rice, coffee, corn flakes, fertilizer, and bearing balls. Research into granular materials is thus directly applicable and goes back at least to Charles-Augustin de Coulomb, whose law of friction was originally stated for granular materials. Granular materials are commercially important in applications as diverse as pharmaceutical industry, agriculture, and energy production. Powders are a special class of granular material due to their small particle size, which makes them more cohesive and more easily suspended in a gas. The soldier/physicist Brigadier Ralph Alger Bagnold was an early pioneer of the physics of granular matter and whose book ''The Physics of Blown Sand and Desert Dunes'' remains an important reference to this day. According to material scientist Patrick Richard, "Granular materials are ubiquitous in nature and are the second-most manipulated material in industry (the first one is water)". In some sense, granular materials do not constitute a single phase of matter but have characteristics reminiscent of solids, liquids, or gases depending on the average energy per grain. However, in each of these states, granular materials also exhibit properties that are unique. Granular materials also exhibit a wide range of pattern forming behaviors when excited (e.g. vibrated or allowed to flow). As such granular materials under excitation can be thought of as an example of a complex system.


Granular Matter is a system composed of many macroscopic particles. Microscopic particles (atoms\molecules) are described (in classical mechanics) by all DOF of the system. Macroscopic particles are described only by DOF of the motion of each particle as a rigid body. In each particle are a lot of internal DOF. Consider inelastic collision between two particles - the energy from velocity as rigid body is transferred to microscopic internal DOF. We get “Dissipation” - irreversible heat generation. The result is that without external driving, eventually all particles will stop moving. In macroscopic particles thermals fluctuations are irrelevant. When a matter is dilute and dynamic (driven) then it is called granular gas and dissipation phenomenon dominates. When a matter is dense and static, then it is called granular solid and jamming phenomenon dominates. When the density is intermediate, then it is called granular liquid.

Static behaviors

Coulomb friction Law

Coulomb regarded internal forces between granular particles as a friction process, and proposed the friction law, that the force of friction of solid particles is proportional to the normal pressure between them and the static friction coefficient is greater than the kinetic friction coefficient. He studied the collapse of piles of sand and found empirically two critical angles: the maximal stable angle \theta_m and the minimum angle of repose \theta_r. When the sandpile slope reaches the maximum stable angle, the sand particles on the surface of the pile begin to fall. The process stops when the surface inclination angle is equal to the angle of repose. The difference between these two angles, \Delta \theta=\theta_m - \theta_r, is the Bagnold angle, which is a measure of the hysteresis of granular materials. This phenomenon is due to the force chains: stress in a granular solid is not distributed uniformly but is conducted away along so-called force chains which are networks of grains resting on one another. Between these chains are regions of low stress whose grains are shielded for the effects of the grains above by vaulting and arching. When the shear stress reaches a certain value, the force chains can break and the particles at the end of the chains on the surface begin to slide. Then, new force chains form until the shear stress is less than the critical value, and so the sandpile maintains a constant angle of repose.

Janssen Effect

In 1895, H. A. Janssen discovered that in a vertical cylinder filled with particles, the pressure measured at the base of the cylinder does not depend on the height of the filling, unlike Newtonian fluids at rest which follow Stevin's law. Janssen suggested a simplified model with the following assumptions: 1) The vertical pressure, \sigma_, is constant in the horizontal plane; 2) The horizontal pressure, \sigma_, is proportional to the vertical pressure \sigma_, where K=\frac is constant in space; 3) The wall friction static coefficient \mu = \frac sustains the vertical load at the contact with the wall; 4) The density of the material is constant over all depths. The pressure in the granular material is then described in a different law, which accounts for saturation:p(z)=p_\infin -\exp(-z/\lambda)/math>where \lambda = \frac and R is the radius of the cylinder, and at the top of the silo z=0. The given pressure equation does not account for boundary conditions, such as the ratio between the particle size to the radius of the silo. Since the internal stress of the material cannot be measured, Janssen's speculations have not been verified by any direct experiment.

Rowe Stress - Dilatancy Relation

In the early 1960s, Rowe studied dilatancy effect on shear strength in shear tests and proposed a relation between them. The mechanical properties of assembly of mono-dispersed particles in 2D can by analyzed based on the representative elementary volume, with typical lengths, \ell_1, \ell_2, in vertical and horizontal directions respectively. The geometric characteristics of the system is described by \alpha=\arctan(\frac) and the variable \beta, which describes the angle when the contact points begin the process of sliding. Denote by \sigma_ the vertical direction, which is the direction of the major principal stress, and by \sigma_ the horizontal direction, which is the direction of the minor principal stress. Then stress on the boundary can be expressed as the concentrated force borne by individual particles. Under biaxial loading with uniform stress \sigma_=\sigma_=0 and therefore F_=F_=0. At equilibrium state: \frac=\frac=\tan(\theta +\beta) where \theta , the friction angle, is the angle between the contact force and the contact normal direction. \theta_, which describes the angle that if the tangential force falls within the friction cone the particles would still remain steady. It is determined by the coefficient of friction \mu=tg\phi_u, so \theta \leq \theta_\mu. Once stress is applied to the system then \theta gradually increases while \alpha,\beta remains unchanged. When \theta \geq \theta_ then the particles will begin sliding, resulting in changing the structure of the system and creating new force chains. \Delta_1,\Delta_2,the horizontal and vertical displacements respectively satisfies: \frac=\frac=-\tan\beta

Granular gases

If the granular material is driven harder such that contacts between the grains become highly infrequent, the material enters a gaseous state. Correspondingly, one can define a granular temperature equal to the root mean square of grain velocity fluctuations that is analogous to thermodynamic temperature. Unlike conventional gases, granular materials will tend to cluster and clump due to the dissipative nature of the collisions between grains. This clustering has some interesting consequences. For example, if a partially partitioned box of granular materials is vigorously shaken then grains will over time tend to collect in one of the partitions rather than spread evenly into both partitions as would happen in a conventional gas. This effect, known as the granular Maxwell's demon, does not violate any thermodynamics principles since energy is constantly being lost from the system in the process.

Ulam Model

Consider N particles each has energy.\varepsilon_. at some constant rate per unit time, randomly choose two particles with energies \varepsilon_,\varepsilon_ and compute the sum \varepsilon_+\varepsilon_. Now, randomly distribute the total energy between the two particles: choose randomly z\in\left,1\right/math> so that the first particle, after the collision, has energy z\left(\varepsilon_+\varepsilon_\right), and the second \left(1-z\right)\left(\varepsilon_+\varepsilon_\right). the stochastic evolution equation:\varepsilon_(t+dt)=\begin \varepsilon_(t) & probability:\,1-\Gamma dt\\ z\left(\varepsilon_(t)+\varepsilon_(t)\right) & probability:\,\Gamma dt \endwhere \Gamma is the collision rate, z is randomly picked from \left,1\right/math> (uniform distribution) and j is an index also randomly chosen from a uniform distribution. The average energy per particle: \begin \left\langle \varepsilon(t+dt)\right\rangle & =\left(1-\Gamma dt\right)\left\langle \varepsilon(t)\right\rangle +\Gamma dt\cdot\left\langle z\right\rangle \left(\left\langle \varepsilon_\right\rangle +\left\langle \varepsilon_\right\rangle \right)\\ & =\left(1-\Gamma dt\right)\left\langle \varepsilon(t)\right\rangle +\Gamma dt\cdot\dfrac\left(\left\langle \varepsilon(t)\right\rangle +\left\langle \varepsilon(t)\right\rangle \right)\\ & =\left\langle \varepsilon(t)\right\rangle \end The second moment: \begin \left\langle \varepsilon^(t+dt)\right\rangle & =\left(1-\Gamma dt\right)\left\langle \varepsilon^(t)\right\rangle +\Gamma dt\cdot\left\langle z^\right\rangle \left\langle \varepsilon_^+2\varepsilon_\varepsilon_+\varepsilon_^\right\rangle \\ & =\left(1-\Gamma dt\right)\left\langle \varepsilon^(t)\right\rangle +\Gamma dt\cdot\dfrac\left(2\left\langle \varepsilon^(t)\right\rangle +2\left\langle \varepsilon(t)\right\rangle ^\right) \end Now the time derivative of the second moment: \dfrac=lim_\dfrac=-\dfrac\left\langle \varepsilon^\right\rangle +\dfrac\left\langle \varepsilon\right\rangle ^ In steady state: \dfrac=0\Rightarrow\left\langle \varepsilon^\right\rangle =2\left\langle \varepsilon\right\rangle ^ Solving the differential equation for the second moment: \left\langle \varepsilon^\right\rangle -2\left\langle \varepsilon\right\rangle ^=\left(\left\langle \varepsilon^(0)\right\rangle -2\left\langle \varepsilon(0)\right\rangle ^\right)e^ However, instead of characterizing the moments, we can analytically solve the energy distribution, from the moment generating function. Consider the Laplace transform: g(\lambda)=\left\langle e^\right\rangle =\int_^e^\rho(\varepsilon)d\varepsilon . Where g(0)=1 , and \dfrac=-\int_^\varepsilon e^\rho(\varepsilon)d\varepsilon=-\left\langle \varepsilon\right\rangle the n derivative: \dfrac=\left(-1\right)^\int_^\varepsilon^e^\rho(\varepsilon)d\varepsilon=\left\langle \varepsilon^\right\rangle now: e^=\begin e^ & 1-\Gamma t\\ e^ & \Gamma t \end \left\langle e^\right\rangle =\left(1-\Gamma dt\right)\left\langle e^\right\rangle +\Gamma dt\left\langle e^\right\rangle g\left(\lambda,t+dt\right)=\left(1-\Gamma dt\right)g\left(\lambda,t\right)+\Gamma dt\int_^\undersetdz Solving for g(\lambda) with change of variables \delta=\lambda z : \lambda g(\lambda)=\int_^g^(\delta)d\delta\Rightarrow\lambda g'(\lambda)+g(\lambda)=g^(\lambda)\Rightarrow g(\lambda)=\dfrac We will show that \rho(\varepsilon)=\dfrace^ (Boltzmann Distribution) by taking its Laplace transform and calculate the generating function: \int_^\dfrace^\cdot e^d\varepsilon=\dfrac\int_^e^d\varepsilon=-\dfrace^|_^=\dfrac=g(\lambda)

Jamming transition

Granular systems are known to exhibit jamming and undergo a jamming transition which is thought of as a thermodynamic phase transition to a jammed state. The transition is from fluid-like phase to a solid-like phase and it is controlled by temperature, T, volume fraction, \phi, and shear stress, \Sigma. The normal phase diagram of glass transition is in the \phi ^-T plane and it is divided into a jammed state region and unjammed liquid state by a transition line. The phase diagram for granular matter lies in the \phi^-\Sigma plane, and the critical stress curve \Sigma(\phi) divides the state phase to jammed\unjammed region, which corresponds to granular solids\liquids respectively. For isotropically jammed granular system, when \phi is reduced around a certain point, J,the bulk and shear moduli approach 0. The J point corresponds to the critical volume fraction \phi_c. Define the distance to point J, the critical volume fraction, \Delta\phi\equiv\phi-\phi_c. The behavior of granular systems near the J point was empirically found to resemble second-order transition: the bulk modulus shows a power law scaling with \Delta\phi and there are some divergent characteristics lengths when \Delta\phi approaches zero. While \phi_c is constant for an infinite system, for a finite system boundary effects result in a distribution of \phi_c over some range. The Lubachevsky-Stillinger algorithm of jamming allows one to produce simulated jammed granular configurations.

Pattern formation

Excited granular matter is a rich pattern-forming system. Some of the pattern-forming behaviours seen in granular materials are: * The un-mixing or segregation of unlike grains under vibration and flow. An example of this is the so-called Brazil nut effect where Brazil nuts rise to the top of a packet of mixed nuts when shaken. The cause of this effect is that when shaken, granular (and some other) materials move in a circular pattern. some larger materials (Brazil nuts) get stuck while going down the circle and therefore stay on the top. * The formation of structured surface or bulk patterns in vibrated granular layers. These patterns include but are not limited to stripes, squares and hexagons. These patterns are thought to be formed by fundamental excitations of the surface known as oscillons. The formation of ordered volumetric structures in granular materials is known as Granular Crystallisation, and involves a transition from a random packing of particles to an ordered packing such as hexagonal close-packed or body-centred cubic. This is most commonly observed in granular materials with narrow size distributions and uniform grain morphology. * The formation of sand ripples, dunes, and sandsheets Some of the pattern-forming behaviours have been possible to reproduce in computer simulations. There are two main computational approaches to such simulations, time-stepped and event-driven, the former being the most efficient for a higher density of the material and the motions of a lower intensity, and the latter for a lower density of the material and the motions of a higher intensity.

Acoustic effects

Some beach sands, such as those of the aptly named Squeaky Beach, exhibit squeaking when walked upon. Some desert dunes are known to exhibit booming during avalanching or when their surface is otherwise disturbed. Granular materials discharged from silos produce loud acoustic emissions in a process known as silo honking.


Granulation is the act or process in which primary powder particles are made to adhere to form larger, multiparticle entities called ''granules.''


When water is cooled sufficiently slowly, the randomly positioned molecules rearrange and ice crystallites emerge and grow. Unlike removing energy by cooling, crystallization in granular material can be achieved by external driving. Crystallites have been observed to form and grow in periodically sheared granular matter, where, in contrast to molecular systems, the positions of the individual particles can be tracked in the experiment. Computer simulations for a system of spherical grains reveal that homogeneous crystallization emerges at a volume fraction \phi = 0.646 \pm 0.001. The computer simulations identify the minimal ingredients necessary for granular crystallization. In particular, gravity and friction are not necessary.

Computational modeling of granular materials

Several methods are available for modeling of granular materials. Most of these methods consist of the statistical methods by which various statistical properties, derived from either point data or an image, are extracted and used to generate stochastic models of the granular medium. A recent and comprehensive review of such methods is available i
Tahmasebi and other (2017)
Another alternative for building a pack of granular particles that recently has bee
is based on the level-set algorithm by which the real shape of the particle can be captured and reproduced through the extracted statistics for particles' morphology.

See also

* Aggregate (composite) * Fragile matter * Random close pack * Soil liquefaction * Metal powder * Particulates * Paste (rheology)


External links

Fundamentals of Particle Technology – free book
* * Mester, L.
The new physical-mechanical theory of granular materials
2009, Homonnai, * Pareschi, L., Russo, G., Toscani, G.
Modelling and Numerics of Kinetic Dissipative Systems
Nova Science Publishers, New York, 2006. {{Authority control Category:Discrete-phase flow