Lie Symmetry Analysis and Exact Solutions of Generalized Fractional Zakharov-Kuznetsov Equations

This paper considers the Lie symmetry analysis of a class of fractional Zakharov-Kuznetsov equations. We systematically show the procedure to obtain the Lie point symmetries for the equation. Accordingly, we study the vector fields of this equation. Meantime, the symmetry reductions of this equation are performed. Finally, by employing the obtained symmetry properties, we can get some new exact solutions to this fractional Zakharov-Kuznetsov equation.


Introduction
As one of the important high dimensional nonlinear evolution equations, the Zakharov-Kuznetsov (ZK) equation was first used to discuss the evolution of the propagation of plane waves in magnetized plasma containing cold ions and isothermal electrons [1].The dimensionless form in renormalized variables of the ZK equation is as follows: where u = u(t, x, y) is the normalized electric potential.a, b, c are all normalized constants with respect to different physical meanings.For details of these coefficients, see [1,2].Due to its wide application in maths and physics, the study of this equation was a wide range of theoretical and practical significance.
In paper [12], Blaha et al. considered the following modified version of ZK equation.
where the signs of a represents different physical phenomenon.In addition, Wazwaz [13] did a further study about a nonlinear dispersive modified ZK equation(mZK) as follows: where n ≥ 3 is odd and the sign is either positive or negative.
To include as many physical applications as possible, lots of papers discuss the generalized Zakharov-Kuznetsov equation of the following form [14,15].
Please note that the paper [15] systematically illustrated the detailed group classification algorithm and the process of reduction by discussing Equation (4).
Recently, many important phenomena in various areas of science were well described by fractional order differential equations [16].Due to the realistic senses, a lot of attention has been given to seek solutions of fractional differential equations (FDEs).Various methods such as the homotopy perturbation method (HPM) [17], the variational iteration method (VIM) [18] and the homotopy analysis method (HAM) [19] have been applied for fractional PDEs.Unfortunately, up to now, there are no general methods effective enough to solve the fractional order systems.
Using a new extended trial equation method, the authors of [20] considered the exact solutions of the generalized fractional Zakharov-Kuznetsov equations(FZK(p, q, r)) as follows: where u = u(t, x, y) represents the electrostatic wave potential in plasmas, 0 < α ≤ 1 is the order of fractional derivative.a, b, c are arbitrary constants and the coefficient of a is the nonlinear term, the coefficients of b and c characterize the spatial dispersions in multi-dimensions.p, q, r = 0 are integers.They successfully constructed some new exact solutions, i.e., the elliptic integral function F, Π solutions to (5).In particular, it is important to note that the symmetry analysis of ( 5) is not considered yet in [20].
As one of the most efficient and important methods of studying differential equations, the symmetry group theory has been extensively used to consider the symmetry properties of Zakharov-Kuznetsov equations, see [21][22][23][24][25][26][27] for example.Furthermore, the authors of [28] provided an interesting Appendix of how to proceed in the symmetry analysis of PDEs.Moreover, there have been excellent books about the symmetry analysis, one can be referred to [29][30][31].
However, not the same as it has been done in PDEs, the symmetry Lie group method is not so efficiently used in fractional differential equations (FDEs).To the best of our knowledge, there are lots of studies on the group properties of FDEs.In [32], the authors considered the Lie symmetries of a class of fractional order differential equations with an arbitrary number of independent variables.In [33], analysis of Lie symmetries with conservation laws of a (3+1)-dimensional fractional KdV-Zakharov-Kuznetsov (mKdV-ZK) equation has been considered.Some other results can be found in [34][35][36][37].
When it comes to fractional ZK equations, Lie symmetry analysis, conservation laws and exact solutions for a modified fractional (2+1)-mZK equation was considered in paper [38].The equation reads as Please note that Equation ( 6) is the special case of (5).
In view of the above discussion, we meant to study the group invariant properties of the fractional ZK Equation (5) studied in [20].For the convenience of discussion, system (5) can be rewritten as follows The organization of the rest of this paper is as follows: In Section 2, some preliminary results needed in later sections and some notations, definitions are stated.In Section 3, we establish our main results about Lie symmetries for the fractional ZK Equation (5).In Section 4, examples of group transformations of solutions are considered and new exact solutions are constructed.In Section 5, we study symmetry reductions to the time fractional ZK equation.Finally, we present discussions and conclusions.

