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.

Arbitrary order spline representation of cohomology generators for isogeometric analysis of eddy current problems

A Low-Frequency-Stable Higher-Order Isogeometric Discretization of the Augmented Electric Field Integral Equation

Isogeometric Analysis of Electric Field Integral Equation on Multipatch NURBS Surfaces with Discontinuous Galerkin

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

Low-Frequency Stabilization for the B-Spline Based Isogeometric Discretization of the Electric Field Integral Equation

A Novel Approach to Deal with Rotationally Symmetrical Conditions for 3 D Eddy Current Field Problems

Isogeometric Analysis for the Scattering of NURBS Surfaces with Coinciding Knots

Branched splines

Recent Advances of Isogeometric Analysis in Computational Electromagnetics

An arbitrary high-order Spectral Difference method for the induction equation

High order approximation of Hodge Laplace problems with local coderivatives on cubical meshes

A general class of $C^1$ smooth rational splines: Application to construction of exact ellipses and ellipsoids

Spectral Analysis of High Order Continuous FEM for Hyperbolic PDEs on Triangular Meshes: Influence of Approximation, Stabilization, and Time-Stepping

A high order discontinuous Galerkin method for the symmetric form of the anisotropic viscoelastic wave equation

Isogeometric analysis for solving discontinuous two-phase engineering problems with precise and explicit interface representation

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

A new spline method on graded mesh for fourth-order time-dependent PDEs: Application to Kuramoto-Sivashinsky and extended Fisher- Kolmogorov equations

Edge Functions for Spectral Element Methods

Multiply Periodic Splines

Discontinuous Galerkin Based Isogeometric Analysis for Geometric Flows and Applications in Geometric Modeling.

Discrete Differential Geometry for $C^{1,1}$ Hyperbolic Surfaces of Non-Constant Curvature