一. 计算电磁学的重要性
在现代科学研究中,“科学试验,理论分析,高性能计算”已经成为三种重要的研究手段。在电磁学领域中,经典电磁理论只能在11 种可分离变量坐标系中求解麦克斯韦方程组或者其退化形式,最后得到解析解。解析解的优点在于:
①可将解答表示为己知函数的显式,从而可计算出精确的数值结果;
②可以作为近似解和数值解的检验标准;
③在解析过程中和在解的显式中可以观察到问题的内在联系和各个参数对数值结果所起的作用。
这种方法可以得到问题的准确解,而且效率也比较高,但是适用范围太窄,只能求解具有规则边界的简单问题。当遇到不规则形状或者任意形状边界问题时,则需要比较复杂的数学技巧,甚至无法求得解析解。20 世纪60 年代以来,随着电子计算机技术的发展,一些电磁场的数值计算方法也迅速发展起来,并在实际工程问题中得到了广泛地应用,形成了计算电磁学研究领域,已经成为现代电磁理论研究的主流。简而言之,计算电磁学是在电磁场与微波技术学科中发展起来的,建立在电磁场理论基础上,以高性能计算机技术为工具,运用计算数学方法,专门解决复杂电磁场与微波工程问题的应用科学。相对于经典电磁理论分析而言,应用计算电磁学来解决电磁学问题时受边界约束大为减少,可以解决各种类型的复杂问题。原则上来讲,从直流到光的宽广频率范围都属于该学科的研究范围。近几年来,电磁场工程在以电磁能量或信息的传输、转换过程为核心的强电与弱电领域中显示了重要作用。
二. 电磁问题的分析过程
电磁工程问题分析时所经历的一般过程为:
三. 计算电磁学的分类
(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)计算电磁学常用方法汇总
(4) 几种主要方法之间的比较
这里对计算电磁学中几种主要的数值方法进行简单的比较,即时域有限差分法(FDTD)、有限元(FEM)、矩量法(MoM)、多极子法(MMP)、几何光学绕射法(GTD)、物理光学绕射法(PTD)和传输线法(TLM)。
性能 | MoM | GTD/PTD | MMP | FDTD | FEM | TLM |
使用求解的问题 | 天线建模、线建模和表面结构、导线结构的问题 | 大电尺寸结构的范围的应用 | 直接计算,不需要中间步骤 | 可以直接求解麦克斯韦方程 | 电的和物体几何尺寸的特性可分开定义和处理 | 所有的场分量可以在同一点进行计算 |
数值建模特点 | 可以对任意结构形状的物体上的电流结构建模 | 在高频散射问题中非常有效,例如雷达散射截面问题 | 不需要存储空间形状参数 | 可以克服FDTD中必需的阶梯建模空间问题 | 可用于非均匀煤质建模和分析 | |
适于计算电磁场的区域 | 辐射条件允许求解在辐射物体外的任何地点的E和H场 | 满足远区平面波近似的空间,节省计算机资源 | 很容易对非均匀煤质的场问题建模 | 适于分析复杂结构,对内部EM问题建模有效 | 适于分析复杂结构,对表面域建模很有效 | |
适于研究的问题 | 计算天线参数、输入阻抗、增益、雷达问题 | 对内部复杂煤质问题可以有效地建模 | 可以对非均匀煤质问题建模 | 比FDTD有较小的数值色散误差 | ||
数值建模中存在的问题 | 对内部区域建模问题困难大 | 几乎不提供有关天线参数的信息 | 场强以外的其它参数必须进行计算 | 对无边界问题需要吸收边界条件处理 | 对无边界问题需要对边界进行建模 | 比FDTD使用更多的计算资源 |
计算机实现遇到的问题 | 在非均匀煤质中会遇到困难,要用大量的内部资源,所以,通常只用于低频问题 | 只在高频有效,不能提供任何电流分布的情况 | 计算密集型,占用的计算量和内存都很大,使用者必须熟悉多极子理论 | 计算密集型,有数值色散误差,内存量大 | 计算密集型,处理开放区域内的封闭面上的未知场点问题难 | 带宽受色散误差限制,不能解围绕散射体和需要大空间的问题 |
计算场强以外的其它物理量的能力 | 只能计算远区场 | 计算场传播和电流分布等参数很难 | 同FDTD |
由于实际问题的多样性,单独使用以上介绍的方法可能并不能满足需要,比如涂敷介质的目标、印刷电路板及微带天线的辐射散射/EMC分析、带复杂腔体和缝隙结构的目标的散射等等。因此工程界常常将各种方法搭配起来使用,形成各种混合方法。常见的混合方法包括边界积分方程与体积分方程/微分方法混合、高频近似方法与低频精确方法的混合、解析方法与数值方法的混合等。
高频方法与低频方法的混合技术一般针对含有复杂细节的电大尺寸目标而提出的。由于完全使用低频的精确方法来处理电大尺寸部分往往超出了目前计算机的能力,而单纯使用高频方法又得不到足够精确的近场,所以这种分而治之的折中方案就出现了。常用的混合方法包括弹跳射线法/矩量法混合(SBR/MoM)、物理绕射理论/矩量法混合(PTD/MoM)、几何绕射理论/矩量法混合(GTD/MOM)等等。当然,引入了高频近似,赢得了速度和空间,同时在一定程度上也损失了精度。
除了上述几种混合方法之外,将解析方法和数值方法混合也是一种非常有用的方法。比如二维非均匀介质电磁问题中将二维的数值计算转化为径向本征模式展开与纵向的解析递推的数值模式匹配法(NMM)以及对于n维偏微分方程先使用(n一l)维数值离散转化为常微分方程后再用解析方法求其通解的直线法都是很好的例子。
(6) 算法的快速求解
快速算法:快速算法是为了解决矩量法求解过程中存储量和计算量过大的问题而出现的。近年来,许多学者致力于精确方法的快速求解以满足工程中日益增长的对电大尺寸复杂物体精确模拟之需要。由于矩量法产生的是一个满阵,存储量为O( N2),采用直接求解的计算复杂度为O (N3),采用迭代求解的计算复杂度为O( N2),当未知量N增大的时候,存储量和计算量都会快速增加,这极大的限制了其求解能力。而某些基于矩量法的快速算法,如多层快速多极子算法,可以成功得将存储量和计算复杂度分别降到O (N)和O (N logN)量级,极大的扩大了其求解能力。这些方法主要有基于分组思想的快速多极子方法(FMM),多层快速多极子算法(MLFMA),快速非均匀平面波算法(FIPWA),自适应积分方法(AIM),共轭梯度快速傅立叶变换(CG-FFT)等方法。
并行计算,也称之为高性能计算,则是在现有的算法基础上,增加计算资源等硬件设施,把待求解的问题分解为许多小问题,分别在不同的处理器上求解,通过网络等方式实现进程间的通信,最后得到需要的解,从而实现联合求解大问题。并行计算机从上世纪中期出现以来,出现了很多种不同的体系,主要有并行向量机(PVP),对称多处理机(SMP),大规模并行处理机(MPP),集群(Cluster),分布式共享存储多处理机(DSM)等。