Commun. Comput. Phys., 11 (2012), pp. 1643-1672.

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

Kenneth Duru 1*, Gunilla Kreiss 1

1 Division of Scientific Computing, Department of Information Technology, Uppsala University, Box 337, SE-751 05 Uppsala, Sweden.

Received 12 February 2010; Accepted (in revised version) 24 May 2011
Available online 12 January 2012


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.

AMS subject classifications: 35B35, 35L05, 35L15, 37C75

Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: Perfectly matched layer, well-posedness, stability, hyperbolicity, elastic waves.

*Corresponding author.
Email: (K. Duru), (G. Kreiss)

The Global Science Journal