修正Benjamin⁃Bona⁃Mahony方程的精确行波解

杨春飞 ,  刘小华

内蒙古师范大学学报(自然科学版) ›› 2024, Vol. 53 ›› Issue (05) : 532 -541.

PDF (1629KB)
内蒙古师范大学学报(自然科学版) ›› 2024, Vol. 53 ›› Issue (05) : 532 -541. DOI: 10.3969/j.issn.1001-8735.2024.05.013

修正Benjamin⁃Bona⁃Mahony方程的精确行波解

作者信息 +

The Exact Traveling Wave Solution of the Modified Benjamin-Bona-Mahony Equation

Author information +
文章历史 +
PDF (1667K)

摘要

利用平面动力系统理论,研究修正Benjamin⁃Bona⁃Mahony (mBBM)方程行波解的存在性,得到不同参数条件下的相图和行波解的存在条件。运用广义Riccati辅助方程法给出mBBM方程的双曲解、三角解和有理函数解。并利用Maple软件,分析扭状孤波解、周期解的性状,数值上验证了精确行波解的表达式。

Abstract

The existence of the traveling wave solution of the modified Benjamin-Bona-Mahony (mBBM) equation is studied by using the theory of plane dynamical systems. The phase diagram and conditions for the existence of the traveling wave solution are obtained under different parameters. The generalized Riccati auxiliary equation method is used to derive double perversion, trigonometric, and rational function solutions for the mBBM equation. The properties of torsional solitary wave and periodic wave solutions are analyzed by using Maple software, and numerical verification confirms the expression of exact traveling wave solution.

Graphical abstract

关键词

mBBM方程 / 解的存在性 / 广义Riccati辅助方程法 / 行波解

Key words

mBBM equation / the existence of solutions / generalized Riccati auxiliary equation method / traveling wave solution

引用本文

引用格式 ▾
杨春飞,刘小华. 修正Benjamin⁃Bona⁃Mahony方程的精确行波解[J]. 内蒙古师范大学学报(自然科学版), 2024, 53(05): 532-541 DOI:10.3969/j.issn.1001-8735.2024.05.013

登录浏览全文

4963

注册一个新账户 忘记密码

构造非线性发展方程的精确行波解是孤波理论的研究课题之一,它为合理地解释相关非线性自然现象提供了理论依据,许多求解非线性发展方程的方法也日益建立和发展,比如:扩展的tanh函数法1、一阶积分法(FIM)2、广义G'/G展开法3、修正的辅助方程法4等。
本文主要探讨修正Benjamin⁃Bona⁃Mahony(mBBM)方程
ut+ux-αu2ux+uxxx=0
的精确行波解,文献[5]利用exp(-ϕ(ξ))展开法构造了方程(1)的精确行波解,包括双曲解、三角解和有理函数解;文献[6]分别利用tanh和sine⁃cosine方法构造了方程(1)的精确行波解;文献[7]使用截断Painlevé展开方法构造了方程(1)的孤立波解和椭圆函数周期波解;文献[8]利用双曲函数展开法构造了方程(1)的一类扭结-反扭结状的双扭结孤立波解。
广义Riccati辅助方程法实现了扩展tanh函数展开法和广义Riccati方程的结合,是研究非线性发展方程的先进数学工具和构造行波解的一种有效的方法,它可以用于其他非线性数学物理方程。例如:文献[9]利用广义Riccati辅助方程法构造了非线性传输线方程的精确行波解,包括双曲解、三角解和有理函数解,文献[10]利用广义Riccati辅助方程法构造了非线性Bogoyavlenskii方程的精确解族,包括孤子和其他解。通过应用该方法,可以以统一的方式为mBBM方程构造了丰富的行波解,丰富其解系。该方法简单、简洁,是一种常用于求解非线性发展方程精确解的有效方法。

1 方程(1)行波解的存在性

对方程(1)进行行波变换,设u(x,t)=u(ξ), ξ=x-ct,方程(1)可化作常微分方程

βu'-αu2u'+u=0

其中β=1-cc为波速,对方程(2)关于ξ积分一次,并令积分常数为零可得

βu-α3u3+u=0

