[論文レビュー] Convergence analysis and a novel Lagrange multiplier partitioned method for fluid-poroelastic interaction

Amy de Castro, Hyesuk Lee|arXiv (Cornell University)|Jan 20, 2026

Advanced Numerical Methods in Computational Mathematics被引用数 0

ひとこと要約

要約：本論文は、界面条件を課すために3つのラグランジュ乗数を用いた强結合・非反復的な分割法を用いた流体-孔弾性相互作用を提案し、シュール補完ソルバーと前処理を用いた完全離散系の空間収束を証明する。

ABSTRACT

We propose a partitioned method for the monolithic formulation of the Stokes-Biot system that incorporates Lagrange multipliers enforcing the interface conditions. The monolithic system is discretized using finite elements, and we establish convergence of the resulting approximation. A Schur complement based algorithm is developed together with an efficient preconditioner, enabling the fluid and poroelastic structure subproblems to be decoupled and solved independently at each time step. The Lagrange multipliers approximate the interface fluxes and act as Neumann boundary conditions for the subproblems, yielding parallel solution of the Stokes and Biot equations. Numerical experiments demonstrate the effectiveness of the proposed algorithm and validate the theoretical error estimate.

研究の動機と目的

• 接着されたStokes-Biot FPSI系をラグランジュ乗数を用いて界面条件を課す分割アプローチを開発する。
• 完全離散スキームの良定性と空間収束を確立する。
• 分解可能な副問題を並列計算可能とする効率的なシュール補完ベースのアルゴリズムを導出する。
• 結果として得られるシュール系を効率的に解くための堅牢な前処理を提供する。

提案手法

• 完全動的な二場Biot方程式をStokes方程式と結合させたStokes-Biot FPSI問題を定式化する。
• 質量保存、界面応力つり合い、Beavers-Joseph-Saffman接線条件を界面に課すために3つのラグランジュ乗数を導入する。
• Backward Eulerで時間離散化し、適合有限要素で空間離散化して、完全な離散鞍点問題を得る。
• 界面量のシュール補完方程式を導出し、各時間ステップで副問題を独立に解く効率的なアルゴリズムを開発する。
• シュール系の反復解法の収束を加速する前処理を提案する。
• 離散解の有限要素空間収束を証明し、数値実験で検証する。

Figure 1 : Manufactured solution: relative residual versus iteration number for BiCGstab(l).

実験結果

リサーチクエスチョン

• RQ1 副問題間のサブ反復を行わずに、完全動的なStokes-Biot FPSI系を効率的に解くことができるか。
• RQ2 界面条件を課すラグランジュ乗数は、安定かつ収束するモノリシック離散化を与えるか。
• RQ3 シュール補完に基づく分割アルゴリズムは、精度を保ちながら流体と孔弾性副問題を分離できるか。
• RQ4 各時間ステップでのシュール系の解法を効率化する前処理戦略は何か。
• RQ5 提案手法の数値的な収束性と頑健性はどうか。

主な発見

• 分割シュール補完法は流体と孔弾性副問題を分離し、各時間ステップで独立に解くことを可能にする。
• ラグランジュ乗数は界面フラックスを近似し、サブ問題のNeumannデータとして機能し、シュール系を介して結合を維持する。
• 提案FE空間と射影を用いた完全離散FPSI定式化に対して空間離散の収束性を確立する。
• シュール補完方程式の効率的な前処理を開発し、解法の性能を向上させる。
• 数値実験は理論的誤差推定を裏付け、アルゴリズムの有効性を示す。

(a) Velocities: $\bm{u}$ and $\frac{\partial\bm{\eta}}{\partial t}-\kappa\nabla p_{p}$ (arrows); vertical components (color)