[论文解读] Analysis of parallel Schwarz algorithms for time-harmonic problems using block Toeplitz matrices. ETNA - Electronic Transactions on Numerical Analysis

Alexander Kyriakis, Jennifer Pestana+2 more|arXiv (Cornell University)|Nov 9, 2021

Matrix Theory and Algorithms参考文献 34被引用 5

一句话总结

本论文通过块托普利茨矩阵分析了一类带罗宾传输条件的单重并行 Schwarz 方法在时 harmonic Helmholtz 方程和 Maxwell 方程中的收敛性和可扩展性。结果表明，在吸收条件（σ > 0）下，该方法在连续层面对于条带型域分解实现了弱可扩展性，当子域数量 N 较大时，收敛因子与 N 无关——该结论对一维和二维问题均成立。

ABSTRACT

In this work we study the convergence properties of the one-level parallel Schwarz method with Robin transmission conditions applied to the one-dimensional and two-dimensional Helmholtz and Maxwell's equations. One-level methods are not scalable in general. However, it has recently been proven that when impedance transmission conditions are used in the case of the algorithm being applied to the equations with absorption, then, under certain assumptions, scalability can be achieved and no coarse space is required. We show here that this result is also true for the iterative version of the method at the continuous level for strip-wise decompositions into subdomains that are typically encountered when solving wave-guide problems. The convergence proof relies on the particular block Toeplitz structure of the global iteration matrix. Although non-Hermitian, we prove that its limiting spectrum has a near identical form to that of a Hermitian matrix of the same structure. We illustrate our results with numerical experiments.

研究动机与目标

• 研究带罗宾传输条件的单重并行 Schwarz 方法在时 harmonic 波问题中的可扩展性。
• 解决在波传播背景下，特别是连续层面，Schwarz 方法缺乏理论分析的问题。
• 为一维和二维 Helmholtz 与 Maxwell 方程中的条带型分解建立弱可扩展性（收敛率与子域数量 N 无关）。
• 证明当存在吸收时，在特定参数条件下，无需粗网格空间即可保持可扩展性。
• 将谱分析技术扩展至由迭代 Schwarz 格式产生的非厄米特块托普利茨矩阵。

提出的方法

• 将 Schwarz 迭代形式化为作用于界面数据的块托普利茨矩阵，利用连续层面的公式化方法。
• 分析当 N → ∞ 时全局迭代矩阵的极限谱，证明其尽管非厄米特，却具有近似厄米特结构。
• 利用傅里叶分析与谱理论，刻画在子域数量众多时的收敛行为。
• 将文献 [6] 中的迭代矩阵形式化方法应用于条带型分解，适配带有吸收的波方程。
• 证明在一维问题中，对任意 σ > 0 均收敛；在二维问题中，在 σ、δ 和 L 满足特定条件时收敛。
• 将理论结果与实际预条件化方法联系起来，将分析与一重优化受限加性 Schwarz（ORAS）方法关联。

实验结果

研究问题

• RQ1 带罗宾条件的单重 Schwarz 方法能否在无粗网格空间的情况下，对时 harmonic 波问题实现弱可扩展性？
• RQ2 在 Helmholtz 与 Maxwell 方程中，条带型分解的 Schwarz 迭代矩阵的极限谱行为如何？
• RQ3 吸收（σ > 0）如何影响连续公式化下的收敛性与可扩展性？
• RQ4 在何种条件下方法具有 k-鲁棒性，这与域分解参数的选择有何关系？
• RQ5 能否利用非厄米特块托普利茨矩阵的谱分析推导出收敛性与可扩展性结果？

主要发现

• 尽管非厄米特，Schwarz 迭代矩阵的极限谱结构与同种块托普利茨形式的厄米特矩阵几乎完全相同。
• 在一维情况下，对任意 σ > 0 均收敛，且当 N 足够大时收敛因子与 N 无关，表明实现了弱可扩展性。
• 在二维情况下，对倏逝模态（˜k > k）或当 σ、δ 和 L 足够大时，可扩展性成立，且当 σ ≥ k 时收敛性可保证。
• 数值实验表明，当 N 超过某一临界值后，GMRES 迭代次数趋于稳定，表明其与 N 和波数 k 均无关（当 σ = 1 时）。
• 当吸收参数 σ 与 k 成正比时，方法实现 k-鲁棒性，与文献中已有结果一致。
• 该理论框架支持设计鲁棒的预条件子（如 ORAS 方法），并可在连续与离散层面进行分析。