A modified invariant subspace method for solving partial differential equations with non-singular kernel fractional derivatives

In this work, the well known invariant subspace method has been modified and extended to solve some partial differential equations involving Caputo-Fabrizio (CF) or Atangana-Baleanu (AB) fractional derivatives. The exact solutions are obtained by solving the reduced systems of constructed fractional differential equations. The results show that this method is very simple and effective for constructing explicit exact solutions for partial differential equations involving new fractional derivatives with nonlocal and non-singular kernels, such solutions are very useful to validate new numerical methods constructed for solving partial differential equations with CF and AB fractional derivatives.


Introduction
Fractional calculus provides an important characteristic to describe the complicated physical phenomena with memory effects. For this reason, the fractional calculus is becoming increasingly used as a modeling tool in physics, engineering and control processing in various fields of sciences such as fluid dynamics, plasma physics, mathematical biology and chemical kinetics, diffusion, etc [1][2][3][4]. Due to their properties, fractional derivatives and integrals make this kind of calculus a good candidate to describe such phenomena. Some fundamental definitions of fractional derivatives were given by Riemann-Liouville and Liouville-Caputo [5][6][7][8][9]. Recently, Caputo and Fabrizio defined a new fractional derivative without singular kernel [10] named Caputo-Fabrizio derivative with specific properties, the derivative of a constant is zero and the initial conditions used in the fractional differential equations having a physical interpretation. Later, Atangana and Baleanu proposed another fractional derivative with non-local and non-singular kernel named Atangana-Baleanu derivative [11]. Besides, seeking exact solutions of fractional partial differential equations is not an easy task, and it's remain a relevant problem. Therefore, many powerful methods have been proposed for solving analytically the fractional partial differential equations. Such methods include; Homotopy Perturbation Method [12], Homotopy Perturbation coupled with Sumudu Transform [13], Adomian Mecomposition Method [14], Variational Iteration Method [15], Fractional Iteration Method [16], etc. On the author hand, recent investigations show that the invariant subspace method, developed by V.A. Galaktionov and S.R. Svirshchevski [17], is an effective tool to construct exact solutions of some fractional partial differential equations with Caputo fractional derivative. R.Sahadevan and P.Prakash [18] used invariant subspace method to derive exact solutions of certain time fractional nonlinear partial differential equations, Hashemi [19] also adopted the same method to solve partial differential equations with conformable derivatives, Choudhary et al. [20] used this technique to explore solutions of some fractional differential equations, etc. In the present paper, we present a modified version of the invariant subspace method which does not require any use of the Laplace transformation. We then make use of this novel technique to solve some fractional partial differential equations using fractional operators of Caputo-Fabrizio and also Atangana-Baleanu type. The exact solutions of these equations are obtained by solving the reduced systems constructed from the studied equations.
The laout of the paper is organized as follows: In section 2, we present some basic definitions of fractional derivatives and integrals. Section 3 describes the modified invariant subspace method. Construction of exact solutions to some partial differential equations with Caputo-Fabrizio and Atangan-Baleanu derivatives is presented in section 4. Finally, concluding remarks are given in section 5.

Fractional Calculus tools
In this section, we present some important defnitions and mathematical concepts on fractional derivatives with nonsingular kernels and related tools. Definition 1. The Mittag-Leffler function E α [21], is defined as where z is a complex variable, α ∈ C and ℜ(α) > 0.
This function arises naturally in the solution of fractional order integral equations or fractional order differential equations. It interpolate between a purely exponential law and power-law like behavior of phenomena governed by ordinary kinetic equations and their fractional counterparts. On the other hand, Caputo and Fabrizio [10] developed a new fractional derivative as follows Definition 2. Let u be a function in H 1 (a, b), b > a et 0 < α < 1 then, the new Caputo-Fabrizio derivative of fractional order α is defined as [10] where M(α) is a normalization function satisfying M(0) = M(1) = 1.
From [10], we recall that if the function u does not belong to H 1 (a; b) then, CF derivative can be writted as The fractional integral operator associated to the CF fractional derivative is expressed as It's clear that the Caputo-Fabrizio derivative has no singular kernel, since the kernel is based on exponential function.
Recently, Atangana and Baleanu proposed a new fractional derivative wich has non-local and non-singular kernel based on the generalized Mittag-Leffler function. More recently, they claimed that there is two general definitions of their derivative in the Riemann-Liouville and Caputo sense. Moreover, this fractional derivative has a fractional integral as an anti-derivative of their operators. The Atangana-Baleanu fractional derivative in Caputo sense (ABC) is given by The AB fractional derivative of order α in Caputo sense is given by [11] where B(α) is a normalization function and B(0) = B(1) = 1 and 0 < α < 1.
The AB fractional integral operator of order α is given by [11] Definition 4. The Atangana-Baleanu fractional integral of order α is defined as [11]

Description of the Modified Method
This section is devoted to descrive the invariant subspace method. Such method has been firstly used in [17] to construct particular exact solutions for partial differential equations of the form where u = u(x,t), u ix = ∂ i u ∂ x i is the i th order derivative of u with respect to the space variable x and F is a nonlinear differential operator. Recently, Gazizov and Kasatkin [22] showed that the invariant subspace method can be applied also to equations with time fractional derivative. In fact, consider the time fractional partial differential equation of the form where F[u] = F(u, u 1x , u 2x , ..., u kx ) and D α t is the time fractional derivative. The modified invariant subspace method is based on the following basic definitions and results [22].
.., f n (x) be an n linearly independent functions and W n is the n-dimensional linear space namely W n = f 1 (x), ..., f n (x) . W n is said to be invariant under the given operator F if F[u] ∈ W n whenever u ∈ W n .
is a solution of equation (8) if and only if the expansion coeffcients u i (t) satisfy the following system of fractional ordinary differential equations Remark 1. The important question concerning the modified invariant subspace method was how to obtain the corresponding invariant subspace of a given differential operator. The answer of this question is given by the following proposition, for more details we refer the reader to [22].

