Peakon, Cuspon, Compacton, and Loop Solutions of a Three-Dimensional 3DKP(3, 2) Equation with Nonlinear Dispersion

Abstract

We study peakon, cuspon, compacton, and loop solutions for the three-dimensional Kadomtsev-Petviashvili equation (3DKP ( 3,2 ) equation) with nonlinear dispersion. Based on the method of dynamical systems, the 3DKP ( 3,2 ) equation is shown to have the parametric representations of the solitary wave solutions such as peakon, cuspon, compacton, and loop solutions. As a result, the conditions under which peakon, cuspon, compacton, and loop solutions appear are also given.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Academic Editor:Mohammad T. Darvishi

1, School of Mathematics and Computing Science, Guilin University of Electronic Technology, Guilin, Guangxi 541004, China

Received 15 April 2014; Accepted 13 July 2014; 13 August 2014

This 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.

1. Introduction

Nonlinear phenomena play a crucial role in applied mathematics and physics. Studies of various physical structures of nonlinear dispersive equations had attracted much attention in connection with the important problems that arise in scientific applications. Mathematically, these physical structures have been studied by using various powerful and efficient methods, such as inverse scattering method [1], Darboux transformation method [2, 3], Hirota bilinear method [4], Lie group method [5, 6], bifurcation method of dynamic systems [7, 8], tanh function method [9-12], Fan-expansion method [13, 14], and homogenous balance method [15]. Practically, there is no unified technique that can be employed to handle all types of nonlinear dispersive equations.

Recently, Xie and Yan [16] considered the following three-dimensional Kadomtsev-Petviashvili equation with nonlinear dispersion (the 3DKP ( m , n ) equations in short): [figure omitted; refer to PDF] where a is a nonzero real number.

When n = 1 , m = 2 , 3DKP ( 2,1 ) equation becomes the 3DKP equation [17, 18] [figure omitted; refer to PDF]

When n = 1 , m = 3 , 3DKP ( 3,1 ) equation becomes the 3DmKP equation [figure omitted; refer to PDF]

When a = 1 / 2 , m = n = 2 , Rosenau and Hyman [19] have given a compacton solution to 3DKP ( 2,2 ) equation [figure omitted; refer to PDF] By using some transformations, some compactons and solitary patterns are obtained [16]. Recently, by using the ansatz method and the Exp-function method, Inc [20] considered some compact and noncompact solutions for (1) and obtained a new traveling wave solution for the 3DKP ( 2,1 ) equation. The authors did not study the bifurcation behavior of the traveling wave solutions of the corresponding traveling wave equations in its parameter space. It is important to understand the dynamical behavior for the traveling wave solutions governed by a singular traveling wave equation. We emphasize that peakon, cuspon, and loop solutions have not been available as yet. To answer this question, we will consider the bifurcations of traveling wave solutions of 3DKP ( 3,2 ) equation in the three-parameter space ( α , r , β ) .

We now assume that m = 3 , n = 2 and make the transformations u ( x , y , z , t ) = [varphi] ( k ( x + l y + s z - λ t ) ) = [varphi] ( ξ ) , where k , l , s , λ ...0; 0 . Thus, (1) becomes [figure omitted; refer to PDF] where " [variant prime] " is the derivative with respect to ξ . Integrating (5) twice, we have [figure omitted; refer to PDF] where g is an integration constant, α = a / ^{k 2} , and β = ( 1 / a ) ( ^{l 2} + ^{s 2} - λ ) . We further make the transformation [figure omitted; refer to PDF] where r satisfies the equation ^{[varphi] 3} + β [varphi] + g = 0 . Then (6) is equivalent to the following two-dimensional systems: [figure omitted; refer to PDF] which admits the following first integral: [figure omitted; refer to PDF]

Systems (8) are planar dynamical systems defined in the 3-parameter space ( α , r , β ) . For a fixed α , we will investigate the bifurcations of phase portraits of (8) in the phase plane ( ψ , y ) as the parameters r , β are changed.

We emphasize that when ψ = r , the right-hand side of the second equation of systems (8) is discontinuous. We call such systems the singular traveling wave systems. The straight line ψ = r in the ψ - y -phase plane is called a singular straight line. There are some general theories and methods for investigating this type of singular traveling wave systems [21-25]. It is now well known that the existence of the singular straight lines implies the occurrence of some nonsmooth dynamical behaviors and curve breaking phenomena of the traveling wave solutions of such system, more precisely the so-called peakon and cuspon, and so forth.

