コルテウェグ・ド・フリース方程式

ヤコビ楕円関数cnの2 乗による Korteweg–De Vries 方程式のクノイド波解(パラメータmの値は0.9 )。
KdV方程式u t + u u x + δ 2 u x x x = 0 ( δ = 0.022 ) の数値解。初期条件はu ( x , 0) = cos(π x )である。時間発展はザブスキー・クラスカル法によって行われた。[ 1 ]初期の余弦波は孤立波の列へと発展する。
KdV方程式の2つのソリトン解が相互作用し(紫色)、それらが互いを通過する際に生じる位相シフトを強調している。赤と青の解は、一方のソリトンが存在しない状態での個々のソリトンの運動を示している[ 2 ]

数学において、コルテヴェク・ド・フリース(KdV)方程式は、浅い水面の波の数学的モデルとして機能する偏微分方程式(PDE)である。この方程式は積分可能なPDEの典型的な例として特に注目されており、PDEを一般的に扱いにくくする非線形性にもかかわらず、多数の明示的解(特にソリトン解)や無限数の保存量などの典型的な挙動を示す。KdVは逆散乱法(ISM)によって解くことができる。 [ 3 ]実際、クリフォード・ガードナージョン・M・グリーン、マーティン・クラスカル、ロバート・ミウラはKdV方程式を解くために古典的な逆散乱法を開発した。

KdV方程式は、 ジョセフ・ヴァレンティン・ブシネスク (1877年、360ページの脚注)によって初めて導入され、 1895年にディーデリック・コルテウェググスタフ・デ・フリースによって再発見され、最も単純な解である1ソリトン解が発見されました。[ 4 ] [ 5 ]方程式と解の挙動の理解は、1965年のノーマン・ザブスキーおよびクラスカルのコンピュータシミュレーションと、1967年の逆散乱変換の開発によって大きく進歩しました。

1972年にT.Kawaharaは、KdV方程式の係数が非常に小さくなるかゼロになる場合の分散波を記述するKawahara方程式として知られる5次のKdV型方程式を提案した。[ 6 ]

意味

KdV方程式は、(空間的に)1次元の非線形分散非散逸波をモデル化する偏微分方程式であり、次式に従う関数によって記述される: [ 7 ]ϕ×t{\displaystyle \phi (x,t)}

tϕ+×3ϕ6ϕ×ϕ0×Rt0{\displaystyle \partial _{t}\phi +\partial _{x}^{3}\phi -6\,\phi \,\partial _{x}\phi =0\,\quad x\in \mathbb {R} ,\;t\geq 0,}

ここで は分散を表し、非線形要素は移流項です。 ×3ϕ{\displaystyle \partial _{x}^{3}\phi }ϕ×ϕ{\displaystyle \phi \partial _{x}\phi }

浅瀬の波をモデル化する場合、水面の高さの平衡高さからの変位です。 ϕ{\displaystyle \phi }

最後の項の前の定数は慣例的なものですが、大きな意味はありません。 、、およびを定数で乗算すると、3 つの項のいずれかの係数を任意のゼロ以外の定数に等しくすることができます。 6{\displaystyle 6}t{\displaystyle t}×{\displaystyle x}ϕ{\displaystyle \phi }

ソリトンソリューション

1ソリトン解

で与えられる固定波形が、位相速度で右方向に移動する際にその形状を維持する解を考える。このような解は で与えられる。これをKdV方程式に代入すると、常微分方程式が得られる。fX{\displaystyle f(X)}c{\displaystyle c}φ×tf×ct1つのfX{\displaystyle \varphi (x,t)=f(x-ct-a)=f(X)}

cdfdX+d3fdX36fdfdX0{\displaystyle -c{\frac {df}{dX}}+{\frac {d^{3}f}{dX^{3}}}-6f{\frac {df}{dX}}=0,}

または、 について積分すると、 X{\displaystyle X}

