DOI:10.12677/AAM.2018.712186,PDF,HTML,XML,被引量下载: 1,253浏览: 1,621
作者:陈 静：西南科技大学理学院，四川 绵阳
关键词:广义非线性KdV方程；紧孤子；孤立子；Generalized KdV-Type Equations；Compactons；Solitons

文章引用：陈静. 一类广义色散方程的行波解[J]. 应用数学进展, 2018, 7(12): 1593-1599. https://doi.org/10.12677/AAM.2018.712186

1. 引言

近几十年来，对刻画自然界中非线性波现象的偏微分方程的研究已引起众多学者的高度关注，特别是为平衡弱非线性和色散而产生的孤立子解更是被广泛地研究。为了获得非线性方程的精确解，大量方法已经被用于处理相应的偏微分方程(参见文 [1] - [10] )。在文 [1] 中，作者就运用推广的tanh方法研究了具有任意阶非线性项的广义ZK方程，得到了包含了几种孤子解在内的大量行波解。

通过对非线性色散方程的研究，Rosenau和Hyman在文 [2] 中第一次提出了紧孤子的概念，将具有紧支集的孤立子叫做紧孤子。为了搞清楚非线性色散项在流体中形成的模型，Rosenau和Hyman研究了一类特殊的非线性KdV问题，即如下定义的K(m,n)方程

$u_t+{\left(u^m\right)}_x+{\left(u^m\right)}_{xxx}=0,m>1,1<n\le 3.$(1)

在文 [3] 中，Wazwaz采用Adomian分解法(见文 [11] )求解非线性色散方程。他指出，对于所有m = n的K(m,n)型方程都能用该方法求精确解，为了说明这一方法，他充分地讨论了K(2,2)和K(3,3)两个方程。

Wazwaz在文 [4] 中研究了两个如下形式的(2 + 1)维Boussinesq方程

$u_{tt}-u_{xx}-u_{yy}-a{\left(u^{2n}\right)}_{xx}-b{\left[u^n{\left(u^n\right)}_{xx}\right]}_{xx}=0,n>1,$(2)

$u_{tt}-u_{xx}-u_{yy}-a{\left(u^{-2n}\right)}_{xx}-b{\left[u^{-n}{\left(u^{-n}\right)}_{xx}\right]}_{xx}=0,n>1.$(3)

Wazwaz总结了正余弦拟设法的主要步骤，并用它获得了方程(2) (3)的紧孤子、孤立子、孤立波相似解和周期解。他进一步指出，方程中函数的指数和系数比a/b、以及正负号的变化都会导致解的物理结构产生质的变化。在文 [3] [4] [5] [6] 中，Wazwaz还运用该方法彻底地讨论了一维及更高维的KdV、mKdV、KP方程。

本文的主要工作是寻求另外的方法去扩展Wazwaz在文 [3] 和 [4] 中的工作以获得非线性色散方程的精确解。为了实现这个这一目标，我们将讨论下列具有正负指数的广义KdV型方程

$u_t+a{u}_x+b{\left(u^n\right)}_x+b{\left(u^n\right)}_y+{\left(u^n\right)}_{xxx}+{\left(u^n\right)}_{yyy}=0,n>0,$(4)

$u_t+a{u}_x+b{\left(u^{-n}\right)}_x+b{\left(u^{-n}\right)}_y+{\left(u^{-n}\right)}_{xxx}+{\left(u^{-n}\right)}_{yyy}=0,n>1,$(5)

其中a，b为非零常数。方程(4) (5)可看作是K(n,n)方程添加了新的一项 $a{u}_x$并将它们推广到(2 + 1)维。运用变量 $\xi =\mu x+\eta y-ct$代换，我们把非线性偏微分方程(4) (5)转化成更容易求解的关于 $\xi$的ODES，并因此得到方程(4) (5)行波解的解析表达式。

本文中我们引入了如下q型双曲函数

${\text{sinh}}_q\left(x\right)=\frac{\exp \left(x\right)-q\exp \left(-x\right)}{2},{\text{cosh}}_q\left(x\right)=\frac{\exp \left(x\right)+q\exp \left(-x\right)}{2},$

${\text{csch}}_q\left(x\right)=\frac{1}{{\text{sinh}}_q\left(x\right)},{\text{sech}}_q\left(x\right)=\frac{1}{{\text{cosh}}_q\left(x\right)},$

