一篇文章入门计算电磁学

计算电磁学CEM（computational electromagnetics）是笔者在研发过程中认为最复杂的物理场，难度在CFD和计算材料学之上。计算电磁学的复杂主要表现在物理场抽象，计算规模大，同时求解方法众多，涉及到大量的底层技术知识。

求解的偏微分方程是麦克斯韦方程组，麦克斯韦在奥斯特，法拉利等前人试验基础上通过数学推理得到了完整的方程组，在该方程组的理论支持下，有了后来的电磁学的飞速发展。该方程组完整的描述了电，磁，材料，频率，时间之间的关系。

求解电磁学可分为三类：解析法，数值法，以及半解析半数值。

(1) 时域方法与谱域方法

电磁学的数值计算方法可以分为时域方法(Time Domain或TD)和频域方法(Frequeney Domain或FD)两大类。

时域方法对Maxwell方程按时间步进后求解有关场量。最著名的时域方法是时域有限差分法(Finite Difference Time Domain或FDTD)。这种方法通常适用于求解在外界激励下场的瞬态变化过程。若使用脉冲激励源，一次求解可以得到一个很宽频带范围内的响应。时域方 法具有可靠的精度，更快的计算速度，并能够真实地反应电磁现象的本质，特别是在诸如短脉冲雷达目标识别、时域测量、宽带无线电通讯等研究领域更是具有不可 估量的作用。

频域方法是基于时谐微分、积分方程，通过对N个均匀频率采样值的傅立叶逆变换得到所需的脉冲响应，即研究时谐(Time Harmonic)激励条件下经过无限长时间后的稳态场分布的情况，使用这种方法，每次计算只能求得一个频率点上的响应。过去这种方法被大量使用，多半是 因为信号、雷达一般工作在窄带。当要获取复杂结构时域超宽带响应时，如果采用频域方法，则需要在很大带宽内的不同频率点上的进行多次计算， 然后利用傅立叶变换来获得时域响应数据，计算量较大；如果直接采用时域方法，则可以一次性获得时域超宽带响应数据，大大提高计算效率。特别是时域方法还能 直接处理非线性媒质和时变媒质问题，具有很大的优越性。时域方法使电磁场的理论与计算从处理稳态问题发展到能够处理瞬态问题，使人们处理电磁现象的范围得 到了极大的扩展。

频域方法可以分成基于射线的方法(Ray-based)和基于电流的方法(Current-based)。前者包括几何光 学法(GO)、几何绕射理论(GTD)和一致性绕射理论(UTD)等等。后者主要包括矩量法(MoM)和物理光学法(PO)等等。基于射线的方法通常用光 的传播方式来近似电磁波的行为，考虑射向平面后的反射、经过边缘、尖劈和曲面后的绕射。当然这些方法都是高频近似方法，主要适用于那些目标表面光滑，其细 节对于工作频率而言可以忽略的情况。同时，它们对于近场的模拟也不够精确。另一方面，基于电流的方法一般通过求解目标在外界激励下的感应电流进而再求解感 应电流产生的散射场，而真实的场为激励场与散射场之和。基于电流的方法中最著名的是矩量法。矩量法严格建立在积分方程基础上，在数字上是精确的。其实，我 们并不能判断它是一种低频方法或者是高频方法，只是矩量法所需要的存储空间和计算时间随未知元数的快速增长阻止了其对高频情况的应用，因而它只好被限定在 低频至中频的应用上。物理光学法可以认为是矩量法的一种近似，它忽略了各子散射元间的相互耦合作用，这种近似对大而平滑的目标是适用的，但是目标上含有边 缘、尖劈和拐角等外形的部件时，它就失效了。当然，对于简单形状的物体，PO法还是一个常用的方法，毕竟，它的求解过程很迅速，并且所需的存储空间也非常 少(O(N))。

(2)积分方程法与微分方程法

从求解的方程形式又可以分成积分方程法(IF)和微分方程法(DE)。IE法 与DE法相比，特点如下：(1)IE法的求解区域维数比DE法少一维，误差仅限于求解区域的边界，故精度高；(2)IE法适宜于求解无限域问题，而DE法 用于无限域问题的求解时则要遇到网格截断问题；(3)IE法产生的矩阵是满的，阶数小，DE法所产生的矩阵是稀疏的，但阶数大；(4)IE法难处理非均 匀、非线性和时变媒质问题，而DE法则可以直接用于这类问题。因此，求解电磁场工程问题的出发点有四种方式：频域积分方程(FDIE)、频域微分方程 (FDDE)、时域微分方程(TDDE)和时域积分方程(TDIE)。

