[論文レビュー] Error Estimates for Gaussian Beam Superpositions

Hailiang Liu, Olof Runborg+1 more|arXiv (Cornell University)|Aug 7, 2010

Advanced Mathematical Physics Problems参考文献 32被引用数 7

ひとこと要約

本稿は、厳密な双曲型PDEおよびシュレーディンガー方程式に対して、高周波初期データを伴うガウスビーム重ね合わせの誤差推定を確立する。波動方程式ではエネルギーノルム、シュレーディンガー方程式では$L^2$ノルムにおいて、$k$次ビーム重ね合わせが$mathcal{O}(\varepsilon^{k/2})$の収束率を示すことを証明した。これは曲率点が存在する場合でも成り立つ。数値結果は理論的収束率を確認しており、奇数次ビームでは誤差相殺により収束が改善される可能性が示唆されている。

ABSTRACT

Gaussian beams are asymptotically valid high frequency solutions to hyperbolic partial differential equations, concentrated on a single curve through the physical domain. They can also be extended to some dispersive wave equations, such as the Schrödinger equation. Superpositions of Gaussian beams provide a powerful tool to generate more general high frequency solutions that are not necessarily concentrated on a single curve. This work is concerned with the accuracy of Gaussian beam superpositions in terms of the wavelength $ε$. We present a systematic construction of Gaussian beam superpositions for all strictly hyperbolic and Schrödinger equations subject to highly oscillatory initial data of the form $Ae^{iΦ/ε}$. Through a careful estimate of an oscillatory integral operator, we prove that the $k$-th order Gaussian beam superposition converges to the original wave field at a rate proportional to $ε^{k/2}$ in the appropriate norm dictated by the well-posedness estimate. In particular, we prove that the Gaussian beam superposition converges at this rate for the acoustic wave equation in the standard, $ε$-scaled, energy norm and for the Schrödinger equation in the $L^2$ norm. The obtained results are valid for any number of spatial dimensions and are unaffected by the presence of caustics. We present a numerical study of convergence for the constant coefficient acoustic wave equation in $\Real^2$ to analyze the sharpness of the theoretical results.

研究の動機と目的

• 厳密な双曲型PDEおよびシュレーディンガー方程式に従う高周波数波動伝播におけるガウスビーム重ね合わせの収束率を厳密に確立すること。
• 従来の幾何光学が失敗する曲率点が存在する状況における重ね合わせの精度を分析すること。
• 収束率が空間次元に依存せず、曲率点の形成に対しても影響を受けないことを証明すること。
• $bb{R}^2$における音響波動方程式に対して、理論的誤差推定を数値的に検証すること。
• 奇数次ビームにおいて、隣接するビーム間の誤差相殺が原因で収束率が向上するかどうかを調査すること。

提案手法

• 光線に沿った位相および振幅のテイラー展開を用いた$k$次ビーム重ね合わせの体系的構築。
• 誤差が適切なノルムにおいて真の解と重ね合わせの間に抑えられることを示すために、振動積分の推定を用いる。
• 適切なエネルギーおよび$L^2$ノルムの定義のため、適切な初期値問題の適切性推定を適用する。
• 漸近的収束を検証するため、$k=1,2,3$に対して複数の時刻における$u_k - u$のスケーリング済みエネルギーノルムを数値的に計算する。
• ビームの局在化と数値実験における不自然な振動の低減のため、パrameter $eta = 1/10$ を持つカットオフ関数を組み込む。
• 異なる$ varepsilon$領域における収束挙動の分析により、漸近的減衰とカットオフに起因する誤差を区別する。

実験結果

リサーチクエスチョン

• RQ1 厳密な双曲型PDEに対して、$k$次ガウスビーム重ね合わせの収束率は、波長$ varepsilon$に対してどのように変化するか？
• RQ2 曲率点がガウスビーム重ね合わせの収束率に影響を与えるか？
• RQ3 2次元空間における音響波動方程式に対して、理論的誤差推定$mathcal{O}(\varepsilon^{k/2})$は数値的に検証可能か？
• RQ4 3次ビームの誤差が$mathcal{O}(\varepsilon^2)$に減少する一方で、2次ビームも$mathcal{O}(\varepsilon^1)$の挙動を示すのはなぜか？（高次であるにもかかわらず）
• RQ5 隣接するビーム間の誤差相殺により、奇数次重ね合わせにおいて収束率が向上する可能性はあるか？

主な発見

• 音響波動方程式において、$k$次ガウスビーム重ね合わせは、$ varepsilon$スケーリング済みエネルギーノルムにおいて真の解へ$ mathcal{O}(\varepsilon^{k/2})$の速度で収束する。
• シュレーディンガー方程式においては、$L^2$ノルムにおいて同様に$ mathcal{O}(\varepsilon^{k/2})$の収束率を示す。
• 数値結果により、$k=1,2,3$においてエネルギーノルムにおける理論的収束率$ mathcal{O}(\varepsilon^{k/2})$が確認された。$k=1,2$では$ mathcal{O}(\varepsilon^1)$、$k=3$では$ mathcal{O}(\varepsilon^2)$の収束が観察された。
• 数値例において$t=0.5$および$t>0.5$で、カスプおよびフォールド曲率点が存在する場合でも、収束率に影響がないことが示された。
• 中程度の$ varepsilon$領域では、2次ビームがカットオフ関数に起因する偶然の誤差相殺を経験し、予測よりも速い収束を示した。
• 3次ビームは、理論的推定値$ mathcal{O}(\varepsilon^{1.5})$よりも速い$ mathcal{O}(\varepsilon^2)$の減少を示し、ビーム間の誤差相殺により理論値を上回る収束が可能である可能性が示唆された。