On the Accuracy and Stability of the Perfectly Matched Layer in Transient Waveguides

Energy transmitted along a waveguide decays less rapidly than in an unbounded medium. In this paper we study the efficiency of a PML in a time-dependent waveguide governed by the scalar wave equation. A straight forward application of a Neumann boundary condition can degrade accuracy in computations. To ensure accuracy, we propose extensions of the boundary condition to an auxiliary variable in the PML. We also present analysis proving stability of the constant coefficient PML, and energy estimates for the variable coefficients case. In the discrete setting, the modified boundary conditions are crucial in deriving discrete energy estimates analogous to the continuous energy estimates. Numerical stability and convergence of our numerical method follows. Finally we give a number of numerical examples, illustrating the stability of the layer and the high order accuracy of our proposed boundary conditions.

[1]  Jean-Pierre Bérenger,et al.  Perfectly Matched Layer (PML) for Computational Electromagnetics , 2007, PML for Computational Electromagnetics.

[2]  D. Givoli.  High-order local non-reflecting boundary conditions: a review☆ , 2004 .

[3]  Tim Colonius,et al.  A high-order super-grid-scale absorbing layer and its application to linear hyperbolic systems , 2009, J. Comput. Phys..

[4]  Dan Givoli,et al.  Radiation boundary conditions for time-dependent waves based on complete plane wave expansions , 2010, J. Comput. Appl. Math..

[5]  T. Rylander,et al.  Perfectly matched layer in three dimensions for the time-domain finite element method applied to radiation problems , 2005, IEEE Transactions on Antennas and Propagation.

[6]  M. Grote,et al.  Efficient PML for the Wave Equation , 2010, 1001.0319.

[7]  Björn Engquist,et al.  High order difference methods for wave propagation in discontinuous media , 2006 .

[8]  Dan Givoli,et al.  High-order local absorbing conditions for the wave equation: Extensions and improvements , 2008, J. Comput. Phys..

[9]  Ya Yan Lu,et al.  Propagating modes in optical waveguides terminated by perfectly matched layers , 2005 .

[10]  L. Greengard,et al.  Nonreflecting Boundary Conditions for the Time-Dependent Wave Equation , 2002 .

[11]  B. Gustafsson.  High Order Difference Methods for Time Dependent PDE , 2008 .

[12]  Heinz-Otto Kreiss,et al.  Hyperbolic Initial Boundary Value Problems which are not Boundary Stable , 2008 .

[13]  Marcus,et al.  On Local Nonreflecting for Time Dependent Boundary Conditions Wave Propagation , 2009 .

[14]  S. Gedney.  An anisotropic perfectly matched layer-absorbing medium for the truncation of FDTD lattices , 1996 .

[15]  N. Anders Petersson,et al.  Perfectly matched layers for Maxwell's equations in second order formulation , 2005 .

[16]  Kenneth Duru,et al.  A Well-Posed and Discretely Stable Perfectly Matched Layer for Elastic Wave Equations in Second Order Formulation , 2012 .

[17]  A. Majda,et al.  Absorbing boundary conditions for the numerical simulation of waves , 1977 .

[18]  Patrick Joly,et al.  Mathematical Modelling and Numerical Analysis on the Analysis of B ´ Erenger's Perfectly Matched Layers for Maxwell's Equations , 2022 .

[19]  G. Kreiss,et al.  Stable and conservative time propagators for second order hyperbolic systems , 2011 .

[20]  Gianluca Iaccarino,et al.  Stable Boundary Treatment for the Wave Equation on Second-Order Form , 2009, J. Sci. Comput..

[21]  Gunilla Kreiss,et al.  Application of a perfectly matched layer to the nonlinear wave equation , 2007 .

[22]  Gianluca Iaccarino,et al.  Stable and accurate wave-propagation in discontinuous media , 2008, J. Comput. Phys..

[23]  Raj Mittra,et al.  Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers , 1996 .

[24]  Patrick Joly,et al.  Stability of perfectly matched layers, group velocities and anisotropic waves , 2003 .

[25]  Jan Nordström,et al.  High order finite difference methods for wave propagation in discontinuous media , 2006, J. Comput. Phys..

[26]  Marcus J. Grote,et al.  Exact Nonreflecting Boundary Conditions for the Time Dependent Wave Equation , 1995, SIAM J. Appl. Math..

[27]  Dan Givoli,et al.  Comparison of high‐order absorbing boundary conditions and perfectly matched layers in the frequency domain , 2010 .

[28]  E. A. Skelton,et al.  Guided elastic waves and perfectly matched layers , 2007 .

[29]  Kenneth Duru,et al.  Perfectly matched layers for second order wave equations , 2010 .

[30]  D. Gottlieb,et al.  Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes , 1994 .

[31]  J. Bérenger,et al.  Application of the CFS PML to the absorption of evanescent waves in waveguides , 2002, IEEE Microwave and Wireless Components Letters.

[32]  S. Gedney,et al.  On the long-time behavior of unsplit perfectly matched layers , 2004, IEEE Transactions on Antennas and Propagation.

[33]  Magnus Svärd,et al.  On the order of accuracy for difference approximations of initial-boundary value problems , 2006, J. Comput. Phys..

[34]  Marcus J. Grote,et al.  On local nonreflecting boundary conditions for time dependent wave propagation , 2009 .

[35]  Christer Johansson,et al.  Waveguide Truncation Using UPML in the Finite-Element Time-Domain Method , 2005 .

[36]  Weng Cho Chew,et al.  A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates , 1994 .