arxiv: v1 [nlin.ps] 3 Sep 2009

Transcription

1 Soliton, kink and antikink solutions o a 2-component o the Degasperis-Procesi equation arxiv: v1 [nlin.ps] 3 Sep 2009 Jiangbo Zhou, Liin Tian, Xinghua Fan Nonlinear Scientiic Research Center, Facult o Science, Jiangsu Universit, Zhenjiang, Jiangsu , China Abstract In this paper, we emplo the biurcation theor o planar dnamical sstems to investigate the traveling wave solutions o a 2-component o the Degasperis-Procesi equation. The epressions or smooth soliton, kink and antikink solutions are obtained. Ke words: 2-component o the Degasperis-Procesi equation, biurcation method, soliton, kink solution, antikink solution 2000 MSC: 35Q51, 34C23, 37G10, 35Q35 1. Introduction Since the theor o solitons has ver wide applications in luid dnamics, nonlinear optics, biochemistr, microbiolog, phsics and man other ields, the stud o soliton solutions has become one o the important issues o nonlinear partial dierential equations [1]-[11]. In 1999, Degasperis and Procesi [12] derived a nonlinear dispersive equa- Corresponding author address: zhoujiangbo@ahoo.cn (Jiangbo Zhou) Ma 28, 2018

2 tion u t u t +4uu = 3u u +uu, (1.1) which is called the Degasperis-Procesi equation. Here u(t, ) represents the luid velocit at time t in the direction in appropriate nondimensional units (or, equivalentl the height o the water s ree surace above a lat bottom). The nonlinear convection term uu in Eq.(1.1) causes the steepening o wave orm, whereas the nonlinear dispersion eect term 3u u +uu = ( 1 2 u2 ) in Eq.(1.1) makes the wave ormspread. Eq.(1.1) can be regarded as a model or nonlinear shallow water dnamics [13]. An important issue regarding the Degasperis-Procesi equation is to ind its traveling wave solutions. Vakhnenko and Parkes [14] derived periodic and solitar wave solutions o Eq.(1.1). Matsuno [15]-[17] obtained multisoliton, N-soliton, cusp and loop soliton solutions o Eq.(1.1). Lundmark and Szmigielski [18] investigated multi-peakon solutions o Eq.(1.1). Shock wave solutions o Eq.(1.1) were obtained in [19]. Lenells [20] classiied all its weak traveling wave solutions. Chen and Tang [21] showed that Eq.(1.1) has kink-like and antikink-like wave solutions. Qiao [8] obtained the peakons, dehisced solitons, cuspons and new 1-peak solitons o Eq.(1.1). It is known that the Degasperis-Procesi equation has solitons, but has no kink or antikink solutions. In this paper, we generalize the Degasperis-Procesi equation to the ollowing 2-component o Degasperis-Procesi equation m t +um +3u m = v, (1.2) v t = (uv), where m = u u. This 2-component o the Degasperis-Procesi equation can be applied to describe shallow water waves. In Eqs.(1.2), u(, t) denotes 2

3 the height o the water surace above a horizontal bottom and v(,t) is related to the horizontal velocit ield. Obviousl, i v 0, then Eqs.(1.2) is reduced to the Degasperis-Procesi equation (1.1). Using the biurcation theor o planar dnamical sstems (see [22]), we will show that, in addition to smooth solitons, Eqs.(1.2) has kink and antikink solutions. The solitons denote the nonlinear localized waves on the shallow water s ree surace that retain their individualit under interaction and eventuall travel with their original shapes and speeds. The kink and antikink waves, which rise or descend rom one asmptotic state to another, orm another tpe. The balance between the nonlinear steepening and dispersion eect under Eqs.(1.2) gives rise to these solutions. The solutions presented in this paper ma help people to know deepl the described phsical process o Eqs.(1.2). The remainder o the paper is organized as ollows. In Section 2, using the traveling wave transormation, we transorm Eqs.(1.2) into a planar dnamical sstem and then discuss biurcations o phase portraits o this sstem. In Section 3, we obtain the epressions or the smooth soliton, kink and antikink solutions o Eqs.(1.2). A short conclusion is given in Section Biurcation and phase portraits o the traveling wave sstem Eqs.(1.2) also takes the orm u t u t +4uu = 3u u +uu +v, v t = (uv). (2.1) Let ξ = +ct, where c is the wave speed. B using the traveling wave transormation u(,t) = ϕ( + ct) = ϕ(ξ), v(,t) = ψ( + ct) = ψ(ξ), we 3

