volume | PIER Journals

Abstract

Modeling wave propagation often requires a truncation of the computational domain to a smaller subdomain to keep computational cost reasonable. The mere volume of papers on absorbing boundary conditions indicates that a perfect solution is not available. A method is proposed that is exact, at least in the case of a time-domain finite-difference scheme for the scalar wave equation. The word `exact' is used in the sense that there is no difference between a computation on the truncated domain with this method and one on an enlarged domain with reflecting boundaries that are placed so far away that their reflections cannot reach the original domain within the modeled time span. Numerical tests in 1D produce stable results with central difference schemes from order 2 to 24 for the spatial discretization. The difference with a reference solution computed on an enlarged domain with the boundary moved sufficiently far away only contains accumulated numerical round-off errors. Generalization to more than one space dimension is feasible if there is a single non-reflecting boundary on one side of a rectangular domain or two non-reflecting boundaries at opposing sides, but not for a corner connecting non-reflecting boundaries. The reason is that the method involves recursion based on translation invariance in the direction perpendicular to the boundary, which does not hold in the last case. This limits the applicability of the method to, for instance, modeling waveguides.

1. Givoli, D., "Non-reflecting boundary conditions," Journal of Computational Physics, Vol. 94, No. 1, 1-29, 1991.
doi:10.1016/0021-9991(91)90135-8 Google Scholar

2. Mulder, W. A., "Experiments with Higdon’s absorbing boundary conditions for a number of wave equations," Computational Geosciences, Vol. 1, No. 1, 85-108, 1997.
doi:10.1023/A:1011556926362 Google Scholar

3. Tsynkov, S. V., "Numerical solution of problems on unbounded domains. A review," Applied Numerical Mathematics, Vol. 27, No. 4, 465-532, Special Issue on Absorbing Boundary Conditions, 1998.
doi:10.1016/S0168-9274(98)00025-7 Google Scholar

4. Tourrette, L. and L. Halpern, Absorbing Boundaries and Layers, Domain Decomposition Methods: Applications to Large Scale Computers, Nova Science Publishers, Inc., 2001.

5. Givoli, D., "High-order local non-reflecting boundary conditions: A review," Wave Motion, Vol. 39, No. 4, 319-326, 2004.
doi:10.1016/j.wavemoti.2003.12.004 Google Scholar

6. Antoine, X., E. Lorin, and Q. Tang, "A friendly review of absorbing boundary conditions and perfectly matched layers for classical and relativistic quantum waves equations," Molecular Physics, Vol. 115, No. 15–16, 1861-1879, 2017.
doi:10.1080/00268976.2017.1290834 Google Scholar

7. Gao, Y., H. Song, J. Zhang, and Z. Yao, "Comparison of artificial absorbing boundaries for acoustic wave equation modelling," Exploration Geophysics, Vol. 48, No. 1, 76-93, 2017.
doi:10.1071/EG15068 Google Scholar

8. Berenger, J.-P., "A perfectly matched layer for the absorption of electromagnetic waves," Journal of Computational Physics, Vol. 114, No. 2, 185-200, 1994.
doi:10.1006/jcph.1994.1159 Google Scholar

9. Mulder, W. A. and B. Hak, "An ambiguity in attenuation scattering imaging," Geophysical Journal International, Vol. 178, No. 3, 1614-1624, 2009.
doi:10.1111/j.1365-246X.2009.04253.x Google Scholar

10. Ting, L. and M. J. Miksis, "Exact boundary conditions for scattering problems," The Journal of the Acoustical Society of America, Vol. 80, No. 6, 1825-1827, 1986.
doi:10.1121/1.394297 Google Scholar

11. Givoli, D. and D. Cohen, "Nonreflecting boundary conditions based on Kirchhoff-type formulae," Journal of Computational Physics, Vol. 117, No. 1, 102-113, 1995.
doi:10.1006/jcph.1995.1048 Google Scholar

12. Teng, Z.-H., "Exact boundary condition for time-dependent wave equation based on boundary integral," Journal of Computational Physics, Vol. 190, No. 2, 398-418, 2003.
doi:10.1016/S0021-9991(03)00281-X Google Scholar

13. Engquist, B. and A. Majda, "Radiation boundary conditions for acoustic and elastic wave calculations," Communications on Pure and Applied Mathematics, Vol. 32, No. 3, 313-357, 1979.
doi:10.1002/cpa.3160320303 Google Scholar

14. Higdon, R. L., "Absorbing boundary conditions for difference approximations to the multi-dimensional wave equation," Mathematics of Computation, Vol. 47, No. 176, 437-459, 1986. Google Scholar

15. Mur, G., "Total-field absorbing boundary conditions for the time-domain electromagnetic field equations," IEEE Transactions on Electromagnetic Compatibility, Vol. 40, No. 2, 100-102, 1998.
doi:10.1109/15.673614 Google Scholar

16. Komatitsch, D. and R. Martin, "An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation," Geophysics, Vol. 72, No. 5, SM155-SM167, 2007.
doi:10.1190/1.2757586 Google Scholar

17. Zhebel, E., S. Minisini, A. Kononov, and W. A. Mulder, "A comparison of continuous mass-lumped finite elements with finite differences for 3-D wave propagation," Geophysical Prospecting, Vol. 62, No. 5, 1111-1125, 2014.
doi:10.1111/1365-2478.12138 Google Scholar

18. Fornberg, B., "Generation of finite difference formulas on arbitrarily spaced grids," Mathematics of Computation, Vol. 51, No. 184, 699-706, 1988.
doi:10.1090/S0025-5718-1988-0935077-0 Google Scholar

19. Higdon, R. L., "Numerical absorbing boundary conditions for the wave equation," Mathematics of Computation, Vol. 49, No. 179, 65-90, 1987.
doi:10.1090/S0025-5718-1987-0890254-1 Google Scholar

20. Mulder, W. A., "Working around the corner problem in numerically exact non-reflecting boundary conditions for the wave equation," Conference Proceedings, 82nd EAGE Conference and Exhibition 2020, Amsterdam, The Netherlands, 2020. Google Scholar

Mr. William Alexander Mulder