Bifurcation, Traveling Wave Solutions, and Stability Analysis of the Fractional Generalized Hirota–Satsuma Coupled KdV Equations

In this paper, the bifurcation, phase portraits, traveling wave solutions, and stability analysis of the fractional generalized Hirota–Satsuma coupled KdV equations are investigated by utilizing the bifurcation theory. Firstly, the fractional generalized Hirota–Satsuma coupled KdV equations are transformed into two-dimensional Hamiltonian system by traveling wave transformation and the bifurcation theory. )en, the traveling wave solutions of the fractional generalized Hirota–Satsuma coupled KdV equations corresponding to phase orbits are easily obtained by applying the method of planar dynamical systems; these solutions include not only the bell solitary wave solutions, kink solitary wave solutions, anti-kink solitary wave solutions, and periodic wave solutions but also Jacobian elliptic function solutions. Finally, the stability criteria of the generalized Hirota–Satsuma coupled KdV equations are given.


Introduction
In this paper, we consider the fractional generalized Hirota-Satsuma coupled KdV (FGHSCKdV) equations as follows [1][2][3][4]: where u � u(t, x), v � v(t, x), and w � w(t, x) represent unknown functions about time t and space x, respectively, and D α t u and D β x u represent the conformable fractional derivatives of u with respect to time t of order α and space x of order β, respectively. e FGHSCKdV equations are related to most types of long waves, acoustic waves, and planetary waves in geophysical hydrodynamics.
In addition, when α � β � 1, equation (1) is reduced to the well-known generalized Hirota-Satsuma coupled KdV (GHSCKdV) equations as follows [5][6][7][8][9][10][11][12]: where u, v, and w are the unknown functions of t and x. System (2) is usually used to explain the interactions between two long waves with different dispersive relations in physics, which was first proposed in Ref. [13]. When w � 0, system (2) is reduced to the coupled KdV equation, which was first introduced by Satsuma and Hirota [14]. Equation (2) is a very important mathematical physical equation. In recent years, many scholars and experts have devoted themselves to study equation (2), especially the property of orbits, the bifurcation of phase portrait, and the exact solution of these equations. It has become a very important subject in the field of physics and engineering technology because it can further help physicists to explain the dynamic behavior of GHSCKdV equations. Recently, many important methods have been established to construct the traveling wave solutions of equation (2), such as the subequation method, the Hirota direct method, the Hirota bilinear method, plane dynamic system analysis method, and so on (see [5][6][7][8][9][10][11][12][13][14] and the references therein).
Fractional partial differential equation (FPDE), usually used to simulate phenomena in natural science, plays a very important role in the field of applied sciences (such as fluid mechanics, nonlinear optics, and so on). e most commonly used definitions of fractional derivatives include Caputo derivative, Riemann-Liouville derivative, conformable derivative, etc. But, for FGHSCKdV equations, the research progress is very slow due to the complexity of fractional derivative. e main reason is that most of the fractional derivatives do not obey the chain rule. For example, in Ref. [1] and Ref. [2], the researchers applied Jumarie's modified Riemann-Liouville derivative to seek the exact solutions of FGHSCKdV. With the development of fractional calculus, in 2014, Khalil et al. [15] gave a new definition of conformable fractional derivative, which satisfies the chain rule and Leibniz formula. We can easily simplify the FPDE into the nonlinear ordinary differential equation (NLODE). erefore, the method of constructing the exact solution of PDE is also suitable for finding the exact solution of FPDE. In Ref. [4], Ali et al. obtained the exact solutions of the time-fractional GHSCKdV in the sense of the conformable derivative by the subequation method and residual power series method. In Ref [3], although Sirisubtawee and his collaborators studied FGHSCKdV by the (G′/G, 1/G)-expansion method, the obtained solutions only include the hyperbolic function solutions, rational function solutions, and trigonometric function solutions. Based on the lack of research on the FGHSCKdV equations, the main purpose of this paper is to study the property of orbits, the bifurcation of phase portrait, the exact solution, and stability of FGHSCKdV equations.
is paper is organized as follows. In Section 2, phase portraits of system (1) are drawn. In Section 3, traveling wave solutions of system (1) are constructed. In Section 4, the stability criteria of system (1) are given. Finally, some conclusions are discussed in Section 5.

Phase Portraits of System (1)
We first do the following traveling wave transformation: where α and β are nonzero arbitrary constants and c stands for the wave speed.
Substituting (3) into (1), we get the NLODE as follows: In order to eliminate the coupling relationship of the above equations, we assume that where a 1 , b 1 , c 1 , a 2 , and b 2 are constants.
Substituting (5) into (4) and then integrating once, we find that equation (4) has the same result as follows: where d 1 is the integral constant.
Multiplying equation (6) by V ′ and then integrating equation (6), we obtain where d 2 is the integral constant.
us, by the bifurcation theory of planar dynamical system (see [17][18][19][20][21]), we can easily construct the exact solution of equation (1) depending on the parameters A and B.

Remark 1.
In the above discussion, we have obtained the solution v(t, x) of system (1). Obviously, we can get the solutions u(t, x), w(t, x) of system (1) from coupling relation (5).

Conclusion
In this paper, we have drawn the bifurcations of phase portraits of FGHSCKdV equations, and we have constructed the traveling wave solutions of FGHSCKdV equations which include implicit analytical solutions, hyperbolic function solutions, and Jacobian elliptic function solutions. rough a fractional traveling wave transformation, the FGHSCKdV equations can be simplified into coupled nonlinear ordinary differential equations. In order to eliminate the coupling relationship, we assume that the solutions U, V, and W of equation (1) satisfy equation (5). erefore, we can easily obtain the Hamiltonian system of (14) and construct the solutions of equation (1). Compared with the reported literature, the solution of Jacobian functions obtained in the paper has not been published. More importantly, we also give the stability condition of GHSCKdV equations. erefore, the study of the paper has very important value.

Data Availability
No data were used to support this study.

Conflicts of Interest
e authors declare that they have no conflicts of interest. Discrete Dynamics in Nature and Society 5