The complementary derivatives method: a second-order accurate interpolation scheme for non-uniform grid in FDTD simulation

Non-uniform grid in finite-difference time-domain methods, which is typically used to resolve fine structures, can reduce the computational domain and therefore lead to a reduction of the computational cost. However, for high-accuracy problems, such as partially-filled parallel plate waveguide and resonators, using different grid size increases the truncation error at the boundary of domains having different grid size. To address this problem, in this work, we introduce the complementary derivatives method (CDM). Theoretical discussion and numerical results will be presented to show that the CDM can maintain second-order accuracy throughout the computational domain.