TY - GEN
T1 - Optimization and local time stepping of an ader-dg scheme for fully anisotropic wave propagation in complex geometries
AU - Wolf, Sebastian
AU - Gabriel, Alice Agnes
AU - Bader, Michael
N1 - KAUST 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.
PY - 2020/6/15
Y1 - 2020/6/15
N2 - 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.
AB - 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.
UR - http://hdl.handle.net/10754/679497
UR - http://link.springer.com/10.1007/978-3-030-50420-5_3
UR - http://www.scopus.com/inward/record.url?scp=85087274573&partnerID=8YFLogxK
U2 - 10.1007/978-3-030-50420-5_3
DO - 10.1007/978-3-030-50420-5_3
M3 - Conference contribution
SN - 9783030504199
SP - 32
EP - 45
BT - Lecture Notes in Computer Science
PB - Springer International Publishing
ER -