cf+d2fdX23f2{\displaystyle -cf+{\frac {d^{2}f}{dX^{2}}}-3f^{2}=A}

ここで、積分定数である。上記の独立変数を仮想時間変数として解釈すると、これは ニュートンの立方体ポテンシャルにおける単位質量の粒子の 運動方程式を満たすことを意味する。A{\displaystyle A}X{\displaystyle X}f{\displaystyle f}

V(f)=(f3+12cf2+Af){\displaystyle V(f)=-\left(f^{3}+{\frac {1}{2}}cf^{2}+Af\right)}

もし

A=0,c>0{\displaystyle A=0,\,c>0}

とすると、ポテンシャル関数はで極大値をとります。 には、 の解が存在し、 は「仮想時間」 にこの点から始まり 、最終的に極小値まで滑り落ち、その後反対側へ上昇して同じ高さに達し、その後方向を反転して の時点で再び極大値に達します。つまり、は に近づくにつれて に近づきます。これが孤立波解の特徴的な形状です。 V(f){\displaystyle V(f)}f=0{\displaystyle f=0}f(X){\displaystyle f(X)}{\displaystyle -\infty }{\displaystyle \infty }f(X){\displaystyle f(X)}0{\displaystyle 0}X{\displaystyle X\to -\infty }

より正確に言えば、解決策は

ϕ(x,t)=12csech2[c2(xcta)]{\displaystyle \phi (x,t)=-{\frac {1}{2}}\,c\,\operatorname {sech} ^{2}\left[{{\sqrt {c}} \over 2}(x-c\,t-a)\right]}

ここで は双曲正割を表し、は任意の定数である。[ 8 ]これは速度 で右方向に動くソリトンを表す。 sech{\displaystyle \operatorname {sech} }a{\displaystyle a}c{\displaystyle c}

Nソリトン解

解を表す式が知られている。これは-ソリトン解であり、後期には別々の単一ソリトンに分解される。[ 9 ]この解は、減少する正のパラメータの集合と非ゼロのパラメータの集合に依存する。この解は、 行列の成分が N{\displaystyle N}N{\displaystyle N}χ1>>χN>0{\displaystyle \chi _{1}>\cdots >\chi _{N}>0}β1,,βN{\displaystyle \beta _{1},\cdots ,\beta _{N}}ϕ(x,t)=22x2log[detA(x,t)]{\displaystyle \phi (x,t)=-2{\frac {\partial ^{2}}{\partial x^{2}}}\mathrm {log} [\mathrm {det} A(x,t)]}A(x,t){\displaystyle A(x,t)}Anm(x,t)=δnm+βne8χn3te(χn+χm)xχn+χm.{\displaystyle A_{nm}(x,t)=\delta _{nm}+{\frac {\beta _{n}e^{8\chi _{n}^{3}t}e^{-(\chi _{n}+\chi _{m})x}}{\chi _{n}+\chi _{m}}}.}

これは逆散乱法を使用して導出されます。

運動の積分

KdV方程式には、時間とともに変化しない解の関数である運動の積分が無限に存在する。 [ 10 ]これらは次のように明示的に表すことができる。 ϕ(t){\displaystyle \phi (t)}

+P2n1(ϕ,xϕ,x2ϕ,)dx{\displaystyle \int _{-\infty }^{+\infty }P_{2n-1}(\phi ,\,\partial _{x}\phi ,\,\partial _{x}^{2}\phi ,\,\ldots )\,{\text{d}}x\,}

ここで多項式は再帰的に次のように定義される。 Pn{\displaystyle P_{n}}

P1=ϕ,Pn=dPn1dx+i=1n2PiPn1i for n2.{\displaystyle {\begin{aligned}P_{1}&=\phi ,\\P_{n}&=-{\frac {dP_{n-1}}{dx}}+\sum _{i=1}^{n-2}\,P_{i}\,P_{n-1-i}\quad {\text{ for }}n\geq 2.\end{aligned}}}