Preliminaries
From the viewpoint of Brown motion, the modified Riemann-Liouville(RL) derivative is defined as where n ∈ N, I ν f (t) is the RL fractional integral of order ν, i.e., where Γ(•) is the Gamma function.
For the function u(t, x), a Riemann-Liouville time fractional partial derivative of order µ can be defined as below [35]: In this paper, we also need the following generalized definitions.According to [32,39], the generalized Erd élyi − Kober fractional differential operator is defined as follows: where and the generalized Erd élyi − Kober fractional integral operator reads In the following, we will give a short introduction about how to find Lie point symmetries of fractional order systems.For further information, one can refer to [30,31].
We now consider the following FDE with t, x, y independent variables, u dependent variables.
According to the Lie theory, for some group parameter ε, we need to determine a one parameter Lie group of infinitesimal transformations The associated infinitesimal generator is defined by where According to the infinitesimal invariance criterion [31], Equation ( 13) admits transformation group (14) if and only if the following equation holds By keeping the essential terms, the operator pr (α,3) X takes the following form where Here, the total derivative operator D i is defined by and (x 1 , x 2 , x 3 ) = (t, x, y), (u 1 ) = (u).We see that the explicit expression for the above ones can be obtained in a standard procedure [31].In addition, according to [40], after similar calculation as [35], we can have the explicit expression for φ α,t : where Because we will use FracSym [40] in our paper, as request, hereafter, we only consider symmetries where φ is linear in u (assume µ α t = 0), i.e., φ(t, x, y, u) = uG(t, x, y) + H(t, x, y).
According to (15), by applying the operators pr (α,3) X to Equation ( 13), after splitting the obtained relations by independent variables, we obtain a system of linear PDEs and FDEs by equating these coefficients to zero.Finally, by solving this over-determined system, we can obtain the vector fields X admitted by FDE (13).
Hence, the final symmetry for the time fractional ZK equation is: Finally, the symmetry group of the time fractional ZK equation is given by the follwing vector fields From ( 26), we can find the symmetry generators.They form a closed Lie algebra as shown in Table 1.(26).

Table 1. Commutator table of Lie algebra
Here, the entry in row i and column j means [X i , X j ].It is the commutator for the Lie algebra defined by [X i , X j ] = X i X j − X j X i .

Examples of Group Transformations of Solutions
In this part, by solving the following initial problems, we can get the Lie symmetry group from the related symmetries to get some new exact solutions from the known ones.
Therefore, for the infinitesimal generator X 1 = ∂ ∂x , the corresponding Lie symmetry group are translation along the x-axis g 1 : (t , x , y , u ) → (t, x + ε 1 , y, u), where ε 1 is an arbitrary real number.The group g 1 shows the space-invariance of the equation along the x-axis.Hence, if u = f (t, x, y) is a solution of ( 5), by group g 1 , we can obtain the corresponding new solutions of (5), i.e., For X 2 = ∂ ∂y , the corresponding Lie symmetry group are translation along the y-axis where ε 2 is an arbitrary real number.The group g 2 shows the space-invariance of the equation along the y-axis.Hence, if u = f (t, x, y) is a solution of ( 5), by group g 2 , we can obtain the corresponding new solutions of (5), i.e., In addition, αy ∂ ∂y corresponds to the nonhomogeneous scaling group where ε 3 is an arbitrary real number.The group g 3 is the well-known scaling transformations.Hence, if u = f (t, x, y) is any solution of ( 5), by group g 3 , we can obtain the corresponding new solutions of ( 5), that is, Moreover, the above scaling transformations either expand or contract the size of not only the independent variables but also the dependent ones.In addition, scaling transformation can provide a way to associate the behavior of the solution from different perspectives, for example, a short time solution with large initial values can be rescaled to a longer time solution with small initial values.
To illustrate this, we consider the following example discussed in [20]: Consider p = q = r = n in the Equation ( 5), i.e., The author of [20] obtained two exact solutions of (34): Then, one can obtain the following two new exact solutions of (34) by applying the group (32). where Compared with the existed solutions, we see that these new solutions are the size dilations of not only the independent variables but also the dependent ones.
For the other two symmetry groups, new invariant solutions can be found through existed solutions for the time fractional ZK equation.Maybe more interesting solutions from the physics point of view can be found by applying the full group to the equation.Thus we enrich the former results in [20].

