IMPLEMENTATION OF A NON-SPLITTING FORMULATION OF PERFECTLY MATCHED LAYER IN A 3D – 4TH ORDER STAGGERED-GRID VELOCITY-STRESS FINITE-DIFFERENCE SCHEME


Published: Jan 25, 2010
Keywords:
Finite-Difference Perfectly Matched Layer Non-Splitting Formulation
A.A Skarlatoudis
C.B. Papazachos
Abstract
One of the major problems in numerical simulations of wave propagation in elastodynamics using grid-point methods is the truncation of the computational space by artificial boundaries. These boundaries produce spurious reflections, polluting the results with artificial noise, therefore several efforts have been realized in order to achieve transparent or non-reflecting boundaries for truncating the computational space. Perfectly Matched Layer (PML) has been one of the most efficient methods for implementing artificial boundaries at the edges of the computational models. However the application of PML requires the “tuning” of several variables, for which very limited work has been presented. In the present study we employ the Non-Splitting formulation of PML (NPML) technique presented by Wang & Tang (2003), based on the introduction of small perturbations in the wavefield. The NPML formulation shows the same accuracy with the originally introduced PML but is easier to be implemented because it does not require the field splitting. The main scope of this paper is to perform a full analysis of the efficiency of NMPL based on numerical tests using a 3D–4th order staggered-grid velocity-stress finite-difference scheme. The analysis consists of direct, quantified comparisons with reference solution waveforms from a semi-analytical method (Discrete Wavenumber Method) with synthetic waveforms produced using the most popular Absorbing Boundary Conditions (ABCs) and PML, in various canonical models (homogeneous and layer over half-space).
Article Details
  • Section
  • Geophysics
Downloads
Download data is not yet available.
References
Berenger, J., (1994). A perfectly matched layer for the absorption of electromagnetic waves: J. Comput.
Phys., 114, 185–200.
Berenger, J. P. (1996). Three-dimensional perfectly matched layer for the absorption of electromagnetic
waves, J. Comp. Phys., 127, 363–379.
XLIII, No 4 – 1939
Bouchon, M. (1981). A simple method to calculate Green’s functions for elastic layered media, Bull.
Seism. Soc. Am., 71, 959–971.
Chew, W., and Q. H. Liu, (1996). Perfectly matched layers for elastodynamics: a new absorbing boundary
condition, J. Comp. Acoust., 4, 341–359.
Clayton, R., and B. Engquist, (1977). Absorbing boundary conditions for acoustic and elastic wave equations,
Bull. Seismol. Soc. Am., 67, 1529 – 1540.
Collino, F., and P., Monk, (1998a). Optimizing the perfectly matched layer, Comput. Methods Appl. Mech.
Eng., 164, 157–171.
Collino, F., and C. Tsogka, (2001). Application of the PML absorbing layer model to the linear elastodynamic
problem in anisotropic heterogeneous media, Geophysics, 66, 294–307.
Coutant, O. (1989). Program of numerical simulation AXITRA. Res. Rep. LGIT (in French), Universite
Joseph Fourier, Grenoble.
Engquist, B., and Majda, A., (1997). Absorbing boundary conditions for the numerical simulation of
Waves: Math. Comp., 31, 629–651.
Festa, G. and S. Nielsen, (2003). PML Absorbing Boundaries, Bull. Seismol. Soc. Am., 93, 891–903.
Frankel, A. (1993). Three-dimensional simulations of ground motions in the San Bernardino Valley, California,
for hypothetical earthquakes on the San Andreas fault, Bull. Seism. Soc. Am., 83, 1020–1041.
Graves, R. W. (1996). Simulating seismic wave propagation in 3D elastic media using staggered-grid finite
differences, Bull. Seism. Soc. Am., 86, 1091–1106.
Hastings, F., Schneider, J. B., and Broschat, S. L., (1996). Application of the perfectly matched layer (PML)
absorbing boundary condition to elastic wave propagation, J. Acoust. Soc. Am., 100, 3061– 3069.
Higdon, R. L., (1991). Absorbing boundary conditions for elastic waves, Geophysics, 56, 231-241.
Komatitsch D. and R. Martin (2007). An unsplit convolutional perfectly matched layer improved at grazing
incidence for the seismic wave equation, Geophysics, 72, 155–167.
Kristek, J., P. Moczo, and R. J. Archuleta, (2002). Efficient methods to simulate planar free surface in the
D 4th-order staggered-grid finite-difference schemes, Studia Geophys. Geodet., 46, 355–381.
Kristek, J., and P. Moczo, (2003). Seismic wave propagation in viscoelastic media with material discontinuities—A
D 4th-order staggered-grid finite-difference modeling, Bull. Seismol. Soc. Am., 93,
–2280.
Kristekova, M., J. Kristek, P. Moczo and S.M. Day, (2006). Misfit criteria for quantitative comparison of
seismograms, Bull. Seismol. Soc. Am., 96, 1836-1850.
Marcinkovich C. and K. Olsen, (2003). On the implementation of perfectly matched layers in a three-dimensional
fourth-order velocity-stress finite difference scheme, J. Geophys. Res., 108, 2276 – 2292.
Martin R., D. Komatitsch, and A. Ezziani, (2008). An unsplit convolutional perfectly matched layer improved
at grazing incidence for seismic wave propagation in poroelastic media, Geophysics, 73, 51–61.
Moczo, P., J. Kristek, and L. Halada (2000). 3D fourth-order staggered grid finite-difference schemes: stability
and grid dispersion, Bull. Seism. Soc. Am., 90, 587–603.
Moczo, P., J. Kristek, V. Vavrycuk, R. J. Archuleta, and L. Halada, (2002). 3D heterogeneous staggeredgrid
finite-difference modeling of seismic motion with volume harmonic and arithmetic averaging of
elastic moduli and densities, Bull. Seism. Soc. Am., 92, 3042–3066.
Peng, C. and M. N. Toksöz, (1994). An optimal absorbing boundary condition for finite difference modeling
of acoustic and elastic wave propagation, J. Acoust. Soc. Am., 95, 733-745.
Peng, C., and Toksoz, M., (1995). An optimal absorbing boundary condition for elastic wave modeling,
Geophysics, 60, 296–301.
Wang, T., and X. Tang, (2003). Finite-difference modeling of elastic wave propagation: A nonsplitting perfectly
matched layer approach, Geophysics, 68, 1749–1755.
Most read articles by the same author(s)