Approximate Analytical Solutions of Fractional Perturbed Diffusion Equation by Reduced Differential Transform Method and the Homotopy Perturbation Method

The approximate analytical solutions of differential equations with fractional time derivative are obtained with the help of a general framework of the reduced differential transform method (RDTM) and the homotopy perturbation method (HPM). RDTM technique does not require any discretization, linearization, or small perturbations and therefore it reduces significantly the numerical computation. Comparing the methodology (RDTM) with some known technique (HPM) shows that the present approach is effective and powerful. The numerical calculations are carried out when the initial conditions in the form of periodic functions and the results are depicted through graphs.The two different cases have studied and proved that themethod is extremely effective due to its simplistic approach and performance.


Introduction
Various fields of science and engineering deal with dynamical systems (see [1,2]), some of which can be described by fractional-order equations (see [3][4][5] and the references therein).The last two decades have witnessed a great progress in fractional calculus and fractional-order dynamical systems.It has been found that fractional calculus is a mathematical tool that works adequately for anomalous social and physical systems with nonlocal, frequency-and historydependent properties, and for intermediate states such as soft materials, which are neither idea solid nor idea fluid (see [3,4,6]).Differential equations with fractional-order derivatives/integrals are called fractional differential equations, and they have found many successful applications in viscoelasticity, heat conduction, electromagnetic wave, diffusion wave, control theory, and so on (see [7][8][9] and the references therein).
Bistable systems play an important role in the study of spatial patterns.A typical example, which appears in population dynamics, leads to referred to as the classical Fisher-Kolmogorov (FK) equation.
In investigating such spatial patterns, a key role is played by a model equation, which is simpler than the full equation describing the process.Recently, interest has turned to higher-order model equations involving bistable dynamics, such as the extended Fisher-Kolmogorov (EFK) equation proposed by Coullet et al. [10] in 1987 as well as by Dee and Van Saarloos [11] in 1988.The EFK equation has appeared in studies of phase transitions, for instance, near a Lifshitz point (cf.Zimmerman [12]).Another well-known equation of this type is the Swift-Hohenberg equation proposed in 1977.
Whereas the EFK equation and the SH equation involve fourth-order spatial derivatives, certain phase-field models lead to even-higher-order spatial gradients.The following equation studied by Gardner and Jones [13] and Caginalp and Fife [14], in which  and  are constants.In this paper, we consider the initial value problem (IVP) of perturbed diffusion equation with fractional time derivative where the fractional derivative in ( 5) is in the sense of Caputo and 0 ≤  < 1.
Historically, comparatively little was known about the extraordinary range of behavior presented by the solutions of nonlinear partial differential equations.Except for a few special cases, it is impossible to find a closed form solution for a fractional differential equation.Effective methods for solving such equations are needed.So approximate and numerical techniques must be used.
The homotopy perturbation method (HPM) is relatively new approach to provide an analytical approximation to linear or nonlinear problem.This method was first presented by He [15,16] and applied to various nonlinear problems [17,18].The basic difference of this method from the other perturbation techniques is that it does not require small parameters in the equation, which overcomes the limitations of the traditional perturbation techniques.Recently, the application of this method is extended for fractional differential equations [5,[19][20][21].Zhou [22] was the first one to use differential transform method (DTM) in engineering applications.He employed DTM in solution of initial boundary value problems in electric circuit analysis.This method constructs an analytical solution in the form of a polynomial.It is different from the traditional high order Taylor series method, which requires symbolic computation of necessary derivatives of the data functions.The differential transform is an iterative procedure for obtaining analytical differential equations.The concept of DTM has broadened to problems involving partial differential equations and systems of differential equations [23,24].Some researchers have lately applied DTM for analysis of uniform and non-uniform beams [25][26][27].Recently, the application of reduced differential transform method is successfully extended to obtain analytical approximate solutions to linear and nonlinear ordinary differential equations of fractional order [28][29][30].
In this paper, we try to find an approximate analytical solution of (5) for fractional time derivatives with the help of powerful analytical method.We use the reduced differential transform method (RDTM) and homotopy perturbation method (HPM) to obtain the solutions and compare them with each other.We know that the HPM method is based on the use of homotopy parameter for classification of most favorable values of parameters in between [0, 1].While, RDTM technique does not require any parameter, discretization, linearization, or small perturbations and therefore it reduces significantly the numerical computation.For the standard cases, comparing the methodology with some known techniques shows that the present approach is effective and powerful.

Preliminaries
There are several approaches to define the fractional calculus, for example, Riemann-Liouville, Gruunwald-Letnikow, Caputo, and Generalized Functions approach.For the readers' convenience, definitions of fractional integral/derivative and some preliminary results are given in this section.
Definition 1 (see [6]).The fractional integral of order  ≥ 0 of a function (, ) is given by provided that the right side is point-wise defined on (0, +∞), where Γ(⋅) is the well-known gamma function.
Riemann-Liouville fractional derivative is mostly used by mathematicians but this approach is not suitable for real world physical problems since it requires the definition of fractional order initial conditions, which have no physically meaningful explanation yet.Caputo introduced an alternative definition, which has the advantage of defining integer order initial conditions for fractional order differential equations.We have chosen the Caputo fractional derivative because it allows traditional initial and boundary conditions to be included in the formulation of the problem.And some other properties of fractional derivative can be found in [4,6].
The relation between the Riemann-Liouville operator and Caputo operator is given by By using (9), from (8), we obtain fractional derivative in the Caputo sense as follows: Let us expand the analytical and continuous function (, ) in terms of a fractional power series as follows: where  is the order of fraction and   () is the fractional differential transform of (, ).
Since the initial conditions are implemented to the integer order derivatives, the transformations of the initial conditions are defined as follows: for  = 0, 1, 2, . . ., ( − 1) The following theorems that can be deduced from ( 7) and ( 8) are given below.

Solution of the Problem by the RDTM and HPM
To illustrate the basic ideas of HPM [15], we consider the following general nonlinear differential equation: with boundary conditions where  is a general differential operator,  is a boundary operator, () is a known analytic function, and Γ is the boundary of the domain Ω.
The operator  can be generally divided into two parts  and , where  is linear, while  is nonlinear.Therefore ( 14) can be written as follows: We construct a homotopy of ( 14) (, ) : Ω × [0, 1] → R which satisfies which is equivalent to  (, ) =  () −  ( 0 ) +  ( 0 ) +  [ () −  ()] = 0, where  ∈ [0,1] is an embedding parameter, and  0 is an initial guess approximation of ( 14) which satisfies the boundary conditions.It follows from ( 17) and ( 18) that Thus, the changing process of  from 0 to 1 is just that of (, ) from  0 () to ().In topology this is called deformation and () − ( 0 ) and () − () are called homotopic.Here the embedding parameter is introduced much more naturally, unaffected by artificial factors; further it can be considered as a small parameter for 0 ≤  ≤ 1.So it is very natural to assume that the solution of ( 18) and ( 19) can be expressed as The approximate solutions of the original equations can be obtained by setting  = 1; that is, The convergence of series (21) has been proved by He in his paper [15].

Mathematical Problems in Engineering
Because of the knowledge of various perturbation methods that low-order approximate solution leads to high accuracy, there requires no infinite series.Then after a series of recurrent calculation by using Mathematica software, we will get approximate solutions of fractional model (5).
According to HPM, we construct the following homotopy of ( 5) where the homotopy parameter  is considered to be small, 0 ≤  ≤ 1.
Assuming the solution of ( 22) has the form when  → 1, (23) becomes the approximate solution of (5).Substituting ( 23) into ( 22) and equating the terms with identical powers of , we obtain the following set of linear differential equations: and so on.The method is based on applying the operators    (the inverse operators of the Caputo derivative    ) on both sides of the above linear differential equations.Using this selection we obtain the successive approximations and the solution may be obtained.
Finally, we approximate the analytical solutions of (, ) by the truncated series For the illustration of the methodology of the reduced differential transform method, we write (5) in the standard operator form with initial condition where ((, )) = (  /  )(, ) is the fractional time derivative operator and ((, )) = From initial condition, we write  0 () =  () .

Numerical Results and Discussion
Here, we take three different values of () and comparing the results of RDTM and HPM in the form of three-and two dimensional figures for each case, we would see that RDTM and HPM solutions are in excellent agreement.In this section we assume that  =  = 1.

Conclusions
In this paper, RDTM and HPM are successfully applied to find the solution of the fractional differential equation and standard motion with different initial conditions.Unlike the traditional methods, the solutions here are given in series 9 form.The approximate solution to the equation was computed without any need for special transformations, linearization, or discretization.In addition, we compare these two methods and show that the results of the RDTM method are in excellent agreement with results of the HPM method and the obtained numerical solutions are shown graphically.We use the mathematical simulation to calculate the functions obtained from the RDTM and HPM.It was shown that RDTM and HPM methods are powerful tools for solving analytically of nonlinear equations.0.05 0 0 0.0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5