KdV方程式

数学の解析学ないし物理学の流体力学関連分野に於ける「KdV方程式」とはソリトン解を持つような非線形の偏微分方程式である。名称はオランダの数学者コルトヴェーグとド・フリースに因む。

概要

KdV方程式は以下の如き3階非線形偏微分方程式で表される。；

• $\frac{\partial u}{\partial t}+6u\frac{\partial u}{\partial x}+\frac{{\partial }^{3}u}{\partial x^3}=0$

上記方程式は略式に

• $u_t+6u{u}_x+u_{xxx}=0$

や

• ${\partial }_{t}u+6u{\partial }_{x}u+\underset{x}{\overset{3}{\partial }}u=0$

などと書かれる事もある。

解法

この節では ちょっぴり大変な計算になるけどKdV方程式の一番簡単な解(1ソリトン解)を導いてみる事にしよう。まず上記方程式の解が

$u(x,t)=\varphi (x-ct)$

とゆー性質(進行波解☆)を満たすとしよう。これを偏微分したら

${\partial }_{t}u={\varphi }^{\prime }\cdot (-c),{\partial }_{x}u={\varphi }^{\prime }\cdot 1,\underset{x}{\overset{3}{\partial }}u={\varphi }^{‴}$

となるが、これらをKdV方程式に代入すれば次式が成り立つ。；

$-c{\varphi }^{\prime }+6\varphi {\varphi }^{\prime }+{\varphi }^{‴}=0$

ここで合成関数の微分法の公式

$\{(y^n){\}}^{\prime }=n{y}^{n-1}\cdot y^{\prime }$

から

$(-c\varphi +3{\varphi }^2+{\varphi }^{″}{)}^{\prime }=0$

が得られる。ここで極限

$\varphi ,{\varphi }^{\prime },{\varphi }^{″}\to 0,(x\to \pm \mathrm{\infty })$

を仮定して上述の方程式をxから+∞まで積分したら

$\begin{array}{rl}\underset{x}{\overset{\mathrm{\infty }}{\int }}(-c\varphi +3{\varphi }^2+{\varphi }^{″}{)}^{\prime }dx& =[-c\varphi +3{\varphi }^2+{\varphi }^{″}{]}_x^{\mathrm{\infty }}\\ & =c\varphi -3{\varphi }^2-{\varphi }^{″}\\ & =0\end{array}$

となる。ここで両辺×φ'計算したら