这里 $q>0$，且为常数。

2. 主要结论

2.1. 具有正指数的KdV型方程

我们首先考虑具有正指数的方程(4)。假设方程(4)具有形如 $u\left(x,y,t\right)=u\left(\xi \right)$的解，其中

$\xi =\mu x+\eta y-ct$， $\mu ,\eta ,c$是不为0的常数，且 $\mu +\eta \ne 0$。易知

$\frac{\partial }{\partial t}=-c\frac{\text{d}}{\text{d}\xi },\frac{\partial }{\partial x}=\mu \frac{\text{d}}{\text{d}\xi },\frac{\partial }{\partial x^3}={\mu }^3\frac{{\text{d}}^3}{\text{d}{\xi }^3},\frac{\partial }{\partial y}=\sigma \frac{\text{d}}{\text{d}\xi },\frac{\partial }{\partial y^3}={\sigma }^3\frac{{\text{d}}^3}{\text{d}{\xi }^3}.$(6)

将(6)式代入方程(4)中得到下面这个非线性ODE

$\left(a\mu -c\right)\frac{\text{d}u}{\text{d}\xi }+b\left(\mu +\eta \right)\frac{\text{d}{u}^n}{\text{d}\xi }+\left({\mu }^3+{\eta }^3\right)\frac{{\text{d}}^3{u}^n}{\text{d}{\xi }^3}=0$(7)

为方便起见，在本文中我们令 $\rho ={\mu }^2-\mu \eta +{\eta }^2$，且对任意的 $\mu \ne 0,\eta \ne 0$都有 $\rho >0$，下面不再赘述。

1) 当 $n=1$时，方程(7)变形为

$\frac{{\text{d}}^{3}u}{\text{d}{\xi }^3}+\frac{1}{\rho }\left(b+\frac{a\mu -c}{\mu +\eta }\right)\frac{\text{d}u}{\text{d}\xi }=0,$(8)

若 $b+\frac{a\mu -c}{\mu +\eta }\ge 0$，解得紧孤子解

$u\left(x,y,t\right)=\{\begin{array}{ll}{\left({\lambda }_1^2+{\lambda }_2^2\right)}^{-\frac{1}{2}}\sin \left(\sqrt{\frac{1}{\rho }\left(b+\frac{a\mu -c}{\mu +\eta }\right)}\xi +\theta \right)+A,\hfill & \left|\phi \right|\le \text{π}\hfill \\ 0,\hfill & 其它\hfill \end{array}$

这里 $\phi =\sqrt{\frac{1}{\rho }\left(b+\frac{a\mu -c}{\mu +\eta }\right)}\xi +\theta$， ${\lambda }_1,{\lambda }_2,A$为任意常数，且当 ${\lambda }_1={\lambda }_2=0$时， $\theta =0$；否则 $\theta$由式子 $\sin \theta ={\lambda }_1{\left({\lambda }_1^2+{\lambda }_2^2\right)}^{-\frac{1}{2}},\cos \theta ={\lambda }_2{\left({\lambda }_1^2+{\lambda }_2^2\right)}^{-\frac{1}{2}}$确定。

若 $b+\frac{a\mu -c}{\mu +\eta }>0$，我们得到下列形式的孤立波相似解(下式中 ${\lambda }_3,{\lambda }_4$都是常数)

$\begin{array}{l}u\left(x,y,t\right)={\lambda }_3{\cosh }_{\left(\frac{{\lambda }_4}{{\lambda }_3}\right)}\left(\sqrt{-\frac{1}{\rho }\left(b+\frac{a\mu -c}{\mu +\eta }\right)}\xi +A\right),当\frac{k_4}{k_3}>0;\\ 或u\left(x,y,t\right)={\lambda }_3{\sinh }_{\left(-\frac{{\lambda }_4}{{\lambda }_3}\right)}\left(\sqrt{-\frac{1}{\rho }\left(b+\frac{a\mu -c}{\mu +\eta }\right)}\xi +A\right),当\frac{k_4}{k_3}<0.\end{array}$

2) 当 $n\ge 2$时，我们考虑下面两种情况：

情况I、 $\mu \ne \frac{c}{a}$时，对方程(7)积分一次并令积分常数为0得

