Fragility and hysteretic creep in frictional granular jamming

The granular jamming transition is experimentally investigated in a two-dimensional system of frictional, bi-dispersed disks subject to quasi-static, uniaxial compression at zero granular temperature. Currently accepted results show the jamming transition occurs at a critical packing fraction $\phi_c$. In contrast, we observe the first compression cycle exhibits {\it fragility} - metastable configuration with simultaneous jammed and un-jammed clusters - over a small interval in packing fraction ($\phi_1 < \phi < \phi_2$). The fragile state separates the two conditions that define $\phi_c$ with an exponential rise in pressure starting at $\phi_1$ and an exponential fall in disk displacements ending at $\phi_2$. The results are explained through a percolation mechanism of stressed contacts where cluster growth exhibits strong spatial correlation with disk displacements. Measurements with several disk materials of varying elastic moduli $E$ and friction coefficients $\mu$, show friction directly controls the start of the fragile state, but indirectly controls the exponential slope. Additionally, we experimentally confirm recent predictions relating the dependence of $\phi_c$ on $\mu$. Under repetitive loading (compression), the system exhibits hysteresis in pressure, and the onset $\phi_c$ increases slowly with repetition number. This friction induced hysteretic creep is interpreted as the granular pack's evolution from a metastable to an eventual structurally stable configuration. It is shown to depend upon the quasi-static step size $\Delta \phi$ which provides the only perturbative mechanism in the experimental protocol, and the friction coefficient $\mu$ which acts to stabilize the pack.Comment: 12 pages, 10 figure

Partially fluidized shear granular flows: Continuum theory and MD simulations

The continuum theory of partially fluidized shear granular flows is tested and calibrated using two dimensional soft particle molecular dynamics simulations. The theory is based on the relaxational dynamics of the order parameter that describes the transition between static and flowing regimes of granular material. We define the order parameter as a fraction of static contacts among all contacts between particles. We also propose and verify by direct simulations the constitutive relation based on the splitting of the shear stress tensor into a``fluid part'' proportional to the strain rate tensor, and a remaining ``solid part''. The ratio of these two parts is a function of the order parameter. The rheology of the fluid component agrees well with the kinetic theory of granular fluids even in the dense regime. Based on the hysteretic bifurcation diagram for a thin shear granular layer obtained in simulations, we construct the ``free energy'' for the order parameter. The theory calibrated using numerical experiments with the thin granular layer is applied to the surface-driven stationary two dimensional granular flows in a thick granular layer under gravity.Comment: 20 pages, 19 figures, submitted to Phys. Rev.

Stress response inside perturbed particle assemblies

