Unconditionally stable D-H FDTD formulation with anisotropic PML boundary conditions

An unconditionally stable finite difference time domain (FDTD) method based on a D-H formulation and the recently proposed alternating-direction-implicit (ADI) marching scheme is presented. The advantage of the D-H algorithm over the conventional E-H is the possibility to easily implement an unsplit field components formulation of the PML absorbing boundary condition that is independent from the background material used in the FDTD grid. The method allows therefore immersing any dielectric in the PML layers without any special consideration, and is amenable for models truncation often used in biomedical simulations. Furthermore, the proposed scheme can be extended to account for frequency dispersive dielectrics.