Abstract

An integrable 2-component Camassa-Holm (2-CH) shallow water system is studied by using integral bifurcation method together with a translation-dilation transformation. Many traveling wave solutions of nonsingular type and singular type, such as solitary wave solutions, kink wave solutions, loop soliton solutions, compacton solutions, smooth periodic wave solutions, periodic kink wave solution, singular wave solution, and singular periodic wave solution are obtained. Further more, their dynamic behaviors are investigated. It is found that the waveforms of some traveling wave solutions vary with the changes of parameter, that is to say, the dynamic behavior of these waves partly depends on the relation of the amplitude of wave and the level of water.

1. Introduction

In this paper, employing the integral bifurcation method together with a translation-dilation transformation, we will study an integrable 2-component Camassa-Holm (2-CH) shallow water system [1] as follows: which is a nonlinear dispersive wave equation that models the propagation of unidirectional irrotational shallow water waves over a flat bed [2], as well as water waves moving over an underlying shear flow [3]. Equation (1.1) also arises in the study of a certain non-Newtonian fluids [4] and also models finite length, small amplitude radial deformation waves in cylindrical hyperelastic rods [5], where and are two dimensionless parameters. The interpretation of , respectively, describes the horizontal fluid velocity and the density in the shallow water regime, where the variable describes the horizontal velocity of the fluid in direction at time , and the variable is related to the free surface elevation from equilibrium position (or scalar density) with the boundary assumptions. The parameter denotes the level of water, the parameter denotes the typical amplitude of the water wave, and the parameter denotes the typical wavelength of the water wave. The constant denotes the speed of the water current which is related to the shallow water wave speed. The case , respectively, corresponds to the two situations in which the gravity acceleration points downwards and upwards. Especially, when the speed of the water current and the parameter , (1.1) becomes the following form where . The system (1.2) appeared in [1], which was first derived by Constantin and Ivanov from the Green-Naghdi equations [6, 7] via the hydrodynamical point of view. Under the scaling , (1.1) can be reduced to the following two-component generalization of the well-known 2-component Camassa-Holm (2-CH) system [1, 810]: where . The equation (1.3) attracts much interest since it appears. Attention was more paid on the local well-posedness, blow-up phenomenon, global existence, and so forth. When , (1.3) has been studied by many authors, see [1126] and references cited therein. Especially, under the parametric conditions , (1.3) can be reduced to the celebrated Camassa-Holm equation [27, 28], In 2006, the integrability of (1.3) for was proved and some peakon and multikink solutions of this system were presented by Chen et al. in [9]. In 2008, the Lax pair of (1.3) for any value of was given by Constantin and Ivanov in [1]. In [29], by using the method of dynamical systems, under the traveling wave transformation , some explicit parametric representations of exact traveling wave solutions of (1.3) were obtained. But the loop solitons were not obtained in [29]. In [30], under , one-loop and two-loop soliton solutions and multisoliton solutions of (1.3) are obtained by using the Darboux transformations.

Although (1.1) can be reduced to (1.3) by the scaling transformation , the dynamic properties of some traveling wave solutions for these two equations are very different. In fact, the dynamic behaviors of some traveling waves of system (1.1) partly depends on the relation () of the amplitude of wave and the level (deepness) of water. In the other words, their dynamic behavior vary with the changes of parameter , that is, the changes of ratio for the amplitude of wave and the deepness of water . In addition, compared with the research results of (1.3), the research results for (1.1) are few in the existing literatures. Thus, (1.1) is very necessary to be further studied.

It is worthy to mention that the solutions obtained by us in this paper are different from others in existing references, such as [9, 1217, 29]. On the other hand, under different transformations, by using different methods, different results will be presented. In this paper, by using the integral bifurcation method [31], we will investigate different kinds of new traveling wave solutions of (1.1) and their dynamic properties under the translation-dilation transformation . By the way, the integral bifurcation method possessed some advantages of the bifurcation theory of the planar dynamic system [32] and auxiliary equation method (see [33, 34] and references cited therein), it is easily combined with computer method [35] and useful for many nonlinear partial diffential equations (PDEs) including some PDEs with high power terms, such as equation [36]. So, by using this method, we will obtain some new traveling wave solutions of (1.1). Some interesting phenomena will be presented.

