3D Acoustic-Elastic Coupling with Gravity: The Dynamics of the 2018 Palu, Sulawesi Earthquake and Tsunami

by   Lukas Krenz, et al.

We present a highly scalable 3D fully-coupled Earth ocean model of earthquake rupture and tsunami generation. We model seismic, acoustic and surface gravity wave propagation in elastic (Earth) and acoustic (ocean) materials sourced by physics-based non-linear earthquake dynamic rupture. Complicated geometries, including high-resolution bathymetry, coastlines and segmented earthquake faults are discretized by adaptive unstructured tetrahedral meshes. A Discontinuous Galerkin discretization with ADER local time-stepping (ADER-DG) yields petascale computational efficiency and high-order accuracy in time and space. We compare the 3D fully-coupled approach to a benchmark problem for 3D-2D linked models that use 2D shallow-water modeling. We present a large-scale fully-coupled model of the 2018 Sulawesi events that links the dynamics from supershear earthquake faulting to elastic and acoustic waves in Earth and ocean to tsunami gravity wave propagation in the narrow Palu Bay. And we demonstrate scalability and performance of the MPI+OpenMP parallelization on three petascale supercomputers.


page 1

page 10


Simulation of 3D elasto-acoustic wave propagation based on a Discontinuous Galerkin Spectral Element method

In this paper we present a numerical discretization of the coupled elast...

Shifted Laplacian multigrid for the elastic Helmholtz equation

The shifted Laplacian multigrid method is a well known approach for prec...

RBF-FD analysis of 2D time-domain acoustic wave propagation in heterogeneous Earth's subsurface

Radial Basis Function-generated Finite Differences (RBF-FD) is a popular...

A high-order discontinuous Galerkin method for the poro-elasto-acoustic problem on polygonal and polyhedral grids

The aim of this work is to introduce and analyze a finite element discon...

Two-grid method on unstructured tetrahedra: Applying computational geometry to staggered solution of coupled flow and mechanics problems

We develop a computational framework that leverages the features of soph...

Please sign up or login with your details

Forgot password? Click here to reset