[论文解读] Recent Results on Domain Decomposition Preconditioning for the High-frequency Helmholtz Equation using Absorption

I. G. Graham, Euan A. Spence+1 more|arXiv (Cornell University)|Jun 23, 2016

Advanced Numerical Methods in Computational Mathematics被引用 5

一句话总结

本文提出了一种用于高频Helmholtz方程的域分解预条件方法，通过引入吸收项来稳定系统并改善收敛性。通过求解带有吸收项的位移Helmholtz问题，并应用具有阻抗或PML型子域求解的多水平域分解方法，该方法在串行机器上实现了迭代次数与$\mathcal{O}(n^{0.2})$成比例，求解时间与$\mathcal{O}(n^{\alpha})$成比例，其中$\alpha \in [1.3, 1.4]$，即使在$ h \sim k^{-3/2} $的污染自由区域也表现良好。

ABSTRACT

In this paper we present an overview of recent progress on the development and analysis of domain decomposition preconditioners for discretised Helmholtz problems, where the preconditioner is constructed from the corresponding problem with added absorption. Our preconditioners incorporate local subproblems that can have various boundary conditions, and include the possibility of a global coarse mesh. While the rigorous analysis describes preconditioners for the Helmholtz problem with added absorption, this theory also informs the development of efficient multilevel solvers for the "pure" Helmholtz problem without absorption. For this case, 2D experiments for problems containing up to about $50$ wavelengths are presented. The experiments show iteration counts of order about $\mathcal{O}(n^{0.2})$ and times (on a serial machine) of order about $\mathcal{O}(n^α)$, { with $α\in [1.3,1.4]$} for solving systems of dimension $n$. This holds both in the pollution-free case corresponding to meshes with grid size $\mathcal{O}(k^{-3/2})$ (as the wavenumber $k$ increases), and also for discretisations with a fixed number of grid points per wavelength, commonly used in applications. Parallelisation of the algorithms is also briefly discussed.

研究动机与目标

• 开发适用于高频Helmholtz问题的鲁棒、可扩展的迭代求解器，此类问题因污染效应和缺乏强制性而长期难以求解。
• 分析并设计基于添加吸收项的Helmholtz方程的域分解预条件器，确保问题的适定性和收敛性。
• 在$ h \sim k^{-3/2} $的污染自由区域，实现与波数$ k $无关的最优收敛速率。
• 将该方法扩展至可变波速和不连续介质，展示其在真实声学建模场景中的鲁棒性。
• 提供一个结合严格分析与数值实验的框架，用于调节吸收参数、粗网格和子域边界条件。

提出的方法

• 通过求解位移Helmholtz问题$-(\Delta + (k^2 + i\varepsilon))u = f$，其中$\varepsilon > 0$，构建预条件器，以确保强制性和适定性。
• 采用域分解方法，使局部子问题可使用狄利克雷、阻抗或PML型边界条件，以改善条件数。
• 实现两水平或多层次预条件器，引入全局粗网格以捕捉低频模态并提升收敛性。
• 采用内外迭代方法，其中粗网格求解由使用阻抗子域求解的一水平方法进行预条件化。
• 调节吸收参数$\varepsilon$，以在逼近原问题（$\varepsilon \to 0$）与稳定位移问题（$\varepsilon$较大）之间取得平衡。
• 使用特征值分布分析，严格界定预条件系统在GMRES方法中的收敛速率。

实验结果

研究问题

• RQ1 何种吸收参数$\varepsilon$的范围可确保位移问题作为原Helmholtz问题的有效预条件器？
• RQ2 何种$\varepsilon$的范围可确保在位移问题上应用的域分解方法本身能快速收敛？
• RQ3 子域边界条件的选择（如狄利克雷与阻抗）如何影响预条件器的收敛性？
• RQ4 该方法是否能在高频问题中实现$\mathcal{O}(n^{\alpha})$的求解时间，且$\alpha < 1.5$，即使在污染自由区域也成立？
• RQ5 该方法在具有不连续波速的问题上表现如何？未解析或解析不连续性的影响是什么？

主要发现

• 对于$ h \sim k^{-3/2} $的问题，该方法在串行机器上实现了迭代次数与$\mathcal{O}(n^{0.2})$成比例，求解时间与$\mathcal{O}(n^{\alpha})$成比例，其中$\alpha \in [1.3, 1.4]$，即使在污染自由区域也表现良好。
• 对于固定每波长网格点数，采用阻抗局部求解的高并行性一水平方法实现了相似的时间复杂度，支持高效并行化。
• 采用较细粗网格和使用阻抗子域求解的内外预条件的两水平方法在总求解时间方面优于其他配置。
• 该方法对可变波速具有鲁棒性：在$ c^* = 1.5 $、$ c^* = 1.0 $和$ c^* = 0.66 $的实验中，迭代次数和时间复杂度在不连续介质中保持一致。
• 当波速不连续性未被解析时，该方法仍表现良好，仅导致迭代次数和求解时间略有增加，与解析情况相比增幅较小。
• 初步的并行实验表明一水平方法具有可扩展性，显示出在大规模应用中的巨大潜力。