It is important to understand the dynamical behavior for the traveling wave solutions governed by a singular traveling wave equation.

The rest of this paper is organized as follows. In Section 2, we discuss the bifurcations of phase portraits of (8), where explicit parametric conditions will be derived. In Section 3, we give explicit parametric representations for peakon, cuspon, compacton, and loop solutions of (1). Section 4 contains the concluding remarks.

2. Bifurcation Set and All Phase Portraits of System (8)

Throughout we assume that α < 0 . Otherwise, we can make a transformation [figure omitted; refer to PDF] to reduce (8) to this case. Based on ψ = r being a straight line solution to the system [figure omitted; refer to PDF] where d ξ = ( ψ - r ) d ζ , for ψ ...0; r , we say that the system (12) is the associated regular system of (8) [25].

Thus, system (12) has five equilibrium points O ( 0,0 ) , _{A ±} ( _{ψ ±} , 0 ) , and _{S ±} ( r , ± Y ) , where [figure omitted; refer to PDF]

Let M ( _{ψ e} , _{y e} ) be the coefficient matrix of the linearized system of (12) at an equilibrium point ( _{ψ e} , _{y e} ) and let J ( _{ψ e} , _{y e} ) be its Jacobin determinant. Then, we have [figure omitted; refer to PDF]

By the theory of planar dynamical systems, we know that for an equilibrium point ( _{ψ e} , _{y e} ) of a planar integrable system, if J < 0 , then the equilibrium point is a saddle point; if J > 0 and Trace ( M ( _{ψ e} , _{y e} ) ) = 0 , then it is a center point; if J > 0 and ( Trace ( M ( _{ψ e} , _{y e} ) ) ^{) 2} - 4 J ( _{ψ e} , _{y e} ) > 0 , then it is a node; and if J = 0 and the Poincare index of the equilibrium point is zero, then it is a cusp.

For the Hamiltonian defined by (10), we write that [figure omitted; refer to PDF]

The relations Δ = 0 , H ( r , ± Y ) = H ( 0,0 ) , and Y = 0 give the following three bifurcation curves: [figure omitted; refer to PDF] Thus, in the ( r , β ) -plane, we have 8 different parameter regions partitioned by the curves ( _{L i} ) , i = 1,2 , 3 , and r = 0 , which are shown in Figure 1.

Figure 1: The bifurcation set of system (8) in ( r , β ) -parameter plane for α < 0 .

[figure omitted; refer to PDF]

We use Figure 2 to show the bifurcations of the phase portraits of (12) for α < 0 .

The phase portraits of (8) for α < 0 .

(a) ( r , β ) ∈ _{A 1}

[figure omitted; refer to PDF]

(b) ( r , β ) ∈ _{A 2}

[figure omitted; refer to PDF]

(c) ( r , β ) ∈ _{A 3}

[figure omitted; refer to PDF]

(d) ( r , β ) ∈ _{A 4}

[figure omitted; refer to PDF]

(e) ( r , β ) ∈ _{A 5}

[figure omitted; refer to PDF]

(f) ( r , β ) ∈ _{A 6}

[figure omitted; refer to PDF]

(g) ( r , β ) ∈ _{A 7}

[figure omitted; refer to PDF]

(h) ( r , β ) ∈ _{A 8}

[figure omitted; refer to PDF]

(i) ( r , β ) ∈ _{L 1 +}

[figure omitted; refer to PDF]

(j) ( r , β ) ∈ _{L 1 -}

[figure omitted; refer to PDF]

(k) ( r , β ) ∈ _{L 2 +}

[figure omitted; refer to PDF]

(l) ( r , β ) ∈ _{L 2 -}

[figure omitted; refer to PDF]

(m) ( r , β ) ∈ _{L 3 +}

[figure omitted; refer to PDF]

(n) ( r , β ) ∈ _{L 3 -}

[figure omitted; refer to PDF]

3. Explicit Parametric Representations of Peakon, Cuspon, Compacton, and Loop Solutions of (8)

In this section, we give some parametric representations of peakon, cuspon, compacton, and loop solutions. To discuss the existence of peakon and cuspon solutions, we need to use the following two lemmas related to the singular straight line (see [25]).

Lemma 1 (the rapid-jump property of the derivative near the singular straight line).

Suppose that in a left (or right) neighborhood of a singular straight line there exists a family of periodic orbits. Then, along a segment of every orbit near the straight line, the derivative of the wave function jumps down rapidly on a very short time interval.