4 reduce sstem (1.2) to the ollowing ordinar dierential equations: cϕ cϕ +4ϕϕ = 3ϕ ϕ +ϕϕ +ψ, cψ = (ϕψ). Integrating (2.2) once with respect to ξ, we have cϕ cϕ = 2ϕ 2 (ϕ ) 2 ϕϕ +ψ +g, cψ = ϕψ +h, (2.2) (2.3) where g, h are two integral constants. From the second epression in sstem (2.3), we can obtain that ψ = h ϕ c. (2.4) Substituting (2.4) into the irst epression in sstem (2.3) ields ϕ = (2ϕ2 cϕ (ϕ ) 2 +g)(ϕ c) h (ϕ c) 2. (2.5) Let = ϕ, then we get the ollowing planar dnamical sstem: dϕ dξ =, d dξ = (2ϕ2 cϕ 2 +g)(ϕ c) h (ϕ c) 2. This sstem has a irst integral (2.6) H(ϕ,) = ϕ 4 2cϕ 3 +(c 2 +g)ϕ 2 2(cg +h)ϕ 2 (ϕ c) 2 = k,(2.7) where k is a constant. Note that (2.6) has a singular line ϕ = c. To avoid the line temporaril, wemaketransormationdξ = (ϕ c) 2 dζ. Underthistransormation, Eq.(2.6) becomes dϕ dζ = (ϕ c)2, d dζ = (2ϕ2 cϕ 2 +g)(ϕ c) h. 4 (2.8)

5 Sstem (2.6) and sstem (2.8) have the same irst integral as (2.7). Consequentl, sstem (2.8) has the same topological phase portraits as sstem (2.6) ecept or the straight line ϕ = c. Foraied k, (2.7)determines aset oinvariant curves osstem (2.8). As k is varied, (2.7) determines dierent amilies o orbits o sstem (2.8) having dierent dnamical behaviors. Let M(ϕ e, e ) be the coeicient matri o the linearized version o sstem (2.8) at the equilibrium point (ϕ e, e ), then 2 M(ϕ e, e ) = e (ϕ e c) (ϕ e c) 2 (2.9) 6ϕ 2 e 6cϕ e +c 2 +g 2 2 e (ϕ e c) and at this equilibrium point, we have J(ϕ e, e ) = detm(ϕ e, e ) = 3e(ϕ 2 e c) 2 (6ϕ 2 e 6cϕ e +c 2 +g)(ϕ e c) 2,(2.10) p(ϕ e, e ) = trace(m(ϕ e, e )) = 0. (2.11) B the qualitative theor o dierential equations (see [22]), or an equilibrium point o a planar dnamical sstem, i J < 0, then this equilibrium point is a saddle point; it is a center point i J > 0 and p = 0; i J = 0 and the Poincaré inde o the equilibrium point is 0, then it is a cusp. Noting that it is impossible to igure out the equilibrium points o sstem (2.8) or arbitrar constants c, g and h, so here or the sake o convenience, we set h = 3c 32 (c2 8g) 0. B using the irst integral value and properties o equilibrium points, we obtain the biurcation curves as ollows: g 1 (c) = (16+7 7)c 2, (2.12) 8 g 2 (c) = c2 8, (2.13) 5

6 g 3 (c) = 3c2 8, (2.14) g 4 (c) = 13c2 32. (2.15) Obviousl, the above our curves have no intersection point and g 1 (c) < g 2 (c) < g 3 (c) < g 4 (c) or arbitrar constant c 0. Using biurcation method o vector ields(e.g., [20]), we have the ollowing result which describes the locations and properties o the equilibrium points o sstem (2.8). Theorem 2.1. For given an constant wave speed c 0, let ϕ 1± = 1 8 (5c± 13c 2 32g) or g g 4 (c). (2.16) Then we have: (1) I g < g 2 (c), then sstem (2.8) has three equilibrium points (ϕ 1,0), (ϕ 1+,0) and ( c 4,0). (ϕ 1,0) and (ϕ 1+,0) are two saddle points, while ( c 4,0) is a center point. Speciall, when g 1 (c) < g < g 2 (c), sstem (2.8) has a homoclinic orbit, which connects with the saddle point (ϕ 1,0). (2) I g 2 (c) < g < g 4 (c), then sstem (2.8) has three equilibrium points (ϕ 1,0), (ϕ 1+,0) and ( c 4,0). (c 4,0) and (ϕ 1+,0) are two saddle points, while (ϕ 1,0) is a center point. Speciall, (i) When g 2 (c) < g < g 3 (c), sstem (2.8) has a homoclinic orbit, which connects with the saddle point ( c 4,0). (ii) When g = g 3 (c), sstem (2.8) has two heteroclinic orbits. The connect with the saddle points ( c,0) and (3c,0). 4 4 (iii) When g 3 (c) < g < g 4 (c), sstem (2.8) has a homoclinic orbit, which connects with the saddle point (ϕ 1+,0). 6

