Convergence of parallel overlapping domain decomposition methods with impedance boundary conditions for time-harmonic Maxwell equations in heterogeneous media

技术摘要：用于时谐麦克斯韦方程组的并行重叠区域分解方法的收敛性

问题陈述
本文研究了非均匀介质中时谐麦克斯韦方程组的数值解问题，该问题在雷达、无线通信和医学成像等应用中至关重要。在高频情况下，这些问题面临显著挑战：为了解析振荡，需要精细的网格，这会导致大规模、非埃尔米特（non-Hermitian）且高度不定（indefinite）的线性系统。虽然两层区域分解（DD）方法已经建立了收敛理论（通常依赖于能够解析波传播的粗空间），但对于麦克斯韦方程组的一层重叠 DD 方法，尤其是在非均匀设置下的理论研究仍然缺乏严谨的收敛理论。作者指出，用于两层方法的标准值域（field-of-values）论证对于本文所讨论的一层方法并不适用，因为其值域包含了原点。

方法论
作者分析了一种在子域界面上采用阻抗边界条件的并行重叠 Schwarz 方法。该方法论分为两个主要阶段：连续层分析和离散层分析。

1. 连续层分析：

   • 适定性（Well-posedness）： 本文证明了并行 Schwarz 迭代在适当函数空间（$X_{imp}$ 和 $U$）中对于一般 Lipschitz 子域的适定性。这依赖于关于麦克斯韦解的调和分析结果。解决的一个关键技术障碍是单位分解函数（partition of unity functions）的正则性及其产生的误差项，这需要关于条带分解几何结构的特定假设（假设 2.11）。
   • 误差传播： 误差传播算子通过“阻抗到阻抗映射”（impedance-to-impedance maps）进行表征。这些映射描述了阻抗数据（切向旋度和场本身）如何在相邻界面之间传播。
   • 条带分解： 对于条带域分解，作者推导出了显式的收敛估计。他们表明，$N$ 次误差传播算子的收敛性取决于特定阻抗到阻记映射（记作 $\rho$ 和 $\gamma$）的范数，这些映射分别代表“左到右”以及“左到左”（或“右到右”）的相互作用。

2. 离散层分析：

   • 离散化： 连续方法使用 Nédélec 有限元进行离散化。离散并行 Schwarz 方法被识别为预处理 Richardson 迭代的限制加性 Schwarz（RAS-imp）方法。
   • 离散误差传播： 作者利用离散阻抗到阻抗映射构建了离散误差传播算子。
   • 收敛继承： 在假设离散阻抗到阻抗映射随网格尺寸 $h \to 0$ 收敛于其连续对应物的情况下，本文证明了离散方法继承了连续方法的收敛行为。

核心贡献

• PDE 层级框架： 开发了一个适用于变系数麦克斯韦方程组的并行重叠 Schwarz 方法与阻抗传输条件的框架。这包括证明了对于一般 Lipschitz 子域的适定性，并推导出了通过阻抗到阻抗映射表达的误差方程。
• 条带收敛估计： 推导了条带分解下的收敛充分条件。其分析将限制误差传播算子范数的任务简化为对界面映射（$\rho$ 和 $\gamma$）范数的估计。
• 离散收敛理论： 证明了 Nédélec 离散化（RAS-imp）在离散阻抗到阻抗映射近似连续映射的前提下，能够继承连续方法的收敛特性。
• 数值验证： 通过广泛的数值实验验证了理论。这些实验包括：
  • 验证离散误差传播算子的幂收缩性（power contractivity）。
  • 计算阻抗到阻抗映射的范数，确认其是有界的且与理论预测一致。
  • 在均匀和非均匀介质中的条带分解和棋盘格分解上的性能测试，显示出 RAS-imp 既可以作为有效的平稳迭代法，也可以作为 Krylov 预条件子（用于 GMRES）。

结果

• 理论界限： 本文确立了对于条带分解，如果阻抗到阻抗映射的范数足够小，则误差传播算子 $T^N$ 具有收缩性。具体而言，$N$ 次算子的范数受限于 $\rho$ 和 $\gamma$ 的函数。
• 数值观察：
  • 在条带分解中，残差历史呈现出“阶梯式衰减”，对应于幂收敛性，即误差范数在大约 $N$ 次迭代后（其中 $N$ 是子域数量）显著下降。
  • 数值计算了阻抗到阻抗映射的范数（$\rho$ 和 $\gamma$）。发现 $\rho$（代表跨越整个区域的传输）远低于 1，而 $\gamma$（代表同向的反射/相互作用）接近 1 但保持有界。
  • 该方法对于波数 $k$ 表现出鲁棒性（对于固定的 $N$），并且对子域数量具有可扩展性。
  • 对于棋盘格分解（超出了条带分析的严格理论范围），该方法仍然表现出幂收敛性，尽管收缩阈值移动到了大约 $2\sqrt{N}$ 次迭代。
  • 在非均匀介质中，由于空间变化的介电常数引起内部反射，收敛速度通常比均匀情况慢，但该方法仍然有效。

意义与主张
本文声称提供了关于时谐麦克斯韦方程组在一层重叠区域分解方法在连续层和离散层层面的首个收敛理论。作者强调，这项工作填补了现有两层理论留下的空白，因为两层理论往往依赖于会掩盖一层方法中关键传播行为的粗空间。

作者对其研究范围相对于类似工作（引用自 [20, 21, 27] 关于 Helmholtz 方程的研究）持谦虚态度。他们明确指出了三个局限性：

1. 连续迭代的适定性结果需要一个正则性假设（假设 2.11），该假设在麦克斯韦方程组的通用 Lipschitz 子域中尚未得到证明，这与 Helmholtz 方程的情况不同。
2. 离散阻抗到阻抗映射向连续映射收敛的假设（假设 6.5）是假设而非证明，这是由于麦克斯韦解的正则性较低以及必要的对偶论证的复杂性。
3. 本文未能在连续层级提供阻抗到阻抗映射的高频界限，因为所需的半经典缺陷测度理论（semiclassical defect measure theory）尚未在时谐麦克斯韦方程组中建立。

尽管存在这些局限性，本文仍断言所开发的框架成功地通过界面映射表征了误差传播，且数值实验有力地支持了理论预测，验证了 RAS-imp 是解决这些挑战性电磁问题的鲁棒求解器。