[论文解读] Domain Decomposition with local impedance conditions for the Helmholtz equation with absorption

Ivan G. Graham, Euan A. Spence+1 more|arXiv (Cornell University)|Jun 10, 2018

Electromagnetic Scattering and Analysis参考文献 52被引用 5

一句话总结

本文针对带有吸收项的Helmholtz方程提出了一种单重叠加Schwarz预条件子，采用具有局部阻抗条件的重叠子域以及基于划分单位的延拓算子。数值实验表明其在GMRES方法中实现了与k无关的收敛性，理论分析进一步证明了当ε/k为常数或缓慢增长时，该方法对吸收型问题具有鲁棒性，从而为通过预条件子处理纯Helmholtz问题提供了有效支持。

ABSTRACT

We consider one-level additive Schwarz preconditioners for a family of Helmholtz problems with absorption and increasing wavenumber $k$. These problems are discretized using the Galerkin method with nodal conforming finite elements of any (fixed) order on meshes with diameter $h = h(k)$, chosen to maintain accuracy as $k$ increases. The action of the preconditioner requires solution of independent (parallel) subproblems (with impedance boundary conditions) on overlapping subdomains of diameter $H$ and overlap $δ\leq H$. The solutions of these subproblems are linked together using prolongation/restriction operators defined using a partition of unity. In numerical experiments (with $δ\sim H$) for a model interior impedance problem, we observe robust (i.e. $k-$independent) GMRES convergence as $k$ increases. This provides a highly-parallel, $k-$robust one-level domain decomposition method. We provide supporting theory by studying the preconditioner applied to a range of absorptive problems, $k^2\mapsto k^2+ \mathrm{i} \varepsilon$, with absorption parameter $\varepsilon$. Working in the Helmholtz ``energy'' inner product, and using the underlying theory of Helmholtz boundary-value problems, we prove a $k-$independent upper bound on the norm of the preconditioned matrix, valid for all $\vert \varepsilon\vert \lesssim k^2$. We also prove a strictly-positive lower bound on the distance of the field of values of the preconditioned matrix from the origin which holds when $\varepsilon/k$ is constant or growing arbitrarily slowly with $k$. These results imply robustness of the preconditioner for the corresponding absorptive problem as k increases and give theoretical support for the observed robustness of the preconditioner for the pure Helmholtz problem.

研究动机与目标

• 为高频Helmholtz问题的吸收型情形开发一种并行且与k无关的预条件子。
• 分析在重叠子域上采用局部阻抗条件的单重叠加Schwarz方法的收敛性质。
• 通过界定预条件矩阵的范数和数值域，为数值实验中观察到的与k无关的鲁棒性提供理论依据。
• 支持将吸收型Helmholtz问题作为纯Helmholtz问题有效预条件子的使用。

提出的方法

• 该方法采用直径为H、重叠量δ ≤ H的重叠子域，并在子域界面处施加阻抗边界条件。
• 子问题通过有限元离散化并行求解，采用节点型协调有限元并设置网格尺寸h = h(k)，以在k增大时保持精度。
• 通过划分单位定义延拓与限制算子，确保在叠加Schwarz框架中的一致性与稳定性。
• 在Helmholtz能量内积下分析预条件子，利用已知的带吸收项Helmholtz边值问题的理论。
• 在|ε| ≤ Ck²且ε/k有界或缓慢增长的条件下，推导出预条件矩阵范数和数值域的理论界。
• 数值实验测试了不同k、H、ε和子域尺寸下的收敛性，包括固定H和采用Dirichlet子域条件的情形。

实验结果

研究问题

• RQ1 基于局部阻抗条件的区域分解预条件子是否能实现高频波数Helmholtz问题的与k无关的GMRES收敛？
• RQ2 在吸收型Helmholtz设置下，预条件矩阵的范数和数值域可建立何种理论界？
• RQ3 当k增大时，预条件子中吸收参数ε的选择如何影响收敛鲁棒性？
• RQ4 与Dirichlet条件相比，子域上采用阻抗条件是否显著提升预条件子的效率？
• RQ5 当ε ~ k时，吸收型问题能否作为纯Helmholtz问题的鲁棒预条件子？该结论是否得到理论支持？

主要发现

• 数值实验表明，当H ~ k^(-α)且α ∈ [0, 0.4]时，GMRES实现与k无关的收敛性，表明对吸收型问题具有鲁棒性。
• 证明了在|ε| ≤ Ck²条件下，预条件矩阵范数存在与k无关的上界，适用于所研究的吸收型问题范围。
• 当ε/k为常数或增长极其缓慢时，建立了预条件矩阵数值域到原点距离的严格正下界，从而确保收敛鲁棒性。
• 即使在预条件子中移除吸收项（ε_prec = 0），预条件子仍保持有效，尽管其收敛鲁棒性略低于ε_prec = k的情形。
• 若在子域上使用Dirichlet条件而非阻抗条件，收敛性能显著变差，尤其在大k时，表明阻抗条件具有关键作用。
• 出人意料的是，GMRES中使用随机初始猜测比零初始猜测迭代次数更少，尽管未提供解释。