常用数值计算方法应用研究及仿真软件简介a.docVIP

常用数值计算方法 自1864年麦克斯韦建立电磁场基本方程以来，电磁波理论与应用的发展已经过了100多年的历史。对电磁分布边值问题的求解从图解、模拟、解析到目前所采用的数值计算方法，经历了四个过程。解析方法只能解决一些经典问题，具体到复杂的实际环境，往往需要通过数值解得到具体环境中的电磁波特性。随着高速和大容量计算机技术的飞速发展，电磁数值计算已经发展成为一门新兴的重要学科，已提出多种实用有效的求解麦克斯韦方程的数值方法，主要有矩量法（MM）（FEM）（FIT）（FDTD）基于这些数值计算方法开发出了许多优秀的电磁仿真软件。表 13所示。针对本论文所应用到的方法，下面简要叙述常用的几种数值方法及相应的商业软件。 表 13 数值算法 频域方法 积分形式 矩量法(MOM) 微分形式 有限差分法(FDM) 有限元法(FEM) 时域方法 积分形式 有限积分法(FIT) 微分形式 时域有限差分法(FDTD) 高频近似 绕射几何理论(GTD) 一致性绕射(UTD) 一致性渐进(UAT) 绕射物理理论(PTD) 绕射谱理论(STD) 有限元法 基于有限元方法(FEM)计算电磁问题，其基本构想是将由偏微分方程表征的连续函数所在的封闭场域划分成有限个小区域，每个小区域用一个选定的近似函数来代替，于是整个场域上的函数被离散化，由此获得一组近似的代数方程，并联立求解，以获得该场域中函数的近似数值。 广义的来说，三维麦克斯韦方程是三维电磁问题的三维支配方程，但是，一般情况下为了方便求解和建模，大多选取由麦克斯韦方程组的前两个旋度方程导出的电场强度满足矢量亥姆赫兹方程作为支配方程。如Ansoft HFSS软件[]的支配方程为： (11) 由变分原理，上式的泛函可以写为： (12) 将这一个三维问题的泛函通过多面体离散成单元小矩阵，矩形块、四面体和六面体等都可以被选用做基本的离散单元，但是，不同离散单元对于有限元运算的精度、速度和内存需求都有不同。Ansoft HFSS软件采用四面体作为基本离散单元，如图 1.4180 年代以后才被应用于电磁学中的棱边元作为矢量基函数。 假设图 1.41能够近似为 (13) 图 1.41()来表示： (14) 式中插值函数为 (15) 而有下列等式获得： (16) (17) (18) (19) 其中 (110) 在利用变分原理和离散化方法建立了有限元矩阵方程后，我们就面临着求解以结点值为未知数的矩阵方程。将方程写为 (111) 式中系数矩阵A是一个n×n方阵，x是待求解的未知量，b表示已知向量。求得这个矩阵方程得到问题空间的电磁场解，并进而求得所需参数，例如散射参数等。 有限积分法 有限积分法(FIT)法早在1977年由托马斯· 魏兰特教授 (Prof. Thomas Weiland) []提出，进而成为其后在电磁仿真领域中一个重要算法的基石。由FIT所导出的矩阵方程保持了解析麦克斯韦方程各种固有的特性，如：电荷守恒性和能量守恒性。解析下的梯度、散度和旋度算子在FIT下具有一一对应的矩阵。这些矩阵满足解析形式下的算子恒等式。故FIT保证了非常好的数值收敛性。另一个区别于其它算法的关键之处在于FIT可被用于所有频段的电磁仿真问题中。软件CST Microwave Studio是由德国CST公司开发的一种商用电磁仿真软件，主要是是采用这种数值计算方法。 CST空间离散化也是建立在Yee网格基础之上，典型划分方法如图 1.42示[]： 图 1.42 麦克斯韦积分方程离散化图 通过这种离散方法，导出相对应的麦克斯韦网格方程： (112) (113) (114) (115) 经过这些步骤，将积分方程转化为线形方程组来求解，得出问题