Invariant Analysis, Analytical Solutions, and Conservation Laws for Two-Dimensional Time Fractional Fokker- Planck Equation

The main purpose of this paper is to apply the Lie symmetry analysis method for the two-dimensional time fractional FokkerPlanck (FP) equation in the sense of Riemann–Liouville fractional derivative. The Lie point symmetries are derived to obtain the similarity reductions and explicit solutions of the governing equation. By using the new conservation theorem, the new conserved vectors for the two-dimensional time fractional Fokker-Planck equation have been constructed with a detailed derivation. Finally, we obtain its explicit analytic solutions with the aid of the power series expansion method.


Introduction
Fractional calculus has attracted more attention of many researches in various scientific areas including biology, physics, financial theory, gas dynamics, engineering, fluid mechanics, and other areas of science, see for example [1][2][3][4][5][6][7][8]. The theory of fractional calculus is considered as a generalization of classical differential and integral calculus; it is an excellent tool for describing the memory effect and hereditary properties of various processes and viscoelastic materials. Due to its realistic senses, many researchers have tried to look for exact, analytical, and numerical solutions of fractional partial differential equations using different powerful methods such us ðG′/GÞ-expansion method [9], the Variational iteration method [10,11], functional variable method [12], subequation method [13], Finite difference method [14], Exp function method [15], Homotopy analysis method [16], Adomian decomposition method [17], the First integral method [18], Laplace transform method [19], Sumudu transform method [20], and so many other approaches.
The Lie symmetry method was firstly advocated by the Norwegian mathematician Sophus Lie [21,22], who has made great achievements in the theories of continuous groups and differential equations. It is an efficient approach and widely employed for solving ordinary differential equations (ODEs), partial differential equations (PDEs), and fractional partial differential equations (FPDEs). This popularity is due to its utility in determining the explicit solutions of both ODEs and PDEs, linearization of some nonlinear equations, reducing the order of independent variables, and so on. Many papers focused on constructing symmetries of different fractional differential equations [23,[24][25][26][27]. Furthermore, the concept of conservation laws is fundamental and widely used in the study of the resolution of PDEs. Moreover, they convey a large deal of information about the studied physical system. The new conservation laws were introduced by Ibragimov [28], based on the notion of Lie symmetry generators without Lagrangian for solving FPDEs. Therefore, this new conservation law plays an increasingly important role in solving the conservation laws of FPDEs. More details about conservation laws can be found in [29][30][31][32].
In this paper, we consider the following two-dimensional time fractional Fokker-Planck equation It is well known that the Markovian diffusion process can be described with the Fokker-Planck equation. The Fokker-Planck equation is a partial differential equation for the probability density and the transition probability of these stochastic processes. It plays an important role in control theory, fluid mechanics, astrophysics, and quantum [33]. Moreover, it has been applied in various natural science fields such as quantum optics, solid-state physics, chemical physics, theoretical biology, and circuit theory. It is firstly proposed by Fokker and Planck to characterize the Brownian motion of particles [34]. Many researchers have solved the Fokker-Planck equation using various powerful methods, for more details see [35][36][37].
The fractional derivatives described here are in the Riemann-Liouville sense of order αðα > 0Þ, see [38,39], which is defined by where ΓðzÞ is the Gamma function defined by The main motivation behind this article is to make use of the Lie symmetry method to get the infinitesimal generators, group invariant solutions for the time fractional twodimensional Fokker-Planck equation (FP), and to construct conservation laws given by Ibragimov [28]. Therefore, the new conserved vectors have been obtained using the new conservation theorem. Based on the power series method [32], the explicit power series solutions of the two-dimensional time fractional Fokker-Planck equation are derived. The rest of this paper is organized as follows: in Section 2, we review some basic definitions of the Lie Symmetry method for fractional partial differential equations (FPDEs) and its properties. By employing the proposed method, Lie point symmetries of the Eq. (1) are obtained; by using similarity variables, the reduced equations are obtained; solving some of them, then the similarity solutions of Eq. (1) are deduced in Section 3. In Section 4, the conservation laws of Eq. (1) are obtained. Section 5 is devoted to constructing the explicit analytical power series solutions. Some conclusions and discussions are given in Section 6.

Method of Lie Symmetry Analysis for FPDEs
In this section, we briefly review the main points about Lie symmetry analysis of FPDEs [39][40][41] of the following form We assume that the Eq. (4) is invariant under a oneparameter ε Lie group of infinitesimal transformations which are given asx where ε < <1 is a group parameter and ξ, ζ, η, and τ are infinitesimals and η 0 α is extended infinitesimal. The explicit expressions of η x , η y , η xx , η yy are given by D x , D y , and D t are the total derivatives with respect to x, y , and t, respectively, which are defined as

Journal of Function Spaces
The infinitesimal generator X is given by the following expression The infinitesimal generator X satisfy the following invariance condition of Eq. (4): where The structure of the Riemann-Liouville derivative must be invariant under transformations (5), because the lower limit of the integral (2) is fixed. The invariance condition yields The explicit form of the α th extended infinitesimal can be obtained as follows: where It is worth noting that μ = 0 if the infinitesimal η is linear in u, due to the presence of ∂ k η/∂u k .

Symmetry Analysis and Similarity Reductions of the Two-Dimensional Time Fractional Fokker-Planck Equation
In the present section, the Lie symmetry analysis method has been applied for deriving the infinitesimal generators of the two-dimensional time fractional Fokker-Planck (1). By using the third prolongation [43,44], the symmetry determining equation for Eq. (1) has been obtained as By substituting the expressions η 0 α given in Eq. (6) and Eq. (12) into Eq. (15), and equating various powers of derivatives of u to zero, we obtain an overdetermined system of linear equations; by solving this system, we obtain the following infinitesimals where C i , i = 1, ::, 6 are arbitrary constants. So, the associated vector fields of Eq. (1) are given by

Journal of Function Spaces
Case 2. For X 1 = ∂/∂y, the characteristic equation is By solving the above characteristic equation, we obtain the solution u = f ðx, tÞ. Substituting it into Eq. (1), we derive the following reduced fractional ordinary Fokker-Planck equation: Using the symmetry method, we obtain the following infinitesimals: where M i , i = 1, ::, 4 are arbitrary constants. Then, the Lie algebra of infinitesimal symmetries of Eq. (19) is given by Case 3. The characteristic equation for the infinitesimal generator X 11 can be expressed symbolically as follows: By solving the above characteristic equation, we obtain the solution u = f ðtÞ. Substituting it into Eq. (19), we derive the reduced fractional ordinary equation: The above can be solved through the Laplace transform method Since the Laplace transform of the Riemann-Liouville deriv-ative is defined by the following form then, According to Eq. (25), we have By using the inverse Laplace transform, it gives where is the Mittag-Leffler function.
Case 4. For X 11 + γX 12 , the similarity transformation corresponding to this generator can be derived by solving the associated characteristic equation which take the form replacing it in Eq. (19) yields the following reduced FODE: By using the Laplace transform, we obtain the following solution

Conservation Laws of the Two-Dimensional Time Fractional Fokker-Planck Equation
In this section, the conservation laws of the two-dimensional time fractional Fokker-Planck equation have been investigated by using a new conservation theorem [28]. The conserved vectors C t , C x , C y have been obtained, and it satisfies 4 Journal of Function Spaces the following conservation equation: The formal Lagrangian for Eq. (1) can be written as follows: where ωðx, y, tÞ is the new dependent variable. Based on the definition of the Lagrangian, the action integral of Eq. (35) is given by The Euler-Lagrange operator is defined as The adjoint operator ðD α t Þ * of is defined by where R n−α T is the right-sided operator of fractional integration of order ðn − αÞ that is defined by So, the adjoint equation of Eq. (1) as the Euler-Lagrange equation, given by For the case of three independent variables x, y, t and one dependent variable uðx, y, tÞ, we get where I is the identity operator and δ/δu is denoted as the Euler-Lagrange operator. So X is presented as and the Lie characteristic function W is given by Using the Lie symmetries V 1 , V 2 , V 3 , we have Based on the fractional generalizations of the Noether operators, the components of conserved vectors can be presented as follows: where Jð:Þ is defined by And the other components C i are defined as where ξ 1 = ξ, ξ 2 = ζ, θ = 1, 2. Using Eqs. (45) and (47), we obtain the following components of conserved vectors.

Power Series and Analytical Solutions for Eq. (1)
In this section, based on the power series method [45], the exact analytic solutions are a kind of exact power series solutions for Eq. (1), constructed with a detailed derivation.

Conclusion
In this paper, the invariance properties of the twodimensional time fractional Fokker-Planck equation with the Riemann-Liouville fractional derivative have been investigated in the sense of Lie point symmetries. Then, the power series method has been applied to get an explicit solution for the two-dimensional time fractional Fokker-Planck equation.
For obtaining new components of conserved vectors, a new theorem of conservation law has been employed along with the formal Lagrangian, which allows us to construct conservation laws for the two-dimensional time fractional Fokker-Planck equation. Our results show that the extended Lie group analysis approach and the power series method provide powerful mathematical tools to investigate other FDEs in different fields of applied mathematics. In addition, it shows that the proposed analysis is very efficient to construct conservation laws of the two-dimensional time fractional Fokker-Planck equation. Moreover, we can employ symmetry analysis to the time-space fractional Fokker-Planck equation; it will be valuable as future subject works.

Data Availability
No data were used to support this study.

Conflicts of Interest
The authors declare that they have no conflicts of interest.