New hybrid FDTD algorithm for electromagnetic problem analysis*

Abstract

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.

Photocurrent improvement of an ultra-thin silicon solar cell using the localized surface plasmonic effect of clustering nanoparticles

“…The absorbed power, short-circuit current density and electric field distribution are calculated using the finitedifference time-domain (FDTD) method. [50] A perfectly matched layer (PML) was used for upper and lower boundary conditions and periodic boundary conditions for lateral boundaries. A plane wave source was used in the wavelength range from 0.3 µm to 1.1 µm.…”

Applicability of the Effective Index Method for the Simulation of X-Cut LiNbO3 Waveguides

“…The above methods directly solve Maxwell's equations without any approximations and can therefore provide the most accurate result. The main disadvantages of these methods are their long calculation times and their programming complexity [26,27]. During the simulation process of basic integrated optical elements, hundreds of calculation iterations are performed in order to determine the optimum parameters.…”

Two-dimensional finite element mesh generation algorithm for electromagnetic field calculation*

“…With its favorable ability to adapt the geometric shape of the field boundary and the variation of the physical properties of the medium, the finite element method has become one of the powerful numerical calculation methods for numerical analysis and optimal design of electromagnetic field. [1,2] Generally speaking, the finite element analysis is divided into three stages: modeling, calculation and solution, result processing and evaluation. Their proportions of time spent in each stage are 40%-60%, 5%-10%, 30%-50%, respectively.…”