[論文レビュー] Parameter-robust convergence analysis of fixed-stress split iterative method for multiple-permeability poroelasticity systems

Qingguo Hong, Johannes Kraus+2 more|arXiv (Cornell University)|Dec 31, 2018

Advanced Numerical Methods in Computational Mathematics参考文献 23被引用数 5

ひとこと要約

本稿では、多孔透性ポアレラスティシティ（MPET）系に対する固定応力分割反復法を提案し、すべての物理的パラメータに依存しない収束率を有する線形収束を証明した。この手法は、全流体圧力の和を用いた安定化を組み込むことで、パラメータに頑健な収束を実現し、プリコンディショニング付きMinResソルバと比較して、広範なパラメータ変動において優れた効率性と頑健性を示した。

ABSTRACT

We consider flux-based multiple-porosity/multiple-permeability poroelasticity systems describing multiple-network flow and deformation in a poro-elastic medium, sometimes also referred to as MPET models. The focus of the paper is on the convergence analysis of the fixed-stress split iteration, a commonly used coupling technique for the flow and mechanics equations in poromechanics. We formulate the fixed-stress split method in the present context and prove its linear convergence. The contraction rate of this fixed-point iteration does not depend on any of the physical parameters appearing in the model. The theory is confirmed by numerical results which further demonstrate the advantage of the fixed-stress split scheme over a fully implicit method relying on norm-equivalent preconditioning.

研究の動機と目的

• 広範な物理的パラメータの変動に対しても効果を保つ、多孔透性ポアレラスティシティ（MPET）系のための頑健な反復的結合手法の開発。
• MPET文脈における固定応力分割法の線形収束を証明し、収束率がモデルパラメータに依存しないことの確立。
• 固定応力スキームが、ノルム同等のプリコンディショニングを施した完全陰解法よりも優れているという理論的・計算的証拠の構築。
• バイオトのモデルからより複雑なMPET系（複数の流体ネットワークを有する）への固定応力分割フレームワークの拡張。
• メッシュの細分化やパラメータ領域を変化させた状況下でも、理論的収束率の数値的検証の実施。

提案手法

• 固定応力分割法は、MPET系に対して固定点反復として定式化され、流れと力学の式が分離される。
• 全流体圧力の和を用いた安定化が導入され、頑健性を確保するとともに、パラメータに頑健な収束解析を可能にする。
• 収束解析は、フラックスベースのMPET系に関する先行研究の安定性結果に依拠し、エネルギーノルムを用いて収縮推定を導出する。
• 変分式を用いてMPET系の弱形式を導出し、有限要素離散化と反復的解法を可能にする。
• 変位と圧力には標準的な有限要素が使用され、フラックスはダルシーの法則により計算される。
• 数値実験では、固定応力スキームとプリコンディショニング付きMinResソルバを比較し、残差を10^8倍に減少させるまでの反復回数を測定した。

実験結果

リサーチクエスチョン

• RQ1 固定応力分割法は、複数の流体ネットワークを有するMPET系へ拡張可能であり、物理的パラメータに依存しない収束性が証明可能か？
• RQ2 固定応力反復におけるパラメータに頑健な収束を保証する安定化戦略は何か？
• RQ3 固定応力スキームの収束率は、ノルム同等のプリコンディショニングを施した完全陰解法と比べてどのように異なるか？
• RQ4 固定応力法は、透水率、剛性、貯留係数の極端な変動に対しても、効率性と頑健性を維持できるか？
• RQ5 理論的収縮率は、メッシュの細分化や異なるパラメータ領域において、数値的に検証可能か？

主な発見

• 固定応力分割法は、MPETモデルのすべての物理的パラメータ（透水率、剛性、貯留係数など）に依存しない収束率を有する線形収束を示した。
• 4ネットワークMPETモデルにおいて、固定応力スキームは、透水率や剛性が極端な値（λ·10^8 や K₃·10^10）であっても、残差を10^8倍に減少させるために常に10〜14回の反復で十分であった。
• プリコンディショニング付きMinResソルバは、特に透水率や剛性が広範囲に変動する場合に、最大で45回の反復を要し、これは収束性の低さを示している。
• λ（ベース値の10^8倍）、K₃（ベース値の10^10倍）、K（ベース値の10^2倍）といった極端な値を含むすべてのテストパラメータ範囲で、この手法は安定的かつ効率的であった。
• 数値結果は理論的収縮率推定を確認しており、すべてのパラメータ変動およびメッシュ細分化において反復回数が低値（10〜14）に安定した。
• すべてのテスト状況において、固定応力スキームはノルム同等のプリコンディショニングを施したMinRes法を上回り、優れた頑健性と効率性を示した。