# Exact Solutions of a High-Order Nonlinear Wave Equation of Korteweg-de Vries Type under Newly Solvable Conditions

Authors: Weiguo Rui
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

By using the integral bifurcation method together with factoring technique, we study a water wave model, a high-order nonlinear wave equation of KdV type under some newly solvable conditions. Based on our previous research works, some exact traveling wave solutions such as broken-soliton solutions, periodic wave solutions of blow-up type, smooth solitary wave solutions, and nonsmooth peakon solutions within more extensive parameter ranges are obtained. In particular, a series of smooth solitary wave solutions and nonsmooth peakon solutions are obtained. In order to show the properties of these exact solutions visually, we plot the graphs of some representative traveling wave solutions.


---

## Body

## 1. Introduction

In this work, we will study the following high-order nonlinear wave equation of Korteweg-de Vries type:(1)
η
t
+
η
x
+
α
η
η
x
+
β
η
x
x
x
+
ρ
1
α
2
η
2
η
x
+
α
β
(
ρ
2
η
η
x
x
x
+
ρ
3
η
x
η
x
x
)
+
ρ
4
α
3
η
3
η
x
+
α
2
β
(
ρ
5
η
2
η
x
x
x
+
ρ
6
η
η
x
η
x
x
+
ρ
7
η
x
3
)
=
0
.
This is an important model of water wave derived by Fokas [1] in 1995, where α
=
3
A
/
2, β
=
B
/
6, and 0
<
α, β
≪
1.Obviously, (1) is a very complex partial deferential equation, it has nine parameters α, β, ρ
i, (i
=
1,2
,
…
,
7), and contains both high-order derivative terms and multinonlinear terms. It is very different from the original KdV equation. Regarding the ρ
1, ρ
2, ρ
3, ρ
4, ρ
5, ρ
6, ρ
7 as free parameters, Tzirtzilakis et al. [2] investigated solitary wave solutions of (1) and they called (1) high-order wave equation of Korteweg-de Vries type. Just as Tzirtzilakis et al. [2] said that investigations of solitary wave solutions of (1) are more physically and practically meaningful. The motion described by the model (1) is a 2-dimensional, inviscid, and incompressible fluid (water) lying above a horizontal flat bottom located at y
¯
=
-
ħ
0 (ħ
0 is a constant) and letting the air above the water. It turns out that, for such a system if the vorticity is zero initially, it remains zero. The fluids (waters) analyzed by Fokas are only irrotational flows. This system is characterized by two parameters α
=
3
A
/
2 and β
=
B
/
6 with A
=
a
/
ħ
0 and B
=
ħ
0
2
/
ℓ
2, where a and ℓ are two typical values of the amplitude and of the wavelength of the waves. The parameters a, ħ
0, and ℓ satisfy the condition a
≪
ħ
0
<
ℓ because the system is a model of short amplitude and long wavelength.Whenρ
1
=
ρ
2
=
ρ
3
=
ρ
4
=
ρ
5
=
ρ
6
=
ρ
7
=
0, (1) can be reduced to the classical (original) KdV equation:
(2)
η
t
+
η
x
+
α
η
η
x
+
β
η
x
x
x
=
0
.
In [1], Fokas assumed that O
(
β
)
<
O
(
α
) and β
≈
α
2. According to this assumption, we easily know that O
(
α
2
β
)
<
O
(
α
3
) and O
(
α
2
β
)
<
O
(
α
β
). Neglecting two high-order infinitesimal terms of O
(
α
3
,
α
2
β
), (1) can be reduced to another high-order wave equation of KdV type [2–5] as follows:
(3)
η
t
+
η
x
+
α
η
η
x
+
β
η
x
x
x
+
ρ
1
α
2
η
2
η
x
+
α
β
(
ρ
2
η
η
x
x
x
+
ρ
3
η
x
η
x
x
)
=
0
.
Equation (3) can be regarded as a special case of (1) for ρ
4
=
ρ
5
=
ρ
6
=
ρ
7
=
0. In [1], it was observed that (3) can be reduced by the local transformation of coordinates
(4)
η
=
v
-
α
ρ
1
v
2
-
β
(
3
ρ
1
+
7
4
ρ
2
-
1
2
ρ
3
)
v
x
x
to a completely integrable PDE as follows:
(5)
v
t
-
3
2
β
ρ
2
v
x
x
t
+
β
(
1
-
3
2
ρ
2
)
v
x
x
x
+
α
v
v
x
-
1
2
α
β
ρ
2
(
v
v
x
x
x
+
2
v
x
v
x
x
)
.
Equation (5) was first derived in [6] by using the method of bi-Hamiltonian systems and its Lax pair was given in [7].Neglecting the highest-order infinitesimal term ofO
(
α
2
β
), (1) can be reduced to a new generalized KdV equation as follows:
(6)
η
t
+
η
x
+
α
η
η
x
+
β
η
x
x
x
+
ρ
1
α
2
η
2
η
x
+
α
β
(
ρ
2
η
η
x
x
x
+
ρ
3
η
x
η
x
x
)
+
ρ
4
α
3
η
3
η
x
=
0
.
We call it a generalized KdV equation of neglecting the highest-order infinitesimal term [8]. In fact, (6) can be regarded as another special case of (1) for ρ
5
=
ρ
6
=
ρ
7
=
0; it is also third-order approximate equation of higher-order KdV type.From the above references and the references cited therein, we know that (1) is a very important model of water wave. However, (1) is too complex to obtain its exact solution under universal conditions. Only under some very special parametric conditions, its exact solutions were obtained in existing literatures [2–5, 8–12]. In addition, in [13], under different kinds of parametric conditions, Marinakis discussed two integrable cases for the third-order approximation model (1). In [14], Marinakis proved that (1) and its some special cases are integrable. Generally, a system is regarded as that it is integrable if it has Darboux transformation, Lax pair, bilinear structure and multilinear structure, Hamilton function, first integral function (equation), symmetrical structure (i.e., symmetry), conservation law, and so forth. In [15], Gandarias and Bruzon proved that (1) is self-adjoint if and only if ρ
3
=
2
ρ
2, and ρ
7
=
ρ
6
-
3
ρ
5.From the above research backgrounds of (1), we can see that its exact solutions under universal conditions are hard to obtain because it is highly nonlinear equation and most probably it is not integrable equation in general. Thus, large numbers of research results are still concentrated on the generalized KdV, mKdV equations [16, 17] and other some high-order equations with KdV type, such as KdV-Burgers equation [18, 19] and KdV-Burgers-Kuramoto equation [20]. Therefore, studying solvability and finding exact solutions of (1) within more extensive parameter ranges are very important and necessary. In this paper, based on the works in [8], by using the integral bifurcation method [21], we will investigate solvable conditions and exact traveling wave solutions of (1) within more extensive parameter ranges. It is different from those special cases which were considered by authors in existing literatures; we will discuss a newly solvable case ρ
3
=
(
m
+
1
)
ρ
2, ρ
6
=
2
(
m
+
1
)
ρ
5, and ρ
7
=
m
ρ
5, m
∈
N for (1); indeed, (1) is solvable under these parametric conditions. Though this is still a special case for (1), it is more general than the cases on the assumptions of parametric conditions which appeared in existing literatures. In particular, the results which will be obtained in this paper are very new and different from those in [8, 9].The rest of this paper is organized as follows. In Section2, we will derive two-dimensional dynamical system which is equivalent to (1) and give its first integrals. In Section 3, by using the integral bifurcation method, we will investigate different kinds of exact traveling wave solutions of (1) within more extensive parameter ranges and discuss their dynamic properties.


