Preavalanche instabilities in a granular pile

Evolution of force networks during stick-slip motion of an intruder in a granular material: Topological measures extracted from experimental data

In quasi-two-dimensional experiments with photoelastic particles confined to an annular region, an intruder constrained to move in a circular path halfway between the annular walls experiences stick-slip dynamics. We discuss the response of the granular medium to the driven intruder, focusing on the evolution of the force network during sticking periods. Because the available experimental data do not include precise information about individual contact forces, we use an approach developed in our previous work [Basak et al., J. Eng. Mech. 147, 04021100 (2021)0733-939910.1061/(ASCE)EM.1943-7889.0002003] based on networks constructed from measurements of the integrated strain magnitude on each particle. These networks are analyzed using topological measures based on persistence diagrams, revealing that force networks evolve smoothly but in a nontrivial manner throughout each sticking period, even though the intruder and granular particles are stationary. Characteristic features of persistence diagrams show identifiable slip precursors. In particular, the number of generators describing the structure and complexity of force networks increases consistently before slips. Key features of the dynamics are similar for granular materials composed of disks or pentagons, but some details are consistently different. In particular, we find significantly larger fluctuations of the measures computed based on persistence diagrams and, therefore, of the underlying networks, for systems of pentagonal particles.

Stochastic methods for slip prediction in a sheared granular system

We consider a sheared granular system experiencing intermittent dynamics of stick-slip type via discrete element simulations. The considered setup consists of a two-dimensional system of soft frictional particles sandwiched between solid walls, one of which is exposed to a shearing force. The slip events are detected using stochastic state space models applied to various measures describing the system. The amplitudes of the events spread over more than four decades and present two distinctive peaks, one for the microslips and the other for the slips. We show that the measures describing the forces between the particles provide earlier detection of an upcoming slip event than the measures based solely on the wall movement. By comparing the detection times obtained from the considered measures, we observe that a typical slip event starts with a local change in the force network. However, some local changes do not spread globally over the force network. For the changes that become global, we find that their size strongly influences the further behavior of the system. If the size of a global change is large enough, then it triggers a slip event; if it is not, then a much weaker microslip follows. Quantification of the changes in the force network is made possible by formulating clear and precise measures describing their static and dynamic properties.

On intermittency in sheared granular systems

We consider a system of granular particles, modeled by two dimensional frictional soft elastic disks, that is exposed to externally applied time-dependent shear stress in a planar Couette geometry. We concentrate on the external forcing that produces intermittent dynamics of stick-slip type. In this regime, the top wall remains almost at rest until the applied stress becomes sufficiently large, and then it slips. We focus on the evolution of the system as it approaches a slip event. Our main finding is that there are two distinct groups of measures describing system behavior before a slip event. The first group consists of global measures defined as system-wide averages at a fixed time. Typical examples of measures in this group are averages of the normal or tangent forces acting between the particles, system size and number of contacts between the particles. These measures do not seem to be sensitive to an approaching slip event. On average, they tend to increase linearly with the force pulling the spring. The second group consists of the time-dependent measures that quantify the evolution of the system on a micro (particle) or mesoscale. Measures in this group first quantify the temporal differences between two states and only then aggregate them to a single number. For example, Wasserstein distance quantitatively measures the changes of the force network as it evolves in time while the number of broken contacts quantifies the evolution of the contact network. The behavior of the measures in the second group changes dramatically before a slip event starts. They increase rapidly as a slip event approaches, indicating a significant increase in fluctuations of the system before a slip event is triggered.

Power-Law Scaling in Granular Rheology across Flow Geometries

Based on discrete element method simulations, we propose a new form of the constitutive equation for granular flows independent of packing fraction. Rescaling the stress ratio μ by a power of dimensionless temperature Θ makes the data from a wide set of flow geometries collapse to a master curve depending only on the inertial number I. The basic power-law structure appears robust to varying particle properties (e.g., surface friction) in both 2D and 3D systems. We show how this rheology fits and extends frameworks such as kinetic theory and the nonlocal granular fluidity model.

