[論文レビュー] Numerical study of blow-up in solutions to generalized Kadomtsev-Petviashvili equations

Christian Klein, Ralf Peter|arXiv (Cornell University)|Oct 19, 2013

Nonlinear Waves and Solitons参考文献 13被引用数 4

ひとこと要約

本稿では、動的スケーリングを用いて特異性の型を分類するため、臨界および超臨界非線形性を有する一般化されたKadomtsev-Petviashvili（gKP）方程式の解における爆発現象について数値的検討を行う。爆発の状況は明確に異なる：gKP Iにおいて $ n = 4/3 $ のとき、爆発は $ L^2 $-臨界ソリトン型プロファイルに従い、$ \|u\|_{\infty} \propto (t^* - t)^{-3/4} $ となる。一方、$ n > 4/3 $ の場合、爆発は指数的スケーリング則を示し、$ \|u\|_{\infty} \propto (t^* - t)^{-2/(3n)} $ となる。これは理論的予想を裏付けるものである。

ABSTRACT

We present a numerical study of solutions to the generalized Kadomtsev-Petviashvili equations with critical and supercritical nonlinearity for localized initial data with a single minimum and single maximum. In the cases with blow-up, we use a dynamic rescaling to identify the type of the singularity. We present a discussion of the observed blow-up scenarios.

研究の動機と目的

• 一般化されたKadomtsev-Petviashvili（gKP）方程式の解が有限時間にわたって爆発する条件を調査すること。
• 動的スケーリング技術を用いて、爆発状況における特異性の型を分類すること。
• gKP II方程式において $ n > 2 $ の場合に爆発が発生するかどうかを特定し、その場合に限ってその漸近的挙動を特徴づけること。
• 初期データの質量およびエネルギーが、特にソリトン解と関連して爆発を引き起こす役割を調査すること。
• gKP IおよびgKP II方程式における爆発に関する予想された漸近的スケーリング則を数値的に検証すること。

提案手法

• 空間的領域では擬スペクトル法を用い、時間領域では分割ステップフォーリエ法を用いた適応的時間刻み法によりgKP方程式を数値的に解く。
• 爆発問題を定常問題に変換するため、動的スケーリングを適用し、爆発の型を特定可能にする。
• スケーリング変換 $ u(x,y,t) \to L(t)^{-2/n} u(x/L(t), y/L^2(t), t^*) $ を用い、$ L(t) \propto (t^* - t)^{\alpha} $ とすることでスケーリング指数を抽出する。
• $ \ln \|u\|_{\infty} $ および $ \ln \|u_y\|_2^2 $ を $ c \ln(t^* - t) + C $ にフィッティングし、爆発速度を特定する。
• 初期データを $ u_0(x,y) = \beta \partial_{xx} \exp(-x^2 - y^2) $ として構築し、$ \partial_x^{-1} u \in \mathcal{S}(\mathbb{R}^2) $ を満たし、制約条件 $ \int \partial_{yy} u \, dx = 0 $ を満たすようにする。
• $ u_y $ の $ L^2 $-ノルムの分析により爆発を確認し、$ t^* $ で $ \|u_y\|_2 \to \infty $ となることから、理論的予想と整合的であることを示す。

実験結果

リサーチクエスチョン

• RQ1 gKP IおよびgKP II方程式において、それぞれ $ n > 4/3 $ および $ n > 2 $ の場合に、初期データのどのような条件下で爆発が発生するか。
• RQ2 爆発時刻 $ t^* $ に近い領域における $ \|u\|_{\infty} $ および $ \|u_y\|_2 $ の漸近的スケーリング挙動は何か。
• RQ3 gKP Iにおいて $ n = 4/3 $ のとき、爆発プロファイルはスケーリングされたソリトン解と整合的か。
• RQ4 gKP Iにおいて $ n > 4/3 $ およびgKP IIにおいて $ n > 2 $ の場合、爆発が予想された指数的スケーリングに従うか。
• RQ5 初期データの質量およびエネルギーは、特にソリトン解と比較して爆発の発生にどのように関与するか。

主な発見

• gKP Iにおいて $ n = 4/3 $ のとき、初期データのエネルギーが十分に小さく、質量が大きい場合に爆発が発生し、$ \|u\|_{\infty} \propto (t^* - t)^{-3/4} $ となる。これは $ L^2 $-臨界ソリトン型プロファイルと一致する。
• $ u_y $ の $ L^2 $-ノルムは $ \|u_y\|_2 \propto (t^* - t)^{-1} $ と発散し、$ L^2 $-臨界爆発に対する理論的予測と整合的である。
• gKP Iにおいて $ n > 4/3 $ の場合、爆発速度は $ \|u\|_{\infty} \propto (t^* - t)^{-2/(3n)} $ と表され、非ソリトン的で自己相似的なプロファイルが指数的スケーリングを示すことを示唆する。
• $ \ln \|u\|_{\infty} $ を $ c_1 \ln(t^* - t) + C_1 $ にフィッティングした結果、$ c_1 = -0.1623 $ が得られ、期待される $ -1/6 $ に非常に近い。これは $ n = 4/3 $ の場合に $ -2/(3n} $ スケーリングが正当化されることを確認する。
• gKP IIにおいて $ n > 2 $ の場合、大きな質量と小さなエネルギーを有する初期データに対して爆発が観測され、$ \|u\|_{\infty} \propto (t^* - t)^{-2/(3n)} $ とならず、両gKP型に共通する普遍的な爆発メカニズムを示唆する。
• $ n > 4/3 $ の場合の $ u_y $ の $ L^2 $-ノルムは $ \|u_y\|_2 \propto (t^* - t)^{-(1+4/n)/6} $ とスケーリングし、フィッティングにより $ c_2 = -0.858 $ が得られ、期待される $ -2/3 $ に近い。これにより、予想されたスケーリングが支持される。