Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media

We present and analyze a perfectly matched, absorbing layer model for the velocity-stress formulation of elastodynamics. The principal idea of this method consists of introducing an absorbing layer in which we decompose each component of the unknown into two auxiliary components: a component orthogonal to the boundary and a component parallel to it. A system of equations governing these new unknowns then is constructed. A damping term finally is introduced for the component orthogonal to the boundary. This layer model has the property of generating no reflection at the interface between the free medium and the artificial absorbing medium. In practice, both the boundary condition introduced at the outer boundary of the layer and the dispersion resulting from the numerical scheme produce a small reflection which can be controlled even with very thin layers. As we will show with several experiments, this model gives very satisfactory results; namely, the reflection coefficient, even in the case of heterogeneous, anisotropic media, is about 1% for a layer thickness of five space discretization steps.

[1]  B. Auld,et al.  Accoustic Fields And Waves In Solids Vol-2 , 1973 .

[2]  A. Majda,et al.  Absorbing boundary conditions for the numerical simulation of waves , 1977 .

[3]  Albert C. Reynolds,et al.  Boundary conditions for the numerical solution of wave propagation problems , 1978 .

[4]  S. Orszag,et al.  Approximation of radiation boundary conditions , 1981 .

[5]  J. Virieux P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method , 1986 .

[6]  James S. Sochacki,et al.  Absorbing boundary conditions and surface waves , 1987 .

[7]  Robert L. Higdon,et al.  Radiation boundary conditions for elastic wave propagation , 1990 .

[8]  E. Bécache Resolution par une methode d'equations integrales d'un probleme de diffraction d'ondes elastiques transitoires par une fissure , 1991 .

[9]  R. Higdon Absorbing boundary conditions for elastic waves , 1991 .

[10]  Robert L. Higdon,et al.  Absorbing boundary conditions for acoustic and elastic waves in stratified media , 1992 .

[11]  An Optimal Absorbing Boundary Condition For Elastic Wave Modeling , 1993 .

[12]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[13]  Carey M. Rappaport,et al.  Perfectly matched absorbing boundary conditions based on anisotropic lossy mapping of space , 1995 .

[14]  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, IMS 1996.

[15]  John B. Schneider,et al.  Application of the perfectly matched layer (PML) absorbing boundary condition to elastic wave propagation , 1996 .

[16]  Andreas C. Cangellaris,et al.  A general approach for the development of unsplit-field time-domain implementations of perfectly matched layers for FDTD grid truncation , 1996 .

[17]  J. Bérenger Three-Dimensional Perfectly Matched Layer for the Absorption of Electromagnetic Waves , 1996 .

[18]  Eliane Bécache,et al.  Éléments finis mixtes et condensation de masse en élastodynamique linéaire. (I) Construction , 1997 .

[19]  F. Collino Perfectly Matched Absorbing Layers for the Paraxial Equations , 1997 .

[20]  Thomas Hagstrom,et al.  On High-Order Radiation Boundary Conditions , 1997 .

[21]  Jean-Pierre Berenger,et al.  Improved PML for the FDTD solution of wave-structure interaction problems , 1997 .

[22]  Stig Hestholm,et al.  Instabilities in applying absorbing boundary conditions to high‐order seismic modeling algorithms , 1998 .

[23]  P. Monk,et al.  Optimizing the Perfectly Matched Layer , 1998 .

[24]  E. Turkel,et al.  Absorbing PML boundary layers for wave-like equations , 1998 .

[25]  Weng Cho Chew,et al.  Analytical derivation of a conformal perfectly matched absorber for electromagnetic waves , 1998 .