# Semianalytic Solution of Space-Time Fractional Diffusion Equation

Authors: A. Elsaid; S. Shamseldeen; S. Madkour
Journal: International Journal of Differential Equations
 (2016)

---

## Abstract

We study the space-time fractional diffusion equation with spatial Riesz-Feller fractional derivative and Caputo fractional time derivative. The continuation of the solution of this fractional equation to the solution of the corresponding integer order equation is proved. The series solution of this problem is obtained via the optimal homotopy analysis method (OHAM). Numerical simulations are presented to validate the method and to show the effect of changing the fractional derivative parameters on the solution behavior.


---

## Body

## 1. Introduction

Fractional derivatives have been utilized for describing numerous models in different branches of applied science. Anomalous diffusion is one process that has been successfully modeled by fractional derivatives. This type of diffusion is characterized by the nonlinear dependence of the mean square displacementx
(
t
) of a diffusing particle over time t: x
2
(
t
)
∝
k
α
t
α, and it can be interpreted as the Lévy stable densities. On the other hand, in the case of classical diffusion, linear dependence x
2
(
t
)
∝
k
t occurs and it follows Gaussian statistics and Fick’s second law for running processes at time t [1].Anomalous diffusion is described by fractional partial differential equations (FPDEs) by replacing classical derivatives with derivatives of fractional order. A type of anomalous diffusion that received great attention is the fractional diffusion equation with spatial Riesz and Riesz-Feller fractional derivatives [2, 3]. These fractional derivative operators are used to model phenomena as anomalous diffusion [1], discrete random walk models [3, 4], and continuous random walk models [5] in the form of FPDEs.General solutions are obtained for some classes of linear FPDEs, using integral transforms. For example, in [6], a general solution is given for a time fractional diffusion-wave equation defined in a bounded space domain. The fractional time derivative is described in the Caputo sense. The finite sine Laplace transform technique is used to obtain the solution. In [7], analytic techniques are introduced to solve certain types of multiterm time-space Caputo-Riesz fractional advection-diffusion equations on a finite domain using the equivalent relationship between the Laplacian operator and the Riesz fractional derivative. In [8], a fundamental solution is investigated for the space-time fractional diffusion-wave equation, which is obtained from the standard diffusion equation by replacing the second-order space derivative with a Riesz-Feller derivative and the first-order time derivative with a Caputo derivative. The paper [9] presents a solution of a unified reaction-diffusion equation with the Riemann-Liouville fractional derivative as the time derivative and Riesz-Feller derivative as the space derivative in terms of H-function. In [10], a closed form solution in terms of Mittag-Leffler function is obtained for a space-time fractional telegraph equation with Riesz-Feller derivative as the space derivative. In [11], a solution of a space-time fractional advection-dispersion equation is derived in terms of the Green functions.The proposed abovementioned methods work in case of linear FPDEs. However, for the nonlinear case, the numerical methods and the semianalytic approaches have the advantage. Several numerical techniques are introduced in the literature, and we state some of the known methods here. In [12], the matrix transform method with the method of lines is used to obtain the numerical solution of a fractional partial differential equation with Riesz space fractional derivatives on a finite domain. An explicit finite-difference approximation is used to solve the variable-order nonlinear space fractional diffusion equation [13]. Crank-Nicolson method for the fractional diffusion equation with the Riesz fractional derivative is used in [14], while the improved matrix transform method with Pade approximation is used in [15]. Yet, semianalytic solutions for this type of problems include only the series solution by the variational iteration method approach (see [16, 17]) and by the homotopy analysis method (HAM) (see, e.g., [2, 18]). In [19], the homotopy analysis method with the Fourier transform is used to obtain and analytically approximate solutions of Riesz fractional diffusion equation and Riesz fractional advection-dispersion equation.As seen in the literature review, only few articles dealt with applying iterative techniques to these types of FPDEs due to the difficulty of repeated application of Riesz or Riesz-Feller fractional derivatives to the solution components. In this work, we consider the space-time fractional diffusion equation of the following form:(1)
D
t
β
C
u
x
,
t
=
k
u
R
x
α
,
θ
u
x
,
t
+
g
u
,
-
∞
<
x
<
∞
,
t
>
0
,
u
x
,
0
=
f
x
,where D
C
t
β denotes the Caputo fractional derivative (in time) of order β and R
x
α
,
θ denotes the Riesz-Feller fractional derivative (in space) of order α and skewness θ. Parameters α
,
θ, and β are real numbers restricted to(2)
0
<
α
≤
2
,
α
≠
1
,
θ
≤
min
⁡
α
,
2
-
α
,
0
<
β
≤
1
, and the two functions k and g are continuous functions in u. We aim to establish the continuation of the solution of (1) to the exact solution of the corresponding equation in Riesz fractional derivative as the skewness parameter approaches zero. This objective is carried out theoretically then via approximate series solution obtained iteratively by applying the OHAM.This paper is organized as follows. In Section2, basic definitions of fractional derivative operators involved are presented. Proof of solution continuation is presented in Section 3. The OHAM is illustrated in Section 4. The results of numerical experiments are presented in Section 5. Section 6 contains the conclusion of this work.


