The Exp-function method is generalized to construct N-soliton solutions of a new generalization of the associated Camassa-Holm equation. As a result, one-soliton, two-soliton, and three-soliton solutions are obtained, from which the uniform formulae of N-soliton solutions are derived. It is shown that the Exp-function method may provide us with a straightforward, effective, and alternative mathematical tool for generating N-soliton solutions of nonlinear evolution equations in mathematical physics.
1. Introduction
The investigation of the traveling wave solutions to nonlinear evolution equations (NLEEs) plays an important role in mathematical physics. A lot of physical models have supported a wide variety of solitary wave solutions. In the recent years, much efforts have been spent on this task and many significant methods have been established such as inverse scattering transform [1], Bäcklund and Darboux transform [2], Hirota bilinear method [3], homogeneous balance method [4], Jacobi elliptic function method [5], tanh-function method [6], Exp-function method [7], simple equation method [8], F-expansion method [9, 10], improved F-expansion method [11], and extended F-expansion method [12]. Here, we study a new generalization of the associated Camassa-Holm equation.
The Camassa-Holm (CH) equation
(1)UT+2k2UX-UXXT+3UUX=2UXUXX+UUXXX,
where k is a nonzero real constant, was derived as a model for shallow water waves by Camassa and Holm in 1993 [13]. This equation is integrable with the following Lax pair:
(2)ψxx=λ(u-uxx+k2)ψ+14ψ,ψt=(12λ-u)ψx+12uxψ.
Considerable interest was paid to the CH equation in recent decades about its integrability and exact solutions [14–20]. Schiff and Fisher showed that the Camassa-Holm equation possessed the Bäcklund transformations and an infinite number of local conserved quantities by using the loop group approach [15, 16]. Parker gave explicit multisoliton solutions for the CH equation by taking the Hirota bilinear method and a coordinate transformation [18]. Its structure and dynamics were investigated in the different parameter regimes. According to [18], there is a reciprocal transformation, (T,X)→(t,x), such that
(3)dx=RdX-URdT,dt=dT,R=U-UXX+k2.
Let us apply the reciprocal transformation to the Lax pair (2) and define the following potential function u(x,t):
(4)u=12RxxR-14(RxR)2+14R2-14k2;
then (1) is transformed into the following associated Camassa-Holm (ACH) equation:
(5)ut+2k3ux+4k2uut+2k2ux∂x-1ut-k2uxxt=0.
Hone showed in [21] how the ACH equation (5) is related to Schrodinger operators and the KdV equation and described how to construct solutions of the ACH equation from tau-functions of the KdV hierarchy, including rational, N-soliton, and elliptic solutions.
Recently, integrable negative order flows, mixed equations, and the relationship of different hierarchies attracted much attention, including continuous and discrete cases, such as the negative KdV, mixed KdV, and Volterra lattice equations. In [22], Luo et al. introduced a new generalization of the associated Camassa-Holm equation; namely,
(6)ut+α(uxxx-6uux)+β(2k3ux+4k2uut+2k2ux∂x-1ut-k2uxxt)=0,
where α and β are two arbitrary constants.
Apparently, (6) is reduced to the ACH equation (5) when we take α=0, β=1. For α=1, β=0, (6) gives the KdV equation. So, (6) may be called the ACH-KdV equation. In [22], Luo et al. show that (6) is integrable in the sense of Lax pair and construct some exact solutions of (6) by Darboux transformation through Lax pair.
The Exp-function method [7] proposed by He and Wu in 2006 provides us with a straightforward and effective method for obtaining generalized solitary wave solutions and periodic solutions of NLEEs. The method was used by many researchers to study various NLEEs. More recently, Marinakis [23] did very interesting work to generalize the Exp-function method for constructing N-soliton solutions of NLEEs. Marinakis chose the famous Korteweg-de Vries (KdV) equation to illustrate the generalized work and successfully obtained the known 2-soliton and 3-soliton solutions in a simple and straightforward way.
In the present paper, we would like to generalize the Exp-function method for constructing N-soliton solutions of ACH-KdV equation (6).
The rest of this paper is organized as follows. In Section 2, we give the description of the Exp-function method for constructing N-soliton solutions of NLEEs. In Section 3, we apply the method to (6). In Section 4, some conclusions and discussions are given.
2. Basic Idea of the Exp-Function Method for N-Soliton Solutions of NLEEs
In this section, we recall the Exp-function method [23] for N-soliton solutions of NLEEs. For a given NLEE, say, in two variables x and t,
(7)F(u,ux,ut,uxx,uxt,utt,…)=0,
the Exp-function method for one-soliton solution is based on the assumption
(8)u(x,t)=∑i1=0p1ai1ei1ξ1∑j1=0q1aj1ej1ξ1,ξ1=w1x-c1t+λ1,
where ai1, bj1, k1, c1, λ1, and w1 are unknown constants and the values of p1 and q1 can be determined by balancing the linear term of highest order in (7) with the highest order nonlinear term.
In order to seek N-soliton solutions for integer N>1, we generalize (8) to the following form:
(9)u(x,t)=∑i1=0p1∑i2=0p2⋯∑iN=0pNai1i2⋯iNe∑g=1Nigξg∑j1=0q1∑j2=0q2⋯∑jN=0qNaj1j2⋯jNe∑g=1Njgξg,hhhhhhhhhhhhhhhhξg=wgx-cgt+λg;
given the value of N=2, it becomes
(10)u(x,t)=∑i1=0p1∑i2=0p2ai1i2e∑g=12igξg∑j1=0q1∑j2=0q2aj1j2e∑g=12jgξg,
which can be used to construct two-soliton solution.
When N=3, (9) changes into
(11)u(x,t)=∑i1=0p1∑i2=0p2∑i3=0p3ai1i2i3e∑g=13igξg∑j1=0q1∑j2=0q2∑j3=0q3aj1j2j3e∑g=13jgξg,
which can be used to obtain three-soliton solution.
Substituting (10) into (7) and using Mathematica, then equating to zero each coefficient of the same order power of the exponential functions yields a set of equations. Solving the set of equations, we can determine the 2-soliton solution and the following 3-soliton solution by means of (11), provided they exist. If possible, we may conclude with the uniform formula of N-soliton solutions for any N⩾1.
3. Multisoliton Solutions of the ACH-KdV Equation
In this section we apply Exp-function method to the ACH-KdV equation (6). We first remove the integral term in (6) by introducing the potential
(12)u(x,t)=qx(x,t);
then substituting (12) into (6) yields
(13)qxt+α(qxxxx-6qxqxx)+β(2k3qxx+4k2qxqxt+2k2qxxqt-k2qxxxt)=0.
Suppose that (13) admits the one-soliton solution of the form
(14)q(x,t)=a1eξ11+b1eξ1,
where ξ1=w1x-c1t+λ1, w1, c1 are undetermined constants, and λ1 is an arbitrary constant. Obviously, (14) is included in the same form as (7). Substituting (14) into (13) and then equating to zero each coefficient of the same order power of eiξ1 (i=1,2,3,4) yield a set of equations for a1, b1, w1, and c1 as follows:
(15)βk2a1c1w13+2βk3a1w12+αa1w14-a1c1w1=0,-11βk2a1b1c1w13+2βk3a1b1w12-6βk2a12c1w12-11αa1b1w14-6αa12w13-a1b1c1w1=0,11βk2a1b12c1w13-2βk3a1b12w12+6βk2a12b1c1w12+11αa1b12w14+6αa12b1w13+a1b12c1w1=0,-βk2a1b13c1w13-2βk3a1b13w12-αa1b13w14+a1b13c1w1=0.
Solving these equations by Maple, one has
(16)a1=-2b1w1,b1=b1,w1=w1,c1=w1(2βk3+αw12)1-βk2w12.
Substituting (16) into (14), we have
(17)q(x,t)=-2b1w1eξ11+b1eξ1.
Using potential (12), we can get the one-soliton solutions of (6) as follows:
(18)u(x,t)=qx(x,t)=-2b1w12eξ1(1+b1eξ1)2,
where ξ1=w1x-(w1(2βk3+αw12)/(1-βk2w12))t+λ1 and w1, λ1, and b1 are arbitrary constants. The one-soliton solution (18) is shown in Figure 1.
Figures of solution (18) and with α=1, β=-2, k=1, w1=1, b1=1, λ1=0. (a) Spatial plots in the intervals x∈[-5,5] and t∈[-5,5]; (b) plan plots with t=0, in the interval x∈[-10,10].
Next, we suppose that (13) has the 2-soliton solution in the form
(19)q(x,t)=a1eξ1+a2eξ2+a3eξ1+ξ21+b1eξ1+b2eξ2+b3eξ1+ξ2,
where ξ1=w1x-c1t+λ1; ξ2=w2x-c2t+λ2; a1, a2, a3, b1, b2, b3, w1, and w2 are constants to be determined; and λ1 and λ2 are arbitrary constants. Obviously, (19) has the same form as (10). Substituting (19) into (13) and using manipulations similar to those illustrated above, we obtain
(20)a1=-2b1w1,a2=-2b2w2,a3=-2b1b2(w1+w2)A12,b3=b1b2A12,c1=w1(2βk2+αw12)1-βk2w12,c2=w2(2βk2+αw22)1-βk2w22,
where
(21)A12=(w1-w2)2(w1+w2)2.
Substituting (20) into (19), we have
(22)q(x,t)=-2b1w1eξ1+b2w2eξ2+b1b2(w1+w2)A12eξ1+ξ21+b1eξ1+b2eξ2+b1b2A12eξ1+ξ2.
Using potential (12), we can get the two-soliton solutions of (6) as follows:
(23)u(x,t)=-2((1+b1eξ1+b2eξ2+b1b2A12eξ1+ξ2)-2(b12b2w22A12e2ξ1+ξ2+b1b22w12A12eξ1+2ξ2hhhhihh+2b1b2(w1-w2)2eξ1+ξ2+b1w12eξ1+b2w22eξ2)hhhhh×(1+b1eξ1+b2eξ2+b1b2A12eξ1+ξ2)-2),
where ξ1=w1x-(w1(2βk3+αw12)/(1-βk2w12))t+λ1; ξ2=w2x-(w2(2βk3+αw22)/(1-βk2w22))t+λ2; and w1, w2, λ1, λ2, b1, and b2 are arbitrary constants. The two-soliton solution (18) is shown in Figure 2.
Figures of solution (23) and with α=1, β=-1, w1=1, w2=1.1, b1=1, b2=1, k=1, λ1=0, λ2=0. (a) Spatial plots in the intervals x∈[-10,15] and t∈[-10,10]; (b)–(d) plan plots and with t=-10,0,10, x∈[-5,20], x∈[-15,15], x∈[-15,15], respectively.
From Figure 2 we can find the wave chase phenomenon. The speed of the first wave w1=1; the speed of the second wave w2=1.1. As time increases, the second wave exceeds the first one.
In what follows, we now suppose that the three-soliton solution of (13) can be expressed as follows:
(24)q(x,t)=(a1eξ1+a2eξ2+a3eξ3+a4eξ1+ξ2hhh+a5eξ1+ξ3+a6eξ2+ξ3+a7eξ1+ξ2+ξ3)×(1+b1eξ1+b2eξ2+b3eξ3+b4eξ1+ξ2hhhih+b5eξ1+ξ3+b6eξ2+ξ3+b7eξ1+ξ2+ξ3)-1,
where, ξ1=w1x-c1t+λ1; ξ2=w2x-c2t+λ2; ξ3=w3x-c3t+λ3; ai,bi (i=1,2,…,7), w1, w2, and w3 are constants to be determined; and λ2,λ2, and λ3 are arbitrary constants. Obviously, (24) has the same form as (11). Substituting (24) into (13) and using manipulations similar to those illustrated above, we obtain
(25)a1=-2w1b1,a2=-2w2b2,a3=-2w3b3,a4=-2b1b2(w1+w2)A12,a5=-2b1b3(w1+w3)A13,a6=-2b2b3(w2+w3)A23,a7=-2b1b2b3(w1+w2+w3)A12A13A23,b4=b1b2A12,b5=b1b3A13,b6=b2b3A23,b7=b1b2b3A12A13A23,c1=w1(2βk2+αw12)1-βk2w12,c2=w2(2βk2+αw22)1-βk2w22,c3=w3(2βk2+αw32)1-βk2w32,
where
(26)Aij=(wi-wj)2(wi+wj)2,1⩽i<j⩽3.
Substituting (25) into (24) and using potential (12), we can obtain the three-soliton solution of (6). The expression is so complicated; therefore we do not give out it in detail. The three-soliton solution is shown in Figure 3, from which wave chase phenomenon also can be found.
Three-soliton waves and with α=1, β=-1, w1=1, w2=1.1, w3=0.9, b1=1, b2=1, b3=1, k=1, λ1=0, λ2=0, λ3=0. (a) Spatial plots in the intervals x∈[-10,20] and t∈[-15,15]; (b)–(d) plan plots and with t=-15,0,15, x∈[-5,30], x∈[-10,20], x∈[-20,15].
When N⩾4, similar computational work becomes more and more complicated since the coefficients of the exponential functions become a highly nonlinear system as shown in [23]. Fortunately, we can find a uniform formula of the N-soliton solutions by analyzing the obtained solutions (17), (22), and (24). We rewrite solutions (17), (22), and (24) in an alternative form:
(27)u(x,t)=-2[ln(1+b1eξ1)]xx,
where ξ1=w1x-c1t+λ1;
(28)u(x,t)=-2[ln(1+b1eξ1+b2eξ2+b1b2A12eξ1+ξ2)]xx,
where ξi=wix-cit+λi (i=1,2), A12=(w1-w2)2/(w1+w2)2; and
(29)u(x,t)=-2[hhhhiihh+b1b2b3A12A13A23eξ1+ξ2+ξ3)ln(1+b1eξ1+b2eξ2+b3eξ3hhhhhhhi+b1b2A12eξ1+ξ2+b1b3A13eξ1+ξ3hhhhhhhi+b2b3A23eξ2+ξ3hhhhiihh+b1b2b3A12A13A23eξ1+ξ2+ξ3)]xx,
where ξi=wix-cit+λi (i=1,2,3), Aij=((wi-wj)2/(wi+wj)2) (1⩽i<j⩽3).
The uniform formula of the N-soliton solutions can be constructed as follows:
(30)u(x,t)=-2[ln(1+∑i=1Nbieξi+∑1⩽i<j⩽NbibjAijeξi+ξjhhhhhhihhh+⋯+∏i=1Nbi∏1⩽i<j⩽NAijeξ1+ξ2+⋯ξN)]xx,
where ξi=wix-cit+λi (i=1,2,…,N) and Aij=((wi-wj)2/(wi+wj)2) (1⩽i<j⩽N).
4. Conclusions
In this paper, one-soliton, two-soliton, and three-soliton solutions of the ACH-KdV equation have been successfully obtained, from which the uniform formula of N-soliton solutions is derived. This is due to the generalization of the Exp-function method. Figures 1–3 imply that these obtained solutions have rich local structures. It may be important to explain some physical phenomena. The method with the help of mathematical software for generating 1-soliton, 2-soliton, and 3-soliton solutions is more simple and straightforward than Hirota bilinear method, without employing the bilinear operator defined in the Hirota bilinear method. The paper shows that the Exp-function method may provide us with a straightforward and effective mathematical tool for generating N-soliton solutions or testing its existence and can be extended to other NLEEs in mathematical physics.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research is supported by the Natural Science Foundation of China (nos. 11161020, 11361023) as well as the Young and Middle-Aged Academic Backbone of Honghe University (no. 2014GG0105).
AblowitzM. J.SegurH.Solitons, nonlinear evolution equations and inverse scattering1992244721725MatveevV. B.SalleM. A.1991Berlin, GermanySpringer10.1007/978-3-662-00922-2MR1146435HirotaR.Exact solution of the modified Korteweg-de Vries equation for multiple collisions of solitons19723351456145810.1143/JPSJ.33.14562-s2.0-0001508211MingliangW.ZhouY.LiZ.Application of a homogeneous balance method to exact solutions of nonlinear equations in mathematical physics19962161–567752-s2.0-30244524644LiuS.FuZ.ZhaoQ.Jacobi elliptic function expansion method and periodic wave solutions of nonlinear wave equations20012891-2697410.1016/S0375-9601(01)00580-1MR18620822-s2.0-0035828886ParkesE. J.DuffyB. R.An automated tanh-function method for finding solitary wave solutions to non-linear evolution equations199698328830010.1016/0010-4655(96)00104-X2-s2.0-0030284901HeJ.-H.WuX.-H.Exp-function method for nonlinear wave equations200630370070810.1016/j.chaos.2006.03.020MR22386952-s2.0-33745177020Mohamad JawadA. J.PetkovićM. D.BiswasA.Modified simple equation method for nonlinear evolution equations2010217286987710.1016/j.amc.2010.06.030MR26786032-s2.0-77955664465FanE.Uniformly constructing a series of explicit exact solutions to nonlinear equations in mathematical physics200316581983910.1016/S0960-0779(02)00472-1MR1964991ZBL1030.351362-s2.0-0037408898WangM.LiX.Applications of F-expansion to periodic wave solutions for a new Hamiltonian amplitude equation20052451257126810.1016/j.chaos.2004.09.044MR2123272ZBL1092.370542-s2.0-13444310842WangD.ZhangH.-Q.Further improved F-expansion method and new exact solutions of Konopelchenko-Dubrovsky equation200525360161010.1016/j.chaos.2004.11.026MR21323592-s2.0-14944374499YombaE.The extended Fan's sub-equation method and its application to KdV-MKdV, BKK and variant Boussinesq equations2005336646347610.1016/j.physleta.2005.01.027MR21187752-s2.0-14144251108CamassaR.HolmD. D.An integrable shallow water equation with peaked solitons1993711116611664CamassaR.HolmD. D.HymanJ. M.A new integrable shallow water equation19943113310.1016/S0065-2156(08)70254-02-s2.0-56349110516SchiffJ.The Camassa-Holm equation: a loop group approach19981211-2244310.1016/S0167-2789(98)00099-2MR16443862-s2.0-0002632382FisherM.SchiffJ.The Camassa Holm equation: conserved quantities and the initial value problem1999259537137610.1016/S0375-9601(99)00466-1MR17096862-s2.0-0001500560ParkerA.On the Camassa-Holm equation and a direct method of solution. II . Soliton solutions200546120633611363210.1098/rspa.2005.1536MR21712802-s2.0-28544450796ParkerA.On the Camassa-Holm equation and a direct method of solution. III. N-soliton solutions200546120643893391110.1098/rspa.2005.1537MR21860112-s2.0-28544439875QiaoZ. J.ZhangG. P.On peaked and smooth solitons for the Camassa-Holm equation200673565766310.1209/epl/i2005-10453-yMR22164352-s2.0-33644926968HolmD. D.IvanovR. I.Smooth and peaked solitons of the CH equation201043431843400310.1088/1751-8113/43/43/434003MR27277772-s2.0-78649678927HoneA. N. W.The associated Camassa-Holm equation and the KdV equation19993227L307L31410.1088/0305-4470/32/27/103MR17233012-s2.0-0033538373LuoL.QiaoZ.LopezJ.Integrable generalization of the associated Camassa-Holm equation2014378967768310.1016/j.physleta.2013.11.029MR3158737MarinakisV.The exp-function method and n-soliton solutions20086310-116536562-s2.0-58449099032