The rest of this paper is organized as follows. In Section 2, we will derive the two-dimensional planar system of (1.1) and its first integral equations. In Section 3, by using the integral bifurcation method, we will obtain some new traveling wave solutions of nonsingular type and singular type and investigate their dynamic behaviors.

2. The Two-Dimensional Planar System of (1.1) and Its First Integral Equations

Obviously, (1.1) can be rewritten as the following form In order to change the PDE (2.1) into an ordinary differential equation, we make a transformation where and is an arbitrary nonzero constant. In fact, (2.2) is a translation-dilation transformation, it has been used extensively in many literatures. Its idea came from many existing references. For example, in [37], the expression is a dilation transformation. In [38], the expression is a translation transformation. In [39, 40], the expression which was first used by Parkes and Vakhnenko is also a translation-dilation transformation.

After substituting (2.2) into (2.1), integrating them once yields where the integral constant and denotes . From the second equation of (2.3), we easily obtain where . Substituting (2.4) into the first equation of (2.3), we obtain where denotes .

Let . Thus (2.5) can be reduced to 2-dimensional planar system as follows: where denotes .

Making a transformation (2.6) becomes where is a parameter. From the point of view of the geometric theory, the parameter is a fast variable, but the parameter is a slow one. The system (2.8) is still a singular system through the system (2.6) becomes (2.8) by the transformation (2.7). This case is not same as that in [29]. Of course, as in [29], we can also change the system (2.6) into a regular system under another transformation . But we do not want to do that due to the following two reasons: (i) we need to keep the singularity of the original system; (ii) we need not to make any analysis of phase portraits as in [29].

Obviously, systems (2.6) and (2.8) have the same first integral as follows: where is an integral constant.

3. Traveling Wave Solutions of Nonsingular Type and Singular Type for (1.1) and Their Dynamic Behaviors

In this section, we will investigate different kinds of exact traveling wave solutions for (1.1) and their dynamic behaviors.

It is easy to know that (2.9) can be reduced to four kinds of simple equations when the parametric conditions satisfy the following four cases:

Case 1. and .

Case 2. and .

Case 3. and .

Case 4. .

We mainly aim to consider the new results of (1.1), so we only discuss the first two typical cases in this section. The other two cases can be similarly discussed, here we omit them.

3.1. The Exact Traveling Wave Solutions under Case 1

Under the parametric conditions and , (2.9) can be reduced to where . Substituting (3.1) into the first expression in (2.8) yields Write The if only if . By using the exact solutions of (3.2), we can obtain different kinds of exact traveling wave solutions of parametric type of (1.1), see the following discussion.(i)If , then (3.2) has one exact solution as follows: Substituting (3.4) into (2.7), and then integrating it yields Substituting (3.4) into (2.2) and (2.4), then combining with (3.5), we obtain a couple of soliton-like solutions of (1.1) as follows: (ii)If , then (3.2) has one exact solution as follows: Similarly, by using (3.8), (2.7), (2.2), and (2.4), we obtain two couples of soliton-like solutions of (1.1) as follows: (iii)If and , then (1.1) has one couple of soliton-like solutions as follows: (iv)If and (i.e., ), then (1.1) has one couple of soliton-like solutions as follows: (v)If and (i.e. ), then (1.1) has one couple of kink and antikink wave solutions as follows: (vi)If and (i.e. ), then (1.1) has one couple of soliton-like solutions as follows: In order to show the dynamic properties of above soliton-like solutions and kink and antikink wave solutions intuitively, as examples, we plot their graphs of some solutions, see Figures 1, 2, 3, and 4. Figures 1(a)1(h) show the profiles of multiwaveform to solution the first solution of (3.6) for fixed parameters and different -values. Figures 2(a)2(d) show the profiles of multiwaveform to solution (3.7) for fixed parameters and different -values. Figures 3(a)3(d) show the profiles of multiwaveform to solution (3.9) for fixed parameters and different -values. Figures 4(a)-4(b) show the profiles of kink wave solution (the first formula of (3.12)) for fixed parameters and different -values.

