An analytical study of the discrete perfectly matched layer for the time-domain maxwell equations in cylindrical coordinates

Abstract

We present an analysis of the perfectly matched layer in cylindrical coordinates discretized with a staggered second-order accurate finite difference time domain method. For fixed discretization parameters, layer width, and a quadratic loss function, we find the numerical reflection produced by the discrete layer is accurately predicted by the infinite resolution reflection coefficient for /spl sigma//sub max//spl isin/[0,/spl sigma//sub max//sup c/], where /spl sigma//sub max/ is the maximum value of the absorption parameter in the layer. We also find that the finite resolution reflection coefficient achieves its minimum value at a /spl sigma//sup m//sub max/>/spl sigma//sub max//sup c/. Numerical experiments validate the analysis.