Kim E-K, Ha S-G, Lee J, Park Y B, Jung K-Y
Opt Express. 2015 Jan 26;23(2):873-81. doi: 10.1364/OE.23.000873.
Efficient unconditionally stable FDTD method is developed for the electromagnetic analysis of dispersive media. Toward this purpose, a quadratic complex rational function (QCRF) dispersion model is applied to the alternating-direction-implicit finite-difference time-domain (ADI-FDTD) method. The 3-D update equations of QCRF-ADI-FDTD are derived using Maxwell's curl equations and the constitutive relation. The periodic boundary condition of QCRF-ADI-FDTD is discussed in detail. A 3-D numerical example shows that the time-step size can be increased by the proposed QCRF-ADI-FDTD beyond the Courant-Friedrich-Levy (CFL) number, without numerical instability. It is observed that, for refined computational cells, the computational time of QCRF-ADI-FDTD is reduced to 28.08 % of QCRF-FDTD, while the L relative error norm of a field distribution is 6.92 %.
开发了一种高效的无条件稳定时域有限差分(FDTD)方法用于色散介质的电磁分析。为此,将二次复有理函数(QCRF)色散模型应用于交替方向隐式时域有限差分(ADI-FDTD)方法。利用麦克斯韦旋度方程和本构关系推导了QCRF-ADI-FDTD的三维更新方程。详细讨论了QCRF-ADI-FDTD的周期边界条件。一个三维数值例子表明,所提出的QCRF-ADI-FDTD方法能够使时间步长超过柯朗-弗里德里希-列维(CFL)数,而不会出现数值不稳定。可以观察到,对于细化的计算单元,QCRF-ADI-FDTD的计算时间减少到QCRF-FDTD的28.08%,而场分布的L相对误差范数为6.92%。