Optimization and local time stepping of an ader-dg scheme for fully anisotropic wave propagation in complex geometries

Sebastian Wolf, Alice Agnes Gabriel, Michael Bader

Research output: Chapter in Book/Report/Conference proceedingConference contribution

8 Scopus citations


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 languageEnglish (US)
Title of host publicationLecture Notes in Computer Science
PublisherSpringer International Publishing
Number of pages14
ISBN (Print)9783030504199
StatePublished - Jun 15 2020
Externally publishedYes

Bibliographical note

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.


Dive into the research topics of 'Optimization and local time stepping of an ader-dg scheme for fully anisotropic wave propagation in complex geometries'. Together they form a unique fingerprint.

Cite this