Upwind Summation By Parts Finite Difference Methods for Large Scale Elastic Wave Simulations In Complex Geometries

11/05/2020
by   Kenneth Duru, et al.
0

High-order accurate summation-by-parts (SBP) finite difference (FD) methods constitute efficient numerical methods for simulating large-scale hyperbolic wave propagation problems. Traditional SBP FD operators that approximate first-order spatial derivatives with central-difference stencils often have spurious unresolved wave-modes in their numerical solutions. On marginally resolved computational grids, these spurious wave-modes have the potential to destroy the accuracy of numerical solutions for a first-order hyperbolic partial differential equation, such as the elastic wave equation. To ensure the accuracy of numerical solutions of the three space dimensional (3D) elastic wave equation in complex geometries, we discretise the 3D elastic wave equation with a pair of non-central (upwind) finite-difference stencils, on boundary-conforming curvilinear meshes. Using the energy method we prove that the numerical method is numerically stable, and energy conserving. Furthermore, computational results show the robustness of the scheme. We present numerical simulations of the 3D elastic wave equation in heterogeneous media with complex non-planar free surface topography, including numerical simulations of community developed seismological benchmark problems. Our results show that the upwind SBP operators are more robust and less prone to numerical dispersion errors on marginally resolved meshes when compared to traditional SBP operators, thereby increasing efficiency.

READ FULL TEXT
research
07/16/2022

Accurate simulations of nonlinear dynamic shear ruptures on pre-existing faults in 3D elastic solids with dual-pairing SBP methods

In this paper we derive and analyse efficient and stable numerical metho...
research
12/13/2020

A non-balanced staggered-grid finite-difference scheme for the first-order elastic wave-equation modeling

We introduce an efficient and accurate staggered-grid finite-difference ...
research
07/13/2021

Spatially-optimized finite-difference schemes for numerical dispersion suppression in seismic applications

Propagation characteristics of a wave are defined by the dispersion rela...
research
10/11/2021

Provably Stable Full-Spectrum Dispersion Relation Preserving Schemes

The dispersion error is often the dominant error for computed solutions ...
research
10/04/2021

An Unconditionally Stable Conformal LOD-FDTD Method For Curved PEC Objects and Its Application to EMC Problems

The traditional finite-difference time-domain (FDTD) method is constrain...
research
10/01/2022

Stability analysis of the perfectly matched layer for the elastic wave equation in layered media

In this paper, we present the stability analysis of the perfectly matche...
research
04/22/2020

Scaling through abstractions – high-performance vectorial wave simulations for seismic inversion with Devito

[Devito] is an open-source Python project based on domain-specific langu...

Please sign up or login with your details

Forgot password? Click here to reset