Lemma 2 (existence of finite time interval of solution with respect to wave variable in the positive or negative direction).

For a singular nonlinear traveling wave system of the first class with possible change of the wave variable, if an orbit transversely intersects with a singular straight line at a point or it approaches a singular straight line, but the derivative tends to infinity, then it only takes a finite time interval to make moved point of the orbit arrive on the singular straight line.

In the following, we give explicit parametric representations of peakon, cuspon, compacton, and loop solutions.

3.1. Peakon

(1) Suppose that α < 0 , ( r , β ) ∈ _{L 1 +} . In this case, we have the phase portrait of (8) shown in Figure 2(i). Notice that H ( _{ψ ±} , 0 ) = H ( - ( 3 / 2 ) r , 0 ) = H ( r , ± Y ) = ( 27 / 40 ) ^{r 5} . We see from (10) that two arch curves connecting _{S ±} in the left side of the straight line ψ = r have the algebraic equation [figure omitted; refer to PDF] Thus, by Lemma 2, we can take initial value as ψ ( 0 ) = 0 . Then, we have [figure omitted; refer to PDF] which is a a solitary peaked wave solution to (1) (so called "peakon" [26]). The profile of peakon soliton solution is shown in Figure 3(a).

The peakon solutions of (8) for α < 0 .

(a) ( r , β ) ∈ _{L 1 +}

[figure omitted; refer to PDF]

(b) ( r , β ) ∈ _{L 2 +}

[figure omitted; refer to PDF]

(2) Suppose that α < 0 , ( r , β ) ∈ _{L 2 +} . In this case, we have the phase portrait of (8) shown in Figure 2(k). Notice that H ( _{ψ ±} , 0 ) = H ( 0,0 ) = 0 . We see from (10) that two arch curves connecting _{S ±} in the left side of the straight line ψ = r have the algebraic equation [figure omitted; refer to PDF] Thus, by Lemma 2, we can take initial value as ψ ( 0 ) = r . Then, we have [figure omitted; refer to PDF] which is a peakon soliton solution to (1). The profile of peakon soliton solution is shown in Figure 3(b).

Therefore, we have the following.

Theorem 3.

(1) When the parameter groups α , r , β of system (8) satisfy the condition α < 0 with r > 0 , β = - ( 3 / 4 ) ^{r 2} , there exists a heteroclinic loop of system (12) given by three branches of the curves H ( ψ , y ) = ( 27 / 40 ) ^{r 5} ( = _{h s} ) . As the limit curves of a family of periodic orbits of system (8), the curve triangle (i.e., heteroclinic loop) in Figure 2(i) gives rise to a solitary peaked wave solution (a peakon) of (1), which has the exact parametric representation given by (18).

(2) When the parameter groups α , r , β of system (8) satisfy the condition α < 0 with r > 0 , β = - ( 24 / 5 ) ^{r 2} , there exists a heteroclinic loop of system (12) given by three branches of the curves H ( ψ , y ) = 0 . As the limit curves of a family of periodic orbits of system (8), the curve triangle (i.e., heteroclinic loop) in Figure 2(k) gives rise to a peakon soliton solution to (1), which has the exact parametric representation given by (20).

Remark 4.

To the best of our knowledge, solutions (18) and (20) obtained for (1) have not been reported in the literature.

3.2. Cuspon

(3) Suppose that α < 0 , ( r , β ) ∈ _{L 3 +} . In this case, we have the phase portrait of (8) shown in Figure 2(m). Notice that H ( 0,0 ) = 0 . We see from (10) that two arch curves connecting O ( 0,0 ) in the right side of the straight line ψ = r have the algebraic equation [figure omitted; refer to PDF] where _{ψ 0} < _{ψ -} and ^{δ 2} = ( - 1 / 4 _{ψ 0} ) ( ( 5 / 2 ) r ^{ψ 0 2} + ( 215 / 12 ) ^{r 2} _{ψ 0} + 30 ^{r 3} ) .