The effect of structural disorder on the stress response inside three dimensional particle assemblies is studied using computer simulations of frictionless sphere packings. Upon applying a localised, perturbative force within the packings, the resulting {\it Green's} function response is mapped inside the different assemblies, thus providing an explicit view as to how the imposed perturbation is transmitted through the packing. In weakly disordered arrays, the resulting transmission of forces is of the double-peak variety, but with peak widths scaling linearly with distance from the source of the perturbation. This behaviour is consistent with an anisotropic elasticity response profile. Increasing the disorder distorts the response function until a single-peak response is obtained for fully disordered packings consistent with an isotropic description.Comment: 8 pages, 7 figure captions To appear in Granular Matte

Effect of boundaries on the force distributions in granular media

The effect of boundaries on the force distributions in granular media is illustrated by simulations of 2D packings of frictionless, Hertzian spheres. To elucidate discrepancies between experimental observations and theoretical predictions, we distinguish between the weight distribution {\cal P} (w) measured in experiments and analyzed in the q-model, and the distribution of interparticle forces P(f). The latter one is robust, while {\cal P}(w) can be obtained once the local packing geometry and P(f) are known. By manipulating the (boundary) geometry, we show that {\cal P}(w) can be varied drastically.Comment: 4 pages, 4 figure

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

A constitutive law for dense granular flows

A continuum description of granular flows would be of considerable help in predicting natural geophysical hazards or in designing industrial processes. However, the constitutive equations for dry granular flows, which govern how the material moves under shear, are still a matter of debate. One difficulty is that grains can behave like a solid (in a sand pile), a liquid (when poured from a silo) or a gas (when strongly agitated). For the two extreme regimes, constitutive equations have been proposed based on kinetic theory for collisional rapid flows, and soil mechanics for slow plastic flows. However, the intermediate dense regime, where the granular material flows like a liquid, still lacks a unified view and has motivated many studies over the past decade. The main characteristics of granular liquids are: a yield criterion (a critical shear stress below which flow is not possible) and a complex dependence on shear rate when flowing. In this sense, granular matter shares similarities with classical visco-plastic fluids such as Bingham fluids. Here we propose a new constitutive relation for dense granular flows, inspired by this analogy and recent numerical and experimental work. We then test our three-dimensional (3D) model through experiments on granular flows on a pile between rough sidewalls, in which a complex 3D flow pattern develops. We show that, without any fitting parameter, the model gives quantitative predictions for the flow shape and velocity profiles. Our results support the idea that a simple visco-plastic approach can quantitatively capture granular flow properties, and could serve as a basic tool for modelling more complex flows in geophysical or industrial applications.Comment: http://www.nature.com/nature/journal/v441/n7094/abs/nature04801.htm

Long-ranged attraction between disordered heterogeneous surfaces

Long-ranged attractions across water between two surfaces that are randomly covered with (mobile) positive and negative charge domains have been attributed to induced correlation of the charges (positive lining up with negative) as the surfaces approach. Here we show, by directly measuring normal forces under a rapid shear field, that these attractions may not in fact be due to such correlations. It is rather the inherent interaction-asymmetry between equally- and between oppositely-charged domains that results in the long-ranged attraction even in the complete absence of any charge correlation

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]

The Role of Friction in Compaction and Segregation of Granular Materials

We investigate the role of friction in compaction and segregation of granular materials by combining Edwards' thermodynamic hypothesis with a simple mechanical model and mean-field based geometrical calculations. Systems of single species with large friction coefficients are found to compact less. Binary mixtures of grains differing in frictional properties are found to segregate at high compactivities, in contrary to granular mixtures differing in size, which segregate at low compactivities. A phase diagram for segregation vs. friction coefficients of the two species is generated. Finally, the characteristics of segregation are related directly to the volume fraction without the explicit use of the yet unclear notion of compactivity.Comment: 9 pages, 6 figures, submitted to Phys. Rev.

Jamming at Zero Temperature and Zero Applied Stress: the Epitome of Disorder

We have studied how 2- and 3- dimensional systems made up of particles interacting with finite range, repulsive potentials jam (i.e., develop a yield stress in a disordered state) at zero temperature and applied stress. For each configuration, there is a unique jamming threshold, $\phi_c$, at which particles can no longer avoid each other and the bulk and shear moduli simultaneously become non-zero. The distribution of $\phi_c$ values becomes narrower as the system size increases, so that essentially all configurations jam at the same $\phi$ in the thermodynamic limit. This packing fraction corresponds to the previously measured value for random close-packing. In fact, our results provide a well-defined meaning for "random close-packing" in terms of the fraction of all phase space with inherent structures that jam. The jamming threshold, Point J, occurring at zero temperature and applied stress and at the random close-packing density, has properties reminiscent of an ordinary critical point. As Point J is approached from higher packing fractions, power-law scaling is found for many quantities. Moreover, near Point J, certain quantities no longer self-average, suggesting the existence of a length scale that diverges at J. However, Point J also differs from an ordinary critical point: the scaling exponents do not depend on dimension but do depend on the interparticle potential. Finally, as Point J is approached from high packing fractions, the density of vibrational states develops a large excess of low-frequency modes. All of these results suggest that Point J may control behavior in its vicinity-perhaps even at the glass transition.Comment: 21 pages, 20 figure