卡東穆塞夫-彼得韋亞斯維利方程

卡東穆塞夫-彼得韋亞斯維利方程(Kadomtsev-Petviashvili equation)，簡稱KP方程，是1970年蘇聯物理學家波里斯·卡東穆塞夫 和弗拉基米爾-彼得韋亞斯維利創立以模擬非線性波動的非線性偏微分方程^[1]：

${\displaystyle \displaystyle \partial _{x}(\partial _{t}u+u\partial _{x}u+\epsilon ^{2}\partial _{xxx}u)+\lambda \partial _{yy}u=0}$

這張在法國大西洋岸雷島（RHE）鯨魚燈塔拍攝的照片，顯示淺海上田字形的橢圓餘弦波列。這種淺水中的孤波可以由卡東穆夫-彼得韋亞斯維利方程模擬。

其中 ${\displaystyle \lambda =\pm 1}$.

解析解

卡東穆塞夫-彼得韋亞斯維利方程有解析解^[2]

行波解

${\displaystyle u(x,y,t)=C5+12.*_{C}2*\tanh(_{C}1+_{C}2*x+_{C}3*y-(.50000000000000000000*(8.*_{C}2^{4}+_{C}3^{2}))*t/_{C}2)}$

代人參數： C5 = 1, _C1 = 0, _C2 = 1, _C3 = 3

得：${\displaystyle u=1+12.*tanh(x+3*y-8.5000000000000000000*t)}$

• 
  Kadomtsev-Petviahivili equation 3D plot
• 
  Kadomtsev-Petviashivili pde animation

Sech 函數亮孤立子解

利用sech函數展開法可得卡東穆塞夫-彼得韋亞斯維利方程的sech函數解和tanh函數解^[3]。

${\displaystyle u:=a*sech(a*x+b*y+c*z-(a^{4}+3*b^{2}+3*c^{2})/a)*t}$

參數：a = -2 .. 2, b = -2 .. 2, c = 0

• 
  Kadomtsev Petviashivili pde sech solution 3d plot
• 
  卡東穆塞夫-彼得韋亞斯維利方程sech函數亮孤立子

tanh 函數解

${\displaystyle u:=2*a^{2}*tanh(a*x+b*y+(8*a^{4}-3*b^{2})/a)^{2}*t}$^[4]。

參數：a = 2, b = -2；

• 
  KP方程暗tanh函數暗孤立子

雅可比橢圓函數解

通過朗斯基行列式展開法可得卡東塞穆夫-彼得韋亞斯維利方程多個雅可比橢圓函數解^[5]。

${\displaystyle u4:={\frac {(-4*m^{2}*k[1]^{2}*g)}{({\sqrt {1-m^{2}}}*sn(\xi [1],k)*sin(\xi [2])+dn(\xi [1],k)*cos(\xi [2])*cn(\xi [1],k))^{2})}}}$

其中：

${\displaystyle g=(m^{2}-1)*sn(\xi [1],k)^{2}+(2-2*m^{2})*sn(\xi [1],k)^{4}+cos(\xi [2])^{2};-2*sn(\xi [1],k)^{2}*cos(\xi [2])^{2}+m^{2}*sn(\xi [1],k)^{4}*cos(\xi [2])^{2}}$

${\displaystyle \xi [1]=k[1]*x+\lambda [1]*y+(4*m^{2}*k[1]^{3}+16*k[1]^{3}-3*\sigma ^{2}*\lambda [1]^{2}/k[1])*t+\gamma [1]}$

${\displaystyle \xi [2]={\sqrt {1-m^{2}}}*(k[1]*x+\lambda [1]*y+(-4*m^{2}*k[1]^{3}+16*k[1]^{3}-3*\sigma ^{2}*\lambda [1]^{2}/k[1])*t)-\gamma [2]}$

代入後得：

${\displaystyle f4:=-4*m^{2}*k[1]^{2}*((m^{2}-1)*JacobiSN(k[1]*x+\lambda [1]*y+(4*m^{2}*k[1]^{3}+16*k[1]^{3}}$

${\displaystyle -3*\sigma ^{2}*\lambda [1]^{2}/k[1])*t+\gamma [1],k)^{2}+(2-2*m^{2})*JacobiSN(k[1]*x+\lambda [1]*y}$

${\displaystyle +(4*m^{2}*k[1]^{3}+16*k[1]^{3}-3*\sigma ^{2}*\lambda [1]^{2}/k[1])*t+\gamma [1],k)^{4}+}$

${\displaystyle cos({\sqrt {(1-m^{2})}}*(k[1]*x+\lambda [1]*y+(-4*m^{2}*k[1]^{3}+16*k[1]^{3}-3*\sigma ^{2}*\lambda [1]^{2}/k[1])*t)}$

${\displaystyle -\gamma [2])^{2})/(sqrt(1-m^{2})*JacobiSN(k[1]*x+\lambda [1]*y+(4*m^{2}*k[1]^{3}+16*k[1]^{3}-}$

${\displaystyle 3*\sigma ^{2}*\lambda [1]^{2}/k[1])*t+\gamma [1],k)*sin({\sqrt {(}}1-m^{2})*(k[1]*x+\lambda [1]*y+}$

${\displaystyle (-4*m^{2}*k[1]^{3}+16*k[1]^{3}-3*\sigma ^{2}*\lambda [1]^{2}/k[1])*t)-\gamma [2])+JacobiDN(k[1]*x+\lambda [1]*y+}$

${\displaystyle (4*m^{2}*k[1]^{3}+16*k[1]^{3}-3*\sigma ^{2}*\lambda [1]^{2}/k[1])*t+\gamma [1],k)*cos({\sqrt {1-m^{2}}}*(k[1]*x+\lambda [1]*y+}$

${\displaystyle (-4*m^{2}*k[1]^{3}+16*k[1]^{3}-3*\sigma ^{2}*\lambda [1]^{2}/k[1])*t)-\gamma [2])*JacobiCN(k[1]*x+}$

${\displaystyle \lambda [1]*y+(4*m^{2}*k[1]^{3}+16*k[1]^{3}-3*\sigma ^{2}*\lambda [1]^{2}/k[1])*t+\gamma [1],k))^{2}}$

• 
  Kadomtsev Petviashivili pde elliptic function solution 3d plot