$c\varphi {\varphi }^{\prime }-3{\varphi }^2{\varphi }^{\prime }-{\varphi }^{\prime }{\varphi }^{″}=(\frac{1}{2}c{\varphi }^2-{\varphi }^3-\frac{1}{2}\varphi {\text{'}}^2)=0$

を得る。この両辺を零からxまで積分すれば

$c{\varphi }^2/2-{\varphi }^3-\varphi {\text{'}}^2/2=0$

$\therefore c{\varphi }^2-2{\varphi }^3-\varphi {\text{'}}^2=0$

が求まる。これをφ'に関して解けば

${\varphi }^{\prime }=\pm \sqrt{c}\varphi \sqrt{1-2\varphi /c}$

となる。ここでφ'(0)=0とすれば簡単な計算により

$\varphi (0)=c/2$

が得られる(初期条件みたいなもんです☆)。

さて、ここで関数φは元々変数(x,t)の2変数関数であったが合成関数の偏微分法により或る1変数関数と見做しても差し支え無い事が分かるであろう。従って便宜上以下ではしばらくの間φをφ=φ(x)と書く事にする。

ゆえに以下の如き変数分離形の1階常微分方程式が得られた事になる。；

$\frac{d\varphi }{dx}=\pm \sqrt{c}\varphi \sqrt{1-\frac{2\varphi }{c}},\varphi (0)=\frac{c}{2}$

で、これを

$\frac{1}{\varphi \sqrt{1-2\varphi /c}}\frac{d\varphi }{dx}=\pm \sqrt{c}$

と変形して両辺を零からxまで積分すれば

$\underset{0}{\overset{x}{\int }}\frac{1}{\varphi \sqrt{1-2\varphi /c}}\frac{d\varphi }{dx}dx=\pm \sqrt{c}x$

が成り立つ。ここで

$y=\varphi (x),dy/dx={\varphi }^{\prime }(x)$

とおけば上記積分は

$\pm \sqrt{c}x=\underset{\varphi (0)}{\overset{\varphi (x)}{\int }}\frac{1}{y\sqrt{1-2y/c}}dy$

と変形できる(定積分の置換積分公式より)。

で、更に

$y=z^{-1},dy/dz=-z^{-2}$

と置くと再び置換積分により

$\begin{array}{rl}\pm \sqrt{c}x& =\underset{1/\varphi (0)}{\overset{1/\varphi (x)}{\int }}\frac{1}{z^{-1}\sqrt{1-2z^{-1}/c}}(-z^{-2})dz\\ & =-\underset{1/\varphi (0)}{\overset{1/\varphi (x)}{\int }}\frac{1}{\sqrt{z^2-2z/c}}dz\\ & =-\underset{2/c}{\overset{1/\varphi (x)}{\int }}\frac{1}{\sqrt{(z-1/c{)}^2-1/c^2}}dz\end{array}$

となる。(※√の中を平方完成した。また定積分の下端に1/φ(0)=2/cを適用した。)

次に

$\zeta =z-1/c,d\zeta /dz=1$

と置いて再度置換積分したら(※上端・下端の値の変化に注意！) $\begin{array}{rl}\pm \sqrt{c}x& =-\underset{1/c}{\overset{1/\varphi (x)-1/c}{\int }}\frac{1}{\sqrt{{\zeta }^2-1/c^2}}d\zeta \\ & =-{[\log (\zeta +\sqrt{{\zeta }^2-\frac{1}{c^2}})]}_{1/c}^{1/\varphi (x)-1/c}\\ & =\log \frac{1}{c}-\log (\frac{1}{\varphi }-\frac{1}{c}+\sqrt{{(\frac{1}{\varphi }-\frac{1}{c})}^2-\frac{1}{c^2}})\\ & =-\log (\frac{c}{\varphi }-1+\sqrt{{(\frac{c}{\varphi }-1)}^2-1})\end{array}$

$\therefore e^{\mp \sqrt{c}x}=\frac{c}{\varphi }-1+\sqrt{{(\frac{c}{\varphi }-1)}^2-1}$

が得られる。この辺の式変形なかなか難しいんだけど途中で以下の公式使ってるんでこれ見て頑張って計算してみてね☆

• $\int \frac{1}{\sqrt{x^2-a^2}}dx=\log (x+\sqrt{x^2-a^2})$
• $\log M+\log N=\log MN$
• $t=\log x\iff e^t=x$

書くのが大変なんで以下ψ=c/φ-1と略記する事にしよう。すると

$e^{\mp \sqrt{c}x}=\psi +\sqrt{{\psi }^2-1}$

$\therefore e^{\mp \sqrt{c}x}-\psi =\sqrt{{\psi }^2-1}$

となる。上式の両辺を2乗したら

$e^{\mp 2\sqrt{c}x}-2e^{\mp \sqrt{c}x}\psi +{\psi }^2={\psi }^2-1$

より以下の数式が成り立つ。；

$2e^{\mp \sqrt{c}x}\psi =1+e^{\mp 2\sqrt{c}x}$

$\begin{array}{rl}\therefore \psi & =\frac{1+e^{\mp 2\sqrt{c}x}}{2e^{\mp \sqrt{c}x}}=\frac{e^{\pm \sqrt{c}x}(1+e^{\mp 2\sqrt{c}x})}{2}\\ & =\frac{e^{\pm \sqrt{c}x}+e^{\mp \sqrt{c}x}}{2}=\cosh (\sqrt{c}x)\end{array}$

(※ちなみに上記のcoshは双曲線関数の1つ(双曲線余弦関数)である。)

従って

$\frac{c}{\varphi }=1+\cosh (\sqrt{c}x)$

が得られる。

ここで双曲線関数についても数多ある三角関数の公式と類似の公式が知られていて以下の如き半角公式(もどき)みたいな等式が成立する。；

${\cosh }^2\left(\frac{\gamma }{2}\right)=\frac{1+\cosh \gamma }{2}$

これにより

$\frac{c}{\varphi }=2{\cosh }^2\left(\frac{\sqrt{c}}{2}x\right)$

$\varphi =\frac{c}{2{\cosh }^2(\sqrt{c}x/2)}=\frac{c}{2}{\mathrm{sech}}^2\left(\frac{\sqrt{c}}{2}x\right)$

が言える。(※ただしsechθ=1/coshθ)

ここで$\kappa =\sqrt{c}/2$とおけば

${\kappa }^2=c/4,2{\kappa }^2=c/2$

と書けるから件の関数φは

$\varphi (x)=2{\kappa }^2{\mathrm{sech}}^2(\kappa x)$

と表される事になる。で、これを元の進行波解の関数形

$u(x,t)=\varphi (x-ct)$

に戻したら

• $u(x,t)=2{\kappa }^2{\mathrm{sech}}^2(\kappa (x-ct))$

となって遂に漸くやっと晴れてKdV方程式の解が導かれた事になる(☆お疲れ様でした♪)

…で、このu=u(x,t)がKdV方程式を満たす事を証明したいと思った時はこの関数をそのまま偏微分するのも良いんだけど、以下の双曲線関数の性質＆微分法

$1-{\tanh }^2\theta =\frac{1}{{\cosh }^2\theta }(={\mathrm{sech}}^2\theta )$

$(\tanh \theta {)}^{\prime }=1-{\tanh }^{}\theta$

を使って(双曲線正接関数の微分法で)検算した方が簡単かも知んない★

関連項目

• 変形KdV方程式
• KP方程式
• サインゴルドン方程式
• 非線形シュレディンガー方程式
• ブジネ方程式
• ベンジャミン–オノ方程式
• 戸田格子方程式