## 2. The First Integrals of (1) under Newly Integrable Conditions

Making a transformationη
(
x
,
t
)
=
ϕ
(
ξ
) with ξ
=
x
-
c
t, (1) can be reduced to the following ODE:
(7)
(
1
-
c
)
ϕ
′
+
α
ϕ
ϕ
′
+
β
ϕ
′′′
+
ρ
1
α
2
ϕ
2
ϕ
′
+
α
β
(
ρ
2
ϕ
ϕ
′′′
+
ρ
3
ϕ
′
ϕ
′′
)
+
ρ
4
α
3
ϕ
3
ϕ
′
+
α
2
β
[
ρ
5
ϕ
2
ϕ
′′′
+
ρ
6
ϕ
ϕ
′
ϕ
′′
+
ρ
7
(
ϕ
′
)
3
]
=
0
,
where ϕ
′
=
d
ϕ
/
d
ξ, c is a wave velocity which moves along the direction of x-axis, and c
≠
0. It is easy to find that (7) can be integrated once under the conditions ρ
3
=
(
m
+
1
)
ρ
2, ρ
6
=
2
(
m
+
1
)
ρ
5, and ρ
7
=
m
ρ
5, where m
∈
N. These are newly solvable conditions and they are different from others in [8, 9, 13–15], and the ranges of parametric values are more extensive than others in existing references. Thus, under these solvable conditions, integrating (7) once and setting the integral constant as zero yield
(8)
(
1
-
c
)
ϕ
+
1
2
α
ϕ
2
+
β
ϕ
′′
+
1
3
ρ
1
α
2
ϕ
3
+
1
4
ρ
4
α
3
ϕ
4
+
ρ
2
α
β
[
ϕ
ϕ
′′
+
1
2
m
(
ϕ
′
)
2
]
+
ρ
5
α
2
β
[
ϕ
2
ϕ
′′
+
m
ϕ
(
ϕ
′
)
2
]
=
0
.When1
+
α
ρ
2
ϕ
+
α
2
ρ
5
ϕ
2
≠
0, (8) can be reduced to a singular two-dimensional system as follows:
(9)
d
ϕ
d
ξ
=
y
,
d
y
d
ξ
=
(
12
(
c
-
1
)
ϕ
-
6
α
ϕ
2
-
4
ρ
1
α
2
ϕ
3
-
3
ρ
4
α
3
ϕ
4
i
-
6
m
α
β
ρ
2
y
2
-
12
m
α
2
β
ρ
5
ϕ
y
2
)
×
(
12
β
(
1
+
α
ρ
2
ϕ
+
α
2
ρ
5
ϕ
2
)
)
-
1
,
where ϕ
′
=
d
ϕ
/
d
ξ
=
y. However, when 1
+
α
ρ
2
ϕ
+
α
2
ρ
5
ϕ
2
=
0, (9) is not equivalent to (8). In order to obtain a system which is equivalent to (8), we make the following scalar transformation
(10)
d
ξ
=
12
β
(
1
+
α
ρ
2
ϕ
+
α
2
ρ
5
ϕ
2
)
d
τ
.
Equation (9) can be changed into a regular two-dimensional system as follows:
(11)
d
ϕ
d
τ
=
12
β
(
1
+
α
ρ
2
ϕ
+
α
2
ρ
5
ϕ
2
)
y
,
d
y
d
τ
=
12
(
c
-
1
)
ϕ
-
6
α
ϕ
2
-
4
ρ
1
α
2
ϕ
3
-
3
ρ
4
α
3
ϕ
4
-
6
m
α
β
ρ
2
y
2
-
12
m
α
2
β
ρ
5
ϕ
y
2
,
where τ is a parameter. Thus, (11) is equivalent to (8). Obviously, (11) and (9) have the same one first integral as follows:
(12)
(
1
+
α
ρ
2
ϕ
+
α
2
ρ
5
ϕ
2
)
y
2
=
1
6
β
∫
[
12
(
c
-
1
)
ϕ
+
6
α
ϕ
2
-
4
ρ
1
α
2
ϕ
3
hihhhhhh
-
3
ρ
4
α
3
ϕ
4
]
(
1
+
α
ρ
2
ϕ
+
α
2
ρ
5
ϕ
2
)
m
-
1
d
ϕ
.
It is easy to find that the integral of right side of (12) is always integrable once the m is given by an idiographic integer. In fact, the first integrals obtained by (12) have certain determinate orderliness; see the following discussions.(i)
Whenm
=
1, (12) can be reduced to
(13)
y
2
=
(
ϕ
2
[
(
c
-
1
)
-
(
1
3
)
α
ϕ
-
(
1
6
)
ρ
1
α
2
ϕ
2
i
22222
-
(
1
10
)
ρ
4
α
3
ϕ
3
]
+
h
1
)
×
(
β
(
1
+
α
ρ
2
ϕ
+
α
2
ρ
5
ϕ
2
)
)
-
1
, where h
1 is an integral constant.(ii)
Whenm
=
2, (12) can be reduced to
(14)
y
2
=
ϕ
2
[
(
c
-
1
)
+
a
1
ϕ
+
a
2
ϕ
2
+
a
3
ϕ
3
+
a
4
ϕ
4
+
a
5
ϕ
5
]
+
h
2
β
(
1
+
α
ρ
2
ϕ
+
α
2
ρ
5
ϕ
2
)
2
, where h
2 is an integral constant and a
1
=
(
1
/
3
)
α
(
2
c
ρ
2
-
1
-
2
ρ
2
), a
2
=
(
1
/
12
)
α
2
(
-
2
ρ
1
-
6
ρ
5
+
6
c
ρ
5
-
3
ρ
2
), a
3
=
-
(
1
/
30
)
α
3
(
6
ρ
5
+
4
ρ
1
ρ
2
+
3
ρ
4
), a
4
=
-
(
1
/
36
)
α
4
(
4
ρ
1
ρ
5
+
3
ρ
4
ρ
2
), a
5
=
-
(
1
/
14
)
ρ
4
α
5
ρ
5.(iii)
Whenm
=
3, (12) can be reduced to
(15)
y
2
=
(
ϕ
2
[
(
c
-
1
)
+
b
1
ϕ
+
b
2
ϕ
2
+
b
3
ϕ
3
iii
22
ii
+
b
4
ϕ
4
+
b
5
ϕ
5
+
b
6
ϕ
6
+
b
7
ϕ
7
]
+
h
3
)
×
(
β
(
1
+
α
ρ
2
ϕ
+
α
2
ρ
5
ϕ
2
)
3
)
-
1
,where h
3 is an integral constant and b
1
=
(
1
/
3
)
α
(
-
1
-
4
ρ
2
+
4
c
ρ
2
), b
2
=
(
1
/
6
)
α
2
(
-
3
ρ
2
-
ρ
1
+
3
c
ρ
2
2
-
3
ρ
2
2
-
6
ρ
5
+
6
c
ρ
5
), b
3
=
(
1
/
30
)
α
3
(
-
24
ρ
2
ρ
5
+
24
c
ρ
2
ρ
5
-
3
ρ
4
-
12
ρ
5
-
6
ρ
2
2
-
8
ρ
1
ρ
2
), b
4
=
(
1
/
18
)
α
4
(
-
6
ρ
5
2
-
6
ρ
2
ρ
5
+
6
c
ρ
5
2
-
2
ρ
2
2
ρ
1
-
3
ρ
4
ρ
2
-
4
ρ
1
ρ
5
), b
5
=
-
(
1
/
42
)
α
5
(
8
ρ
1
ρ
2
ρ
5
+
6
ρ
5
2
+
6
ρ
4
ρ
5
+
3
ρ
4
ρ
2
2
), b
6
=
-
(
1
/
24
)
α
6
ρ
5
(
2
ρ
1
ρ
5
+
3
ρ
4
ρ
2
), b
7
=
-
(
1
/
18
)
ρ
4
α
7
ρ
5
7.(iv)
Whenm is an integer and m
≥
4, (12) can be reduced to the following form:
(16)
y
2
=
ϕ
2
[
(
c
-
1
)
+
γ
1
ϕ
+
γ
2
ϕ
2
+
⋯
+
γ
2
m
+
1
ϕ
2
m
+
1
]
+
h
m
β
(
1
+
α
ρ
2
ϕ
+
α
2
ρ
5
ϕ
2
)
m
, where h
m (m
∈
N and m
≥
4) are integral constants and γ
1
,
γ
2
,
…
,
γ
2
m
+
1 are certain expressions of the parameters ρ
1, ρ
2, ρ
4, ρ
5, α, and β; these expressions are always determined from (12) once the m is given by an idiographic integer; we omit them here because their expressions are tediously long.


