ARTICLE

A symmetric discretization of the perfectly matched layer for the 2D Helmholtz equation

HYUNSEO PARK1 HYEONJUN SONG2 YOONSEO PARK3 CHANGSOO SHIN2
Show Less
1 Department of Mechanical and Aerospace Engineering, Seoul National University, South Korea,
2 Department of Energy Systems Engineering, Seoul National University, South Korea,
3 Computational Science and Technology Program, Seoul National University, South Korea,
JSE 2017, 26(6), 541–560;
Submitted: 8 April 2017 | Accepted: 18 September 2017 | Published: 1 December 2017
© 2017 by the Authors. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution -Noncommercial 4.0 International License (CC-by the license) ( https://creativecommons.org/licenses/by-nc/4.0/ )
Abstract

Park, H., Song, H., Park, Y. and Shin, C., 2017. Journal of Seismic Exploration, 26: 541- 560. A symmetric discretization of the Perfectly Matched Layer (PML) for the 2D Helmholtz equation is introduced. The PML is an efficient method to suppress spurious reflections at the boundaries of the computational domain, so that the Sommerfeld radiation condition in unbounded medium is effectively achieved. The PML can be formulated in the symmetric form, which has not been used with dispersion-minimizing finite difference methods in the exploration geophysics community. We suggest a simple symmetrization of the discretized matrix that can be used with a dispersion-minimizing method. The symmetric discretization of the PML enables us to utilize the LDLT (LDL) decomposition with the Bunch-Kaufman pivoting, which considerably reduces not only the number of arithmetic operations but also storage requirement for numerical factorization of a sparse matrix, compared to the LU decomposition. Some numerical experiments are shown to demonstrate the efficiency of the suggested scheme.

Keywords
acoustic wave
Helmholtz equation
symmetric matrix
finite difference
absorbing boundary
numerical dispersion minimization.
References
  1. Aminzadeh, F., Brac, J. and Kunz, T., 1997. 3D Salt and Overthrust Model. ModelingSeries i, SEG, Tulsa, OK.
  2. Baysal, E., Kosloff, D.D. and Sherwood, J.W.C., 1983. Reverse time migration.Geophysics, 48: 1514-1524.
  3. Bérenger, J.-P., 1994. A perfectly matched layer for the absorption of electromagneticwaves: J. Computat. Phys., 114: 185-200.
  4. Bérenger, J.-P., 2007. Perfectly matched layer (pml) for computationalelectromagnetics. Synthesis Lectures on Computat. Electromagnet., 2: 1-117.
  5. Bunch, J.R., and Kaufman, L., 1977. Some stable methods for calculating inertia andsolving symmetric linear systems. Mathemat. Computat., 31: 163-179.
  6. Chen, J.-B., 2013. A generalized optimal 9-point scheme for frequency-domain scalarwave equation. J. Appl. Geophys., 92: 1-7.
  7. Chen, J.-B., 2014. A 27-point scheme for a 3D frequency-domain scalar wave equationbased on an average-derivative method. Geophys. Prosp., 62: 258-277.
  8. Chew, W.C. and Weedon, W.H., 1994. A 3D perfectly matched medium from modified
  9. Maxwell's equations with stretched coordinates. Microwave Optic. Technol.Lett., 7: 599-604.
  10. George, A., 1973,. Nested dissection of a regular finite element mesh. SIAM J. Numer.Analys., 10: 345-363.
  11. Hustedt, B., Operto, S. and Virieux, J., 2004. Mixed-grid and staggered-grid finite-difference methods for frequency-domain acoustic wave modeling. Geophys. J.Internat., 157: 1269-1296.
  12. Jo, C., Shin, C. and Suh, J., 1996. An optimal 9-point, finite-difference, frequency- space,2-D scalar wave extrapolator. Geophysics, 61: 529-537.
  13. Kalinkin, A., Anders, A. and Anders, R., 2014. Intel math Kernel library parallel directsparse olver for clusters. EAGE Workshop on High Performance Computingfor Upstream, Chania, Crete, Greece.
  14. Karypis, G. and Kumar, V., 1998. A fast and high quality multilevel scheme forpartitioning irregular graphs. SIAM J. Scientif. Comput., 20: 359-392.
  15. Marfurt, K., 1984. Accuracy of finite-difference and finite-element modeling of thescalar and elastic wave equations. Geophysics, 49: 533-549.
  16. Operto, S., Virieux, J., Amestoy, P., L'Excellent, J.-Y., Giraud, L. and Ali, H.B.H.,3D finite-difference frequency-domain modeling of visco-acoustic wavepropagation using a massively parallel direct solver: A feasibility study.Geophysics, 72(5): SM195-SM211.
  17. Pratt, R.G., Shin, C. and Hick, G., 1998 Gauss-Newton and full Newton methods infrequency-space seismic waveform inversion. Geophys. J. Internat., 133: 341-
  18. Sacks, Z., Kingsland, D., Lee, R. and Lee, J.-F., 1995. A perfectly matched anisotropicabsorber for use as an absorbing boundary condition. [EEE Transact. Antenn.Propagat., 43: 1460-1463.
  19. Shin, C., 1988. Nonlinear Elastic Wave Inversion by Blocky Parameterization. Ph.D.thesis, University of Tulsa, OK.
  20. Tarantola, A., 1984. Inversion of seismic reflection data in the acoustic approximation.Geophysics, 49: 1259-1266.
  21. Turkel, E. and Yefet, A., 1998. Absorbing PML boundary layers for wave-likeequations. Appl. Numer. Mathemat., 27:533-557.
Share
Back to top
Journal of Seismic Exploration, Electronic ISSN: 0963-0651 Print ISSN: 0963-0651, Published by AccScience Publishing