[論文レビュー] Simulation of Deformation and Flow in Fractured, Poroelastic Materials

Katja K. Hanowski, Oliver Sander|arXiv (Cornell University)|Jun 18, 2016

Advanced Numerical Methods in Computational Mathematics参考文献 3被引用数 8

ひとこと要約

本稿では、単一の有限長の亀裂を有する破壊的多孔質媒体における連成流体流れと多孔質弾性変形をシミュレートするための新しい変分形式と数値法を提示する。マトリックスの力学にはXFEMを、亀裂内の流れには低次元のダルシー模型を組み合わせ、亀裂先端における亀裂幅の消失にもかかわらず重み付きソボレフ空間内での解の存在を証明し、2次元および3次元のシミュレーションで最適収束率を示し、高速なサブストラクチャリングソルバー収束を実現する。

ABSTRACT

We introduce a coupled system of PDEs for the modeling of the fluid-fluid and fluid-solid interaction in a poroelastic material with a single static fracture. The fluid flow in the fracture is modeled by a lower-dimensional Darcy equation, which interacts with the surrounding rock matrix and the fluid it contains. We explicitly allow the fracture to end within the domain, and the fracture width is an unknown of the problem. The resulting weak problem is nonlinear, elliptic and symmetric, and can be given the structure of a fixed-point problem. We show that the coupled fluid-fluid problem has a solution in a specially crafted Sobolev space, even though the fracture width cannot be bounded away from zero near the crack tip. For numerical simulations, we combine XFEM discretizations for the rock matrix deformation and pore pressure with a standard lower-dimensional finite element method for the fracture flow problem. The resulting coupled discrete system consists of linear subdomain problems coupled by nonlinear coupling conditions. We solve the coupled system with a substructuring solver and observe very fast convergence. We also observe optimal mesh dependence of the discretization errors even in the presence of crack tips.

研究の動機と目的

• マトリックスの周囲の多孔質弾性マトリックスにおける流体の流れと変形の非線形結合をモデル化すること。流体-流体および流体-固体の相互作用を含む。
• 領域内に終端する亀裂、特に亀裂先端における亀裂幅の消失により解の特異性が生じるという課題に対処すること。
• 亀裂幅がゼロから離れていない場合でも、特別に構築された重み付きソボレフ空間内での連成系の弱解の存在を確立すること。
• マトリックス変位および圧力にXFEMを、亀裂内流れに標準FEMを用いた、モノリシックで適合的な離散化法を開発し、特異性が存在しても最適収束率を保証すること。
• 2次元および3次元設定において、非線形連成系を効率的に解くサブス トラクチャリングソルバーの実装と検証を行うこと。

提案手法

• マトリックス内に3次元の多孔質弾性方程式と、亀裂内に1次元のダルシー方程式を組み合わせたPDE連成系を定式化し、亀裂幅を未知変数として扱う。
• 非線形結合を扱うために固定点反復フレームワークを採用し、亀裂の透水率が未知の亀裂幅に逆比例することを考慮する。
• 亀裂先端における応力および流体圧力の特異性を、メッシュの整列を必要としない特異関数（例：$ r^{1/2} $）によるXFEMの拡張を用いる。
• マトリックス問題を、変位および空隙圧力にXFEMを用いた適合有限要素で離散化し、亀裂部分は標準的な低次元FEMを用いる。
• 得られた連成系を、サブドメイン問題を別々に処理し、非線形結合条件を効率的に強制するサブス トラクチャリングソルバーで解く。
• 亀裂幅が亀裂先端でゼロに近づくという退化にもかかわらず、解の存在を厳密に証明するための重み付きソボレフ空間を構築する。

実験結果

リサーチクエスチョン

• RQ1 亀裂幅が亀裂先端でゼロに近づき、ゼロから離れていない場合でも、重み付きソボレフ空間内に弱解が存在するか？
• RQ2 亀裂先端における特異性が存在する中で、逆透水率に起因する非線形結合（亀裂内流れとマトリックス変形の間）を的確に離散化する方法は何か？
• RQ3 特異性が存在し、亀裂が領域内に終端する場合、数値スキームの収束特性はどのように振る舞うか？
• RQ4 モノリシックな連成問題の定式化から生じる非線形に結合された系を、サブス トラクチャリングソルバーが効率的に処理できるか？
• RQ5 適切な拡張関数を用いたXFEMを用いることで、亀裂先端付近においても変位、圧力、亀裂内流れの最適収束率が保持されるか？

主な発見

• 亀裂幅が亀裂先端でゼロに近づく場合でも、特別に構築された重み付きソボレフ空間内に、連成多孔質弾性-亀裂内流れ系の弱解が存在する。
• マトリックス圧力の$ L^2 $誤差は$ ext{O}(h^3) $で減少し、$ H^1 $誤差は$ ext{O}(h^{3/2}) $で減少する。これは、圧力勾配特異性の正確な再現によるスーパー収束の可能性がある。
• 亀裂内圧力の$ L^2 $誤差は$ ext{O}(h^2) $で減少し、予想を上回る性能を示す。これは、マトリックス圧力からの結合項に高められた正則性が寄与している可能性がある。
• 変位および体積圧力の$ L^2 $誤差と$ H^1 $誤差は、それぞれ$ ext{O}(h^2) $および$ ext{O}(h) $の最適収束率を示し、XFEM離散化の堅牢性を確認する。
• サブス トラクチャリングソルバーは、すべてのテストケースで非常に少ない反復回数で収束し、非線形に結合された系に対する高い効率性を示す。
• 3次元シミュレーションは、本手法が現実的な幾何形状への適用可能性を示しており、物理的に整合した結果を示している：流体の流出が亀裂の開口を引き起こし、亀裂先端でvon Mises応力のピークが観察される。