A Non-stiff Summation-By-Parts Finite Difference Method for the Wave Equation in Second Order Form: Characteristic Boundary Conditions and Nonlinear Interfaces

06/01/2021
by   Jeremy E. Kozdon, et al.
0

Curvilinear, multiblock summation-by-parts finite difference methods with the simultaneous approximation term method provide a stable and accurate method for solving the wave equation in second order form. That said, the standard method can become arbitrarily stiff when characteristic boundary conditions and nonlinear interface conditions are used. Here we propose a new technique that avoids this stiffness by using characteristic variables to "upwind" the boundary and interface treatment. This is done through the introduction of an additional block boundary displacement variable. Using a unified energy, which expresses both the standard as well as characteristic boundary and interface treatment, we show that the resulting scheme has semidiscrete energy stability for the anistropic wave equation. The theoretical stability results are confirmed with numerical experiments that also demonstrate the accuracy and robustness of the proposed scheme. The numerical results also show that the characteristic scheme has a time step restriction based on standard wave propagation considerations and not the boundary closure.

READ FULL TEXT
research
03/02/2021

An energy-based summation-by-parts finite difference method for the wave equation in second order form

We develop an energy-based finite difference method for the wave equatio...
research
01/08/2020

On some stable boundary closures of finite difference schemes for the transport equation

We explore in this article the possibilities and limitations of the so-c...
research
10/11/2022

An efficient neural-network and finite-difference hybrid method for elliptic interface problems with applications

A new and efficient neural-network and finite-difference hybrid method i...
research
10/24/2022

Boundary and interface methods for energy stable finite difference discretizations of the dynamic beam equation

We consider energy stable summation by parts finite difference methods (...
research
11/13/2021

An explicit predictor/multicorrector time marching with automatic adaptivity for finite-strain elastodynamics

We propose a time-adaptive predictor/multi-corrector method to solve hyp...
research
08/21/2021

High-order accurate schemes for Maxwell's equations with nonlinear active media and material interfaces

We describe a fourth-order accurate finite-difference time-domain scheme...

Please sign up or login with your details

Forgot password? Click here to reset