电磁学计算方法的比较.docxVIP

电磁学计算方法的比较 | [] [] 摘 要：介绍了电磁学计算方法的研究进展和状态，对几种富有 代表性的算法做了介绍，并比较了各自的优势和不足，包括矩量法、 有限元法、时域有限差分方法以及复射线方法等。 关键词：矩量法；有限元法；时域有限差分方法；复射线方法 1 引 言 1864 年 Maxwell 在前 人 的理 论 （ 高斯 定 律、安 培定 律、 法拉 第 定律和自由磁极不存在）和实验的基础上建立了统一的电磁场理论， 并用数学模型揭示了自然界一切 宏观电磁现象所遵循的普遍规律， 这就是 著名 的 Maxwell 方程 。在 11 种可 分离 变量 坐 标系 求解 Maxwel 方程组或者其退化形式，最后得到解析解。这 种方法可以得到问题的准确解，而且效率也比较高，但是适用范围太窄，只能求解具有规则边界的简单问题。对于不规则形状或者任意形状边界则需要比较高的数学 技巧，甚至无法求得解析解。20 世纪 60年代以来，随着电子计算机技术的发展，一些电磁场的数值计算方法发展起来，并得到广泛地应用，相对于经典电磁理论 而言，数值方法受边界形状的约束大为减少，可以解决各种类型的复杂问题。但各种数值计算方法都有优缺点，一个复杂的问题往往难以依靠一种单一方法解决，常 需要将多种方法结合起来，互相取长补短，因此混和方法日益受到人 们的重视。 本文综述了国内外计算电磁学的发展状况，对常用的电磁计算方 法做了分类。 电磁场数值方法的分类 电磁学问题的数值求解方法可分为时域和频域 2大类。频域技术 主要有矩量法、有限差分方法等，频域技术发展得比较早，也比较成 熟。时域法主要有时域差分技 术。时域法的引入是基于计算效率的 考虑，某些问题在时域中讨论起来计算量要小。例如求解目标对冲激 脉冲的早期响应时，频域法必须在很大的带宽内进行多次采 样计算， 然后做傅里叶反变换才能求得解答，计算精度受到采样点的影响。若 有非线性部分随时间变化，采用时域法更加直接。另外还有一些高频 方法，如 GTD， UTD和射线理论。 从求解方程的形式看，可以分为积分方程法（IE）和微分方程法 DE） 。IE 和 DE相比，有如下特点：IE 法的求解区域维数比 DE法 少一维，误差限于求解 区域的边界，故精度高；IE 法适合求无限域问题，DE法此时会遇到网格截断问题；IE 法产生的矩阵是满的，阶数小，DE法所产生的是稀疏矩阵，但阶数大； IE 法难以处理非均匀、非线性和时变媒质问题，DE法可直接用于这类问题〔1〕。 几种典型方法的介绍 有限元方法是在 20 世纪 40 年代被提出，在 50 年代用于飞机设 计。后来这种方法得到发展并被非常广泛地应用于结构分析问题中。 目前，作为广泛应用于工程和数学问题的一种通用方法，有限元法已 非常著名。 有限元法是以变分原理为基础的一种数值计算方法。其定解问题为： 应用变分原理，把所要求解的边值问题转化为相应的变分问题， 利用对区域 D的剖分、插值，离散化变分问题为普通多元函数的极值 问题，进而得到一组多元的代数方程组，求解代数方程组就可以得到 所求边值问题的数值解。 （１）时域有限差分方法 时域有限差分（ FDTD）是电磁场的一种时域计算方法。传 统上电 磁场的计算主要是在频域上进行的，这些年以来，时域计算方法也越 来越受到重视。他已在很多 方面显示出独特的优越性，尤其是在解 决有关非均匀介质、任意形状和复杂结构的散射体以及辐射系统的电 磁问题中更加突出。FDTD法直接求解依赖时间变量的 麦克斯韦旋度 方程，利用二阶精度的中心差分近似把旋度方程中的微分算符直接转 换为差分形式，这样达到在一定体积内和一段时间上对连续电磁场的 数据取样压 缩。电场和磁场分量在空间被交叉放置，这样保证在介 质边界处切向场分量的连续条件自然得到满足。在笛卡儿坐标系电场 和磁场分量在网格单元中的位置是每一磁 场分量由 4个电场分量包 围着，反之亦然。 这种电磁场的空间放置方法符合法拉第定律和安培定律的自然 几何结构。因此 FDTD算法是计算机在数据存储空间中对连续的实际 电磁波的传播过程在时间进程上 进行数字模拟。而在每一个网格点 上各场分量的新值均仅依赖于该点在同一时间步的值及在该点周围 邻近点其他场前半个时间步的值。这正是电磁场的感应原理。这 些 关系构成 FDTD法的基本算式，通过逐个时间步对模拟区域各网格点 的计算，在执行到适当的时间步数后，即可获得所需要的结果。 在 上述 算法 中， 时 间增 量 t 和 空间 增量 x ， y 和 z 不 是 相互独立的，他们的取值必须满足一定的关系，以避免数值不稳定。 这种不稳定表现为在解显式 差分方程时随着时间步的继续计算结果 也将无限制的 67 增加。为了保证数值稳定性必须满足数值稳定条件： 用差分方法对麦克