FUNDAMENTAL SOLUTION TO THE CAUCHY PROBLEM FOR THE TIME-FRACTIONAL ADVECTION-DIFFUSION
EQUATION
Yuriy Povstenko
1, Joanna Klekot
21
Institute of Mathematics and Computer Science, Jan Długosz University in Częstochowa Częstochowa, Poland
2
Institute of Mathematics, Częstochowa University of Technology Częstochowa, Poland
1
j.povstenko@ajd.czest.pl,
2joanna.klekot@im.pcz.pl
Abstract. The one-dimensional time-fractional advection-diffusion equation with the Caputo time derivative is considered. The fundamental solution to the Cauchy problem is obtained using the integral transform technique. The numerical results are illustrated graphically.
Introduction
In the last few decades, a considerable research effort has been expended to the study of fractional differential equations which have many applications in physics, geophysics, geology, chemistry, rheology, engineering and bioengineering (see, for example, [1-8] and references therein). The time-fractional diffusion-wave equa- tion is a mathematical model of important physical phenomena ranging from amorphous, colloid, glassy and porous materials through fractals, random and disordered media to comb structures, dielectrics, semiconductors, polymers and biological systems. The time-fractional advection-diffusion equation can be inter- preted in terms of diffusion or heat conduction with additional velocity field, trans- port processes in porous media and groundwater hydrology. In the case of the frac- tional advection-diffusion equation, as a rule, different numerical methods have been used to find the solution: the implicit difference method based on the shifted Grünwald approximation [9] and the explicit difference method [10], the fractional variational method [11], the finite volume method [12], etc. In [13, 14] the solution to one-dimensional time-fractional advection-diffusion equation was obtained in terms of the H-function.
In this paper, we get the fundamental solution to the Cauchy problem for the time-fractional advection-diffusion equation in terms of the Mittag-Leffler function.
The Laplace transform with respect to time and the exponential Fourier transform
with respect to the space coordinate are used. The results of numerical calculations
are illustrated graphically.
1. Formulation of the problem
Consider the time-fractional advection-diffusion equation
( ) ( ) ( )
x t x v T x
t x a T t
t x T
∂
− ∂
∂
= ∂
∂
∂ , , ,
2 2 α
α
, (1)
where: − ∞ < x < ∞ , 0 < t < ∞ , 0 < α ≤ 1 , a > 0 , v > 0 . In equation (1) ∂ α T / ∂ t α is the Caputo fractional derivative of the order α :
( ) ( t ) T ( ) d n n
n t
T t
n
n n − < <
∂
− ∂
− Γ
=
∂
∂ α α ∫ τ − α − τ τ τ α
α
1 1 ,
0
1 . (2)
Here Γ ( ) x is the gamma function.
In what follows we will study the fundamental solution to the Cauchy problem for the time-fractional advection-diffusion equation (1) with the initial condition
( x ) p ( ) x
T δ
0 0
, = , (3)
where δ ( ) x is the Dirac delta function. In the initial condition (3) we have intro- duced the constant multiplier p 0 to obtain the nondimensional quantity displayed in figures.
The zero condition at infinity is also assumed:
( ) , 0
lim =
∞
±
→
t x T
x
. (4)
2. Method of solution
To solve the Cauchy problem under consideration we use the Laplace transform with respect to time t . For the function T (t ) , 0 < t < ∞ , this transform is defined as
{ } = = ∫ ∞ −
0
* ( ) ( )
)
( t T s e T t dt
T
L s t (5)
with the inverse carrying out according to the Fourier-Mellin formula
{ } ( ) + ∫ ∞
∞
−
− = =
i c
i c
t
s T s ds e
i t T s T
L ( ) ,
2 ) 1
( *
* 1
π (6)
where c is a fixed positive number.
For the Laplace transform rule the Caputo fractional derivative requires the knowledge of the initial values of the function and its integer derivatives of the order k = 1 , 2 , K , n − 1 :
. 1
, ) 0 ( )
) (
( 1
0
) ( 1
* s s T n n
T s t
d t T L d
n
k
k
k − < ≤
−
=
∑ −
=
+
−
− α
α α
α α
(7)
Applying the Laplace transform to equation (1) while taking into account the initial condition (3) and the rule (7) gives
( ) ( )
x s x v T x
s x a T x p s s x T
s ∂
− ∂
∂
= ∂
− − , ,
) ( )
,
( 2
* 2 0
1
* α δ
α . (8)
Next, we use the exponential Fourier transform with respect to the spatial coordi- nate x:
{ } ∫ ∞
∞
−
=
= T e T x dx
x T
F ix ( )
2 ) 1
~ ( )
( ξ
ξ π , (9)
{ } ( ) ∫ ∞
∞
−
−
− = = ~ ( ) .
2 ) 1
~ (
1 ξ ξ
π
ξ T x e ξ T d
T
F ix (10)
The Fourier transform of the derivative of a function is defined by the following relation:
)
~ ( ) ) (
( i ξ T ξ
t d
x T
F d n
n
n = −
. (11)
Application of the exponential Fourier transform (9) to equation (8) using (11) leads to
) ,
~ ( )
,
~ ( ) 2
,
~ * ( 1 0 2 * *
s T v i s T p a
s s T
s ξ ξ ξ ξ
ξ α π
α − − = − + , (12)
and in the transform domain we get
ξ π ξ
ξ α
α
v i a s
s s p
T
−
= +
− 2 1
* 0
) 2 ,
~ (
. (13)
Inversion of the Laplace and Fourier transforms results in the solution:
( )
[ ]
∞ ∫
∞
−
− −
−
= ξ ξ ξ
π
ξ α
α a i v t e d
p E t x
T 0 2 i x
) 2 ,
( . (14)
To obtain (14) the following formula [1]
( α )
α α α
bt E b s
L s = −
+
−
−
1 1 (15)
has been used, where E α (z ) is the Mittag-Leffler function in one parameter α having the series representation
( 1 ) , 0 .
) (
0
∑ ∞
=
+ >
= Γ
n
n
n z z
E α
α
α (16)
For negative real values of the argument the Mittag-Leffler function E α ( x − ) can be calculated as the integral [15]:
( )
( )
∞ ∫ −
− < <
+
= +
−
0
2
1
1 0 1 , cos
2 ) sin
(
/ 1
α π
π α π α
α α
α α
α
du u
u e u x
E x u . (17)
Next, we recast the argument of the Mittag-Leffler function in the solution (14) as
∞ ∫
∞
−
−
+
−
−
= ξ ξ
π
ξ α
α t e d
a v a i v a
p E t x
T 2 i x
2 2 0
4 2
) 2 ,
( (18)
and use the substitution
a i v
− 2
η = ξ (19)
which allows us to rewrite the solution (18) in the form
∫ ∞
∞
−
+
−
+
−
= η η η
π
α
α d
a i v x i a t
a v p E
t x
T exp 2
4 ) 2
,
( 2
2
0 2 (20)
and consequently as
∫ ( )
∞
∞
−
+
−
= η η η
π
α
α t x d
a a v
a E x v t p
x
T cos
2 4 2 exp ) ,
( 2
2 2
0 . (21)
In the particular case α = 1 the Mittag-Leffler function E 1 ( ) z = e z , hence for the
classical advection-diffusion equation we have
∫ ( )
∞
∞
−
+
−
= η η η
π
d x a t
a v a
x v t p
x
T cos
exp 4 exp 2
) 2 ,
( 2
2 2
0 (22)
or
∫ ( )
∞
∞
−
−
−
= cos .
exp 4 exp 2
) 2 ,
(
22
0 η η
π
η x d
a e t v a
x v t p
x
T a t (23)
Since [16]
( )
−
∫ =
∞
∞
−
−
t a x t
a d x e a t
exp 4 cos
2