一篇高效电磁计算新方法的学术报告
本报告旨在向国内研究者介绍一篇发表于2025年6月《IEEE Transactions on Microwave Theory and Techniques》第73卷第6期、题为《一种用于具有内外球形边界的电磁问题的高效三维混合有限元谱积分方法》的原创性研究论文。该研究由厦门大学电磁声学研究院的孙泽飞、庄明伟（IEEE会员）以及宁波东方理工大学（筹）电子科学与技术学院、同时兼任厦门大学电磁声学研究院的刘清虎（IEEE Fellow）共同完成。
一、 研究背景与目标
本研究属于计算电磁学领域，核心目标是解决无界区域（如开放空间）中电磁波传播与散射问题的数值模拟难题。在使用有限元法进行此类模拟时，一个关键挑战是如何有效地“截断”计算域，即在有限的计算区域内模拟无限空间，并精确地施加辐射边界条件，以避免非物理的数值反射。
传统的解决方案主要包括吸收边界条件和完美匹配层。尽管这些方法应用广泛，但存在一些固有局限，例如PML需要额外的网格层，增加了计算量和内存消耗，且其精度可能受到波入射角的影响。另一种精确的途径是将有限元法与边界积分法耦合，形成FEM-BIM（或称FEBI）方法，该方法能精确满足无穷远处的辐射条件，但往往带来计算复杂度高、内存需求大的问题，尤其对于大规模三维问题。
在此背景下，谱积分法作为一种高效的边界积分方法，通过与有限元或谱元法结合，已在处理特定结构（如多层圆柱、周期性结构）的电磁散射问题上显示出优势。然而，针对三维球形或部分球形边界的通用、高效混合方法尚未见报道。因此，本研究旨在开发并验证一种新型的三维混合方法，将有限元法的区域建模灵活性与谱积分法在球形边界上处理开放问题的精确高效性相结合，以超越传统FEM-PML和现有FEM-BIM方法的性能。
二、 研究方法与流程
本研究提出并详细阐述了名为“有限元谱积分法”（FESI）的新方法。其核心工作流程可概括为以下几个关键步骤：
1. 问题建模与有限元弱形式推导： 研究针对一个包含非均匀介质、并由内部和外部两个同心球形表面界定的有限元计算域。该域之外是均匀的背景介质。研究从描述电场行为的矢量亥姆霍兹方程出发，利用伽辽金方法推导出标准的有限元弱形式。在此过程中，边界上的磁场切向分量被明确地作为未知量（等效为表面电流密度 J）引入方程，为后续与外部场方程的耦合奠定了基础。这部分工作采用了基于四面体网格的一阶Nédélec矢量基函数进行离散化，形成了包含内部电场、边界切向电场和边界等效电流密度的方程组。
2. 球形边界上的谱积分法公式化： 这是本研究的主要创新点之一。对于内部和/或外部的均匀球形区域，研究摒弃了在边界上直接进行空间域积分的传统边界元法，转而采用谱积分法。其核心是利用矢量球谐函数作为一组在球面上完备且正交的基函数，将边界上的切向电磁场（等效电流 J 和等效磁流 M）以及并矢格林函数在此基底下展开。 具体流程包括： * 场与格林函数展开： 将边界上的未知电流/磁流、已知入射场以及描述均匀介质中场传播的并矢格林函数，全部用矢量球谐函数及其相关联的球贝塞尔函数进行展开表示。 * 建立谱域方程： 将上述展开式代入电场积分方程，并利用矢量球谐函数的正交性，将在空间域上复杂的积分方程，转化为在谱域（即球谐模系数空间）上的一系列解耦的、标量的代数方程。这是SIM效率的关键：原本稠密耦合的边界积分矩阵，在谱域中变成了对角阵。 * 处理不完整球面： 研究还扩展了方法的应用范围，提出了处理不完整球形边界（如半球）的策略。虽然此时无法直接使用SIM计算外场，但通过结合数值积分和快速矢量球谐变换，仍能精确求解边界场并进而计算外部辐射场。
3. FEM与SIM的对称耦合方案： 有限元域产生的是基于局部网格基函数的离散方程，而SIM产生的是基于全局球谐基函数的谱域方程。为了将两者结合成一个统一的系统，需要建立两种表示之间的“桥梁”。本研究设计了一个对称耦合方案： * 投影矩阵： 构造投影矩阵 Y，其作用是将谱域中的球谐系数向量，转换（投影）到有限元边界上的矢量基函数空间，反之亦然。这通过计算球谐函数与有限元边界基函数在球面上的内积实现。 * 构建整体系统矩阵： 将有限元矩阵、谱积分对角矩阵以及上述投影矩阵块，按照内部未知量、边界电场未知量、边界电流未知量、谱域系数的顺序组装成一个大型线性系统。 * 对称化预处理： 由于投影操作，初始的整体矩阵是非方阵且非对称的。研究巧妙地设计了一个预条件子，作用于该系统后，得到一个近似对称的方阵系统，从而可以使用高效的标准稀疏/稠密混合求解器进行求解。这种对称耦合策略是FESI方法得以实用化的关键步骤。
4. 近场与远场计算： 一旦从耦合系统中解出所有未知量，即可高效计算任意位置的场： * 近场计算： 在有限元域内，直接利用有限元解进行插值。在内部或外部的均匀球域内，利用已求得的谱域系数和球谐函数/球贝塞尔函数的解析表达式，通过求和快速计算出任意点的散射场。这避免了传统FEM-PML方法中需要通过数值微分计算边界电流带来的误差和复杂性。 * 远场计算： 对于外部区域，利用球谐函数和球贝塞尔函数在远区的渐进形式，可以直接从边界谱域系数解析地导出远场方向图和雷达散射截面。研究集成了快速矢量球谐变换算法来加速这一过程，使得远场计算时间几乎与有限元网格密度无关，极具效率。
三、 主要结果与验证
研究通过四个详尽的数值算例，系统地验证了FESI方法的精度、效率和优越性。
1. 算例一：与解析解及COMSOL的对比验证 以一个介电球体的平面波散射问题为例，与经典的Mie级数解析解以及商业软件COMSOL的FEM-PML方法进行对比。结果显示： * 精度： 当谱采样密度设置为4点每波长时，FESI在有限元域内计算得到的近场相对均方根误差即可达到约1%，与充分收敛的FEM-PML精度相当甚至更优。在远场计算中，FESI表现出更高的精度。 * 效率： 在达到相同精度的前提下，FESI的计算时间显著低于FEM-PML。特别是在远场计算中，得益于FAVEST算法，FESI的计算时间恒定且极短，而传统方法时间随精度要求增加。 * 内存与自由度： 通过移除内部均匀球形区域的网格，FESI减少了总的计算自由度，从而降低了内存使用量。
2. 算例二：雷达罩散射分析 模拟一个包含复杂介质雷达罩和内部抛物面天线的实际工程问题。对比FESI、COMSOL (FEM-PML) 以及FEKO (FEBI和MoM-BEM) 的结果。 * 性能优势： 在相同网格密度下，FESI比COMSOL减少了约34.7%的自由度，求解时间缩短约73.5%，内存使用降低约40.1%。与FEKO的FEBI方法相比，FESI在保证精度的同时，在求解时间和内存消耗上展现出巨大优势。 * 场连续性： FESI计算的由SIM区域过渡到FEM区域的场显示出完美的连续性，而FEKO的FEBI #2（内外均使用BIM）结果在交接处出现了可见的不连续，表明了FESI耦合方案的有效性。
3. 算例三：Ku波段波纹喇叭天线 此算例展示了FESI处理不完整球形边界（半球）的能力。天线辐射边界是一个半球面。 * 有效性验证： 尽管半球边界无法直接使用SIM公式计算外场，但通过研究提出的“谱域校正”流程，FESI仍然能精确求解半球边界上的电流分布，并进而计算远场。其近场、边界电流及远场结果与COMSOL和FEKO均吻合良好。 * 边界场精度： 研究发现，COMSOL由于需要数值计算磁场以得到边界电流，在边界某些区域产生了非物理的振荡，而FESI直接求解的电流则光滑准确，凸显了其边界求解的精度优势。
4. 算例四：滤波环形耦合器频率响应分析 对一个基于半球形腔体的微波无源器件进行频域扫掠分析，计算其S参数。 * 频域性能： 在10-14 GHz范围内，FESI计算得到的S参数曲线与COMSOL结果高度一致。 * 工程应用价值： 此例证明了FESI在分析具有辐射缝隙的封装结构频率特性方面的适用性，其高效性对于需要大量频率采样点的器件优化设计具有重要意义。
四、 结论与意义
本研究成功提出并验证了一种创新的三维有限元谱积分方法。其主要结论和价值体现在： * 方法学创新： FESI方法首次实现了三维有限元法与基于矢量球谐函数的谱积分法在球形及部分球形边界上的高效、对称耦合，为解决一类常见电磁问题提供了新的精确计算工具。 * 性能卓越： 相较于广泛使用的FEM-PML方法，FESI在保证高精度的同时，能显著降低计算自由度、减少内存占用、并大幅提升求解效率，尤其在近场和远场计算方面优势明显。 * 应用范围广： 方法不仅适用于完整的球形边界（如球体散射、球形天线罩），通过扩展也能有效处理不完整的球形边界（如喇叭天线、半球形谐振器），拓宽了其工程应用场景。 * 为未来发展奠基： 该工作为将谱积分法推广至更一般的几何边界，以及处理多个离散球面的耦合问题奠定了重要的理论基础和技术储备。
五、 研究亮点
六、 其他有价值内容
论文附录提供了标量 tesseral 函数在球极坐标奇异点（θ=0, π）处的解析极限表达式，这对于确保矢量球谐函数数值计算的准确性至关重要，体现了研究的严谨性。此外，文中对采样密度选择、带宽截断准则的讨论，为实际应用该方法提供了明确的参数设置指导。