[論文レビュー] Wave propagation across acoustic / Biot's media: a finite-difference method

Guillaume Chiavassa, Bruno Lombard|arXiv (Cornell University)|Sep 15, 2011

Lattice Boltzmann Simulation Studies参考文献 40被引用数 30

ひとこと要約

本稿では、Biotの式を用いて流体および多孔質弾性媒体を貫通する波動伝播をシミュレートする高次精度有限差分法を提示する。開口、密封、不完全な多孔質界面境界条件を組み込む。4次精度のADER時間積分、時空間メッシュ局所細分化、およびサブセル解像度を有する埋め込み界面法を組み合わせることで、拡散的性質を示す遅い縦波を正確に捉え、すべての界面タイプについて解析解と一致する高い精度を達成する。

ABSTRACT

Numerical methods are developed to simulate the wave propagation in heterogeneous 2D fluid / poroelastic media. Wave propagation is described by the usual acoustics equations (in the fluid medium) and by the low-frequency Biot's equations (in the porous medium). Interface conditions are introduced to model various hydraulic contacts between the two media: open pores, sealed pores, and imperfect pores. Well-possedness of the initial-boundary value problem is proven. Cartesian grid numerical methods previously developed in porous heterogeneous media are adapted to the present context: a fourth-order ADER scheme with Strang splitting for time-marching; a space-time mesh-refinement to capture the slow compressional wave predicted by Biot's theory; and an immersed interface method to discretize the interface conditions and to introduce a subcell resolution. Numerical experiments and comparisons with exact solutions are proposed for the three types of interface conditions, demonstrating the accuracy of the approach.

研究の動機と目的

• 流体／多孔質異種媒質における波動伝播を記述する力学的・Biotの式に従う堅牢な数値法の開発。
• 開放的多孔質、密封的多孔質、有限透過率を有する不完全多孔質の3種類の水力的界面境界条件の正確なモデル化。
• 拡散的性質を示す遅い縦波をシミュレートする課題に取り組み、界面近傍で細かい空間解像度を要する。
• 結合された流体／多孔質弾性系における初期境界値問題の適切な定式化を保証する。
• 非構造メッシュを用いずに、カルテシアングリッド上で複雑な幾何形状において高次精度と安定性を達成する。

提案手法

• 4次精度のADERスキームにStrang分割を適用し、最適なCFL安定性条件を満たす。
• 流体（音響）および多孔質弾性（Biotの式）領域に対して別々のソルバーを用い、埋め込み界面法により結合する。
• 水力的接触（開放的、密封的、不完全多孔質）の界面条件は、非整合な界面に対処可能なサブセル解像度技術を用いて離散化する。
• 遅い縦波は局在的かつ拡散的であるため、界面周辺に時空間メッシュ局所細分化戦略を適用して解像度を高める。
• 特異値分解を用いたテイラー展開に基づく再構成法により、界面点における修正された有限差分スタencilを計算し、一貫性と精度を確保する。
• 数値実験において正規化因子 𝒩 = 1000 を用いて、計算に適した正規化物理パラメータを用いる。

実験結果

リサーチクエスチョン

• RQ1 複雑な界面条件を有する流体および多孔質弾性媒質を貫通する波動伝播シミュレーションにおいて、どのように高次精度を達成できるか？
• RQ2 開放的、密封的、不完全多孔質の異なる水力的接触タイプが、波動の透過および反射に与える影響は何か？
• RQ3 非構造メッシュを用いずに、カルテシアングリッド上の有限差分法が、Biot理論が予測する遅い縦波を正確に解像できるか？
• RQ4 ADER時間積分、時空間メッシュ局所細分化、埋め込み界面法の組み合わせが、安定性および精度をどのように向上させるか？
• RQ5 提案手法が、任意の形状の界面および変化する材料特性に対しても、どの程度の精度を維持できるか？

主な発見

• 本手法は、流体および多孔質弾性媒質を貫通する波動伝播のシミュレーションにおいて高い精度を達成しており、3種類の界面条件すべてについて数値結果が解析解とよく一致する。
• 埋め込み界面法により、サブセル解像度を有する界面条件を正確に捉え、非整合な界面に対しても正確な取り扱いが可能である。
• 時空間メッシュ局所細分化は、界面近傍に局在し拡散的性質を示す遅い縦波を効果的に解像する。
• Strang分割を用いた4次精度ADERスキームは、高粘性を有する剛性系においても最適なCFL条件を満たし、安定性を維持する。
• 正規化変数の使用により数値的条件数が改善され、現実的な物理パラメータを用いた安定なシミュレーションが可能になる。
• 初期境界値問題の適切な定式化が厳密に証明されており、モデルの数学的整合性が保証される。