We present an extension of the earthquake simulation software SeisSol to support seismic wave propagation in fully triclinic anisotropic materials. To our best knowledge, SeisSol is one of the few open-source codes that offer this feature for simulations at petascale performance and beyond. We employ a Discontinuous Galerkin (DG) method with arbitrary high-order derivative (ADER) time stepping. Here, we present a novel implementation of fully physical anisotropy with a two-sided Godunov flux and local time stepping. We validate our implementation on various benchmarks and present convergence analysis with respect to analytic solutions. An application example of seismic waves scattering around the Zugspitze in the Bavarian Alps demonstrates the capabilities of our implementation to solve geophysics problems fast.
|Original language||English (US)|
|Title of host publication||Lecture Notes in Computer Science|
|Publisher||Springer International Publishing|
|Number of pages||14|
|State||Published - Jun 15 2020|
Bibliographical noteKAUST Repository Item: Exported on 2022-06-30
Acknowledged KAUST grant number(s): ORS-2017-CRG6 3389.02
Acknowledgements: The research leading to these results has received funding from European Union Horizon 2020 research and innovation program (ENERXICO, grant agreement No. 828947), KAUST-CRG (FRAGEN, grant no. ORS-2017-CRG6 3389.02) and the European Research Council (TEAR, ERC Starting grant no. 852992). Computing resources were provided by the Leibniz Supercomputing Centre (project no. pr45fi on SuperMUC-NG).
This publication acknowledges KAUST support, but has no KAUST affiliated authors.