Unsplit-field FDTD implementations of the PML using DSP algorithms

Abstract

A simple algorithm is presented for the Finite Difference Time Domain (FDTD) implementation of the stretched coordinate Perfectly Matched Layer (PML) formulations. The algorithm is based on incorporating the Digital Signal Processing (DSP) algorithms developed for digital filters into the PML formulations. The advantage of this method is that it allows direct FDTD implementation of Maxwell's equations in the PML region without the need for splitting the field components. In addition, the formulations are implemented using D and B fields, rather than the conventional E and H. This allows the PML to be independent from the material properties of the FDTD computational domains. Numerical tests are carried out in two dimensions to validate the proposed formulations.