$\left(a\mu -c\right)u+b\left(\mu +\eta \right)u^n+\left({\mu }^3+{\eta }^3\right)\frac{{\text{d}}^2{u}^n}{\text{d}{\xi }^2}=0.$(9)

令 $\frac{\text{d}{u}^n}{\text{d}\xi }=Z,\frac{{\text{d}}^2{u}^n}{\text{d}{\xi }^2}=Z\frac{\text{d}Z}{\text{d}{u}^n}$，代入式(9)中分离变量得

$\left({\mu }^3+{\eta }^3\right)Z\text{d}Z=n\left[\left(c-a\mu \right)u^n-b\left(\mu +\eta \right)u^{2n-1}\right]\text{d}u.$(10)

进一步对方程 积分并令积分常数为0有

$\frac{{\mu }^3+{\eta }^3}{2}{Z}^2=n\left[\frac{c-a\mu }{n+1}{u}^{n+1}-\frac{b\left(\mu +\eta \right)}{2n}{u}^{2n}\right].$(11)

由此有

$\frac{n{u}^{\frac{n-3}{2}}\text{d}u}{\sqrt{\frac{2n\left(c-a\mu \right)}{\left(n+1\right)\left(\mu +\eta \right)}-b{u}^{n-1}}}=\pm \frac{\text{d}\xi }{\sqrt{\rho }}.$(12)

令 $V=\sqrt{\frac{2n\left(c-a\mu \right)}{\left(n+1\right)\left(\mu +\eta \right)}-b{u}^{n-1}}$，

则当 $b>0$时，方程(12)变形为

$\frac{\text{d}V}{\sqrt{\frac{2n\left(c-a\mu \right)}{\left(n+1\right)\left(\mu +\eta \right)}-V^2}}=\pm \frac{n-1}{2n}\sqrt{\frac{b}{\rho }}\text{d}\xi .$(13)

当 $\frac{a\mu -c}{\mu +\eta }\ne 0$时，对式(13)积分解得如下紧孤子解