On the basis of Lemma 2, we can take initial value ψ ( 0 ) = r . Then, we have [figure omitted; refer to PDF] By introducing a new variable χ , from (22), one obtains the parametric representation of the following cuspon solution to (1) [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the elliptic integral of the third kind and sn ( u , k ) and dn ( u , k ) are the Jacobian elliptic functions [27]. According to (23), we may plot the graph of cuspon solution to (1) shown in Figure 4(a).

The cuspon and compacton solutions of (8) for α < 0 .

(a) ( r , β ) ∈ _{L 3 +}

[figure omitted; refer to PDF]

(b) ( r , β ) ∈ _{A 3}

[figure omitted; refer to PDF]

(4) Suppose that α < 0 , ( r , β ) ∈ _{A 3} . In this case, we have the phase portrait of (8) shown in Figure 2(c). Notice that H ( 0,0 ) = 0 . This case is completely similar to the case of α < 0 , ( r , β ) ∈ _{L 3 +} .

Theorem 5.

When the parameter groups α , r , β of system (8) satisfy the condition α < 0 with r > 0 , β = - 7 ^{r 2} or α < 0 , ( r , β ) ∈ _{A 3} , corresponding to the stable and unstable manifolds in the right phase plane of the equilibrium point 0 ( 0,0 ) in Figure 2(m) or Figure 2(c) defined by H ( ψ , y ) = 0 ; (1) has a cuspon solution given by (23).

Remark 6.

To the best of our knowledge, solutions (23) obtained for (1) have not been reported in the literature.

3.3. Compacton

(5) Suppose that α < 0 , ( r , β ) ∈ _{L 3 -} . In this case, we have the phase portrait of (8) shown in Figure 2(n). Notice that H ( _{ψ -} , 0 ) = H ( r , 0 ) = - ( 11 / 30 ) ^{r 5} . We see from (10) that arch curve connecting _{A +} in the right side of the straight line ψ = r has the algebraic equation [figure omitted; refer to PDF] where _{ψ 1,2} = ( ( - 11 ± 5 11 / 3 ) r ) / 4 . Thus, by using the first equation of (8) and (25), we obtain the parametric representation of the arch curve as follows: [figure omitted; refer to PDF]

Therefore, we have the following.

Theorem 7.

(1) When the parameter groups α , r , β of system (8) satisfy the condition α < 0 with r < 0 , β = - 7 ^{r 2} , there exists a homoclinic loop of system (12) given by branch of the curve H ( ψ , y ) = - ( 11 / 30 ) ^{r 5} .

(2) As the limit curves of a family of periodic orbits of system (8), the periodic curve in Figure 2(n) gives rise to a compacton soliton solution to (1), which has the exact parametric representation given by (26).

The profile of compacton soliton solution is shown in Figure 4(b).

Remark 8.

To the best of our knowledge, solutions (26) obtained for (1) have not been reported in the literature.

3.4. Loop Solution

(6) Suppose that α < 0 , ( r , β ) ∈ _{L 3 -} . We notice that the curves defined by H ( ψ , y ) = _{h +} correspond to different orbits of (8) consisting of two stable manifolds, two unstable manifolds of the saddle point _{A +} ( _{ψ +} , 0 ) , and the open curve passing through the point ( ( 7 / 2 ) r , 0 ) , respectively (see Figure 2(n)).

We next discuss the parametric representation of ψ ( ξ ) for these curves. We see from (10) that the arch curves in the right side of the straight line ψ = r have the algebraic equation [figure omitted; refer to PDF]

We first consider the unstable (or stable) manifold of the saddle point _{A +} ( _{ψ +} , 0 ) . On the basis of Lemma 2, we can take initial value ψ ( 0 ) = - ( 103 r / 50 ) , y ( 0 ) = - ( 97 r / 573750 ) 276735795 α r (or y ( 0 ) = ( 97 r / 573750 ) 276735795 α r ); we have from (27) that [figure omitted; refer to PDF] where Ψ ( ψ ) = ( ψ - r ) / ( ( ψ + 4 r ) ( ψ - ( 7 / 2 ) r ) ( ^{ψ 2} - 2 r ψ + ( 28 / 3 ) ^{r 2} ) ) . By introducing a new variable χ , from (28), one obtains the parametric representation of the unstable manifold: [figure omitted; refer to PDF] where [figure omitted; refer to PDF]

In addition, taking the initial condition ψ ( 0 ) = 7 r / 2 , y ( 0 ) = 0 , for the left open curve defined by a branch of H ( ψ , y ) = _{h +} , we have the parametric representation: [figure omitted; refer to PDF]

Employing the above formulas to draw the graphs of ψ ( ξ ) , we obtain corresponding wave profiles shown in Figure 5.

The profiles of waves when α < 0 , ( r , β ) ∈ _{L 3 -} . (a) shows the kink wave corresponding to the unstable manifold. (b) shows the kink wave corresponding to the stable manifold. (c) shows the breaking wave corresponding to the open curve passing through the point _{A +} ( _{ψ +} , 0 ) . (d) shows the wave profile of three curves--loop solution.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

Remark 9.

The loop solution, that is, the so-called loop soliton solution, is not one real soliton solution (see [28, 29]).

(7) Suppose that α < 0 , ( r , β ) ∈ _{A 4} . We notice that the curves defined by H ( ψ , y ) = 0 correspond to different orbits of (8) consisting of two stable manifolds, two unstable manifolds of the saddle point O ( 0,0 ) , and the open curve passing through the point ( _{ψ 2} , 0 ) , respectively (see Figure 2(d)).

We next discuss the parametric representation of ψ ( ξ ) for these curves. We see from (10) that the arch curves in the left side of the straight line ψ = r have the algebraic equation [figure omitted; refer to PDF] where _{ψ 1} < 0 < _{ψ 2} < _{ψ 3} .

We first consider the unstable (or stable) manifold of the saddle point O ( 0,0 ) . On the basis of Lemma 2, we can take initial value [figure omitted; refer to PDF]

We have from (32) that [figure omitted; refer to PDF] By introducing a new variable χ , from (34), one obtains the parametric representation of the unstable manifold: [figure omitted; refer to PDF] where [figure omitted; refer to PDF]

In addition, taking the initial condition ψ ( 0 ) = _{ψ 2} , y ( 0 ) = 0 , for the right open curve defined by a branch of H ( ψ , y ) = 0 , we have the parametric representation: [figure omitted; refer to PDF]

Employing the above formulas to draw the graphs of ψ ( ξ ) , we obtain corresponding wave profiles shown in Figure 6.

The profiles of waves when α < 0 , ( r , β ) ∈ _{A 4} . (a) shows the kink wave corresponding to the unstable manifold. (b) shows the kink wave corresponding to the stable manifold. (c) shows the breaking wave corresponding to the open curve passing through the point O ( 0,0 ) . (d) shows the wave profile of three curves--loop solution.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

Finally, we note that in Figures 2(e), 2(f), 2(g), 2(h), and 2(l), there is at least one loop solution to (1). Due to space limitations, we omit them.

Remark 10.

To the best of our knowledge, solutions (31) and (37) obtained for (1) have not been reported in the literature.

4. Discussion

In this paper, we used the qualitative analysis methods of a dynamical system to investigate the peakon, cuspon, compacton, and loop solutions of 3DKP ( 3,2 ) equation. Our procedure shows that the 3DKP ( 3,2 ) equation either has peakon, cuspon, compacton, or loop solutions. The phase portrait bifurcation of the traveling wave system corresponding to the equation is given. Particularly, cuspon and loop belong to the compound-type solutions; that is, they consist of two or three branches of nonsmooth solutions due to the existence of the singular straight line ψ = r in the corresponding phase portraits. The approach we used is simple and can be extended to study the soliton solutions of some other equations.

Acknowledgments

The authors wish to thank the anonymous reviewers for their helpful comments and suggestions. This work is supported by the National Natural Science Foundation of China (nos. 11361017 and 11161013) and Natural Science Foundation of Guangxi (nos. 2012GXNSFAA053003 and 2013GXNSFAA019010).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] M. J. Ablowitz, P. A. Clarkson Solitons, Nonlinear Evolution Equations and Inverse Scattering , vol. 149, of London Mathematical Society Lecture Note Series, Cambridge University Press, Cambridge, UK, 1991.