x(ξ)=u(ξ),y(ξ)=u'(ξ),则可将(3)化为如下等价的Hamilton系统:

dxdξ=y ,dydξ=α3x3-βx 

系统(4)具有首次积分

H(x,y)=y22-α12x4+β2x2=h

其中h为哈密尔顿量。

接下来对系统(4)进行奇点分析,假设

f(x)=α3x3-βx=x(α3x2-β)=0

f(x)有三个零点,即x0=0,x1=3βα,x2=-3βα,故系统(4)有三个奇点,即

u0(0,0)u13βα,0,u2-3βα,0

由方程(6)可得f'(x)=αx2-β,故在奇点ui(i=0,1,2)处的Jacobi矩阵为

J(ui,0)=01αxi2-β0

其行列式值为J(ui,0)=DetJ(ui,0)=-αxi2+β,根据平面动力系统理论及方法1112可知,若J(ui,0)<0,则对应的奇点为鞍点;若J(ui,0)>0,则对应的奇点为中心点; 若J(ui,0)=0,则对应的奇点为尖点。因此对系统(4)进行奇点分析有如下结论成立。

(1)当βα>0时,系统(4)有三个奇点,即u0(0,0),u13βα,0,u2-3βα,0;当α>0,β>0时,u0为中心点,u1,2为鞍点(相图见图1(a));当α<0,β<0时,u0为鞍点,u1,2为中心点(相图见图1(b))。

(2)当βα<0时,系统(4)仅有一个奇点,即u0(0,0);当α<0,β>0时,u0为中心点(相图见图1(c));当α>0,β<0时,u0为鞍点(相图见图1(d))。

根据以上定性分析,利用数学软件Maple绘制系统(4)的相图如图1所示。由图1可知,系统(4)存在两条异宿轨,两条同宿轨,无数条闭轨和开轨。根据平面动力系统理论及方法1112,系统(4)的异宿轨、同宿轨、闭轨和开轨分别对应着方程(3)的扭状孤波解、钟状孤波解、周期解和无界行波解。并且根据行波变换,可得如下关于方程(1)行波解的性质。

性质1α>0,β>0时,方程(1)存在两个扭状孤波解和无数个周期解(见图1(a));当α<0,β<0时,方程(1)存在两个钟状孤波解和无数个周期解(见图1(b));当α<0,β>0时,方程(1)存在无数个周期解(见图1(c));当α>0,β<0时,方程(1)存在无界行波解(见图1(d))。

2 方程(1)的精确行波解

采用广义Riccati辅助方程法对方程(1)的行波解进行研究。首先给出广义Riccati辅助方程法的基本步骤。

2.1 广义Riccati辅助方程法

考虑如下(1+1)维的非线性发展方程

H(u,ut,ux,uxx,)=0

其中H为其变元u,ut,ux,uxx,的多项式。

步骤1 对方程(8)作如下的行波变换

u(x,t)=u(ξ),  ξ=x-ct

其中c为波速。将行波变换(9)代入方程(8),可得如下常微分方程

