[論文レビュー] An Energy-Stable, Bound-Preserving and Locally Conservative Numerical Framework for Multicomponent Gas Flow in Poroelastic Media

Huangxin Chen, Yuxiang Chen+2 more|arXiv (Cornell University)|Feb 11, 2026

Advanced Numerical Methods in Computational Mathematics被引用数 0

ひとこと要約

論文は、多孔弾性媒体における多成分気体の流れを模擬するための堅牢でエネルギー安定かつ境界を守る局所保守的な数値フレームワークを開発し、安定化線形スキームと適応時間ステップ実装、および特殊な変数消去戦略を組み合わせている。

ABSTRACT

In this paper, we propose a robust and efficient numerical framework for simulating multicomponent gas flow in poroelastic media, with a focus on preserving fundamental thermodynamic principles and ensuring computational reliability. The model captures the complex nonlinear coupling between multicomponent transport and solid deformation, while addressing critical numerical challenges such as mass conservation, energy stability, and molar density boundedness. To achieve this, we develop a stabilized discretization approach that guarantees the preservation of the original energy dissipation law and ensures the boundedness of each gas component's molar density. Furthermore, the proposed method incorporates an adaptive time-stepping strategy that dynamically adjusts the time step size based on the system's dynamics, significantly enhancing computational efficiency without compromising stability or accuracy. For spatial discretization, a mixed finite element method combined with an upwind scheme is employed for the flow and transport equations to ensure local mass conservation, while a discontinuous Galerkin (DG) method is utilized for discretizing the momentum equation of poroelasticity to effectively overcome numerical locking phenomena. Numerical experiments are presented to demonstrate the performance, robustness, and applicability of the method in simulating multicomponent gas flow under various scenarios.

研究の動機と目的

• 熱力学的一貫性の下で多成分気体流の正確なシミュレーションを動機付ける。
• 数値スキームにおいてエネルギー散逸、成分別モル密度の境界、局所質量保存を保持する。
• 多成分輸送と固体変形を結ぶ頑健で効率的な数値フレームワークを開発する。
• 適切な離散化と適応時間ステップを通じて数値安定性と物理的忠実性を確保する。

提案手法

• 熱力学的一貫性を持つ generalized Maxwell–Stefan–Darcy モデルを導出し、輸送とBiot poroelasticityを結合させる。
• 元のエネルギー散逸則を保持し、モル密度の境界を課す安定化線形スキームを用いる。
• 安定性または精度を損なうことなく効率性を高めるための適応時間ステップ戦略を採用する。
• 局所的質量保存と前方風上安定性を確保するため、流れと輸送には混合有限要素法を用いる。
• 数値ロックを緩和するため、ポロエラスティック運動方程式を不連続ガレージャン（DG）法で離散化する。
• 境界を(primary variables)として (c1,...,c_{M-1},c) を用いる変数変換を導入し、M番目成分密度を除去して境界を課す。

実験結果

リサーチクエスチョン

• RQ1 多成分気体流をどのようにモデル化すれば熱力学原理と質量保存を満たすか。
• RQ2 安定化されエネルギー散逸性を持つ数値スキームは全成分のモル密度を境界内に保てるか。
• RQ3 適応時間ステップ戦略は結合輸送と機械を安定させつつ効率を向上させるか。
• RQ4 どの離散化選択（流れには混合FVEs、力学にはDGなど）が数値ロックを最小化し局所保存を確保するのに最適か。

主な発見

• 提案された半離散スキームは各時間ステップでエネルギー安定性を保証する離散エネルギー散逸不等式を満たす。
• 変数消去戦略（cM = c − sum_{i=1}^{M-1} c_i）は主変数を (c1,...,c_{M-1},c) にした境界保持型の定式化を生む。
• 安定化された化学ポテンシャル形式はスキーム下で離散エネルギーが非増加になることを保証する。
• 完全離散スキームは混合有限要素によって全成分の局所質量保存を達成し、固体力学のロックを回避するDG離散化により実現する。
• 適応時間ステップアルゴリズムは安定性と精度を維持しつつ効率的な時間積分を提供する。
• 数値実験は多孔弾性媒体における多成分ガス流の多様なシナリオに対して頑健性と適用性を示す。