Fractional partial differential equations with Caputo-Fabrizio derivative
In this section, we apply the modified invariant subspace method to construct exact solutions of some partial differential equations with Caputo-Fabrizio derivative in time.
Setting F[u] := CF D α t u(x,t), it is obvious that Eq.(12) admits the following invariant subspace Since Therefore, the exact solution of Eq.(12) can be written as where c 1 (t) and c 2 (t) satisfy the following system of FDEs From the second equation of (16), we find that the function c 2 (t) is a constant and then we assume that c 2 (t) = 1. Thus, the first equation of (16) has the following solution Therefore, Eq. (12) has an exact solution of the form • Example 2: Consider now, the following time-fractional partial differential equation where t > 0, x ∈ R and 0 < α < 1.
It is easy to check that is an invariant subspace of Eq (19), seeing that consequently, an exact solution of Eq. (19) can be written as where c 1 (t) end c 2 (t) are unknown functions to be determined. Substituting Eq. (22) in Eq. (19) yields: The second equation of (23) shows that the function c 2 (t) is a constant and we infer that c 2 (t) = 1 2 . AAccordingly, the first equation of (23) can be expressed as Finally, we obtain an exact solution of Eq. (19) as (25) Fig. 2 Profile of the solution (25) for α = 0.9.
• Example 3: Now we deal with the nonlinear time-fractional partial differential equation where t > 0, x ∈ R and 0 < α < 1. Eq.(26) admits an invariant subspace defined through as As a deduction, an exact solution of Eq.(26) can take the form A modified invariant subspace method for solving partial differential equations with non-singular kernel fractional derivatives 41 Substituting Eq.(29) in Eq.(26) and equating coefficients of different powers of x, we get Solving second equation of (30) gives c 2 (t) = 1.Therefore, the solution of the first equation of (30) is given by: it then follows that an exact solution of Eq.(26) is given by (32) Fig. 3 Profile of the solution (32) for α = 0.9.
• Example 4: Let us consider the following equation where t > 0, x ∈ R, α ∈ [0, 1] and α = 1 2 . It is clear that the above equation Eq.(33) admits an invariant subspace by cause of In an analogous way, the exact solution of Eq.(33) has the form where c 1 (t) and c 2 (t) satisfy the following system of FDEs Similarly, we find that the function c 2 (t) is a constant and then we assume that c 2 (t) = 1. Therefore, the first equation of (37) has the following solution: This is leads eventually to an exact solution to the system Eq.(33) as:

Fractional partial differential equations with Atangana-Baleanu derivative
In what follows, we discuss four examples of getting exact solutions to some partial differential equations with Atangana-Baleanu fractional derivative.
• Example 1: Consider the time-fractional partial differential equation: where t > 0, x ∈ R and 0 < α < 1. which admits an invariant subspace defined through by virtue of It then follows that the form of exact solution for Eq.(40) is Substituting Eq.(42) in Eq.(40) and equating different powers of x to zero yields From the second equation of (43), we find that the function c 2 (t) is a constant, then we assume that c 2 (t) = 1. WWe then conclude that the solution of the first equation of (43) is expressed as Consequently, the exact solution of Eq.(40) reads • Example 2: Let us second consider the following time-fractional partial differential equation where t > 0, x ∈ R, 0 < α < 1.
It is easy to check that Eq.(46) admits an invariant subspace as Since Hence, an exact solution of Eq(46) has the following form In a similar way, substitution of Eq.(49) in Eq. (46) gives From second equation of (50) it comes c 2 (t) is a constant then we assume that c 2 (t) = 1 2 . Therefore, the solution of the first equation of (50) is Finally, we obtain an exact solution of Eq.(46) as (52) Fig. 6 Profile of the solution (52) for α = 0.9.
• Example 3: Consider now the partial differential equation where t > 0, x ∈ R and 0 < α < 1. Equation (26) admits an invariant subspace of the form as far as Then we can form an exact solution of Eq.(26) as where c 1 (t) and c 2 (t) satisfy the following system of FDEs Solving second equation of (30), we get c 2 (t) = 1. Therefore, the solution of the first equation of (57) is constructed as We finally obtain an exact solution of Eq.(53) as • Example 4: WWe finally consider the nonlinear time-fractional partial differential equation where t > 0, x ∈ R, α ∈ [0, 1] and α = 1 2 . It is easy to check that the above Eq.(60) admits an invariant subspace as Since F[c 1 (t) + c 2 (t)x] = c 2 (t)t ∈ W 4 .
The functions c 1 (t) and c 2 (t) satisfy the following system of FDEs CF D 2α t c 1 (t) = c 2 (t)t, CF D 2α t c 2 (t) = 0.
From (64), it can be infered that c 2 (t) is a constant, let us assume that c 2 (t) = 1. The first equation of (64) has then the following solution Accordingly, we get an exact solution of Eq.(60) as

Conclusion
The modifed invariant subspace method was used to seek exact solutions to a class of nonlinear equations with fractional derivatives having nonsingular kernels. Several examples illustrated the effectiveness of the invariant subspace theory for exploring solution of various structures. It is also worth mentionning that the present method does not need any use of laplace transform. Furthermore, some graphical reprensentations are given to show the profiles of the obtained solutions. We stress here that those solutions are very useful to test the efficiency of newly suggested numerical methods for solving partial differential equations with Caputo-Fabrizio or Atangana-Baleanu fractional derivatives.
T h i s p a g e i s i n t e n t i o n a l l y l e f t b l a n k