9:15 AM - 9:30 AM
[STT43-02] Acceleration of dynamic rupture simulation using FDP=LH-matrices
Keywords:Boundary integral equation method, Dynamic rupture simulation
Ando(2016) suggested Fast Domain Partitioning (FDP) Method which partition causality cone the region where elastic waves affect into P and S wave fronts (Domain Fp, FS), the domain between P and S waves (Domain I), and the domain of the static equilibrium (Domain S). Each domain shows the different nature of kernel and this allows for fast computation.
Sato and Ando(2021) developed a new approximation algorithm called FDP=H-matrices and they acheived computational complexity O(NlogNM) by incorporating H matrices into FDPM with the use of travel-time approximation. They solved 2D elastodynamic wave propagation problems and confirmed its utilities.
In this research, we slightly changed the algorithm of Domain I not to use approximation and applied to 3D dynamic rupture simulation. Furthermore, instead of using H matrices, we adopted Lattice H matrices(Ida 2018) as low-rank approximation methods. Lattice H matrices are constructed by blocks of lattice structures which consist of H matrices as submatrices in order to efficiently use a large number of process. For simplicity, we set the fault plane in parallel and investigated the accuracy of calculation and the spped of the simulation. The results indicate the potential to apply this FDP=LH-matrices method to actual earthquake simulations.