## 2. Fractional Derivatives and Integrals

Definition 1.
A real functionf
(
x
), x
>
0, is said to be in the space C
μ, μ
∈
R, if there exists a real number p
>
μ, such that f
(
x
)
=
x
p
f
1
(
x
), where f
1
(
x
)
∈
C
(
0
,
∞
), and it is said to be in the space C
μ
m if f
m
∈
C
μ, m
∈
N
.Definition 2.
The Riemann-Liouville fractional integral operator of orderα
≥
0 of a function f
(
x
)
∈
C
μ, μ
≥
-
1, is defined as(3)
J
α
f
x
=
1
Γ
α
∫
0
x
x
-
τ
α
-
1
f
τ
d
τ
,
α
>
0
,
x
>
0
,
J
0
f
x
=
f
x
.The operatorJ
α satisfies the following properties. For f
∈
C
μ, μ
≥
-
1, α
,
β
≥
0, and γ
>
-
1,(1)
J
α
J
β
f
(
x
)
=
J
α
+
β
f
(
x
),
(2)
J
α
J
β
f
(
x
)
=
J
β
J
α
f
(
x
),
(3)
J
α
x
γ
=
(
Γ
γ
+
1
/
Γ
γ
+
α
+
1
)
x
α
+
γ.Definition 3.
The fractional derivative in Caputo sense off
(
x
)
∈
C
-
1
m, m
∈
N, x
>
0, is defined as(4)
D
C
x
β
f
x
=
J
m
-
β
d
m
d
x
m
f
x
,
m
-
1
<
β
<
m
,
d
m
d
x
m
f
x
,
β
=
m
.Some basic properties of Caputo fractional derivative are listed in the following lemma.Lemma 4.
Ifm
-
1
<
α
≤
m, m
∈
N, and f
∈
C
μ
m, μ
≥
-
1, then,(5)
D
α
c
=
0
,
c
is  a  constant
D
α
x
γ
=
Γ
γ
+
1
Γ
γ
-
α
+
1
x
α
+
γ
,
γ
>
α
-
1
,
0
,
γ
≤
α
-
1
,
D
α
J
α
f
x
=
f
x
,
J
α
D
α
f
x
=
f
x
-
∑
k
=
0
m
-
1
f
k
0
x
k
k
!
,
x
>
0
.Definition 5.
The Riesz-Feller partial fractional derivativeR
x
α
,
θ is defined as [3](6)
R
x
α
,
θ
u
x
=
-
c
+
D
+
α
u
x
+
c
-
D
-
α
u
x
,
0
<
α
<
2
,
α
≠
1
,where D
±
α
u
(
x
) are the Weyl fractional derivatives defined as(7)
D
±
α
u
x
=
±
d
d
x
W
±
1
-
α
u
x
,
0
<
α
<
1
d
2
d
x
2
W
±
2
-
α
u
x
,
1
<
α
<
2
,and W
±
β denote the Weyl fractional integrals of order β
>
0, given by(8)
W
+
β
u
x
=
1
Γ
β
∫
-
∞
x
x
-
z
β
-
1
u
z
d
z
,
W
-
β
u
x
=
1
Γ
β
∫
x
∞
z
-
x
β
-
1
u
z
d
z
.The coefficients c
+ and c
- have the following forms:(9)
c
±
=
s
i
n
α
∓
θ
π
/
2
s
i
n
α
π
,whereθ
≤
min
⁡
(
α
,
2
-
α
). When α
=
0, the Weyl fractional derivative degenerates into the identity operator:(10)
D
±
0
u
x
=
I
u
x
=
u
x
.For continuity, we get(11)
D
±
1
u
x
=
±
d
d
x
u
x
,
D
±
2
u
x
=
d
2
d
x
2
u
x
.Evidently, in the case of α
=
2 and θ
=
0, we define(12)
R
x
α
,
0
u
x
=
d
2
d
x
2
u
x
.For the case of α
=
1 and θ
=
0, we have(13)
R
x
1,0
u
x
=
d
d
x
H
u
x
=
d
d
x
1
π
∫
-
∞
∞
u
z
z
-
x
d
z
,where H is the Hilbert transform and the integral is understood in the Cauchy principal value sense. For the cases α
=
1 and θ
=
±
1, we obtain the first derivative:(14)
R
x
1
,
±
1
u
x
=
±
d
d
x
u
x
.The special case θ
=
0 is known as Riesz fractional derivative.


