A symmetric discretization of the perfectly matched layer for the 2D Helmholtz equation
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.
- Aminzadeh, F., Brac, J. and Kunz, T., 1997. 3D Salt and Overthrust Model. Modeling
- Series i, SEG, Tulsa, OK.
- Baysal, E., Kosloff, D.D. and Sherwood, J.W.C., 1983. Reverse time migration.
- Geophysics, 48: 1514-1524.
- Bérenger, J.-P., 1994. A perfectly matched layer for the absorption of electromagnetic
- waves: J. Computat. Phys., 114: 185-200.
- Bérenger, J.-P., 2007. Perfectly matched layer (pml) for computational
- electromagnetics. Synthesis Lectures on Computat. Electromagnet., 2: 1-117.
- Bunch, J.R., and Kaufman, L., 1977. Some stable methods for calculating inertia and
- solving symmetric linear systems. Mathemat. Computat., 31: 163-179.
- Chen, J.-B., 2013. A generalized optimal 9-point scheme for frequency-domain scalar
- wave equation. J. Appl. Geophys., 92: 1-7.
- Chen, J.-B., 2014. A 27-point scheme for a 3D frequency-domain scalar wave equation
- based on an average-derivative method. Geophys. Prosp., 62: 258-277.
- Chew, W.C. and Weedon, W.H., 1994. A 3D perfectly matched medium from modified
- Maxwell's equations with stretched coordinates. Microwave Optic. Technol.
- Lett., 7: 599-604.
- George, A., 1973,. Nested dissection of a regular finite element mesh. SIAM J. Numer.
- Analys., 10: 345-363.
- 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.
- 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.
- Kalinkin, A., Anders, A. and Anders, R., 2014. Intel math Kernel library parallel direct
- sparse olver for clusters. EAGE Workshop on High Performance Computing
- for Upstream, Chania, Crete, Greece.
- Karypis, G. and Kumar, V., 1998. A fast and high quality multilevel scheme for
- partitioning irregular graphs. SIAM J. Scientif. Comput., 20: 359-392.
- Marfurt, K., 1984. Accuracy of finite-difference and finite-element modeling of the
- scalar and elastic wave equations. Geophysics, 49: 533-549.
- 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 wave
- propagation using a massively parallel direct solver: A feasibility study.
- Geophysics, 72(5): SM195-SM211.
- Pratt, R.G., Shin, C. and Hick, G., 1998 Gauss-Newton and full Newton methods in
- frequency-space seismic waveform inversion. Geophys. J. Internat., 133: 341-
- Sacks, Z., Kingsland, D., Lee, R. and Lee, J.-F., 1995. A perfectly matched anisotropic
- absorber for use as an absorbing boundary condition. [EEE Transact. Antenn.
- Propagat., 43: 1460-1463.
- Shin, C., 1988. Nonlinear Elastic Wave Inversion by Blocky Parameterization. Ph.D.
- thesis, University of Tulsa, OK.
- Tarantola, A., 1984. Inversion of seismic reflection data in the acoustic approximation.
- Geophysics, 49: 1259-1266.
- Turkel, E. and Yefet, A., 1998. Absorbing PML boundary layers for wave-like
- equations. Appl. Numer. Mathemat., 27:533-557.