運動の最初の積分は次のとおりです。

  • 質量ϕdx,{\displaystyle \int \phi \,\mathrm {d} x,}
  • 勢いϕ2dx,{\displaystyle \int \phi ^{2}\,\mathrm {d} x,}
  • エネルギー。[2ϕ3(xϕ)2]dx{\displaystyle \int \left[2\phi ^{3}-\left(\partial _{x}\phi \right)^{2}\right]\,\mathrm {d} x}

奇数項のみが非自明な(つまりゼロでない)運動積分をもたらす。[ 11 ]P2n+1{\displaystyle P_{2n+1}}

ラックスペア

KdV方程式

tϕ=6ϕxϕx3ϕ{\displaystyle \partial _{t}\phi =6\,\phi \,\partial _{x}\phi -\partial _{x}^{3}\phi }

はラックス方程式として再定式化できる。

Lt=[L,A]LAAL{\displaystyle L_{t}=[L,A]\equiv LA-AL\,}

シュトゥルム・リウヴィル演算子を用いると: L{\displaystyle L}

L=x2+ϕ,A=4x36ϕx3[x,ϕ]{\displaystyle {\begin{aligned}L&=-\partial _{x}^{2}+\phi ,\\A&=4\partial _{x}^{3}-6\phi \,\partial _{x}-3[\partial _{x},\phi ]\end{aligned}}}

ここで はとなる交換子である。[ 12 ] Lax 対はKdV 方程式の最初の積分の無限数を説明する。 [ 13 ][x,ϕ]{\displaystyle [\partial _{x},\phi ]}[x,ϕ]f=fxϕ{\displaystyle [\partial _{x},\phi ]f=f\partial _{x}\phi }

実際、は(定数を無視すれば)時間に依存しないシュレーディンガー作用素であり、そのポテンシャルはである。このラックス定式化により、固有値は実際には に依存しないことがわかる。[ 14 ]L{\displaystyle L}ϕ(x,t){\displaystyle \phi (x,t)}t{\displaystyle t}

ゼロ曲率表現

Lax接続の成分をKdV方程式に 設定することは 、Lax接続のゼロ曲率方程式と同等である。 Lx=(01ϕλ0),Lt=(ϕx2ϕ+4λ2ϕ2ϕxx+2ϕλ4λ2ϕx),{\displaystyle L_{x}={\begin{pmatrix}0&1\\\phi -\lambda &0\end{pmatrix}},L_{t}={\begin{pmatrix}-\phi _{x}&2\phi +4\lambda \\2\phi ^{2}-\phi _{xx}+2\phi \lambda -4\lambda ^{2}&\phi _{x}\end{pmatrix}},}tLxxLt+[Lx,Lt]=0.{\displaystyle \partial _{t}L_{x}-\partial _{x}L_{t}+[L_{x},L_{t}]=0.}

最小作用原理

Korteweg-De Vries 方程式

tϕ+6ϕxϕ+x3ϕ=0,{\displaystyle \partial _{t}\phi +6\phi \,\partial _{x}\phi +\partial _{x}^{3}\phi =0,}

はラグランジュ密度から導かれるオイラー・ラグランジュの運動方程式であり、L{\displaystyle {\mathcal {L}}\,}

定義​ ϕ{\displaystyle \phi }

ϕ:=ψx.{\displaystyle \phi :={\frac {\partial \psi }{\partial x}}.}
オイラー・ラグランジュ方程式の導出

ラグランジアン(式(1))は2次導関数を含むので、この場のオイラー・ラグランジュの運動方程式は

ここで、は成分に関する微分です。 {\displaystyle \partial }μ{\displaystyle \mu }

上の和が暗黙的に表されるので、式(2)は実際には次のように表される。 μ{\displaystyle \mu }

式(1)に式(3)の5つの項を代入して評価する。