$u\left(x,y,t\right)=\{\begin{array}{ll}{\left[\frac{2n\left(c-a\mu \right)}{b\left(n+1\right)\left(\mu +\eta \right)}{\cos }^2\left(\frac{n-1}{2n}\sqrt{\frac{b}{\rho }}\xi \pm r\right)\right]}^{\frac{1}{n-1}},\hfill & \left|\varphi \right|\le \frac{\text{π}}{2}\hfill \\ 0,\hfill & 其它\hfill \end{array}$

这里 $\varphi =\frac{n-1}{2n}\sqrt{\frac{b}{\rho }}\xi \pm r$。在本文中，为方便起见我们用r表示任意积分常数。

类似地，在 $b<0$的情况下，我们解得孤立波相似解

$u\left(x,y,t\right)={\left[\frac{2n\left(c-a\mu \right)}{b\left(n+1\right)\left(\mu +\eta \right)}{\cosh }^2\left(\frac{n-1}{2n}\sqrt{-\frac{b}{\rho }}\xi \pm ir\right)\right]}^{\frac{1}{n-1}},i=\sqrt{-1}.$

情况II、 $\mu =\frac{c}{a}$。

方程(7)变形为线性方程

$b\frac{\text{d}{u}^n}{\text{d}\xi }+\left({\mu }^2-\mu \eta +{\eta }^2\right)\frac{{\text{d}}^3{u}^n}{\text{d}{\xi }^3}=0.$(14)

当 $b>0$时，解方程(14)得 $u^n={\left(k_1^2+k_2^2\right)}^{\frac{1}{2}}\text{sin}\left(\sqrt{\frac{b}{\rho }}\xi +\alpha \right)+A$，其中 $k_1,k_2$为任意常数，且当 $k_1=k_2=0$时， $\alpha =0$，否则 $\alpha$满足 $\sin \alpha =k_1{\left(k_1^2+k_2^2\right)}^{-\frac{1}{2}}$， $\cos \alpha =k_2{\left(k_1^2+k_2^2\right)}^{-\frac{1}{2}}$。

由此得到紧孤子解

$u\left(x,y,t\right)=\{\begin{array}{ll}{\left[{\left(k_1^2+k_2^2\right)}^{\frac{1}{2}}\sin \left(\sqrt{\frac{b}{\rho }}\xi +\alpha \right)+A\right]}^{\frac{1}{n}},\hfill & \left|\psi \right|\le \text{π}\hfill \\ 0,\hfill & 其它\hfill \end{array}$

这里 $\psi =\sqrt{\frac{b}{\rho }}\xi +\alpha$。

当 $b<0$，我们有解

$u\left(x,y,t\right)={\left[k_3\exp \left(\sqrt{-\frac{b}{\rho }}\xi \right)+k_4\exp \left(-\sqrt{-\frac{b}{\rho }}\xi \right)+A\right]}^{\frac{1}{n}}$， $k_3,k_4$为任意常数。

综上，我们整理得到下列形式孤立波相似解

$\begin{array}{l}u\left(x,y,t\right)={\left[2k_3{\cosh }_{\left(\frac{k_4}{k_3}\right)}\left(\sqrt{-\frac{b}{\rho }}\xi \right)+A\right]}^{\frac{1}{n}},当\frac{k_4}{k_3}>0;\\ 或u\left(x,y,t\right)={\left[2k_3{\cosh }_{\left(-\frac{k_4}{k_3}\right)}\left(\sqrt{-\frac{b}{\rho }}\xi \right)+A\right]}^{\frac{1}{n}},当\frac{k_4}{k_3}<0.\end{array}$

2.2. 具有负指数的KdV型方程

考虑具有负指数的方程(5)。按照前面所述方法，我们将行波变量 $\xi =\mu x+\eta y-ct$代入方程(5)，使其变为下列ODE

$\left(a\mu -c\right)\frac{\text{d}u}{\text{d}\xi }+b\left(\mu +\eta \right)\frac{\text{d}{u}^{-n}}{\text{d}\xi }+\left({\mu }^3+{\eta }^3\right)\frac{{\text{d}}^3{u}^{-n}}{\text{d}{\xi }^3}=0,$(15)

其中 $\mu ,\eta ,c$为不为0的常数，且 $\mu ,\eta ,c$。

情况1、 $\mu \ne \frac{c}{a}$。

对方程(15)积分并令积分常数为0可得

$\left(a\mu -c\right)u+b\left(\mu +\eta \right)u^{-n}+\left({\mu }^3+{\eta }^3\right)\frac{{\text{d}}^2{u}^{-n}}{\text{d}{\xi }^2}=0.$(16)

进一步变形有

$\frac{-n{u}^{\frac{-n-\text{3}}{2}}\text{d}u}{\sqrt{\frac{2n\left(c-a\mu \right)}{\left(n-1\right)\left(\mu +\eta \right)}-b{u}^{-n-1}}}=\pm \frac{\text{d}\xi }{\sqrt{\rho }},$(17)

当 $b>0$时，解得周期解

$u\left(x,y,t\right)={\left[\frac{b\left(n-1\right)\left(\mu +\eta \right)}{2n\left(c-a\mu \right)}{\sec }^2\left(\frac{n+1}{2n}\sqrt{\frac{b}{\rho }}\xi \pm r\right)\right]}^{\frac{1}{n+1}}.$

当 $b<0$时，解得孤立子解

$u\left(x,y,t\right)={\left[\frac{b\left(n-1\right)\left(\mu +\eta \right)}{2n\left(c-a\mu \right)}{\text{sech}}^2\left(\frac{n+1}{2n}\sqrt{-\frac{b}{\rho }}\xi \pm ir\right)\right]}^{\frac{1}{n+1}}.$

情况2、 $\mu =\frac{c}{a}$。

此时我们考虑如下方程

$b\frac{\text{d}{u}^{-n}}{\text{d}\xi }+\left({\mu }^2-\mu \eta +{\eta }^2\right)\frac{{\text{d}}^3{u}^{-n}}{\text{d}{\xi }^3}=0.$(18)

当 $b>0$时，解之得周期解

$u\left(x,y,t\right)={\left[{\left(l_1^2+l_2^2\right)}^{\frac{1}{2}}\sin \left(\sqrt{\frac{b}{\rho }}\xi +\beta \right)+A\right]}^{-\frac{1}{n}},$

这里 $l_1,l_2,A$为常数， $l_1^2+l_2^2\ne 0$或 $A\ne \text{0}$，且 $\beta$满足等式 $\sin \beta =l_1{\left(l_1^2+l_2^2\right)}^{-\frac{1}{2}}$， $\cos \beta =l_2{\left(l_1^2+l_2^2\right)}^{-\frac{1}{2}}$。

当 $b<0$时，解之得一般解

$u\left(x,y,t\right)={\left[l_3\exp \left(\sqrt{-\frac{b}{\rho }}\xi \right)+l_4\exp \left(-\sqrt{-\frac{b}{\rho }}\xi \right)+A\right]}^{-\frac{1}{n}}.$

当 $A=0$时则有两个特殊的孤立子解

$\begin{array}{l}u\left(x,y,t\right)={\left[\frac{1}{2l_3}{\text{sech}}_{\left(\frac{l_4}{l_3}\right)}\left(\sqrt{-\frac{b}{\rho }}\xi \right)\right]}^{\frac{1}{n}},当\frac{l_4}{l_3}>0;\\ 或u\left(x,y,t\right)={\left[\frac{1}{2l_3}\csc h_{\left(-\frac{l_4}{l_3}\right)}\left(\sqrt{-\frac{b}{\rho }}\xi \right)\right]}^{\frac{1}{n}},当\frac{l_4}{l_3}<0.\end{array}$

3. 结论

本文的主要目的是研究广义非线性具有正负n指数的KDV型方程(4)和(5)的精确解。我们运用变量代换法获得了包括紧孤子、孤立子、孤立波相似解和周期解在内的行波解。在下表1中，我们容易观察

到不同的指数和系数决定了具有不同物理结构的解。

研究偏微分方程以获得其精确解进而研究非线性波现象的运动规律和特点是众多科学工作者一直以来的重要工作，值得指出的是,文中所使用的方法对于获得这些非线性偏微分方程的一些特殊解有一定的优势。

参考文献

[1]	Li, B., Chen, Y. and Zhang, H.Q. (2003) Exact Travelling Wave Solutions for a Generalized Zakharov-Kuznetsov Equation. Applied Mathematics and Computation, 146, 653-666. https://doi.org/10.1016/S0096-3003(02)00610-0
[2]	Rosenau, P. and Hyman, J.M. (1993) Compactons: Solitons with Finite Wavelength. Physical Review Letters, 70, 564-567. https://doi.org/10.1103/PhysRevLett.70.564
[3]	Wazwaz, A.M. (2001) A Study of Nonlinear Dispersive Equa-tions with Solitary-Wave Solutions Having Compact Support. Mathematics and Computers in Simulation, 56, 269-276. https://doi.org/10.1016/S0378-4754(01)00291-9
[4]	Wazwwaz, A.M. (2005) Variants of the Two-Dimensional Boussinesq Equation with Compactons, Solitons, and Periodic Solutions. Computers and Mathematics with Applications, 49, 295-301. https://doi.org/10.1016/j.camwa.2004.06.029
[5]	Wazwwaz, A.M. (2003) A Study on Nonlinear Dispersive Par-tial Differenial Equations of Compact and Nocompact Solutions, and Nocompact Structures. Applied Mathematics and Computation, 135, 399-409. https://doi.org/10.1016/S0096-3003(02)00005-X
[6]	Wazwwaz, A.M. (2004) Dinstinct Variants of the KdV Equation with Compact and Nocompact Structures. Applied Mathematics and Computation, 150, 365-377. https://doi.org/10.1016/S0096-3003(03)00238-8
[7]	Kivshar, Y. (1994) Compactons in Dispersive Lattices. Nonlinear Coherent Structures in Physics and Biology, 329, 255-258. https://doi.org/10.1007/978-1-4899-1343-2_38
[8]	Dinda, P.T. and Remoissenet, M. (1999) Breather Compactons in Nonlinear Klein-Gordon Systems. Physical Review E, 60, 6218-6221. https://doi.org/10.1103/PhysRevE.60.6218
[9]	Ismail, M.S. and Taha, T. (1998) A Numerical Study of Com-pactons. Mathematics and Computers in Simulation, 47, 519-530.
[10]	Rosenau, P. (2000) Compact and Nocompact Dispersive Structures. Physics Letters A, 275, 193-203. https://doi.org/10.1016/S0375-9601(00)00577-6
[11]	Adomian, G. (1994) Solving Frontier Problems of Physics: The Decomposition Method. Kluwer Academic Publishers, Boston. https://doi.org/10.1007/978-94-015-8289-6