## 3. Optimal Homotopy Analysis Method (OHAM)

In this section, we begin by illustrating the basic steps of the classical homotopy analysis method (HAM). Consider the following nonlinear equation:(15)
N
u
x
,
t
=
0
,where N is a nonlinear operator, u
(
x
,
t
) is the unknown function, and x and t denote spatial and temporal independent variables, respectively. Liao [20] generalized the traditional homotopy method to construct the so-called zero-order deformation equation given by(16)
1
-
p
L
ϕ
x
,
t
;
p
-
u
0
x
,
t
=
p
ħ
H
x
,
t
N
ϕ
x
,
t
;
p
,where p
∈
[
0,1
] is an embedding parameter, ħ is a nonzero auxiliary parameter, H
(
x
,
t
) is an auxiliary function, L is an auxiliary linear operator, u
0
(
x
,
t
) is an initial guess of u
(
x
,
t
), and ϕ
(
x
,
t
;
p
) is an unknown function. Obviously, when p
=
0 and p
=
1, we have ϕ
(
x
,
t
;
0
)
=
u
0
(
x
,
t
) and ϕ
(
x
,
t
;
1
)
=
u
(
x
,
t
), respectively. Thus, as p increases from 0 to 1, the solution ϕ
(
x
,
t
;
p
) varies from the initial guess u
0
(
x
,
t
) to the required solution u
(
x
,
t
). By expanding ϕ
(
x
,
t
;
p
) in Taylor series with respect to p, we have(17)
ϕ
x
,
t
;
p
=
u
0
x
,
t
+
∑
m
=
1
∞
u
m
x
,
t
p
m
,where(18)
u
m
x
,
t
=
1
m
!
∂
m
ϕ
x
,
t
;
p
∂
p
m
p
=
0
.If the auxiliary linear operator, the initial guess, and the auxiliary parameter ħ and the auxiliary function are so properly chosen, then, as proved by Liao [20], series (17) converges at p
=
1 and one has(19)
u
x
,
t
=
u
0
x
,
t
+
∑
m
=
1
∞
u
m
x
,
t
,which, as proven by Liao [20], must be one of the solutions of the original nonlinear equation. Using definition (18), the governing equation of the HAM can be deduced from the zero-order deformation equation (16) as follows. Define the vector(20)
u
→
n
=
u
0
x
,
t
,
u
1
x
,
t
,
u
2
x
,
t
,
…
,
u
n
x
,
t
.From (16), the so-called mth-order deformation equation is given by(21)
L
u
m
x
,
t
-
χ
m
u
m
-
1
x
,
t
=
ħ
H
x
,
t
R
m
u
→
m
-
1
x
,
t
,where(22)
R
m
u
→
m
-
1
=
1
m
-
1
!
∂
m
-
1
N
ϕ
x
,
t
;
p
∂
p
m
-
1
p
=
0
,
χ
m
=
0
,
m
≤
1
,
1
,
m
>
1
.Applying the inverse operator L
-
1 to both sides of (21), u
m
(
x
,
t
) can be calculated by symbolic computations software. The HAM has been successfully applied to solve various classes of equations and applied problems [21–25].In the classical HAM, choosing the value of parameterħ depends on inspecting the graph of the quantity of interest: the solution or one of its derivatives. Yet, when H
(
x
,
t
) is fixed, it is obvious that u
m
(
x
,
t
) contains only one control parameter ħ. Thus, by constructing a formula for the residual error, the OHAM solution is obtained by choosing the value for parameter ħ that minimizes the error. Here, the averaged residual error defined for ordinary differential equations in [26] is generalized to the case of two variable partial differential equations in the following form:(23)
E
m
ħ
=
1
M
K
∑
i
=
0
M
∑
j
=
0
K
N
∑
n
=
0
m
u
n
i
M
,
j
K
2
,which is a nonlinear algebraic equation of one unknown: the convergence-control parameter ħ. Thus, the optimal value of ħ is determined by the minimum of the averaged residual error E
m to ensure the fast convergence of the homotopy series.


