Degasperis-Procesi方程是偏微分方程中非常重要的一个方程。本文主要介绍扩展的Degasperis-Procesi方程,利用行波变换将偏微分方程转换成常微分方程,再对方程中的参数进行分析得到相图分支,根据相图轨道构建出方程精确行波解,求出其参数表达式。同时证明了在某些参数条件下,扩展的Degasperis-Procesi方程具有光滑的孤立波解和周期波解。 Degasperis-Procesi equation is one of the most important partial differential equations. The presented paper mainly introduces the extended Degasperis-Procesi equation. The partial differential equation is transformed into ordinary differential equation by traveling wave transformation, and then the parameters in the formula are analyzed to obtain the phase diagram branch. According to the phase diagram orbit, the exact traveling wave solutions for the equation are constructed. It is proved that the extended Degasperis-Procesi equation has smooth solitary wave solution and periodic wave solution under some certain parametric conditions.
Degasperis-Procesi方程是偏微分方程中非常重要的一个方程。本文主要介绍扩展的Degasperis-Procesi方程,利用行波变换将偏微分方程转换成常微分方程,再对方程中的参数进行分析得到相图分支,根据相图轨道构建出方程精确行波解,求出其参数表达式。同时证明了在某些参数条件下,扩展的Degasperis-Procesi方程具有光滑的孤立波解和周期波解。
扩展的Degasperis-Procesi方程,动力系统分支理论,同宿轨,周期轨,行波解
Minzhi Wei1, Fengling Yang2*
1Guangxi University of Finance and Economics, Nanning Guangxi
2Bank of Communications Co., Ltd. Liuzhou Branch, Liuzhou Guangxi
Received: Jan. 23rd, 2021; accepted: Feb. 17th, 2021; published: Feb. 25th, 2021
Degasperis-Procesi equation is one of the most important partial differential equations. The presented paper mainly introduces the extended Degasperis-Procesi equation. The partial differential equation is transformed into ordinary differential equation by traveling wave transformation, and then the parameters in the formula are analyzed to obtain the phase diagram branch. According to the phase diagram orbit, the exact traveling wave solutions for the equation are constructed. It is proved that the extended Degasperis-Procesi equation has smooth solitary wave solution and periodic wave solution under some certain parametric conditions.
Keywords:Extended Degasperis-Procesi Equation, Bifurcation Theory of Dynamical Systems, Homoclinic Orbit, Periodic Orbit, Traveling Wave Solution
Copyright © 2021 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
孤立波最初由英国科学家罗素(Russell)于1842年发现并于1844年提出 [
1965年Zabusky和Kruskal发现,假定FPU问题的震动是连续的,就可以用KdV方程近似描述。对于KdV方程中波速不同的两个孤波,相互碰撞后,只发生位移,而波形波速不改变,类似于粒子的碰撞,所以他们将这样的孤波命名为孤立子。Cammassa和Holm于 [
与孤立波理论研究相关的方程有很多,其中就包括Cammassa-Holm (CH)方程和Degasperis-Procesi (DP)方程。Degasperis和Procesi [
u t − u t x x + ( a + 1 ) u u x = a u x u x x + u u x x x , (1)
当且仅当 a = 2 或 a = 3 时,方程(1)完全可积。当 a = 2 时,方程(1)是Cammassa-Holm (CH)方程,可积且具有多孤子尖波峰解;当 a = 3 时,方程(1)是无色散形式的Degasperis-Procesi (DP)方程。
最近Helal [
u t − u t x x + ( a + 1 ) u 2 u x = ( b + 1 ) u x u x x + u u x x x . (2)
情况1: a = 2 , b = 1 ,方程(2)是扩展的Camassa-Holm (mCH)浅水波方程(见 [
情况2: a = 3 , b = 2 ,方程(2)是扩展的Degasperis-Procesi (mDP)方程(见 [
本文研究 b = 2 , a ≠ 3 的情况,即以下形式的扩展Degasperis-Procesi方程:
u t − u t x x + ( a + 1 ) u 2 u x = 3 u x u x x + u u x x x , (3)
利用动力系统分支理论 [
对方程(3)作行波变换 u ( x , t ) = u ( ξ ) , ξ = x − c t ,其中c是波速。代入式(3),积分忽略积分常数,方程(3)转化为如下的常微分方程:
− c u + c u ″ + a + 1 3 u 3 = ( u ′ ) 2 + u u ″ , (4)
将其转化为以下的二维系统:
d u d ξ = y ; d y d ξ = c u − a + 1 3 u 3 + y 2 c − u , (5)
并得到它的首次积分为:
H ( u , y ) = − c 2 u 2 + 2 3 c u 3 + 1 6 c ( a + 1 ) u 4 − 2 15 ( a + 1 ) u 5 + ( c − u ) 2 y 2 = h , (6)
通过变形可以得到:
y 2 = h + c 2 u 2 − 2 3 c u 3 − 1 6 c ( a + 1 ) u 4 + 2 15 ( a + 1 ) u 5 ( c − u ) 2 , (7)
当 c − u = 0 时,记 h = h s 。
可以看出,方程(4)是一个平面动力系统,定义在参数空间 ( c , a , n ) 上。由于在实际模型中,有界解才具有物理意义,同时(4)式的相图轨道决定了的所有行波解,所以我们将研究(4)式在相平面 ( u , y ) 中的相图分支,对于给定的物理模型,仅求出方程(4)的有界解。通过观察在系统参数发生变化时,相图的分支情况,得出方程的行波解。
设 u ( ξ ) 是(5)在 ξ ∈ ( − ∞ , + ∞ ) 上的一个连续解,且 lim ξ → − ∞ u ( ξ ) = a , lim ξ → + ∞ u ( ξ ) = b 。①若 a = b ,称 u ( x , t ) 为一个孤立波解;②若 a ≠ b ,称 u ( x , t ) 为扭结或反扭结波解。式(1)的一个孤立波解对应式(5)的一个同宿轨道,式(1)的一个扭结或反扭结波解对应式(5)的一个异宿轨道。因此想要分析所有式(1)孤立波的相图分支,就得先找到式(5)基于系统参数的同宿轨道。
考虑式(5),首先令 d ξ = ( c − u ) d τ ,式(5)变为
d u d τ = ( c − u ) y ; d y d τ = c u − a + 1 3 u 3 + y 2 , (8)
显然,除直线 c − u = 0 外,式(8)与式(5)具有相同的拓扑相图。且式(7)与式(5)可积,首次积分与式(6)相同。当h不变时,式(7)确定式(8)的一组不变曲线,包含不同的曲线分支;当h变化时,式(7)定义了式(8)不同曲线的动力学行为。
令 F ( u ) = c u − a + 1 3 u 3 ,为了研究系统的临界点,我们需要找到方程 F ( u ) = 0 的所有零点。在 ( u , y ) 相平面中,方程(5)在u轴上的横坐标是 F ( u ) 的零点。
可以发现,若 c ( a + 1 ) > 0 ,则 F ( u ) 有三个平衡点,分别是 E 0 ( 0 , 0 ) 、 G 1 , 2 ( 3 c a + 1 , 0 ) 。若 c ( a + 1 ) < 0 ,则 F ( u ) 有一个平衡点 E 0 ( 0 , 0 ) 。
设 E ( u i , 0 ) 是式(7)的平衡点, M ( u i , 0 ) 是式(7)线性化系统在 E ( u i , 0 ) 处的系数矩阵, J ( u i , 0 ) 是它的雅可比行列式,则
J ( u i , 0 ) = det M ( u i , 0 ) = − ( c − u i ) [ c − ( a + 1 ) u i 2 ]
由平面动力系统分支定理可知,当 J ( u i , 0 ) > 0 时,哈密顿系统平衡点 E ( u i , 0 ) 是中心;当 J ( u i , 0 ) < 0 时, E ( u i , 0 ) 是鞍点;当 J ( u i , 0 ) < 0 ,且Poincaré指数为0时, E ( u i , 0 ) 是尖点。
由(6)定义的函数,有
h ± = H ( u ± , 0 ) = u 2 ( − c 2 + 2 3 c u + 1 6 c ( a + 1 ) u 2 − 2 15 ( a + 1 ) u 3 ) ; h o = H ( ± y , 0 ) = 0 , (9)
由以上分析得到四种不同的参数条件下的相图分支,见图1。
图1. 系统(5)在不同参数条件下的相图
可以发现,(1-1) c > 0 , a + 1 < 0 、(1-2) c < 0 , a + 1 > 0 两种情况下,方程不存在行波解,(1-3)与(1-4)对称。
本章根据第二章所求得的结论,对动力系统进行相图分析,在 取不同的值时,得到不同的轨线分支,通过(7)式的变形,借助文献 [
由图1(1-3)可知,当 h = 0 时,出现同宿轨,则(7)可变形为
y 2 = u 2 [ c 2 − 2 c 3 u − c ( a + 1 ) 6 u 2 + 2 15 ( a + 1 ) u 3 ] ( c − u ) 2 = u 2 ( u − t 1 ) ( u − z 1 ) ( u − z 2 ) ( c − u ) 2 , (10)
其中 t 1 是 c 2 − 2 c 3 u − c ( a + 1 ) 6 u 2 + 2 15 ( a + 1 ) u 3 的实根, z 1 、 z 2 是它的复数根。结合(10)和(5)的第一个式子得
± d ξ = c − u u 1 ( u − t 1 ) ( u − z 1 ) ( u − z 2 ) d u , (11)
由 [
± ξ = g 1 F ( φ 1 , k 1 ) − c g 1 A − t 1 [ − u + A + t 1 2 t 1 ∏ ( φ 1 , α 2 α 2 − 1 , k 1 ) − α f 1 ] , (12)
其中,
A 2 = ( z 1 + z 2 2 − t 1 ) 2 − ( z 1 − z 2 ) 2 4 , g 1 = [ ( z 1 + z 2 2 − t 1 ) 2 − ( z 1 − z 2 ) 2 4 ] − 1 4 ,
k 1 2 = A + z 1 + z 2 2 − t 1 2 A , k ′ 1 = 1 − k 1 2 , φ 1 = cos − 1 ( u − t 1 − A u − t 1 + A ) , α = A − t 1 A + t 1 ,
f 1 = 1 − α 2 k 1 2 + α 2 k ′ 1 2 tan − 1 ( k 1 2 + α 2 k ′ 1 2 1 − α 2 s d u ) .
(12)是方程(3)在 c > 0 , a + 1 > 0 条件下孤立波解的参数表达式。同理,也可求出方程在 c < 0 , a + 1 < 0 条件下孤立波解的参数表达式。
由图1(1-3)可知,当 h = h s 时,出现周期轨,则(7)式可变形为
y 2 = ( c − u ) 2 ( u − t 1 ) ( u − t 2 ) ( u − t 3 ) ( c − u ) 2 = ( u − t 2 ) ( u − t 3 ) ( u − t 4 ) , (13)
其中 t 2 > t 3 > t 4 。结合(13)和(5)的第一个式子得
± d ξ = d u ( u − t 2 ) ( u − t 3 ) ( u − t 4 ) , (14)
由 [
± ξ = g 2 s n − 1 ( sin φ 2 , k 2 ) , (15)
由(14)式计算(15)式,得
u = t 2 − t 3 s n 2 ( ξ g 2 , k 2 ) 1 − s n 2 ( ξ g 2 , k 2 ) , (16)
其中, g 2 = 2 t 2 − t 3 , φ 2 = sin − 1 u − t 2 u − t 3 , k 2 2 = t 3 − t 4 t 2 − t 4 , s n u = sin φ 2 。
本文利用动力系统分支理论研究扩展的Degasperis-Procesi方程,在不同的参数条件下获得四个相图,证明了孤立波解、周期波解的存在,并获得对应精确解的参数表达式。
本文受广西高校中青年教师科研基础能力提升项目(No. 2020KY16019)和广西财经学院青年教师科研发展基金项目(No. 2019QNA03)资助。
韦敏志,杨凤玲. 扩展的Degasperis-Procesi方程的行波解分支Traveling Wave Solutions and Bifurcations for Extended Degasperis-Procesi Equation[J]. 应用数学进展, 2021, 10(02): 531-536. https://doi.org/10.12677/AAM.2021.102057