KdV方程式(KdVほうていしき、)、もしくはコルトヴェーグ・ドフリース方程式とは、非線形波動を記述する非線形偏微分方程式の一つである。ソリトン解を有する可積分系の代表的な例として知られる。方程式の名前は、定式化を行った (D. Korteweg) と (G. de Vries) に因む。
概要 時間変数 と空間変数 をもつ一次元実数値関数 に対して、 をゼロではない任意の実定数として
で与えられる非線形偏微分方程式をKdV方程式というコルテヴェク‐ドフリースの方程式 - コトバンク(法則の事典)。また、各変数と に適当なを施せば、係数を と取りなおすことができる。このとき、各変数に対する偏微分を右下の添え字として表せば、
となる。
KdV方程式は多様な非線形波動、特に孤立波と呼ばれる大振幅波動を説明する偏微分方程式として知られる。例えば水深が波長と比べて十分小さい浅水波はKdV方程式で記述でき、津波はKdVソリトン的な性質を示す。他に、プラズマにはイオン音波や磁気音波などの多様な波動が存在するが、その大振幅波に対応したKdV波動が理論的に知られている。
KdV方程式の一般的な解法としては、や広田の直接法が存在する。
非線形項・分散項 KdV方程式の第二項 を波の立上りの効果を表す非線形項、第三項 を波の広がりの効果を表す分散項という。ここでの分散とは波数空間での分散を意味しており、波動がより幅広い波数を持つ多数の波に分散していくことを現わしている。KdV方程式が支配する物理系では、非線形効果による波の成長と波の分散が釣り合うことで、孤立波(ソリトン)が形成されて形を崩すことなく伝播する。
その他の表示 KdV方程式の係数のとり方はいくつかの流儀が存在するが、いずれも適当な変数変換の下で、互いに移り変われる。例えば、 なる変換による
や なる変換による
もよく用いられる。
歴史的背景 KdV方程式の研究の歴史は、1834年、造船技師がエジンバラ郊外の運河で孤立波を観察したことに遡るWeisstein。彼は艀用の狭い運河のそばで馬に乗っているときに、孤立波が運河を伝播する様子を偶然、目撃した。彼はその時の状況を次のように記している。
造船技師かつ流体力学の研究者であった彼は、実験用の水槽を作り、研究を進めた。そして、次の結果を得た。 浅水波の伝播において、孤立波(永久型の長波)が存在する。 一定水深の水路において、孤立波の速度は で与えられる。ここで は重力加速度、 は静止した状態の流体水面から測った波の高さ、 は静止流体の深さである。この研究結果は、ケルヴィン卿、ストークス、、エアリーといった当時の科学者の間で、孤立波の存在についての大きな論争を起こした。特にエアリーは孤立波の存在に否定的であった。孤立波の存在に最終決着がつくのは、スコット・ラッセルの観測から60年経った後であった。
1895年、 (D. Korteweg) と (G. de Vries) とは適度に振幅の小さい浅水波を記述する方程式として、KdV方程式を導いた。この方程式は孤立波を含む永久波解を持っていた。但し、彼らの研究は注目を浴びることなく、長い間、忘れ去られていた。
right|thumb|ザブスキーとクルスカルの報告結果を基に同条件で計算したKdV方程式の解の時間発展の様子。時刻 での初期条件として与えられた余弦波は、時間とともに形を変え、いくつかの孤立波の集まりとなる。
70年後、KdV方程式はとによる、非線形格子におけるエネルギー伝播の問題(フェルミ・パスタ・ウラムの問題)の研究過程で再発見された。1965年、彼らは非線形格子の連続体モデルの数値計算において、不思議な現象を見い出した。一つは、余弦波で与えた初期状態がいくつかの孤立波に分裂する現象であり、もう一つは二つの孤立波の伝播において、速度が速い孤立波が速度が遅い孤立波を追い越す形で衝突してもそれぞれの波形が壊れず、そのまま伝播する現象である。彼らはこうした粒子性を有する波動現象を、孤立波と粒子を表す接尾語 -on を合わせ、ソリトンと名付けた。このザブスキーとクルスカルによる研究を契機に、こうした可積分系の性質は注目を集め、その後の研究活発化と理論の発展につながった。
KdV方程式の解 KdV方程式の解として、次のものが存在する。
1-ソリトン解 一つの孤立波を表す1ソリトン解は次の形で与えられる。
ここで、 は で与えられる双曲線正割関数を表す。速度での定常伝播を仮定して、KdV方程式を常微分方程式に還元して解くことで導出することができる。この解は という関数形を有しており、 で表される一つのピークを持つ孤立波が形を保ったまま、速度 で伝播する状況に対応している。また、振幅値 は速度 に比例しており、波の高さ(振幅)が高いほど、速く伝播する性質を持つ。
特筆すべきは孤立波の振幅、空間的な幅、速度がすべてないしの1つの自由度に制約されることである。これを満たさない波形を初期条件として与えると、波は条件を満たすソリトンの列に分裂する。
2-ソリトン解 二つの孤立波を表す2ソリトン解は次の形で与えられる。
但し、
である。
この解は次のような行列式による表示を行うことも可能である。
-ソリトン解 個の孤立波を表す ソリトン解は次の形で与えられる。
ここで、 は 次の正方行列で、その 行 列成分 が
で与えられる。但し、 はクロネッカーのデルタを表す。
周期解(クノイダル波) KdV方程式はヤコビの楕円関数 (クノイダル関数)で表される周期解
をもつ。ただし、
である。
保存量 可積分系であるKdV方程式は、時間に対して不変となる無限個の独立な保存量を持つという著しい性質を持つ。1968年、日系人数学者らによって、この性質は見出された。KdV方程式
については、
が保存量となる。ソリトンはソリトン間の相互作用や分裂といった複雑な現象を示すが、これらの保存量を用いることでソリトンの挙動を初期値問題として予測することができる。
脚注
参考文献 論文
書籍
関連項目 可積分系 ソリトン KP方程式 (2次元版KdV方程式) Burgers方程式(KdV方程式と並んで非線形波動を表す偏微分方程式。2階微分を含み衝撃波を表す。)
外部リンク