We observe that some profiles of above soliton-like solutions are very much sensitive to one of parameters, that is, their profiles are transformable (see Figures 13). But the others are not, their waveforms do not vary no matter how the parameters vary (see Figure 4). Some phenomena are very similar to those in [41, 42]. In [41], the waveforms of soliton-like solution of the generalized KdV equation vary with the changes of parameter and depend on the velocity extremely. Similarly, in [42], the properties of some traveling wave solutions of the generalized KdV-Burges equation depend on the dissipation coefficient ; if dissipation coefficient , it appears as a monotonically kink profile solitary wave; if , it appears as a damped oscillatory wave.

From Figures 1(a)1(h), it is easy to know that the profiles of solution (3.6) vary gradually, its properties depend on the parameter . When parametric values of increase from to , the solution (3.6) has eight kinds of waveforms: Figure 1(a) shows a shape of antikink wave when ; Figure 1(b) shows a shape of transmutative antikink wave when ; Figure 1(c) shows a shape of thin and dark solitary wave when ; Figure 1(d) shows a shape of fat and dark solitary wave when ; Figure 1(e) shows a shape of compacton wave when ; Figure 1(f) shows a shape of fat loop soliton when ; the Figure 1(g) shows a shape of thin loop soliton when ; Figure 1(h) shows a shape of oblique loop soliton when .

Similarly, from Figures 2(a)2(d), it is also easy to know that the profiles of solution (3.7) are transformable, but their changes are not gradual, it depends on the parameter extremely. When parametric values of increase from to , the solution (3.7) has four kinds of waveforms: Figure 2(a) shows a shape of smooth kink wave when ; Figure 2(b) shows a shape of fat and bright solitary wave when ; Figure 2(c) shows a shape of thin and bright solitary wave when ; Figure 2(d) shows a shape of singular wave of cracked loop soliton when . Especially, the changes of waveforms from Figure 2(a) to Figure 2(b) and from Figure 2(c) to Figure 2(d) happened abruptly.

As in Figure 1, the profiles of the first solution of (3.9) in Figure 3 vary gradually. When parametric values of increase from to , the waveform becomes a bright compacton from a shape of bright solitary wave, then becomes a shape of fat loop soliton and a shape of thin loop soliton at last.

Different from the properties of solutions (3.6), (3.7) and (3.9), the property of the first solution of (3.12) is stable. The profile of the first solution of (3.12) is not transformable no matter how the parameters vary. Both Figures 4(a) and 4(b) show a shape of smooth kink wave.

3.2. The Traveling Wave Solutions Under Case 2

Under the parametric conditions and , (2.9) becomes where . Substituting (3.14) into the first expression in (2.8) yields We know that the case corresponds to the situation in which the gravity acceleration points upwards. As an example, in this subsection, we only discuss the case . The case can be similarly discussed, but we omit them here. Especially, when , the above values of can be reduced to .

(i) If then the . Under these parametric conditions, (3.15) has a Jacobi elliptic function solution as follows: or where is the model of Jacobi elliptic function and .

Substituting (3.16) and (3.17) into (2.7), respectively, we obtain Substituting (3.16) into (2.2) and (2.4), using (3.18) and the transformation , we obtain a couple of periodic solutions of (1.1) as follows: Substituting (3.17) into (2.2) and (2.4), using (3.19) and the transformation , we also obtain a couple of periodic solutions of (1.1) as follows:

(ii)  If , then the . In the case of these parametric conditions, (3.15) has a Jacobi elliptic function solution as follows:

