A stable discontinuous Galerkin method for the perfectly matched layer for elastodynamics in first order form

DOI10.1007/S00211-020-01160-WMaRDI QIDQ2217869zbMATH OpenOpenAlexWikidataFDO

Authors Kenneth Duru, Leonhard Rannabauer, Alice-Agnes Gabriel, Gunilla Kreiss, Michael Bader

Publication date 11 January 2021

Published in Numerische Mathematik (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1910.06477

zbMATH Keywords

boundary and interface conditionsdiscontinuous Galerkin methodelastic wavesfirst order systemshigh order accuracyLaplace transformsperfectly matched layerstability

Mathematics Subject Classification ID

Initial-boundary value problems for systems of linear first-order PDEs (35F46) Initial value problems for systems of nonlinear first-order PDEs (35F55) PDEs in connection with mechanics of deformable solids (35Q74) Laplace transform (44A10) Finite volume methods for initial value and initial-boundary value problems involving PDEs (65M08) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12) Error bounds for initial value and initial-boundary value problems involving PDEs (65M15) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Linear elasticity with initial stresses (74B10)

Abstract: We present a stable discontinuous Galerkin (DG) method with a perfectly matched layer (PML) for three and two space dimensional linear elastodynamics, in velocity-stress formulation, subject to well-posed linear boundary conditions. First, we consider the elastodynamics equation, in a cuboidal domain, and derive an unsplit PML truncating the domain using complex coordinate stretching. Leveraging the hyperbolic structure of the underlying system, we construct continuous energy estimates, in the time domain for the elastic wave equation, and in the Laplace space for a sequence of PML model problems, with variations in one, two and three space dimensions, respectively. They correspond to PMLs normal to boundary faces, along edges and in corners. Second, we develop a DG numerical method for the linear elastodynamics equation using physically motivated numerical flux and penalty parameters, which are compatible with all well-posed, internal and external, boundary conditions. When the PML damping vanishes, by construction, our choice of penalty parameters yield an upwind scheme and a discrete energy estimate analogous to the continuous energy estimate. Third, to ensure numerical stability of the discretization when PML damping is present, it is necessary to extend the numerical DG fluxes, and the numerical inter-element and boundary procedures, to the PML auxiliary differential equations. This is crucial for deriving discrete energy estimates analogous to the continuous energy estimates. By combining the DG spatial approximation with the high order ADER time stepping scheme and the accuracy of the PML we obtain an arbitrarily accurate wave propagation solver in the time domain. Numerical experiments are presented in two and three space dimensions corroborating the theoretical results.

Review text

This article discusses a a stable discontinuous Galerkin method with a perfectly matched layer (PML) for 2D and 3D space dimensional linear elastodynamics, in velocity-stress formulation, subject to well-posed linear boundary conditions. First, the elastodynamics equation in a cuboidal domain is discusses and an unsplit PML is derived by truncating the domain using complex coordinate stretching. Second, the authors develop a discontinuous Galerkin numerical method for the linear elastodynamics equation using physically motivated numerical flux and penalty parameters, which are compatible with all well-posed, internal and external, boundary conditions. Third, to ensure numerical stability of the discretization in the presence of PML damping, the numerical DG fluxes and the numerical inter-element and boundary procedures are extended to the PML auxiliary differential equations. This step turns out to be crucial for deriving discrete energy estimates analogous to the continuous energy case. Numerical experiments included in the paper support the theoretical findings.

Recommendations

• A well-posed and discretely stable perfectly matched layer for elastic wave equations in second order formulation
• Discontinuous Galerkin method based on peridynamic theory for linear elasticity
• Stability analysis of discontinuous Galerkin approximations to the elastodynamics problem
• Discontinuous Galerkin methods with staggered hybridization for linear elastodynamics
• Higher order accurate discontinuous and continuous \(p\)-Galerkin methods for linear elastodynamics
• On energy stable discontinuous Galerkin spectral element approximations of the perfectly matched layer for the wave equation
• Convergence of a stabilized discontinuous Galerkin method for incompressible nonlinear elasticity
• A stable discontinuous Galerkin method for linear elastodynamics in 3D geometrically complex elastic solids using physics based numerical fluxes
• On perfectly matched layers for discontinuous Petrov-Galerkin methods
• Discontinuous Galerkin Method for the Solution of Elasto-Dynamic and Fluid-Structure Interaction Problems

