[論文レビュー] 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

ひとこと要約

本稿では、吸収を用いて安定化させ、収束を向上させる高周波数ヘルムホルツ方程式に対するドメイン分解プリコンディショナの手法を提示する。吸収を追加したシフトされたヘルムホルツ問題を解き、インピーダンスまたは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.

研究の動機と目的

• 汚染と強制性の欠如のため、とりわけ困難とされる高周波数ヘルムホルツ問題に対する、強固でスケーラブルな反復解法の開発。
• 追加された吸収に基づくドメイン分解プリコンディショナの分析と設計を行い、適切な定式化と収束性を保証すること。
• 波数 $k$ に依存しない最適な収束レートを達成すること、特に $h \sim k^{-3/2}$ の汚染なし領域でも同様に成立させること。
• 可変波動速度および不連続な媒質への拡張を行い、現実的な音響モデリングの状況でも強健性を示すこと。
• 吸収パラメータ、粗いグリッド、部分領域境界条件のチューニングを可能にする、厳密な解析と数値実験を統合したフレームワークの提供。

提案手法

• 吸収パrameter $\varepsilon > 0$ を追加したシフトされたヘルムホルツ問題 $-(\Delta + (k^2 + i\varepsilon))u = f$ を解くことでプリコンディショナを構築し、強制性と適切な定式化を保証する。
• 局所的な部分問題にディリクレ、インピーダンス、またはPML型の境界条件を適用することで、条件数の改善を図るドメイン分解を用いる。
• 低周波数モードを捉えるためにグローバルな粗いグリッドを用いた二層または多層プリコンディショナを実装する。
• 粗いグリッドの解法をインナーオーバーサイド反復法でプリコンディショニングし、インピーダンス部分領域解法を用いる。
• 元の問題に近づける（$\varepsilon \to 0$）ことと、シフトされた問題の安定化（$\varepsilon$ を大きくすること）のトレードオフを最適化するために、吸収パラメータ $\varepsilon$ をチューニングする。
• 場の値解析を用いて、プリコンディショニングされた系のGMRES収束レートを厳密に評価する。

実験結果

リサーチクエスチョン

• RQ1 どの範囲の吸収 $\varepsilon$ が、元のヘルムホルツ問題に対する有効なプリコンディショナとしてのシフト問題を保証するか？
• 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]$）する。これは汚染なし領域でも成立する。
• 1波長あたりのグリッド点数を固定した場合、インピーダンス局所解法を用いた高並列性のある1層法は、類似した時間計算量を達成し、効率的な並列化が可能である。
• 比較的細かい粗いグリッドと、インナーオーバーサイドプリコンディショニングを用いた2層法が、他の構成よりも総合的な解法時間において優れている。
• 可変波動速度に対しても強健性を示す。$c^* = 1.5$、$c^* = 1.0$、$c^* = 0.66$ の実験では、不連続媒質においても一貫した反復回数と時間計算量のスケーリングが得られた。
• 波動速度の不連続性が解像されていなくても、反復回数と解法時間はわずかに増加するにとどまり、解像済みのケースと比較して著しい劣化は見られない。
• 初期の並列実験により、1層法のスケーラビリティが示され、大規模応用への強い可能性が示唆された。