DOI:10.12677/AAM.2019.83056,PDF,HTML,XML,被引量下载: 887浏览: 1,204
作者:林 桐：华南理工大学数学学院，广东 广州
关键词:行波系统；分支曲线；数值模拟；分支相图；The Traveling Wave System；Branch Curve；Numerical Simulation；Branch Phase Diagram

文章引用：林桐. 一类广义KdV行波系统的相图研究[J]. 应用数学进展, 2019, 8(3): 503-511. https://doi.org/10.12677/AAM.2019.83056

1. 引言

对非线性微分方程的研究，比线性微分方程的研究难度大得多，大多只能依靠数值解法。数学家们提出反散射法、达布变换法，tanh、雅可比函数展开法等构造性解法，是非线性偏微分方程研究的福音。目前已知的非线性偏微分方程，大约只有3000多种，而有名的，也不过一百多种。多以发现它们的学者名字命名，例如，LAX可积系统下的有：BBM方程、Belousov-Zhabotinsky方程、本杰明–小野方程、Bogoyavlenski-Konoplechenko方程、波恩–英菲尔德方程、Boussinesq方程、博格斯方程、KdV方程、MKdV方程、五阶KdV方程、五阶色散KdV方程等。

关于KdV方程的研究，已有许多著名的结果，B. Dey研究精确的哈密顿密度和守恒定律，给出了方程的两个扭解 [1] [2] 。J. F. Zhang等人使用两种不一样的方法给出了方程的孤立波解和奇异波解，他们针对KdV和mKdV方程，提出了两种构造新的精确孤波解的方法，并且总结出12种类型的孤子解，归结为mKdV方程的孤子解 [3] [4] 。J. Yu采用齐次平衡法得到了一个扭解 [5] 。R. Grimshaw等人为方程研究了大幅度孤子 [6] 。E. G. Fan使用代数方法给出了方程的钟形孤立子、扭形孤立子和雅克比周期解 [7] ，继而改进了这个方法，又构造了一般非线性方程的一系列显式精确解 [8] 。M. Y. Tang运用动力系统分支法给方程研究了孤立波和它们的分支 [9] 。Y. Z. Peng用改进的隐射方法得到了一些包含双曲型函数的孤立波解、包含雅各比椭圆函数的周期波解和三角函数的周期波解 [10] 。刘正荣将研究微分方程数值解的基本方法编撰在 [11] 中，宋明运用了这本书的方法，研究了广义KP-BBM方程的周期波解 [12] 。

本文研究的方程形式如下：

$u_t+\left(a+bu\right)u{u}_x+u_{xxx}=0$(1)

我们将讨论方程在由参数a和b构成的平面上的分支相图。

2. 预备知识

2.1. 行波系统的奇点

假设我们有系统：

$\begin{array}{l}\frac{\text{d}x}{\text{d}t}=p\left(x,y\right)\\ \frac{\text{d}y}{\text{d}t}=q\left(x,y\right)\end{array}$(2)

使得 $p\left(x,y\right)=0,q\left(x,y\right)=0$的特殊点 $\left(x_0,y_0\right)$称为奇点，计算雅各比矩阵

$J=\left(\begin{array}{cc}{p}_x& p_y\\ q_x& q_y\end{array}\right)$

令行列式 $\left|J-\lambda I\right|=0$的 $\lambda$就是该奇点的特征值。当奇点的特征值都不为零时，对应的情况如下：

① 如果 ${\lambda }_1$和 ${\lambda }_2$都是实数，且 ${\lambda }_1$不等于 ${\lambda }_2$， ${\lambda }_1{\lambda }_2>0$，那么 $\left(x_0,y_0\right)$是一个结点，并且， ${\lambda }_1$和 ${\lambda }_2$都小于0时是稳定结点，都大于0时是非稳定结点；

② 如果 ${\lambda }_1$和 ${\lambda }_2$都是实数，且 ${\lambda }_1{\lambda }_2<0$，那么 $\left(x_0,y_0\right)$是一个鞍点；

③ 如果 ${\lambda }_1={\lambda }_2=\lambda$是实数，且 $p_y\left(x_0,y_0\right)$和 $q_x\left(x_0,y_0\right)$中至少有一个非零，那么 $\left(x_0,y_0\right)$是一个退化结点，当 $\lambda <0$时是稳定的， $\lambda >0$时是不稳定的；

④ 如果 ${\lambda }_1={\lambda }_2=\lambda$是实数，且 $p_y\left(x_0,y_0\right)=q_x\left(x_0,y_0\right)=0$，那么 $\left(x_0,y_0\right)$是一个临界结点，当 $\lambda <0$时是稳定的， $\lambda >0$时是不稳定的；

⑤ 如果 ${\lambda }_{1,2}=\alpha \pm i\beta$， $\alpha \ne 0$，那么 $\left(x_0,y_0\right)$是一个焦点，当 $\alpha <0$时是稳定的， $\alpha >0$时是不稳定的；

