Numerical solution of convected wave equation in free field using artificial boundary method

Abstract

In this article, we propose two procedures focusing on the computation of the time‐dependent convected wave equation in a free field with a uniform background flow. Both procedures are based on a framework, expended from Du et al. (SIAM J. Sci. Comput. 40 (2018), A1430–A1445.), of constructing the Dirichlet‐to‐Dirichlet (DtD)‐type discrete absorbing boundary conditions (ABCs). The first procedure is dedicated to reducing the infinite problem into a finite problem by a direct application of the framework on the finite difference discretization of the convected wave equation. However, the presence of convection terms makes the stability analysis hard to implement, which motivates us to develop the second procedure. First, the convected wave equation is transformed into a standard wave equation by using the Prandtl‐Glauert‐Lorentz transformation. After that, we obtain the DtD‐type ABC by using the above framework, and on this basis, derive an equivalent Dirichlet‐to‐Neumann‐type ABCs, which makes stability and convergence analysis easy to be obtained by the classical energy method. The effectiveness and comparison of these two procedures are investigated through numerical experiments.