论文部分内容阅读
计算电磁学是电磁场理论、数学和计算机技术相结合的产物。计算电磁学正向着高精度、高效能、高速度的目标快速发展,很多以前无法解决的疑难电磁场问题找到了很好的解决方法,得到了精确的计算结果。越来越多的实际工程电磁场疑难问题摆在人们面前,这更促进了计算电磁学向更高水平发展。基于麦克斯韦方程组的解析方法和数值方法求解是计算电磁学的主要任务。通常,只有经典的电磁场问题才有解析解,数值计算逐渐成为解决复杂电磁场问题的主要手段,有时甚至是唯一手段。FDTD方法是典型的时域全波分析方法,应用范围非常广泛,是近些年最受关注、发展最迅速的数值方法之一。麦克斯韦方程是描述电磁现象的基本方程,FDTD方法从其旋度方程,即麦克斯韦微分形式出发,对时间和空间的一阶偏导数采取中心差分近似,直接转换为显式差分运算,这样可以在时间和空间上对连续电磁场数据实现抽样离散。FDTD方法能够描述时域电磁场的传播特性,只要给出问题的初始条件和边界条件,即可应用FDTD方法迭代递推得到各个时间步和空间步的电磁场分布。随着计算机硬件条件的发展,与FDTD方法相关的创新研究不断涌现,FDTD方法将会赢得越来越多的计算电磁学领域专家的关注和青睐。然而,FDTD方法也有其自身不可忽视的不足:一方面,采用差分法对麦克斯韦方程近似求解时,会在计算网格中引起数值色散,这种关系随数值波模的传播方向以及离散化程度不同而改变。因此,FDTD方法受限于数值色散条件,一般空间步长不大于波长的十分之一,当数值模拟对象的电尺寸较大时,将导致所需内存剧增;另一方面,FDTD方法的时间步长和空间步长不独立,Courant-Friedrich-Levy (CFL)稳定条件限制了时间离散,使得时间步长必须随空间步长变化而变化,如果空间步长变小,则需要增加时间迭代步数以实现收敛。因此,对于目前普通的PC机而言,FDTD方法在计算效率方面存在着不足。围绕着电磁场数值方法的计算效率问题,近年来出现了FDTD的多种改进算法:针对内存问题的R-FDTD方法,针对计算时间步长问题的ADI-FDTD和LOD-FDTD方法都是较热门的数值计算方法。当然,在实际应用中,会发现这些算法都或多或少存在一些问题,对这些FDTD方法的改进算法作进一步的研究,以节约内存使用量和计算时间,提高计算准确程度为目的,实现更加高效、更加精确的数值计算,具有理论和应用意义。论文首先对课题的研究背景及意义进行了阐述,包括计算电磁学进展、发展时域数值计算的必要性以及时域有限差分方法的介绍。本文涉及到的FDTD改进算法包括R-FDTD方法、ADI-FDTD方法和LOD-FDTD方法。它们能够很好的解决内存使用量大,计算时间长的问题。介绍了几种FDTD改进算法的研究现状。最后说明了本文所做的主要工作。针对具有对称结构的计算模型,从理论上分析了采用PEC边界和PMC边界截断的对称边界条件,提出了对称截断FDTD方法,利用该方法能够确定截断边界以外场分量的值,以实现截断边界处的FDTD方法计算。数值计算验证了对称截断FDTD方法的正确性和可行性。论述了R-FDTD方法中对暂存场分量边值进行补充计算的必要性。论证了三维R-FDTD求解感应电荷密度处理导体问题的方法与传统FDTD方法等价。提出了在内存使用量和计算时间上都具有明显优势的周期对称结构R-FDTD方法,该方法结合了R-FDTD方法和第二章对称截断方法的优点,将计算区域的内存使用量最多降为FDTD算法的1/6。由于每个时间步迭代计算的复杂度降低,需要计算的网格数明显减少,总的迭代计算时间也大幅缩减。数值计算验证了该方法的正确性和有效性。为了准确求解ADI-FDTD方法实现PEC边界和PMC边界的待求场分量系数,通过在获得该系数前应用理想导体边界条件,推导出了相应的修正系数。计算了单个金属立方体和对称的两个金属立方体的双站RCS。结果表明,理想导体边界作为理想导体表面,采用修正系数的计算结果与FDTD方法计算结果更为吻合;理想导体边界作为截断计算空间对称面,采用修正系数的计算结果与ADI-FDTD方法计算结果相同,与理论推导结论一致。证明了LOD-FDTD方法实现PEC边界和PMC边界时的待求场分量系数与传统的LOD-FDTD方法系数不同。通过在获得该系数前应用理想导体边界条件,得到对应的修正系数。针对将理想导体边界条件作为理想导体表面和截断计算空间对称面的不同情况,讨论了修正系数与传统LOD-FDTD系数的区别。具有统一的表达式的修正系数LOD-FDTD方法计算理想导体表面较传统LOD-FDTD方法误差更小,并对其进行了数值验证。