计算电磁学也可以分成基于微分方程的方法(Differential Equation)和基于积分方程的方法(Integral Equation)两类。前者包括FDTD、时域有限体积法FVTD、频域有限差分法FDFD、有限元法FEM。在微分方程类数值方法中，其未知数理论上 讲应定义在整个自由空间以满足电磁场在无限远处的辐射条件。但是由于计算机只有有限的存贮量，人们引入了吸收边界条件来等效无限远处的辐射条件，使未知数 局限于有限空间内。即便如此，其所涉及的未知数数目依然庞大(相比于边界积分方程而言)。同时，由于偏微分方程的局域性，使得场在数值网格的传播过程中形 成色散误差。所研究的区域越大，色散的积累越大。数目庞大的未知数和数值耗散问题使得微分方程类方法在分析电大尺寸目标时遇到了困难。对于FEM方法，早 期基于节点(Node-based)的处理方式非常有可能由于插值函数的导数不满足连续性而导致不可预知的伪解问题，使得这种在工程力学中非常成功的方法 在电磁学领域内无法大展身手，直到一种基于棱边(Edge-based)的处理方式的出现后，这个问题才得以解决。

积分方程类方法主要包括各类基于边界积分方程(Boundary Integral Equation)与体积分方程(Volume Integral Equation)的方法。与微分类方法不同，其未知元通常定义在源区，比如对于完全导电体(金属)未知元仅存在于表面，显然比微分方程类方法少很多;而 格林函数(Green’s Function)的引入，使得电磁场在无限远处的辐射条件己解析地包含在方程之中。场的传播过程可由格林函数精确地描述，因而不存在色散误差的积累效 应。

(3) 几种主要方法之间的比较

这里对计算电磁学中几种主要的数值方法进行简单的比较，即时域有限差分法(FDTD)、有限元(FEM)、矩量法(MoM)、多极子法(MMP)、几何光学绕射法(GTD)、物理光学绕射法(PTD)和传输线法(TLM)。

(4) 多种方法的混合使用

由于实际问题的多样性，单独使用以上介绍的方法可能并不能满足需要，比如涂敷介质的目标、印刷电路板及微 带天线的辐射散射/EMC分析、带复杂腔体和缝隙结构的目标的散射等等。因此工程界常常将各种方法搭配起来使用，形成各种混合方法。常见的混合方法包括边 界积分方程与体积分方程/微分方法混合、高频近似方法与低频精确方法的混合、解析方法与数值方法的混合等。

高频方法与低频方法的混合技术一 般针对含有复杂细节的电大尺寸目标而提出的。由于完全使用低频的精确方法来处理电大尺寸部分往往超出了目前计算机的能力，而单纯使用高频方法又得不到足够 精确的近场，所以这种分而治之的折中方案就出现了。常用的混合方法包括弹跳射线法/矩量法混合(SBR/MoM)、物理绕射理论/矩量法混合(PTD /MoM)、几何绕射理论/矩量法混合(GTD/MOM)等等。当然，引入了高频近似，赢得了速度和空间，同时在一定程度上也损失了精度。

除 了上述几种混合方法之外，将解析方法和数值方法混合也是一种非常有用的方法。比如二维非均匀介质电磁问题中将二维的数值计算转化为径向本征模式展开与纵向 的解析递推的数值模式匹配法(NMM)以及对于n维偏微分方程先使用(n一l)维数值离散转化为常微分方程后再用解析方法求其通解的直线法都是很好的例 子。

(5) 算法的快速求解

快速算法：快速算法是为了解决矩量法求解过程中存储量和计算量过大的问题而出现的。近年来，许多 学者致力于精确方法的快速求解以满足工程中日益增长的对电大尺寸复杂物体精确模拟之需要。由于矩量法产生的是一个满阵，存储量为O( N2)，采用直接求解的计算复杂度为O (N3)，采用迭代求解的计算复杂度为O( N2)，当未知量N增大的时候，存储量和计算量都会快速增加，这极大的限制了其求解能力。而某些基于矩量法的快速算法，如多层快速多极子算法，可以成功的将存储量和计算复杂度分别降到O (N)和O (N logN)量级，极大的扩大了其求解能力。这些方法主要有基于分组思想的快速多极子方法（FMM），多层快速多极子算法（MLFMA），快速非均匀平面波 算法（FIPWA），自适应积分方法（AIM），共轭梯度快速傅立叶变换（CG-FFT）等方法。