Application of Spatial Filtering High Order FDTD Method in the Multi-physics Field Coupled Equations

The traditional finite-difference time-domain method (FDTD(2, 2)) is widely employed to solve the coupled equations of quantum-electromagnetic multi-physics. However, the time step of the traditional explicit FDTD(2, 2) method is limited by the Courant-Friedrichs-Lewy (CFL) condition. Besides, the FDTD(2, 2) method only achieves the second-order accuracy, which can result in significant error accumulation. To overcome this problem, we propose the SF-HO-FDTD(2, q) method with scalable time stability conditions for solving the multi-physics field coupled equations (q denotes q-th order spatial difference scheme). This approach integrates spatial filtering (SF) and higher-order FDTD (HO-FDTD) methods. Moreover, the proposed method requires no further derivation of the iterative equations of the traditional HO-FDTD(2, q) method. It only needs to introduce the SF operation in each numerical iteration process to filter out the spatial high-frequency components arising from the use of time steps that do not satisfy the CFL stability condition. This condition is applied to ensure the stability of the numerical method. Therefore, the proposed method has a high compatibility with the traditional HO-FDTD(2, q) method. Then, the numerical dispersion analysis of the SF-HO-FDTD(2, q) method is theoretically studied. Finally, numerical examples are carried out to confirm the accuracy and efficiency of the strategy suggested in this study.