tt(L(ttψ))=0{\displaystyle \partial _{tt}\left({\frac {\partial {\mathcal {L}}}{\partial (\partial _{tt}\psi )}}\right)=0}
xx(L(xxψ))=xx(xxψ){\displaystyle \partial _{xx}\left({\frac {\partial {\mathcal {L}}}{\partial (\partial _{xx}\psi )}}\right)=\partial _{xx}\left(-\partial _{xx}\psi \right)}
t(L(tψ))=t(12xψ){\displaystyle \partial _{t}\left({\frac {\partial {\mathcal {L}}}{\partial (\partial _{t}\psi )}}\right)=\partial _{t}\left({\frac {1}{2}}\partial _{x}\psi \right)}
x(L(xψ))=x(12tψ+3(xψ)2){\displaystyle \partial _{x}\left({\frac {\partial {\mathcal {L}}}{\partial (\partial _{x}\psi )}}\right)=\partial _{x}\left({\frac {1}{2}}\partial _{t}\psi +3(\partial _{x}\psi )^{2}\right)\,}
Lψ=0{\displaystyle {\frac {\partial {\mathcal {L}}}{\partial \psi }}=0}

定義を覚えておき、それを使って上記の用語を簡略化すると、 ϕ=xψ{\displaystyle \phi =\partial _{x}\psi }

xx(xxψ)=xxxϕ{\displaystyle \partial _{xx}\left(-\partial _{xx}\psi \right)=-\partial _{xxx}\phi }
t(12xψ)=12tϕ{\displaystyle \partial _{t}\left({\frac {1}{2}}\partial _{x}\psi \right)={\frac {1}{2}}\partial _{t}\phi }
x(12tψ+3(xψ)2)=12tϕ+3x(ϕ)2=12tϕ+6ϕxϕ{\displaystyle \partial _{x}\left({\frac {1}{2}}\partial _{t}\psi +3(\partial _{x}\psi )^{2}\right)={\frac {1}{2}}\partial _{t}\phi +3\partial _{x}(\phi )^{2}={\frac {1}{2}}\partial _{t}\phi +6\phi \partial _{x}\phi }

最後に、これら3つの非ゼロ項を式(3)に代入して、

(xxxϕ)(12tϕ)(12tϕ+6ϕxϕ)=0,{\displaystyle \left(-\partial _{xxx}\phi \right)-\left({\frac {1}{2}}\partial _{t}\phi \right)-\left({\frac {1}{2}}\partial _{t}\phi +6\phi \partial _{x}\phi \right)=0,}

これはまさにKdV方程式である

tϕ+6ϕxϕ+x3ϕ=0.{\displaystyle \partial _{t}\phi +6\phi \,\partial _{x}\phi +\partial _{x}^{3}\phi =0.}

長期漸近解析

十分に速く減衰する滑らかな解は、最終的には右方向に移動するソリトンの有限の重ね合わせと、左方向に減衰する分散成分に分裂することが示される。これはZabuskyとKruskal (1965)によって初めて観察され、振動型リーマン・ヒルベルト問題に対する非線形最急降下法解析を用いて厳密に証明できる。[ 15 ]

歴史

KdV 方程式の歴史は、1834 年のジョン・スコット・ラッセルの実験から始まり、その後 1870 年頃にレイリー卿ジョセフ・ブシネスクによる理論的研究が行われ、最終的に 1895 年にコルテウェグとド・フリースによって研究されました。

KdV方程式はその後、ZabuskyとKruskal (1965) が数値的にその解が長時間にわたって「ソリトン」、すなわち十分に離れた孤立波の集合に分解するように見えることを発見するまで、あまり研究されていませんでした。さらに、ソリトンは互いに通過しても形状にほとんど影響を及ぼさないようです(ただし、これにより位置が変化する可能性はあります)。彼らはまた、 KdV方程式がFPUT系の連続極限であることを示し、フェルミ、パスタ、ウラム、ツィンゴウによる以前の数値実験との関連を示しました。逆散乱変換を用いた解析解の開発は、1967年にガードナー、グリーン、Kruskal、三浦によって行われました。[ 3 ] [ 16 ]

