Three-dimensional discrete element modeling of triggered slip in sheared granular media

A Discrete Elements Study of the Frictional Behavior of Fault Gouges

A series of discrete elements simulations is presented for the study of fault gouges' frictional response. The gouge is considered to have previously undergone ultra-cataclastic flow and long-time consolidation loading. We explore the effect of different particle characteristics such as size, polydispersity, and also shearing velocities on gouge's response under the conditions met in the seismogenic zone. Monte-Carlo analyses suggest that the local stick-slip events disappear when averaging over a large number of numerical samples. Moreover, the apparent material frictional response remains almost unaffected by the spatial randomness of particles' position and by the particle's size distribution. On the contrary, the mean particle size controls the formation and thickness of the observed shear bands, which appear after the peak friction is met. Furthermore, the apparent friction evolution fits well to an exponential decay law with slip, which involves a particle size dependent critical slip distance. For the studied conditions and depth, the shearing velocity is found to play a secondary role on the apparent frictional response of the gouge, which highlights the importance of analyses involving multiphysics for studying the rheology of fault gouges. Besides improving the understanding of the underlying physics of the problem, the above findings are also useful for deriving pertinent constitutive models in the case of modeling with continuum theories.

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.

Shear bands in dense fault gouge

Earthquakes happen with frictional sliding, by releasing all the stresses accumulated in the prestressed surrounding medium. The geological fault gouge, coming from the wear of previous slips, acts on friction stability and plays a key role in this sudden energy release. A large part of slip mechanisms are influenced, if not controlled, by the characteristics and environment of this tribological “third body”. A 2D granular fault (mm scale) is implemented with Discrete Element Modelling (DEM). A displacement-driven model with dry contact is studied to observe kinematics and properties of the slipping zone. Increasing the length of the granular media increases the slip needed to weaken the friction from friction peak to steadystate. Low-angle Riedel shear bands are mostly observed. Their number increases with the inter-particle friction coefficient, which also influences shear bands formation in their orientation angle (higher friction leads to higher angle with the main slip direction).

Induced and endogenous acoustic oscillations in granular faults

The frictional properties of disordered systems are affected by external perturbations. These perturbations usually weaken the system by reducing the macroscopic friction coefficient. This friction reduction is of particular interest in the case of disordered systems composed of granular particles confined between two plates, as this is a simple model of seismic fault. Indeed, in the geophysical context frictional weakening could explain the unexpected weakness of some faults, as well as earthquake remote triggering. In this manuscript, we review recent results concerning the response of confined granular systems to external perturbations, considering the different mechanisms by which the perturbation could weaken a system, the relevance of the frictional reduction to earthquakes, as well as discussing the intriguing scenario whereby the weakening is not monotonic in the perturbation frequency, so that a re-entrant transition is observed, as the system first enters a fluidized state and then returns to a frictional state.This article is part of the theme issue 'Statistical physics of fracture and earthquakes'.

Strain localization in dry sheared granular materials: A compactivity-based approach

Shear banding is widely observed in natural fault zones as well as in laboratory experiments on granular materials. Understanding the dynamics of strain localization under different loading conditions is essential for quantifying strength evolution of fault gouge and energy partitioning during earthquakes and characterizing rheological transitions and fault zone structure changes. To that end, we develop a physics-based continuum model for strain localization in sheared granular materials. The grain-scale dynamics is described by the shear transformation zone (STZ) theory, a nonequilibrium statistical thermodynamic framework for viscoplastic deformation in amorphous materials. Using a finite strain computational framework, we investigate the initiation and growth of complex shear bands under a variety of loading conditions and identify implications for strength evolution and the ductile to brittle transition. Our numerical results show similar localization patterns to field and laboratory observations and suggest that shear zones show more ductile response at higher confining pressures, lower dilatancy, and loose initial conditions. Lower pressures, higher dilatancy, and dense initial conditions favor a brittle response and larger strength drops. These findings shed light on a range of mechanisms for strength evolution in dry sheared granular materials and provide a critical input to physics-based multiscale models of fault zone instabilities.

Controlled Viscosity in Dense Granular Materials

We experimentally investigate the fluidization of a granular material subject to mechanical vibrations by monitoring the angular velocity of a vane suspended in the medium and driven by an external motor. On increasing the frequency, we observe a reentrant transition, as a jammed system first enters a fluidized state, where the vane rotates with high constant velocity, and then returns to a frictional state, where the vane velocity is much lower. While the fluidization frequency is material independent, the viscosity recovery frequency shows a clear dependence on the material that we rationalize by relating this frequency to the balance between dissipative and inertial forces in the system. Molecular dynamics simulations well reproduce the experimental data, confirming the suggested theoretical picture.

Synchronized oscillations and acoustic fluidization in confined granular materials