Rheology of granular materials composed of crushable particles

We investigate sheared granular materials composed of crushable particles by means of contact dynamics simulations and the bonded-cell model for particle breakage. Each particle is paved by irregular cells interacting via cohesive forces. In each simulation, the ratio of the internal cohesion of particles to the confining pressure, the relative cohesion, is kept constant and the packing is subjected to biaxial shearing. The particles can break into two or more fragments when the internal cohesive forces are overcome by the action of compressive force chains between particles. The particle size distribution evolves during shear as the particles continue to break. We find that the breakage process is highly inhomogeneous both in the fragment sizes and their locations inside the packing. In particular, a number of large particles never break whereas a large number of particles are fully shattered. As a result, the packing keeps the memory of its initial particle size distribution, whereas a power-law distribution is observed for particles of intermediate size due to consecutive fragmentation events whereby the memory of the initial state is lost. Due to growing polydispersity, dense shear bands are formed inside the packings and the usual dilatant behavior is reduced or cancelled. Hence, the stress-strain curve no longer passes through a peak stress, and a progressive monotonic evolution towards a pseudo-steady state is observed instead. We find that the crushing rate is controlled by the confining pressure. We also show that the shear strength of the packing is well expressed in terms of contact anisotropies and force anisotropies. The force anisotropy increases while the contact orientation anisotropy declines for increasing internal cohesion of the particles. These two effects compensate each other so that the shear strength is nearly independent of the internal cohesion of particles.

Avalanche precursors in a frictional model

We present a one-dimensional numerical model based on elastically coupled sliders on a frictional incline of variable tilt. This very simple approach makes it possible to study the precursors to the avalanche and to provide a rationalization of different features that have been observed in experiments. We provide a statistical description of the model leading to master equations describing the state of the system as a function of the angle of inclination. Our central results are the reproduction of large-scale regular events preceding the avalanche, on the one hand, and an analytical approach providing an internal threshold for the outbreak of rearrangements before the avalanche in the system, on the other hand.

Relaxation-type nonlocal inertial-number rheology for dry granular flows

We propose a constitutive model to describe the nonlocality, hysteresis, and several flow features of dry granular materials. Taking the well-known inertial number I as a measure of sheared-induced local fluidization, we derive a relaxation model for I according to the evolution of microstructure during avalanche and dissipation processes. The model yields a nonmonotonic flow law for a homogeneous flow, accounting for hysteretic solid-fluid transition and intermittency in quasistatic flows. For an inhomogeneous flow, the model predicts a generalized Bagnold shear stress revealing the interplay of two microscopic nonlocal mechanisms: collisions among correlated structures and the diffusion of fluidization within the structures. In describing a uniform flow down an incline, the model reproduces the hysteretic starting and stopping heights and the Pouliquen flow rule for mean velocity. Moreover, a dimensionless parameter reflecting the nonlocal effect on the flow is discovered, which controls the transition between Bagnold and creeping flow dynamics.

Stress Transmission and Failure in Disordered Porous Media

By means of extensive lattice-element simulations, we investigate stress transmission and its relation with failure properties in increasingly disordered porous systems. We observe a non-Gaussian broadening of stress probability density functions under tensile loading with increasing porosity and disorder, revealing a gradual transition from a state governed by single-pore stress concentration to a state controlled by multipore interactions and metric disorder. This effect is captured by the excess kurtosis of stress distributions and shown to be nicely correlated with the second moment of local porosity fluctuations, which appears thus as a (dis)order parameter for the system. By generating statistical ensembles of porous textures with varying porosity and disorder, we derive a general expression for the fracture stress as a decreasing function of porosity and disorder. Focusing on critical sites where the local stress is above the global fracture threshold, we also analyze the transition to failure in terms of a coarse-graining length. These findings provide a general framework which can also be more generally applied to multiphase and structural heterogeneous materials.