Cites work

• A new absorbing layer for elastic waves
• A perfectly matched layer for the absorption of electromagnetic waves
• A time domain analysis of PML models in acoustics
• Absorbing boundary conditions for numerical simulation of waves
• Boundary conditions and stability of a perfectly matched layer for the elastic wave equation in first order form
• Boundary waves and stability of the perfectly matched layer for the two space dimensional elastic wave equation in second order form
• Discretizing singular point sources in hyperbolic wave propagation problems
• ExaHyPE: an engine for parallel dynamically adaptive simulations of wave problems
• Guided elastic waves and perfectly matched layers
• High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: viscous heat-conducting fluids and elastic solids
• High-order local non-reflecting boundary conditions: a review
• scientific article; zbMATH DE number 3868743 (Why is no real title available?)
• scientific article; zbMATH DE number 2042291 (Why is no real title available?)
• scientific article; zbMATH DE number 872523 (Why is no real title available?)
• scientific article; zbMATH DE number 3423635 (Why is no real title available?)
• Long-time stable high-order absorbing boundary conditions for elastodynamics
• Nodal high-order methods on unstructured grids. I: Time-domain solution of Maxwell's equations
• On energy stable discontinuous Galerkin spectral element approximations of the perfectly matched layer for the wave equation
• Perfectly matched layers for convex truncated domains with discontinuous Galerkin time domain simulations
• Perfectly matched layers for Maxwell's equations in second order formulation
• Stability of perfectly matched layers, group velocities and anisotropic waves.
• The analysis of matched layers
• The role of numerical boundary procedures in the stability of perfectly matched layers

Cited in

(23)

• An energy-based discontinuous Galerkin discretization of the elastic wave equation in second order form
• On energy stable discontinuous Galerkin spectral element approximations of the perfectly matched layer for the wave equation
• A new discontinuous Galerkin method for elastic waves with physically motivated numerical fluxes
• A stable discontinuous Galerkin method for linear elastodynamics in 3D geometrically complex elastic solids using physics based numerical fluxes
• Discontinuous Galerkin discretizations and analysis for the Cohen-Monk PML model
• On mathematical and numerical modelling of multiphysics wave propagation with polytopal discontinuous Galerkin methods: a review
• An immersed discontinuous Galerkin method for wave propagation in acoustic elastic media
• Elasto-acoustic modeling and simulation for the seismic response of structures: the case of the Tahtalı dam in the 2020 İzmir earthquake
• Boundary conditions and stability of a perfectly matched layer for the elastic wave equation in first order form
• On perfectly matched layers for discontinuous Petrov-Galerkin methods
• A two-and-a-half-dimensional displacement-based PML for elastodynamic wave propagation
• An efficient time-domain perfectly matched layers formulation for elastodynamics on spherical domains
• Numerical analysis of a method for solving 2D linear isotropic elastodynamics with traction free boundary condition using potentials and finite elements
• Compact second-order time-domain perfectly matched layer formulation for elastic wave propagation in two dimensions
• Analysis of a high-order space and time discontinuous Galerkin method for elastodynamic equations. Application to 3D wave propagation
• Analysis and FDTD simulation of a perfectly matched layer for the Drude metamaterial
• Developing and analyzing an explicit unconditionally stable finite element scheme for an equivalent Bérenger’s PML model
• Analysis and application of a time-domain finite element method for the Drude metamaterial perfectly matched layer model
• Analysis and Application of Two Novel Finite Element Methods for Solving Ziolkowski’s PML Model in the Integro-Differential Form
• Development and analysis of an equivalent Drude metamaterial perfectly matched layer model
• On the stability analysis of the perfectly matched layer for the elastic wave equation in layered media
• A perfectly matched layer for damping vertically propagating waves in the compressible Boussinesq equations
• A non-uniform basis order for the discontinuous Galerkin method of the 3D dissipative wave equation with perfectly matched layer

Describes a project that uses

Uses Software

• ExaHyPE

This page was built for publication: A stable discontinuous Galerkin method for the perfectly matched layer for elastodynamics in first order form

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2217869)