KdV方程式は現在ではホイヘンスの原理と密接に関係していると考えられています。[ 17 ] [ 18 ]

アプリケーションと接続

KdV方程式は物理問題といくつかの関連があります。連続体極限におけるフェルミ-パスタ-ウラム-ツィンゴウ問題における弦の支配方程式であることに加えて、KdV方程式は、以下のような多くの物理的設定における長い1次元波の発展を近似的に記述します。

KdV 方程式は、非線形シュレーディンガー方程式に適用されるような逆散乱変換を使用して解くこともできます。

KdV方程式とグロス・ピタエフスキー方程式

フォームの簡略化されたソリューションを考慮すると

ϕ(x,t)=ϕ(x±t){\displaystyle \phi (x,t)=\phi (x\pm t)}

KdV方程式は次のようになる。

±xϕ+x3ϕ+6ϕxϕ=0{\displaystyle \pm \partial _{x}\phi +\partial _{x}^{3}\phi +6\,\phi \,\partial _{x}\phi =0\,}

または

±xϕ+x(x2ϕ+3ϕ2)=0{\displaystyle \pm \partial _{x}\phi +\partial _{x}(\partial _{x}^{2}\phi +3\phi ^{2})=0\,}

積分し、積分定数がゼロとなる特殊なケースをとると、次のようになります。

x2ϕ3ϕ2=±ϕ{\displaystyle -\partial _{x}^{2}\phi -3\phi ^{2}=\pm \phi \,}

これは一般化定常グロス・ピタエフスキー方程式(GPE) の特別な場合である。λ=1{\displaystyle \lambda =1}

x2ϕ3ϕλϕ=±ϕ{\displaystyle -\partial _{x}^{2}\phi -3\phi ^{\lambda }\phi =\pm \phi \,}

したがって、一般化GPEの特定の解のクラス(真の1次元凝縮体の場合と、 1次元で3次元方程式を用いた場合)では、2つの方程式が1つになる。さらに、負の符号と実数の場合を例にとると、明るいソリトンを生み出す引力的な自己相互作用が得られる。 λ=4{\displaystyle \lambda =4}λ=2{\displaystyle \lambda =2}λ=3{\displaystyle \lambda =3}ϕ{\displaystyle \phi }

バリエーション

KdV方程式には様々なバリエーションが研究されてきました。いくつかを次の表に示します。