7 (3) I g = g 4 (c), then sstem (2.8) has two equilibrium points ( c,0) and 4 ( 5c,0). 8 (c,0) is a saddle point, while (5c,0) is a cusp. 4 8 (4) I g > g 4 (c), then sstem (2.8) has onl one equilibrium point ( c 4,0). It is a saddle point. Without loss o generalit, we show the phase portraits o sstem (2.8) in each region and on the biurcation curves in Fig.1 or the wave speed c > Soliton, kink and antikink solutions o Eqs.(1.2) Suppose that ϕ(ξ)(ξ = +ct) is a traveling wave solution o the irst component u in Eqs.(1.2) or ξ (,+ ), and lim ϕ(ξ) = A, lim ϕ(ξ) = ξ ξ + B, where A and B are two constants. I A = B, then ϕ(ξ) is called a soliton solution. I A B, then ϕ(ξ) is called a kink (or an antikink) solution. Usuall, a soliton solution o Eqs.(1.2) corresponds to a homoclinic orbit o sstem (2.6) and a kink (or an antikink) solution o Eq.(1.1) corresponds to a heteroclinic orbit (or the so-called connecting orbit) o sstem (2.6). In Fig.1(c), (d), (), the homoclinic orbits o sstem (2.6) or (2.8) can be epressed respectivel as = ± (ϕ ϕ 1 ) ϕ 2 +l 1 ϕ+l 2 ϕ c or ϕ 1 ϕ ϕ 2, (3.1) = ± (ϕ c 4 ) ϕ 2 +l 1 ϕ+l 2 ϕ c or c 4 ϕ ϕ 1, (3.2) = ± (ϕ ϕ 1+) ϕ 2 +l 1+ ϕ+l 2 ϕ c or ϕ 2+ ϕ ϕ 1+, (3.3) where l 1± = 1 4 ( 3c± 13c 2 32g), (3.4) 7

8 l 2± = 1 32 ( 9c2 +16g ±c 13c 2 32g), (3.5) l 1 = 3c 2, (3.6) l 2 = 3 16 c2 g, (3.7) ϕ 1 = 1 4 (3c 6c 2 16g), (3.8) ϕ 2± = 1 8 (3c±( 13c 2 32g +2 c(c+ 13c 2 32g))), (3.9) and ϕ 1± is as in (2.16). Substituting Eq.(3.1)-(3.3) into the irst equation o sstem (2.6), respectivel, and integrating along the corresponding homoclinic orbit, we have ϕ2 ϕ ϕ1 ϕ ϕ ϕ 2+ s c (s ϕ 1 ) s 2 +l 1 s+l 2 ds = ξ, (3.10) s c (s c 4 ) s 2 +l 1 s+l 2 ds = ξ, (3.11) s c (s ϕ 1+ ) s 2 +l 1+ s+l 2+ ds = ξ. (3.12) It ollows rom (3.10)-(3.12) that where β (ϕ 2 ) = β (ϕ)ep( ξ ), (3.13) β(ϕ 1 ) = β(ϕ)ep( ξ ), (3.14) β + (ϕ) = β + (ϕ 2+ )ep( ξ ), (3.15) β ± (ϕ) = (2 ϕ 2 +l 1± ϕ+l 2± +2ϕ+l 1± )(ϕ ϕ 1± ) α 1± (2, (3.16) a 1± ϕ2 +l 1± ϕ+l 2± +b 1± ϕ+l 3± ) α 1± 8

9 β(ϕ) = (2 ϕ 2 +l 1 ϕ+l 2 +2ϕ+l 1 )(ϕ c 4 )α 1 (2, (3.17) a 1 ϕ2 +l 1 ϕ+l 2 +b 1 ϕ+l 3 ) α 1 l 3± = 1 2 c2 2g, (3.18) l 3 = 2g, (3.19) a 1± = 1 16 (13c2 32g ±3c 13c 2 32g), (3.20) a 1 = 1 8 c2 +g, (3.21) b 1± = 1 2 (c± 13c 2 32g), (3.22) b 1 = c, (3.23) α 1± = 5c± 13c 2 32g 2 13c 2 32g 3c, (3.24) 13c 2 32g α 1 = c 2c2 +16g, (3.25) and l 1±, l 2±, l 1, l 2, ϕ 1± and ϕ 2± are as in (3.4), (3.5), (3.6), (3.7), (2.16) and (3.9), respectivel. (3.13)-(3.15) are the implicit epressions or the soliton solutions o the irst component u in Eqs.(1.2). From (2.4), (3.13)-(3.15), we can give the parametric epressions or the soliton solutions o the second component v in Eqs.(1.2) as ollows: ξ = ±(lnβ (ϕ) lnβ (ϕ 2 )), ψ = 3c(c2 8g) 32(ϕ c), (ϕ 1 ϕ ϕ 2 ), (3.26) 9

