A discontinuous Galerkin method for three-dimensional poroelastic wave propagation: forward and adjoint problems

DOI10.1007/S40315-021-00395-4MaRDI QIDQ2129488zbMATHOpenAlexWikidataFDO

Authors Nick Dudley Ward, Timo Lähivaara, Simon P. Eveson

Publication date 22 April 2022

Published in Computational Methods and Function Theory (Search for Journal in Brave)

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

zbMATH Keywords

discontinuous Galerkin methodadjoint methodporoelastic waves

Mathematics Subject Classification ID

Inverse problems for PDEs (35R30) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Finite element methods applied to problems in fluid mechanics (76M10) Computational methods for problems pertaining to geophysics (86-08)

Abstract: We develop a numerical solver for three-dimensional wave propagation in coupled poroelastic-elastic media, based on a high-order discontinuous Galerkin (DG) method, with the Biot poroelastic wave equation formulated as a first order conservative velocity/strain hyperbolic system. To derive an upwind numerical flux, we find an exact solution to the Riemann problem, including the poroelastic-elastic interface; we also consider attenuation mechanisms both in Biot's low- and high-frequency regimes. Using either a low-storage explicit or implicit-explicit (IMEX) Runge-Kutta scheme, according to the stiffness of the problem, we study the convergence properties of the proposed DG scheme and verify its numerical accuracy. In the Biot low frequency case, the wave can be highly dissipative for small permeabilities; here, numerical errors associated with the dissipation terms appear to dominate those arising from discretisation of the main hyperbolic system. We then implement the adjoint method for this formulation of Biot's equation. In contrast with the usual second order formulation of the Biot equation, we are not dealing with a self-adjoint system but, with an appropriate inner product, the adjoint may be identified with a non-conservative velocity/stress formulation of the Biot equation. We derive dual fluxes for the adjoint and present a simple but illuminating example of the application of the adjoint method.

Recommendations

• A discontinuous Galerkin method for poroelastic wave propagation: the two-dimensional case
• A discontinuous Galerkin finite element method for dynamic and wave propagation problems in non-linear solids and saturated porous media
• A nodal discontinuous Galerkin finite element method for the poroelastic wave equation
• A discontinuous Galerkin method for wave propagation in orthotropic poroelastic media with memory terms
• Discontinuous Galerkin method for wave propagation in elastic media with inhomogeneous inclusions
• A weight-adjusted discontinuous Galerkin method for the poroelastic wave equation: penalty fluxes and micro-heterogeneities
• Numerical modeling of wave propagation phenomena in thermo-poroelastic media via discontinuous Galerkin methods
• Multiphysics discontinuous Galerkin method for a poroelasticity model
• Analysis of a high-order space and time discontinuous Galerkin method for elastodynamic equations. Application to 3D wave propagation
• Full-anisotropic poroelastic wave modeling: a discontinuous Galerkin algorithm with a generalized wave impedance

Cites work

• scientific article; zbMATH DE number 5179816 (Why is no real title available?)
• A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs
• A discontinuous Galerkin method for poroelastic wave propagation: the two-dimensional case
• A nodal discontinuous Galerkin finite element method for the poroelastic wave equation
• A weight-adjusted discontinuous Galerkin method for the poroelastic wave equation: penalty fluxes and micro-heterogeneities
• Additive Runge-Kutta schemes for convection-diffusion-reaction equations
• Discretely exact derivatives for hyperbolic PDE-constrained optimization problems discretized by the discontinuous Galerkin method
• Estimation of aquifer dimensions from passive seismic signals with approximate wave propagation models
• Finite Volume Methods for Hyperbolic Problems
• High-resolution finite volume modeling of wave propagation in orthotropic poroelastic media
• Low-storage implicit/explicit Runge-Kutta schemes for the simulation of stiff high-dimensional ODE systems
• Mechanics of Deformation and Acoustic Propagation in Porous Media
• Spectral-element and adjoint methods in seismology
• The adjoint method in seismology
• Unified Riemann solution for multi-physics coupling: anisotropic poroelastic/elastic/fluid interfaces

Cited in

(6)

• A discontinuous Galerkin method for poroelastic wave propagation: the two-dimensional case
• Implementation of hybridizable discontinuous Galerkin method for time-harmonic anisotropic poroelasticity in two dimensions
• Full-anisotropic poroelastic wave modeling: a discontinuous Galerkin algorithm with a generalized wave impedance
• A nodal discontinuous Galerkin finite element method for the poroelastic wave equation
• A weight-adjusted discontinuous Galerkin method for the poroelastic wave equation: penalty fluxes and micro-heterogeneities
• A discontinuous Galerkin method for wave propagation in orthotropic poroelastic media with memory terms

This page was built for publication: A discontinuous Galerkin method for three-dimensional poroelastic wave propagation: forward and adjoint problems

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