Energy conservative SBP discretizations of the acoustic wave equation in covariant form on staggered curvilinear grids

DOI10.1016/J.JCP.2020.109386zbMATH Open1436.65108arXiv1907.01105OpenAlexW2954954758MaRDI QIDQ777587FDOQ777587

Authors: Ossian O'Reilly, N. Anders Petersson 

Publication date: 7 July 2020

Published in: Journal of Computational Physics (Search for Journal in Brave)

Abstract: We develop a numerical method for solving the acoustic wave equation in covariant form on staggered curvilinear grids in an energy conserving manner. The use of a covariant basis decomposition leads to a rotationally invariant scheme that outperforms a Cartesian basis decomposition on rotated grids. The discretization is based on high order Summation-By-Parts (SBP) operators and preserves both symmetry and positive definiteness of the contravariant metric tensor. To improve accuracy and decrease computational cost, we also derive a modified discretization of the metric tensor that leads to a conditionally stable discretization. Bounds are derived that yield a point-wise condition that can be evaluated to check for stability of the modified discretization. This condition shows that the interpolation operators should be constructed such that their norm is close to one.

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

Recommendations

• SBP-SAT finite difference discretization of acoustic wave equations on staggered block-wise uniform grids
• Stormer-Numerov HDG methods for acoustic waves
• Symplectic-mixed finite element approximation of linear acoustic wave equations
• Acoustic wave propagation in complicated geometries and heterogeneous media
• A new mimetic scheme for the acoustic wave equation

zbMATH Keywords

curvilinear coordinateshigh order accuracystaggered gridsacoustic wave equationsummation-by-partscovariant formulation

Mathematics Subject Classification ID

Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06)

Cites Work

• The Theoretical Accuracy of Runge–Kutta Time Discretizations for the Initial Boundary Value Problem: A Study of the Boundary Error
• Title not available (Why is that?)
• Introduction to tensor analysis and the calculus of moving surfaces
• Summation by parts for finite difference approximations for \(d/dx\)
• Finite volume methods, unstructured meshes and strict stability for hyperbolic problems
• A skew-symmetric discontinuous Galerkin spectral element discretization and its relation to SBP-SAT finite difference methods
• An explicit fourth-order orthogonal curvilinear staggered-grid FDTD method for Maxwell's equations
• Convergence of summation-by-parts finite difference methods for the wave equation
• Stable coupling of nonconforming, high-order finite difference methods
• Stable and Accurate Interpolation Operators for High-Order Multiblock Finite Difference Methods
• Stable Grid Refinement and Singular Source Discretization for Seismic Wave Simulations
• A fourth order accurate finite difference scheme for the elastic wave equation in second order formulation
• A stable hybrid method for hyperbolic problems
• Title not available (Why is that?)
• The Convergence Rate for Difference Approximations to Mixed Initial Boundary Value Problems
• Title not available (Why is that?)
• Summation by Parts, Projections, and Stability. I
• Summation by Parts, Projections, and Stability. II
• Wave propagation in anisotropic elastic materials and curvilinear coordinates using a summation-by-parts finite difference method
• Finite-difference strategy for elastic wave modelling on curved staggered grids
• Discretizing singular point sources in hyperbolic wave propagation problems
• Diagonal-norm summation by parts operators for finite difference approximations of third and fourth derivatives
• A hybrid framework for coupling arbitrary summation-by-parts schemes on general meshes
• Treatment of the polar coordinate singularity in axisymmetric wave propagation using high-order summation-by-parts operators on a staggered grid
• Diagonal-norm upwind SBP operators
• Energy stable and high-order-accurate finite difference methods on staggered grids
• High order accurate finite difference modeling of seismo-acoustic wave propagation in a moving atmosphere and a heterogeneous Earth model coupled across a realistic topography
• SBP-SAT finite difference discretization of acoustic wave equations on staggered block-wise uniform grids
• Combining finite element and finite difference methods for isotropic elastic wave simulations in an energy-conserving manner
• Summation-by-parts operators for correction procedure via reconstruction

Cited In (10)

• Summation-by-parts finite-difference shallow water model on the cubed-sphere grid. I: Non-staggered grid
• Full-spectrum dispersion relation preserving summation-by-parts operators
• Elastic Wave Propagation in Curvilinear Coordinates with Mesh Refinement Interfaces by a Fourth Order Finite Difference Method
• Provably non-stiff implementation of weak coupling conditions for hyperbolic problems
• High-order mimetic finite differences for anisotropic elliptic equations
• Dual-pairing summation by parts finite difference methods for large scale elastic wave simulations in 3D complex geometries
• Mimetic finite difference operators and higher order quadratures
• Acoustic shape optimization using energy stable curvilinear finite differences
• Implicit summation by parts operators for finite difference approximations of first and second derivatives
• Boundary-optimized summation-by-parts operators for finite difference approximations of second derivatives with variable coefficients

This page was built for publication: Energy conservative SBP discretizations of the acoustic wave equation in covariant form on staggered curvilinear grids

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