[論文レビュー] Frozen Gaussian approximation for general linear strictly hyperbolic system: formulation and Eulerian methods

Jianfeng Lu, Xu Yang|arXiv (Cornell University)|Oct 10, 2010

Quantum chaos and dynamical systems参考文献 12被引用数 4

ひとこと要約

本稿は、一般の線形厳密な双曲型系に、凍結ガウス近似（FGA）を拡張し、ラグランジュ的FGAに内在する軌道発散問題を克服するためのオイラー的数値手法を開発する。位相空間上のリウヴィル方程式を用いてFGAをオイラー的枠組みで定式化することで、焦点やビーム拡散が生じる状況下でも、安定的かつ高精度に高周波数波動伝搬を計算可能となる。この手法は波動方程式および量子力学におけるエルマン＝クルック propagator への応用を含む。

ABSTRACT

The frozen Gaussian approximation, proposed in [Lu and Yang, [15]], is an efficient computational tool for high frequency wave propagation. We continue in this paper the development of frozen Gaussian approximation. The frozen Gaussian approximation is extended to general linear strictly hyperbolic systems. Eulerian methods based on frozen Gaussian approximation are developed to overcome the divergence problem of Lagrangian methods. The proposed Eulerian methods can also be used for the Herman-Kluk propagator in quantum mechanics. Numerical examples verify the performance of the proposed methods.

研究の動機と目的

• 長時間の高周波数波動伝搬における粒子軌道発散に起因するラグランジュ的凍結ガウス近似の不安定性を解消する。
• 滑らかな係数を有する一般の線形厳密双曲型系へ、凍結ガウス近似フレームワークを一般化する。
• 固定グリッド上での安定計算を可能にするために、位相空間上のリウヴィル方程式に基づくオイラー的数値手法を開発する。
• 高周波数シュレーディンガー方程式に対するエルマン＝クルック propagator へのこの手法の応用を可能にする。
• 波長に比例した細かいメッシュを必要とする従来の手法に対する、数値的に安定した代替手法を提供する。

提案手法

• M×M 線形双曲型系へのFGAの拡張として、幅を凍結したガウス関数を用いた位相空間表現を定義する。
• 位相空間上のリウヴィル方程式を用いて、オイラー的枠組みにおける凍結ガウスビームの時間発展を定式化する。
• 固定位相空間グリッド上で局所的にリウヴィル方程式を解くことで、粒子軌道発散を回避する。
• リウヴィル方程式の時間積分を効率的に行うために、半ラグランジュ法および有限差分法を用いる。
• 適切なメッシュサイズを有する数値積分則を用いて、位相空間積分から物理空間解を再構成する。
• 音響波動方程式およびシュレーディンガー方程式の両方に対して、エルマン＝クルック propagator を通じてこの手法を適用する。

実験結果

リサーチクエスチョン

• RQ1 凍結ガウス近似は、スカラー波動方程式を超えて、一般の線形厳密双曲型系へどのように一般化可能か？
• RQ2 長時間の高周波数波動伝搬におけるラグランジュ的FGAの限界は何か？ それらはどのように克服できるか？
• RQ3 位相空間ダイナミクスに基づくオイラー的手法は、高周波数波動伝搬に対して安定的かつ高精度な解を提供できるか？
• RQ4 オイラー的FGA手法は、エルマン＝クルック propagator のような量子力学的系へ、どの程度応用可能か？
• RQ5 ビーム拡散や焦点の発生下でも、オイラー的FGAはラグランジュ的FGAおよびスペクトル法と比較して、精度と安定性の面でどのように優れるか？

主な発見

• 固定グリッド上での位相空間発展を用いることで、オイラー的FGA手法は、ラグランジュ的FGAの発散問題を効果的に解消した。
• 2次元音響波動系に対する数値結果では、スペクトル解と良好に一致し、ε=1/64、T=1.0 時に圧力および速度の相対誤差が5%未満であった。
• 拡散する波束を有する1次元シュレーディンガー方程式に対しては、オイラー手法はT=10を越えて精度を維持したが、ラグランジュ手法は軌道発散により失敗した。
• 調和ポテンシャルを有する2次元シュレーディンガー方程式において、オイラー的エルマン＝クルック propagator はT=0.5およびT=1において、拡散および局在化ダイナミクスを正確に捉えた。
• 波長分解能グリッドを必要とせず、中程度のメッシュサイズ（例：δx=1/64、δq=δp=1/32）で高精度を達成した。
• 提案されたオイラー的FGAフレームワークは、エルマン＝クルック propagator へ直接適用可能であり、高周波数初期データを有する量子ダイナミクスへの応用を拡張した。