[論文レビュー] On a Robin-Robin domain decomposition method with optimal convergence rate

Wenbin Chen, Xuejun Xu+1 more|arXiv (Cornell University)|Dec 19, 2011

Advanced Numerical Methods in Computational Mathematics参考文献 25被引用数 2

ひとこと要約

本論文は、2パラメータのロビン＝ロビン非重複領域分割法を提案することで、長年の未解決問題を解決した。この手法は、メッシュサイズ h に依存しない最適収束速度を達成する。従来の方法は収束速度が 1 − O(h¹/²) に制限されるが、本手法は h に依存しない収束速度の安定性を保証しており、理論的考察と数値実験の両方で確認された。

ABSTRACT

In this paper, we shall solve a long-standing open problem: Is it possible that the convergence rate of the Lions’ Robin-Robin nonoverlapping domain decomposition(DD) method is independent of the mesh size h? We shall design a two-parameter Robin-Robin domain decomposition method. It is shown that the new DD method is optimal, which means the convergence rate is independent of the mesh size h. The traditional Robin-Robin domain decomposition method converges at a rate of 1 − O(h1/2), even under the optimal parameter. Numerical implementation confirming our theoretical findings shall be given. AMS subject classifications. Primary 65N30; Secondary 65M60

研究の動機と目的

• リヨンのロビン＝ロビン非重複領域分割法の収束速度がメッシュサイズ h に依存しないかどうかという未解決問題を解消すること。
• 最適収束特性を達成する2パラメータのロビン＝ロビンDD法を設計すること。
• 理論的に収束速度が h に依存しないことを証明し、数値実験でその妥当性を検証すること。これにより、古典的手法の 1 − O(h¹/²) の限界を克服すること。

提案手法

• 非重複領域分割フレームワークにおけるロビン伝達条件に対して2パラメータ形式を導入する。
• 最適パラメータ選択下での新手法の収束挙動を導出し、分析する。
• メッシュサイズ h に依存しない収束速度の理論的境界を確立する。
• スペクトル解析とエネルギー推定を用いて、収束速度の最適性を証明する。
• モデル問題に対して数値的に手法を実装し、理論的予測の妥当性を検証する。
• 古典的手法と収束挙動を比較することで、性能向上を明確にする。

実験結果

リサーチクエスチョン

• RQ1 ロビン＝ロビン非重複領域分割法の収束速度をメッシュサイズ h に依存させない形にできるか？
• RQ2 ロビン伝達条件のどのパラメータ化が非重複DD法における最適収束をもたらすか？
• RQ3 古典的手法の収束速度 1 − O(h¹/²) は根本的な限界を示しているか？
• RQ4 2パラメータのロビン＝ロビン形式はメッシュに依存しない収束を達成できるか？
• RQ5 理論的最適性が数値実験によって確認されるか？

主な発見

• 提案された2パラメータのロビン＝ロビン領域分割法は、メッシュサイズ h に依存しない収束速度を達成する。
• この手法は最適性が証明されており、h が小さくなるに従って収束速度が劣化しない。
• 古典的手法の収束速度が、最適パラメータ選択下でも 1 − O(h¹/²) に制限されることを示した。
• 数値実験により、理論的予測であるメッシュに依存しない収束が確認された。
• 収束速度の h 依存性による劣化を排除することで、新手法は古典的手法を上回る性能を示した。
• 理論的解析により、新手法の収束速度は、h にかかわらず1から一様に離れていることが示された。