无条件稳定和有条件稳定混合的时域谱元电磁分析方法
【技术领域】
[0001] 本发明属于有条件稳定和无条件稳定混合的时域方法，特别是一种针对多尺度电 磁问题的快速分析技术。
【背景技术】
[000引1864年，麦克斯韦（Maxwell)在前人的研究基础上，提出了能够对宏观电磁场基 本归路进行概括的数学方程组一著名的麦克斯韦方程，由此奠定了电磁学理论研究的基 础。电磁场理论的研究渐渐渗透到地学、生命科学、医学、材料科学和信息科学等众多技术 科学领域，大大促进了科学技术的发展和人类生活的变化。
[0003] 早期很长一段时间，电磁场理论的研究致力于得到一些问题的解析解，然而完全 用解析方法求解的问题是十分有限的，不能解决什么问题。于是，为了解决科学技术中的电 磁场问题，又发展了一些近似方法和数值方法。但是，限于当时的计算条件，无法充分发挥 该些方法的作用，令某些问题得不到实质性的解决。随着电子计算机技术的飞速发展，W高 性能计算机技术为手段，结合电磁场理论和计算数学提供的各种数值方法，应运而生了一 口交叉学科一计算电磁学。
[0004] 使用计算电磁学对电磁现象进行分析时，首先根据分析对象的特点建立相应的电 磁、数学模型。然后，选择合适的算法并在计算机上实现。当前计算电磁学方法按解域分 可分为频域方法和时域方法。频域方法主要有：W电磁场问题的积分方程为基础的矩量法 (MOM)和基于变分原理的有限元法（FEM)等；时域方法主要有：时域有限差分法（抑TD)、时 域有限元法（FETD)、时域积分法（TDIE)和时域伪谱方法（PSTD)等。
[0005]时域谱元法（Joon-Ho Lee and Qing Huo Liu, "A3-D Spectral-Element Time-Domain Method for Electromagnetic Simulation,"IEEE Transactions on Microwave Theory and "Techniques. , vol. 55, no. 5, pp. 983-991, May2007)可认为是一种特 殊的时域有限元法，由于时域谱元法采用的差分方式为中也差分，系数矩阵只含有质量矩 阵，又由于该方法中所选用的基函数为正交基函数，所W系数矩阵是块对角的，矩阵求逆会 变的很容易，与时域有限元方法比，该将大大减少计算时间。时域谱元法对网格采用的离散 方式为曲六面体离散，该能很好地拟合各种复杂的电磁结构，又因为时域谱元法网格离散 的尺寸可W很大，与时域有限差分方法相比，该将大大减少计算的未知量。
[0006] 不连续迦迂金方法（Discontinuous Galerkin, DG)的发展已有了显著的进展，该 些方法适用于具有复杂结构和不均匀煤质的大尺寸问题。该些方法很大程度上得益于80 年代上半时有里德和希尔提出的中子传输方程的求解。该些方法的最重要的特点是允许基 函数（因此，数值解）在不同单元的交界面上不连续。在每个单元上引入一套局部的基函数 而不是在整个计算区域，且不同类型的单元例如六面体，棱柱体或者四面体可W在模型中 共存，同时可W允许两边的方程使用不同的差分格式，该样计算起来更加灵活方便。
[0007] 现有的分析多尺度电磁问题的时域谱元法主要存在W下两个问题：
[0008] (1)采用六面体跟四面体混合剖分，但是六面体部分采用时域谱元法，四面体部分 采用时域有限元方法，而时域有限元方法生成的质量矩阵为稀疏症，求逆花费大量时间再 加上时间步长设置的很小，所W求解速度较慢。
[0009] (2)如果不同区域均采用六面体网格剖分，使用时域谱元法的时候，为保证算法稳 定性，时间步长必须按照最小的网格尺寸设定，该样整体的时间迭代步数就会很多，导致求 解时间慢。
【发明内容】
[0010] 本发明的目的在于提供一种无条件稳定和有条件稳定混合的时域谱元电磁分析 方法，从而快速分析复杂多尺度电磁问题。
[0011] 实现本发明目的的技术解决方案为：一种无条件稳定和有条件稳定混合的时域谱 元电磁分析方法，步骤如下：
[0012] 第一步，对所要分析的电磁问题进行几何建模，将整体模型采用曲六面体进行剖 分，剖分之后得到各个体单元的顶点编号、坐标W及体的编号；
[0013] 第二步，采用比较边长的方法找出剖分得到的网格中尺寸小于设定值的六面体， 将其标记为小尺寸区域，剩余网格标为大尺寸区域；
[0014] 第H步，将电场值定义在整体模型剖分后的网格中的每个点上，并用时域谱元法 中的化L多项式作为矢量基函数对电场在XYZ H个方向进行展开，代入时域波动方程，并采 用伽迂金测试，即测试基函数与展开基函数相同，得到矩阵方程。
[0015] 第四步，将方程中的时间项用时间差分展开，在标记出来的小尺寸的区域采用具 有无条件稳定的Newmark-目差分格式，其他区域采用有条件稳定的中也差分格式，在进行 时间迭代时整体时间步长按中也差分区域设定，根据总的时间步数每一步直接求解，最终 求得时域电场值。
[0016] 步骤一中剖分采用的曲六面体单元边长为1/10A，A为电磁波波长。
[0017] 步骤二中设定值的选取与剖分后得到网格的尺寸有关，既要保证小于此设定值的 网格的数目占整体网格数目的比例尽可能小，又要保证大尺寸区域中最小网格的尺寸跟整 个区域最小网格的尺寸之比尽可能大。
[0018] 步骤H中使用的化L基函数形式如下：
[002引其中，
，j = 0, 1，LN，LnU)是 N 阶 Legen化e多项式，将S G [-1，1]内的节点U J，j = 0, 1，LN}作为化L积分点，它们是方 程式0 -苗斯.冷）=0的（化1)个根；
[0023] 将电场用基函数展开，代入到时域波动方程
[0024] 采用伽迂金法测试，即测试函数与基函数相同，获得总系数矩阵，求解方程
[00巧]
[002引步骤四中，在小尺寸的区域采用Newmark-目差分格式后方程为：
[0027] ([T]+A *2目[S])e。"=(2 [T] - A *2 (1-2目）[S])e。-([T]+A *2目[S])
[0028] 在大尺寸区域采用中也差分格式后方程变为：
[002引[T] e。"=(2 [T] - A*2[S])e。- [T]
[0030] 求解方程，在每一步时间迭代中，先求解中也差分区域处的电场，再求解 Newmark-目差分区域的电场，最终得到总的时域电场值。
[0031] 本发明与现有技术相比，其显著优点；（1)将模