[論文レビュー] On the Polysymplectic Integrator for the Short Pulse Equation

Monika E. Pietrzyk, Igor V. Kanatchikov|arXiv (Cornell University)|Dec 30, 2015

Numerical methods for differential equations参考文献 13被引用数 4

ひとこと要約

本稿では、短パルス方程式（SPE）の幾何的多symplectic積分法を、De Donder-Weyl ハミルトニアン形式を用いて提案する。空間および時間の両方で中点法を用いて離散化することで、長時間の伝搬において、擬似スペクトル法に比べて1桁の精度向上と2.5倍の高速化を達成し、エネルギー保存則を満たすとともに、安定性を確保する。

ABSTRACT

The polysymplectic analysis of the Short Pulse Equation known in nonlinear optics is used in order to construct a geometric polysymplectic integrator for it. The proposed scheme turns out to be much more effective than other standard integration schemes for nonlinear PDEs, such as the pseudo-spectral integrator. In our numerical experiments the polysymplectic integrator appears to be an order of magnitude more precise and approximately 2.5 times faster at long propagation times than the pseudo-spectral method.

研究の動機と目的

• 短パルス方程式（SPE）の幾何的数値積分法を構築し、その背後にある多symplectic構造を保存すること。
• 緩やかに変化する包絡線近似が破綻するため非線形シュレーディンガー方程式が不適切となる超短パルスのシミュレーションにおいて、数値的安定性と精度を向上させること。
• 精度、速度、エネルギー保存則の観点から、提案された多symplectic積分法と標準的手法（例：擬似スペクトル法）の性能を比較すること。
• 非線形光学分野におけるSPEおよびその一般化の長時間シミュレーションに耐えうる、強固で効率的なスキームを確立すること。

提案手法

• 一次ラグランジアンを用いてSPEを定式化し、多運動量 $ p^t $ と $ p^x $ を用いたDe Donder-Weyl（DW）ハミルトニアン形式を導出する。
• Duffin-Kemmer-Petiau行列 $ \beta^x $ と $ \beta^t $ を用いて、DWハミルトニアン方程式を行列形式で表現し、多symplectic構造を実現する。
• 空間 $ x $ および時間 $ t $ の両方で2次精度の中点法を用いてDW方程式を離散化し、$ \phi_{i,j+1/2} $ や $ p^t_{i+1/2,j} $ のようなステaggerグリッド点を用いる。
• 各グリッド点で得られる非線形代数方程式系を、式 (18a) の立方方程式を用いて $ p^t_{i,j+1} $ について解き、その後 $ \phi_{i,j+1} $ と $ p^x_{i,j+1} + p^x_{i,j} $ を明示的に更新する。
• 再構成された電場振幅は $ u_{i,j+1} = 2p^t_{i,j+1} $ で与えられ、元のSPEの $ u_{xt} = u + \frac{1}{6}(u^3)_{xx} $ と整合性を保つ。
• 周期的境界条件または吸収境界条件を適用し、初期条件は正確な Sakovich 解から導出する。エネルギー関数の保存則を強制的に満たす。

実験結果

リサーチクエスチョン

• RQ1 De Donder-Weyl形式を用いて、短パルス方程式のための幾何的多symplectic積分法を構築可能か？
• RQ2 長時間シミュレーションにおいて、提案された多symplectic積分法は、擬似スペクトル法に比べて精度と計算効率で優れているか？
• RQ3 多symplecticスキームはエネルギー関数を保存し、SPEの数値解において長期間にわたって安定性を維持できるか？
• RQ4 離散化スキームは、元の多symplectic構造および関連する保存則をどの程度正確に保持しているか？

主な発見

• 多symplectic積分法は、擬似スペクトル法に比べて約10倍の高い精度を達成し、$ N = 2^{17} $ および $ \Delta t = 0.0001 $ の条件下で $ t = 100 $ 時に $ \sigma \approx 6.5 \times 10^{-6} $ の誤差を示す。
• 計算時間は約2.5倍高速で、同じ条件下で擬似スペクトル法の100,000秒に対して本手法は40,000秒のCPU時間を要する。
• 本スキームは長時間にわたる伝搬においてもエネルギー関数を保存し、標準的な積分スキームに比べて優れた耐障害性と安定性を示す。
• 式 (18a) の立方方程式により $ p^t_{i,j+1} $ が連続的になることが保証され、根の選択プロセスにより数値的一致性が維持される。
• 数値的結果は、多symplectic積分法が、正確なSakovich解との比較によって検証されたように、長時間にわたり精度と安定性を維持していることを確認する。
• 本手法の有効性は、複数の離散化ステップにおいて確認されており、$ \Delta x $ および $ \Delta t $ を小さくするに従い、誤差が一貫して減少する。