汪春江,舒 级,李 倩,王云肖,杨 袁
(四川师范大学 数学与软件科学学院,四川 成都 610066)
一类(3+1)维KdV方程的有理解及其怪波
汪春江,舒 级*,李 倩,王云肖,杨 袁
(四川师范大学 数学与软件科学学院,四川 成都 610066)
讨论一类经典的(3+1)维KdV方程,该方程在流体动力学、等离子物理、气体动力学等方面有广泛应用.通过一个简单的符号计算方法得到方程的有理解,并讨论了在某些条件下的怪波解.
KdV方程; 精确解; 符号计算方法; 有理解; 怪波
一直以来,非线性现象都是基础数学和应用数学研究的主题,非线性演化方程的精确解研究在数学物理上有着重大作用,常系数方程[1-3]、变系数方程[4-5]、随机方程[6-7]的精确解已经被广泛研究,并产生了很多关于精确解的方法,如逆散映射法[8-9]、Backlünd变换[10-11]、达布变换[12]、齐次平衡法[13]、(G′/G)-展开法[14]等.
近年来,怪波已成为国内外研究的焦点.从怪波解的形式上看,通常是有理分式.在海洋学和其他学科领域,科学家们都发现了怪波现象[15-17],例如,Bose-Einstein凝聚物怪波事件[18-19],在等离子体的空间和表面的异常波[20-21],尤其是在光学领域[22-23],当光脉冲在光子晶体纤维中传输高能量时,怪波就会存在.对于非线性Schrödinger方程[24-25]的怪波与驻波、多怪波和高阶怪波[26-28]、明暗怪波解[29-30]已广泛被讨论.怪波不仅出现在深水中,浅水中也发现了怪波[31].从直观上看,怪波具有超常的波高,因此大多数学者和研究人员只能从波高角度对它进行定义,即认为波高大于有效波高2倍(或2.2倍)的单波可以称为怪波.在浅水中,怪波的产生取决于调制不稳定性:当波高kH<1.363 m[32],非线性聚焦过程停止.以水为介质的波不同于一般的阈值,浅水波高为kH<1/3 m[33].
在本文中,考虑(3+1)维KdV方程
(1)
其中u是关于x、y、z、t的函数,x、y、z、t是独立变量.它是物理学家和数学家感兴趣的方程之一.KdV方程的复合解[35]、变系数KdV方程的精确解[36]、高阶KdV方程的精确解[37-39]、新的KdV-mKdV方程的孤波[40]、广义的Hirota Satsuma耦合方程和偶合的MKdV方程[41-42]的孤立波解已经被研究.
本文的目的是通过一个简单的符号计算方法[43],构建出(3+1)维KdV的怪波和有理解.
一个简单的符号计算方法对于求解非线性偏微分方程是有效的.下面给出求解(3+1)维偏微分方程的主要步骤.
步骤 1 利用截断展开法[44-46],作代换
(2)
其中,α是一个常数,j1,j2,j3,j4≥1(j1,j2,j3,j4∈0,1,2,…).通过上述变换将一个(3+1)维的偏微分方程
(3)
转化为一个双线性方程
(4)
其中F是f,ft,fx,fy,fz,fxx,fxy,fxz,…的多项式.
步骤 2 假设f是一个关于x、y、z、t的2N阶多项式,给定
(5)
其中,系数ai,j,k,l(0≤i,j,k,l≤2N)是常数,满足
(6)
如果j1≥1.
步骤 3 把(5)式代入(4)式,并令关于xiyjzktl(0≤i+j+k+l≤2N)的系数为0,得到一系列的多项式方程.
步骤 4 利用Maple软件,可以算出系数aijkl,i,j,k,l∈(0,1,2,…).
步骤 5 把满足条件(6)的系数aijkl(i,j,k,l∈{0,1,2,…})代入(2)式之后得到非线性偏微分方程(3)的解.
下面将应用上述方法来求解(3+1)维KdV方程的解.根据截断展开法,用变换
(7)
将(3+1)维KdV方程变为
令f成为上述方程的一个解
(8)
为了构建方程(1)的有理解,假定f是一个2阶关于x、y、z、t的多项式
(9)
其中系数aijkl(0≤i,j,k,l≤2)是常数,满足
(10)
为了简单起见,假定方程(9),a2000=a0200=a0020=a0000=1,a0002=2,代入到方程(8)中,并令关于xiyjzktl(0≤i+j+k+l≤2)的系数为0,得到15个多项式方程:
为了后面描述的方便,令a0011=b,a1100=c,a1001=d,a0110=e,a0101=f,a1010=g,a0001=h,a0010=j,a0100=k,a1000=l,解方程组,得到如下解.
情形 1
f=-ge,
情形 2
f=-ge,
情形 3
f=-ge,
情形 4
f=-ge,
情形 5
情形 6
情形 7
情形 8
情形 9
情形 10
情形 11
情形 12
情形 13
情形1~13代入(7)式得到(3+1)维KdV方程的解,情形1~4有2个自由变量a1100、a0010.在这些情形中可以得到实值型的怪波解.情形5~12有一个自由变量a0010.在这些情况中得到复值型的有理解.然而,在情形13中,有一个自由变量a0010.在这个情形中,可以得到实值型的有理解.
1) 从情形1中,挑选出系数:
可以得到(3+1)维KdV方程的实值型怪波解
其中
(12)
2) 从情形5中,挑选出系数:
可以得到(3+1)维KdV方程的复值型有理解
(13)
其中
(14)
3) 从情形13中,挑选出系数:
可以得到(3+1)维KdV方程的实值型有理解
其中
(16)
本文运用符号计算方法得到了(3+1)维KdV方程的有理解和怪波解,这些怪波和有理解是非奇异的.这些解对理解怪波的产生机制有一定帮助.下一步将研究如何用简单的符号计算方法构造非线性演化方程的高阶怪波解.
[1] JEFFREY A.Exact solutions tp the KdV-Burgers’ equation[J].Wave Motion,1991,14:369-375.
[2] KHALFALLAH M.New exact traveling wave solution of the (3+1) dimensional Kadomtsev-Petviashvili(KP) equation[J].Commun Nonlinear Sci Numer Simulat,2009,14:1169-1175.
[3] MALLORY K,VAN GORDER R A,VAJRAVELU K.Several classes of exact solutions to the (1+1) Born-Infeld equation[J].Commun Nonlinear Sci Numer Simulat,2014,19:1669-1674.
[4] ABDEL LATIF M S.Some exact solutions of KdV equation with variable coeffcients[J].Commun Nonlinear Sic Numer Simulat,2011,16:1783-1786.
[5] LI J B.Exact traveling wave solution and dynamical behavior for the (n+1)-dimensional multiple sine-Gordon equation[J].Science in China,2007,A50:153-164.
[6] KIM H S.Exact solutions of Wick-Type stochastic equations with variable coefficients[J].Reports on Mathematical Physics,2011,67:415-429.
[7] XIE Y.Exact solutions for stochastic KdV equations[J].Phys Lett,2003,A310:161-167.
[8] HEREMAN W,NUSEIR A.Symbolic methods to construct exact solutions of nonlinear partial differential equations[J].Math Comput Simul,1997,43:13-27.
[9] SACKS P,SHIN J.Computational methods for some inverse scattering problems[J].Appl Math Comput,2009,207:111-123.
[10] MA W X,ABDELJABBAR A.A bilinear Backlund transformation of a (3+1)-dimensional generalized KP equation[J].Applied Mathematics Letters,2012,25:1500-1504.
[11] QUINTINO A.Spectral deformation and Backlund transformation of constrained Willmore surfaces[J].Differential Geometry and its Applications,2011,29:S261-S270.
[12] BAGROV V G,SAMSONOV B F.Darboux transformation of the Schrödinger equation[J].Physics of Particles and Nuclei,1997,28:374-397.
[13] ESLAMI M,FATHIVAJARGAH B.Exact solutions of modified Zakharov-Kuznetsov equation by the homogeneous balance method[J].Ain Shams Engineering J,2014,5:221-225.
[14] WANG M,LI X,ZHANG J.The (G′/G)-expansion method and traveling wave solutions of nonlinear evolution equations in mathematical physics[J].Phys Lett,2008,A372:417-423.
[15] DYSTHE K,KROGSTAD H E.Oceanic rogue waves[J].Annu Rev Fluid Mech,2008,40:287-310.
[16] KIM H S.Exact solutions of Wick-type stochastic equation with variable coefficients[J].Reports on Mathematical Physics,2011,67:415-429.
[17] PELINOVSKY K E,SLUNYAEV A.Rogue Waves in the Ocean: Observation,Theories and Modeling[M].New York:Springer,2009.
[18] LIA J Z,LIU H Y,WANG Q.Enhanced multilevel linear sampling methods for inverse scattering problems[J].J Comput Phys,2014,257:554-571.
[19] SACKS P,SHIN J.Computatioal methods for some inverse scattering problems[J].Appl Math Comput,2009,207:111-123.
[20] RUDERMAN M S.Freak waves in laboratory and space plasmas[J].Eur Phys J:Special Topics,2010,185:57-66.
[21] MOSLEM W M,SHUKLA P K,ELIASSON B.Surface plasma rogue waves[J].Eur Phys Lett,2011,96:25002.
[22] SOLLI D R,ROPERS C,KOONATH P,et al.Optical rogue waves[J].Nature,2007,450:1054-1058.
[23] KIBLER B,FATOME J,FINOT C.The Peregrine soliton in nonlinear fibre optics[J].Nature Physics,2010,6:790-795.
[24] SHU J,ZHANG J.Sharp criterion of global existence for a class of nonlinear Schrödinger equation with critical exponent[J].Appl Math Comput,2006,182:1482-1487.
[25] HE J S,CHARALAMPADIS E G,KEVREKIDIS P G.Rogue waves in nonlinear Schrödinger models with variable coefficients: application to Bose-Einstein condensates[J].Phys Lett,2014,A378:557-583.
[26] YU F J.Multi-rogue waves for a higher-order nonlinear Schrodinger equation in optical fibers[J].Appl Math Comput,2013,220:176-184.
[27] DAI C Q,HUANG W H.Multi-rogue wave and multi-breather solutions in PT-symmetric coupled waveguides[J].Applied Mathematics Letters,2014,32:35-40.
[28] LI L J,WU Z W,WANG L H,et al.High-order rogue waves for the Hirota equatoin[J].Ann Phys,2013,343:198-211.
[29] YU F J,YAN Z Y.New rogue waves and dark-bright soliton solutions for a coupled nonlinear Schrödinger equation with variable coefficients[J].Appl Math Comput,2014,233:351-358.
[30] MA W X.Complexiton solutions to the Korteweg-de Vries equation[J].Phys Lett,2002,A301:35-44.
[31] SOOMERE T.Rogue waves in shallow water[J].Eur Phys J:Special Topics,2010,185:81-96.
[32] PETERSON P,SOOMERE T,ENGELBRECHT J.Soliton interaction as apossible model for extreme waves in shallow water[J].Nonlinear Processes in Geophysics,2003,10:503-510.
[33] RUDERMAN M S.Freak waves in laboratory and space plasmas[J].Eur Phys J:Special Topics,2010,185:57-66.
[34] BOUNTIS T C.Singularities and Dynamical Systems[M].New York:North-Holland,1985.
[35] MA H C,DENG A P,WANG Y.Exact solution of a KdV equation with variable coefficients[J].Appl Math Comput,2011,61:2278-2280.
[36] TRIKI H,WAZWAZ A-M.Traveling wave solutions for 5th-order KdV type equations with time-dependent coefficients[J].Commun Nonlinear Sic Numer Simulat,2014,19:404-408.
[37] OSBORNE A R,ONORATO M,SERIO M.The nonlinear dynamics of rogue waves and holes in deep-water gravity wave trains[J].Physics Letters,2000,A275:386-393.
[38] NOBUHIT M.Statistical properties of freak waves observed in the sea of Japan[C]//Proc International of Shore and Polar Engineering Conference,2000,3:109-115.
[39] DAVIDL K.Simulation of extreme waves in a background random sea[C]//Proc International of Shore and Polar Engineering Conference,2000,3:31-37.
[40] TAKASHI Y.Characteristics of giant freak waves observed in the sea of Japan[C]//Proc International 5-Symposium on Ocean Wave Measurement and Analysis,1997,1:316-328.
[41] HAVER S.A possible freak wave event measured at the draupner,jacket January 1 1995[C]//Actes de colloques-IFREMER,2004:1-8.
[42] OSBORNE A R,ONORATO M,SERIO M.The nonlinear dynamics of rogue waves and holes in deep water gravity wave trains[J].Physics Letters,2000,A275:386-393.
[43] OSBORNE A R.The random and deterministic dynamics of “rogue waves” in unidirectional deep-water wave trains[J].Marine Structures,2001,14:275-293.
[44] SLUNYAEV A,KHARIF C,PELINOVSKY E.Nonlinear wave focusing on water of finite depth[J].Physica D,2002,173:77-96.
[45] DYSTHE K,KROGSTAD H E,MULLER P.Oceanic rogue waves[J].Annu Rev Fluid Mech,2008,40:287-310.
[46] FAN E G.Solitons for a generalized Hirota-Satsuma coupled KdV equatioan and a coupled MKdV equation[J].Phys Lett,2001,A282:18-22.
2010 MSC:35Q55
(编辑 周 俊)
Rogue Waves and Rational Solutions for a Class of (3+1)-dimensional KdV Equation
WANG Chunjiang,SHU Ji,LI Qian,WANG Yunxiao,YANG Yuan
(CollegeofMathematicsandSoftwareScience,SichuanNormalUniversity,Chengdu610066,Sichuan)
This paper discusses a classical (3+1)-dimensional KdV equation,which has broad applications in hydrodynamics,plasma physics,gas dynamics.We obtain rational solutions of this equation by a simple symbolic computation approach.Under some conditions,we find that some of rational solutions are rogue waves.
KdV equation; exact solution; symbolic computation approach; rational solution; rogue wave
2016-03-30
四川省科技厅应用基础计划项目(2016JY0204)和四川省教育厅自然科学重点基金(14ZA0031)
O175.27
A
1001-8395(2017)02-0157-06
10.3969/j.issn.1001-8395.2017.02.003
*通信作者简介:舒 级(1976—),男,教授,主要从事随机动力系统和偏微分方程的研究,E-mail:shuji2008@hotmail.com