Bifurcation and exact solutions for the (2+1$2+1$)-dimensional conformable time-fractional Zoomeron equation

In this paper, the bifurcation and new exact solutions for the (2 + 1)-dimensional conformable time-fractional Zoomeron equation are investigated by utilizing two reliable methods, which are generalized (G′/G)-expansion method and the integral bifurcation method. The exact solutions of the (2 + 1)-dimensional conformable time-fractional Zoomeron equation are obtained by utilizing the generalized (G′/G)-expansion method, these solutions are classified as hyperbolic function solutions, trigonometric function solutions, and rational function solutions. Giving different parameter conditions, many integral bifurcations, phase portraits, and traveling wave solutions for the equation are obtained via the integral bifurcation method. Graphical representations of different kinds of the exact solutions reveal that the two methods are of significance for constructing the exact solutions of fractional partial differential equation.


Introduction
It is well known that the fractional partial differential equations (FPDEs) [1][2][3][4][5][6][7][8][9][10][11] have received considerable attention [12][13][14][15][16][17][18][19] due to their wide use to describe various complex physical phenomena in the domain of science and engineering. Among the research of the FPDEs, analyzing the bifurcations and the exact traveling wave solutions of FPDEs have been widely investigated as an important subject. Recently, many effective methods have been established and developed to analyze the dynamical behavior of the FPDEs. These methods include the (G /G)-expansion method, the integral bifurcations, the Lie symmetry analysis method, the exp-function method, the Kudryashov method, and so on.
It is worth noting that the (G /G)-expansion method, which was first introduced in [20], has made significant achievements in searching for the exact traveling wave solutions of partial differential equations (PDEs). But the exact solutions of FPDEs have been developed very slowly compared to the exact traveling solutions of PDEs. Most of the methods directly transform an FPDE into an ordinary differential equation by a fractional complex transformation. But the Jumarie's fractional chain rule does not hold. Therefore, the fractional complex transformation cannot be used to obtain the exact traveling solutions of FPDEs when the Riemann-Liouville derivative is used. Recently, Khalil and coworkers [21] introduced the conformable fractional derivative. After that, some scholars [22][23][24][25] have begun to discuss the exact solutions of FPDEs in the sense of the conformable fractional derivative. In this paper, we will introduce the procedure of the generalized (G /G)expansion method for FPDEs, and will discuss the exact traveling wave solutions of the (2 + 1)-dimensional conformable time-fractional Zoomeron equation by the generalized (G /G)-expansion method together with conformable fractional derivative.
The bifurcation method first proposed by Liu and Li [26] is one of the most powerful tools to study the dynamic behavior of PDEs, especially in the analysis of the bifurcation and exact traveling wave solutions [27][28][29][30]. As far as we know, the bifurcation method has not been used to investigate the exact traveling wave solutions of FPDEs in the sense of the conformable fractional derivative. In the paper, we will introduce the procedure of bifurcation approach for constructing the exact traveling wave solutions of FPDEs. By using this method, we will analyze the bifurcation and exact solutions of the (2 + 1)-dimensional conformable time-fractional Zoomeron equation.
The Zoomeron equation is a very convenient model which displays the novel phenomena related with boomerons and trappons, this equation is usually used to describe the evolution of a single scalar field. Recently, Odabasi [31] studied the following (2 + 1)dimensional conformable time-fractional Zoomeron equation: where ∂ α u ∂t α is the conformable fractional derivative of u depending on the variable t. Odabasi applied the modified trial equation method to obtain the exact solutions of the (2 + 1)-dimensional conformable time-fractional Zoomeron equation. Kumar and Kaplan [32] applied the extended exp(-(ξ ))-expansion technique and the exponential rational functional technique to find the explicit and exact solutions of the (2+1)-dimensional conformable time-fractional Zoomeron equation. Hosseini et al. [33] adopted the exp(-(ξ ))expansion approach and modified Kudryashov method to search for the exact solutions of the (2 + 1)-dimensional conformable time-fractional Zoomeron equation.
The main objective of the paper is to employ the generalized (G /G)-expansion method and bifurcation method to construct exact traveling wave solutions of the (2 + 1)dimensional conformable time-fractional Zoomeron equation. The remainder of the article is structured as follows: In Sect. 2, we review the definition of the conformable fractional derivative, and introduce two effective methods for constructing the exact traveling wave solutions of FPDEs. Then in Sect. 3, we discuss the exact solutions of the (2 + 1)dimensional conformable time-fractional Zoomeron equation by using the generalized (G /G)-expansion method and bifurcation method, respectively. Moreover, we obtain the bifurcation and phase portraits of this equation. Finally, we give a brief conclusion in Sect. 4.

The conformable fractional derivative
The definition of the conformable fractional derivative is defined as in [34].
Remark 2.1 The conformable fractional derivative possesses the following properties: The conformable fractional derivative has many important properties. The detailed proof is given in the Appendix.

Description of the methods
Consider the following conformable FPDE: where t, x, y ∈ R, u = u(t, x, y) ∈ R, F is a polynomial in u and its partial fractional-order derivatives. Introduce a traveling wave transformation where k, m and l are arbitrary constants. Equation (2.2) is reduced to the following integer-order ordinary differential equation: where P is a polynomial in u and its derivatives, notation ( ) means the derivative with respect to ξ . If it is possible, we should integrate several times equation (2.4) and take the integral constants as zero.

The generalized (G /G)-expansion method
Step 1. Assume that the solution of equation (2.4) can be expressed as where d is an arbitrary constant, while a i (i = 0, 1, 2, . . . , N) and b i (i = 1, 2, . . . , N) are to be determined later. Then G = G(ξ ) satisfies the following nonlinear ordinary differential equation: where A, B, C, and D are real parameters.
Step 2. The positive integer N can be computed by balancing the highest order derivative and nonlinear term appearing in equation (2.4).

The bifurcation method
Step 1. Let du dξ = y. Equation (2.4) can be transformed into the following two-dimensional system: where R(u, y) is an integral expression.
Step 2. Solve system (2.7) is an integral system, which has the first integral where h is an integral constant.
Step 3. By employing the first integral H(u, y) and analyzing the orbit properties in the phase plane, we can obtain the exact solutions of equation (2.2).

Applications
By employing the transformation u(t, x, y) = u(ξ ), ξ = kx + my + lt α α , the (2 + 1)-dimensional conformable time-fractional Zoomeron equation can be reduced to an ordinary differential equation having form Integrating equation (3.1) twice with respect to ξ , we obtain where ρ is a constant of integration.

Conclusion
By using the generalized (G /G)-expansion method and bifurcation theory method, we obtain exact traveling wave solutions, bifurcation, and phase portraits for the (2 + 1)dimensional conformable time-fractional Zoomeron equation under the given parameter conditions. Many exact solutions have been obtained, which include hyperbolic function solutions, Jacobi elliptic function solutions, trigonometric function solutions, and rational function solutions. Compared with the previous work, the solution method obtained in the paper has not been reported. Furthermore, two methods we employ here can be used to analyze the exact solutions and bifurcation for other FPDE.