A simple diffuse interface approach on adaptive Cartesian grids for the linear elastic wave equations with complex topography

Abstract

In most classical approaches of computational geophysics for seismic wave propagation problems, complex surface topography is either accounted for by boundary-fitted unstructured meshes, or, where possible, by mapping the complex computational domain from physical space to a topologically simple domain in a reference coordinate system. In this paper we propose a completely different strategy. We address the problem of geometrically complex free surface boundary conditions with a novel diffuse interface method on adaptive Cartesian meshes that consists in the introduction of a characteristic function $0\le \alpha \le 1$ which identifies the location of the solid medium and the surrounding air and thus implicitly defines the location of the free surface boundary. Our new approach completely avoids the problem of mesh generation, since all that is needed for the definition of the complex surface topography is to set a scalar color function to unity inside the regions covered by the solid and to zero outside. An analysis of the eigenvalues of the PDE system shows that the complexity of the geometry has no influence on the admissible time step size due to the CFL condition. The model reduces to the classical linear elasticity equations inside the solid medium where the gradients of $\alpha$ are zero, while in the diffuse interface zone at the free surface boundary the governing PDE system becomes nonlinear. We can prove that the solution of the Riemann problem with arbitrary data and a jump in $\alpha$ from unity to zero yields a Godunov-state at the interface that satisfies the free-surface boundary condition exactly. In order to reduce numerical dissipation, we use high order DG finite element schemes on adaptive AMR grids together with a high resolution shock capturing subcell finite volume (FV) limiter in the diffuse interface region.

References in corpus (6)

• A pressure-based semi-implicit space-time discontinuous Galerkin method on staggered unstructured meshes for the solution of the compressible Navier-Stokes equations at all Mach numbers
• Free and smooth boundaries in 2-D finite-difference schemes for transient elastic waves
• High order ADER schemes for a unified first order hyperbolic formulation of Newtonian continuum mechanics coupled with electro-dynamics
• Space-time adaptive ADER-DG schemes for dissipative flows: compressible Navier-Stokes and resistive MHD equations
• Semi-implicit discontinuous Galerkin methods for the incompressible Navier-Stokes equations on adaptive staggered Cartesian grids
• Arbitrary high order accurate space-time discontinuous Galerkin finite element schemes on staggered unstructured meshes for linear elasticity

Cited by in corpus (11)

• High order ADER schemes for continuum mechanics
• ExaHyPE: An Engine for Parallel Dynamically Adaptive Simulations of Wave Problems
• A structure-preserving staggered semi-implicit finite volume scheme for continuum mechanics
• A unified first order hyperbolic model for nonlinear dynamic rupture processes in diffuse fracture zones
• Space-time adaptive ADER discontinuous Galerkin schemes for nonlinear hyperelasticity with material failure
• Simulation of non-Newtonian viscoplastic flows with a unified first order hyperbolic model and a structure-preserving semi-implicit scheme
• Modelling wave propagation in elastic solids via high-order accurate implicit-mesh discontinuous Galerkin methods
• Enclave Tasking for Discontinuous Galerkin Methods on Dynamically Adaptive Meshes
• TeaMPI -- Replication-based Resilience without the (Performance) Pain
• A cell-centered implicit-explicit Lagrangian scheme for a unified model of nonlinear continuum mechanics on unstructured meshes
• Efficient implementation of ADER discontinuous Galerkin schemes for a scalable hyperbolic PDE engine