FD3D_TSN: A Fast and Simple Code for Dynamic Rupture Simulations with GPU Acceleration

Premus, Jan; Gallovic, Frantisek; Hanyk, Ladislav and Gabriel, Alice-Agnes (2020): FD3D_TSN: A Fast and Simple Code for Dynamic Rupture Simulations with GPU Acceleration. In: Seismological Research Letters, Vol. 91, No. 5: pp. 2881-2889

Abstract

We introduce FD3D_TSN-an open-source Fortran code for 3D dynamic earthquake rupture modeling based on the staggered grid fourth-order finite-difference method employing a regular cubical spatial discretization. Slip-weakening and fast-velocity weakening rate-and-state fault friction laws are combined with vertical planar fault geometry, orthogonal to a planar free surface. FD3D_TSN demonstrates good agreement with other methods in a range of benchmark exercises of the Southern California Earthquake Center and U.S. Geological Survey dynamic rupture code verification project. Efficient graphic processing units (GPU) acceleration using the OpenACC framework yields a factor of 10 speed-up in terms of time to solution compared to a single core solution for current hardware (Intel i9-9900K and Nvidia RTX 2070). The software is fast and easy-to-use and suitable explicitly for data-driven applications requiring a large number of forward simulations such as dynamic source inversion or probabilistic ground-motion modeling. The code is freely available for the scientific community and may be incorporated in physics-based earthquake source imaging and seismic hazard assessment, or for teaching purposes.

• Item Type: Journal article
• Faculties: Geosciences > Department of Earth and Environmental Sciences
• Subjects: 500 Science > 550 Earth sciences and geology
• ISSN: 0895-0695
• Language: English