Meshing scheme in the computation of spontaneous fault rupture process

Abstract

The choice of spatial grid size has been being a crucial issue in all kinds of numerical algorithms. By using BIEM (Boundary Integral Equation Method) to calculate the rupture process of a planar fault embedded in an isotropic and homogeneous full space with simple discretization scheme, this paper focuses on what grid size should be applied to control the error as well as maintaining the computing efficiency for different parameter combinations of (Dc, Te), where Dc is the critical slip-weakening distance and Te is the initial stress on the fault plane. We have preliminarily found the way of properly choosing the spatial grid size, which is of great significance in the computation of seismic source rupture process with BIEM.