Bonded-cell model for particle fracture

Particle degradation and fracture play an important role in natural granular flows and in many applications of granular materials. We analyze the fracture properties of two-dimensional disklike particles modeled as aggregates of rigid cells bonded along their sides by a cohesive Mohr-Coulomb law and simulated by the contact dynamics method. We show that the compressive strength scales with tensile strength between cells but depends also on the friction coefficient and a parameter describing cell shape distribution. The statistical scatter of compressive strength is well described by the Weibull distribution function with a shape parameter varying from 6 to 10 depending on cell shape distribution. We show that this distribution may be understood in terms of percolating critical intercellular contacts. We propose a random-walk model of critical contacts that leads to particle size dependence of the compressive strength in good agreement with our simulation data.

Effect of size polydispersity versus particle shape in dense granular media

We present a detailed analysis of the morphology of granular systems composed of frictionless pentagonal particles by varying systematically both the size span and particle shape irregularity, which represent two polydispersity parameters of the system. The microstructure is characterized in terms of various statistical descriptors such as global and local packing fractions, radial distribution functions, coordination number, and fraction of floating particles. We find that the packing fraction increases with the two parameters of polydispersity, but the effect of shape polydispersity for all the investigated structural properties is significant only at low size polydispersity where the positional and/or orientational ordering of the particles prevail. We focus in more detail on the class of side/side contacts, which is the interesting feature of our system as compared to a packing of disks. We show that the proportion of such contacts has weak dependence on the polydispersity parameters. The side- side contacts do not percolate but they define clusters of increasing size as a function of size polydispersity and decreasing size as a function of shape polydispersity. The clusters have anisotropic shapes but with a decreasing aspect ratio as polydispersity increases. This feature is argued to be a consequence of strong force chains (forces above the mean), which are mainly captured by side-side contacts. Finally, the force transmission is intrinsically multiscale, with a mean force increasing linearly with particle size.

Slow relaxation and aging phenomena at the nanoscale in granular materials

Granular matter exhibits a rich variety of dynamic behaviors, for which the role of thermal fluctuations is usually ignored. Here we show that thermal fluctuations can pronouncedly affect contacting nanoscale asperities at grain interfaces and brightly manifest themselves through the influence on nonlinear-acoustic effects. The proposed mechanism based on intrinsic bistability of nanoscale contacts comprises a wealth of slow-dynamics regimes including slow relaxations and aging as universal properties of a wide class of systems with metastable states.

Granular avalanches in a two-dimensional rotating drum with imposed vertical vibration

We present statistics on granular avalanches in a rotating drum with and without imposed vertical vibration. The experiment consists of a quasi-two-dimensional, vertical drum containing pentagonal particles and rotated at a constant angular velocity. The drum rests on an electromagnetic shaker to allow vibration of the assembly as it rotates. We measure time series of the slope of the interface and find that the critical angle for slope failure θ(c) and the resulting angle of repose θ(r) are broadly distributed with an approximate power-law distribution of avalanches θ(c)-θ(r) for large avalanches. The faceted pentagonal grains used lead to significant interlocking with critical and repose angles (θ(c)≈45° and θ(r)≈39°) larger than experiments using spherical grains, even with vibration, and avalanche magnitudes correlated with the prior build-up and anti-correlated with the prior avalanche. We find that the stability of the assembly increases with small vibrations and is destabilized at vibration amplitudes above a dimensionless acceleration (peak acceleration divided by acceleration due to gravity) of Γ=0.2. We also study history dependence of the avalanches by periodically oscillating the drum to compare the initial avalanche upon reversal of shear to steady-state distributions for avalanches during continuous rotation. We observe history dependence as an initial decrease in critical angle upon reversal of the drum rotation direction, indicating that a texture is induced to resist continued shear such that the surface is weaker to reversals in shear direction. Memory of this history is removed by sufficient external vibration (Γ≥0.8), which leads to compaction and relaxation of the surface layer grains responsible for avalanching dynamics, as initial and steady-state avalanche distributions become indistinguishable.

