Phase Portraits and Traveling Wave Solutions of a Fractional Generalized Reaction Duffing Equation ()
1. Introduction
Some famous nonlinear fractional wave equations, such as the fractional Klein-Gordon equation, Landau-Ginzburg-Higgs equation, the fractional
equation, the fractional Duffing equation and the fractional Sine-Gordon equation, can be summarized as the fractional generalized reaction Duffing model. A lot of authors have done a lot of research on the exact solutions of this equation. By using the new ansatz method, the solitary wave solutions and periodic solutions of gRDM were obtained in [1]. Furthermore, the exact soliton solutions of gRDM have been obtained by using the generalized hyperbolic function method, the Bäcklund transformation obtained by the homogeneous balance method, the first integration method of the fractional derivative in the sense of the improved Riemann-Liouville derivative, and the compatible fractional complex transformation method, respectively in [2] [3] [4] [5]. Based on an extended first-type elliptic sub-equation method and its algorithm, the new bell-shaped and kink-shaped solitary wave solutions, triangular periodic wave solutions and singular solutions of gRDM were solved in [6]. The accurate soliton solutions were obtained by using Bäcklund transformation of fractional Riccati equation, function variable method, and general projective Riccati equation [7] [8] [9]. In addition, other authors have used auxiliary function methods, Hermite transformation and Riccati equations, fractional sub-equations and other methods to study the exact solutions of gRDM in [10] [11] [12]. Recently, some new traveling wave solutions of the (2+1)-dimensional time-fractional Zoomeron equation and the superfield gardner equation have been obtained in [13] [14]. The fractional derivatives and fractional derivative equations have been deeply studied in [15] [16] [17] [18]. In this paper, we consider the following fractional order generalized reaction Duffing equation
(1)
where
and s are all real constants,
, and
is defined in Section 2. The following equations are special cases of Equation (1), for example
1) Fractional Klein-Gordon equation
2) Fractional Landau-Ginzburg-Higgs equation
3) Fractional
equation
4) Fractional Duffing equation
5) Fractional Sine-Gordon equation
In this paper, we use the dynamic system approach [19] [20] to study the phase portraits and traveling wave solutions of the Equation (1), and try to construct all possible exact traveling wave solutions of this equation.
The rest of this paper is organized as follows. In Section 2, we introduce some basic definitions and important properties of the fractional derivative. In Section 3, by applying the dynamic system approach [19] [20], we give the phase portraits of the Equation (1). In Section 4, we give all possible exact traveling wave solutions of the Equation (1) under different parameters. In Section 5, we state the main conclusions of this paper.
2. Definition and Properties of the Fractional Derivative
The idea of fractional derivatives originated from the semi-derivative discussed by Leibniz and Lopida in 1695. Subsequently, many authors studied fractional derivatives and formed several different definitions, such as Riemann-Liouville, Caputo and other fractional derivatives. In this section, we introduce the common fractional derivatives proposed by Khalil et al. [21]. Let
. Then, the conformal fractional derivative of f of order
is defined as
(2)
for all
,
. And the conformal fractional derivative has the following properties. Let
, and
be
-differentiable at a point
. Then
(3)
In addition, if f is differentiable, then
(4)
3. Phase Portraits of Equation (1)
Inspired by [22], we introduce the following fractional transformation
(5)
where
are all arbitrary constants. According to (3)-(4), it infers
(6)
By (6), substituting Equation (5) into Equation (1), we get
(7)
where ' is the derivative with respect to
. Furthermore, it follows from [20] [23] that (7) is equivalent to the plane Hamiltonian system
(8)
with the Hamiltonian
where
.
In order to study the phase pictures of the system (8), it is necessary to study the equilibrium points of the system (8). Let
. When
, the system (8) has two equilibrium points
. When
,
the system has three equilibrium points
and
. When
, the system has only one equilibrium point
. Let
be the coefficient matrix of the linearized system of the system (8) at an equilibrium point
. Let
. We have
By the planar dynamical theory [20], the above analysis and Maple, we obtain the following results and the phase portraits.
Case 1.
.
When
,
is a saddle point and
is a cusp point.
When
,
is a center point and
is a cusp point.
The corresponding phase portraits of the system (8) are shown in Figure 1.
Case 2.
.
When
,
and
,
and
are saddle points and
is a center point.
Figure 1. The phase portraits of the system (8). (a)
; (b)
.
When
,
and
,
and
are saddle points and
is a center point.
When
,
,
is a saddle point, and
and
are center points.
When
,
,
is a center point, and
and
are saddle points.
When
,
and
,
and
are center points and
is a saddle point.
When
,
and
,
and
are center points and
is a saddle point.
When
and
,
and
are cusp points and
is a saddle point.
When
and
,
and
are cusp points and
is a center point.
The corresponding phase portraits of the system (8) are shown in Figures 2-6.
Case 3.
When
,
is a saddle point.
When
,
is a center point.
The corresponding phase portraits of the system (8) are shown in Figure 7.
4. Exact Solutions of Equation (1)
We use the elliptic integral theory and direct integration method to give all possible explicit parameter representations of the traveling wave solution of Equation (1). We first denote
(9)
where
,
.
Figure 2. The phase portraits of the system (8). (a)
; (b)
.
Figure 3. The phase portraits of the system (8). (a)
; (b)
.
Figure 4. The phase portraits of the system (8). (a)
; (b)
.
Figure 5. The phase portraits of the system (8). (a)
; (b)
.
4.1. Consider Case 1 in Section 3
By
, it obtains
.
(1) If
,
, corresponding to the homoclinic orbit
defined by
from (9), the Equation (1) has a solution of shown in Figure 1(a). By
, it gets
Figure 6. The phase portraits of the system (8). (a)
; (b)
.
Figure 7. The phase portraits of the system (8). (a)
; (b)
.
(10)
Using the first equation of system (8) and equation (10), we get the following parameter expression
(11)
where
,
,
,
, and
. The solution (11) is shown in Figure 1(a) with
,
,
.
2) If
,
, corresponding to the homoclinic orbit
defined by
, the Equation (1) has a solution of shown in Figure 1(b). By
, it gets
(12)
By (8) and Equation (12), we get
(13)
where
,
,
,
, and
. The solution (13) is shown in Figure 1(b) with
,
,
.
4.2. Consider Case 2 in Section 3
1) If
,
, corresponding to the homoclinic orbit
defined by
, Equation (1) has a solution of shown in Figure 6(a). By
, it infers
(14)
Combining the first equation of system (8) and Equation (14), we have
(15)
where
,
, and
. The solution (15) is shown in Figure 6(a) with
,
.
2) If
,
, corresponding to the homoclinic orbit
defined by
, the Equation (1) has a solution of shown in Figure 6(b). By
, it gets
(16)
Appling (8) and (16), we get
(17)
where
,
, and
. The solution (17) is shown in Figure 6(b) with
,
.
3) If
,
and
or
,
and
, corresponding to the homoclinic orbit
defined by
, Equation (1) has the solution of shown in Figure 2 or Figure 5. By
, it obtains
(18)
It follows from (8) and (18), we obtain
where
,
,
,
, and
.
4) If
,
or
,
corresponding to the homoclinic orbit
defined by
, Equation (1) has the solution of shown in Figure 3 or Figure 4. By
, it gets
(19)
Using the first equation of system (8) and Equation (19), we get the following parameter expression:
where
,
,
,
, and
.
5) If
,
or
,
corresponding to the homoclinic orbit
defined by
has the solution of shown in Figure 2 or Figure 4. It follows from
that
(20)
The relation
holds on the U-axis. Therefore, by using the first equation of system (8) and equation (20), we get the following parameter expression
where
,
and
.
4.3. Consider Case 3 in Section 3
1) If
,
, corresponding to the homoclinic orbit
defined by
, Equation (1) has the solution of shown in Figure 7(b). By
, it gets
(21)
Using the first equation of system (8) and Equation (21), we get the following parameter expression
(22)
where
,
,
,
, and
. The solution (22) is shown in Figure 7(b) with
,
,
.
(2) If
,
, corresponding to the homoclinic orbit
defined by
, Equation (1) has the solution of shown in Figure 7(a). By
, it gets
(23)
Using the first equation of system (8) and Equation (23), we get the following parameter expression
(24)
where
,
,
,
, and
. The solution (24) is shown in Figure 7(a) with
,
,
.
5. Conclusion
In conclusion, we obtained the phase portraits of the traveling wave system by using the fractional complex transformation and the dynamical system method [19] [20]. Moreover, we construct all possible accurate traveling wave solutions of Equation (1) under different parameter conditions.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.