According to the acoustic fluidization hypothesis, elastic waves at a characteristic frequency form inside seismic faults even in the absence of an external perturbation. These waves are able to generate a normal stress which contrasts the confining pressure and promotes failure. Here, we study the mechanisms responsible for this wave activation via numerical simulations of a granular fault model. We observe the particles belonging to the percolating backbone, which sustains the stress, to perform synchronized oscillations over ellipticlike trajectories in the fault plane. These oscillations occur at the characteristic frequency of acoustic fluidization. As the applied shear stress increases, these oscillations become perpendicular to the fault plane just before the system fails, opposing the confining pressure, consistently with the acoustic fluidization scenario. The same change of orientation can be induced by external perturbations at the acoustic fluidization frequency.

Rattler-induced aging dynamics in jammed granular systems

Granular materials jam when developing a network of contact forces able to resist the applied stresses. Through numerical simulations of the dynamics of the jamming process, we show that the jamming transition does not occur when the kinetic energy vanishes. Rather, as the system jams, the kinetic energy becomes dominated by rattler particles, which scatter within their cages. The relaxation of the kinetic energy in the jammed configuration exhibits a double power-law decay, which we interpret in terms of the interplay between backbone and rattler particles.

Cascading elastic perturbation in Japan due to the 2012 M w 8.6 Indian Ocean earthquake

Since the discovery of extensive earthquake triggering occurring in response to the 1992 M w (moment magnitude) 7.3 Landers earthquake, it is now well established that seismic waves from earthquakes can trigger other earthquakes, tremor, slow slip, and pore pressure changes. Our contention is that earthquake triggering is one manifestation of a more widespread elastic disturbance that reveals information about Earth's stress state. Earth's stress state is central to our understanding of both natural and anthropogenic-induced crustal processes. We show that seismic waves from distant earthquakes may perturb stresses and frictional properties on faults and elastic moduli of the crust in cascading fashion. Transient dynamic stresses place crustal material into a metastable state during which the material recovers through a process termed slow dynamics. This observation of widespread, dynamically induced elastic perturbation, including systematic migration of offshore seismicity, strain transients, and velocity transients, presents a new characterization of Earth's elastic system that will advance our understanding of plate tectonics, seismicity, and seismic hazards.

Dynamic Weakening by Acoustic Fluidization during Stick-Slip Motion

The unexpected weakness of some faults has been attributed to the emergence of acoustic waves that promote failure by reducing the confining pressure through a mechanism known as acoustic fluidization, also proposed to explain earthquake remote triggering. Here we validate this mechanism via the numerical investigation of a granular fault model system. We find that the stick-slip dynamics is affected only by perturbations applied at a characteristic frequency corresponding to oscillations normal to the fault, leading to gradual dynamical weakening as failure is approaching. Acoustic waves at the same frequency spontaneously emerge at the onset of failure in the absence of perturbations, supporting the relevance of acoustic fluidization in earthquake triggering.

Stick-slip instabilities in sheared granular flow: The role of friction and acoustic vibrations

We propose a theory of shear flow in dense granular materials. A key ingredient of the theory is an effective temperature that determines how the material responds to external driving forces such as shear stresses and vibrations. We show that, within our model, friction between grains produces stick-slip behavior at intermediate shear rates, even if the material is rate strengthening at larger rates. In addition, externally generated acoustic vibrations alter the stick-slip amplitude, or suppress stick-slip altogether, depending on the pressure and shear rate. We construct a phase diagram that indicates the parameter regimes for which stick-slip occurs in the presence and absence of acoustic vibrations of a fixed amplitude and frequency. These results connect the microscopic physics to macroscopic dynamics and thus produce useful information about a variety of granular phenomena, including rupture and slip along earthquake faults, the remote triggering of instabilities, and the control of friction in material processing.

Softening of stressed granular packings with resonant sound waves

We perform numerical simulations of a two-dimensional bidisperse granular packing subjected to both a static confining pressure and a sinusoidal dynamic forcing applied by a wall on one edge of the packing. We measure the response experienced by a wall on the opposite edge of the packing and obtain the resonant frequency of the packing as the static or dynamic pressures are varied. Under increasing static pressure, the resonant frequency increases, indicating a velocity increase of elastic waves propagating through the packing. In contrast, when the dynamic amplitude is increased for fixed static pressure, the resonant frequency decreases, indicating a decrease in the wave velocity. This occurs both for compressional and for shear dynamic forcing and is in agreement with experimental results. We find that the average contact number Zc at the resonant frequency decreases with increasing dynamic amplitude, indicating that the elastic softening of the packing is associated with a reduced number of grain-grain contacts through which the elastic waves can travel. We image the excitations created in the packing and show that there are localized disturbances or soft spots that become more prevalent with increasing dynamic amplitude. Our results are in agreement with experiments on glass bead packings and earth materials such as sandstone and granite and may be relevant to the decrease in elastic wave velocities that has been observed to occur near fault zones after strong earthquakes, in surficial sediments during strong ground motion, and in structures during earthquake excitation.