Analysis of a dislocation model for earthquakes

Jing Liu,Xin Yang Lu,Noel J Walkington

2024-09-24

Abstract:Approximation of problems in linear elasticity having small shear modulus in a thin region is considered. Problems of this type arise when modeling ground motion due to earthquakes where rupture occurs in a thin fault. It is shown that, under appropriate scaling, solutions of these problems can be approximated by solutions of a limit problem where the fault region is represented by a surface. In a numerical context this eliminates the need to resolve the large deformations in the fault; a numerical example is presented to illustrate efficacy of this strategy.

Numerical Analysis

What problem does this paper attempt to address?

### What problems does this paper attempt to solve?

This paper aims to solve the problem of mathematical modelling for crustal movement simulation in earthquakes, especially the challenges faced by direct numerical simulation when the fault zone (i.e., the thin zone where large deformations occur) is very small. Specifically, the authors explore how to simplify these problems into a limit problem as the width of the thin fault zone tends to zero, thus avoiding the need for fine meshing of the fault zone.

#### Main problems and solutions

1. **Problem description**:
- In earthquake simulation, ground motion is usually represented by underground structures made of linearly elastic materials, with thin fault zones between these materials where large deformations (ruptures) occur.
- Directly simulating these models requires very fine meshes to resolve the large deformations within the fault zone, which is computationally very expensive and complex.

2. **Solution**:
- Through appropriate scaling, the fault zone can be approximated as a surface, and the rupture can be represented as a discontinuity in certain components of the solution.
- The key to this method lies in proving that as the fault width tends to zero, the solution of the original problem converges to the solution of this simplified model.

3. **Theoretical basis**:
- The paper uses Gamma - convergence techniques to prove that as the fault width approaches zero, the energy functional of the original problem converges to the energy functional of a simplified limit problem.
- Specifically, for a given shear modulus \( \mu \) and bulk modulus \( \lambda \), as well as the permanent deformation \( \gamma \) within the fault zone, the energy functional \( I_\epsilon(u) \) converges to the limit energy functional \( I(u) \), where:

\[
I(u) = \frac{1}{2} \int_{\Omega \setminus S_0} \left( 2\mu |D(u)|^2 + \lambda (\text{div}(u))^2 \right) dx + \frac{1}{2} \int_{S_0} \mu ([u_1] - \gamma)^2 dx
\]

- Here \( [u_1(x)] \) represents the jump of the first (horizontal) component of \( u \) on the line \( S_0 = (-1,1) \times \{0\} \).

4. **Numerical verification**:
- The article also provides numerical examples to demonstrate the effectiveness of this simplified model. For example, in Figure 1, it is unnecessary to use a fine mesh to resolve the deformation; and in Figure 2, it is more efficient to represent the shear deformation on both sides of the fault by the discontinuity in the horizontal displacement.

Through this method, the paper provides an effective way to handle the fault zone in earthquake simulation, greatly reducing the computational cost and improving the efficiency of numerical simulation.