⑥ 如果 ${\lambda }_{1,2}=\pm i\beta$，并且可以找到

$\begin{array}{l}\frac{\text{d}x}{\text{d}t}=p\left(x,y\right)\\ \frac{\text{d}y}{\text{d}t}=q\left(x,y\right)\end{array}$

的首次积分，那么 $\left(x_0,y_0\right)$是一个中心。

假设方程(1)存在行波解 $\xi =x-ct$，代入 $u=\phi \left(\xi \right)$，得到

$u_t=-c{\phi }^{\prime }(\xi )$

$u_x={\phi }^{\prime }(\xi )$

$u_{xxx}={\phi }^{‴}(\xi )$

代入方程(1)，得到

$-c{\phi }^{\prime }\left(\xi \right)+a\phi \left(\xi \right){\phi }^{\prime }\left(\xi \right)+b{\phi }^2\left(\xi \right){\phi }^{\prime }\left(\xi \right)+{\phi }^{‴}\left(\xi \right)=0$

对上式进行积分，并且取积分常数为零，得到

$-c\phi \left(\xi \right)+\frac{a}{2}{\phi }^2\left(\xi \right)+\frac{b}{3}{\phi }^3\left(\xi \right)+{\phi }^{″}\left(\xi \right)=0$

下面，令 ${\phi }^{\prime }\left(\xi \right)=y$

代入方程，得到微分方程组

$\begin{array}{l}\frac{\text{d}\phi }{\text{d}\xi }=y\\ \frac{\text{d}y}{\text{d}\xi }=c\phi -\frac{a}{2}{\phi }^2-\frac{b}{3}{\phi }^3\end{array}$(3)

于是我们就把一个非线性偏微分高阶方程化简为一个二阶的常微分方程组。

接下来，令

$y=0,c\phi -\frac{a}{2}{\phi }^2-\frac{b}{3}{\phi }^3=0$

我们可以得到本系统的奇点。

2.2. 行波系统的首次积分

由微分方程组(3)，我们令

$\frac{\text{d}\phi }{\text{d}y}=\frac{y}{c\phi -\frac{a}{2}{\phi }^2-\frac{b}{3}{\phi }^3}$

解得

$y\left(\xi \right)=\pm \frac{\sqrt{6c{\phi }^2-2a{\phi }^3-b{\phi }^4+12C}}{\sqrt{6}}$

于是，我们得到积分函数

$H\left(\phi ,y\right)=6y^2+2a{\phi }^3+b{\phi }^4-6c{\phi }^2$

当H取值为固定的常数时，可以得到 $\phi$和y的关系，得到一条轨线。我们取经过某一点的H值，便得到了经过该点的轨线。取不同的值，我们得到许多条轨线，从而获得相图的形状。

3. 主要结果

3.1. 参数不同取值下的奇点情况

下面，我们需要研究方程：

$c\phi -\frac{a}{2}{\phi }^2-\frac{b}{3}{\phi }^3=0$

当参数c不为零时，由求根公式判别法可知，当 $b=-\left(3a^2\right)/16c$时，方程有一个根；当 $b<-\left(3a^2\right)/16c$且 $c>0$时，无实根；当 $b>-\left(3a^2\right)/16c$且 $c>0$时，有两个不同的根；当 $b>-\left(3a^2\right)/16c$且 $c<0$时，无实根；当 $b<-\left(3a^2\right)/16c$且 $c<0$时，有两个不同的根。

当 $c>0$时，系统奇点判断如下：

① $a=0,b=0$时，有 $\left(0,0\right)$一个鞍点；

② $a\ne 0,b=0$时，有一个鞍点 $\left(0,0\right)$，一个中心；

③ $a=0,b>0$时，有一个鞍点 $\left(0,0\right)$，两个中心；

④ $a\ne 0,b>0$时，有一个鞍点 $\left(0,0\right)$，两个中心；

⑤ $-\left(3a^2\right)/16c<b<0$时，有一个中心，两个鞍点，其中一个鞍点是 $\left(0,0\right)$；

⑥ $b<-\left(3a^2\right)/16c$时，只有一个鞍点 $\left(0,0\right)$；

⑦ $b=-\left(3a^2\right)/16c$时，有一个鞍点 $\left(0,0\right)$，和一个退化奇点。

$c>0$时的积分曲线图和数值模拟图见图1。

当 $c<0$时，奇点判断如下：

① $a=0,b=0$时，有 $\left(0,0\right)$一个中心；

② $a\ne 0,b=0$时，有一个中心 $\left(0,0\right)$，一个鞍点；

③ $a=0,b<0$时，有一个中心 $\left(0,0\right)$，两个鞍点；

④ $a\ne 0,b>0$时，有一个中心 $\left(0,0\right)$，两个鞍点；