## 4. Continuation of the Solution

In this section, we define Riesz-Feller fractional derivative on unbounded domain in Caputo sense. This definition is considered as a generalization of the definition given in [27] for Riesz fractional derivative on bounded domains.Definition 6.
The Riesz-Feller fractional derivative in Caputo sense is defined by(24)
R
C
t
α
,
θ
f
t
=
-
c
+
D
C
+
α
+
c
-
D
C
-
α
,
0
<
α
<
2
,
α
≠
1
,where D
C
+
α and D
C
-
α are Weyl fractional derivatives defined in Caputo sense for α
∈
(
0,2
)
,
α
≠
1, as(25)
D
C
±
α
f
t
=
±
W
±
1
-
α
d
d
t
f
t
,
0
<
α
<
1
,
W
±
2
-
α
d
2
d
t
2
f
t
,
1
<
α
<
2
,and defined for α
=
0, 1, and 2 as(26)
D
C
±
0
f
t
=
f
t
,
D
C
±
1
f
t
=
±
d
d
t
f
t
,
D
C
±
2
f
t
=
d
2
d
t
2
f
t
.Lemma 7.
Letf
(
x
) belong to the class of “good functions” [28]. Then, for α
∈
(
0,2
), α
≠
1, (27)
R
C
x
α
,
θ
f
x
=
R
x
α
,
θ
f
x
.Proof.
Consider the caseα
∈
(
0,1
). Riesz-Feller fractional derivative can be written in the following form:(28)
R
x
α
,
θ
f
x
=
-
d
d
x
∫
0
∞
c
+
f
x
-
λ
-
c
-
f
x
+
λ
λ
α
d
λ
,which can be written as(29)
R
x
α
,
θ
f
x
=
∫
0
∞
c
-
f
′
x
+
λ
-
c
+
f
′
x
-
λ
λ
α
d
λ
.For R
C
x
α
,
θ, we have(30)
R
C
x
α
,
θ
f
x
=
c
+
∫
-
∞
x
f
′
z
x
-
z
α
d
z
-
c
-
∫
x
∞
f
′
y
y
-
x
α
d
y
,which by substitution yields(31)
R
C
x
α
,
θ
f
x
=
∫
0
∞
c
-
f
′
x
+
λ
-
c
+
f
′
x
-
λ
λ
α
d
λ
,which is the same as (29). A similar argument holds for the case α
∈
(
1,2
).As the equivalence between the Riemann and Caputo definitions of Riesz-Feller fractional derivative is deduced, we discuss the continuation of the solution obtained.Theorem 8.
Iff
(
x
) is a function in L
1
(
-
∞
,
∞
), then the exact solution u
β
α
,
θ of the space-time fractional diffusion equation (1) is given by(32)
u
β
α
,
θ
x
,
t
=
1
2
π
∬
-
∞
∞
E
β
-
ω
α
t
β
e
-
i
θ
π
/
2
f
v
e
i
ω
x
-
v
d
ω
d
v
,where E
β
(
z
) is the one parameter Mittag-Leffler function defined as [29](33)
E
β
z
=
∑
k
=
0
∞
z
k
Γ
1
+
k
β
,
β
>
0
.Theorem 9.
Letα
∈
(
1,2
), β
∈
(
0,1
), θ
≤
min
⁡
(
α
,
2
-
α
), and f
(
x
) be a function in L
1
(
-
∞
,
∞
), and let u
β
α
,
θ be solution (32) of the space-time fractional problem (1); then,(34)
lim
θ
→
0
⁡
u
β
α
,
θ
x
,
t
=
u
β
α
x
,
t
,where u
β
α
(
x
,
t
) is the exact solution of Riesz fractional diffusion equation:(35)
D
C
t
β
u
x
,
t
=
R
x
α
u
x
,
t
,
-
∞
<
x
<
∞
,
t
>
0
,
u
x
,
0
=
f
x
,with the same limits for parameters α and β.Proof.
Consider the set of functions:(36)
φ
n
ω
=
1
2
π
E
β
-
ω
α
t
β
e
-
i
π
/
2
n
∫
-
∞
∞
f
v
e
i
ω
x
-
v
d
v
,
ω
∈
-
∞
,
∞
,
n
=
1,2
,
…Then,(37)
φ
n
ω
≤
1
2
π
E
β
-
ω
α
t
β
e
-
i
π
/
2
n
∫
-
∞
∞
f
v
d
v
,and as f
(
x
)
∈
L
1
(
-
∞
,
∞
), then there exists M
>
0 such that(38)
φ
n
ω
≤
M
2
π
E
β
-
ω
α
t
β
e
-
i
π
/
2
n
.From [30, Theorem (1.6)], there exists K
1
>
0 such that(39)
E
β
-
z
≤
K
1
1
+
z
.Then,(40)
φ
n
ω
≤
1
2
π
M
K
1
1
+
ω
α
t
β
,
ω
∈
-
∞
,
∞
,
n
=
1,2
,
…
.By Lebesgue dominated convergence theorem and setting θ
=
1
/
n, we have(41)
lim
θ
→
0
⁡
u
β
α
,
θ
x
,
t
=
lim
n
→
∞
⁡
∫
-
∞
∞
φ
n
ω
d
ω
=
∫
-
∞
∞
l
i
m
n
→
∞
φ
n
ω
d
ω
=
1
2
π
∬
-
∞
∞
E
β
-
ω
α
t
β
f
v
e
i
ω
x
-
v
d
ω
d
v
=
u
β
α
x
,
t
,which is the exact solution of the Riesz fractional diffusion equation (35).Similar arguments hold to prove the continuation of the solution when both parametersα and β tend to their integer limits.