Experimental investigation of plastic deformations before a granular avalanche

We present an experimental study of the deformation inside a granular material that is progressively tilted. We investigate the deformation before the avalanche with a spatially resolved diffusive wave spectroscopy setup. At the beginning of the inclination process, we first observe localized and isolated events in the bulk, with a density which decreases with the depth. As the angle of inclination increases, series of microfailures occur periodically in the bulk, and finally a granular avalanche takes place. The microfailures are observed only when the tilt angles are larger than a threshold angle much smaller than the granular avalanche angle. We have characterized the density of reorganizations and the localization of microfailures. We have also explored the effect of the nature of the grains, the relative humidity conditions, and the packing fraction of the sample. We discuss those observations in the framework of the plasticity of granular matter. Microfailures may then be viewed as the result of the accumulation of numerous plastic events.

Probing the shear-band formation in granular media with sound waves

We investigate the mechanical responses of dense granular materials, using a direct shear box combined with simultaneous acoustic measurements. Measured shear wave speeds evidence the structural change of the material under shear, from the jammed state to the flowing state. There is a clear acoustic signature when the shear band is formed. Subjected to cyclic shear, both shear stress and wave speed show the strong hysteretic dependence on the shear strain, likely associated with the geometry change in the packing structure. Moreover, the correlation function of configuration-specific multiply scattered waves reveals an intermittent behavior before the failure of material.

Simulation of cemented granular materials. II. Micromechanical description and strength mobilization at the onset of macroscopic yielding

This is the second of two papers investigating the mechanical response of cemented granular materials by means of contact dynamics simulations. In this paper, a two-dimensional polydisperse sample with high void ratio is sheared in a load-controlled simple shear numerical device until the stress state of the sample reaches the yield stress. We first study the stress transmission properties of the granular material in terms of the fabric of different subsets of contacts characterized by the magnitude of their normal forces. This analysis highlights the existence of a peculiar force carrying structure in the cemented material, which is reminiscent of the bimodal stress transmission reported for cohesionless granular media. Then, the evolution of contact forces and torques is investigated trying to identify the micromechanical conditions that trigger macroscopic yielding. It is shown that global failure can be associated to the apparition of a group of particles whose contacts fulfill at least one of the local rupture conditions. In particular, these particles form a large region that percolates through the sample at the moment of failure, evidencing the relationship between macroscopic yielding and the emergence of large-scale correlations in the system.

Annular shear of cohesionless granular materials: from the inertial to quasistatic regime

Using discrete simulations, we investigate the behavior of a model granular material within an annular shear cell. Specifically, two-dimensional assemblies of disks are placed between two circular walls, the inner one rotating with prescribed angular velocity, while the outer one may expand or shrink and maintains a constant radial pressure. Focusing on steady state flows, we delineate in parameter space the range of applicability of the recently introduced constitutive laws for sheared granular materials (based on the inertial number). We discuss the two origins of the stronger strain rates observed near the inner boundary, the vicinity of the wall and the heteregeneous stress field in a Couette cell. Above a certain velocity, an inertial region develops near the inner wall, to which the known constitutive laws apply, with suitable corrections due to wall slip, for small enough stress gradients. Away from the inner wall, slow, apparently unbounded creep takes place in the nominally solid material, although its density and shear to normal stress ratio are on the jammed side of the critical values. In addition to rheological characterizations, our simulations provide microscopic information on the contact network and velocity fluctuations that is potentially useful to assess theoretical approaches.

Computer simulation of model cohesive powders: plastic consolidation, structural changes, and elasticity under isotropic loads