名前 方程式
コルテヴェク・デ・フリース(KdV) tu+x3u+6uxu=0{\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u+6u\partial _{x}u=0}
KdV(円筒形) tu+x3u6uxu+12tu=0{\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u-6u\partial _{x}u+{\tfrac {1}{2t}}u=0}
KdV(変形) tu+x(x2u2ηu33u(xu)22(η+u2))=0{\displaystyle \displaystyle \partial _{t}u+\partial _{x}\left({\frac {\partial _{x}^{2}u-2\eta u^{3}-3u(\partial _{x}u)^{2}}{2(\eta +u^{2})}}\right)=0}
KdV(一般化) tu+x3u=x5u{\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u=\partial _{x}^{5}u}
KdV(一般化)tu+x3u+xf(u)=0{\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u+\partial _{x}f(u)=0}
KdV(修正版)tu+x3u±6u2xu=0{\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u\pm 6u^{2}\partial _{x}u=0}
ガードナー方程式tu+x3u(6ε2u2+6u)xu=0{\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u-(6\varepsilon ^{2}u^{2}+6u)\partial _{x}u=0}
KdV(修正版) tu+x3u18(xu)3+(xu)(Aeau+B+Ceau)=0{\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u-{\tfrac {1}{8}}(\partial _{x}u)^{3}+(\partial _{x}u)(Ae^{au}+B+Ce^{-au})=0}
KdV(球状) tu+x3u6uxu+1tu=0{\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u-6u\partial _{x}u+{\tfrac {1}{t}}u=0}
広田・薩摩方程式{ut12uxxx+3uux3(vw)x=0vt+vxxx3uvx=0wt+wxxx3uwx=0{\displaystyle \displaystyle {\begin{cases}u_{t}-{\frac {1}{2}}u_{xxx}+3uu_{x}-3(vw)_{x}=0\\v_{t}+v_{xxx}-3uv_{x}=0\\w_{t}+w_{xxx}-3uw_{x}=0\end{cases}}}
KdV(スーパー){tu=6uxux3u+3wx2wtw=3(xu)w+6uxw4x3w{\displaystyle \displaystyle {\begin{cases}\partial _{t}u=6u\partial _{x}u-\partial _{x}^{3}u+3w\partial _{x}^{2}w\\\partial _{t}w=3(\partial _{x}u)w+6u\partial _{x}w-4\partial _{x}^{3}w\end{cases}}}
KdV(移行型) tu+x3u6f(t)uxu=0{\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u-6f(t)u\partial _{x}u=0}
KdV(変数係数) tu+βtnx3u+αtnuxu=0{\displaystyle \displaystyle \partial _{t}u+\beta t^{n}\partial _{x}^{3}u+\alpha t^{n}u\partial _{x}u=0}
KdV-バーガース方程式tu+μx3u+uxuνx2u=0{\displaystyle \displaystyle \partial _{t}u+\mu \partial _{x}^{3}u+u\partial _{x}u-\nu \partial _{x}^{2}u=0}
KdV-バーガース-フィッシャー方程式(コチャック方程式)tu+μx3u+ϵuxuνx2u=ru(1u){\displaystyle \displaystyle \partial _{t}u+\mu \partial _{x}^{3}u+\epsilon u\partial _{x}u-\nu \partial _{x}^{2}u=ru(1-u)}
川原方程式tu+αuxu+βx3uγx5u=0,{\displaystyle \partial _{t}u+\alpha u\partial _{x}u+\beta \partial _{x}^{3}u-\gamma \,\partial _{x}^{5}u=0,}
非均質KdV tu+αu+βxu+γx2u=Ai(x),u(x,0)=f(x){\displaystyle \partial _{t}u+\alpha u+\beta \partial _{x}u+\gamma \partial _{x}^{2}u=Ai(x),\quad u(x,0)=f(x)}

参照

注記

  1. ^ザブスキー&クラスカル 1965 .
  2. ^ Drazin, PG (1990). 『ソリトン入門』ケンブリッジ、ニューヨーク、メルボルン:ケンブリッジ大学出版局. ISBN 0-521-33655-4
  3. ^ a bガードナーら 1967 .
  4. ^ダリゴル 2005、84ページ。
  5. ^ Korteweg & de Vries 1895 .
  6. ^河原 徹也 (1972). 「分散媒質中の振動孤立波」.日本物理学会誌. 33 (1): 260– 264. Bibcode : 1972JPSJ...33..260K . doi : 10.1143/JPSJ.33.260 .
  7. ^ポリアニンとザイツェフ、2003 年、第 9.1.1 章。 Korteweg-de Vries 方程式。
  8. ^ Vakakis 2002、105–108 ページ。
  9. ^ Dunajski 2009 .
  10. ^三浦、ガードナー、クラスカル 1968 .
  11. ^ディンゲマンス 1997、733ページ。
  12. ^ Polyanin & Zaitsev 2003、第S.10.1章 Laxペア法。
  13. ^ラックス 1968
  14. ^ Dunajski 2009、3​​1–32 ページ。
  15. ^グルネルト&テシュル 2009 .
  16. ^ Dauxois&Peyrard 2006 .
  17. ^チャルブ&ズベリ 2006 .
  18. ^ベレスト&ロウツェンコ 1997 .

参考文献