[論文レビュー] Cross-Points in Domain Decomposition Methods with a Finite Element Discretization

Martin J. Gander, Kévin Santugini-Repiquet|arXiv (Cornell University)|Apr 18, 2014

Advanced Numerical Methods in Computational Mathematics参考文献 15被引用数 16

ひとこと要約

本稿は、非重複領域分割法における有限要素法を用いた際の、2つ以上の部分領域が交差する交差点におけるネウマン境界条件の一貫性のない離散化の課題に取り組む。2つの手法を提案する：安定した収束を実現するための追加界面変数を導入する補助変数法、および、交差点を共有するすべての部分領域を含めて通信を行う完全通信法。両手法とも収束の証明と最適パラメータ選定を可能にし、質量のルーミングにより代数的に高次のヴェンツェル型伝達条件が得られることを示している。

ABSTRACT

Non-overlapping domain decomposition methods necessarily have to exchange Dirichlet and Neumann traces at interfaces in order to be able to converge to the underlying mono-domain solution. Well known such non-overlapping methods are the Dirichlet-Neumann method, the FETI and Neumann-Neumann methods, and optimized Schwarz methods. For all these methods, cross-points in the domain decomposition configuration where more than two subdomains meet do not pose any problem at the continuous level, but care must be taken when the methods are discretized. We show in this paper two possible approaches for the consistent discretization of Neumann conditions at cross-points in a Finite Element setting.

研究の動機と目的

• 有限要素法を用いた領域分割法における、複数の部分領域が交差する交差点におけるネウマン境界条件の離散化の不一致を解消すること。
• 複数の部分領域が交差する点におけるネウマントレースの数学的に一貫性があり収束性のある離散化戦略を提供すること。
• 交差点特異性が存在する状況下でも、有限要素離散化を用いた非重複最適化シュバルツ法の収束を証明すること。
• 質量のルーミングにより、接線微分の離散化を回避して代数的に高次のヴェンツェル型伝達条件を生成できることを示すこと。
• 交差点が存在する状況下で収束の振動を最小化する最適なロビンパラメータおよびオーバーランディング要因を同定すること。

提案手法

• 補助変数法を導入し、部分領域が追加の界面変数を保持することで、交差点におけるネウマンデータの一貫性のある離散化を可能にする。
• 完全通信法を採用し、交差点を共有するすべての部分領域（点で接続されるものでさえも）と情報交換を行うことで、一貫性を維持する。
• エネルギー推定を用いて、補助変数を用いた場合の領域分割アルゴリズムの収束を証明する。これは、交差点における標準的な有限要素離散化では達成できない結果である。
• 交差点における複数の可能なネウマントレース分割のうち、反復過程での振動を最小化する方法を提案する。
• ロビン条件の質量のルーミングが、接線微分の直接離散化を回避して代数的に高次のヴェンツェル型伝達条件と同等の結果をもたらすことを示す。
• 2×2の部分領域と$Q_1$有限要素を用いたラプラス方程式の数値実験を行い、さまざまなロビンパラメータおよびオーバーランディングレベルにおける収束因子を評価する。

実験結果

リサーチクエスチョン

• RQ1 有限要素法を用いた領域分割法における、交差点におけるネウマン条件の一貫性のある離散化はどのように実現できるか？
• RQ2 交差点が存在する状況下でも、有限要素離散化を用いた非重複最適化シュバルツ法の収束を証明できるか？
• RQ3 反復プロセスにおける振動を最小化するために、交差点におけるネウマントレースをどのように分割すべきか？
• RQ4 有限要素の文脈において、ロビン条件の質量のルーミングによって代数的に高次のヴェンツェル型伝達条件を生成できるか？
• RQ5 交差点が存在する状況下で、最適なロビンパラメータおよびオーバーランディング要因はメッシュの細分化にどのようにスケーリングされるか？

主な発見

• 補助変数法により、エネルギー推定を用いた厳密な収束証明が可能となり、これは交差点における標準的な有限要素離散化では達成できない。
• 完全通信法により、交差点を共有するすべての部分領域（点で接続されるものでさえも）を含めた通信が保証され、ネウマン条件の離散化の一貫性が確保される。
• $100\times100$メッシュ（1部分領域あたり）の場合、完全通信法を用いて$\kappa=0.6013464$の最良収束因子が$ p=2.0 $および$ \omega=122.5 $の条件下で達成された。
• 補助変数法を用いた場合、$100\times100$セルにおいて$ p=2.0 $および$ \omega=122.0 $の条件下で最良収束因子$ \kappa=0.6006753 $が得られた。
• 最適なロビンパラメータ$ p $は、交差点において$ O(1/h) $のスケーリングを示すことが確認され、従来の理論的推定と一致する。これは、通常の界面点における$ O(1/\sqrt{h}) $とは異なり、交差点特有の性質を反映している。
• ロビン条件の質量のルーミングにより、接線微分の離散化を回避して代数的に高次のヴェンツェル型伝達条件が得られ、交差点においても同様に有効である。