[2] V. B. Matveev, M. A. Salle Darboux Transformation and Solitons , Springer, Berlin, Germany, 1991.

[3] C. H. Gu, H. S. Hu, Z. X. Zhou Darboux Transformations in Soliton Theory and Its Geometric Applications , Shanghai Science and Technology, Shanghai, China, 1999.

[4] R. Hirota, J. Satsuma, "Soliton solutions of a coupled Korteweg-de Vries equation," Physics Letters A , vol. 85, no. 8-9, pp. 407-408, 1981.

[5] P. J. Olver Applications of Lie Groups to Differential Equations , Springer, New York, NY, USA, 1993.

[6] G. W. Bluman, S. Kumei Symmetries and Differential Equations , Springer, Berlin, Germany, 1989.

[7] J. Li, Z. Liu, "Traveling wave solutions for a class of nonlinear dispersive equations," Chinese Annals of Mathematics B , vol. 23, no. 3, pp. 397-418, 2002.

[8] M. Z. Wei, S. Q. Tang, H. L. Fu, G. X. Chen, "Single peak solitary wave solutions for the generalized KP-MEW ( 2 , 2 ) equation under boundary condition," Applied Mathematics and Computation , vol. 219, no. 17, pp. 8979-8990, 2013.

[9] A. Wazwaz, "A reliable treatment of the physical structure for the nonlinear equation K ( m , n ) ," Applied Mathematics and Computation , vol. 163, no. 3, pp. 1081-1095, 2005.

