Abstract:High order accurate and explicit time-stable solvers are well suited for hyperbolic wave propagation problems. As a result of the complexities of real geometries, internal interfaces and nonlinear boundary and interface conditions, discontinuities and sharp wave fronts may become fundamental features of the solution. Thus, geometrically flexible and adaptive numerical algorithms are critical for high fidelity and efficient simulations of wave phenomena in many applications. Adaptive curvilinear meshes hold promise to minimise the effort to represent complicated geometries or heterogeneous material data avoiding the bottleneck of feature-preserving meshing. To enable the design of stable DG methods on three space dimensional (3D) curvilinear elements we construct a structure preserving anti-symmetric coordinate transformation motivated by the underlying physics. Using a physics-based numerical penalty-flux, we develop a 3D provably energy-stable discontinuous Galerkin finite element approximation of the elastic wave equation in geometrically complex and heterogenous media. By construction, our numerical flux is upwind and yields a discrete energy estimate analogous to the continuous energy estimate. The ability to treat conforming and non-conforming curvilinear elements allows for flexible adaptive mesh refinement strategies. The numerical scheme has been implemented in ExaHyPE, a simulation engine for parallel dynamically adaptive simulations of wave problems on adaptive Cartesian meshes. We present 3D numerical experiments of wave propagation in heterogeneous isotropic and anisotropic elastic solids demonstrating stability and high order accuracy. We demonstrate the potential of our approach for computational seismology in a regional wave propagation scenario in a geologically constrained 3D model including the geometrically complex free-surface topography of Mount Zugspitze, Germany.

Domain of Dependence stabilization for the acoustic wave equation on 2D cut-cell meshes

Error analysis of a first-order DoD cut cell method for 2D unsteady advection

An Energy Stable High-Order Cut Cell Discontinuous Galerkin Method with State Redistribution for Wave Propagation

Weight-adjusted discontinuous Galerkin methods: matrix-valued weights and elastic wave propagation in heterogeneous media

An Optimal Nearly Analytic Discrete-Weighted Runge-Kutta Discontinuous Galerkin Hybrid Method for Acoustic Wavefield Modeling

Bifurcation tracking on moving meshes and with consideration of azimuthal symmetry breaking instabilities

A space-time domain RBF method for 2D wave equations

A Parallel Weighted Runge-Kutta Discontinuous Galerkin Method For Solving Acousitc Wave Equations In 3d D'Alembert Media On Unstructured Meshes

An Effective Discontinuous Galerkin Method for Solving Acoustic Wave Equations in Heterogeneous Media

A stable discontinuous Galerkin method for linear elastodynamics in 3D geometrically complex media using physics based numerical fluxes

A Cut Discontinuous Galerkin Method for Coupled Bulk-Surface Problems

A Moving Boundary Flux Stabilization Method for Cartesian Cut-Cell Grids using Directional Operator Splitting

Domain decomposition for quasi-periodic scattering by layered media via robust boundary-integral equations at all frequencies

Stability estimates of Nyström discretizations of Helmholtz decomposition boundary integral equation formulations for the solution of Navier scattering problems in two dimensions with Dirichlet boundary conditions

Numerical Boundary Control of Multi-Dimensional Hyperbolic Equations

An Additive Overlapping Domain Decomposition Method for the Helmholtz Equation

An integral equation method for the advection-diffusion equation on time-dependent domains in the plane

An Energy Based Discontinuous Galerkin Method for Coupled Elasto-Acoustic Wave Equations in Second Order Form

Convergence of finite difference methods for the wave equation in two space dimensions

An iterative method for Helmholtz boundary value problems arising in wave propagation