[論文レビュー] Multi-scale methods for wave propagation in heterogeneous media

Björn Engquist, Henrik Holst+1 more|ArXiv.org|Nov 13, 2009

Advanced Mathematical Modeling in Engineering参考文献 11被引用数 4

ひとこと要約

本稿では、急激に振動する係数を有する不均質な媒質における波動伝播問題を解くための非一様マルチスケール法（HMM）を提示する。マクロスケールの波動シミュレーションと局所的なマイクロスケールの解法を結合することで、小さなスケールεにほぼ依存しない計算効率を達成し、2次精度を維持するとともに、長時間にわたる分散効果を捉えることができる。1次元、2次元、3次元の周期的および非周期的問題において検証済みである。

ABSTRACT

Multi-scale wave propagation problems are computationally costly to solve by traditional techniques because the smallest scales must be represented over a domain determined by the largest scales of the problem. We have developed and analyzed new numerical methods for multi-scale wave propagation in the framework of heterogeneous multi-scale method. The numerical methods couples simulations on macro- and micro-scales for problems with rapidly oscillating coefficients. We show that the complexity of the new method is significantly lower than that of traditional techniques with a computational cost that is essentially independent of the micro-scale. A convergence proof is given and numerical results are presented for periodic problems in one, two and three dimensions. The method is also successfully applied to non-periodic problems and for long time integration where dispersive effects occur.

研究の動機と目的

• 高周波数成分を有する係数を伴う波動方程式を解く際の従来手法の高い計算コストを低減すること。
• マクロスケールとマイクロスケールのシミュレーションを結合することで、計算複雑度を低減するマルチスケール数値フレームワークの構築。
• 周期的係数を有する多次元設定におけるHMMの収束性を厳密に証明すること。
• 非周期的問題および長時間積分における分散効果の捉え込みの観点から、手法の頑健性を示すこと。
• 均質化方程式が未知または適用不可能な状況において、マルチスケール波動ソルバーの実用的応用を可能にすること。

提案手法

• マクロスケールの波動解を、局所的なマイクロスケールの解法から得たフラックスによって進める、非一様マルチスケール法（HMM）フレームワークを採用する。
• 空間的および時間的精度を確保するため、マクロスケールおよびマイクロスケールの両方で2次精度の有限差分スキームを用いる。
• 各マクロスケールグリッド点の周囲の小さなサブドメイン上でマイクロスケール波動方程式を解き、マクロスケールグリッド点における有効フラックスを計算する。
• 多次元ステンシルにおける2次精度を維持するために、マイクロスケールでの空間微分をステaggered平均で近似する。
• 周期的および非周期的係数構造を有する1次元、2次元、3次元問題に本手法を適用する。
• マクロスケールおよびマイクロスケールの両方で、Newmarkに類似したループフロッグ法に基づく時間積分スキームを用い、安定性と精度を確保する。

実験結果

リサーチクエスチョン

• RQ1 HMMフレームワークは、多次元における周期的係数を有するマルチスケール波動方程式に厳密に適用可能か？
• RQ2 HMMの収束速度および計算複雑度は、不均質媒質における波動伝播に対してどのように評価できるか？
• RQ3 均質化解が非自明な振る舞いを示す長時間波動シミュレーションにおいて、HMMは分散効果をどの程度正確に捉えられるか？
• RQ4 古典的均質化理論が適用できない非周期的係数構造へとHMMを拡張可能か？
• RQ5 標準的手法と比較して、HMMの計算コストは小スケールパラメータεに対してどのようにスケーリングするか？

主な発見

• HMMは、周期的係数を有する1次元、2次元、3次元の波動伝播問題において、空間的および時間的両方で2次精度を達成する。
• 本手法の計算コストは、小スケールεにほぼ依存せず、εをグローバルに解像する必要がある標準的手法に比べて顕著な利点を示す。
• 多次元周期的ケースにおけるHMMの収束性について、均質化理論の理論的期待値と整合する誤差境界を示す厳密な証明が確立された。
• 本手法は、均質化解に明示的にモデル化されていない項を含む長時間分散効果を効果的に捉えることができた。
• 数値結果により、非周期的係数問題に対しても頑健な性能を示しており、理論的仮定を越えた応用可能性が確認された。
• 分散項をマクロスケール方程式に明示的にモデル化する必要なく、長時間積分においても安定性と精度を維持した。