[論文レビュー] A finite element method with strong mass conservation for Biot's linear consolidation model

Béatrice Rivière, Guido Kanschat|arXiv (Cornell University)|Dec 20, 2017

Advanced Numerical Methods in Computational Mathematics参考文献 10被引用数 46

ひとこと要約

本稿では、混合有限要素法と内部ペナルティ不連続ガラーキン項を組み合わせることで、強い質量保存を保証するBiotの線形圧縮モデルに対するH(div)-適合有限要素法を提案する。このスキームは、流体速度および変位の両方で最適なL²収束率を達成し、流体および固体の体積変化率が変化しても安定した性能を示し、メッシュの細分化パラメータに依存せずに正確な質量バランスを維持する。

ABSTRACT

An H(div) conforming finite element method for solving the linear Biot equations is analyzed. Formulations for the standard mixed method are combined with formulation of interior penalty discontinuous Galerkin method to obtain a consistent scheme. Optimal convergence rates are obtained.

研究の動機と目的

• ポーラス弾性解析における強い質量保存を保証する有限要素法の開発。
• 従来のTaylor-Hood混合法が抱える、流体速度と変位の収束率の差を解消すること。
• 非圧縮性流体（cₛ = 0）を含む一般の材料パラメータ下でも、速度および変位の両方のL²収束率を最適化すること。
• Biot-Willis定数αおよび貯留係数cₛに対してロバストであること、特に非圧縮性流体の極限（cₛ = 0）においても同様。
• 射影を用いない一貫性のある流体および固体速度の結合を実現し、パラメータ領域全体にわたる安定性と精度を確保すること。

提案手法

• 流体速度および固体変位にH(div)-適合Raviart-Thomas有限要素を用い、局所的質量保存を保証する。
• 流体圧力と速度の標準的混合式と、内部ペナルティ不連続ガラーキン項を組み合わせ、変位場のH¹連続性を強制する。
• 速度および変位の両方に同一のベクトル有限要素空間を用いることで、L²内での等価次数近似を可能にする。
• 最適な誤差推定を導出するために、変位場の発散に関するスーパーアプロキマション仮定を用いる。
• 時間積分にはθ=0.501を用いたθスキームを採用し、強いA安定性を確保するとともに、時間積分誤差を低減する。
• 離散的質量バランス残差Δm(t)を用いて質量保存を検証し、粗いメッシュ上でも機械精度（≈10⁻¹⁷）以内に保たれていることが示された。

実験結果

リサーチクエスチョン

• RQ1 射影に基づく結合に依存せずに、Biotの圧縮モデルにおける強い質量保存を保証する有限要素法を構築可能か？
• RQ2 H(div)-適合混合法と内部ペナルティDG項を組み合わせることで、速度および変位の両方のフィールドに最適な収束率が得られるか？
• RQ3 Biot-Willis定数αおよび貯留係数cₛに対して、特に非圧縮性流体の極限（cₛ = 0）においてもロバストか？
• RQ4 L²内での等価次数近似を維持しつつ、安定性と一貫性を保ちながら、最適な収束率を達成できるか？
• RQ5 メッシュサイズおよび時間ステップに依存せずに、離散レベルで正確な質量バランスが維持されるか？

主な発見

• RT₁/Q₁要素を用いた場合、数値実験により流体速度および変位の両方で最適なL²収束率（h²）を達成した。
• RT₂/Q₂要素を用いた場合、L²収束率はh³に達し、高次近似における最適収束性を示した。
• 変位場の発散誤差はh¹の収束率を示し、理論的仮定および数値的検証と整合的であった。
• 離散的質量バランス残差Δm(t)は、cₛ = 0およびα ≠ 1を含むすべてのテストパラメータ条件下でも機械精度（≈10⁻¹⁷）以内に保たれた。
• 従来の方法が性能を低下または最適収束を達成できなかったcₛ = 0（非圧縮性流体）の場合でも、本手法は安定した性能を示した。
• 数値的検証により、本スキームはさまざまな材料パラメータ下で安定かつ正確であり、メッシュの細分化に伴い期待通りの速度で誤差が減少することが確認された。