10 ξ = ±(lnβ(ϕ) lnβ (ϕ 1 )), ψ = 3c(c2 8g) c 32(ϕ c), ( 4 ϕ ϕ 1), (3.27) ξ = ±(lnβ + (ϕ 2+ ) lnβ + (ϕ)), ψ = 3c(c2 8g) 32(ϕ c), (ϕ 2+ ϕ ϕ 1+), (3.28) with ϕ as the parameter. Now we take a set o data and emplo Maple to displa the graphs o the above obtained soliton solutions in Figs In Fig.1(e), the heteroclinic orbits o sstem (2.6) or (2.8) can be epressed as = ± (ϕ c 3c )(ϕ ) 4 4 ϕ c or c 4 ϕ 3c 4, (3.29) Substituting Eq.(3.29) into the irst equation o sstem (2.6), and integrating along the heteroclinic orbits, we have (ϕ ± c 4 )3 2 = ( 3c 4 ϕ ±) 1 2 ep(±ξ) or c 4 ϕ ± 3c 4, (3.30) where ϕ + = ϕ + (ξ) is a kink solution, and ϕ = ϕ (ξ) is an antikink solution. The corresponding solutions o the second component v in Eqs.(1.2) can be epressed in parametric orm as ollows. ξ = ±( 3 ln(ϕ 2 ± c) ln(3c 4 ±)), ψ ± = c 3 16(ϕ ± c), c ( 4 ϕ ± 3c ), 4 (3.31) with ϕ ± as the parameter. Here ψ + = ψ + (ξ) is an antikink solution, and ψ = ψ (ξ) is a kink solution. Also, we take a set o data and emplo Maple to displa the graphs o the above obtained kink and antikink solutions in Fig

11 Remark 3.1. (1) The Degasperis-Procesi equation (1.1) does not have kink and antikink solutions because there is no heteroclinic orbit in the corresponding phase portraits (see [23]). (2) In[24], Vakhnenko and Parkes also suggested a generalized Degasperis- Procesi equation. Using direct integration, the obtain its kink and antikink solutions. So the connection between our suggested Eqs.(1.2) and the generalized Degasperis-Procesi equation considered in [24] needs urther studing. (3) The biurcation method onl enables us to analze the restricted class o solutions o Eqs.(1.2), namel, the one-valued solutions. However the Degasperis-Procesi equation has intriguing solutions such as loop-like solitar wave solutions (see [14]). So whether the 2-component o Degasperis-Procesi equation (1.2) has loop-like solutions needs other mathematical methods. 4. Conclusion In this work, we propose a generalization o the Degasperis-Procesi equation, that is, a 2-component o the Degasperis-Procesi equation (1.2) or the irst time. Emploing the biurcation method, we obtain the analtic epressions or smooth soliton, kink and antikink solutions o the 2-component o the Degasperis-Procesi equation (1.2). However, its integrabilit needs urther investigation. Reerences [1] A. M. Wazwaz, Compact and noncompact structures or a variant o KdV equation in higher dimensions, Appl. Math. Comput. 132 (2002)

12 [2] Z. Y. Yan, Modiied nonlinearl dispersive mk(m, n, k) equations: I. New compacton solutions and solitar pattern solutions, Comput. Phs. Comm. 1 (2003) [3] J. W. Shen, W. Xu, Biurcations o smooth and non-smooth traveling wave solutions in the generalized Camassa-Holm equation, Chaos, Solitons and Fractals 26 (2005) [4] J. H. He, X. H. Wu, Construction o solitar solution and compacton-like solution b variational iteration method, Chaos, Solitons and Fractals 29 (2006) [5] A. Parker, Cusped solitons o the Camassa-Holm equation. I. Cuspon solitar wave and antipeakon limit, Chaos, Solitons and Fractals 34 (2007) [6] Z. R. Liu, Z. R. Ouang, A note on solitar waves or modiied orms o Camassa-Holm and Degasperis-Procesi equations, Phs. Lett. A 366 (2007) [7] A. Biswas, 1-soliton solution o the K(m, n) equation with generalized evolution, Phs. Lett. A 372 (2008) [8] Z. J. Qiao, M-Shape peakons, dehisced solitons, cuspons and new 1- peak solitons or the Degasperis-Procesi equation, Chaos, Solitons and Fractals 37 (2008) [9] J. B. Zhou, L. X. Tian, Soliton solution o the osmosis K(2,2) equation, Phs. Lett. A 372 (2008)

