[論文レビュー] Numerical study of a multiscale expansion of KdV and Camassa-Holm equation

Тамара Грава, Christian Klein|ArXiv.org|Feb 12, 2007

Nonlinear Waves and Solitons参考文献 22被引用数 23

ひとこと要約

本稿では、勾配崩壊付近における双曲型方程式のハミルトニアン摂動に関する普遍性予想を数値的に検証し、KdV方程式およびCamassa-Holm方程式の解が崩壊に近い領域で、Painlevé I階層の第二形（PI2）によって漸近的に記述されることを示している。PI2に基づくマルチスケール解は、古典的漸近理論と比較して臨界点付近で優れた近似を与える。KdVでは振動が先行側に、Camassa-Holmでは振動が後行側に現れるという明確な差異が観察される。

ABSTRACT

We study numerically solutions to the Korteweg-de Vries and Camassa-Holm equation close to the breakup of the corresponding solution to the dispersionless equation. The solutions are compared with the properly rescaled numerical solution to a fourth order ordinary differential equation, the second member of the Painlevé I hierarchy. It is shown that this solution gives a valid asymptotic description of the solutions close to breakup. We present a detailed analysis of the situation and compare the Korteweg-de Vries solution quantitatively with asymptotic solutions obtained via the solution of the Hopf and the Whitham equations. We give a qualitative analysis for the Camassa-Holm equation

研究の動機と目的

• ハミルトニアン摂動を受ける双曲型方程式の解が勾配崩壊付近でPainlevé I階層によって漸近的に記述されることを示す普遍性予想を数値的に検証すること。
• 分散なし方程式の崩壊点付近において、マルチスケール漸近解（PI2）と古典的Lax-LevermoreおよびDeift-Venakides-Zhou理論の精度を比較すること。
• KdVおよびCamassa-Holm方程式の解が勾配崩壊付近で示す振るまいの定性的および定量的差異、特に振動の位置と近似精度の違いを調査すること。
• コロケーション法を用いて、所定の漸近的境界条件を満たす4階常微分方程式（PI2）を数値的に解き、残差チェックにより解の精度を検証すること。

提案手法

• 空間的・時間的精度をスペクトル的に保つ擬似スペクトル法を用いて、KdVおよびCamassa-Holm方程式の数値解を求める。
• Painlevé I階層の第二形（PI2）である4階常微分方程式を用いてマルチスケール漸近解を構築する：$ X = 6T U - \left[ U^3 + \frac{1}{2}U_X^2 + U U_{XX} + \frac{1}{10}U_{XXXX} \right] $。
• MATLABのbvp4cソルバーを用いて、有限区間$[X_l, X_r]$上でPI2常微分方程式をコロケーション法と立方スプライン補間を用いて数値的に解く。
• 無限遠における漸近的境界条件は、$ Y = X^{1/3} $に関するローラン級数を用いて近似し、再帰的に係数を計算することで一貫性を確保する。
• Chebyshevスペクトル微分法を用いた残差チェックにより、$ |T| < 1 $、$ |X| < 10 $ の範囲で解の精度が$ 10^{-4}$以内であることを検証する。
• 臨界点$ (x_c, t_c) $付近の複数の時間ステップにおいて、KdVおよびCH方程式の全数値解とマルチスケール解を比較する。

実験結果

リサーチクエスチョン

• RQ1 Painlevé I階層（PI2）に基づくマルチスケール解は、KdVおよびCamassa-Holm方程式の解が勾配崩壊点付近で有効な漸近的記述を提供するか？
• RQ2 PI2に基づくマルチスケール解の近似精度は、崩壊時刻付近において古典的Lax-LevermoreおよびDeift-Venakides-Zhou漸近理論と比較してどうなるか？
• RQ3 KdVおよびCamassa-Holm方程式の解が臨界点付近で振動を形成する様子に、空間的分布（先行側対後行側）においてどのような定性的な差異が生じるか？
• RQ4 PI2常微分方程式の数値解が、普遍性予想が予言する普遍的振るまいをどの程度正確に再現するか？
• RQ5 $ \epsilon $ を小さくするに従って、マルチスケール近似の精度はどのように変化するか？また、ダブルスケーリング極限において改善が見られるか？

主な発見

• PI2常微分方程式に基づくマルチスケール解は、Lax-Levermore理論およびDeift-Venakides-Zhou理論と比較して、勾配崩壊付近におけるKdV方程式の解をより優れた漸近的記述で表している。
• KdV方程式では、振動が臨界点の先行側に現れ、マルチスケール解はそれらをよく近似するが、後行側ではやや精度が劣る。
• Camassa-Holm方程式では、振動が臨界点の後行側に現れ、マルチスケール解は先行部分を後行の振動よりもよく近似する。これはKdVの場合と同様の傾向を示すが、空間的構造が逆転している。
• Chebyshevスペクトル微分法による残差チェックにより、PI2常微分方程式の数値解が$ 10^{-4}$の精度で方程式を満たしていることが確認され、漸近的アンザッツの信頼性が裏付けられた。
• CH方程式の解が崩壊付近で示す定性的な振るまい—伝搬方向の左側に振動がなく、後行側に振動が現れる—は、KdVの場合とは顕著に異なり、分散的ダイナミクスの違いを示している。
• $ \epsilon $ が小さくなるに従って、CH方程式の数値解とマルチスケール解の一致度が向上する。これはダブルスケーリング極限における期待通りであり、PI2解の漸近的普遍性が確認された。