As in the first case (i) using the same method, we obtain a couple of periodic solutions of (1.1) as follows: or Especially, when (i.e., ), . From (3.23) and (3.24), we obtain a couple of kink-like solutions of (1.1) as follows:

(iii)  If , then the . In the case of these parametric conditions, (3.15) has a Jacobi elliptic function solution as follows: As in the first case (i), similarly, we obtain a couple of periodic solutions of (1.1) as follows: where .

(iv)  If , then the . In the case of these parametric conditions, (3.15) has a Jacobi elliptic function solution as follows: Similarly, we obtain a couple of periodic solutions of (1.1) as follows: where .

(v)  If , then the . In the case of these parametric conditions, (3.15) has a Jacobi elliptic function solution as follows:

As in the first case (i) using the same method, we obtain a couple of periodic solutions of (1.1) as follows:

In order to show the dynamic properties of above periodic solutions intuitively, as examples, we plot their graphs of the solutions (3.20), (3.23) and (3.24), see Figures 5 and 6.

Figures 5(a) and 5(b) show two shapes of smooth and continuous periodic waves, all of them are nonsingular type. Figure 6(a) shows a shape of periodic kink wave. Figure 6(b) shows a shape of singular periodic wave. Both Figures 6(a) and 6(b) show two discontinuous periodic waves, all of them are singular type.

From the above illustrations, we find that the waveforms of some solutions partly depend on wave parameters. Indeed, in 2006, Vakhnenko and Parkes’s work [43] successfully explained similar phenomena. In [43], the graphical interpretation of the solution for gDPE is presented. In this analysis, the 3D-spiral (whether one loop from a spiral or a half loop of a spiral) has the different projections that is the essence of the possible solutions. Of caurse, this approach also can be employed to the 2-component Camassa-Holm shallow water system (1.1). By using the Vakhnenko and Parkes’s theory, the phenomena which appeared in this work are easily understood, so we omit this analysis at here. However, it is necessary to say something about the cracked loop soliton (Figure 2(d)) and the singular periodic wave (Figure 6(b)), what are the 3D-curves with projections associated with two peculiar solutions (3.7) and (3.24) which are shown in Figure 2(d) and Figure 6(b)? In order to answer this question, by using the Vakhnenko and Parkes’s approach, as an example we give the 3D-curves of solution (3.24) and their projection curve in the -plane, which are shown in Figure 7. For convenience to distinguish them, we colour the 3D-curves red and colour their projection curve green. From Figure 7, we can see that the 3D-curves are not intersected, but their projection curve is intersected in -plane.

4. Conclusion

In this work, by using the integral bifurcation method together with a translation-dilation transformation, we have obtained some new traveling wave solutions of nonsingular type and singular type of 2-component Camassa-Holm equation. These new solutions include soliton solutions, kink wave solutions, loop soliton solutions, compacton solutions, smooth periodic wave solutions, periodic kink wave solution, singular wave solution, and singular periodic wave solution. By investigating these new exact solutions of parametric type, we found some new traveling wave phenomena, that is, the waveforms of some solutions partly depend on wave parameters. For example, the waveforms of solutions (3.6), (3.7), and (3.9) vary with the changes of parameter. These are three peculiar solutions. The solution (3.6) has five kinds of waveforms, which contain antikink wave, transmutative antikink wave, dark soliton, compacton, and loop soliton according as the parameter varies. The solution (3.7) has three kinds of waveforms, which contain kink wave, bright soliton, and singular wave (cracked loop soliton) according as the parameter varies. The solution (3.9) also has three kinds of waveforms, which contain bright soliton, compacton, and loop soliton according as the parameter varies. These phenomena show that the dynamic behavior of these waves partly depends on the relation of the amplitude of wave and the level of water.

Acknowledgments

The authors thank reviewers very much for their useful comments and helpful suggestions. This work was financially supported by the Natural Science Foundation of China (Grant no. 11161020) and was also supported by the Natural Science Foundation of Yunnan Province (no. 2011FZ193).