[論文レビュー] Computational multiscale methods for parabolic wave approximations in heterogeneous media

Eric T. Chung, Yalchin Efendiev+2 more|arXiv (Cornell University)|Apr 6, 2021

Advanced Mathematical Modeling in Engineering参考文献 47被引用数 9

ひとこと要約

本稿では、時間方向に後退オイラー差分法、空間方向に多スケール還元のための制約エネルギー最小化一般化マルチスケール有限要素法（CEM-GMsFEM）、および時間発展方向の次元削減のための固有直交分解（POD）を組み合わせることで、極めて不均質な媒体における放物型波近似を解くための計算マルチスケール手法を提案する。数値実験では、相対誤差が10⁻³〜10⁻⁴程度の安定的で高精度な解を得ており、変動する媒体および層状媒体において、効率的で頑健な性能を示している。

ABSTRACT

In this paper, we develop a computational multiscale to solve the parabolic wave approximation with heterogeneous and variable media. Parabolic wave approximation is a technique to approximate the full wave equation. One benefit of the method is that: one wave propagation direction can be taken as an evolution direction, and we then can discretize it using a classical scheme like Backward Euler. Consequently, we obtain a set of quasi-gas-dynamic (QGD) models with different heterogeneous permeability fields. Then, we employ constraint energy minimization generalized multiscale finite element method (CEM-GMsFEM) to perform spatial discretization for the problem. The resulting system can be solved by combining the central difference in time evolution. Due to the variable media, we apply the technique of proper orthogonal decomposition (POD) to further the dimension of the problem and solve the corresponding model problem in the POD space instead of in the whole multiscale space spanned by all possible multiscale basis functions. We prove the stability of the full discretization scheme and give the convergence analysis of the proposed approximation scheme. Numerical results verify the effectiveness of the proposed method.

研究の動機と目的

• 高不均一で変動する媒体における放物型波近似を解くための効率的な計算マルチスケールフレームワークの開発。
• 空間的不均一性とマルチスケール特徴により、各zレベルで準ガス動力学的（QGD）モデルを解く際の計算課題の解決。
• 空間方向およびz発展方向における問題の次元削減を、それぞれCEM-GMsFEMおよびPODを用いて実現。
• 放物型波モデルの完全な差分スキームの安定性と収束性の確立。
• 層状およびMarmousi型の不均一な媒体における、手法の精度と効率の数値的妥当性の検証。

提案手法

• z微分を後退オイラースキームで離散化し、各zレベルで2次元QGDモデルの系列に変換する。
• 粗メッシュ上での空間離散化にCEM-GMsFEMを用い、局所的な不均一性を捉えるマルチスケール基底関数を構築する。
• 得られたQGDモデルにおける時間微分 ∂ₜₜv および ∂ₜv を中心差分法で離散化する。
• エネルギーノルムにおける射影誤差を最小化するように、組み合わせたCEM基底関数上に固有直交分解（POD）を適用し、低次元モデルを構築する。
• POD空間内で低次元系を解き、自由度を著しく削減しつつも精度を保持する。
• 後退オイラースキームの無条件安定性を証明し、CEM-GMsFEMの粗メッシュ収束結果に基づいて収束解析を提供する。

実験結果

リサーチクエスチョン

• RQ1 マルチスケール有限要素法は、変動する媒体における放物型波近似の空間的不均一性を効果的に扱えるか？
• RQ2 時間依存のマルチスケール問題において、精度を損なわずにz発展方向の次元をどのように削減できるか？
• RQ3 CEM-GMsFEMとPODを組み合わせることで、完全な差分スキームの安定性と収束性にどのような影響を与えるか？
• RQ4 層状および複雑な不均一な媒体（例：Marmousi領域）において、相対誤差はどのようにスケーリングするか？
• RQ5 提案手法は、細かいメッシュソルバーと比較して、著しく計算コストを低減しつつも高精度を達成できるか？

主な発見

• 最初の実験では、層構造の不均一な媒体において、相対誤差が約10⁻³に達している。
• Marmousi実験では、最終時間ステップでの相対誤差は約3.9×10⁻⁴であり、複雑で現実的な媒体においても高精度であることが示された。
• 最初の実験では、300個のCEM基底関数から抽出した100個のPOD基底関数が、主要モードを最小限の射影誤差で捉えていた。
• Marmousiケースでは、300個のCEM基底関数から50個のPOD基底関数を選択したため、効果的な次元削減が達成された。
• 後退オイラースキームは無条件安定性を保証しており、提案されたフレームワーク下で完全スキームの安定性と収束性が証明された。
• 数値結果により、構造的および非構造的な不均一な媒体、特に現実的な地質的配置を含む場合にも、手法の効率性と頑健性が確認された。