Daftardar-Jafari Method for Fractional Heat-Like and Wave-Like Equations with Variable Coefficients

In this paper, Daftardar-Gejji and Jafari method is applied to solve fractional heat-like and wave-like models with variable coefficients. The method is proved for a variety of problems in one, two and three dimensional spaces where analytical approximate solutions are obtained. The examples are presented to show the efficiency and simplicity of this method.


Introduction
Various phenomena in engineering physics, chemistry, other sciences can be described very successfully by models using mathematical tools from fractional calculus, i.e., the theory of derivatives and integrals of fractional non-integer order [1] [2] [3] [4].A survey of some applications of fractional calculus in continuum and statistical mechanics is given by Mainardi [5].Oldham and Spanier [6], Miller and Ross [7] provide the history and a comprehensive treatment of this subject.
The fractional calculus has a three-centuary and two decades long history.The idea appeared in a letter by Leibniz to L'Hôpital in 1695.The subject of fractional calculus has gained importance during the past three decades due mainly to its demonstrated applications in different area of physics and engineering.Indeed it provides several potentially useful tools for solving differential and integral equations.It is important to solve time fractional partial differential equations.It was found that fractional time derivatives arise generally as infinitesimal gene-rators of the time evolution when taking a long-time scaling limit.Hence, the importance of investigating fractional equations arises from the necessity to sharpen the concepts of equilibrium, stability states, and time evolution in the longtime limit (for more details, see [8] [9] [10]).
In this work, we will consider the fractional heat-like and wave-like equations of the form xx yy zz u f x y z u g x y z u h x y z u t x a y b z c t subject to the Neumann boundary conditions

y z t f y z t u a y z t f y z t u x z t g x z t u x b z t g x z t u x y t h x y t u x y c t h x y t
and the initial conditions where α is a parameter describing the fractional derivative and t u is the rate of change of temperature at a point over time, ( ) , , h x y z are any functions in x, y, and z.
In the case of 0 1 α < ≤ , then Equation (1) reduces to a fractional heat-like equation with variable coefficients.And in the case of 1 2 α < ≤ , then Equation Recently, Molliq et al. [15] used the Variational iteration method (VIM) for solving Equation (1).The same equation was solved by Momani [16] utilizing the Adomian decomposition method (ADM).Xu and Cang [17], and Xu et al. [18] applied the homotopy analysis method (HAM) to solve it.
The Daftardar-Jafari method (DJM) developed in 2006 has been extensively used by many researchers for the treatment of linear and nonlinear ordinary and partial differential equations of integer and fractional order [19] [20] [21] [22] [23].The method converges to the exact solution if it exists through successive approximations.For concrete problems, a few number of approximations can be used for numerical purposes with high degree of accuracy.The DJM does not require any restrictive assumptions for nonlinear terms as required by some existing techniques.
The main objective of this paper is to apply DJM to obtain fractional solutions for different models of Equation (1).While the VIM [15] requires the determination of Lagrange multiplier in its computational algorithm, DJM is independent of any such requirements.Moreover, unlike the ADM [16], where the cal-culation of the tedious Adomian polynomials is needed to deal with nonlienar terms, DJM handles linear and nonlinear terms in a simple and straightforward manner without any additional requirements.

Fractional Calculus
In this section we will introduce some definitions and properties of the fractional calculus to enable us to follow the solutions of the problem given in this paper.
≤ ≤ , is called the Riemann-Liouville fractional integral operator of or- der α .Properties of the operator J α can be found in [3] [4], we mention the fol- lowing: For The Riemann-Liouville derivative has certain disadvantages when trying to model real-world phenomena with fractional differential equations.Therefore, we shall introduce now a modified fractional differential operator D α proposed by M. Caputo in his work on the theory of viscoelasticity [24].For more information about Caputo definition and properties see [3] [5] [24].Definition 2. The fractional derivative of ( ) The Caputo fractional derivative is considered here because it allows traditional initial and boundary conditions to be included in the formulation of the problem.In this paper, we consider the three-dimensional time fractional heatlike and wave-like Equation (1.1), where the unknown function

( )
, , , u u x y z t = is a assumed to be a causal function of time, i.e., vanishing for 0 t < , and the fractional derivative is taken in Caputo sense to be: Definition 3.For m to be the smallest integer that exceeds α , the Caputo fractional derivative of order 0 α > is defined as For mathematical properties of fractional derivatives and integrals one can consult the above mentioned references.

The DJ Method
Consider the following general functional equation where N is a nonlinear operator from a Banach space B B → and f is a known function, , , , n x x x x =  .We are looking for a solution u of Equation (4) having the series form ( ) ( ) The nonlinear operator N can be decomposed as ( ) From Equations ( 5) and ( 6), Equation ( 4) is equivalent to ( ) We define the recurrence relation ( ) Then ( ) .
The n -term approximate solution of Equation ( 4) is given by In what follows, we apply DJM to six physical models to demonstrate the strength of the method and to establish exact solutions of these models.

Fractional Heat-Like Equations
In this section, we illustrate our analysis by examining the following three fractional heat-like equations.
Example 1. Firstly, let us consider the one-dimensional initial boundary value problems (IBVP) subject to the boundary conditions and the initial condition Operating with m = and using the initial condition (13), we set ( ) ( ) ( ) ( ) Following the algorithm (8), the successive approximations are , . 1 Thus, the approximate solution in a series form is given by So, the solution for the standard heat-like equation ( 1 α = ) is given by ( ) and the initial condition Operating with ( ) ( ) ⋅ on both sides of Equation ( Thus, the approximate solution in a series form is given by ( ) So, the solution for the standard heat-like equation ( 1 This series has the closed form and the initial condition ( ) , , , 0 0.


Thus, the approximate solution in a series form is given by ( So, the solution for the standard heat-like equation ( 1 , , , This series has the closed form

Fractional Wave-Like Equations
In this section, we illustrate our analysis by examining the following three fractional wave-like equations.
Example 1.We first consider the one-dimensional IBVP subject to the boundary conditions ( ) ( ) and the initial conditions , 0 , , 0 .
Following the algorithm (8), the successive approximations are Thus, the approximate solution in a series form is given by So, the solution for the standard heat-like equation ( 2 α = ) is given by ( ) ( ) and the initial conditions , , 0 , , , 0 .


Thus, the approximate solution in a series form is given by This series has the closed form n ∈  .Thus, the approximate solution in a series form is given by ( ) ( )

Conclusion
In this work, DJM has been successfully used to solve fractional heat-like and wave-like equations with variable coefficients giving it a wider applicability.The proposed scheme was applied directly without any need for transformation formulae or restrictive assumptions.Results have shown that the analytical approximate solution process of DJM is compatible with those methods in the literature providing analytical approximation such as VIM, ADM and HAM.The results obtained in all studied cases demonstrate the reliability and the efficiency of this method.

∑
Using the algorithm (8), the successive approximations are

1 m
= and using the initial condition(19), we set

∑
Utilizing the algorithm (8), the successive approximations are exact solution of the problem (17)-(19) compatible with VIM, ADM and HAM.
2 m = and using the initial conditions(22), we set

(
on both sides of Equation (23) yields is the exact solution of the problem (23)-(25) compatible with VIM, ADM and HAM.Example 3. Finally, we consider the three-dimensional inhomogeneous IBVP ,