Symmetry Reductions and Exact Solutions to the Time Fractional ZK Equation
In this section, we mainly consider the symmetry reductions to the time fractional ZK Equation.
(i) For the generator X 1 = ∂ ∂x , we have the invariant where τ = t, η = y are the group-invariant.Substituting it into (5) yields the following reduced fractional ODE To solve (38), we need the Laplace transform After taking the inverse transform to F(s), we have the solution, which is where f 4 (η) is an arbitrary function about η, C 0 is a constant.
Hence, for Equation ( 5), we give the following group-invariant solution: Please note that this solution can be viewed as a kind of standing wave solution and it is independent of the space variable x.In addition, due to 0 < α < 1, it decays in time.This solution has not appeared in previous papers.(ii) For the generator X 2 = ∂ ∂y , we have the invariant u = g(τ, ξ), (42) where τ = t, ξ = x are the group-invariant.Substituting it into (5), we obtain the following reduced fractional differential equation D α τ g + apg p−1 g ξ + b (q 3 − 3q 2 + 2q)g q−3 g 3 ξ + (3q 2 − 3q)g q−2 g ξ g ξξ + qg q−1 g ξξξ = 0. (43) (iii) For the generator where ξ = xt 3p−q−2 are the group-invariant.
In the following, we will do further study on symmetry reductions and exact solutions of (43) in detail.
If Equation ( 43) is invariant under the point transformations with the group parameter , the associated Lie algebra is spanned by in which τ(ξ, τ, g), ξ(ξ, τ, g), ḡ(ξ, τ, g) are to be determined.If the vector fields above generates a symmetry of (43), we obtain the following Lie symmetry condition pr where As with the similar discussion given in previous sections, again by using the Maple package [31,40,41], we can have the symmetry algebra of (43), which is spanned by the following vector fields For V 1 , we obtain the group invariant Inserting it into (43), we have the following reduced fractional ODE which indicates that f 5 (τ) = C 1 τ α−1 , where C 1 is a constant.As a result, we obtain a group-invariant solution of (5) of the form This solution is only related to the time variable.In addition, it also decays in time.For V 2 , we have the group invariant where ξ = ξτ a 0 , τ = t, ξ = x.

Conclusions
In this paper, by using the Lie symmetry analysis method, we have considered the invariance properties of a class of generalized fractional Zakharov-Kuznetsov equations.Lie point symmetries to this equation is performed.The Lie algebra and the symmetry reductions of this fractional Zakharov-Kuznetsov equation are obtained.Finally, some new exact solutions were constructed to the fractional Zakharov-Kuznetsov equation.Although there have been lots of symmetry results for the time fractional Zakharov-Kuznetsov equations, all those models considered can be viewed as the special cases of the one we considered in this paper.Hence, we have extended some existing results.However, as we see, the coefficients a, b and c are all normalized constants according to different physical meanings.To stay closely to the former research and to keep its specific physical meaning of parameters, therefore, we only consider the fractional Zakharov-Kuznetsov equation with constant coefficients.Furthermore, note that in [24], a class of generalized Zakharov-Kuznetsov equations with variable coefficients was considered.Perhaps we will pay our attention to this more generalized model for future studies.