Search

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly to simulate a sheared granular fault system. We focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements measured at locations on the upper and lower plates just adjacent to the gouge in the direction parallel to the shear direction (x direction). The simulations show that at slips the plate velocities are proportional to the normal load and may be inversely proportional to the square root of the plate's Young's modulus; whereas the driving shear velocity does not show distinct influence on the plate velocities. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than the half of the driving shear velocity, and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory scale fault gouge system and may help reveal the complexities of earthquake frictional dynamics.

When a dynamic earthquake rupture propagates on a fault in the Earth's crust, the medium around the fault is dynamically damaged due to stress concentrations around the rupture tip. Recent field observations, laboratory experiments and canonical numerical models show the coseismic off-fault damage is essential to describe the coseismic off-fault deformation, rupture dynamics, radiation and overall energy budget. However, the numerical modeling of "localized" off-fault fractures remains a challenge mainly because of computational limitations and model formulation shortcomings. We thus developed a numerical framework for modeling coseismic off-fault fracture networks using the combined finite-discrete element method (FDEM) and we applied it to simulate dynamic ruptures with coseismic off-fault damage on various fault configurations. This paper addresses the role of coseismic off-fault damage on rupture dynamics associated with a planar fault, as a base case, and with a number of first-order geometrical complexities, such as fault kink, step-over and roughness., Comment: arXiv admin note: text overlap with arXiv:1901.01771

Most earthquake ruptures propagate at speeds below the shear wave velocity within the crust, but in some rare cases, ruptures reach supershear speeds. The physics underlying the transition of natural subshear earthquakes to supershear ones is currently not fully understood. Most observational studies of supershear earthquakes have focused on determining which fault segments sustain fully-grown supershear ruptures. Experimentally cross-validated numerical models have identified some of the key ingredients required to trigger a transition to supershear speed. However, the conditions for such a transition in nature are still unclear, including the precise location of this transition. In this work, we provide theoretical and numerical insights to identify the precise location of such a transition in nature. We use fracture mechanics arguments with multiple numerical models to identify the signature of supershear transition in coseismic off-fault damage. We then cross-validate this signature with high-resolution observations of fault zone width and early aftershock distributions. We confirm that the location of the transition from subshear to supershear speed is characterized by a decrease in the width of the coseismic off-fault damage zone. We thus help refine the precise location of such a transition for natural supershear earthquakes.

A numerical scheme using the combined finite-discrete element methods (FDEM) is employed to study a model of an earthquake system comprising a granular layer embedded in a formation. When the formation is driven so as to shear the granular layer, a system of stress chains emerges. The stress chains endow the layer with resistance to shear and on failure launch broadcasts into the formation. These broadcasts, received as acoustic emission, provide a remote monitor of the state of the granular layer, of the earthquake system.

Earthquake ruptures dynamically activate coseismic off-fault damage around fault cores. Systematic field observation efforts have shown the distribution of off-fault damage around main faults, while numerical modeling using elastic-plastic off-fault material models has demonstrated the evolution of coseismic off-fault damage during earthquake ruptures. Laboratory scale micro-earthquake experiments have pointed out the enhanced high-frequency radiation due to the coseismic off-fault damage. However, the detailed off-fault fracturing mechanisms, subsequent radiation and its contribution to the overall energy budget remain to be fully understood because of limitations of current observational techniques and model formulations. Here, we constructed a new physics-based dynamic earthquake rupture modeling framework, based on the combined finite-discrete element method (FDEM), to investigate the fundamental mechanisms of coseismic off-fault damage, and its effect on the rupture dynamics, the radiation and the overall energy budget. We conducted a 2-D systematic case study with depth and showed the mechanisms of dynamic activation of the coseismic off-fault damage. We found the decrease in rupture velocity and the enhanced high-frequency radiation in near-field due to the coseismic off-fault damage. We then evaluated the overall energy budget, which shows a significant contribution of the coseismic off-fault damage to the overall energy budget even at depth, where the damage zone width becomes narrower. The present numerical framework for the dynamic earthquake rupture modeling thus provides the insight into the earthquake rupture dynamics with the coseismic off-fault damage.