## 3. Exact Traveling Wave Solutions of (1) and Their Dynamic Properties

In this section, we will investigate exact traveling wave solutions of (1) and discuss their dynamic properties under different kinds of parametric conditions in the greatly possible parameter regions; see the following discussions. For the convenience of discussion, we always consider the cases of all the integral constants as zero (i.e., h
i
=
0, i
=
1,2
,
…
,
m) in the next discussions.Case 1.
Under the parametric conditionsh
1
=
0, m
=
1, and ρ
5
=
(
1
/
4
)
ρ
2
2, (13) can be reduced to
(17)
y
=
±
ϕ
(
c
-
1
)
-
(
1
/
3
)
α
ϕ
-
(
1
/
6
)
ρ
1
α
2
ϕ
2
-
(
1
/
10
)
ρ
4
α
3
ϕ
3
β
(
1
+
(
1
/
2
)
α
ρ
2
ϕ
)
.(a) Whenρ
4
<
0, c
=
1, and ρ
5
=
(
1
/
4
)
ρ
2
2, substituting (17) into the first equation of (9) yields
(18)
d
ϕ
X
11
(
ϕ
)
+
1
2
α
ρ
2
d
ϕ
X
12
(
ϕ
)
=
±
α
-
α
ρ
4
10
β
d
ξ
,
where X
11
(
ϕ
)
=
(
ϕ
-
ϕ
1
)
(
ϕ
-
0
)
3
(
ϕ
-
ϕ
2
) and X
12
(
ϕ
)
=
(
ϕ
-
ϕ
1
)
(
ϕ
-
0
)
(
ϕ
-
ϕ
2
) with ϕ
1,2
=
(
5
ρ
1
±
25
ρ
1
2
-
120
ρ
4
)
/
6
(
-
ρ
4
)
α and ϕ
2
<
0
<
ϕ
1
<
ϕ. Taking ϕ
(
0
)
=
ϕ
1 as the initial constant and then integrating (18), we obtain an implicit solution of (1) as follows:
(19)
2
ϕ
1
-
ϕ
2
-
ϕ
1
ϕ
2
E
(
θ
1
,
k
1
)
+
2
+
α
ρ
2
ϕ
2
ϕ
2
ϕ
1
-
ϕ
2
F
(
θ
1
,
k
1
)
=
α
-
α
ρ
4
10
β
|
ξ
|
,
where F
(
θ
1
,
k
1
) is an incomplete elliptic integral of the first kind and E
(
θ
1
,
k
1
) is a Legendre’s incomplete elliptic integral of the second kind with θ
1
=
arcsin
⁡
(
η
-
ϕ
1
)
/
η, k
1
=
-
ϕ
2
/
(
ϕ
1
-
ϕ
2
), and η
(
x
,
t
)
=
ϕ
(
ξ
) with ξ
=
x
-
t.(b) Whenρ
4
>
0, ρ
1
>
2
6
ρ
4
/
5, c
=
1, and ρ
5
=
(
1
/
4
)
ρ
2
2, substituting (17) into the first equation of (9) yields
(20)
d
ϕ
X
21
(
ϕ
)
+
1
2
α
ρ
2
d
ϕ
X
22
(
ϕ
)
=
±
α
α
ρ
4
10
β
d
ξ
,
where X
21
(
ϕ
)
=
(
0
-
ϕ
)
3
(
ϕ
-
ϕ
~
2
)
(
ϕ
-
ϕ
~
1
) and X
22
(
ϕ
)
=
(
0
-
ϕ
)
(
ϕ
-
ϕ
~
2
)
(
ϕ
-
ϕ
~
1
) with ϕ
~
1,2
=
(
-
5
ρ
1
∓
5
ρ
1
2
-
24
ρ
4
/
5
)
/
6
ρ
4
α and ϕ
~
1
<
ϕ
~
2
<
ϕ
<
0. Taking ϕ
(
0
)
=
ϕ
~
2 as the initial constant and then integrating (20), we obtain an implicit solution of (1) as follows:
(21)
-
2
ϕ
~
2
+
α
ρ
2
ϕ
~
2
-
ϕ
~
1
F
(
θ
2
,
k
2
)
+
2
ϕ
~
2
-
ϕ
~
1
E
(
θ
2
,
k
2
)
-
2
ϕ
~
2
η
-
ϕ
~
2
η
(
ϕ
~
1
-
η
)
=
α
α
ρ
4
10
β
|
ξ
|
,
where θ
2
=
arcsin
ϕ
~
1
(
η
-
ϕ
~
2
)
/
ϕ
~
2
(
η
-
ϕ
~
1
), k
2
=
ϕ
~
2
/
ϕ
~
1, and η
(
x
,
t
)
=
ϕ
(
ξ
) with ξ
=
x
-
t.(c) Whenρ
4
>
0, ρ
1
<
-
2
6
ρ
4
/
5, c
=
1, and ρ
5
=
(
1
/
4
)
ρ
2
2, substituting (17) into the first equation of (9) yields
(22)
d
ϕ
X
31
(
ϕ
)
+
1
2
α
ρ
2
d
ϕ
X
32
(
ϕ
)
=
±
α
α
ρ
4
10
β
d
ξ
,
where X
31
(
ϕ
)
=
(
ϕ
~
2
-
ϕ
)
(
ϕ
-
ϕ
~
1
)
(
ϕ
-
0
)
3 and X
32
(
ϕ
)
=
(
ϕ
~
2
-
ϕ
)
(
ϕ
-
ϕ
~
1
)
(
ϕ
-
0
) with ϕ
~
1,2
=
(
-
5
ρ
1
∓
5
ρ
1
2
-
24
ρ
4
/
5
)
/
6
ρ
4
α and 0
<
ϕ
~
1
<
ϕ
<
ϕ
~
2. Taking ϕ
(
0
)
=
ϕ
~
1 as the initial constant and then integrating (22), we obtain an implicit solution of (1) as follows:
(23)
2
ϕ
~
1
ϕ
~
2
E
(
θ
2
,
k
2
)
+
α
ρ
2
ϕ
~
2
F
(
θ
2
,
k
2
)
=
α
α
ρ
4
10
β
|
ξ
|
,
where θ
2 and k
2 are given above and η
(
x
,
t
)
=
ϕ
(
ξ
) with ξ
=
x
-
t.(d) Whenρ
4
<
0, c
≠
1, and equation (
c
-
1
)
-
(
1
/
3
)
α
ϕ
-
(
1
/
6
)
ρ
1
α
2
ϕ
2
-
(
1
/
10
)
ρ
4
α
3
ϕ
3
=
0 has three real roots, substituting (17) into the first equation of (9) yields
(24)
d
ϕ
ϕ
X
4
(
ϕ
)
+
1
2
α
ρ
2
d
ϕ
X
4
(
ϕ
)
=
±
α
-
α
ρ
4
10
β
d
ξ
,
where X
4
(
ϕ
)
=
(
ϕ
-
a
1
)
(
ϕ
-
b
1
)
(
ϕ
-
c
1
) and the three real roots a
1, b
1, and c
1 are defined by 2
-
p
/
3
cos
⁡
(
θ
~
/
3
)
-
(
5
ρ
1
/
9
α
ρ
4
), 2
-
p
/
3
cos
⁡
(
θ
~
/
3
+
2
π
/
3
)
-
(
5
ρ
1
/
9
α
ρ
4
), 2
-
p
/
3
cos
⁡
(
θ
~
/
3
+
4
π
/
3
)
-
(
5
ρ
1
/
9
α
ρ
4
) with p
=
-
(
5
/
27
)
(
-
18
ρ
4
+
5
ρ
1
2
)
/
(
ρ
4
2
α
2
)
<
0, q
=
-
(
10
/
729
)
(
135
ρ
1
ρ
4
+
729
c
ρ
4
2
-
25
ρ
1
3
-
729
ρ
4
2
)
/
(
ρ
4
3
α
3
), θ
~
=
arccos
[
-
(
q
/
2
)
(
-
p
3
/
27
)
]
-
1
/
2. We always write a
1
>
b
1
>
c
1; that is, the largest root is denoted by a
1 among the three real roots. Taking ϕ
(
0
)
=
a
1 as initial value and then integrating (24) once, we obtain a traveling wave solution as follows:
(25)
1
α
1
2
[
F
(
ψ
,
k
1
)
-
(
α
1
2
-
1
)
Π
(
ψ
,
α
1
2
,
k
1
)
]
+
α
ρ
2
a
1
-
c
1
F
(
ψ
,
k
1
)
=
α
-
α
ρ
4
10
β
|
ξ
|
,
where ψ
=
arcsin
(
η
-
a
1
)
/
(
η
-
b
1
), α
1
=
b
1
/
a
1, k
1
=
(
b
1
-
c
1
)
/
(
a
1
-
c
1
), and η
(
x
,
t
)
=
ϕ
(
ξ
) with ξ
=
x
-
c
t.(e) Whenρ
4
>
0, c
≠
1, and equation (
c
-
1
)
-
(
1
/
3
)
α
ϕ
-
(
1
/
6
)
ρ
1
α
2
ϕ
2
-
(
1
/
10
)
ρ
4
α
3
ϕ
3
=
0 has three real roots a
1, b
1, and c
1, substituting (17) into the first equation of (9) yields
(26)
d
ϕ
ϕ
X
5
(
ϕ
)
+
1
2
α
ρ
2
d
ϕ
X
5
(
ϕ
)
=
±
α
α
ρ
4
10
β
d
ξ
,
where X
5
(
ϕ
)
=
(
a
1
-
ϕ
)
(
ϕ
-
b
1
)
(
ϕ
-
c
1
) and a
1, b
1, and c
1 are given above. Taking ϕ
(
0
)
=
a
1 as initial value and then integrating (26) once, we obtain a traveling wave solution as follows:
(27)
2
a
1
a
1
-
c
1
Π
(
φ
,
α
2
2
,
k
2
)
+
α
ρ
2
a
1
-
c
1
F
(
φ
,
k
2
)
=
α
α
ρ
4
10
β
|
ξ
|
,
where φ
=
arcsin
(
a
1
-
η
)
/
(
a
1
-
b
1
), α
2
=
(
a
1
-
b
1
)
/
a
1, k
2
=
(
a
1
-
b
1
)
/
(
a
1
-
c
1
), and η
(
x
,
t
)
=
ϕ
(
ξ
) with ξ
=
x
-
c
t.Case 2.
Under the parametric conditionsm
=
1, h
1
=
0, ρ
5
=
(
1
/
4
)
ρ
2
2, and ρ
4
=
0, (13) can be reduced to
(28)
y
=
±
ϕ
(
c
-
1
)
-
(
1
/
3
)
α
ϕ
-
(
1
/
6
)
ρ
1
α
2
ϕ
2
β
(
1
+
(
1
/
2
)
α
ρ
2
ϕ
)
.
Substituting (28) and ρ
5
=
(
1
/
4
)
ρ
2
2 into the first equation of (9) yields
(29)
d
ϕ
ϕ
X
6
(
ϕ
)
+
1
2
α
ρ
2
d
ϕ
X
6
(
ϕ
)
=
±
1
β
d
ξ
,
where X
6
(
ϕ
)
=
a
~
+
b
~
ϕ
+
c
~
ϕ
2 with a
~
=
(
c
-
1
), b
~
=
-
(
1
/
3
)
α, and c
~
=
-
(
1
/
6
)
ρ
1
α
2. Write Δ
=
4
a
~
c
~
-
b
~
2
=
(
1
/
9
)
α
2
[
6
ρ
1
(
c
-
1
)
+
1
].(a) Whena
~
>
0, c
~
>
0, Δ
>
0 (i.e., c
>
1, -
(
1
/
6
(
c
-
1
)
)
<
ρ
1
<
0), integrating (29) and setting the integral constant as zero, we obtain two exact traveling wave solutions of implicit function type as follows:
(30)
-
1
a
~
Arsh
(
2
a
~
+
b
~
η
η
Δ
)
+
α
ρ
2
2
c
~
Arsh
(
2
c
~
η
+
b
~
Δ
)
=
±
ξ
β
,
where a
~, b
~, c
~, and Δ are given above and η
(
x
,
t
)
=
ϕ
(
ξ
) with ξ
=
x
-
c
t; the “Arsh” is an inverse function of the hyperbolic-sine function, that is, sinh
⁡
-
1.(b) Whena
~
=
0, c
~
<
0, Δ
>
0 (i.e., c
=
1, ρ
1
>
0), integrating (29) and setting the integral constant as zero, we obtain two exact traveling wave solutions of implicit function type as follows:
(31)
2
c
~
η
+
b
~
b
~
=
±
sin
(
-
c
~
β
ξ
±
α
ρ
2
-
c
~
b
~
η
+
c
~
η
2
b
~
η
)
,
where η
(
x
,
t
)
=
ϕ
(
ξ
) with ξ
=
x
-
t.Case 3.
Under the parametric conditionsm
=
1, h
1
=
0, c
=
(
ρ
4
/
10
ρ
2
ρ
5
)
+
1
-
(
1
/
3
ρ
2
), and ρ
1
=
(
2
ρ
5
/
ρ
2
)
-
(
3
ρ
4
/
5
ρ
2
)
+
(
3
ρ
2
ρ
4
/
5
ρ
5
), (13) can be reduced to
(32)
y
=
±
ϕ
3
ρ
4
-
10
ρ
5
30
β
ρ
2
ρ
5
-
α
ρ
4
10
β
ρ
5
ϕ
.(a) When(
(
3
ρ
4
-
10
ρ
5
)
/
ρ
2
ρ
5
)
>
0 and ρ
4
≠
0, substituting (32) into the first equation of (9) to integrate, we obtain a smooth solitary wave solution as follows:
(33)
η
(
x
,
t
)
≡
ϕ
(
ξ
)
=
-
3
ρ
4
-
10
ρ
5
3
α
ρ
2
ρ
4
sec
h
2
(
1
2
3
ρ
4
-
10
ρ
5
30
β
ρ
2
ρ
5
ξ
)
,
where ξ
=
x
-
[
(
ρ
4
/
10
ρ
2
ρ
5
)
+
1
-
(
1
/
3
ρ
2
)
]
t.(b) When(
(
3
ρ
4
-
10
ρ
5
)
/
ρ
2
ρ
5
)
<
0 and ρ
4
≠
0, substituting (32) into the first equation of (9) to integrate, we obtain a periodic wave solution as follows:
(34)
η
(
x
,
t
)
≡
ϕ
(
ξ
)
=
-
3
ρ
4
-
10
ρ
5
3
α
ρ
2
ρ
4
sec
2
(
1
2
10
ρ
5
-
3
ρ
4
30
β
ρ
2
ρ
5
ξ
)
,
where ξ
=
x
-
[
(
ρ
4
/
10
ρ
2
ρ
5
)
+
1
-
(
1
/
3
ρ
2
)
]
t.(c) Whenρ
2
<
0 and ρ
4
=
0, (32) becomes
(35)
y
=
±
1
-
3
β
ρ
2
ϕ
,
(
ϕ
>
0
)
.
Taking ϕ
(
0
)
=
1 as initial value, substituting (35) into the first equation of (9) to integrate, we obtain a peakon solution as follows:
(36)
η
(
x
,
t
)
=
e
-
(
1
/
-
3
β
ρ
2
)
|
x
-
(
1
-
(
1
/
3
ρ
2
)
)
t
|
.Case 4.
Under the parametric conditionsm
=
2, h
2
=
0, and ρ
4
=
0, (14) can be reduced to
(37)
y
=
±
ϕ
(
c
-
1
)
+
a
1
ϕ
+
a
2
ϕ
2
+
a
3
ϕ
3
+
a
4
ϕ
4
β
(
1
+
α
ρ
2
ϕ
+
α
2
ρ
5
ϕ
2
)
,
where a
1, a
2, a
3, and a
4 are given above. Substituting (37) into the first equation of (9) yields
(38)
d
ϕ
ϕ
R
1
(
ϕ
)
+
α
ρ
2
d
ϕ
R
1
(
ϕ
)
+
α
2
ρ
5
ϕ
d
ϕ
R
1
(
ϕ
)
=
±
d
ξ
β
,
where R
1
(
ϕ
)
=
(
c
-
1
)
+
a
1
ϕ
+
a
2
ϕ
2
+
a
3
ϕ
3
+
a
4
ϕ
4. In fact, the cases of roots of the equation R
1
(
ϕ
)
=
0 determine the forms of solutions of (38); different kinds of cases of roots correspond to different kinds of solutions of (38). However, the expressions of the roots of the equation R
1
(
ϕ
)
=
0 are very complex, so we omit these expressions. Of course, the roots of the equation R
1
(
ϕ
)
=
0 can be solved once the parameters α, ρ
1, ρ
2, and ρ
5 are fixed concretely. For example, we can obtain four real roots 2.919208358, -
1.821216553, -
3.236247427, and -
10.46174438 of the equation R
1
(
ϕ
)
=
0 when α
=
0.5, ρ
1
=
6.0, ρ
2
=
2.0, ρ
5
=
0.4, and c
=
4 by using computer.(1) Whenρ
1
ρ
5
>
0 and the equation R
1
(
ϕ
)
=
0 has four real roots ϕ
1, ϕ
2, ϕ
3, and ϕ
4, respectively, taking ϕ
(
0
)
=
ϕ
1, ϕ
2, ϕ
3, ϕ
4 as initial value and then integrating (38), we obtain four kinds of exact traveling wave solutions of implicit function type as follows:
(39)
-
Ω
1
Π
(
μ
,
κ
1
,
r
)
+
ω
(
1
ϕ
4
+
α
ρ
2
+
α
2
ρ
5
ϕ
4
)
F
(
μ
,
r
)
-
Ω
2
Π
(
μ
,
κ
2
,
r
)
=
1
3
α
2
ρ
1
ρ
5
β
|
ξ
|
,
(40)
Ω
3
Π
(
λ
,
κ
3
,
r
)
+
ω
(
1
ϕ
3
+
α
ρ
2
+
α
2
ρ
5
ϕ
3
)
F
(
λ
,
r
)
+
Ω
4
Π
(
λ
,
κ
4
,
r
)
=
1
3
α
2
ρ
1
ρ
5
β
|
ξ
|
,
(41)
Ω
1
Π
(
γ
,
κ
5
,
r
)
+
ω
(
1
ϕ
2
+
α
ρ
2
+
α
2
ρ
5
ϕ
2
)
F
(
γ
,
r
)
+
Ω
2
Π
(
γ
,
κ
6
,
r
)
=
1
3
α
2
ρ
1
ρ
5
β
|
ξ
|
,
(42)
-
Ω
3
Π
(
ϑ
,
κ
7
,
r
)
+
ω
(
1
ϕ
1
+
α
ρ
2
+
α
2
ρ
5
ϕ
1
)
F
(
ϑ
,
r
)
-
Ω
4
Π
(
ϑ
,
κ
8
,
r
)
=
1
3
α
2
ρ
1
ρ
5
β
|
ξ
|
,
where ω
=
2
/
(
ϕ
1
-
ϕ
3
)
(
ϕ
2
-
ϕ
4
), Ω
1
=
ω
(
ϕ
1
-
ϕ
4
)
/
ϕ
1
ϕ
4, Ω
2
=
ω
α
2
ρ
5
(
ϕ
4
-
ϕ
1
), Ω
3
=
ω
(
ϕ
3
-
ϕ
2
)
/
ϕ
2
ϕ
3, Ω
4
=
ω
α
2
ρ
5
(
ϕ
2
-
ϕ
3
), μ
=
arcsin
(
ϕ
2
-
ϕ
4
)
(
ϕ
1
-
η
)
/
(
ϕ
1
-
ϕ
2
)
(
η
-
ϕ
4
), and η is limited by ϕ
1
≥
η
>
ϕ
2
>
ϕ
3
>
ϕ
4; λ
=
arcsin
⁡
(
ϕ
1
-
ϕ
3
)
(
η
-
ϕ
2
)
/
(
ϕ
1
-
ϕ
2
)
(
η
-
ϕ
4
) and η is limited by ϕ
1
≥
η
>
ϕ
2
>
ϕ
3
>
ϕ
4; γ
=
arcsin
(
ϕ
2
-
ϕ
4
)
(
ϕ
3
-
η
)
/
(
ϕ
3
-
ϕ
4
)
(
ϕ
2
-
η
) and η is limited by ϕ
1
>
ϕ
2
>
ϕ
3
>
η
≥
ϕ
4; ϑ
=
arcsin
(
ϕ
1
-
ϕ
3
)
(
η
-
ϕ
4
)
/
(
ϕ
3
-
ϕ
4
)
(
ϕ
1
-
η
) and η is limited by ϕ
1
>
ϕ
2
>
ϕ
3
≥
η
>
ϕ
4, and κ
1
=
ϕ
4
(
ϕ
2
-
ϕ
1
)
/
ϕ
1
(
ϕ
2
-
ϕ
4
), κ
2
=
(
ϕ
2
-
ϕ
1
)
/
(
ϕ
2
-
ϕ
4
), κ
3
=
ϕ
3
(
ϕ
1
-
ϕ
2
)
/
ϕ
2
(
ϕ
1
-
ϕ
3
), κ
4
=
(
ϕ
1
-
ϕ
2
)
/
(
ϕ
1
-
ϕ
3
), κ
5
=
ϕ
2
(
ϕ
3
-
ϕ
4
)
/
ϕ
3
(
ϕ
2
-
ϕ
4
), κ
6
=
(
ϕ
3
-
ϕ
4
)
/
(
ϕ
2
-
ϕ
4
), κ
7
=
ϕ
1
(
ϕ
4
-
ϕ
3
)
/
ϕ
4
(
ϕ
1
-
ϕ
3
), κ
8
=
(
ϕ
4
-
ϕ
3
)
/
(
ϕ
1
-
ϕ
3
), and r
=
(
ϕ
1
-
ϕ
2
)
(
ϕ
3
-
ϕ
4
)
/
(
ϕ
1
-
ϕ
3
)
(
ϕ
2
-
ϕ
4
).(2) Whenρ
1
ρ
5
<
0 and the equation R
1
(
ϕ
)
=
0 has four real roots φ
1, φ
2, φ
3, and φ
4, respectively, taking ϕ
(
0
)
=
φ
1, φ
2, φ
3, φ
4 as initial value and then integrating (38), we obtain another four kinds of exact traveling wave solutions of implicit function type as follows:
(43)
Q
1
Π
(
ν
,
m
1
,
q
)
+
ω
(
1
φ
2
+
α
ρ
2
+
α
2
ρ
5
φ
2
)
F
(
ν
,
q
)
+
Q
2
Π
(
ν
,
m
2
,
q
)
=
1
3
α
2
-
ρ
1
ρ
5
β
|
ξ
|
,
-
Q
1
Π
(
ϵ
,
m
3
,
q
)
+
ω
(
1
φ
1
+
α
ρ
2
+
α
2
ρ
5
φ
1
)
F
(
ϵ
,
q
)
-
Q
2
Π
(
ϵ
,
m
4
,
q
)
=
1
3
α
2
-
ρ
1
ρ
5
β
|
ξ
|
,
Q
3
Π
(
δ
,
m
5
,
q
)
+
ω
(
1
φ
4
+
α
ρ
2
+
α
2
ρ
5
φ
4
)
F
(
δ
,
q
)
+
Q
4
Π
(
δ
,
m
6
,
q
)
=
1
3
α
2
-
ρ
1
ρ
5
β
|
ξ
|
,
-
Q
3
Π
(
ζ
,
m
7
,
q
)
+
ω
(
1
φ
3
+
α
ρ
2
+
α
2
ρ
5
φ
3
)
F
(
ζ
,
q
)
-
Q
4
Π
(
ζ
,
m
8
,
q
)
=
1
3
α
2
-
ρ
1
ρ
5
β
|
ξ
|
,
where Q
1
=
ω
(
φ
2
-
φ
1
)
/
φ
1
φ
2, Q
2
=
ω
α
2
ρ
5
(
φ
1
-
φ
2
), Q
3
=
ω
(
φ
4
-
φ
3
)
/
φ
3
φ
4, Q
4
=
ω
α
2
ρ
5
(
φ
3
-
φ
4
), ν
=
arcsin
(
φ
2
-
φ
4
)
(
η
-
φ
1
)
/
(
φ
1
-
φ
4
)
(
η
-
φ
2
), and η is defined by φ
1
≥
η
>
φ
2
>
φ
3
>
φ
4; ϵ
=
arcsin
(
φ
1
-
φ
3
)
(
φ
2
-
η
)
/
(
φ
2
-
φ
3
)
(
φ
1
-
η
) and η is defined by φ
1
>
φ
2
>
η
≥
φ
3
>
φ
4; δ
=
arcsin
(
φ
2
-
φ
4
)
(
η
-
φ
3
)
/
(
φ
2
-
φ
3
)
(
η
-
φ
4
) and ηis defined by φ
1
>
φ
2
≥
η
>
φ
3
>
φ
4; ζ
=
arcsin
(
φ
1
-
φ
3
)
(
φ
4
-
η
)
/
(
φ
1
-
φ
4
)
(
φ
3
-
η
) and η is defined by φ
1
>
φ
2
>
φ
3
>
φ
4
>
η, and m
1
=
φ
2
(
φ
1
-
φ
4
)
/
φ
1
(
φ
2
-
φ
4
), m
2
=
(
φ
1
-
φ
4
)
/
(
φ
2
-
φ
4
), m
3
=
φ
1
(
φ
2
-
φ
3
)
/
φ
2
(
φ
1
-
φ
3
), m
4
=
(
φ
2
-
φ
3
)
/
(
φ
1
-
φ
3
), m
5
=
φ
4
(
φ
2
-
φ
3
)
/
φ
3
(
φ
2
-
φ
4
), m
6
=
(
φ
2
-
φ
3
)
/
(
φ
2
-
φ
4
), m
7
=
φ
3
(
φ
1
-
φ
4
)
/
φ
4
(
φ
1
-
φ
3
), m
8
=
(
φ
1
-
φ
4
)
/
(
φ
1
-
φ
3
), and q
=
(
φ
2
-
φ
3
)
(
φ
1
-
φ
4
)
/
(
φ
1
-
φ
3
)
(
φ
2
-
φ
4
).Case 5.
Under the parametric conditionsm
=
2, h
2
=
0, c
=
(
3
ρ
4
/
56
ρ
2
ρ
5
)
+
1
-
(
1
/
4
ρ
2
), and ρ
1
=
(
9
ρ
5
/
4
ρ
2
)
-
(
27
ρ
4
/
56
ρ
2
)
+
(
15
ρ
2
ρ
4
/
28
ρ
5
), (14) can be reduced to
(44)
y
=
±
ϕ
3
ρ
4
-
14
ρ
5
56
β
ρ
2
ρ
5
-
α
ρ
4
14
β
ρ
5
ϕ
.(a) When(
3
ρ
4
-
14
ρ
5
)
/
ρ
2
ρ
5
>
0 and ρ
4
≠
0, substituting (44) into the first equation of (9) to integrate, we obtain a smooth solitary wave solution as follows:
(45)
η
(
x
,
t
)
≡
ϕ
(
ξ
)
=
-
3
ρ
4
-
14
ρ
5
4
α
ρ
2
ρ
4
sec
h
2
(
1
2
3
ρ
4
-
14
ρ
5
56
β
ρ
2
ρ
5
ξ
)
,
where ξ
=
x
-
[
(
3
ρ
4
/
56
ρ
2
ρ
5
)
+
1
-
(
1
/
4
ρ
2
)
]
t.(b) When(
3
ρ
4
-
14
ρ
5
/
ρ
2
ρ
5
)
<
0 and ρ
4
≠
0, substituting (44) into the first equation of (9) to integrate, we obtain a periodic wave solution as follows:
(46)
η
(
x
,
t
)
≡
ϕ
(
ξ
)
=
-
3
ρ
4
-
14
ρ
5
4
α
ρ
2
ρ
4
sec
2
(
1
2
14
ρ
5
-
3
ρ
4
56
β
ρ
2
ρ
5
ξ
)
,
where ξ
=
x
-
[
(
3
ρ
4
/
56
ρ
2
ρ
5
)
+
1
-
(
1
/
4
ρ
2
)
]
t.(c) Whenρ
2
<
0 and ρ
4
=
0, (44) becomes
(47)
y
=
±
1
2
-
β
ρ
2
ϕ
,
(
ϕ
>
0
)
.
Taking ϕ
(
0
)
=
1 as initial value, substituting (47) into the first equation of (9) to integrate, we obtain a peakon solution as follows:
(48)
η
(
x
,
t
)
=
e
-
(
1
/
2
-
β
ρ
2
)
|
x
-
(
1
-
(
1
/
4
ρ
2
)
)
t
|
.Case 6.
Under the parametric conditionsm
=
3, h
3
=
0, c
=
(
ρ
4
/
30
ρ
2
ρ
5
)
+
1
-
(
1
/
5
ρ
2
), and ρ
1
=
(
12
ρ
5
/
5
ρ
2
)
-
(
2
ρ
4
/
5
ρ
2
)
+
(
ρ
2
ρ
4
/
2
ρ
5
), (15) can be reduced to
(49)
y
=
±
ϕ
3
ρ
4
-
18
ρ
5
90
β
ρ
2
ρ
5
-
α
ρ
4
18
β
ρ
5
ϕ
.(a) When(
(
3
ρ
4
-
18
ρ
5
)
/
ρ
2
ρ
5
)
>
0 and ρ
4
≠
0, substituting (49) into the first equation of (9) to integrate, we obtain a smooth solitary wave solution as follows:
(50)
η
(
x
,
t
)
≡
ϕ
(
ξ
)
=
-
3
ρ
4
-
18
ρ
5
5
α
ρ
2
ρ
4
sech
2
⁡
(
1
2
3
ρ
4
-
18
ρ
5
90
β
ρ
2
ρ
5
ξ
)
,
where ξ
=
x
-
[
(
ρ
4
/
30
ρ
2
ρ
5
)
+
1
-
(
1
/
5
ρ
2
)
]
t.(b) When(
(
3
ρ
4
-
18
ρ
5
)
/
ρ
2
ρ
5
)
<
0 and ρ
4
≠
0, substituting (49) into the first equation of (9) to integrate, we obtain a periodic wave solution as follows:
(51)
η
(
x
,
t
)
≡
ϕ
(
ξ
)
=
-
3
ρ
4
-
18
ρ
5
5
α
ρ
2
ρ
4
sec
2
⁡
(
1
2
18
ρ
5
-
3
ρ
4
90
β
ρ
2
ρ
5
ξ
)
,
where ξ
=
x
-
[
(
ρ
4
/
30
ρ
2
ρ
5
)
+
1
-
(
1
/
5
ρ
2
)
]
t.(c) Whenρ
2
<
0 and ρ
4
=
0, (49) becomes
(52)
y
=
±
1
-
5
β
ρ
2
ϕ
,
(
ϕ
>
0
)
.
Taking ϕ
(
0
)
=
1 as initial value, substituting (52) into the first equation of (9) to integrate, we obtain a peakon solution as follows:
(53)
η
(
x
,
t
)
=
e
-
(
1
/
-
5
β
ρ
2
)
|
x
-
(
1
-
(
1
/
5
ρ
2
)
)
t
|
.Case 7.
Under the parametric conditionsm
=
4, h
4
=
0, c
=
(
3
ρ
4
/
132
ρ
2
ρ
5
)
+
1
-
(
1
/
6
ρ
2
), and ρ
1
=
(
5
ρ
5
/
2
ρ
2
)
-
(
15
ρ
4
/
44
ρ
2
)
+
(
21
ρ
2
ρ
4
/
44
ρ
5
), (16) can be reduced to
(54)
y
=
±
ϕ
3
ρ
4
-
22
ρ
5
132
β
ρ
2
ρ
5
-
α
ρ
4
22
β
ρ
5
ϕ
.(a) When(
(
3
ρ
4
-
22
ρ
5
)
/
ρ
2
ρ
5
)
>
0 and ρ
4
≠
0, substituting (54) into the first equation of (9) to integrate, we obtain a smooth solitary wave solution as follows:
(55)
η
(
x
,
t
)
≡
ϕ
(
ξ
)
=
-
3
ρ
4
-
22
ρ
5
6
α
ρ
2
ρ
4
sec
h
2
(
1
2
3
ρ
4
-
22
ρ
5
132
β
ρ
2
ρ
5
ξ
)
,
where ξ
=
x
-
[
(
3
ρ
4
/
132
ρ
2
ρ
5
)
+
1
-
(
1
/
6
ρ
2
)
]
t.(b) When(
(
3
ρ
4
-
22
ρ
5
)
/
ρ
2
ρ
5
)
<
0 and ρ
4
≠
0, substituting (54) into the first equation of (9) to integrate, we obtain a periodic wave solution as follows:
(56)
η
(
x
,
t
)
≡
ϕ
(
ξ
)
=
-
3
ρ
4
-
22
ρ
5
6
α
ρ
2
ρ
4
sec
2
(
1
2
22
ρ
5
-
3
ρ
4
132
β
ρ
2
ρ
5
ξ
)
,
where ξ
=
x
-
[
(
3
ρ
4
/
132
ρ
2
ρ
5
)
+
1
-
(
1
/
6
ρ
2
)
]
t.(c) Whenρ
2
<
0 and ρ
4
=
0, (54) becomes
(57)
y
=
±
1
-
6
β
ρ
2
ϕ
,
(
ϕ
>
0
)
.
Taking ϕ
(
0
)
=
1 as initial value, substituting (57) into the first equation of (9) to integrate, we obtain a peakon solution as follows:
(58)
η
(
x
,
t
)
=
e
-
(
1
/
-
6
β
ρ
2
)
|
x
-
(
1
-
(
1
/
6
ρ
2
)
)
t
|
.When the value ofm goes on increasing, the expression (16) becomes more and more complex; thus we cannot obtain exact solutions of (1) as in Cases 1, 2, and 4 under general parameter conditions by integrating this expression. But the exact smooth solitary wave solutions and nonsmooth peakon solutions can always be obtained under certain special parameter conditions. From Cases 3, 6, and 7, by using the mathematical induction, we easily obtain the following results.Under the parametric conditionsh
m
=
0, c
=
(
3
ρ
4
/
(
2
m
+
3
)
(
2
m
+
4
)
ρ
2
ρ
5
)
+
1
-
(
1
/
(
m
+
2
)
ρ
2
), and ρ
1
=
(
3
(
m
+
1
)
ρ
5
/
(
m
+
2
)
ρ
2
)
-
(
9
(
m
+
1
)
ρ
4
/
(
2
m
+
3
)
(
2
m
+
4
)
ρ
2
)
+
(
3
(
m
+
3
)
ρ
2
ρ
4
/
4
(
2
m
+
3
)
ρ
5
) and all-in positive integers m, (16) can be reduced to
(59)
y
=
±
ϕ
3
ρ
4
-
(
4
m
+
6
)
ρ
5
(
2
m
+
3
)
(
2
m
+
4
)
β
ρ
2
ρ
5
-
α
ρ
4
(
4
m
+
6
)
β
ρ
5
ϕ
.(a) When(
3
ρ
4
-
(
4
m
+
6
)
ρ
5
/
ρ
2
ρ
5
)
>
0 and ρ
4
≠
0, substituting (59) into the first equation of (9) to integrate, for all-in positive integers m, we obtain a series of smooth solitary wave solutions as follows:
(60)
η
(
x
,
t
)
≡
ϕ
(
ξ
)
=
-
3
ρ
4
-
(
4
m
+
6
)
ρ
5
(
m
+
2
)
α
ρ
2
ρ
4
×
sec
h
2
(
1
2
3
ρ
4
-
(
4
m
+
6
)
ρ
5
(
2
m
+
3
)
(
2
m
+
4
)
β
ρ
2
ρ
5
ξ
)
,
where ξ
=
x
-
[
(
3
ρ
4
/
(
2
m
+
3
)
(
2
m
+
4
)
ρ
2
ρ
5
)
+
1
-
(
1
/
(
m
+
2
)
ρ
2
)
]
t. Obviously, when m
=
1
,
2
,
3
,
4, the solution (60), respectively, becomes the solutions (33), (45), (50), and (55).(b) When(
(
3
ρ
4
-
(
4
m
+
6
)
ρ
5
)
/
ρ
2
ρ
5
)
<
0 and ρ
4
≠
0, substituting (59) into the first equation of (9) to integrate, for all-in positive integers m, we obtain a series of periodic wave solutions as follows:
(61)
η
(
x
,
t
)
≡
ϕ
(
ξ
)
=
-
3
ρ
4
-
(
4
m
+
6
)
ρ
5
(
m
+
2
)
α
ρ
2
ρ
4
×
sec
2
(
1
2
(
4
m
+
6
)
ρ
5
-
3
ρ
4
(
2
m
+
3
)
(
2
m
+
4
)
β
ρ
2
ρ
5
ξ
)
,
where ξ
=
x
-
[
(
3
ρ
4
/
(
2
m
+
3
)
(
2
m
+
4
)
ρ
2
ρ
5
)
+
1
-
(
1
/
(
m
+
2
)
ρ
2
)
]
t. In particular, when m
=
1
,
2
,
3
,
4, the solution (61), respectively, becomes the solutions (34), (46), (51), and (56).(c) Whenρ
2
<
0 and ρ
4
=
0, (59) becomes
(62)
y
=
±
1
-
(
2
m
+
3
)
(
2
m
+
4
)
β
ρ
2
ϕ
,
(
ϕ
>
0
)
.
Taking ϕ
(
0
)
=
1 as initial value, substituting (62) into the first equation of (9) to integrate, for all-in positive integers m, we obtain a series of peakon solutions as follows:
(63)
η
(
x
,
t
)
=
e
-
(
1
/
-
(
m
+
2
)
β
ρ
2
)
|
x
-
(
1
-
(
1
/
(
m
+
2
)
ρ
2
)
)
t
|
.
Similarly, when m
=
1
,
2
,
3
,
4, the solution (63), respectively, becomes the solutions (36), (48), (53), and (58).In order to intuitively describe and expediently discuss the dynamic properties of the above exact traveling wave solutions of implicit type and explicit type, we draw their profile graphs; see Figures1, 2, 3, and 4 and the discussions below them.The graphs of solutions (19) and (39) under some fixed parametric values.
(a)
Dark broken-soliton
(b)
Bright broken-solitonThe graphs of solution (60) under some fixed parametric values.
(a)
m
=
1
(b)
m
=
4
(c)
m
=
40
(d)
m
=
100The graphs of solution (61) under some fixed parametric values.
(a)
m
=
1
(b)
m
=
4The graphs of solution (63) under some fixed parametric values.
(a)
m
=
1
(b)
m
=
100Whenα
=
0.5, β
=
0.25, ρ
1
=
3, ρ
2
=
2, ρ
4
=
-
2.5, t
=
0.1, x
∈
[
-
3,3
], and η
∈
[
0,10
], we draw graph of solution (19); see Figure 1(a). When α
=
0.5, β
=
0.25, ρ
1
=
6, ρ
2
=
2, ρ
5
=
0.4, c
=
4, t
=
0.01, x
∈
[
-
3,3
], and η
∈
[
0,2.8
], we draw profile graph of solution (39); see Figure 1(b).Figure1(a) shows a shape of dark broken-soliton; it seems that a dark soliton is broken. Figure 1(b) shows a shape of bright broken-soliton; it seems that a bright soliton is broken.Whenα
=
0.5, β
=
0.25, ρ
2
=
-
0.8, ρ
4
=
8, ρ
5
=
-
0.5, t
=
0.01, and m is taking different values, we draw graphs of solution (60); see Figures 2(a)–2(d).Figure2 shows four smooth solitary waves; it is easy to find that their amplitudes and width of waves are different; the amplitude and width of wave decrease as the value of m increases under the same parametric conditions.Whenα
=
0.5, β
=
0.25, ρ
2
=
-
0.8, ρ
4
=
8, ρ
5
=
0.5, t
=
0.01, and m is taking different values, we draw graphs of solution (61); see Figures 3(a) and 3(b).Figure3 shows two periodic waves of blow-up type; it is easy to find that their width of waves is also different; the width of periodic wave increases as the value of m increases under the same parametric conditions.Whenβ
=
0.2, ρ
2
=
-
2, t
=
0.5, and m is taking different values, we draw graphs of solution (63); see Figures 4(a) and 4(b).Figure4 shows two peakon waves; it is easy to find that their amplitudes are the same, but the width of waves and velocities of waves are different; the width of wave increases as the value of m increases; the velocities of waves decrease as the value of m increases under the same parametric condition.


## 4. Conclusion

In this work, by using the integral bifurcation method together with factoring technique, we investigate exact traveling wave solutions of (1) within more extensive parameter ranges. Some exact traveling wave solutions such as broken-soliton solutions (19), (21), (23), (25) (27), (39)–(42), and (43), periodic wave solutions of blow-up type (34), (46), (51), (56), and (61), smooth solitary wave solutions (33), (45), (50), (55), and (60), and nonsmooth peakon solutions (35), (48), (53), (58), and (63) are obtained. It is very worthy to mention solutions (60) and (63); solution (60) denotes a series of smooth solitary waves; solution (63) denotes a series of nonsmooth peakon waves according as the positive integer m increase. Compared with the literature [9], exact solutions obtained in this paper are much more than those in [9]; only one exact solitary wave solution is obtained in [9]. Indeed, we enrich the results of exact solutions based on the works in [8, 9].In addition, the results obtained in this work are very different from those in existing literatures. However, the exact solutions which we obtained in this paper are not much yet because the model studied in this work is very complex equation. Therefore, we hope more and more researchers pay attention to investigations for exact solutions of (1) in future work.



---

*Source: 714214-2014-03-04.xml*