论文部分内容阅读
近年来,为了使时域有限差分(finite-difference time-domain,FDTD)算法的时间步长不再受Courant-Friedrichs-Lewy(CFL)稳定性条件的限制,提出了无条件稳定FDTD算法;可以减少计算时间,提高计算效率。近年来随着通信行业的飞速发展,通信器件和设备也趋于小型化,无条件稳定FDTD算法在数值计算中的作用也越来越重要。本论文主要研究无条件稳定FDTD算法的分析理论和数值方法,并提出了几种新型具有低数值色散和高阶精度的无条件稳定FDTD算法。同时,对新型无条件稳定FDTD算法的数值特性进行了研究。本研究课题主要创新工作内容有:(一)本论文基于split-step (SS)方案和Crank-Nicolson (CN)方案,提出了一种新型三维分步无条件稳定FDTD算法。该算法具有简洁的表达形式,沿着x、y和z三个坐标轴方向,将Maxwell矩阵分解成三个子矩阵。同时,将一个时间步长分解成三个均等子时间分步,从而将一个复杂的三维问题转换成三个简单的一维问题进行处理,减少了算法的计算复杂度,该算法称为SSCN3-FDTD。另外,将Maxwell矩阵沿着三个坐标轴的正方向和负方向,进一步分解成六个子矩阵。同时,将一个时间步长分解成六个均等子时间分步,本论文提出了一种新型无条件稳定SSCN6-FDTD算法。SSCN3-FDTD算法和SSCN6-FDTD算法同时具有二阶的时间精度和空间精度。并且SSCN3-FDTD算法和SSCN6-FDTD算法的算术因子数比ADI-FDTD算法要少,从而可以减少计算时间,提高计算效率。另外,基于高阶空间中心差分格式,本论文提出了具有高阶空间精度的SSCN6-FDTD算法。为了进一步减小数值色散误差,在矩阵分解的同时,引入了三个色散控制参数,本论文提出了一种改进SSCN6-FDTD算法。(二)本论文提出了一种具有高阶精度的二维分步分解无条件稳定FDTD算法,该算法具有四个子时间分步,记为SS4-FDTD-2D。SS4-FDTD-2D算法具有新的矩阵分解形式,与SS1-FDTD算法不同。SS4-FDTD-2D算法的归一化数值相位速度误差和归一化数值相位速度各向异性误差均比ADI-FDTD-2D算法、SS1-FDTD-2D算法和SS2-FDTD-2D算法要小。为了应用于开域的电磁问题,本论文提出了一种加载NPML吸收边界的无条件稳定SS4-FDTD-2D算法。另外,基于高阶空间中心差分格式,本论文提出了具有高阶空间精度的SS4-FDTD-2D算法。为了进一步减小数值色散误差,在矩阵分解的同时,引入了两个色散控制参数,提出了一种改进SS4-FDTD-2D算法。(三)本论文基于split-step (SS)方案,提出了三种新型具有低数值色散和高阶精度的二维无条件稳定FDTD算法。在第一种算法中,采用对称阵子和均匀分解的方法,将Maxwell矩阵分解成六个子矩阵。同时,将一个时间步长分解成六个均等的子时间分步。通过变换第一种算法中子矩阵的排列次序,得到第二、三种算法。因此,三种新型算法具有类似的结构形式,并且其归一化数值相位速度误差和归一化数值相位速度各向异性误差比ADI-FDTD-2D算法和SS1-FDTD-2D算法要小。(四)本论文将SS4-FDTD-2D算法成功扩展到三维空间,提出了一种新型具有高阶精度的三维无条件稳定FDTD算法,记为SS4-FDTD-3D。采用对称阵子和均匀分解的方法,将Maxwell矩阵分解成四个子矩阵。同时,将一个时间步长分解成四个均等子时间分步。SS4-FDTD-3D算法的归一化数值相位速度误差和归一化数值相位速度各向异性误差均比ADI-FDTD-3D算法要小。并且,SS4-FDTD-3D算法采用粗网格尺寸与ADI-FDTD-3D算法采用细网格尺寸,具有相同水平的精度。另外,基于高阶空间中心差分格式,提出了具有高阶空间精度的SS4-FDTD-3D算法。