Rapidly Convergent Finite-Element Domain Decomposition Method With Two-Channel Transmission Conditions

技术摘要：具有双通道传输条件的快速收敛有限元领域分解法

问题陈述
有限元法（FEM）广泛用于求解由麦克斯韦方程组控制的电磁（EM）边值问题。然而，对于电尺寸较大的问题，FEM 会产生包含数百万或数十亿个未知数的稀疏线性系统，这给直接解法和迭代解法都带来了巨大压力。领域分解法（DDM），特别是有限元撕裂与互连（FETI）系列，通过将定义域划分为若干子域并在子域界面上迭代求解全局界面方程来解决这一问题。

FETI-DP（对偶-原初）变体中的一个关键瓶颈是施加在子域界面上的传输条件（TC）。目前事实上的标准——Robin TC（FETI-DP-2λ）——能有效阻尼传播误差模态，但无法阻尼瞬态（evanescent）模态。因此，在处理涉及大量子域数量、深瞬态波、共振和多尺度特征的问题时，收敛速度会变慢或停滞。虽然二阶传输条件（SOTCs）试图解决瞬态模态问题，但它们会在角点边缘引入冗余（需要惩罚项），并且在深瞬态机制下阻尼性能会恶化。这些局限性促使从业者转向直接解法，从而牺牲了 DDM 所固有的巨大并行性和可扩展性。目前的开放性挑战在于开发一种能够在不牺牲并行性或准确性的前提下，针对多样化问题设置保持鲁棒迭代收敛性的 FETI-DDM。

方法论
作者提出了一种新型的对偶-原初 FETI-DP 方法（FETI-DP-4λ），该方法基于直接从混合 E-B（电场-磁通） 公式中的一阶旋度方程（法拉第和安培-麦克斯韦方程）推导出的双通道传输条件。

1. 混合公式： 该方法并未将系统简化为二阶矢量波方程，而是保留 $E$ 和 $B$ 作为主要未知数，并使用 Whitney 棱边和面元进行扩展。这导致局部块系统中的界面数据由两组切向配对组成：$\hat{n} \times E$ 和 $\hat{n} \times H$。
2. 双通道 TC： 所提方法在每个界面构建了两个不同的通道：
   • 法拉第通道： 将 $\hat{n} \times E$ 与法向通量 $\hat{n}(\hat{n} \cdot B)$ 配对。
   • 安培-麦克斯韦通道： 将 $\hat{n} \times H$ 与法向通量 $\hat{n}(\hat{n} \cdot D)$ 配对。
   • 每个通道分配有各自的拉格朗日乘子（对偶变量），导致每个界面有四个拉格朗日乘子（每侧两个）。
3. 机制： 通过为每个旋度方程分配专用通道，该方法消除了 Robin TC 的“单系数障碍”。这使得能够对 TE 型和 TM 型瞬态族进行单调阻尼，包括在 Robin TC 停滞且 SOTC 回弹的深瞬态机制下的模态。
4. 系统组装： 通过消除子域未知数和原初角点未知数来组装全局界面方程。所得系统使用未预处理的 GMRES 进行求解。该方法每个子域仅需一次稀疏分解和一个用于角点未知数的小规模粗网格求解，保留了标准 FETI-DP 的计算结构。

核心贡献

• 新型传输条件： 引入了双通道 TC，同时强制执行切向场连续性（$\hat{n} \times E$）和法向通量连续性（$\hat{n} \cdot B$）。
• 鲁棒收敛： 该方法对于具有深瞬态波、多尺度特征和复杂划分特征的问题，实现了快速、无停滞的收敛，克服了 Robin TC 和 SOTC 的局限性。
• 保持可扩展性： 该方法保留了 DDM 的大规模并行性和计算可扩展性，同时显著减少了所需的迭代次数。
• 增强准确性： 通过恢复场与通量的连续性，该方法将界面跳跃降低了数个数量级，从而提高了求解精度，尤其是在松弛求解器容差的情况下。

数值结果
在采用共形四面体网格的 3-D 波导几何结构（包括 WR-90、脊加载、L 型和 X 波段带通滤波器）上进行了广泛的数值实验。将所提方法与基于 Robin TC 和 SOTC 的方法进行了对比。

• 迭代次数减少：
  • 在均匀划分的 WR-90 波导中，与 Robin TC 相比，所提 TC 在 $10^{-6}$ 容差下减少了 3.42 倍的 GMRES 迭代次数，在 $10^{-10}$ 容差下减少了 4.42 倍。
  • 在由薄板耦合瞬态模态的不规则划分波导中，所提 TC 的收敛速度比 Robin TC 快 9 倍。
  • 对于由 METIS 划分为 110 个子域的复杂 X 波段带通滤波器，所提 TC 在 180 次迭代后达到 $10^{-10}$ 容差，而 Robin TC 需要 1,020 次（接近六倍的减少）。
• 界面连续性：
  • 该方法大幅降低了界面跳跃。对于 METIS 划分的波导，最差情况下的法向通量跳跃 $|\Delta(\hat{n} \cdot B)|$ 从 $9.1 \times 10^{-7}$（Robin）降低至 $2.2 \times 10^{-12}$（所提方法），降幅接近六个数量级。
  • 最差情况下的切向场跳跃 $|\Delta(\hat{n} \times E)|$ 从 $1.4 \times 10^{-8}$ 降低至 $1.2 \times 10^{-10}$（超过两个数量级）。
• 准确性： 结果保持了高精度，所提方法与 Robin TC 解之间的相对 $L_2$ 范数差异极小（均匀情况下的 $3.2 \times 10^{-7}$）。

意义
论文声称，所提的双通道 FETI-DP-4λ 方法成功解决了大规模电磁领域分解问题中长期存在的收敛停滞挑战。通过直接利用基础麦克斯韦方程组，通过双通道方法强制执行切向场和法向通量的连续性，该方法实现了一种鲁棒、精确且大规模并行的迭代求解器。它证明了在不牺牲 DDM 的核心优势（并行性和可扩展性）或诉诸直接解法的前提下，可以在多样化的问题设置中（包括具有深瞬态波和复杂几何结构的场景）实现快速收敛。结果表明，对于传统 Robin TC 失效或 SOTC 不稳定的挑战性问题，该方法具有显著优势。