Daftardar-Jafari Method for Fractional Heat-Like and Wave-Like Equations with Variable Coefficients ()
1. 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 generators 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 long- time 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
(1)
subject to the Neumann boundary conditions
(2)
and the initial conditions
(3)
where
is a parameter describing the fractional derivative and
is the rate of change of temperature at a point over time,
is temperature as a function of time and space, while
and
are the second spatial derivatives (thermal conductions) of temperature in x, y, and z directions, respectively. Finally,
and
are any functions in x, y, and z.
In the case of
, then Equation (1) reduces to a fractional heat-like equation with variable coefficients. And in the case of
, then Equation (1) reduces to a fractional wave-like equation which models anomalous diffusive and subdiffusive systems, description of fractional random walk, unification of diffusion and wave propagation phenomena [11] [12] [13] [14] .
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 calculation 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.
2. 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. These definitions include, Riemann-Liouville, Weyl, Reize, Compos, Caputo, and Nashimoto fractional operators.
Definition 1. Let
The operator
defined on the usual Lebesque space
by
for
, is called the Riemann-Liouville fractional integral operator of order
.
Properties of the operator
can be found in [3] [4] , we mention the following:
For
and
1.
exists for almost every
2.
3.
4.
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
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
in the Caputo sense is defined as
for
Also, we need here two of its basic properties.
Lemma 1. If
and
then
and
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 heat- like and wave-like Equation (1.1), where the unknown function
is a assumed to be a causal function of time, i.e., vanishing for
, and the fractional derivative is taken in Caputo sense to be:
Definition 3. For
to be the smallest integer that exceeds
, the Caputo fractional derivative of order
is defined as
For mathematical properties of fractional derivatives and integrals one can consult the above mentioned references.
3. The DJ Method
Consider the following general functional equation
(4)
where
is a nonlinear operator from a Banach space
and
is a known function,
. We are looking for a solution
of Equation (4) having the series form
(5)
The nonlinear operator
can be decomposed as
(6)
From Equations (5) and (6), Equation (4) is equivalent to
(7)
We define the recurrence relation
(8)
Then
(9)
and
(10)
The
-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.
4. Fractional Heat-Like Equations
In this section, we illustrate our analysis by examining the following three frac- tional heat-like equations.
Example 1. Firstly, let us consider the one-dimensional initial boundary value problems (IBVP)
(11)
subject to the boundary conditions
(12)
and the initial condition
(13)
Operating with
on both sides of Equation (11) yields
For
with
and using the initial condition (13), we set
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 (
) is given by
This series has the closed form
which is the exact solution of the problem (11)-(13) compatible with VIM, ADM and HAM.
Example 2. Now, let us consider the two-dimensional IBVP
(14)
subject to the Neumann boundary conditions
(15)
and the initial condition
(16)
Operating with
on both sides of Equation (14) yields
For
with
and using the initial condition (16), we set
Using 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 (
) is given by
This series has the closed form
which is the exact solution of the problem (14)-(16) compatible with VIM, ADM and HAM.
Example 3. Let us consider the three-dimensional inhomogeneous IBVP
(17)
subject to the boundary conditions
(18)
and the initial condition
(19)
Operating with
on both sides of Equation (17) yields
For
with
and using the initial condition (19), we set
Utilizing 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 (
) is given by
This series has the closed form
which is the exact solution of the problem (17)-(19) compatible with VIM, ADM and HAM.
5. 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
(20)
subject to the boundary conditions
(21)
and the initial conditions
(22)
Operating with
on both sides of Equation (20) yields
For
with
and using the initial conditions (22), we set
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 (
) is given by
This series has the closed form
which is the exact solution of the problem (20)-(22) compatible with VIM, ADM and HAM.
Example 2. We next consider the two-dimensional IBVP
(23)
subject to the Neumann boundary conditions
(24)
and the initial conditions
(25)
Operating with
on both sides of Equation (23) yields
For
with
and using the initial conditions (25), we set
Using 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 (
) is given by
This series has the closed form
which 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
(26)
subject to the boundary conditions
(27)
and the initial conditions
(28)
Operating with
on both sides of Equation (26) yields
For
with
and using the initial conditions (28), we set
Utilizing 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 (
) is given by
This series has the closed form
which is the exact solution of the problem (26)-(28) compatible with VIM, ADM and HAM.
6. 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.