[论文解读] Recent Advances of Isogeometric Analysis in Computational Electromagnetics

Zeger Bontinck, Jacopo Corno+9 more|arXiv (Cornell University)|Sep 18, 2017

Advanced Numerical Analysis Techniques被引用 9

一句话总结

本文提出等几何分析（IGA）作为一种高精度、几何精确的计算电磁学方法，采用B样条和NURBS基函数，将CAD几何与有限元离散化统一。该方法在形状优化与不确定性量化工作流中展现出优越的收敛性与鲁棒性，关键结果表明：在2D与3D基准测试中，平均磁感应强度梯度提升17.5%，不均匀性因子降低74.4%。

ABSTRACT

In this communication the advantages and drawbacks of the isogeometric analysis (IGA) are reviewed in the context of electromagnetic simulations. IGA extends the set of polynomial basis functions, commonly employed by the classical Finite Element Method (FEM). While identical to FEM with N\\'ed\\'elec's basis functions in the lowest order case, it is based on B-spline and Non-Uniform Rational B-spline basis functions. The main benefit of this is the exact representation of the geometry in the language of computer aided design (CAD) tools. This simplifies the meshing as the computational mesh is implicitly created by the engineer using the CAD tool. The curl- and div-conforming spline function spaces are recapitulated and the available software is discussed. Finally, several non-academic benchmark examples in two and three dimensions are shown which are used in optimization and uncertainty quantification workflows.

研究动机与目标

• 解决电磁模拟中计算成本高、网格生成复杂的问题，特别是在工业应用中，80%的分析时间耗费在网格生成上。
• 通过利用原生CAD基函数，克服经典FEM在处理复杂几何与高阶离散化时的局限性。
• 通过控制点同时变形几何与计算网格，实现形状优化与不确定性量化工作流的无缝集成。
• 证明将等几何分析与边界元方法结合以降低问题维度并消除体网格化需求的可行性与优势。
• 全面回顾旋度与散度相容的样条函数空间及其在麦克斯韦方程组中的应用，确保谱正确离散化。

提出的方法

• 采用B样条与NURBS基函数——这些函数原生存在于CAD工具中——来表示几何与场变量，实现对曲线与曲面（包括圆锥曲线）的精确参数化。
• 实现旋度与散度相容的样条函数空间，以确保与麦克斯韦方程组中微分形式的离散相容性，保持重要的物理特性。
• 在等几何框架内采用伽辽金变分格式，仅需修改基函数例程即可直接复用现有FEM代码库。
• 采用仅边界离散化的等几何边界元方法（IGABEM），降低问题维度，避免三重体参数化需求。
• 通过全局C1连续的样条空间构建高阶、光滑的解场，提升单位自由度的收敛速率与精度。
• 利用高效的压缩技术处理边界元公式中产生的密集线性系统，实现在标准硬件上的高效求解。

实验结果

研究问题

• RQ1 等几何分析能否通过直接使用CAD几何显著减少电磁模拟中的网格生成工作量？
• RQ2 与经典多项式FEM相比，基于NURBS的基函数在电磁问题中如何提升精度与收敛性？
• RQ3 IGA在通过几何与网格同步变形支持形状优化与不确定性量化工作流方面，其适用范围与程度如何？
• RQ4 将等几何分析与边界元方法结合在电磁问题中具有哪些优势与挑战？
• RQ5 等几何方法如何确保与麦克斯韦方程组微分形式结构的离散相容性？

主要发现

• 由于B样条基函数具有更高的整体正则性，IGA在单位自由度上的收敛速度优于经典FEM。
• 该方法可精确表示复杂几何，如圆锥曲线，直接从CAD导入，无需几何近似。
• 在3D形状优化基准测试中，使用IGA优化后，平均磁感应强度梯度提升12.2%，不均匀性因子降低60.0%。
• 在2D基准测试中，平均磁感应强度梯度提升17.5%，不均匀性因子降低74.4%，表明在优化工作流中性能更优。
• 等几何边界元方法（IGABEM）生成更小、更密集的线性系统，可在标准台式计算机上求解，且具备高效的压缩技术以支持可扩展性。
• IGA通过控制点同时变形几何与网格，实现稳健的形状敏感性分析与优化，无需重新生成网格。