MPEMathematical Problems in Engineering1563-51471024-123XHindawi Publishing Corporation24936110.1155/2009/249361249361Research ArticleSolitons, Peakons, and Periodic Cuspons of a Generalized Degasperis-Procesi EquationZhouJiangboTianLixinMacauElbert E. NeherNonlinear Scientific Research CenterFaculty of ScienceJiangsu UniversityZhenjiang, Jiangsu 212013Chinaujs.edu.cn200930032009200924112008230220092009Copyright © 2009This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

We employ the bifurcation theory of planar dynamical systems to investigate the exact travelling wave solutions of a generalized Degasperis-Procesi equation utuxxt+4uux+γ(uuxx)x=3uxuxx+uuxxx. The implicit expression of smooth soliton solutions is given. The explicit expressions of peaked soliton solutions and periodic cuspon solutions are also obtained. Further, we show the relationship among the smooth soliton solutions, the peaked soliton solutions, and the periodic cuspon solutions. The physical relevance of the found solutions and the reason why these solutions can exist in this equation are also given.

1. Introduction

Recently, Degasperis and Procesi  derived a nonlinear dispersive equation ut-uxxt+4uux=3uxuxx+uuxxx which is called the Degasperis-Procesi equation. Here u(t,x) represents the fluid velocity at time t in the x direction in appropriate nondimensional units (or, equivalently the height of the water's free surface above a flat bottom). The nonlinear convection term uux in (1.1) causes the steepening of wave form, whereas the nonlinear dispersion effect term 3uxuxx+uuxxx=((1/2)u2)xxx in (1.1) makes the wave form spread. Equation (1.1) can be regarded as a model for nonlinear shallow water dynamics. Degasperis et al.  showed that the (1.1) is integrable by deriving a Lax pair and a bi-Hamiltonian structure for it. Yin proved local well posedness to (1.1) with initial data u0Hs(), s>3 on the line  and on the circle . The global existence of strong solutions and weak solutions to (1.1) is investigated in . The solution to Cauchy problem of (1.1) can also blow up in finite time when the initial data satisfies certain sign condition. Vakhnenko and Parkes  obtained periodic and solitary-wave solutions of (1.1). Matsuno [12, 13] obtained multisoliton, cusp and loop soliton solutions of (1.1). Lundmark and Szmigielski  investigated multipeakon solutions of (1.1). Lenells  classified all weak travelling wave solutions. The shock wave solutions of (1.1) are investigated in [16, 17].

Yu and Tian  investigated the following generalized Degasperis-Procesi equation: ut-uxxt+4uux=3uxuxx+uuxxx-γuxxx, where γ is a real constant, and the term uxxx denotes the linear dispersive effect. They obtained peaked soliton solutions and period cuspon solutions of (1.2). Unfortunately, they did not obtain smooth soliton solutions of (1.2).

In this paper, we are interesting in the following generalized Degasperis-Procesi equation: ut-uxxt+4uux+γ(u-uxx)x=3uxuxx+uuxxx, where γ is a real constant, the term ux denotes the dissipative effect and the term uxxx represents the linear dispersive effect. Employing the bifurcation theory of planar dynamical systems, we obtain the analytic expressions of smooth solitons, peaked solitons, and period cuspons of (1.3). Our work covers and supplements the results obtained in .

The remainder of the paper is organized as follows. In Section 2, using the travelling wave transformation, we transform (1.3) into the planar dynamical system (2.3) and then discuss bifurcations of phase portraits of system (2.3). In Section 3, we obtain the implicit expression of smooth solitons, the explicit expressions of peaked solitons and periodic cuspon solutions. At the same time, we show that the limits of smooth solitons and periodic cusp waves are peaked solitons. In Section 4, we discuss the physical relevance of the found solutions and give the reason why these solutions can exist in (1.3).

2. Bifurcations of Phase Portraits of System (2.3)

We look for travelling wave solutions of (1.3) in the form of u(x,t)=φ(x-ct)=φ(ξ), where c is the wave speed and ξ=x-ct. Substituting u=φ(ξ) into (1.3), we obtain -cφ+cφ+4φφ-3φφ-φφ+γφ-γφ=0.

By integrating (2.1) once we have φ(φ-c+γ)=g-(c-γ)φ+2φ2-(φ)2, where g is the integral constant.

Let y=φ, then we get the following planar dynamical system: dφdξ=y,dydξ=g-(c-γ)φ+2φ2-y2φ-c+γ, with a first integral H(φ,y)=(φ-c+γ)2(y2-φ2-g)=h, where h is a constant.

Note that (2.3) has a singular line φ=c-γ. To avoid the line temporarily we make transformation dξ=(φ-c+γ)dζ. Under this transformation, (2.3) becomes dφdζ=(φ-c+γ)y,dydζ=g-(c-γ)φ+2φ2-y2.

System (2.3) and system (2.5) have the same first integral as (2.4). Consequently, system (2.5) has the same topological phase portraits as system (2.3) except for the straight line φ=c-γ. Obviously, φ=c-γ is an invariant straight-line solution for system (2.5).

For a fixed h, (2.4) determines a set of invariant curves of system (2.5). As h is varied, (2.4) determines different families of orbits of system (2.5) having different dynamical behaviors. Let M(φe,ye) be the coefficient matrix of the linearized system of (2.5) at the equilibrium point (φe,ye), then

M(φe,ye)=(yeφe-c+γ4φe-c+γ-2ye), and at this equilibrium point, we have J(φe,ye)=detM(φe,ye)=-2ye2-(φe-c+γ)(4φe-c+γ),p(φe,ye)=trace(M(φe,ye))=-ye. By the qualitative theory of differential equations (see ), for an equilibrium point of a planar dynamical system, if J<0, then this equilibrium point is a saddle point; it is a center point if J>0 and p=0; if J=0 and the Poincaré index of the equilibrium point is 0, then it is a cusp.

By using the first integral value and properties of equilibrium points, we obtain the bifurcation curves as follows: g1(c)=(c-γ)28,g2(c)=-(c-γ)2. Obviously, the two curves have no intersection point and g2(c)<0<g1(c) for arbitrary constants cγ.

Using bifurcation method of vector fields (e.g., ), we have the following result which describes the locations and properties of the singular points of system (2.5).

Theorem 2.1.

For a given constant wave speed c0, let φ0±=c-γ±(c-γ)2-8g4for  gg1(c),y0±=±(c-γ)2+gfor  gg2(c). When c=γ,

if g<0, then system (2.5) has two equilibrium points (--g/2,0) and (-g/2,0), which are saddle points;

if g=0, then system (2.5) has only one equilibrium point (0,0), which is a cusp;

if g>0, then system (2.5) has two equilibrium points (0,-g) and (0,g), which are saddle points.

When cγ,

if g<g2(c), then system (2.5) has two equilibrium points (φ0-,0) and (φ0+,0). They are saddle points:

if c>γ, then φ0-<-1/2(c-γ)<1/4(c-γ)<c-γ<φ0+,

if c<γ, then φ0-<c-γ<1/4(c-γ)<-1/2(c-γ)<φ0+;

if g=g2(c), then system (2.5) has three equilibrium points (φ0-,0), (φ0+,0), and (c-γ,0). (c-γ,0) is a cusp:

if c>γ, then φ0-=-(1/2)(c-γ)<(1/4)(c-γ)<c-γ=φ0+. (φ0-,0) is a saddle point, while (φ0+,0) is a degenerate center point,

if c<γ, then φ0-=c-γ<1/4(c-γ)<-1/2(c-γ)=φ0+. (φ0-,0) is a degenerate center point, while (φ0+,0) is a saddle point;

if g2(c)<g<g1(c), then system (2.5) has four equilibrium points (φ0-,0), (φ0+,0), (c-γ,y0-), and (c-γ,y0+). (c-γ,y0-) and (c-γ,y0+) are two saddle points:

if c>γ, then φ0-<1/4(c-γ)<φ0+c-γ. (φ0-,0) is a saddle point, while (φ0+,0) is a center point,

if c<γ, then c-γφ0-<1/4(c-γ)<φ0+. (φ0-,0) is a center point, while (φ0+,0) is a saddle point.

Specially, when g=0,

if c>γ, then the three saddle points (φ0-,0), (c-γ,y0-), and (c-γ,y0+) form a triangular orbit which encloses the center point (φ0+,0),

if c<γ, then the three saddle points (φ0+,0), (c-γ,y0-), and (c-γ,y0+) form a triangular orbit which encloses the center point (φ0-,0);

if g=g1(c), then system (2.5) has three equilibrium points ((c-γ)/4,0), (c-γ,y0-), and (c-γ,y0+). ((c-γ)/4,0) is a degenerate center point, while (c-γ,y0-) and (c-γ,y0+) are two saddle points;

if g>g1(c), then system (2.5) has two equilibrium points (c-γ,y0-) and (c-γ,y0+), which are saddle points.

Corresponding to the case c=γ and the case cγ, we show the phase portraits of system (2.5) in Figures 1 and 2, respectively.

The phase portraits of system (2.5)(c=γ). (a) g<0; (b) g=0; (c) g>0.

The phase portraits of system (2.5) (cγ). (a) g<g2(c),c>γ; (b) g<g2(c),c<γ; (c) g=g2(c),c>γ; (d)g=g2(c),c<γ; (e) g2(c)<g<0,c>γ; (f) g2(c)<g<0,c<γ; (g) g=0,c>γ; (h) g=0,c<γ; (i) 0<g<g1(c),c>γ; (j) 0<g<g1(c),c<γ; (k) g=g1(c),c>γ; (l) g=g1(c),c<γ; (m)g>g1(c),c>γ; (n)g>g1(c),c<γ.

3. Solitons, Peakons, and Periodic Cusp Wave SolutionsTheorem 3.1.

Given arbitrary constant cγ, let ξ=x-ct, then

when 0<g<g1(c),

if c>γ, then (1.3) has the following smooth hump-like soliton solutions: β1(φ1+)=β1(φ)e-|ξ|for  φ0-<φ<φ1+,

if c<γ, then (1.3) has the following smooth valley-like soliton solutions: β2(φ)=β2(φ1-)e-|ξ|for  φ1-<φ<φ0+,

when g=0, (1.3) has the following peaked soliton solutions: φ=(c-γ)e-|ξ|,

when g2(c)<g<0,

if c>γ, then (1.3) has the following periodic cusp wave solutions: u(x,t)=φ3(x-ct-2nT)for  (2n-1)T<x-ct<(2n+1)T,

if c<γ, then (1.3) has the following periodic cusp wave solutions: u(x,t)=φ4(x-ct-2nT)for  (2n-1)T<x-ct<(2n+1)T,

where β1(φ)=(2φ2+l1φ+l2+2φ+l1)(φ-φ0-)α1(2a1φ2+l1φ+l2+b1φ+l3)α1,β2(φ)=(2φ2+m1φ+m2+2φ+m1)(φ-φ0+)α2(2a2φ2+m1φ+m2+b2φ+m3)α2,l1=-3(c-γ)+(c-γ)2-8g2,l2=3(c-γ)2-4g+5(c-γ)(c-γ)2-8g8,l3=(c-γ)2-4g+3(c-γ)(c-γ)2-8g2,m1=-3(c-γ)-(c-γ)2-8g2,m2=3(c-γ)2-4g-5(c-γ)(c-γ)2-8g8,m3=(c-γ)2-4g-3(c-γ)(c-γ)2-8g2,a1=(c-γ)2-8g+3(c-γ)(c-γ)2-8g4,a2=(c-γ)2-8g-3(c-γ)(c-γ)2-8g4,b1=-(c-γ)-(c-γ)2-8g,b2=-(c-γ)+(c-γ)2-8g,α1=-3(c-γ)+(c-γ)2-8g2(c-γ)2-8g+3(c-γ)(c-γ)2-8g,α2=-3(c-γ)+(c-γ)2-8g2(c-γ)2-8g-3(c-γ)(c-γ)2-8g,φ3(ξ)=l+e-|ξ|+l-e|ξ|for  -gφ3c-γ,φ4(ξ)=l+e|ξ|+l-e-|ξ|for  c-γφ4--g,l±=c-γ±(c-γ)2+g2,T=|ln(-g+-2g)-ln(2l-)|.φ1±=34(c-γ)±14(c-γ)2-8g12(c-γ)2(c-γ)(c-γ)2-8g,φ0+ and φ0- are as in (2.9).

Before proving this theorem, we take a set of data and employ Maple to display the graphs of smooth solion, peaked soliton and periodic cuspon solutions of (1.3), see Figures 3, 4, 5, 6 and 7.

Smooth hump-like soliton solutions of (1.3)(c=2, γ=1). (a) g=0.12; (b) g=0.1; (c) g=0.075; (d) g=0.05.

Smooth valley-like soliton solutions of (1.3) (c=1, γ=2). (a) g=0.12; (b) g=0.1; (c) g=0.075; (d) g=0.05.

Peaked soliton solutions of (1.3). (a) c=2, γ=1; (b) c=1, γ=2.

Period cuspon solutions of (1.3) (c=2, γ=1). (a) g=-0.4; (b) g=-0.1; (c) g=-0.01; (d) g=-0.0000001.

Period cuspon solutions of (1.3) (c=1, γ=2). (a) g =-0.4; (b) g=-0.1; (c) g=-0.01; (d) g=-0.0000001.

Proof.

Usually, a solion solution of (1.3) corresponds to a homoclinic orbit of system (2.5), and a periodic travelling wave solution of (1.3) corresponds to a periodic orbit of system (2.5). The graphs of homoclinic orbit, periodic orbit of system (2.5), and their limit curves are shown in Figure 8.

(1) When 0<g<g1(c),c>γ, system (2.5) has a homoclinic orbit (see Figure 8(a)). This homoclinic orbit can be expressed as

y=±(φ-φ0-)φ2+l1φ+l2c-γ-φfor  φ0-<φ<φ1+. Substituting (3.7) into the first equation of system (2.3) and integrating along this homoclinic orbit, we obtain (3.1).

When 0<g<g1(c),c<γ, we can obtain (3.2) in similar way.

(2) When g=0,c>γ, system (2.5) has a homoclinic orbit that consists of the following three line segments (see Figure 8(c)):

y=±φfor  φ0-φφ1+,φ=c-γfor  -(c-γ)2+gy(c-γ)2+g. Substituting (3.8) into the first equation of system (2.3) and integrating along this orbit, we obtain (3.3).

When g=0,c<γ, we can also obtain (3.3).

(3) When g2(c)<g<0,c>γ, system (2.5) has a periodic orbit (see Figure 8(e)). This periodic orbit can be expressed as

y=±φ2+gfor  -gφc-γ,φ=c-γfor  -(c-γ)2+gy(c-γ)2+g. Substituting (3.10) into the first equation of system (2.3) and integrating along this periodic orbit, we obtain (3.4).

When g2(c)<g<0,c<γ, we can obtain (3.5).

The orbits of system (2.5) connecting with the saddle points. (a) 0<g<g1(c),c>γ; (b) 0<g<g1(c),c<γ; (c) g=0,c>γ; (d) g=0,c<γ; (e) g2(c)<g<0,c>γ; (f) g2(c)<g<0,c<γ.

Remark 3.2.

From the above discussion, we can see that when g<0,g0, the period of the periodic cusp wave solution becomes bigger and bigger, and the periodic cuspon solutions (3.4) and (3.5) tend to the peaked soliton solutions (3.3). When g>0,g0, the smooth hump-like soliton solutions (3.1) and the smooth valley-like soliton solutions (3.2) lose their smoothness and tend to the peaked soliton solutions (3.3).

4. Discussion

In this paper, we obtain the solitons, peakons, and periodic cuspons of a generalized Degasperis-Procesi equation (1.3). These solitons denote the nonlinear localized waves on the shallow water's free surface that retain their individuality under interaction and eventually travel with their original shapes and speeds. The balance between the nonlinear steepening and dispersion effect under (1.3) gives rise to these solitons.

The peakon travels with speed equal to its peak amplitude. This solution is nonanalytic, having a jump in derivative at its peak. Peakons are true solitons that interact via elastic collisions under (1.3). We claim that the existence of a singular straight line for the planar dynamical system (2.3) is the original reason why the travelling waves lose their smoothness.

Also, the periodic cuspon solution is nonanalytic, having a jump in derivative at its each cusp.

Acknowledgments

The authors are deeply grateful to the anonymous referee for the careful reading of our manuscript and many constructive comments and valuable suggestions, which have helped them to improve it. This research was supported by NSF of China (No. 10771088) and Startup Fund for Advanced Talents of Jiangsu University (No. 09JDG013).

DegasperisA.ProcesiM.DegasperisA.GaetaG.Asymptotic integrabilitySymmetry and Perturbation Theory (Rome, 1998)1999River Edge, NJ, USAWorld Scientific2337MR1844104ZBL0963.35167DegasperisA.KholmD. D.KhonA. N. I.A new integrable equation with peakon solutionsTheoretical and Mathematical Physics2002133214631474MR200153110.1023/A:1021186408422YinZ.On the Cauchy problem for an integrable equation with peakon solutionsIllinois Journal of Mathematics2003473649666MR2007229ZBL1061.35142YinZ.Global weak solutions for a new periodic integrable equation with peakon solutionsJournal of Functional Analysis20042121182194MR206524110.1016/j.jfa.2003.07.010ZBL1059.35149YinZ.Global solutions to a new integrable equation with peakonsIndiana University Mathematics Journal200453411891209MR209545410.1512/iumj.2004.53.2479ZBL1062.35094YinZ.Global existence for a new periodic integrable equationJournal of Mathematical Analysis and Applications20032831129139MR199417710.1016/S0022-247X(03)00250-6ZBL1033.35121ZhouY.Blow-up phenomenon for the integrable Degasperis-Procesi equationPhysics Letters A20043282-3157162MR207665310.1016/j.physleta.2004.06.027ZBL1134.37361EscherJ.LiuY.YinZ.Global weak solutions and blow-up structure for the Degasperis-Procesi equationJournal of Functional Analysis20062412457485MR227192710.1016/j.jfa.2006.03.022ZBL1126.35053YinZ.Well-posedness, blowup, and global existence for an integrable shallow water equationDiscrete and Continuous Dynamical Systems. Series A2004112-3393411MR2083424ZBL1061.3512310.3934/dcds.2004.11.393LiuY.YinZ.Global existence and blow-up phenomena for the Degasperis-Procesi equationCommunications in Mathematical Physics20062673801820MR224979210.1007/s00220-006-0082-5ZBL1131.35074VakhnenkoV. O.ParkesE. J.Periodic and solitary-wave solutions of the Degasperis-Procesi equationChaos, Solitons & Fractals200420510591073MR202937910.1016/j.chaos.2003.09.043ZBL1049.35162MatsunoY.Multisoliton solutions of the Degasperis-Procesi equation and their peakon limitInverse Problems200521515531570MR2173410ZBL1086.35095MatsunoY.Cusp and loop soliton solutions of short-wave models for the Camassa-Holm and Degasperis-Procesi equationsPhysics Letters A20063595451457MR225351910.1016/j.physleta.2006.06.065LundmarkH.SzmigielskiJ.Multi-peakon solutions of the Degasperis-Procesi equationInverse Problems200319612411245MR2036528ZBL1041.35090LenellsJ.Traveling wave solutions of the Degasperis-Procesi equationJournal of Mathematical Analysis and Applications200530617282MR213288910.1016/j.jmaa.2004.11.038ZBL1068.35163LundmarkH.Formation and dynamics of shock waves in the Degasperis-Procesi equationJournal of Nonlinear Science2007173169198MR231484710.1007/s00332-006-0803-3EscherJ.LiuY.YinZ. Y.Shock waves and blow-up phenomena for the periodic Degasperis-Procesi equationIndiana University Mathematics Journal20075618711710.1512/iumj.2007.56.3040ZBL1124.35041MR2305931YuL.TianL.WangX.The bifurcation and peakon for Degasperis-Procesi equationChaos, Solitons & Fractals2006304956966MR225026610.1016/j.chaos.2005.08.152ZBL1142.35589LuoD.WangX.ZhuD.HanM.Bifurcation Theory and Methods of Dynamical Systems199715River Edge, NJ, USAWorld Scientificxiii+461Advanced Series in Dynamical SystemsMR2001522ZBL0961.37015