Efficient PML for the wave equation
arXiv:1001.0319
Abstract
In the last decade, the perfectly matched layer (PML) approach has proved a flexible and accurate method for the simulation of waves in unbounded media. Most PML formulations, however, usually require wave equations stated in their standard second-order form to be reformulated as first-order systems, thereby introducing many additional unknowns. To circumvent this cumbersome and somewhat expensive step, we instead propose a simple PML formulation directly for the wave equation in its second-order form. Inside the absorbing layer, our formulation requires only two auxiliary variables in two space dimensions and four auxiliary variables in three space dimensions; hence it is cheap to implement. Since our formulation requires no higher derivatives, it is also easily coupled with standard finite difference or finite element methods. Strong stability is proved while numerical examples in two and three space dimensions illustrate the accuracy and long time stability of our PML formulation.
16 pages, 6 figures
Cited by in corpus (6)
- A General Approach to Seismic Inversion with Automatic Differentiation
- Review and Recent Developments on the Perfectly Matched Layer (PML) Method for the Numerical Modeling and Simulation of Elastic Wave Propagation in Unbounded Domains
- Distributed Machine Learning for Computational Engineering using MPI
- Learning Hidden Dynamics using Intelligent Automatic Differentiation
- Fractional Buffer Layers: Absorbing Boundary Conditions for Wave Propagation
- Full waveform inversion using triangular waveform adapted meshes