基于非零散度关系的交替方向隐式减缩FDTD算法

摘要

该文证明了即使在无源区域,交替方向隐式时域有限差分法(ADI-FDTD)所给出的电磁场量不满足零散度关系,同时推导出了该散度关系的具体表达式.基于该非零散度关系,将不受Courant稳定条件限制的ADI-FDTD法和能节约最多达1/3内存的减缩时域有限差分(R-FDTD)法结合,提出了一种新的交替方向隐式减缩FDTD算法.该算法保留了ADI-FDTD能增大时间步长,缩短计算时间的优点,同时与ADI-FDTD相比节约了最多达1/3(三维)或2/5(二维)的内存.与基于零散度关系的ADI/R-FDTD相比,该算法避免了采用长时间步长计算时的发散现象.应用所提出的ADI/R-FDTD算法计算了二维自由空间波的传播及一维频率选择表面垂直入射的问题,计算结果与ADI-FDTD计算结果完全一致,验证了ADI/R-FDTD的正确性和有效性.