F(u,u',u,)=0

其中u'=dudξ,u=d2udξ2,

步骤2 假设方程(10)具有如下形式的解

u(ξ)=i=0maiϕi

其中系数ai(i=0,1,2,,m)为待定参数。ϕ=ϕ(ξ)满足如下广义Riccati方程

ϕ'(ξ)=qϕ2(ξ)+pϕ(ξ)+r

其中q,p,r是任意非零参数。

步骤3 利用齐次平衡原则,平衡方程(10)中的非线性项和最高阶导数项,由此可确定m的值。将式(11)以及式(12)代入方程(10),可得到一个关于ϕi(ξ)的多项式函数,合并ϕi(ξ)的同类项并令其系数等于零,可得到一个关于q,p,r,c,ai(i=0,1,2,,m)的代数方程组。

步骤4 利用数学软件Maple求解该代数方程组并结合方程(12)以下情况的解13(13)-(39),可得给出方程(10)的精确行波解,再根据行波变换(9),即可得出方程(8)的行波解的精确表达式。

根据文献[13],方程(12)具有以下几种情况的精确解。

(1)当Δ=p2-4qr>0pq0(或qr0)时,方程(12)有解

ϕ1(ξ)=-12qp+ΔtanhΔ2ξ,
ϕ2(ξ)=-12qp+ΔcothΔ2ξ,
ϕ3(ξ)=-12qp+ΔtanhΔξ±isechΔξ,
ϕ4(ξ)=-12qp+ΔcothΔξ±cschΔξ,
ϕ5(ξ)=-14q2p+ΔtanhΔ4ξ±cothΔ4ξ,
ϕ6(ξ)=12q-p+Δ(A2+B2)-AΔcoshΔξAsinhΔξ+B,
ϕ7(ξ)=12q-p-Δ(B2-A2)+AΔcoshΔξAsinhΔξ+B,

其中A,B是两个非零实数,且B2-A2>0

ϕ8(ξ)=2rcoshΔ2ξΔsinhΔ2ξ-pcoshΔ2ξ,
ϕ9(ξ)=-2rsinhΔ2ξpsinhΔ2ξ-ΔcoshΔ2ξ,
ϕ10(ξ)=2rcoshΔ2ξΔsinhΔξ-pcoshΔξ±iΔ,
ϕ11(ξ)=2rsinhΔ2ξ-psinhΔξ+ΔcoshΔξ±Δ,
ϕ12(ξ)=4rsinhΔ4ξcoshΔ4ξ-2psinhΔ4ξcoshΔ4ξ+2Δcosh2Δ4ξ-Δ

(2)当Δ=p2-4qr<0pq0(或qr0)时,方程(12)有解

ϕ13(ξ)=12q-p+-Δtan-Δ2ξ,
ϕ14(ξ)=-12qp+-Δcot-Δ2ξ,
ϕ15(ξ)=12q-p+-Δtan-Δξ±sec-Δξ,
ϕ16(ξ)=-12qp+-Δcot-Δξ±csc-Δξ,
ϕ17(ξ)=14q-2p+-Δtan-Δ4ξ-cot-Δ4ξ,
ϕ18(ξ)=12q-p+±-Δ(A2-B2)-A-Δcos-ΔξAsin-Δξ+B,
ϕ19(ξ)=12q-p-±-Δ(A2-B2)+A-Δcos-ΔξAsin-Δξ+B,

其中A,B是两个非零实数,且A2-B2>0

ϕ20(ξ)=2rcos-Δ2ξ-Δsin-Δ2ξ+pcos-Δ2ξ,
ϕ21(ξ)=2rsin-Δ2ξ-psin-Δ2ξ+-Δcos-Δ2ξ,
ϕ22(ξ)=-2rcos-Δ2ξ-Δsin-Δξ-pcos-Δξ±-Δ,
ϕ23(ξ)=2rsin-Δ2ξ-psin-Δξ+-Δcos-Δξ±-Δ,
ϕ24(ξ)=4rsin-Δ4ξcos-Δ4ξ-2psin-Δ4ξcos-Δ4ξ+2-Δcos2-Δ4ξ--Δ

(3)当r=0pq0时,方程(12)有解

ϕ25(ξ)=-pdqd+coshpξ-sinhpξ,
ϕ26(ξ)=pcoshpξ+sinhpξqd+coshpξ+sinhpξ

(4)当q0r=p=0时,方程(12)有解

ϕ27(ξ)=-1qξ+c1,

其中c1是任意常数。

2.2 方程(1)有界行波解的精确表达式

利用广义Riccati辅助方程法构造方程(1)的精确行波解,通过平衡方程(3)中的非线性项u3和最高阶导数项u,可得m=1。因此,假设方程(1)具有如下形式的解

u(ξ)=a0+a1ϕ(ξ)

其中a0,a1为待确定参数,ϕ(ξ)满足式(12)。由方程(40)和(12)可推得

u3(ξ)=a03+3a02a1ϕ(ξ)+3a0a12ϕ2(ξ)+a13ϕ3(ξ)
u(ξ)=pra1+(p2+2qr)a1ϕ(ξ)+3pqa1ϕ2(ξ)+2q2a1ϕ3(ξ)

将方程(40)、(41)和(42)代入方程(3),合并ϕi(ξ)(i=0,1,2,3)的同类项并令其系数等于零,可得到一个关于q,p,r,c,a0,a1的代数方程组

ϕ0(ξ):-αa033+pra1+(1-c)a0=0,ϕ1(ξ):-a1(αa02-p2-2qr+c-1)=0,ϕ2(ξ):-αa0a12+3pqa1=0,ϕ3(ξ):-13αa13+2q2a1=0

利用数学软件Maple求解该代数方程组,得到a0,a1,c的值为

a0=±p26α, a1=±q6α, c=-p22+2qr+1

其中α>0

由方程(40)、(44)以及方程(12)的解(13)-(39),可得方程(1)有界行波解的精确表达式。

(1)当Δ=p2-4qr>0pq0(或qr0)时,方程(1)有解

u1(x,t)=6αΔ2tanhΔ2ξ,
u2(x,t)=6αΔ2cothΔ2ξ,
u3(x,t)=6αΔ2tanhΔξ±isechΔξ,
u4(x,t)=6αΔ2cothΔξ±cschΔξ,
u5(x,t)=6αΔ4tanhΔ4ξ±cothΔ4ξ,
u6(x,t)=±6αΔ(A2+B2)-AΔcoshΔξ2AsinhΔξ+B,
u7(x,t)=6αΔ(B2-A2)+AΔcoshΔξ2AsinhΔξ+B,

其中A,B是两个非零实数,且B2-A2>0

u8(x,t)=±6αp2+2qrcoshΔ2ξΔsinhΔ2ξ-pcoshΔ2ξ,
u9(x,t)=±6αp2-2qrsinhΔ2ξpsinhΔ2ξ-ΔcoshΔ2ξ,
u10(x,t)=±6αp2+2qrcoshΔ2ξΔsinhΔξ-pcoshΔξ±iΔ,
u11(x,t)=±6αp2+2qrsinhΔ2ξ-psinhΔξ+ΔcoshΔξ±Δ,
u12(x,t)=±6αp2+4qrsinhΔ4ξcoshΔ4ξ-2psinhΔ4ξcoshΔ4ξ+2Δcosh2Δ4ξ-Δ

(2)当Δ=p2-4qr<0pq0(或qr0)时,方程(1)有解

u13(x,t)=±6α-Δ2tan-Δ2ξ,
u14(x,t)=6α-Δ2cot-Δ2ξ,
u15(x,t)=±6α-Δ2tan-Δξ±sec-Δξ,
u16(x,t)=6α-Δ2cot-Δξ±csc-Δξ,
u17(x,t)=±6α-Δ4tan-Δ4ξ-cot-Δ4ξ,
u18(x,t)=±6α±-Δ(A2-B2)-A-Δcos-Δξ2Asin-Δξ+B,
u19(x,t)=±6α±-Δ(A2-B2)+A-Δcos-Δξ2Asin-Δξ+B,

其中A,B是两个非零实数,且A2-B2>0

u20(x,t)=±6αp2-2qrcos-Δ2ξ-Δsin-Δ2ξ+pcos-Δ2ξ,
u21(x,t)=±6αp2+2qrsin-Δ2ξ-psin-Δ2ξ+-Δcos-Δ2ξ,
u22(x,t)=±6αp2-2qrcos-Δ2ξ-Δsin-Δξ-pcos-Δξ±-Δ,
u23(x,t)=±6αp2+2qrsin-Δ2ξ-psin-Δξ+-Δcos-Δξ±-Δ,
u24(x,t)=±6αp2+4qrsin-Δ4ξcos-Δ4ξ-2psin-Δ4ξcos-Δ4ξ+2-Δcos2-Δ4ξ--Δ

(3)当r=0pq0时,方程(1)有解

u25(x,t)=±6αp2-pdd+coshpξ-sinhpξ,
u26(x,t)=±6αp2-pcoshpξ+sinhpξd+coshpξ+sinhpξ

(4)当q0r=p=0时,方程(1)有解

u27(x,t)=±6αp2-1qξ+c1,

其中ξ=x--p22+2qr+1tc1是任意常数。

3 解的性态图

方程(1)的部分解在特定参数值下的性态3D图和2D图,如图2-4所示。图2是方程(1)的解u1在区间-10<x,t<10内的3D图(参数为α=6,p=q=2,r=0.25)和2D图(参数为t=1)。由图2可以看出,u1具有单调递减性,为扭状孤立波解。图3是方程(1)的解u2在区间-10<x,t<10内的3D图(参数为α=6,p=q=2,r=0.25)和2D图(参数为t=1)。由图3可以看出,此解为奇异扭状孤立波解。图4是方程(1)的解u13在区间-10<x,t<10内的3D图(参数为α=6,p=q=2,r=1)和2D图(参数为t=1)。由图4可以看出,此解为周期解。

4 结论

本文主要讨论了修正Benjamin⁃Bona⁃Mahony (mBBM)方程行波解的存在性及精确表达式,并且数值上验证了精确解的性态,这些解的构造不仅丰富了mBBM方程的解系,更为深入理解mBBM方程解的性质奠定了基础,对更全面地理解由mBBM方程所刻画的物理现象提供参考。

参考文献

[1]

ALI K KMEHANNA M SWAZWAZ A M. Analytical and numerical treatment to the (2+1)-dimensional Date-Jimbo-Kashiwara-Miwa equation[J]. Nonlinear Engineering202110(1): 187-200.

[2]

JAVEED SSAIF SWAHEED Aet al. Exact solutions of fractional mBBM equation and coupled system of fractional Boussinesq-Burgers[J]. Results in Physics20189: 1275-1281.

[3]

MOHANTY S KKRAVCHENKO O VDEKA M Ket al. The exact solutions of the (2+1)-dimensional Kadomtsev-Petviashvili equation with variable coefficients by extended generalized G′/G-expansion method[J]. Journal of King Saud University-Science202335(1): 102358.

[4]

ALTHOBAITI A. Travelling waves solutions of the KP equation in weakly dispersive media[J]. Open Physics202220(1): 715-723.

[5]

KHAN KAKBAR M ALL. Application of exp-expansion method to find the exact solutions of modified Benjamin-Bona-Mahony equation[J]. World Applied Sciences Journal201324(10): 1373-1377.

[6]

YUSUFOĞLU EBEKIR A. The tanh and the sine-cosine methods for exact solutions of the MBBM and the Vakhnenko equations[J]. Chaos, Solitons & Fractals, 200838(4): 1126-1133.

[7]

YAN X WTIAN S FDONG M Jet al. Nonlocal symmetries, conservation laws and interaction solutions of the generalised dispersive modified Benjamin-Bona-Mahony equation[J]. Zeitschrift Für Naturforschung A201873(5): 399-405.

[8]

王林雪,宗谨,王雪玲,.mBBM方程的双扭结孤立波及其动力学稳定性[J].计算物理201633(2):212-220.

[9]

MALWE B HBETCHEWE GDOKA S Yet al. Travelling wave solutions and soliton solutions for the nonlinear transmission line using the generalized Riccati equation mapping method[J]. Nonlinear Dynamics201684(1): 171-177.

[10]

ZAYED E M EAL-NOWEHY A G. Solitons and other solutions to the nonlinear Bogoyavlenskii equations using the generalized Riccati equation mapping method[J]. Optical and Quantum Electronics201749(11): 359.

[11]

马知恩,周义仓.常微分方程定性与稳定性方法[M].北京:科学出版社, 2001.

[12]

张芷芬.微分方程定性理论[M].北京:科学出版社,1985.

[13]

AHMED NBABER M ZIQBAL M Set al. Analytical study of reaction diffusion Lengyel-Epstein system by generalized Riccati equation mapping method[J]. Scientific Reports202313(1): 20033.

基金资助

贵州省教育厅自然科学研究资助项目“贵州省高等学校光通讯系统中孤子的数学理论与计算协作创新团队”(黔教技〔2023〕062号)

AI Summary AI Mindmap
PDF (1629KB)

87

访问

0

被引

详细

导航
相关文章

AI思维导图

/