The quasistatic behavior of a simple two-dimensional model of a cohesive powder under isotropic loads is investigated by discrete element simulations. We ignore contact plasticity and focus on the effect of geometry and collective rearrangements on the material behavior. The loose packing states, as assembled and characterized in a previous numerical study [Gilabert, Roux, and Castellanos, Phys. Rev. E 75, 011303 (2007)], are observed, under growing confining pressure P , to undergo important structural changes, while solid fraction Phi irreversibly increases (typically, from 0.4-0.5 to 0.75-0.8). The system state goes through three stages, with different forms of the plastic consolidation curve, i.e., Phi as a function of the growing reduced pressure P;{*}=PaF_{0} , defined with adhesion force F0 and grain diameter a . In the low-confinement regime (I), the system undergoes negligible plastic compaction, and its structure is influenced by the assembling process. In regime II the material state is independent of initial conditions, and the void ratio varies linearly with lnP [i.e., Delta(1Phi)=lambdaDelta(lnP;{*}) ], as described in the engineering literature. Plasticity index lambda is reduced in the presence of a small rolling resistance (RR). In the last stage of compaction (III), Phi approaches an asymptotic, maximum solid fraction Phi_{max} , as a power law Phi_{max}-Phi proportional, variant(P;{*});{-alpha} , with alpha approximately 1 , and properties of cohesionless granular packs are gradually retrieved. Under consolidation, while the range xi of fractal density correlations decreases, force patterns reorganize from self-balanced clusters to force chains, with correlative evolutions of force distributions, and elastic moduli increase by a large amount. Plastic deformation events correspond to very small changes in the network topology, while the denser regions tend to move like rigid bodies. Elastic properties are dominated by the bending of thin junctions in loose systems. For growing RR those tend to form particle chains, the folding of which, rather than tensile ruptures, controls plastic compaction.

Evolving loop structure in gradually tilted two-dimensional granular packings

Granular packings, especially near the jamming transition, form fragile networks where small perturbations can lead to destabilization and large scale rearrangements. A key stabilizing element in two dimensions is the contact loop, yet surprisingly little is known about contact loop statistics in realistic granular networks. In this paper, we use particle dynamics to study the evolution of contact loop structure in a gradually tilted two-dimensional granular bed. We find that the resulting contact loop distributions (1) are sensitive to material properties, (2) deviate from the expected structure of a randomly wired lattice, and (3) are uniquely dependent on tilting angle. Also, we introduce a quantitative measure of loop stability xi and show that increased tilting results in a gradual destabilization of individual loops. We briefly discuss the considerations for extending our approach to three dimensions.

Internal states of model isotropic granular packings. III. Elastic properties

In this third and final paper of a series, elastic properties of numerically simulated isotropic packings of spherical beads assembled by different procedures, as described in the first companion paper, and then subjected to a varying confining pressure, as reported in the second companion paper, are investigated. In addition to the pressure, which determines the stiffness of contacts because of Hertz's law, elastic moduli are chiefly sensitive to the coordination number z , which should not be regarded as a function of the packing density. Comparisons of numerical and experimental results for glass beads in the 10 kPa-10 MPa pressure range reveal similar differences between dry samples prepared in a dense state by vibrations and lubricated packings, so that the greater stiffness of the latter, in spite of their lower density, can be attributed to a larger coordination number. Effective medium type approaches, or Voigt and Reuss bounds, provide good estimates of bulk modulus B , which can be accurately bracketed, but badly fail for shear modulus G , especially in low z configurations under low pressure. This is due to the different response of tenuous, fragile networks to changes in load direction, as compared to load intensity. In poorly coordinated packings, the shear modulus, normalized by the average contact stiffness, tends to vary proportionally to the degree of force indeterminacy per unit volume, even though this quantity does not vanish in the rigid limit. The elastic range extends to small strain intervals and compares well with experimental observations on sands. The origins of nonelastic response are discussed. We conclude that elastic moduli provide access to mechanically important information about coordination numbers, which escape direct measurement techniques, and indicate further perspectives.

Stress transmission in wet granular materials

