An Efﬁcient Mechanism to Solve Fractional Differential Equations Using Fractional Decomposition Method

: We present some new results that deal with the fractional decomposition method (FDM). This method is suitable to handle fractional calculus applications. We also explore exact and approximate solutions to fractional differential equations. The Caputo derivative is used because it allows traditional initial and boundary conditions to be included in the formulation of the problem. This is of great signiﬁcance for large-scale problems. The study outlines the signiﬁcant features of the FDM. The relation between the natural transform and Laplace transform is a symmetrical one. Our work can be considered as an alternative to existing techniques, and will have wide applications in science and engineering ﬁelds.


Introduction
For many years the subject of fractional calculus has been studied by many research scholars. This is an ongoing process and one can recognizes that within the field for this study of fractional calculus new techniques and mechanisms show up, which in turn make it possible, to find important challenging insights and unknown correlations between many areas of physics. Recently, an increase in the interest of scientist in nonlocal field theories. In order for us to deal with problems in high-energy and particle physics which only up to this point can be done with local field theories, there is a valid reason for these late developments. Fractional derivatives have proven their capability to describe several phenomena associated with memory and after effects due to their nonlocality property [1,2]. Such phenomena are commonplace in physical processes, biological structures, and cosmological problems. For example, the fractional rheological models have been employed to test the low applied force frequencies [3][4][5][6]. For this reason, it became necessary to illuminate the solutions of the models that describe these phenomena. Several analytical techniques are presented to achieve their objectives. Actually, all these approaches are accommodation for the existing methods to handle the integer case models, which is natural since the fractional derivative generalizes the classical derivative to an arbitrary order.
Recently, fractional calculus and their applications have been treated by many researchers, see [5][6][7][8][9][10][11][12][13] and the references therein. Even though fractional derivatives have existed as long as their integer order counterparts, only in recent decades have fractional derivative models become exciting new tools in the study of practical problems in disciplines as diverse as physics [14][15][16], finance [14,17], biology [6,7] and hydrology [4,10,18,19]. As fractional derivative models are becoming increasingly popular among the wider scientific community, it becomes the main motivation to study numerical schemes for fractional differential equations. Lately, many techniques have discussed the means of exploring approximate solutions of fractional differential equations, such as finite difference and finite element method (FDFEM) [20], spectral method [21] the fractional sub-equation method [10], the fast and parallel computational methods (FPCM) [15,16], the fractional variation method (FVM) [5], the fractional complex transform (FCT) [22], fractional Laplace method [23] and the fractional Adomian decomposition method (FADM) [19,24]. We will use fractional natural decomposition method (FNDM) to solve FODEs and fractional system of ODEs. The method was proven to converge rapidly, and the existence and uniqueness of the method was proven in [25].
Because of the nonstop needs for another new schemes that exist in the literature, e.g. the Mellin transform, and its relation to transforms namely the Sumudu transform and the Al Zaki transform, we developed a valid mechanism which we called the fractional decomposition method (FDM) to find solutions to fractional nonlinear ODEs and PDEs including systems of FDEs. Our scheme avoids the computational difficulties we usually face in other methods exist in the literature, which usually use linearization, discretization and perturbation schemes. Moreover, we have given detailed proofs to some theorems using the proposed method FDM. Also, we presented some important examples of systems fractional differential equations and showed that the FDM provide solutions, which coincide with the ones done by other researchers. To name one, is the solution of diffusion equation with fractional derivatives. These are appeared in Section 4.
Our work will be presented in the following manner: First, in Section 2, we give the history of the natural transform method, and definitions of fractional derivatives. In Section 3, we present proofs to results related to the natural fractional derivatives. Section 4 is devoted to the applications model of FDEs using the proposed method. In Section 5, we solve fractional systems of ODEs. Finally, our concluding remarks are presented, in Acknowledgment, where we outline what we have accomplished in this research.

Related Materials
We explore some definitions and terminologies of the natural transform that will be needed later on in the proofs of our results, (see, for example, [14,17,26]).

Definition 1.
Suppose h(x) is a real function, where x > 0. Then h(x) is said to be in the space C µ , where µ is areal number, if ∃q ∈ R with q > µ, such that h(x) = x q g(x), where g(x) ∈ C[0, ∞), and it is said to be in the space C m µ if h (m) ∈ C µ , m ∈ N.
The fractional derivative of f in the Caputo sense can be defined as Definition 3. Refs. [27,28] The two-parameter Mittag-Leffler function is given by Next, we define the natural transform (N-transform) along with the definition of exponential order function following [23].

Definition 4.
Let h(t) be a function piecewise continuous, with bounded variation, locally integrable (i.e., the function is absolutely integrable in any real interval [a, b] so that b a |h(t)|dt < ∞), and of exponential order, then there exists the N-transform of h(t). That is, a function of exponential order is the one that does not "grow faster" than given exponentials, as t −→ ±∞. Alternatively, there are real constants K, a > 0 such that |h(t)| < K · e at , when t is large and negative (say, for t < t 1 ∈ R). In addition, there are real constants M, b > 0 such that |h(t)| < M · e bt , when t is large (say, for t > t 1 ∈ R). It also has to be true that b < a.
Then, we define the natural transform where N is the N-transform of h(t) and s and u are the N-transform variables. Note that one can writes Equation (2) as N (h(t)) = R − (s, u) + R + (s, u). Hence, Remark 1. Under these conditions, when R(s, u) is analytic, the integral in Equation (2) converges uniformly and absolutely in the complex plane defined by the strip b < Re(s) < a.
Note that one can obtain the Laplace transform and the Sumudu transform if we plug in u = 1, and s = 1 in the above equations, respectively. Hence, the relation between the natural transform and Laplace transform is a symmetry one.
We shall use the well-known gamma function throughout this paper where Γ(z + 1) = zΓ(z).

