Discrete FEM-BEM coupling with the Generalized Optimized Schwarz Method

技术摘要：基于广义优化 Schwarz 方法的离散 FEM-BEM 耦合

问题陈述
本文研究了对由亥姆霍兹方程（Helmholtz equation）描述的在有界和无界区域内声波传播问题的数值解。目标是开发一种通用且高效的求解器，能够处理非均匀材料和复杂的几何形状，包括具有交叉点（cross-points，即由两个以上子域共享的点）的几何结构。作者致力于将有限元法（FEM，适用于非均匀有界区域）与边界元法（BEM，适用于分段均匀及无界区域）进行耦合。一个特别需要解决的挑战是，经典的 FEM-BEM 形式在“伪共振”（spurious resonances，即变分形式存在非平凡核，但物理问题本身是适定的波数）处可能出现病态问题。

方法论
作者提出了一种基于广义优化 Schwarz 方法（Generalized Optimized Schwarz Method, GOSM）的非重叠子域分解（Domain Decomposition, DD）方法。该方法涉及：

1. 子结构化（Substructuring）： 将定义域划分为若干个非重叠的子域。迭代求解器的未知量与界面（骨架，skeleton）相关，而非体量（volume）。
2. 离散框架： 该分析通过使用 Galerkin 近似（FEM 和 BEM）将 [20] 中提出的连续理论扩展到全离散设置。
3. 算子形式化： 该方法将全局问题重新表述为如下形式的子结构化方程 $(Id + \Pi S)q = \text{rhs}$，其中：
   • $S$ 是一个散射算子，用于求解带有优化阻抗边界条件的局部问题。
   • $\Pi$ 是一个非局部交换算子（用于替代经典优化 Schwarz 方法中使用的标准交换算子），旨在处理交叉点。
   • $q$ 代表阻抗迹（impedance traces）在骨架上的分布。
4. FEM-BEM 耦合： 本研究专门分析了三种经典的 FEM-BEM 耦合方式：Costabel、Johnson-Nédélec 和 Bielak-MacCamy。研究调查了与这些耦合相关的边界算子 $A_\Gamma$ 的“符号假设”（sign assumption，与物理耗散相关）。
5. 迭代求解器： 当符号假设成立时（特别是对于 Costabel 耦合），算子 $\Pi S$ 是一个收缩映射。这允许使用具有保证几何收敛性的 Richardson 迭代方案。对于符号假设不成立的耦合（Johnson-Nédélec、Bielak-MacCamy），则提出了 GMRes 求解器。
6. 传输算子（Transmission Operators）： 数值实验使用了多种传输算子，包括局部 Després 算子、非局部 Yukawa 超奇异算子，以及基于 Schur 补（Schur complement）的、由正定 Dirichlet-to-Neumann 映射导出的算子。

核心贡献

• 离散适定性： 本文确立了经典 FEM-BEM 耦合的子结构化形式在不同于经典伪共振的波数下是适定的。
• 核维度关系： 建立了一个关于初始变分形式、局部问题和子结构化形式的核维度的显式关系（定理 7.3）。
• 对伪共振的鲁棒性： 一个重要的理论结果是，证明了 Costabel FEM-BEM 耦合的子结构化形式即使在伪共振处仍然是适定的。这是因为子结构化算子的核维度与局部问题的核维度一致，使得在特定条件下子结构化系统具有平凡核。
• 收敛理论： 作者证明了对于满足符号假设的 Costabel 耦合，子结构化算子具有“单位阵 + 收缩映射”的形式。这保证了 Richardson 迭代的几何收敛性。
• 交叉点处理： 该方法被证明对交叉点具有鲁棒性，无需在每个交叉点增加额外的未知量或耦合这些点处的局部问题，这是其他 DD 方法中的常见难点。
• $h$-鲁棒性： 理论估计和数值实验表明，只要使用非局部传输算子（如 Yukawa 或 Schur 补类型），即可实现对网格尺寸 $h$ 具有鲁棒性的几何收敛（即 $h$-鲁棒）。

结果
数值实验在 2D 环境下针对均匀和非均匀问题，以及涉及不可渗透障碍物和开放腔体（引入交叉点）的问题进行了测试。

• 迭代次数： 针对网格尺寸（$h$）和波数（$\kappa$）对收敛所需的迭代次数进行了分析。
  • 使用局部传输算子（Després）的配置，其迭代次数随网格加密呈现准二次增长（$O(h^{1.69} \dots O(h^{1.78})$）。
  • 使用非局部传输算子（Yukawa、Schur 补）的配置表现出随网格加密几乎不变的迭代次数，证实了 $h$-鲁棒性。
  • 关于波数，非局部配置通常表现出准常数或亚线性增长，而局部配置则表现出亚线性增长。
• 耦合对比： 虽然针对 Costabel 耦合推导了几何收敛性的理论保证，但数值结果表明，在测试场景中，Johnson-Nédélec 和 Bielak-MacCamy 耦合产生了相似的迭代次数。
• 交叉点鲁棒性： 在涉及交叉点的矩形腔体测试案例中，该方法保持了 $h$-鲁棒性，迭代次数在网格加密过程中保持稳定。
• 传输算子选择： “Y-S”配置（针对 BEM 使用 Yukawa 算子，针对 FEM 使用 Schur 补）通常产生的迭代次数最低。研究强调，将传输算子的类型与子域的性质（稀疏型 vs 稠密型）相匹配，有利于提高实现效率。

意义与主张
本文声称提供了首个具有保证收敛界的 $h$-鲁棒迭代算法用于 FEM-BEM 耦合。通过将经典的 FEM-BEM 耦合重新表述为子结构化形式，作者证明了即使在存在交叉点和伪共振的情况下，也可以使用领域分解技术高效地求解这些经典方法。这项工作弥合了连续领域分解理论与离散 FEM-BEM 实现之间的鸿沟，为耦合这些方法提供了一种系统策略，该策略在理论上是严谨的，在数值上也是高效的，适用于高性能计算应用。作者强调，我们的方法避免了在交叉点处耦合局部问题的复杂性，这是与其他多域 FEM-BEM 方案相比的一个显著优势。