Abstract

The Multi-Transmitting Formula (MTF) proposed by Liao et al. is a local artificial boundary condition widely used in numerical simulations of wave propagation in an infinite medium, while the drift instability is usually caused in its numerical implementation. In view of a physical interpretation of the Gustafsson, Kreiss and Sundström criterion on numerical solutions of initial-boundary value problems in the hyperbolic partial differential equations, the mechanism of the drift instability of MTF was discussed, and a simple measure for eliminating the drift instability was proposed by introducing a modified operator into the MTF. Based on the theory of spherical wave propagation and damping effect of medium, the physical implication on modified operator was interpreted. And the effect of the modified operator on the reflection coefficient of MTF was discussed. Finally, the validity of the proposed stable implementation measure was verified by numerical tests of wave source problem and scattering problem.

Highlights

  • For the numerical simulations of near-field wave motions and the response of geological structures, the control equations of different media should be determined to obtain the reliable wave propagation characteristics [1,2,3,4,5]

  • We proposed to introduce a small positive parameter into the boundary conditions in the numerical implementation of Multi-Transmitting Formula (MTF), and Eq (4) should be rewritten as: 1⁄2ð1 þ gÞB00 À B11 N up0þ1 1⁄4 0

  • The Mechanism of drift instability in the numerical implementation of MTF is theoretically analyzed, and it reveals that drift instability is caused by the reason that the GKS criterion is violated in the MTF with a zero frequency and zero wavenumber

Read more

Summary

OPEN ACCESS

Citation: Su J, Zhou Z, Li Y, Hao B, Dong Q, Li X (2020) A stable implementation measure of multitransmitting formula in the numerical simulation of wave motion. PLoS ONE 15(12): e0243979. https://doi.org/10.1371/journal.pone.0243979 Data Availability Statement: All relevant data are within the paper and its Supporting Information files. Funding: National Natural Science Foundation of China (U1839202, 41374049) and the National Key Research and Development Program of China (2017YFC1500400). Competing interests: The authors have declared that no competing interests exist.

Introduction
Mechanism of and eliminating measure for MTF drift instability
Case A
Case B
Numerical experiments
Wave source problem
Scattering problem
Conclusions
Author Contributions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call