Extension and validation of a perfectly matched layer formulation for the unconditionally stable D-H FDTD method

In this letter, a modification to the recently proposed unconditionally stable D-H ADI FDTD method is presented that considerably reduces the late-time error induced by the corner cells. The PML boundary is derived from the direct discretization of the modified D-H Maxwellʹs equations rather than the superposition of uniaxial PML boundaries. An optimal choice of the PML conductivity profile coefficients is proposed. Results show that the reflection error of the PML is limited for increased time step size beyond the Courant-Friedrichs-Lewy stability bound, and maximum reflection errors are 15 to 20 dB lower than the original formulation.