[論文レビュー] Natural domain decomposition algorithms for the solution of time-harmonic elastic waves

Romain Brunet, Victorita Dolean+1 more|arXiv (Cornell University)|Apr 27, 2019

Advanced Numerical Methods in Computational Mathematics参考文献 27被引用数 7

ひとこと要約

本稿は、ナヴィエ方程式に従う時間調和弾性波問題を解くために最適化されたシュヴァルツ領域分割法を導入し、解析している。古典的シュヴァルツ法がこれらの方程式に対して発散することを証明しているが、物理的に妥当な伝達条件（特に局所的低周波数近似）を用いることで収束が達成され、Krylov法ソルバーの前処理として優れた性能を示す。これは、2次元および多数の部分領域を含む問題、特に不均質な伝達問題において、古典的手法を著しく上回る。

ABSTRACT

We study for the first time Schwarz domain decomposition methods for the solution of the Navier equations modeling the propagation of elastic waves. These equations in the time harmonic regime are difficult to solve by iterative methods, even more so than the Helmholtz equation. We first prove that the classical Schwarz method is not convergent when applied to the Navier equations, and can thus not be used as an iterative solver, only as a preconditioner for a Krylov method. We then introduce more natural transmission conditions between the subdomains, and show that if the overlap is not too small, this new Schwarz method is convergent. We illustrate our results with numerical experiments, both for situations covered by our technical two subdomain analysis, and situations that go far beyond, including many subdomains, cross points, heterogeneous materials in a transmission problem, and Krylov acceleration. Our numerical results show that the Schwarz method with adapted transmission conditions leads systematically to a better solver for the Navier equations than the classical Schwarz method.

研究の動機と目的

• 時間調和弾性波問題（ナヴィエ方程式で記述）に対する古典的シュヴァルツ法の収束特性を調査すること。
• 弾性波問題はスカラー波（ヘルムホルツ問題）よりもベクトル的性質と複数の波型を有するため、より困難であるが、そのような問題に対する効果的な領域分割ソルバーの欠如に取り組むこと。
• ナヴィエ方程式に対するシュヴァルツ法の収束を保証する新たな伝達条件の開発および解析を行うこと。
• 最適伝達条件の局所的低周波数近似を設計し、実装可能性を確保しながら収束性を維持すること。
• 多数の部分領域および不均質材料を含む複雑な設定において、新手法がソルバーや前処理として古典的シュヴァルツ法を上回る性能を示すことを数値的に検証すること。

提案手法

• 時間調和領域における2つの重複部分領域に対して、連続的レベルでの古典的シュヴァルツアルゴリズムを解析し、ナヴィエ方程式に対して収束しないことを証明する。
• より自然で物理的に妥当な伝達条件を部分領域間で導入し、それが有限ステップで収束を実現できることを示す。
• 最適伝達条件の局所的低周波数近似を提案し、非局所的作用素を回避することで実装可能性を高める。
• 部分領域の重なりが小さすぎない限り、新規最適化シュヴァルツ法の収束を証明する。
• 離散問題に対して、制限付き加法的シュヴァルツ（RAS）および最適化RAS（ORAS）を用いて実装し、並列計算を可能にする。
• GMRESをKrylov法ソルバーとして用い、古典的および最適化シュヴァルツ法を前処理として用いて性能を評価する。

実験結果

リサーチクエスチョン

• RQ1 なぜ古典的シュヴァルツ法は、時間調和弾性波のナヴィエ方程式に適用しても収束しないのか？
• RQ2 P波およびS波の物理的挙動を反映する伝達条件は、ナヴィエ方程式に対するシュヴァルツ法の収束を実現できるか？
• RQ3 最適伝達条件の局所的低周波数近似を構築することは可能か？ その場合、ナヴィエ方程式に対する収束性は保たれるか？
• RQ4 複雑な多数の部分領域および不均質な材料を含む問題において、最適化シュヴァルツ法は古典的手法と比較して、ソルバーや前処理としてどれほど優れた性能を示すか？
• RQ5 部分領域の重なりを増加させることで、最適化シュヴァルツ法の収束性は向上するのか？ その条件は何か？

主な発見

• 古典的シュヴァルツ法は、重なりがある場合でさえナヴィエ方程式に対して発散し、反復的ソルバーとして使用できない。
• 非局所的作用素を含む最適伝達条件を用いることで、シュヴァルツ法は2ステップで収束するが、実装には不適切である。
• 最適伝達条件の局所的低周波数近似は、重なりが十分に大きい限り、最適化シュヴァルツ法の収束を保証する。
• 数値実験により、最適化シュヴァルツ法はソルバーとして収束するが、古典的法は発散する。重なりが増加しても同様の傾向が続く。
• GMRESの前処理として用いた場合、最適化シュヴァルツ法は古典的法と比較して、特に多数の部分領域、交差点、材料の不均質性を含む問題で著しく少ない反復回数で収束する。
• 最適化シュヴァルツ法の収束性は重なりの増加とともに向上するが、古典的法の性能は同様の条件下で劣化する。