13 [10] J. B. Li, Y. Zhan, Eact loop solutions, cusp solutions, solitar wave solutions and periodic wave solutions or the special CH-DP equation, Nonlinear Anal.: Real World Appl. 10 (2009) [11] J. B. Zhou, L. X. Tian, Solitons, peakons and periodic cusp wave solutions or the Fornberg-Whitham equation, Nonlinear Anal.: Real World Appl. (2009), in press. [12] A. Degasperis, M. Procesi, Asmptotic integrabilit, in: A. Degasperis, G. Gaeta (Eds.), Smmetr and Perturbation Theor, World Scientiic, Singapore, 1999, pp [13] A. Degasperis, D. D. Holm, A. N. W. Hone, A new integral equation with peakon solutions, Theo. Math. Phs. 133 (2002) [14] V. O. Vakhnenko, E. J. Parkes, Periodic and solitar-wave solutions o the Degasperis-Procesi equation, Chaos Solitons Fractals 20 (2004) [15] Y. Matsuno, Multisoliton solutions o the Degasperis-Procesi equation and their peakon limit, Inverse Problems 21 (2005) [16] Y. Matsuno, The N-soliton solution o the Degasperis-Procesi equation, Inverse Problems 21 (2005) [17] Y. Matsuno, Cusp and loop soliton solutions o short-wave models, or the Camassa-Holm and Degasperis-Procesi equations, Phs. Lett. A. 359 (2006)

14 [18] H. Lundmark, J. Szmigielski, Multi-peakon solutions o the Degasperis- Procesi equation, Inverse Problems 19 (2003) [19] H. Lundmark, Formation and dnamics o shock waves in the Degasperis-Procesi equation, J. Nonlinear Sci. 17 (2007) [20] J. Lenells, Traveling wave solutions o the Degasperis-Procesi equation, J. Math. Anal. Appl. 306 (2005) [21] C. Chen, M. Tang, A new tpe o bounded waves or Degasperis-Procesi equation, Chaos, Solitons and Fractals 27 (2006) [22] D. Luo, et al., Biurcation Theor and Methods o Dnamical Sstems, World Scientiic Publishing Co., London, [23] L. Q. Yu, L. X. Tian, X. D. Wang, The biurcation and peakon or Degasperis-Procesi equation, Chaos, Solitons and Fractals 30(2006) [24] V. O. Vakhnenko, E. J. Parkes, The solutions o a generalized Degasperis-Procesi equation, Rep. NAS Ukr. 8 (2006)

15 K2 K K10 K20 (a) K2 K K K20 (b) K K2 K4 (c) K K1.0 (d) K0.5 K K0.2 K K1 K2 K3 (e) () (g) K2 K4 K6 (h) Figure 1: The phase portraits o sstem (2.8) or the parameter c > 0. (a) g < g 1 (c); (b) g = g 1 (c); (c) g 1 (c) < g < g 2 (c); (d) g 2 (c) < g < g 3 (c); (e)g = g 3 (c); () g 3 (c) < g < g 4 (c); (g) g = g 4 (c); (h) g > g 4 (c). 15

16 K10 K K0.2 K0.4 (a) 2 K8 K6 K4 K (b) Figure 2: The soliton o Eqs.(1.2). (c = 1, g = 2, h = ) K15 K10 K (a) K15 K10 K K0.25 K0.30 K0.35 K0.40 K0.45 K0.50 K0.55 (b) Figure 3: The soliton o Eqs.(1.2). (c = 1, g = 0.37, h = ) K10 K (a) K10 K K0.45 K0.50 K0.55 K0.60 (b) Figure 4: The soliton o Eqs.(1.2). (c = 1, g = 0.4, h = ) 16

17 4 C K8 K6 K4 K (a) K8 K6 K4 K2 0 2 (b) C K0.3 K0.4 K0.5 K0.6 K0.7 K0.8 K0.9 4 K K4 K (c) K K K0.3 K0.4 K0.5 K0.6 K0.7 K0.8 K0.9 (d) Figure 5: The kink and antikink wave solutions o Eqs.(1.2). (c = 1, g = 0.375, h = ) 17