Exact Traveling Wave Solutions and Bifurcation of a Generalized ( 3 + 1 )-Dimensional Time-Fractional Camassa-Holm-Kadomtsev-Petviashvili Equation

In this paper, we study the ( 3 + 1 )-dimensional time-fractional Camassa-Holm-Kadomtsev-Petviashvili equation with a conformable fractional derivative. By the fractional complex transform and the bifurcation method for dynamical systems, we investigate the dynamical behavior and bifurcation of solutions of the traveling wave system and seek all possible exact traveling wave solutions of the equation. Furthermore, the phase portraits of the dynamical system and the remarkable features of the solutions are demonstrated via interesting ﬁ gures.


Introduction
Fractional differential equations have wide applications in physics, biology, geology, and other fields. Fu and Yang [1] derived a (3 + 1)-dimensional time-space fractional ZK equation which describes the fractal process of nonlinear dust acoustic waves and obtained some exact solutions via the extended Kudryashov method. The Lump solution of fractional KP equation in dusty plasma is studied in [2]. Guo et al. [3] derived a time-fractional mZK equation to study gravity solitary waves. And some new model equations were derived to study the dynamics of nonlinear Rossby waves in [4,5]. In recent years, many powerful methods for solving fractional differential equations have been proposed, for instance, the first integral method [6,7], the variational iterative method [8][9][10], the exp-function method [11,12], the fractional complex transform [13], and the (G′/G)-expansion method [14][15][16]. It is rather remarkable that the exact traveling wave solutions of the fractional differential equations were obtained by using the bifurcation theory in [17][18][19].
In this paper, we use the conformable fractional derivative proposed by Khalil et al. [27] which is a natural definition and can be applied in optical fiber, physics, and so on. By using the bifurcation method of dynamical systems [17][18][19][28][29][30][31], we investigate the dynamical behavior and bifurcation of solutions of the traveling wave system and try to construct all possible exact traveling wave solutions of Equation (1). This paper is structured as follows. A brief description on the conformable fractional derivative is given in Section 2. In Section 3, by the fractional complex transformation [13], we transform the time-fractional gCH-KP Equation (1) into an ordinary differential equation. Then, applying the bifurcation method of dynamical systems [17][18][19][28][29][30][31], we give the bifurcations and phase portraits of Equation (1). In Section 4, we give all possible exact traveling wave solutions of Equation (1) under different parameter conditions. In Section 5, we state the main conclusion of this paper.

Overview on Fractional Derivative
Since L'Hospital in 1695 asked what does it mean d n f /dx n if n = ð1/2Þ, a lot of definitions of fractional derivative have been introduced in different senses, such as Riemann-Liouville, Caputo, and Weyl fractional derivatives [21]. In this paper, we consider the conformable fractional derivative proposed by Khalil et al. [27]. Given a function f : ð0, +∞Þ ⟶ R. Then, the conformable fractional derivative of f of order α is defined as for all t > 0, α ∈ ð0, 1. And the conformable fractional derivative has the following properties (see [27]). Let α ∈ ð0, 1 and f , g be α-differentiable at a point t > 0. Then, In addition, if f is differentiable, then

Bifurcations and Phase Portraits of the Time-Fractional gCH-KP Equation
Introduce the following fractional transformation [13]: where k, l, m, and n are arbitrary constants. By (3)-(4), it infers Substituting Equation (5) into Equation (1), we have where ′ is the derivative with respect to ξ. Then, we integrate Equation (7) twice and obtain Following [17][18][19][28][29][30], Equation (8) is equivalent to the planar Hamiltonian system with the Hamiltonian where We shall investigate the bifurcations of phase portraits of the system (9) in the ðU, VÞ-phase plane depending on the parameters a, b, c, c 1 , c 2 , k, l, m, n. It is well known that a smooth homoclinic orbit of the system (9) gives rise to a smooth solitary wave solution of Equation (1). For a continuous solution of Uðkx Usually, a homoclinic orbit of the system (9) corresponds to a solitary wave solution of Equation (1), and a heteroclinic orbit of the system (9) corresponds to a kink (or antikink) wave solution. Similarly, a periodic orbit of the system (9) corresponds to a periodic traveling wave solution of Equation (1). Thus, we try to seek all homoclinic orbits, heteroclinic orbits, and periodic orbits of the system (9), which depend on the parameters a, b, c, c 1 , c 2 , k, l, m, n.
In order to study the phase portraits of the system (9), we first study the equilibrium points of the system (9). It is easy to see that when B ≠ 0, the system (9) has two equilibrium points E 0 ð0, 0Þ and E 1 ð−ðB/AÞ, 0Þ. When B = 0, the system (9) has only one equilibrium point E 0 ð0, 0Þ. Let MðU e , V e Þ be the coefficient matrix of the linearized system of the system (9) at an equilibrium point E j ðj = 0, 1Þ. And let J = det ðMðU e , V e ÞÞ. We have According to the bifurcation method of dynamical systems [28][29][30][31] and the above analysis, it infers the following. Case 1. When B > 0, the origin E 0 ð0, 0Þ is a saddle point and E 1 ð−ðB/AÞ, 0Þ is a center point. The bifurcations of phase portraits of system (9) are shown in Figures 1(a) and 1(b).

Explicit Parametric Expressions of the Solutions of Equation (1)
In this section, by using the elliptic integral theory in Byrd and Fridman [32] and the direct integration method, we give all possible explicit parametric representations of the traveling wave solutions of Equation (1). For convenience, denote Then, using the first equation of the system (9) and Equation (13), we obtain the parametric expression of the smooth solitary wave solution of valley type as follows: where ξ = kx + ly + mz − ðn/Γð1 + αÞÞt α , A = ðb/2cnkÞ, and B = ðak 2 − nk + c 1 l 2 + c 2 m 2 Þ/cnk 3 . The profile of the valley solitary wave solution (14) is shown in Figure 4(a) with A = 1, B = 2, k = 1, n = 1, α = 0:5, and y = z = 0.

Journal of Function Spaces
Then, using the first equation of the system (9) and Equation (20), we obtain the parametric expression of the smooth solitary wave solution of valley type as follows: where ξ = kx + ly + mz − ðn/Γð1 + αÞÞt α .
(2) If B > 0, A < 0, corresponding to the homoclinic orbit to the equilibrium point E 1 ð−ðB/AÞ, 0Þ defined by H ðU, VÞ = h 1 , Equation (1) has a smooth solitary wave solution of peak type shown in Figure 2(b). The expressions of the homoclinic orbits defined by H ðU, VÞ = h 1 can be written: Combining dU/V = dξ, we obtain the exact solitary wave solution of peak type as follows: where ξ = kx + ly + mz − ðn/Γð1 + αÞÞt α . In this case, the parametric expression of the periodic solutions is similar to Equations (18) and (19). Figure 3)

Conclusion
In this study, using the fractional complex transform [13] and the bifurcation method of dynamical systems [17][18][19][28][29][30][31], we obtain the bifurcations and phase portraits of the traveling wave system and construct all possible exact traveling wave solutions of Equation (1) with the conformable fractional derivative under different parameter conditions. The profile of the solitary wave solutions, periodic solutions, and periodic cusp wave solutions is given.

Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.

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