New hybrid FDTD algorithm for electromagnetic problem analysis

Since the time step of the traditional finite-difference time-domain (FDTD) method is limited by the small grid size, it is inefficient when dealing with the electromagnetic problems of multi-scale structures. Therefore, the explicit and unconditionally stable FDTD (US-FDTD) approach has been developed to break through the limitation of Courant-Friedrich-Levy (CFL) condition. However, the eigenvalues and eigenvectors of the system matrix must be calculated before the time iteration in the explicit US-FDTD. Moreover, the eigenvalue decomposition is also time consuming, especially for complex electromagnetic problems in practical application. In addition, compared with the traditional FDTD method, the explicit US-FDTD method is more difficult to introduce the absorbing boundary and plane wave. To solve the drawbacks of the traditional FDTD and the explicit US-FDTD, a new hybrid FDTD algorithm is proposed in this paper. This combines the explicit US-FDTD with the traditional FDTD, which not only overcomes the limitation of CFL condition but also reduces the system matrix dimension, and introduces the plane wave and the perfectly matched layer (PML) absorption boundary conveniently. With the hybrid algorithm, the calculation of the eigenvalues is only required in the fine mesh region and adjacent coarse mesh region. Therefore, the calculation efficiency is greatly enhanced. Furthermore, the plane wave and the absorption boundary introduction of the traditional FDTD method can be directly utilized. Numerical results demonstrate the effectiveness, accuracy, stability, and convenience of this hybrid algorithm.