We construct global artificial boundary conditions (ABCs) for ttle numericM simulation of wave processes on unbounded domains using a special non-deteriorating algorithm that has been developed previously for the long-term computation of wave-radiation solutions. The ABCs are obtained directly for the discrete formulation of the problem; in so doing, neither a rational approximation of "non-reflecting kernels," nor discretization of the continuous boundary conditions is required. The extent of temporal nonlocality of the new ABCs appears fixed and limited; in addition, the ABCs can handle artificial boundaries of irregular shape on regular grids with no fitting/adaptation needed and no accuracy loss induced. The non-deteriorating algorithm, which is the core of the new ABCs, is inherently three-dimensional, it guarantees temporally uniform grid convergence of the solution driven by a continuously operating source on arbitrarily long time intervals, and provides unimprovable linear computational complexity with respect to the grid dimension. The algorithm is based on the presence of lacunae, i.e., aft fronts of the waves, in wave-type solutions in odd-dimension spaces. It can, ill fact, be built as a modification on top of any consistent and stable finite-difference scheme, making its grid convergence uniform in time and at the same time keeping the rate of convergence the same as that of the non-modified scheme. In the paper, we delineate the construction of the global lacunae-based ABCs in the framework of a discretized wave equation. The ABCs are obtained for the most general formulation of the problem that involves radiation of waves by moving sources (e.g., radiation of acoustic waves by a maneuvering aircraft). W_ also present systematic numerical results that corroborate the theoretical design properties of the ABCs' algorithm. Key words, artificial boundary conditions, wave propagation, lacunae, non-deteriorating method Subject classification. Applied and Numerical Mathematics
[1]
J. Bérenger.
Three-Dimensional Perfectly Matched Layer for the Absorption of Electromagnetic Waves
,
1996
.
[2]
S.,et al.
Numerical Solution of Initial Boundary Value Problems Involving Maxwell’s Equations in Isotropic Media
,
1966
.
[3]
Jean-Pierre Berenger,et al.
A perfectly matched layer for the absorption of electromagnetic waves
,
1994
.
[4]
Leslie Greengard,et al.
Rapid Evaluation of Nonreflecting Boundary Kernels for Time-Domain Wave Propagation
,
2000,
SIAM J. Numer. Anal..
[5]
Long-time numerical computation of wave-type solutions driven by moving sources ✩
,
2001
.
[6]
T. Hagstrom.
Radiation boundary conditions for the numerical simulation of waves
,
1999,
Acta Numerica.
[7]
David Gottlieb,et al.
On the construction and analysis of absorbing layers in CEM
,
1998
.
[8]
Andreas C. Cangellaris,et al.
GT-PML: generalized theory of perfectly matched layers and its application to the reflectionless truncation of finite-difference time-domain grids
,
1996,
1996 IEEE MTT-S International Microwave Symposium Digest.
[9]
David Gottlieb,et al.
A Mathematical Analysis of the PML Method
,
1997
.