[10] A.-M. Wazwaz, "Analytic study on nonlinear variants of the RLW and the PHI-four equations," Communications in Nonlinear Science and Numerical Simulation , vol. 12, no. 3, pp. 314-327, 2007.

[11] Z. Yan, H. Zhang, "New explicit and exact travelling wave solutions for a system of variant Boussinesq equations in mathematical physics," Physics Letters A , vol. 252, no. 6, pp. 291-296, 1999.

[12] Z. Yan, "New explicit travelling wave solutions for two new integrable coupled nonlinear evolution equations," Physics Letters A , vol. 292, no. 1-2, pp. 100-106, 2001.

[13] E. Fan, "Uniformly constructing a series of explicit exact solutions to nonlinear equations in mathematical physics," Chaos, Solitons and Fractals , vol. 16, no. 5, pp. 819-839, 2003.

[14] S. Zhang, "New exact solutions of the KdV-Burgers-Kuramoto equation," Physics Letters A , vol. 358, no. 5-6, pp. 414-420, 2006.

[15] M.-L. Wang, "Exact solutions for a compound KdV-Burgers equation," Physics Letters A , vol. 213, no. 5-6, pp. 279-287, 1996.

[16] F. Xie, Z.-Y. Yan, "Compactons and noncompactons to three-dimensional Kadomtsev-Petviashvili equation with nonlinear dispersion," Chaos, Solitons & Fractals , vol. 36, no. 2, pp. 278-282, 2008.

[17] E. A. Kuznetsov, S. K. Turitsyn, "Two-and three-dimensional solitons in weakly dispersive media," Zhurnal Eksperimental'noi i Teoreticheskoi Fiziki , vol. 82, pp. 1457-1463, 1982.

[18] E. A. Kuznietsov, S. L. Musher, "Effect of collapse of sound waves on the structure of collisionless shock waves in a magnetized plasma," Journal of Experimental and Theoretical Physics , vol. 64, no. 5, pp. 947-955, 1986.

[19] P. Rosenau, J. M. Hyman, "Compactons: solitons with finite wavelength," Physical Review Letters , vol. 70, no. 5, pp. 564-567, 1993.

[20] M. Inc, "Compact and noncompact structures of a three-dimensional 3DKP(m ,n ) equation with nonlinear dispersion," Applied Mathematics Letters , vol. 26, no. 4, pp. 437-444, 2013.

[21] S. N. Chow, J. K. Hale Method of Bifurcation Theory , Springer, New York, NY, USA, 1981.

[22] J. Li, Z. Liu, "Smooth and non-smooth traveling waves in a nonlinearly dispersive equation," Applied Mathematical Modelling , vol. 25, no. 1, pp. 41-56, 2000.

[23] Y. Xie, H. Fu, S. Tang, "Peaked and smooth solitons for K * ( 4 , 1 ) equation," Journal of Applied Mathematics , vol. 2013, 2013.

[24] H. Fu, Y. Tang, S. Tang, H. Yan, Q. Liu, "Peaked and smooth solitons for K *(3, 1) equation, Indian," Journal of Physics , vol. 88, pp. 83-91, 2014.

[25] J. Li, Z. Qiao, "Peakon, pseudo-peakon, and cuspon solutions for two generalized Camassa-Holm equations," Journal of Mathematical Physics , vol. 54, no. 12, 2013.

[26] R. Camassa, D. D. Holm, "An integrable shallow water equation with peaked solitons," Physical Review Letters , vol. 71, no. 11, pp. 1661-1664, 1993.

[27] P. F. Byrd, M. D. Friedman Handbook of Elliptic Integrals for Engineers and Scientists , Springer, Berlin, Germany, 1971.

[28] J. B. Li Singular Nonlinear Travelling Wave Equations: Bifurcations and Exact Solutions , Science Press, Beijing, China, 2013.

[29] J. Li, "Dynamical understanding of loop soliton solution for several nonlinear wave equations," Science in China. Series A. Mathematics , vol. 50, no. 6, pp. 773-785, 2007.