The role of numerical boundary procedures in the stability of perfectly matched layers

Kenneth Duru

DOI: https://doi.org/10.48550/arXiv.1405.0536

2014-05-03

Abstract:In this paper we address the temporal energy growth associated with numerical approximations of the perfectly matched layer (PML) for Maxwell's equations in first order form. In the literature, several studies have shown that a numerical method which is stable in the absence of the PML can become unstable when the PML is introduced. We demonstrate in this paper that this instability can be directly related to numerical treatment of boundary conditions in the PML. First, at the continuous level, we establish the stability of the constant coefficient initial boundary value problem for the PML. To enable the construction of stable numerical boundary procedures, we derive energy estimates for the variable coefficient PML. Second, we develop a high order accurate and stable numerical approximation for the PML using summation--by--parts finite difference operators to approximate spatial derivatives and weak enforcement of boundary conditions using penalties. By constructing analogous discrete energy estimates we show discrete stability and convergence of the numerical method. Numerical experiments verify the theoretical results

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the numerical instability and long - time energy growth problems caused by the introduction of the Perfectly Matched Layer (PML) in Maxwell's equations. Specifically:

1. **Instability after PML introduction**:
- In the absence of PML, some numerical methods are stable, but when PML is introduced, these methods may become unstable. This instability is manifested as a linear or exponential growth of energy in the system over time.
- This energy growth will spread into the computational domain and undermine the accuracy of numerical simulations. For example, this phenomenon may be observed in the simulation of electromagnetic waveguides when the wave is trapped in the domain for a long time.

2. **The impact of improper handling of boundary conditions**:
- The paper points out that this instability and long - time energy growth can be directly attributed to the improper handling of numerical boundary conditions in PML.
- The author shows how to ensure the stability and convergence of numerical methods by improving the handling of boundary conditions through continuous and discrete energy estimation methods.

3. **Achievement of high - order accuracy and stability**:
- One of the goals of the paper is to develop high - order accurate and stable numerical approximation methods. In particular, the Summation - by - Parts (SBP) finite - difference operator is used to approximate spatial derivatives in a bounded domain, and the imposition of boundary conditions is weakened by penalty functions.
- By constructing similar discrete energy estimates, the discrete stability and convergence of the numerical method are proved.

In summary, this paper aims to solve the numerical instability and long - time energy growth problems caused by the introduction of PML in Maxwell's equations by improving the numerical boundary treatment method of PML, and to provide high - order accurate and stable numerical approximation methods.