[論文レビュー] A Riemann-Hilbert approach for the Degasperis-Procesi equation

Anne Boutet de Monvel, Dmitry Shepelsky|arXiv (Cornell University)|Jul 29, 2011

Numerical methods in inverse problems被引用数 16

ひとこと要約

本稿では、Degasperis-Procesi (DP) 方程式のためのリーマン・ヒルベルト (RH) 問題の定式化を開発し、逆散乱変換による解の表現を可能にした。この手法は、3×3 行列 RH 問題を用いて3次Lax対を扱い、明確な長時間漸近挙動を導出する。特定の領域では振動的減衰が生じ、他の領域では急速な減衰が観察され、位相変調とスペクトルデータに従う振幅減衰を含む明示的な式が得られる。

ABSTRACT

We present an inverse scattering transform approach to the Cauchy problem on the line for the Degasperis--Procesi equation $u_t-u_{txx}+3ωu_x+4uu_x=3u_xu_{xx}+uu_{xxx}$ in the form of an associated Riemann-Hilbert problem. This approach allows us to give a representation of the solution to the Cauchy problem, which can be efficiently used in studying its long-time behavior.

研究の動機と目的

• Degasperis-Procesi 方程式の逆散乱変換を、リーマン・ヒルベルト問題の定式化を用いて開発すること。
• 従来、Camassa-Holm 方程式に用いられたRH手法を、DP 方程式の3次Lax対構造へ拡張すること。
• 解 $ u(x,t) $ を、特徴的なスペクトルパラメータ点におけるRH問題の解の評価により表現すること。
• 異なる空間的・時間的領域における解の長時間漸近的挙動を分析すること。
• 明示的な位相および振幅依存性を持つ漸近公式を確立し、振動的および減衰的領域を含むこと。

提案手法

• スペクトルパラメータ $ z $ と3×3行列 $ \Phi $ を用いて、DP 方程式のLax対をスカラー形式および行列形式で定式化し、固有関数の微分を表現する。
• 複素 $ z $-平面上にリーマン・ヒルベルト問題を構築し、解 $ \Psi(z) $ が曲線の外で正則であり、散乱データによって定まるジャンプ条件を満たすようにする。
• 解 $ \Psi(z) $ の $ z \to \infty $ における漸近挙動を用いて、解 $ u(x,t) $ を抽出する。この式構造はパrametric multisoliton解と一致する。
• 非線形勾配降下法を用いて曲線を変形し、$ x/t $ の異なる領域における漸近的挙動を分析する。
• 位相関数 $ \Theta $ を用いて $ k $-平面上の臨界点を特定し、符号表を用いて振動的または減衰的挙動を同定する。
• 臨界点が重なる場合（例：$ x/t = -3/8 $ および $ x/t = 3 $）、遷移領域はPainlevé II超越関数を用いて取り扱う。

実験結果

リサーチクエスチョン

• RQ1 3次Lax対を有するDegasperis-Procesi方程式の逆散乱変換を、リーマン・ヒルベルト問題としてどのように定式化できるか？
• RQ2 DP 方程式の解の長時間漸近的挙動は、$ x/t $-平面上の異なる領域でどのように変化するか？
• RQ3 漸近的挙動における振動的および減衰的領域は、スペクトルデータと位相関数の解析からどのように生じるか？
• RQ4 臨界点および曲線変形が、RH 問題の非線形勾配降下解析において果たす役割は何か？
• RQ5 臨界点が重なる遷移領域（例：$ x/t = -3/8 $ および $ x/t = 3 $）では、Painlevé II 解が関与する特殊な漸近的挙動がどのように生じるか？

主な発見

• 解 $ u(x,t) $ は、3×3 行列リーマン・ヒルベルト問題の $ z \to \infty $ における展開により表現され、パラメトリック multisoliton 解と構造が一致する。
• 領域 $ 0 < x/t < 3 $ では、振動的減衰が観察され、$ u(x,t) = \frac{c_1}{\sqrt{t}} \sin(c_2 t + c_3 \log t + c_4)(1 + o(1)) $ と表される。$ c_1, c_2, c_3, c_4 $ は、スペクトル反射係数 $ r(s) $ を用いて明示的に表される。
• 振幅 $ c_1 $ は $ \kappa_0 $, $ h_0 $、および $ p_0^{-1} $ と $ p_0 $ における反射係数の位相に依存し、$ c_1 = -\left(\frac{192h_0\kappa_0}{(3-4\kappa_0^2)(1+4\kappa_0^2)}\right)^{1/2} \cos(\arg r(-p_0^{-1}) - \arg r(p_0)) $ と表される。
• 領域 $ x/t < -3/8 $ では、実数の臨界点が存在しないため、符号表および位相解析により、解は急速な指数的減衰を示す。
• 臨界点が重なる点 $ x/t = -3/8 $ および $ x/t = 3 $ では、漸近的挙動がPainlevé II方程式の解によって支配され、支配的挙動の変化が示される。
• 本手法により、2×2 のCamassa-Holm方程式に対するRH手法が、3×3 のDP方程式へと自然に拡張され、注意深く曲線変形とスペクトル解析を用いることで、詳細な長時間漸近挙動が得られた。