## 5. Numerical Simulation

In this section, we consider linear and nonlinear problems to illustrate the efficiency of the OHAM to this type of problems and to illustrate the continuation of the solution we proved in Section4. Yet, to apply the recursive technique of OHAM to problem (1), a repeated evaluation of Riesz-Feller fractional derivative to solution components is needed. This is accomplished by using a property of Riesz-Feller fractional derivative that is proved in the following lemma.Lemma 10.
Letα
∈
(
0,2
), α
≠
1, θ
≤
min
⁡
(
α
,
2
-
α
), and ω
>
0. Then,(42)
R
x
α
,
θ
e
i
ω
x
=
-
ω
α
e
i
ω
x
-
θ
π
/
2
.Proof.
Case
α
∈
(
0,1
). From the definition of Riesz-Feller fractional derivative and by substitution, we can write(43)
R
x
α
,
θ
e
i
ω
x
=
-
d
d
x
∫
0
∞
c
+
e
i
ω
x
-
τ
-
c
-
e
i
ω
x
+
τ
Γ
1
-
α
τ
α
d
τ
,which reduces to(44)
R
x
α
,
θ
e
i
ω
x
=
-
i
ω
Γ
1
-
α
e
i
ω
x
c
+
Γ
1
-
α
i
ω
1
-
α
-
c
-
Γ
1
-
α
-
i
ω
1
-
α
,
R
x
α
,
θ
e
i
ω
x
=
-
ω
α
e
i
ω
x
c
+
i
α
-
c
-
-
i
α
. Substituting for c
+ and c
-, we have(45)
R
x
α
,
θ
e
i
ω
x
=
-
ω
α
e
i
ω
x
sin
⁡
α
-
θ
π
/
2
e
i
π
α
/
2
+
sin
⁡
α
+
θ
π
/
2
e
-
i
π
α
/
2
sin
⁡
α
π
.Utilizing Euler form and trigonometric identities, we have (46)
R
x
α
,
θ
e
i
ω
x
=
-
ω
α
e
i
ω
x
-
θ
π
/
2
.
Case
α
∈
(
1,2
). Consider(47)
R
x
α
,
θ
e
i
ω
x
=
-
d
2
d
x
2
∫
0
∞
c
+
e
i
ω
x
-
τ
-
c
-
e
i
ω
x
+
τ
Γ
2
-
α
τ
α
-
1
d
τ
,which reduces to(48)
R
x
α
,
θ
e
i
ω
x
=
ω
2
Γ
2
-
α
e
i
ω
x
c
+
Γ
2
-
α
i
ω
2
-
α
-
c
-
Γ
2
-
α
-
i
ω
2
-
α
.Carrying out the same procedure, the lemma follows.Example 11.
Consider problem (1) with k
(
u
)
=
1, g
(
u
)
=
a
u, and f
x
=
ϵ
sin
⁡
(
π
x
/
L
):(49)
D
C
t
β
u
x
,
t
=
R
x
α
,
θ
u
x
,
t
+
a
u
,
-
∞
<
x
<
∞
,
t
>
0
,
u
x
,
0
=
ϵ
sin
⁡
π
x
L
,where a, ϵ, and b are real constants.To solve this problem by means of OHAM, let us consider the linear operator(50)
L
ϕ
x
,
t
;
p
=
D
C
t
β
ϕ
x
,
t
;
pand the operator(51)
N
ϕ
x
,
t
;
p
=
D
C
t
β
ϕ
-
R
x
α
,
θ
ϕ
-
a
ϕ
.Then, one can write themth-order deformation equation, with H
(
x
,
t
)
=
1, for this problem as(52)
L
u
m
x
,
t
-
χ
m
u
m
-
1
x
,
t
=
ħ
R
m
u
m
-
1
x
,
t
,where(53)
R
m
u
m
-
1
=
D
C
t
β
u
m
-
1
-
R
x
α
,
θ
u
m
-
1
-
a
u
m
-
1
.Apply the inverse operatorL
-
1
=
J
t
μ to both sides of (52); we obtain(54)
u
m
x
,
t
=
χ
m
u
m
-
1
x
,
t
+
u
m
x
,
0
-
χ
m
u
m
-
1
x
,
0
+
ħ
J
t
β
R
m
u
m
-
1
x
,
t
,where(55)
u
0
x
,
t
=
u
x
,
0
=
ϵ
sin
⁡
π
x
L
,
u
m
x
,
0
=
0
,
m
>
1
.Referring to (53) and (54) with conditions (55), we get(56)
u
1
x
,
t
=
h
t
β
ϵ
Γ
1
+
β
-
a
sin
⁡
π
x
L
+
π
L
α
sin
⁡
π
x
L
-
π
θ
2
,
u
2
=
h
1
+
h
t
β
ϵ
Γ
1
+
2
β
-
a
sin
⁡
π
x
L
+
π
L
α
sin
⁡
π
x
L
-
π
θ
2
+
h
2
t
2
β
ϵ
Γ
1
+
2
β
a
2
sin
⁡
π
x
L
+
π
L
α
π
L
α
sin
⁡
π
x
L
-
π
θ
-
2
a
sin
⁡
π
x
L
-
π
θ
2
.The series solution is given byu
=
u
0
+
u
1
+
u
2
+
u
3
+
⋯. Figure 1 shows the effect of the skewness parameter θ on the behavior of the solution. The plot represents the sum of the first five terms (u
0 to u
4
) in the OHAM series when a
=
1.0, ϵ
=
1.0, L
=
3.0, t
=
0.2, and 0
<
x
<
π and the fractional parameters β
=
1.0, α
=
1.5, and θ
=
0,0.5
,
-
0.5
. At these different values of θ
=
0,0.5, and -
0.5, the optimal convergence-control parameter ħ is calculated as ħ
=
-
0.995787
,
-
1.049435, and -
1.056491, and the corresponding residual error is obtained as E
=
2.969688
∗
10
-
7
,
4.243949
∗
10
-
6, and 4.453205
∗
10
-
6, respectively.Figure 1
The behavior of the four terms of OHAM solution of (49) at t
=
0.2, 0
<
x
<
π, and the fractional parameters α
=
1.5 and β
=
1 and the skewness θ
=
0,0.25
,
-
0.25.The solution (the first five terms in the OHAM series) behavior as the parametersθ, α, and β change at a fixed time t
=
5.0 and in the interval 0
≤
x
≤
π is shown in Figures 1–3. The amplitude of the solution decreases α and/or β increases. As α reaches 2 and/or β increases up to unity, the series solution approximately coincides with the exact solution of the corresponding integer order equation u
=
e
-
(
π
2
/
L
2
-
a
)
t
sin
⁡
(
π
x
/
L
), represented by the solid line in the graph (refer to Figures 2 and 3), and this numerical result supports the continuation of the solution. The series displayed in plots is the partial sum of the first five terms: n
=
5 (summing u
0 to u
4). The optimal convergence parameter ħ in each case is obtained by minimizing the residual error E
m displayed in (23) in the intervals 0
≤
x
≤
π and 0
≤
t
≤
6.0, and the obtained results are shown in Table 1.Table 1
The estimated optimal convergence parameterħ  and the corresponding residual error E
m  at different fractional parameters used in Example 11.
α
,
θ
,
β
Optimalħ
Residual error (E
4
)
α
=
1.5,  β
=
1.0
θ
=
0.0
-
0.990381137
4.6032
E
-
10
θ
=
0.25
-
1.001711905
2.5426
E
-
5
θ
=
-
0.25
-
1.0023206958
2.54262
E
-
6
β
=
1.0,  θ
=
0
α
=
1.8
-
0.9339901124
2.2375
E
-
7
α
=
1.9
-
0.9303835167
3.4254
E
-
7
α
=
1.95
-
0.9285827510
4.2025
E
-
7
α
=
2.0,  θ
=
0
β
=
0.5
-
0.9099205063
9.7340
E
-
9
β
=
0.7
-
0.9151677201
7.4716
E
-
8
β
=
0.9
-
0.9225579219
3.0301
E
-
8Figure 2
The behavior of the four terms of OHAM solution of (49) at t
=
0.2, 0
<
x
<
π, ⁡and the fractional parameters θ
=
0, β
=
1, and α
=
1.8,1.95, compared with the exact solution.Figure 3
The behavior of the four terms of OHAM solution of (49) at t
=
0.2, 0
<
x
<
π, and the fractional parameters θ
=
0, α
=
2.0, and β
=
0.5,0.7,0.9, compared with the exact solution.Example 12.
Consider the nonlinear fractional problem (1) with k
(
u
)
=
u, g
(
u
)
=
4
u
2
-
u, and f
(
x
)
=
sin
⁡
(
2
x
):(57)
D
C
t
β
u
x
,
t
=
u
R
x
α
,
θ
u
x
,
t
+
4
u
2
-
u
,
-
∞
<
x
<
∞
,
t
>
0
,
u
x
,
0
=
sin
⁡
2
x
.Here, we choose the auxiliary linear operator as(58)
L
ϕ
=
∂
∂
t
ϕ
,and the nonlinear operator N is chosen as(59)
N
ϕ
=
D
C
t
β
ϕ
-
ϕ
R
x
α
,
θ
ϕ
-
4
ϕ
2
+
ϕ
.Then, mth-order deformation equation for this problem is given by(60)
∂
∂
t
u
m
x
,
t
-
χ
m
u
m
-
1
x
,
t
=
ħ
H
x
,
t
R
m
u
→
m
-
1
x
,
t
,where R
m
[
u
→
m
-
1
(
x
,
t
)
] is given by(61)
R
m
u
→
m
-
1
x
,
t
=
D
C
t
β
u
m
-
1
-
∑
k
=
0
m
-
1
u
k
R
x
α
,
θ
u
m
-
1
-
k
-
4
∑
k
=
0
m
-
1
u
k
u
m
-
1
-
k
+
u
m
-
1
.We choose H
(
x
,
t
)
=
1 and u
0
=
s
i
n
(
2
x
). Following the same procedure in the previous example to obtain the series solution terms, (62)
u
0
=
sin
⁡
2
x
,
u
1
=
h
t
β
2
Γ
1
+
β
-
4
+
4
cos
⁡
4
x
+
2
α
cos
⁡
π
θ
2
-
2
α
cos
⁡
4
x
-
π
θ
2
+
2
sin
⁡
2
x
.To find the next terms in the solution series, use the following recursive formula which is valid form
>
1:(63)
u
m
=
1
+
ħ
u
m
-
1
+
ħ
J
t
β
-
∑
k
=
0
m
-
1
u
k
R
x
α
,
θ
u
m
-
1
-
k
-
4
∑
k
=
0
m
-
1
u
k
u
m
-
1
-
k
+
u
m
-
1and the solution is thus obtained as(64)
u
=
u
0
+
u
1
+
u
2
+
u
3
+
⋯
.Figures4–6 show the solution (the first five terms in the OHAM series) behavior as the parameters θ, α, and β change at a fixed time t
=
0.2 and in the interval 0
≤
x
≤
π
/
2. As α and/or β increases, the amplitude of the solution decreases. As α tends to 2 and/or β goes up to unity, the series solution approximately coincides with the exact solution of the corresponding integer order equation u
=
e
-
t
sin
⁡
(
2
x
), represented by the solid line in the graph (see Figures 5 and 6), and this confirms the continuation of the solution numerically. The series displayed in plots is the partial sum of the first five terms: n
=
5 (summing u
0 to u
4). The optimal convergence parameter ħ in each case is obtained by minimizing the residual error E
m displayed in (23) in the intervals 0
≤
x
≤
2 and 0
≤
t
≤
0.4, and the obtained results are tabulated in Table 2.Table 2
The estimated optimal convergence parameterħ  and the corresponding residual error E
m  at different fractional parameters used in Example 12.
α
,
θ
,
β
Optimalħ
Residual error (E
4
)
α
=
1.5,  β
=
1.0
θ
=
0.0
-
0.831733694
2.3498
E
-
8
θ
=
0.5
-
0.954364814
1.0016
E
-
5
θ
=
-
0.5
-
0.637748947
7.4733
E
-
6
β
=
1.0,  θ
=
0
α
=
1.92
-
0.589337116
2.3160
E
-
6
α
=
1.95
-
0.589337116
2.1595
E
-
6
α
=
1.99
-
0.639152091
7.6926
E
-
6
α
=
2.0,  θ
=
0
β
=
0.5
-
0.858637427
2.9308
E
-
9
β
=
0.7
-
0.929790675
7.6259
E
-
12
β
=
0.95
-
0.973269800
2.7929
E
-
12Figure 4
The behavior of the four terms of OHAM solution of (57) at t
=
0.2, 0
<
x
<
π
/
2, and the fractional parameters α
=
1.5 and β
=
1 and the skewness θ
=
0,0.5
,
-
0.5.Figure 5
The behavior of the four terms of OHAM solution of (57) at t
=
0.2, 0
<
x
<
π
/
2, and the fractional parameters θ
=
0, β
=
1, and α
=
1.5,1.7,1.9, compared with the exact solution.Figure 6
The behavior of the four terms of OHAM solution of (57) at t
=
0.2, 0
<
x
<
π
/
2, and the fractional parameters θ
=
0, α
=
2.0, and β
=
0.4,0.8,0.95, compared with the exact solution.


## 6. Conclusion

A definition of the fractional-order Riesz-Feller derivative in the Caputo sense is proposed and its equivalence with the classical definition is proved. Then, we proved the continuation of the solution of the fractional-order anomalous diffusion equation with Riesz-Feller spatial derivative to the corresponding Riesz problem. The iterative series solution for the fractional equation is obtained using the OHAM. The advantage of using this technique is the ability to estimate an approximation to the residual error. The results obtained illustrate graphically the continuation of the solution we proved theoretically.



---

*Source: 2371837-2016-08-08.xml*