While the hydraulic fracturing technology, aka fracking (or fraccing, frac), has become highly developed and astonishingly successful, a consistent formulation of the associated fracture mechanics that would not conflict with some observations is still unavailable. It is attempted here. Classical fracture mechanics, as well as the current commercial softwares, predict vertical cracks to propagate without branching from the perforations of the horizontal well casing, which are typically spaced at 10 m or more. However, to explain the gas production rate at the wellhead, the crack spacing would have to be only about 0.1 m, which would increase the overall gas permeability of shale mass about 10,000$\times$. This permeability increase has generally been attributed to a preexisting system of orthogonal natural cracks, whose spacing is about 0.1 m. But their average age is about 100 million years, and a recent analysis indicated that these cracks must have been completely closed by secondary creep of shale in less than a million years. Here it is considered that the tectonic events that produced the natural cracks in shale must have also created weak layers with nano- or micro-cracking damage. It is numerically demonstrated that a greatly enhanced permeability along the weak layers, with a greatly increased transverse Biot coefficient, must cause the fracking to engender lateral branching and the opening of hydraulic cracks along the weak layers, even if these cracks are initially almost closed. A finite element crack band model, based on recently developed anisotropic spherocylindrical microplane constitutive law, demonstrates these findings.

Numerical simulation of nonlinear elastic wave propagation in solids with cracks is indispensable for decoding the complicated mechanisms associated with the nonlinear ultrasonic techniques in Non-Destructive Testing (NDT). Here, we introduce a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), to explicitly simulate the crack induced nonlinear elasticity in solids with both horizontal and inclined cracks. In the FDEM model, the solid is discretized into finite elements to capture the wave propagation in the bulk material, and the finite elements along the two sides of the crack also behave as discrete elements to track the normal and tangential interactions between crack surfaces. The simulation results show that for cracked models, nonlinear elasticity is generated only when the excitation amplitude is large enough to trigger the contact between crack surfaces, and the nonlinear behavior is very sensitive to the crack surface contact. The simulations reveal the influence of normal and tangential contact on the nonlinear elasticity generation. Moreover, the results demonstrate the capabilities of FDEM for decoding the causality of nonlinear elasticity in cracked solid and its potential to assist in Non-Destructive Testing (NDT).

We propose a machine learning approach to address a key challenge in materials science: predicting how fractures propagate in brittle materials under stress, and how these materials ultimately fail. Our methods use deep learning and train on simulation data from high-fidelity models, emulating the results of these models while avoiding the overwhelming computational demands associated with running a statistically significant sample of simulations. We employ a graph convolutional network that recognizes features of the fracturing material and a recurrent neural network that models the evolution of these features, along with a novel form of data augmentation that compensates for the modest size of our training data. We simultaneously generate predictions for qualitatively distinct material properties. Results on fracture damage and length are within 3% of their simulated values, and results on time to material failure, which is notoriously difficult to predict even with high-fidelity models, are within approximately 15% of simulated values. Once trained, our neural networks generate predictions within seconds, rather than the hours needed to run a single simulation.

Fracture coalescence is a critical phenomenon for creating large fractures from smaller flaws, affecting fracture network flow and seismic energy release potential. In this paper, simulations of fracture coalescence processes in granite specimens with pre-existing cracks are performed. These simulations utilize an in-house implementation of the Combined Finite-Discrete Element method (FDEM) known as the Hybrid Optimization Software Suite (HOSS). The pre-existing cracks within the specimens follow two geometric patterns: 1) a single crack oriented at different angles with respect to the loading direction, and 2) two cracks, where one crack is oriented perpendicular to the loading direction and the other crack is oriented at different angles. The intent of this study is to demonstrate the suitability of FDEM for modeling fracture coalescence processes including: crack initiation and propagation, tensile and shear fracture behavior, and patterns of fracture coalescence. The simulations provide insight into the evolution of fracture tensile and shear fracture behavior as a function of time. The single-crack simulations accurately reproduce experimentally measured peak stresses as a function of crack inclination angle. Both the single- and double-crack simulations exhibit a linear increase in strength with increasing crack angle; the double-crack specimens are systematically weaker than the single-crack specimens.

