Hostname: page-component-78c5997874-4rdpn Total loading time: 0 Render date: 2024-11-16T01:20:29.877Z Has data issue: false hasContentIssue false

A Well-Posed and Discretely Stable Perfectly Matched Layer for Elastic Wave Equations in Second Order Formulation

Published online by Cambridge University Press:  20 August 2015

Kenneth Duru*
Affiliation:
Division of Scientific Computing, Department of Information Technology, Uppsala University, Box 337, SE-751 05 Uppsala, Sweden
Gunilla Kreiss*
Affiliation:
Division of Scientific Computing, Department of Information Technology, Uppsala University, Box 337, SE-751 05 Uppsala, Sweden
*
Corresponding author.Email:[email protected]
Get access

Abstract

We present a well-posed and discretely stable perfectly matched layer for the anisotropic (and isotropic) elastic wave equations without first re-writing the governing equations as a first order system. The new model is derived by the complex coordinate stretching technique. Using standard perturbation methods we show that complex frequency shift together with a chosen real scaling factor ensures the decay of eigen-modes for all relevant frequencies. To buttress the stability properties and the robustness of the proposed model, numerical experiments are presented for anisotropic elastic wave equations. The model is approximated with a stable node-centered finite difference scheme that is second order accurate both in time and space.

Type
Research Article
Copyright
Copyright © Global Science Press Limited 2012

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Sjogreen, B. and Petersson, N. A., Perfectly matched layer for Maxwell’s equations in second order formulation, J. Comput. Phys., 209 (2005), 1946.Google Scholar
[2]Kreiss, H.-O., Petersson, N. A. and Ystrom, J., Difference approximations for second order wave equation, SIAM J. Numer. Anal., 40 (2002), 1940459.CrossRefGoogle Scholar
[3]Berenger, J. P., A perfectly matched layer for the absorption of electromagnetic waves, J. Comput. Phys., 114 (1994), 185200.Google Scholar
[4]Becache, E., Fauqueux, S. and Joly, P., Stability of perfectly matched layers, group velocities and anisotropic waves, J. Comput. Phys., 188 (2003), 399433.CrossRefGoogle Scholar
[5]Appelo, D. and Kreiss, G., A new absorbing layer for elastic waves, J. Comput. Phys., 215 (2006), 642660.CrossRefGoogle Scholar
[6]Taflove, A., Advances in Computational Electrodynamics, The Finite-Difference TimeDomain, Artec House Inc, 1998.Google Scholar
[7]Diaz, J. and Joly, P., Stabilized perfectly matched layer for advective acoustics, Waves, (2003), 115119.Google Scholar
[8]Skelton, E. A., Adams, S. D. M. and Craster, R. V., Guided elastic waves and perfectly matched layers, Wave. Motion, 44 (2007), 573592.Google Scholar
[9]Abarbanel, S., Gottlieb, D. and Hesthaven, J. S., Long time behaviour of the perfectly matched layer equations in computational electromagnetics, J. Sci. Comput., 17(1-4) (2002), 405422.CrossRefGoogle Scholar
[10]Appelo, D. and Petersson, A. N., A stable finite difference method for the elastic wave equation on complex domains with free surfaces, Commun. Comput. Phys., 5 (2009), 84107.Google Scholar
[11]Abarbanel, S. and Gottlieb, D., A mathematicalanalysis of the PML method, J.Comput. Phys., 134 (1997), 357363.Google Scholar
[12]Gustafsson, B., Kreiss, H.-O. and Oliger, J., Time Dependent Problems and Difference Methods, John Wileys and Sons, 1995.Google Scholar
[13]Nilsson, S., Petersson, N. A., Sjogreen, B. and Kreiss, H.-O., Stable difference approximations for the elastic wave equation in second order formulation, SIAM J. Numer. Anal., 42 (2004), 12921323.Google Scholar
[14]Kuzuoglu, M. and Mittra, R., Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers, IEEE Microw. Guided. W., 6(12) (1996), 447449.Google Scholar
[15]Roden, J. A. and Gedney, D. S., Convolutional PML (CPML): an efficient fdtd implementation of the CFS-PML for arbitrary media, Microw. Opt. Tech. Lett., 27(5) (2000), 334339.Google Scholar
[16]Meza-Fajardo, K. C. and Papageogiou, A. S., A nonconvolutional, split-field, perfectly matched layer for wave propagation in isotropic and anisotropic elastic media: stability analysis, B. Seismol. Soc. Am., 98(4) (2008), 18111836.Google Scholar
[17]Chew, W. and Weedon, W., A 3-D perfectly matched medium from modified Maxwell’s equations with stretched coordinates, Micro. Opt. Tech. Lett., 7(13) (1994), 599604.Google Scholar
[18]Komatitsch, D. and Tromp, J., A perfectly matched layer absorbing boundary condition for the second order seismic wave equation, Geophys. J. Int., 154 (2003), 146153.CrossRefGoogle Scholar
[19]Becache, E., Petropoulos, P. G. and Gedney, S. D., On the long-time behaviour of unplit perfectly matched layers, IEEE T. Antenn. Propag., 52(5) (2004), 13351342.Google Scholar
[20]Gallina, P., Effect of damping on assymetric systems, J. Vibrat. Acoust., 125 (2003), 359365.Google Scholar
[21]Appelo, D. and Colonius, T., A high-order super-grid absorbing layer and its application to linear hyperbolic systems, J. Comput. Phys., 228 (2009), 42004217.Google Scholar
[22]Thomsen, L., Weak elastic anisotropy, Geophysics, 51(10) (1986), 19541966.Google Scholar
[23]Collino, F. and Tsogka, C., Application of the PML absorbing layer model to the linear elasto-dynamic problem in anisotropic heterogeneous media, Geophysics, 66 (2001), 294307.Google Scholar
[24]Kato, T., Perturbation Theory for Linear Operators, Springer Verlag, 1966.Google Scholar
[25]Lai, W. M., Rubin, D. and Krempl, E., Introduction to Continuum Mechanics, Butterworth-Heimann Ltd., US, 1993.Google Scholar
[26]Duru, K., Perfectly Matched Layer For Second Order Wave Equations, Licentiate Thesis, Div Sc. Comp., Dept. of Infor Tech., Uppsala University, ISSN 1404-5117; 2010-004, 2010.Google Scholar