Emergent fault friction and supershear in a continuum model of geophysical rupture

Abhishek Arora,Amit Acharya
2023-12-15
Abstract:Important physical observations in rupture dynamics such as static fault friction, short-slip, self-healing, and supershear phenomenon in cracks are studied. A continuum model of rupture dynamics is developed using the field dislocation mechanics (FDM) theory. The energy density function in our model encodes accepted and simple physical facts related to rocks and granular materials under compression. We work within a 2-dimensional ansatz of FDM where the rupture front is allowed to move only in a horizontal fault layer sandwiched between elastic blocks. Damage via the degradation of elastic modulus is allowed to occur only in the fault layer, characterized by the amount of plastic slip. The theory dictates the evolution equation of the plastic shear strain to be a Hamilton-Jacobi (H-J) equation, resulting in the representation of a propagating rupture front, which is fully coupled to elastodynamics in the whole domain. Our simulations recover static friction laws as emergent features of our model, without putting in by hand any such discontinuous criteria. Estimates of material parameters of cohesion and friction angle are deduced. Short-slip and slip-weakening (crack-like) behaviors are also reproduced as a function of the degree of damage behind the rupture front. The long-time behavior of a moving rupture front is probed, and it is deduced that the equilibrium profiles under no shear stress are not traveling wave profiles under non-zero shear load. However, it is shown that a traveling wave structure is likely attained in the limit of long times. Finally, a crack-like damage front is driven by an initial impact loading, and it is observed in our simulations that an upper bound to the crack speed is the dilatational wave speed of the material unless the material is put under pre-stressed conditions, in which case supersonic motion can be obtained.
Materials Science,Mesoscale and Nanoscale Physics
What problem does this paper attempt to address?
The paper aims to study important physical phenomena in the dynamics of geological fault rupture, such as static fault friction, short slip, self-healing, and supershear phenomena, by constructing a continuum model. Specifically, the research team developed a continuum model using the Field Dislocation Mechanics (FDM) theory to simulate the dynamic process of rupture. In this model, the energy density function encodes the fundamental physical properties of rocks and granular materials under compression. The study assumes that the rupture front can only move within the horizontal fault plane and that elastic modulus degradation, characterized by the degree of plastic slip, occurs only on this plane. According to theoretical derivation, the evolution equation of plastic shear strain is a Hamilton-Jacobi (H-J) equation, which allows the representation of the propagating rupture front. To solve the H-J equation, a Central-Upwind scheme was adopted. Through numerical simulations, the model can reproduce the static friction law as an intrinsic feature without artificially setting any discontinuous criteria. Additionally, the model can predict short slip and slip-weakening (crack-like) behavior, which depends on the variation in the degree of damage behind the rupture front. The long-term behavior of the moving rupture front was explored, and it was found that in the absence of shear stress, the equilibrium profile is not a traveling wave profile under non-zero shear stress. However, in the long-term limit, the model suggests that traveling wave structures may form. Finally, by driving the damage front with an initial impact load, it was observed that the upper limit of crack speed is the material's dilatational wave speed, unless the material is in a prestressed state, in which case supersonic motion can be achieved. In the absence of prestress, the supershear state can be restored. In summary, this study constructs a continuum model based on FDM theory, which not only explains and simulates key phenomena such as static friction but also explores other complex behaviors in rupture dynamics, providing a new perspective for understanding geological fault rupture.