In this paper, the Combined Finite-Discrete Element Method (FDEM) has been applied to analyze the deformation of anisotropic geomaterials. In the most general case geomaterials are both non-homogeneous and non-isotropic. With the aim of addressing anisotropic material problems, improved 2D FDEM formulations have been developed. These formulations feature the unified hypo-hyper elastic approach combined with a multiplicative decomposition-based selective integration for volumetric and shear deformation modes. This approach is significantly different from the co-rotational formulations typically encountered in finite element codes. Unlike the co-rotational formulation, the multiplicative decomposition-based formulation naturally decomposes deformation into translation, rotation, plastic stretches, elastic stretches, volumetric stretches, shear stretches, etc. This approach can be implemented for a whole family of finite elements from solids to shells and membranes. This novel 2D FDEM based material formulation was designed in such a way that the anisotropic properties of the solid can be specified in a cell by cell basis, therefore enabling the user to seed these anisotropic properties following any type of spatial variation, for example, following a curvilinear path. In addition, due to the selective integration, there are no problems with volumetric or shear locking with any type of finite element employed.

Large continental earthquakes activate multiple faults in a complex fault system, dynamically inducing co-seismic damage around them. The 2016 Mw 7.8 Kaikoura earthquake in the northern South Island of New Zealand has been reported as one of the most complex continental earthquakes ever documented1, which resulted in a distinctive on and off-fault deformation pattern. Previous geophysical studies confirm that the rupture globally propagated northward from epicenter. However, the exact rupture- propagation path is still not well understood because of the geometrical complexity, partly at sea, and the possibility of a blind thrust. Here we use a combination of state-of- the-art observation of surface deformation, provided by optical image correlation, and first principle physics-based numerical modeling to determine the most likely rupture path. We quantify in detail the observed horizontal co-seismic deformation and identify specific off-fault damage zones in the area of the triple junction between the Jordan, the Kekerengu and the Papatea fault segments. We also model dynamic rupture propagation, including the activation of off-fault damage, for two alternative rupture scenarios through the fault triple junction. Comparing our observations with the results from the above two modeled scenarios we show that only one of the scenarios best explains both the on and off-fault deformation fields. Our results provide a unique insight into the rupture pathway, by observing, and modeling, both on and off-fault deformation. We propose this combined approach here to narrow down the possible rupture scenarios for large continental earthquakes accompanied by co-seismic off-fault damage., Comment: 32 pages, 11 figures

Sheared granular layers undergoing stick slip behavior are broadly employed to study the physics and dynamics of earthquakes. Here, a two dimensional implementation of the combined finite discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular fault system including both gouge and plate, and to investigate the influence of different normal loads on seismic moment, macroscopic friction coefficient, kinetic energy, gouge layer thickness, and recurrence time between slips. In the FDEM model, the deformation of plates and particles is simulated using the FEM formulation while particle particle and particle plate interactions are modeled using DEM derived techniques. The simulated seismic moment distributions are generally consistent with those obtained from the laboratory experiments. In addition, the simulation results demonstrate that with increasing normal load, (i) the kinetic energy of the granular fault system increases; (ii) the gouge layer thickness shows a decreasing trend; and (iii) the macroscopic friction coefficient does not experience much change. Analyses of the slip events reveal that, as the normal load increases, more slip events with large kinetic energy release and longer recurrence time occur, and the magnitude of gouge layer thickness decrease also tends to be larger; while the macroscopic friction coefficient drop decreases. The simulations not only reveal the influence of normal loads on the dynamics of sheared granular fault gouge, but also demonstrate the capabilities of FDEM for studying stick slip dynamic behavior of granular fault systems.