[論文レビュー] Robust preconditioning of monolithically coupled multiphysics problems

Karl Erik Holter, Miroslav Kuchta+1 more|arXiv (Cornell University)|Jan 15, 2020

Advanced Numerical Methods in Computational Mathematics参考文献 53被引用数 6

ひとこと要約

本稿では、界面における分数階および重み付きラプラシアン作用素を用いた、モノリシックに結合されたマルチフィジックス問題に対するロバストでパrameterに頑健な予めんづけ枠組みを提示する。分数ソボレフ空間における作用素予めんづけを活用することで、Darcy-Stokes問題や流体-構造連成問題といった困難な問題に対してもメッシュに依存しない収束性を達成し、数値実験により最適収束率およびさまざまな材料パラメータとメッシュの細分化に対して有界な条件数が確認された。

ABSTRACT

In many applications, one wants to model physical systems consisting of two different physical processes in two different domains that are coupled across a common interface. A crucial challenge is then that the solutions of the two different domains often depend critically on the interaction at the interface and therefore the problem cannot be easily decoupled into its subproblems. Here, we present a framework for finding robust preconditioners for a fairly general class of such problems by exploiting operators representing fractional and weighted Laplacians at the interface. Furthermore, we show feasibility of the framework for two common multiphysics problems; namely the Darcy-Stokes problem and a fluid--structure interaction problem. Numerical experiments that demonstrate the effectiveness of the approach are included.

研究の動機と目的

• 連成マルチフィジックス問題のモノリシックな定式化において、界面で部分問題が強く結合されている場合に、その予めんづけ戦略を構築すること。
• 粘度、透過率、およびBeavers-Joseph-Saffman係数などの物理的パラメータの変動に対し、パラメータに頑健な予めんづけを保証すること。
• 重み付きおよび分数階ソボレフ空間を用いて、界面結合問題への作用素予めんづけ理論を拡張すること。
• Darcy-Stokes流れと流体-構造連成問題という2つの代表的問題に対して、フレームワークを数値的に検証すること。
• 予めんづけられた系が、メッシュの細分化において有界な条件数と最適収束率を維持することを示すこと。

提案手法

• 界面連続性を強制するためにラグランジュ乗数を用いたモノリシックな弱形式で連成マルチフィジックス問題を定式化する。
• 界面変数の正しいスケーリングを捉えるために、分数階および重み付きソボレフ空間におけるリーマン写像に基づく予めんづけを構築する。
• 界面における分数階ラプラシアン作用素を用いて、伝達条件（Beavers-Joseph-Saffman条件を含む）をモデル化する。
• 速度、圧力、およびラグランジュ乗数のノルムを、物理的スケーリングを反映させ、連成系の適切な定式化を保証するように導出する。
• MINRES反復解法内に予めんづけを適用し、メッシュ収束性およびパラメータに頑健な収束性を達成する。
• メッシュサイズおよび材料パラメータを変化させた2次元Darcy-Stokesおよび流体-構造連成問題における数値実験により、手法を検証する。

実験結果

リサーチクエスチョン

• RQ1 粘度や透過率などの物理的パラメータの変動に対しても、依然として頑健であるようなモノリシックな予めんづけ枠組みを構築できるか？
• RQ2 重み付きおよび分数階ソボレフ空間を用いることで、連成系における界面変数の正しいスケーリングを反映する予めんづけを設計できるか？
• RQ3 提案された予めんづけは、Darcy-Stokes問題や流体-構造連成問題において、メッシュの細分化に対して有界な条件数と最適収束率を維持するか？
• RQ4 解法時間と反復回数の観点から、ドメイン分割法と比較して、モノリシックな予めんづけの相対的性能は如何ほどか？
• RQ5 適切な仮定の下で、部分問題の適切な定式化が古典的に確立されていない問題に対しても、このフレームワークを適用できるか？

主な発見

• 提案された予めんづけは、メッシュの細分化にわたって有界な条件数（例：Stokes-Navier問題ではh = 2⁻⁸で12.72）と最適収束率を達成した。
• 数値実験により、モノリシックな予めんづけはドメイン分割法を上回る性能を示し、Darcy-Stokes問題では1反復あたり約47秒、Stokes-Navier問題では約116秒の解法時間であったのに対し、DD法では同程度の性能で60～150秒を要した。
• Beavers-Joseph係数の極端な値（k = 10⁻⁶から1）の変動に対しても収束が安定し、パラメータに頑健であることが確認された。
• MINRESソルバの収束はメッシュに依存せず、物理的パラメータに対しても頑健であり、例としてStokes-Navier問題ではh = 2⁻⁸で47回の反復が一貫して観察された。
• 文献[24]の分析から導出された理論的ノルムが、ラグランジュ乗数を含む形に成功裏に拡張され、ロバストな予めんづけが可能になった。
• Darcy-Stokes問題と流体-構造連成問題という2つのベンチマーク問題においてフレームワークが検証され、誤差の収束率が理論的推定値と一致した。