We analyze stress transmission in wet granular media in the pendular state by means of three-dimensional molecular-dynamics simulations. We show that the tensile action of capillary bonds induces a self-stressed particle network organized in two percolating "phases" of positive and negative particle pressures. Various statistical descriptors of the microstructure and bond force network are used to characterize this partition. Two basic properties emerge: 1) the highest particle pressure is located in the bulk of each phase; 2) the lowest pressure level occurs at the interface between the two phases, involving also the largest connectivity of the particles via tensile and compressive bonds. When a confining pressure is applied, the number of tensile bonds falls off and the negative phase breaks into aggregates and isolated sites.

Failure and strengthening of granular slopes under horizontal vibration

We present experimental measurements of a granular slope under horizontal vibration. We use optical particle tracking to measure the motion of surface beads as the slope fails. We find that for all but the largest inclination angles, initial bead motion leads to strengthening rather than an avalanche. The initial motion of the beads is usually intermittent and evolves differently for different preparations, slope angles, and rates of increase in the vibration amplitude. When a specific criterion is chosen to define failure, the Coulomb friction model adequately describes the average acceleration required to produce failure, as long as slope preparation and experimental protocol are constant. However, the observed intermittent motion and rate dependence indicate that strengthening microrearrangements are important features that affect failure of slopes under external perturbations.

Precursors to avalanches in a granular monolayer

We investigate the stability of a granular monolayer composed of spherical grains on an inclined plate. When the tilt angle alpha increases, some reorganizations are observed throughout the pile. The packing fraction rho of the packing evolves by successive jumps. Those discontinuous events precede the collapse of the pile at a critical angle alphac. The occurrence of precursors before avalanches is modeled by stop-and-go motions of blocks due to the competition between sliding friction and the Janssen effect [J. Durand, (Springer-Verlag, New York, 2000)].

Multi-scale analysis of the stress state in a granular slope in transition to failure

By means of contact dynamics simulations, we analyze the stress state in a granular bed slowly tilted toward its angle of repose. An increasingly large number of grains are overloaded in the sense that they are found to carry a stress ratio above the Coulomb yield threshold of the whole packing. Using this property, we introduce a coarse-graining length scale at which all stress ratios are below the packing yield threshold. We show that this length increases with the slope angle and jumps to a length comparable to the depth of the granular bed at an angle below the angle of repose. This transition coincides with the onset of dilation in the packing. We map this transition into a percolation transition of the overloaded grains, and discuss it in terms of long-range correlations and granular slope metastability.

Memory of the unjamming transition during cyclic tiltings of a granular pile

Discrete numerical simulations are performed to study the evolution of the microstructure and the response of a granular packing during successive loading-unloading cycles, consisting of quasistatic rotations in the gravity field between opposite inclination angles. We show that internal variables--e.g., stress and fabric of the pile--exhibit hysteresis during these cycles due to the exploration of different metastable configurations. Interestingly, the hysteretic behavior of the pile strongly depends on the maximal inclination of the cycles, giving evidence of the irreversible modifications of the pile state occurring close to the unjamming transition. More specifically, we show that for cycles with maximal inclination larger than the repose angle, the weak-contact network carries the memory of the unjamming transition. These results demonstrate the relevance of a two-phase description--strong- and weak-contact networks--for a granular system, as soon as it has approached the unjamming transition.

Friction versus texture at the approach of a granular avalanche

We perform an analysis of the granular texture of a granular bed close to stability limit. Our analysis is based on a unique criterion of friction mobilization in a simulated two-dimensional packing. In this way, we recover the bimodal character of granular texture and the coexistence of weak and strong phases in the sense of distinct contacts populations. Moreover, we show the existence of a well-defined subset of contacts within the weak contact network. These contacts are characterized by their important friction and form a highly coherent population in terms of fabric. They play an antagonistic role with respect to force chains. Thus, we are able to discriminate between incoherent contacts and coherent contacts in the weak phase and to specify the role that the latter plays in the destabilization process.

