A Stabilized Multilevel B-Spline-Based Fast Integral Method for the Solution of the Electric Field Integral Equation

技术摘要：一种用于 EFIE 的稳定多层级 B-样条基快速积分法

问题陈述
电场积分方程（EFIE）是求解电磁散射和辐射问题的标准工具，它仅需进行表面离散化。然而，标准的边界元法（BEM）离散化会导致矩阵组装、存储和矩阵-向量乘法（MVP）具有 $O(N^2)$ 的复杂度，其中 $N$ 是未知数数量。为了将复杂度降低到 $O(N)$ 或接近 $O(N)$，研究人员开发了各种快速积分方法，例如多层级快速多极子算法（MLFMM）以及基于 FFT 加速的自适应积分法（AIM）。

在依赖全局拉格朗日插值的 FFT 加速方法中存在一个特定的局限性。这些方法需要等距插值点以利用 FFT 结构。然而，等距拉格朗日插值在高阶多项式时会遭受 Runge 型不稳定性，导致精度崩溃，从而阻碍了稳健的高阶压缩。虽然 B-样条可以解决等距网格上的 Runge 不稳定性问题，但将其应用于多层级快速积分方法并非易事。与拉格朗日多项式不同，B-样条系数并不与插值点处的函数值重合，因此需要对采样矩阵进行求逆。这些矩阵在高阶情况下会变得严重病态，从而损害数值稳定性。此外，在确保稳定性的同时维持精确的层间传递（这是多层级效率所必需的）一直是一个重大挑战。

方法论
作者提出了一种稳定的多层级 B-样条基快速积分法，该方法结合了兼容 FFT 的核函数插值与稳健的高阶插值。其核心方法论包含三个主要部分：

1. B-样条核函数插值：
   该方法采用定义在等距节点向量上的 B-样条，而非全局拉格朗日多项式。这一选择既减轻了 Runge 型不稳定性，又保持了与 FFT 加速所需的等距采样结构的兼容性。格林函数使用张量积多项式进行分解。对于 B-样条，插值权重并非格林函数的直接采样，而是通过采样矩阵 $B$ 推导而得。B-样条的张量积结构使得对采样矩阵的求逆可以通过 $O(s^3)$ 次运算（其中 $s$ 是样条数量）高效完成，而非朴素的 $O(s^6)$，从而保持了 MVP 的整体 $O(s^3 \log s)$ 复杂度。

2. 通过节点移除实现的稳定性：
   为了解决高阶情况下 B-样条采样矩阵的病态问题，作者引入了一种“节点移除”（knot-removal）稳定策略。经验分析表明，移除等距节点向量边界附近的特定节点，可以在不破坏多层级方案所需的等距结构的前提下，显著降低稳定性常数（即逆采样矩阵的范数）。

   • 缓冲节点（Buffer Knots）： 为了处理超出其相关单元（box）的基函数（这是 BEM 中的常见问题），该方法通过增加“缓冲节点”来扩展节点向量。这确保了所有基函数的作用范围都包含在插值区域内。
   • 精确层间传递： 多层级方法的一个关键要求是在父级和子级单元之间精确地传递矩（moments）。所提出的稳定策略保留了节点向量的嵌套结构。通过使用“节点插入”（一种标准的样条细分技术），该方法构建了层级间的精确传递矩阵。这确保了在稳定化（移除边界节点）后的网格上计算的矩，可以精确地映射到用于核函数插值的原始等距网格上，从而在各层级间保持机器精度。

3. 分解策略：
   本文介绍了一种用于离散化 EFIE 的新型“基于星型”（star-based）分解法，同时也回顾了现有的直接分解和偶极分解法。

   • 基于星型的分解： 该方法利用循环-星型（loop-star）分解概念来表示标量势算符。它将标量势算符分解为一个稀疏的星型矩阵和一个单层算符。这种策略可以将标量势算符的零空间（null space）在机器精度范围内予以保留，这对于低频预处理技术（如 Calderón 预处理）的有效性至关重要。该方法需要同时计算 RWG 和分片常数矩。
   • 直接与偶极分解： 这些作为替代方案被进行了回顾。直接分解直接对核函数进行插值，而偶极形式则使用偶极核。文中强调了基于星型的分解在处理零空间方面具有更优越的稳定性。

主要贡献

• 稳定的 B-样条插值： 引入了节点移除策略，消除了与高阶等距 B-样条插值相关的 Runge 现象和病态问题，从而实现了稳健的高精度压缩。
• 精确的多层级传递： 开发了一种利用节点插入实现精确层间传递的机制，即使在应用了稳定性处理（节点移除）的情况下也是如此。这解决了稳定插值与精确层级传递之间的冲突。
• 基于星型的分解： 提出了一种新的 EFIE 分解策略，能够将标量势算符的零空间保持在机器精度水平，从而增强了对低频稳定化和预处理的兼容性。
• 复杂度保持： 证明了该方法在 MVP、组装和存储方面均能保持 $O(N)$ 的复杂度，这得益于 B-样条的紧支撑特性和张量积结构。

结果
针对典型几何体（球体）和真实几何体（飞机），在各种电尺寸（从 $0.2\lambda$ 到 $4\lambda$）下进行了数值实验。

• 稳定性与精度： 基于 B-样条的方法展示了稳定的高阶插值，避免了拉格朗日方法中观察到的崩溃现象。误差分析显示，随着多项式阶数 $p$ 和样条数 $s$ 的增加，误差收敛符合理论界限。稳定性策略成功控制了高阶情况（$p=5, 6$）下的误差，而拉格朗日方法则表现出不稳定性。
• 分解性能： 在三种测试策略中，基于星型的分解表现出了最稳健的误差行为。
• 可扩展性： 性能分析确认了对于球体和飞机模型，存储、组装时间和 MVP 时间均呈现 $O(N)$ 缩放特性。对于拥有约 300 万个未知数的球体，内存消耗在 91 GB 到 182 GB 之间（取决于 $p$），MVP 时间在 7.6 秒到 54 秒之间。对于高达约 400 万个未知数的飞机模型，也观察到了类似的线性缩放。

意义
本文声称，所提出的稳定多层级 B-样条方法为快速积分方程求解器提供了一种实用的、高精度的拉格朗日插值替代方案。通过克服 Runge 现象和采样矩阵的病态问题，该方法实现了此前在等距 FFT 加速方案中无法实现的稳健高阶压缩。通过基于星型的分解保留标量势的零空间，进一步提升了该方法在低频问题中的实用价值。作者指出，初步结果已在会议上发表，该方法可以与高频技术（如 MLFMM）结合使用，以适用于宽带表述，但其主要贡献在于稳定了低中频核函数插值的框架。