Important Properties
Here are some of the main properties for the N-transforms which can be found in [26]. , υ > −1, Re(s) > 0.

Natural Caputo Fractional Derivatives
Here, we give detailed proofs to some theorems of N-transform of Caputo fractional derivative. The proof of theorem 1 was given in another published paper by the first author.

Caputo Fractional Derivative
For the sake of readers, we give just some of the natural transform properties. For more properties, we direct the reader to see, for example, [15,16,26].

Proof. First note that
One can conclude from Equation (5), Theorem 5. The natural transform of the Caputo derivative for ϕ(t) = e at is N + c D υ t e a t = a n u n s υ−n u υ (s − au) .
Proof. First note that

Applications of FDM for Fractional ODEs and PDEs
Here we employ the new scheme to solve two non-linear fractional ODEs and we present solution to the diffusion fractional differential equations. Finally, we develop numerical tables for these examples for multiple values of υ and t.

Consider the general nonlinear (FODE)
where t > 0 and 0 < υ ≤ 1, and along with initial condition where c D υ t y(t) is the Caputo derivative for y(t), L is the linear differential operator and F represents the nonlinear part. Also g(t) is the non-homogeneous part, and y 0 is defined and continuous.
We can conclude by applying Theorem 1 to Equation (6) Substituting Equation (7) into Equation (8) and then taking N −1 , one can conclude where H(t) is arising from the non-homogeneous part together with the initial condition. Assuming a solution exist of y(t) as follows But the nonlinear part in our problem is where the Adomian polynomials of y 0 , y 1 , . . . y n are computed using Thus, Equation (12) can be reduced to The rest of Adomian polynomials can be obtained likewise. Substituting Equation (11) into Equation (9), one arrives at One can arrive, with the help of Equation (14), at Eventually, we have the general recursion formula as y n+1 (t) = −N −1 u υ s υ N + (L (y n (t)) + A n (t)) , n ≥ 1.
Hence, our approximate solution is of the form Test Problems Example 1. Consider the nonlinear FDE in the form together with condition y(0) = 0.
Applying Theorem 1 to Equation (17), one arrives at Taking the N-transform inverse of Equation (19) one concludes Suppose a solution exist for ϕ(t) and the nonlinear term ϕ 2 (t) is given as Note here, Using Equation (21), one notices Equation (20) becomes Looking at both sides of Equation (22), one can conclude . Likewise, .
Thus, the approximate solution for ϕ(t) becomes If we choose υ = 1, then Equation (23) becomes This exact solution agrees with the one appearing in the literature (see Figure 1 and Table 1).
Employ first Theorem 1 to Equation (24) to see that Taking the N-transform inverse of Equation (26), we arrive at Suppose that our solution for φ(x, t) is One sees from Equations (27) and (28), that Equation (26) becomes Looking at Equation (29), one can get We follow this direction to obtain Finally, with the help of Equation (29), one can easily explore the rest of the iterations Likewise, Our approximate solution for φ(x, t) is now of the form Therefore, our exact solution is (24) which exists through out the literature (see Table 2).  Figures 2-5 show that the solution peak is high, and one can see that the peak of the solutions of the diffusion equation becomes more and more smooth as the fractional factor υ increases.

Fractional Systems of Ordinary Differential Equations
Next, let us examine two models of systems of FDEs. Then, we present numerical values in tables for some values of t. We only used 5th order approximate solutions for the two functions.  Table 3) together with value conditions x(0) = 2, y(0) = 1.
It follows that, Note that when υ = η = 1, then the exact solutions are

Discussion and Conclusions
Many techniques were used, prior to this work, to handle FDEs. In this work, we have implemented an efficient integral transform method called the Fractional Decom-position Method (FDM) and given detailed proofs to some theorems using the proposed method. We have successfully employed the FDM to obtain exact solutions for the diffusion fractional differential equation and analytical solutions for nonlinear fractional ODE including approximate solutions to two systems of fractional ODE. The FDM reduces the computational difficulties currently exist in the literature. The FDM has self-efficient properties, which help in solving some application with fractional derivatives problems. Our numerical results show that the FDM is a valid and easy scheme for obtaining exact and approximate solutions to fractional differential equation problems. It is proven that the FDM can be implemented without discretization, linearization and perturbation methods used in solving fractional differential equation applications. Thus, it is considered as another option to existing schemes and can be used for wide applications. The FDM has been implemented to many linear and nonlinear FDEs and we have not faced any difficulties when using the new mechanism. Our aim, in the future, is to explore more properties to help solve applications of the FDM, test reasonable series, which converge very fast, and employ it to other fractional integral equation applications.
Finally, the numerical results showed that the new scheme is accurate and efficient. We were able to explore solutions to physical models when υ = η = 1. The next step for our research is to further apply the new scheme to other FDEs that arise in other areas of scientific fields. Earlier works [4,6,14,21] have developed efficient techniques but specialized only for usage on solving specific type of problems, but our FNDM expands their applicability since it is so general.