Hybridizable discontinuous Galerkin methods for poroelastic wave propagation with symmetric stress approximation

技術的概要：対称応力近似を伴う多孔弾性波動伝播のためのハイブリッド化不連続ガラーキン法

問題提起
本論文は、固体マトリクスと間隙流体の相互作用を通じて流体飽和多孔質媒体における機械的擾乱を記述する多孔弾性波動伝播の数値シミュレーションを取り扱っている。支配方程式は、元々ビオによって開発されたもので、線形運動量バランスと質量保存を含む連成偏微分方程式系を形成する。これらの方程式は、強い固体 - 流体摩擦と波動速度の大きな対比、特にビオの遅い P 波の散逸性に起因する高い剛性によって特徴づけられる。有限差分法、有限体積法、標準的な有限要素法を含む既存の数値手法は、波動解像度のために高次精度が要求される場合、計算コストに関する課題に直面している。不連続ガラーキン（DG）法は柔軟性を提供するが、自由度（DOFs）のグローバルな結合により、許容しがたい計算コストをもたらす。さらに、標準的な定式化は、弾性応力を主要変数として使用する際、ほぼ非圧縮性材料（ロック現象）への対応や、トラクション境界条件の自然的な課定において困難を伴うことが多い。

手法
著者らは、これらの限界を克服するために設計されたハイブリッド化不連続ガラーキン（HDG）法のファミリーを提案する。このアプローチには、いくつかの重要な方法論的ステップが含まれる：

1. 第一階対称定式化：支配方程式は、第一階対称双曲系として書き換えられる。決定的なことに、この定式化は弾性応力単独ではなく、全応力テンソル（$\sigma = C\varepsilon(u_s) - \alpha p I$）を使用する。この選択は、特にトラクション境界と圧力境界が交差する場所で、トラクション境界条件の自然的な課定を容易にする。
2. ハイブリッド化戦略：この手法は、2 つの異なる HDG アプローチを組み合わせる：
   • HDG+ 法（最近の線形弾性における進歩に触発されたもの）は、固体変位と応力成分に適用される。これにより、最適収束性を損なうことなく、内部変数と同等またはそれ以下の多項式次数を持つハイブリッド化変数の使用が可能となる。
   • LDG-H 法は、拡散的な流体方程式に適応される。
3. 静的縮約：HDG アプローチの定義的特徴は、静的縮約を通じて内部要素の DOFs を排除することである。その結果、グローバルな系は、固体変位（$\hat{v}_s$）と流体圧力（$\hat{p}$）の数値トレース変数のみを保持する。これは、標準的な DG 法と比較して、グローバルに結合された線形系のサイズを大幅に削減する。
4. 時間積分：完全離散スキームに対して、クランク・ニコルソン法が採用される。この陰的、2 次精度、非散逸スキームは、陽的スキーム（CFL 条件）が課す厳しい時間ステップ制限なしに系の剛性を処理し、かつ系の離散エネルギーを保持するように選択されている。
5. 誤差解析：本論文は、半離散定式化および完全離散クランク・ニコルソンスキームの厳密な事前誤差解析を提供する。この解析は、最適収束率を導出するために、特定の HDG 射影演算子とエネルギー恒等式に依存している。

主要な貢献

• 新規 HDG 定式化：本論文は、全応力テンソルを利用し、HDG+ と LDG-H 手法を組み合わせる、動的多孔弾性に対する新しい HDG 法を導入する。
• 対称応力近似：この手法は、混合定式化においてしばしば維持が困難な性質である、応力テンソルの数値近似が強く対称であることを保証する。
• ロックフリー特性：この定式化は、標準的な有限要素法で一般的に見られる体積ロックの問題を回避し、ほぼ非圧縮性材料に対して頑健であることが示される。
• 削減された計算コスト：静的縮約を通じて、この手法はグローバルに結合された DOFs の数を削減する。著者らは、高次近似（$k=3, 4$）の場合、標準的な DG アプローチと比較してグローバルな系のサイズが 77% 以上削減されることを実証している。
• 包括的な誤差解析：本論文は、半離散および完全離散スキームに対する最適収束率を確立し、固体変位、流体速度、圧力、および応力に関する具体的な推定を含む。

結果
数値実験は理論的知見を検証する：

• 収束率：製造された解を用いたテストは、応力、流体速度、および圧力に対して $O(h^{k+1})$ の最適収束次数を確認する。固体変位については、理論解析は $O(h^{k+1})$ を予測するが、数値実験は $L^2$ 誤差に対して $O(h^{k+2})$ の最適超収束率を観測する。
• 頑健性：この手法は、ほぼ非圧縮性材料（ポアソン比 $\nu = 0.499$）に対しても最適収束率を維持し、ロックフリー特性を確認する。
• 物理的精度：等方性、異方性、および不均質媒体における波動伝播のシミュレーションは、高速 P 波および散逸性の遅い P 波を含む複雑な波動場を成功裡に解像する。この手法は、人工的な振動なしに、材料界面（例えば、頁岩 - 砂岩境界）における波動の反射と透過を正確に捉える。

意義と主張
本論文は、提案された HDG 法が多孔弾性波動伝播のための効率的かつ正確な枠組みを提供すると主張している。HDG の利点（削減されたグローバル DOFs、高次収束）を活用し、多孔弾性における特定の課題（対称応力、トラクション境界条件、非圧縮性）に対処することにより、この手法は標準的な DG 法および有限要素法に対する実行可能な代替手段を提供する。著者らは、全応力テンソルの使用が境界条件の処理を簡素化し、HDG+ と LDG-H 戦略の組み合わせが、複雑な地球物理学的設定における高解像度波動シミュレーションに適した計算効率的なスキームを生み出すことを強調している。この研究は、高次解像度と計算効率が不可欠な地球物理学および石油工学の応用において、この手法が非常に適しているとの結論で締めくくられる。