[論文レビュー] Frozen Gaussian approximation for high frequency wave propagation in periodic media

Ricardo Delgadillo, Jianfeng Lu+1 more|arXiv (Cornell University)|Apr 30, 2015

Advanced Mathematical Modeling in Engineering参考文献 23被引用数 3

ひとこと要約

本稿では、ブロッホ分解と位相空間漸近解析を組み合わせることで、周期的媒体における高周波数波動伝播のためのフローズンガウス近似（FGA）を提案する。この手法は、小さなパラメータεに依存しない計算コストで真の解に収束することを実現し、従来の方法のメッシュ解像度制限を克服するとともに、ガウスビーム法に内在するビーム拡散問題を回避する。

ABSTRACT

Propagation of high-frequency wave in periodic media is a challenging problem due to the existence of multiscale characterized by short wavelength, small lattice constant and large physical domain size. Conventional computational methods lead to extremely expensive costs, especially in high dimensions. In this paper, based on Bloch decomposition and asymptotic analysis in the phase space, we derive the frozen Gaussian approximation for high-frequency wave propagation in periodic media and establish its converge to the true solution. The formulation leads to efficient numerical algorithms, which are presented in a companion paper [Delgadillo, Lu and Yang, arXiv:1509.05552].

研究の動機と目的

• 短い波長と小さな格子定数のため、周期的媒体における高周波数波動伝播のシミュレーションが計算的に非現実的であるという問題に対処すること。
• 特に高次元において顕著な、従来の数値的手法のメッシュ解像度制限を克服すること。
• フローズンガウス近似（FGA）を周期的媒体に拡張し、特異点（カウスティクス）付近でも効率性と精度を維持すること。
• 周期的で滑らかなポテンシャルを有するシュレーディンガー方程式の真の解へのFGA法の収束性を厳密に確立すること。

提案手法

• 周期的ポテンシャルを有するシュレーディンガー方程式に対して、ブロッホ分解と漸近解析を用いて位相空間積分表現を導出する。
• 各成分が固定幅を持ち、ハミルトニアンフローに従って時間発展するフローズンガウス関数の重ね合わせとしてFGA解を構築する。
• 振幅関数および位相関数の時間発展方程式を導出するために、εのべき級数展開を用いた摂動展開を実施する。
• ブロッホ作用素の固有値スペクトルを用いて、解成分の正則性および有界性を制御する。
• フーリエ積分作用素理論とエネルギー推定を用いて、FGA近似の収束性を証明する。
• ベッセルの不等式とブロッホバンド内のスペクトルギャップを用いて、解成分の一様有界性を確立する。

実験結果

リサーチクエスチョン

• RQ1 フローズンガウス近似は、周期的ポテンシャルを有する周期的媒体における高周波数波動伝播に拡張可能か？
• RQ2 周期的系において特異点（カウスティクス）が存在する状況下でも、FGA法は精度を保ち、ビーム拡散を回避できるか？
• RQ3 半古典的極限ε → 0における、FGA近似の真の解への収束速度はいかほどか？
• RQ4 周期的媒体における短波長を解像しつつ、計算コストをεに依存させずに抑えることは可能か？

主な発見

• FGA法は、半古典的領域においてシュレーディンガー方程式の真の解に収束し、誤差は近似解の残差にε⁻¹を乗じたもので有界であることが示された。
• 幾何光学やアイルコナル方程式法が破綻するような特異点付近でも、この手法は安定かつ高精度を維持する。
• ブロッホバンド内のスペクトルギャップのおかげで、振幅関数や位相関数を含む解成分が時間的・空間的に一様に有界である。
• FGA定式化は、εに依存しない計算コストを持つ数値アルゴリズムを導き、高次元における効率的シミュレーションを可能にする。
• エネルギー推定とブロッホ作用素の逆作用素の有界性に依拠して、収束性が確立されており、これはスペクトルギャップ条件に依存する。