[論文レビュー] Numerical Study of Blowup in the Davey-Stewartson System

Christian Klein, Benson Muite+1 more|arXiv (Cornell University)|Dec 17, 2011

Advanced Mathematical Physics Problems参考文献 19被引用数 6

ひとこと要約

本稿では、スペクトル法と4次精度の分割スキームを用いて、収束型Davey-Stewartson II方程式における爆発ダイナミクスを数値的に調査する。両方の解、すなわちラムプ解とOzawaの正確解が不安定であることが示された。摂動により有限時間内に爆発が生じ、特にOzawaの爆発は近傍の初期データに対して一般性を持つことが確認され、非線形領域における爆発が安定した特徴であることが裏付けられた。

ABSTRACT

Nonlinear dispersive partial differential equations such as the nonlinear Schrödinger equations can have solutions that blow-up. We numerically study the long time behavior and potential blowup of solutions to the focusing Davey-Stewartson II equation by analyzing perturbations of the lump and the Ozawa exact solutions. It is shown in this way that the lump is unstable to both blowup and dispersion, and that blowup in the Ozawa solution is generic.

研究の動機と目的

• 収束型Davey-Stewartson II方程式における正確解の長期的挙動と安定性を調査すること。
• 既知の正確解に近い初期データに対して爆発が一般現象であるかどうかを特定すること。
• 小さな摂動に対するラムプ解およびOzawaの爆発解の安定性を評価すること。
• エネルギー保存則とフーリエ分解解像度を用いて、爆発検出のための数値手法の妥当性を検証すること。
• 収束型DS II方程式の解は、中間的挙動を示さず、爆発または分散のいずれかに分類されるという予想を検証すること。

提案手法

• 高次元空間分解能を達成するため、並列化されたスペクトルコードを用い、フーリエ擬スペクトル空間離散化を実施する。
• 無条件安定な4次精度の分割ステップ時間積分スキームを用いることで、長時間にわたるシミュレーションの精度と安定性を確保する。
• 数値的に計算されたエネルギーを、爆発の主要な指標として監視し、正則性の喪失を示す急激な上昇を観察する。
• ガウス関数および変形された有理関数形を用いて、ラムプ解およびOzawaの正確解の摂動を初期データとして構築する。
• アリヤシング誤差を回避するため、中間時間におけるフーリエ係数を分析する。
• 大規模計算を処理するため、GENCI、NERSC、NICSなどの高性能コンピューティングリソースを活用してシミュレーションを実行する。

実験結果

リサーチクエスチョン

• RQ1 収束型Davey-Stewartson II方程式において、Ozawa解の爆発挙動は近傍の初期データに対して一般性を持つか？
• RQ2 ラムプ解は小さな摂動を受けると、爆発を示すか、分散を示すか？
• RQ3 非局所項を含む2+1次元非線形分散型PDEに対して、数値的手法は爆発を信頼性高く検出できるか？
• RQ4 既知の正確解に近い初期データに対して、爆発が避けられないエネルギー閾値や条件は存在するか？
• RQ5 DS II方程式は、予想されるように、一般の初期データに対して爆発または分散の二択を示すか？

主な発見

• ラムプ解は不安定である。エネルギーを追加した摂動は有限時間内に爆発を引き起こすが、低エネルギーの摂動では分散が生じる。
• Ozawa解に対して、$ D = 0.1 $ の摂動では $ t_c \sim 0.2332 $ で爆発が発生し、$ D = 0.5 $ の場合はより早く $ t_c \sim 0.1659 $ で爆発が発生する。これは一般の爆発挙動を示している。
• $ \nu = 0.9 $ および $ \nu = 1.1 $ のOzawa初期データの変形では爆発は観測されず、それぞれ $ t = 0.2441 $ および $ t = 0.2254 $ で最大値を示し、エネルギーは保存されている。
• 数値的エネルギー計算は、解の正則性を失うと急激に値が上昇するため、爆発を信頼性高く検出できる。
• フーリエ係数による検証により、空間分解能が高く、爆発時付近でさえもアリヤシング誤差のない正確なシミュレーションが実現されている。
• 結果は、収束型DS II方程式の解が中間的長期的正則挙動を示さず、爆発または分散のいずれかに分類されることを支持する。