Rheophysics of dense granular materials: discrete simulation of plane shear flows

We study the plane shear flow of a dense assembly of dissipative disks using discrete simulation and prescribing the pressure and the shear rate. Those shear states are steady and uniform, and become intermittent in the quasistatic regime. In the limit of rigid grains, the shear state is determined by a single dimensionless number, called the inertial number I , which describes the ratio of inertial to pressure forces. Small values of I correspond to the quasistatic critical state of soil mechanics, while large values of I correspond to the fully collisional regime of kinetic theory. When I increases in the intermediate dense flow regime, we measure an approximately linear decrease of the solid fraction from the maximum packing value, and an approximately linear increase of the effective friction coefficient from the static internal friction value. From those dilatancy and friction laws, we deduce the constitutive law for dense granular flows, with a plastic Coulomb term and a viscous Bagnold term. The mechanical characteristics of the grains (restitution, friction, and elasticity) have a small influence in the dense flow regime. Finally, we show that the evolution of the relative velocity fluctuations and of the contact force anisotropy as a function of I provides a simple explanation of the friction law.

Devil's-staircase-like behavior of the range of random time series with record-breaking fluctuations

We propose insight into the analysis of the record-breaking fluctuations in random time series, which permits to distinguish between the self-organized criticality and the record dynamics (RD) scenarios of system evolution, using a finite time series realization. Performed analysis of the time series associated with the historical prices of different commodities has shown that the evolution of commodity markets is controlled by the record-breaking fluctuations as it is outlined by the RD. Furthermore, we found that the sizes of record-breaking fluctuations follow a fat-tailed distribution and the devil's-staircase-like records of price ranges are multiaffine and persistent, nevertheless, the high moments (q> q(C) >2) of their q-order height-height correlation functions behave logarithmically.

Intermittency and clustering in a system of self-driven particles

Intermittent behavior is shown to appear in a system of self-driven interacting particles. In the ordered phase, most particles move in the same approximate direction, but the system displays a series of intermittent bursts during which the order is temporarily lost. This intermittency is characterized and its statistical properties are found analytically for a reduced system containing only two particles. For large systems, the particles aggregate into clusters that play an essential role in the intermittent dynamics. The study of the cluster statistics shows that both the cluster sizes and the transition probability between them follow power-law distributions. The exchange of particles between clusters is shown to satisfy detailed balance.

Contact dynamics in a gently vibrated granular pile

We use multispeckle diffusive wave spectroscopy to probe the micron-scale dynamics of a water-saturated granular pile submitted to discrete gentle taps. The typical time scale between plastic events is found to increase dramatically with the number of applied taps. Furthermore, this microscopic dynamics weakly depends on the solid fraction of the sample. This process is largely analogous to the aging phenomenon observed in thermal glassy systems. We propose a heuristic model where this slowing-down mechanism is associated with a slow evolution of the distribution of the contact forces between particles. This model accounts for the main features of the observed dynamics.

Granular medium under vertical tapping: change of compaction and convection dynamics around the liftoff threshold

A granular medium submitted to vertical tapping reveals simultaneously compaction and convection. The two phenomena are directly coupled and their dynamics can be quantified by a characteristic compaction time and by an estimation of the convective downhill speed along the wall. A remarkable change of behavior is observed around the liftoff acceleration threshold of the whole packing, with a drastic slowing down of both dynamics below this threshold. Above it, a collective shock wave densifies the packing at each tap, whereas, below it, cumulative localized rearrangements will compact the entire system in the long time range.

Order parameter description of stationary partially fluidized shear granular flows

We carry out a detailed comparison of soft-particle molecular dynamics simulations with the theory of partially fluidized shear granular flows. We verify by direct simulations a constitutive relation based on the separation 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 components is determined by the order parameter. Based on results of the simulations we construct the "free energy" function for the order parameter. We also present the simulations of the stationary deep 2D granular flows driven by an upper wall and compare it with the continuum theory.