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

DOI10.1137/140976443zbMATH Open1339.65112arXiv1405.0536OpenAlexW2343355535MaRDI QIDQ2804997FDOQ2804997

Authors: Kenneth Duru 

Publication date: 9 May 2016

Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)

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

Full work available at URL: https://arxiv.org/abs/1405.0536

Recommendations

• Boundary conditions and stability of a perfectly matched layer for the elastic wave equation in first order form
• Perfectly Matched Layers for Hyperbolic Systems: General Formulation, Well‐posedness, and Stability
• On the analysis of Bérenger's Perfectly Matched Layers for Maxwell's equations
• Efficient and stable perfectly matched layer for CEM
• Studies on some perfectly matched layers for one-dimensional time-dependent systems

zbMATH Keywords

efficiencystabilitynumerical experimentsenergy methodfinite differenceMaxwell's equationswell-posednessnormal mode analysispenalty methodsinitial boundary value problemperfectly matched layershigh-order accuracyguided wavessummation-by-partsboundary waves

Mathematics Subject Classification ID

Maxwell equations (35Q61) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12)

Cites Work

• Title not available (Why is that?)
• A perfectly matched layer for the absorption of electromagnetic waves
• A mathematical analysis of the PML method
• The analysis of matched layers
• Symmetric positive linear differential equations
• Title not available (Why is that?)
• Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: Methodology and application to high-order compact schemes
• A New Method of Imposing Boundary Conditions in Pseudospectral Approximations of Hyperbolic Equations
• Title not available (Why is that?)
• High Order Difference Methods for Time Dependent PDE
• Initial boundary value problems for hyperbolic systems
• Symmetric hyperbolic linear differential equations
• Perfectly Matched Layers for Hyperbolic Systems: General Formulation, Well‐posedness, and Stability
• Summation by parts for finite difference approximations for \(d/dx\)
• A skew-symmetric discontinuous Galerkin spectral element discretization and its relation to SBP-SAT finite difference methods
• Review of summation-by-parts schemes for initial-boundary-value problems
• Stability of perfectly matched layers, group velocities and anisotropic waves.
• Perfectly matched layers for Maxwell's equations in second order formulation
• On the analysis of Bérenger's Perfectly Matched Layers for Maxwell's equations
• On the construction and analysis of absorbing layers in CEM
• Boundary procedures for summation-by-parts operators
• Well-posed absorbing layer for hyperbolic problems
• On the accuracy and stability of the perfectly matched layer in transient waveguides
• Strict stability of high-order compact implicit finite-difference schemes: The role of boundary conditions for hyperbolic PDEs. I
• Strict stability of high-order compact implicit finite-difference schemes: The role of boundary conditions for hyperbolic PDEs. II
• Efficient and stable perfectly matched layer for CEM
• Error Bounded Schemes for Time-dependent Hyperbolic Problems
• A well-posed and discretely stable perfectly matched layer for elastic wave equations in second order formulation
• On the Long-Time Behavior of Unsplit Perfectly Matched Layers
• Long-time performance of unsplit pmls with explicit second order schemes
• Discrete stability of perfectly matched layers for anisotropic wave equations in first- and second-order formulation
• Long time behavior of the perfectly matched layer equations in computational electromagnetics
• Derivation of strictly stable high order difference approximations for variable-coefficient PDE
• Lacunae based stabilization of PMLs
• Title not available (Why is that?)
• Remarks on the stability of Cartesian PMLs in corners
• Boundary waves and stability of the perfectly matched layer for the two space dimensional elastic wave equation in second order form
• On the implementation of boundary conditions for the method of lines

Cited In (14)

• On energy stable discontinuous Galerkin spectral element approximations of the perfectly matched layer for the wave equation
• Perfectly matched layers for convex truncated domains with discontinuous Galerkin time domain simulations
• On the accuracy and stability of the perfectly matched layer in transient waveguides
• Combining finite element and finite difference methods for isotropic elastic wave simulations in an energy-conserving manner
• A stable discontinuous Galerkin method for the perfectly matched layer for elastodynamics in first order form
• Title not available (Why is that?)
• Application of transparent boundary conditions to high-order finite-difference schemes for the wave equation in waveguides
• Energy decay and stability of a perfectly matched layer for the wave equation
• Boundary conditions and stability of a perfectly matched layer for the elastic wave equation in first order form
• Perfectly matched layer absorbing boundary condition for nonlinear two-fluid plasma equations
• Transparent boundary conditions for the wave equation: high-order approximation and coupling with characteristic NRBCs
• A staggered-grid finite-difference scheme optimized in the time-space domain for modeling scalar-wave propagation in geophysical problems
• SBP-SAT finite difference discretization of acoustic wave equations on staggered block-wise uniform grids
• Construction of stable PMLs for general \(2\times 2\) symmetric hyperbolic systems

This page was built for publication: The role of numerical boundary procedures in the stability of perfectly matched layers

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2804997)