⑤ $0<b<-\left(3a^2\right)/16c$时，有一个鞍点，两个中心，其中一个中心是 $\left(0,0\right)$；

⑥ $b>-\left(3a^2\right)/16c$时，只有一个中心 $\left(0,0\right)$；

⑦ $b=-\left(3a^2\right)/16c$时，有一个中心 $\left(0,0\right)$，和一个退化奇点。

$c<0$时的积分曲线图和数值模拟图见图2。

当 $c=0$时， $\left(0,0\right)$是一个退化奇点，奇点判断如下：

① $a=0,b=0$时，有 $\left(0,0\right)$一个退化奇点；

② $b>0$时，有 $\left(0,0\right)$一个退化奇点，还有一个中心；

③ $b<0$时，有 $\left(0,0\right)$一个退化奇点，还有一个鞍点。

$c=0$时的积分曲线图和数值模拟图见图3。

3.2. ab平面上的分支曲线和分支相图

接下来，为了得到分支曲线，我们令

$H\left(-\frac{3a-\sqrt{9a^2+48bc}}{4b},0\right)=H\left(0,0\right)=0$

解得结果为：

$b=-\frac{a^2}{6c}$

于是得到两条分支曲线 $b=-a^2/6c$和 $b=-3a^2/16c$。

下面，我们将ab平面用分支曲线分割成不同区域，见图4：

与之对应的分支相图分别见图5~7。

4. 总结

本文的工作难点在于，方程中带有未知数a、b，而且变量替换时，行波解中也带有未知数c，于是我们需要进行繁杂的分类讨论。本文把c分成大于零、小于零和等于零三种情况来讨论。然后再讨论ab的关系对于奇点坐标的影响，从而得到ab关系的分界线，也就是分支曲线。最后在c的不同取值下，在ab平面上作出分支相图。

通过观察发现，不管c的取值为何，每个ab平面上的分支相图都是关于b轴对称的，也就是说，a的取值只是改变了相图的左右对称的方向。而且每个相邻的区域之间的相图存在缓慢的过渡，在边界上达到定性的变化，然后逐渐展开。

参考文献

[1]	Dey, B. (1986) Domain Wall Solutions of KdV-Like Equations with Higher Order Nonlinearity. Journal of Physics A, 19, L9-L12. https://doi.org/10.1088/0305-4470/19/1/003
[2]	Dey, B. (1988) KdV Like Equations with Domain Wall Solutions and Their Hamiltonians. In: Solitons, Springer, Berlin, 188-194.
[3]	Zhang, J.F. (1998) New Solitary Wave Solution of the Combined KdV and mKdV Equation. International Journal of Theoretical Physics, 37, 1541-1546. https://doi.org/10.1023/A:1026615919186
[4]	Zhang, J.F., Wu, F.M. and Shi, J.Q. (2000) Simple Solition Solution Method for the Combined KdV and mKdV Equation. International Journal of Theoretical Physics, 39, 1697-1702. https://doi.org/10.1023/A:1003648715053
[5]	Yu, J. (2000) Exact Solitary Wave Solutions to a Combined KdV and mKdV Equation. Mathematical Methods in the Applied Sciences, 23, 1667-1670. https://doi.org/10.1002/1099-1476(200012)23:18<1667::AID-MMA180>3.0.CO;2-L
[6]	Grimshaw, R., Peli-novsky, D., Pelinovsky, E. and Slunyaev, A. (2002) Generation of Large-Amplitude Solitons in the Extended Korteweg-de Vries Equation. Chaos, 12, 1070-1076. https://doi.org/10.1063/1.1521391
[7]	Fan, E.G. (2002) Multiple Travelling Wave Solutions of Nonlinear Evolution Equations Using a Unified Algebraic Method. Journal of Physics A, 35, 6853-6872. https://doi.org/10.1088/0305-4470/35/32/306
[8]	Fan, E.G. (2003) Uniformly Con-structing a Series of Explicit Exact Solutions to Nonlinear Equations in Mathematical Physics. Chaos, Solitons and Fractals, 16, 819-839. https://doi.org/10.1016/S0960-0779(02)00472-1
[9]	Tang, M.Y., Wang, R.Q. and Jing, Z.J. (2002) Solitary Waves and Their Bifurcations of KdV Like Equation with Higher Order Nonlinearity. Science in China A, 45, 1255-1267.
[10]	Peng, Y.-Z. (2003) New Exact Solutions to the Combined KdV and mKdV Equation. International Journal of Theoretical Physics, 42, 863-868. https://doi.org/10.1023/A:1024475003055
[11]	刘正荣. 微分方程定性方法与数值模拟[M]. 广州: 华南理工大学, 2013.
[12]	宋明. 几类高次非线性波方程的行波解研究[D]: [博士学位论文]. 广州: 华南理工大学, 2014.