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.

Goodbye Christoffel Symbols: A Flexible and Efficient Approach for Solving Physical Problems in Curved Spaces

A Three-Dimensional Finite Volume Method for Conservation Laws in Conjunction with Modified Solution for Nonlinear Coupled Constitutive Relations

A Geometric Approach to the Navier-Stokes Equations

Numerical approximation of boundary value problems for curvature flow and elastic flow in Riemannian manifolds

Fluid Dynamics in Curvilinear Coordinates without Fictitious Forces

A Conservative Numerical Method for the Cahn-Hilliard Equation with Generalized Mobilities on Curved Surfaces in Three-Dimensional Space

A simple diffuse interface approach for compressible flows around moving solids of arbitrary shape based on a reduced Baer-Nunziato model

Curvature-driven flows: a variational approach

A Genuinely Two-Dimensional Riemann Solver for Compressible Flows in Curvilinear Coordinates

A novel finite element approximation of anisotropic curve shortening flow

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

Algorithmic differentiation of hyperbolic flow problems

Curvilinear Virtual Elements for 2D solid mechanics applications

A Cartesian grid-based boundary integral method for moving interface problems

Conformal approach to physics simulations for thin curved 3D membranes

A high‐order compact scheme for solving the 2D steady incompressible Navier‐Stokes equations in general curvilinear coordinates

Numerical Relativity in Spherical Polar Coordinates: Off-center Simulations

A Tetrad-First Approach to Robust Numerical Algorithms in General Relativity

An Approach for Solving Three-Dimensional Fluid Dynamics Problems with Allowance for Elastic Processes

Aerodynamics Simulations of Three-Dimensional Inviscid Flow Using Curvilinear Discontinuous Galerkin Method on Unstructured Meshes

A Finite Element Method to Compute Three-Dimensional Equilibrium Configurations of Fluid Membranes: Optimal Parameterization, Variational Formulation and Applications.