Next Article in Journal
CIST: Differentiating Concepts and Instances Based on Spatial Transformation for Knowledge Graph Embedding
Next Article in Special Issue
Technology for Obtaining the Approximate Value of Moving Singular Points for a Class of Nonlinear Differential Equations in a Complex Domain
Previous Article in Journal
Entropy Analysis of Sutterby Nanofluid Flow over a Riga Sheet with Gyrotactic Microorganisms and Cattaneo–Christov Double Diffusion
Previous Article in Special Issue
Mathematical Model for Analyzing the Dynamics of Tungro Virus Disease in Rice: A Systematic Literature Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Boundary Value Problem of Space-Time Fractional Advection Diffusion Equation

by
Elsayed I. Mahmoud
1,2,* and
Temirkhan S. Aleroev
3
1
Department of Mathematics, Faculty of Science, Zagazig University, Zagazig 44519, Egypt
2
Nikolskii Mathematical Institute, Peoples Friendship University of Russia, Moscow 117198, Russia
3
Department of Applied Mathematics, Moscow State University of Civil Engineering, Yaroslavskoe Shosse, 26, Moscow 129337, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(17), 3160; https://doi.org/10.3390/math10173160
Submission received: 29 July 2022 / Revised: 13 August 2022 / Accepted: 31 August 2022 / Published: 2 September 2022
(This article belongs to the Special Issue Dynamical Systems and Optimal Control, 2nd Edition)

Abstract

:
In this article, the analytical and numerical solution of a one-dimensional space-time fractional advection diffusion equation is presented. The separation of variables method is used to carry out the analytical solution, the basis of the system eigenfunction and their corresponding eigenvalue for basic equation is determined, and the numerical solution is based on constructing the Crank-Nicolson finite difference scheme of the equivalent partial integro-differential equations. The convergence and unconditional stability of the solution are investigated. Finally, the numerical and analytical experiments are given to verify the theoretical analysis.

1. Introduction

Recently, interest in the topics of fractional calculus has increased [1]. These topics are used by many scientists to simulate a variety of physical, biological, and chemical processes [2,3]. The space-time fractional advection diffusion equation (STFADE) considers a famous model for anomalous diffusion in several practical fields, such as the diffusion processes of pollutants in porous media and the continuous-time random walk. In [4,5], Gorenflo, Mainardi and Luchko developed discrete random walk models and found that the fundamental solution may be regarded as the probability density of a self-similar random process evolving in time, which can be viewed as a generalized diffusion process.
The use of fractional derivatives (FDs) for modeling real physical processes or environments leads to the appearance of equations containing derivatives and integrals of fractional order in addition to the classical ones. Researchers have focused their efforts on fractional-order physical models of non-equilibrium processes [6,7]. As a result, they study the linear elasticity and viscoelastic deformation processes in biomaterials and capillary-porous materials with fractal structures. Recent applications of fractional equations to many biological systems are more adequate than the previously used integer-order biological systems, because fractional-order derivatives and integrals describe the memory and hereditary properties of different substances [8,9,10,11,12].
This paper is concerned with methods for obtaining accurate numerical and analytical approximations to the solutions of fractional sub-diffusion equations (the ultra-slow diffusion model with time fractional derivative 0 < α < 1 ), which is even more slow than the power-law sub-diffusion and is widely observed in a variety of natural and engineering fields. In this paper, we investigate the one-dimensional STFADE with the source term
α z ( ρ , t ) t α = 0 R D ρ β z ( ρ , t ) + s ( ρ , t ) , z ( 0 , t ) = z ( 1 , t ) = 0 , 0 < t < T , z ( ρ , 0 ) = ϕ ( ρ ) , 0 < ρ < 1 ,
in domain ω = ( 0 , 1 ) × ( 0 , T ) ,   0 < T < , where 0 < α < 1 is a parameter describing the order of the fractional derivative in the time direction, and 1 < β < 2 describes the order of the fractional derivative space direction. α z ( ρ , t ) t α is the Caputo fractional derivative of order α , which is defined by
α z ( ρ , t ) t α = 1 Γ ( 1 α ) 0 t ( t σ ) α z ( ρ , σ ) σ d σ , 0 < α < 1 , 0 < σ < t z ( ρ , t ) t α = 1 ,
and the space fractional derivative D ρ β is the Riemann–Liouville fractional derivative operator of order β defined by
D ρ β z ( ρ , t ) = 1 Γ ( 2 β ) d 2 d ρ 2 0 ρ ( ρ σ ) 1 β z ( σ , t ) d σ , 1 < β < 2 , 0 < σ < ρ 2 z ( ρ , t ) ρ 2 β = 2 .
Equation (1) yields the following forms of diffusion equations for various values of the parameters α and β :
For α = 1 and β = 2 ⟶ the classical diffusion equation;
For 0 < α < 1 and β = 2 ⟶ the time-fractional diffusion equation;
For α = 1 and 0 < β < 2 ⟶ the space-fractional diffusion equation,
For 0 < α < 1 and 0 < β < 2 ⟶ the space-time fractional diffusion equation.
This work establishes the fundamental analytical and numerical solution of (1), which is a very general form of STFADE, and many problems that are already studied are special cases of it, as in [13,14,15].
The structure of the paper is organized as follows. In Section 2, some theoretical definitions and proprieties of fractional differentiable classes and Mittage–Leffler functions are presented. In Section 3, we discuss the analytical solution of (1) with the separation of variable method and construct the base system of eigenfunctions and their corresponding eigenvalues of the basic equation. In Section 4 and Section 5, the numerical scheme for solving TSFADE and their convergence and unconditional stability analysis are established. Finally, in Section 6, the numerical experiments are presented.

2. Theoretical Preliminaries

Definition 1
([16]). A real-valued function Ψ ( τ ) can be said to belong to space C α [ 0 , ] ,   α I R if there exists a real number q > α such that Ψ ( τ ) = τ q Ψ ˜ ( τ ) , Ψ ˜ ( τ ) C [ 0 , ] .
Definition 2
([16]). A real valued function Ψ ( τ ) can be said to belong to space C α m [ 0 , ] ,   α I R ,   m I N { 0 } if Ψ m C α [ 0 , ] .
Definition 3
([17]). Let E α , β ( w ) be the two-parameter Mittage–Leffler function on w defined by
E α , β ( w ) = k = 0 w k Γ ( α k + β ) , R e ( α ) > 0 , w , β C .
As an example, it is easy to deduce that the two-parameter Mittag–Leffler function E α , β ( w ) is a generalization of the exponential function ( E 1 , 1 ( w ) = e x p ( w ) ), hyperbolic functions ( E 2 , 1 ( w 2 ) = cosh ( w ) ) and trigonometric functions ( E 2 , 2 ( w 2 ) = sin ( w ) / w , E 2 , 1 ( w 2 ) = cos ( w ) ).
The validity of the following relationships on E α , β ( w ) can be easily tested [15]:
1.
0 t σ α 1 E α , α ( λ σ α ) ( t σ ) β 1 E α , β ( μ ( t σ ) α ) d σ = t β 1 E α , β ( μ t α ) E α , β ( λ t α ) λ μ , λ μ t α + β 1 ( 1 α E α , α + β 1 ( λ t α ) + ( 1 β α ) E α , α + β ( λ t α ) ) λ = μ ,
2.
0 t ( t τ ) ρ 1 τ β 1 E α , β ( λ τ α ) d τ = Γ ( ρ ) t ρ + β 1 E α , β + ρ ( λ t α )
Lemma 1
([16]). If α < 2 , β R , θ is such that π α / 2 < θ < min ( π , π α ) , w C such that | w | 0 , θ | arg ( w ) | π ; then
| E α , β ( w ) | C 1 + | w |

3. Analytical Solution of STFADE

In this section, we proceed with a discussion of the theorem of the existence of aa solution of the boundary value problem for Equation (1) by using the separation of variables method (Fourier method).
Theorem 1.
For every s ( ρ , t ) C α , β ( 1 ) , ( 2 ) ( ω ) and ϕ ( ρ ) C ( 0 , 1 ) , the regular solution z ( ρ , t ) C α , β ( 1 ) , ( 2 ) ( ω ) of the boundary value problems 1 can be represented in the form
z ( ρ , t ) = m = 1 ϕ m E α , 1 ( λ m t α ) + 0 t τ α 1 E α , α ( λ m τ α ) s m ( t τ ) d τ ρ β 1 E β , β ( λ m ρ β ) .
Here, E α , β ( z ) -is the well-known function of two parameters of the Mittag–Leffler type: ϕ m , s m ( t ) —the decomposition coefficients of functions ϕ ( ρ ) and s ( ρ , t ) , respectively, on the basis of functions P m ( ρ , λ m ) = ρ β 1 E β , β ( λ m ρ β ) .
To solve the boundary value problem (1), we represent the function z ( ρ , t ) in the form
z ( ρ , t ) = z 1 ( ρ , t ) + z 2 ( ρ , t )
According to (7), the boundary value problem (1) can be divided into two boundary value problems as follows:
α z 1 ( ρ , t ) t α = 0 R D ρ β z 1 ( ρ , t ) , z 1 ( 0 , t ) = z 1 ( 1 , t ) = 0 , 0 < t < T , z 1 ( ρ , 0 ) = ϕ ( ρ ) , 0 < ρ < 1 ,
and
α z 2 ( ρ , t ) t α = 0 R D ρ β z 2 ( ρ , t ) + s ( ρ , t ) , z 2 ( 0 , t ) = z 2 ( 1 , t ) = 0 , 0 < t < T , z 2 ( ρ , 0 ) = 0 , 0 < ρ < 1 .
In general, as is customary in solving such problems by the Fourier method, we first consider the homogeneous fractional partial differential Equation (8). The essence of this problem is to find a non-trivial solution represented as product z 1 = P ( ρ ) T ( t ) ; therefore, Equation (8) can be valid in the case
c D t α T ( t ) T ( t ) = R D ρ β P ( ρ ) P ( ρ ) = λ ,
Then, we can obtain the following fractional differential equations:
R D ρ β P ( ρ ) + λ P ( ρ ) = 0 , P ( 0 ) = P ( 1 ) = 0 , 1 < β < 2 ,
C D t α T ( t ) + λ T ( t ) , 0 < α < 1 .
Now, we can discuss the existence of a solution of Equation (8) by formulating it in the sense of the eigenfunctions of the Sturm–Liouville eigenvalue problem (10) by using a Fourier series, which is studied in [13,18], with the help of the following lemma.
Lemma 2
([13]). For the Sturm–Liouville eigenvalue problem (10), the following statement holds:
1.
The corresponding eigenfunctions are given by P m ( ρ ) = ρ β 1 E β , β ( λ m ρ β ) , where the eigenvalues { λ m } are zeros of the Mittage–Leffler function E β , β ( λ ) with I m ( λ m ) < 0 , | λ l | | λ l + 1 | , l 1 .
2.
The system of eigenfunctions { ρ β 1 E β , β ( λ m ρ β ) } m = 1 is complete in L 2 ( 0 , 1 )
3.
All eigenvalues { λ m } are in the sector { w C : | a r g w | < π β π 2 } .
Remark 1.
In our work to specify the eigenvalue of problem (10), using Wolfram Mathematica (high-level technical computing language) to calculate the values of λ m at different β, according to Lemma 2, we can rearrange all the zeros of E β , β ( λ m ) as follows: λ 2 , λ 1 , λ 1 , λ 2 , . We observe, as shown in Table 1, that for β 1 , there is only a single real value and corresponding eigenfunction as shown in Figure 1, and with increasing β, we observe an increase in the number of real eigenvalues λ m with no complex when β is close to 2.
The system { P m ( ρ ) } m = 1 of eigenfunctions is complete but not orthogonal [19]. Thus, together with the system { P m ( ρ ) } n = 1 , we will consider the system { Υ m ( ρ ) } m = 1 = { ( 1 ρ ) β 1 E β , β ( λ m ( 1 ρ ) β ) } m = 1 , which represents the eigenfunctions of the adjoint of (10) and is biorthogonal with { P m ( ρ ) } m = 1 .
Equation (11) is solved in [20,21] by using the Laplace transform of the Caputo time fractional differential operator (2), where we obtain the solution in terms of the Mittage–Liffler function as T m = ϕ m E α , 1 ( λ m t α ) , where the Fourier coefficient ϕ m is given by
ϕ m = ϕ ( ρ ) , Υ m ( ρ ) P m ( ρ ) , Υ m ( ρ ) .
Now, we can construct the solution of problem (8), where
z 1 ( ρ , t ) = m = 1 T m ( t ) P m ( ρ ) = m = 1 ϕ m E α , 1 ( λ m t α ) ρ β 1 E β , β ( λ m ρ β ) .
By using Lemma (1) and the fact that for T > 0 ,
| 1 E α , 1 ( λ m T α ) 1 | C 1 ,
C 1 is positive constant, where E α , 1 ( λ m T α ) 1 [17]; it is not difficult to prove the solution of (13) is absolutely convergent.
Now we go to discuss the solution of Equation (9), z 2 ( ρ , t ) , which can be found by using the full basis of the eigenfunctions P m
z 2 ( ρ , t ) = m = 1 u m ( t ) P m ( ρ ) , s ( ρ , t ) = m = 1 s m ( t ) P m ( ρ ) ,
where s m ( t ) = s ( ρ , t ) , Υ m ( ρ ) P m ( ρ ) , Υ m ( ρ ) . Then, from (9), we obtain
C D t α u m ( t ) + ( λ m ) u m ( t ) = s m ( t ) , 0 < α < 1 , u m ( 0 ) = 0 .
Applying the Laplace transform operator to both sides of Equation (16) [1], we have
u m ( t ) = 0 t τ α 1 E α , α ( λ m τ α ) s m ( t τ ) d τ ,
then, from (15), z 2 ( ρ , t ) can formalize in the form
z 2 ( ρ , t ) = m = 1 0 t τ α 1 E α , α ( λ m τ α ) s m ( t τ ) d τ ρ β 1 E β , β ( λ m ρ β ) ,
Thus, from (7), (13) and (18), we obtain (6).

4. Numerical Solution of STFADE

In this section, the Crank—Nicholson finite difference technique is applied to obtain the numerical solution of one-dimensional STFADE (1) and present a numerical approximation for the Riemann–Liouville derivative. We find that if we assume for Equation(1) an equivalent partial integro-differential equation form, the precision of the discrete approximations could be improved. Let us try to take the fractional Riemman–Liouville integral j t α , 0 < α < 1 , on both sides of Equation (1) such that
j t α f ( ρ , t ) = 1 Γ ( α ) 0 t f ( ρ , ξ ) ( t ξ ) 1 α d ξ ,
Thus, we obtain
z ( ρ , t ) = z ( ρ , 0 ) + j t α ( 0 R D ρ β z ( ρ , t ) ) + s ( ρ , t ) ,
where S ( ρ , t ) = j t α ( s ( ρ , t ) ) . In order to carry out discretizations, we introduced the temporal step size τ = T , Z / { 0 } and uniform mesh of interval [ 0 , T ] , ω τ = { t n :   = n τ ; 0 < n < } . Additionally, for a spatial discretization, let = 1 M , M Z / { 0 } and a uniform mesh of interval [ 0 , 1 ] , ω = { ρ i :   = i ;   0 < i < M } . Suppose on ω × ω τ there exists a grid function W = z i n | 0 i M , n 0 , such that for any z , v W , we define the following norms, semi-norm . L 2 and the inner product [22], as follows:
z i n + 1 2 = 1 2 z i n + 1 + z i n , z n , g n = h i = 1 M 1 z i n g i n , z n L 2 2 = z n , z n z n = max 0 i M | z i n | , δ ρ β z n , g n = δ ρ β / 2 z n , δ ρ β / 2 g n .
Additionally, we can represent the Riemann–Lioville fractional operator of any spatially dependent function f ( ρ ) C 4 ( R ) by shifted the Grunwald–Letnikov formula as
0 R D ρ β f ( ρ i ) = 1 Γ ( 4 β ) β l = 0 i + 1 q l β f ( ρ i l + 1 ) + O ( 2 ) ,
where
q l β = 1 , l = 0 2 3 β 4 , l = 1 3 3 β 4 . 2 3 β + 6 , l = 2 ( l + 1 ) 3 β 4 l 3 β + 6 ( l 1 ) 3 β 4 ( l 2 ) 3 β + ( l 3 ) 3 β , l 3 .
For time discretization, the high-order accuracy of the suggested method presented on the first-order approximation of the fractional Riemann–Liouville integral j t α , this approximation established for every f ( t ) C α ( [ 0 , T ] ) that
j t α f ( t n ) = τ α k = 0 n p k α f ( t n k ) + O ( τ ) ,
where p 0 α = 1 , p k α = ( 1 ) k α k .
Lemma 3
([23]). The sequence { q l β } l = 0 , 1 < β < 2 , which is defined in (22), has the following properties
l = 0 q l β = 0 , q 1 β 0 , q 0 β 3 β 0 , q 0 β + q 2 β 0 , q 2 ( β ) > 0 , β ( β 0 , 2 ) 0 , β ( 1 , β 0 ) ,
where β 0 is the zeros of q 2 β .
Lemma 4
([24]). The sequence { p k α } k = 0 , 0 < α < 1 , which is defined in (23), satisfies
p k α > 0 , p k + 1 α < p k α , p k + 1 α + p k 1 α > 2 p k α , k 1 .
Lemma 5
([25]). Consider the sequence { p k α } k = 0 , 0 < α < 1 , which is defined in (22). Then, for any real vector ( Υ 1 , Υ 2 , , Υ j , ) T R j , it holds that
m = 0 j 1 k = 0 m p k α Υ m + 1 k Υ m + 1 0
Therefore, according to discretization, ω and ω τ , a weighted Crank–Nicolson method for Equation (20) at the point ( ρ i , t n + 1 2 ) , produces
z ( ρ i , t n + 1 ) + z ( ρ i , t n ) = 2 z ( ρ i , 0 ) + τ α ( k = 0 n + 1 p k α δ ρ β z ( ρ i , t n + 1 k ) + k = 0 n p k α δ ρ β z ( ρ i , t n k ) ) + S ( ρ i , t n + 1 ) + S ( ρ i , t n ) + R i n + 1 ,
where R i n + 1 < C ( τ + 2 ) , C > 0 , 0 < i < M 1 , 0 < n < 1 .
Let z i n represent the numerical solution of z ( ρ i , t n ) , z i 0 = z ( ρ i , 0 ) and S i n = S ( ρ i , t n ) . Tearranging the above Formula (27), we obtain the following weighted Crank–Nicolson numerical scheme for the one-dimensional STFADE (1):
z i n + 1 + z i n = 2 z i 0 + τ α k = 0 n + 1 p k α δ ρ β z i n + 1 k + k = 0 n p k α δ ρ β z i n k + S i n + 1 + S i n , 1 < i < M 1 , 1 < n < 1 , z 0 n = z M n = 0 , 0 < n < , z i 0 = ϕ i , 0 < i < M

5. Analysis of the Numerical Solution

In this section, we will discuss the convergence and stability of the suggested numerical scheme (28).
Theorem 2.
Consider z ( ρ , t ) C α , 2 ( ω ) , the solution of Equation (1), and consider that { Z i n | 0 i M , 0 n } is the solution of the numerical scheme (28). Then, for small τ and ℏ, it holds that
Λ n C ( τ + 2 ) ,
where Λ i n = z ( ρ i , t n ) Z i n .
Proof. 
Subtracting Equation (28) from Equation (27), we have
Λ i n + 1 + Λ i n = τ α k = 0 n + 1 p k α δ ρ β Λ i n + 1 k + k = 0 n p k α δ ρ β Λ i n k + R n + 1 i , Λ i 0 = 0 , 0 < i < M 1 , 0 < n < 1 .
Equation (30) can be rewritten in a matrix form as
Λ n + 1 + Λ n = k = 0 n p k α A ( Λ n + 1 k + Λ n k ) + R n + 1 ,
where
Λ n = ( Λ 1 n , Λ 2 n , , Λ M 1 n ) T , A = τ α Γ ( 4 β ) β q 1 ( β ) q 0 ( β ) 0 0 q 2 ( β ) q 1 ( β ) q 0 ( β ) q 2 ( β ) q 1 ( β ) q 0 ( β ) q ( β ) q 1 ( β ) q 2 ( β ) q 1 ( β ) .
By the Gershgorin theorem and Lemma 3, we could investigate into the matrix A is negative definite matrix.
Multiplying Equation (31) by ( Λ n + 1 Λ n ) T in both sides, we obtain
( Λ n + 1 Λ n ) T ( Λ n + 1 + Λ n ) = k = 0 n p k α ( Λ n + 1 Λ n ) T A ( Λ n + 1 k + Λ n k ) + ( Λ n + 1 Λ n ) T R n + 1
By using Lemma 5 and the fact that ( Λ n + 1 Λ n ) T A ( Λ n + 1 k + Λ n k ) < 0 , we can deduce that
Λ n + 1 2 Λ n 2 R n + 1 , ( Λ n + 1 Λ n ) ,
By summing over n from 0 to J 1 , it is deduced that
Λ J 2 n = 0 J 1 R n + 1 , ( Λ n + 1 Λ n ) R J , Λ J C ( τ + 2 ) Λ .
Thus, we deduce (27). □
Theorem 3.
The numerical solution { Z i n | 0 i M , 0 n } of scheme (28) is unconditionally stable and holds that
Z J 2 2 Z 0 2 + n = 0 J 1 ( 1 + τ α p n + 1 α δ ρ β ) Z 0 + S n + 1 2 2
Proof. 
Similar to the proof of Theorem 2, by multiplying Equation (28) by ( Z i n + 1 Z i n ) and obtaining the sum over i from 1 M 1 ,
Z n + 1 2 Z n 2 = 2 Z 0 , Z n + 1 Z n + τ α p n + 1 α δ ρ β Z 0 , Z n + 1 Z n + τ α k = 0 n p k α δ ρ β ( Z n + 1 k + Z n k ) , Z n + 1 Z n + 1 2 S n + 1 2 , Z n + 1 Z n .
After applying the Cauchy=-Schwarz inequality and obtaining the sum of Equation (34) over n from 0 to J 1 , we obtain
Z J 2 Z 0 2 = 2 n = 0 J 1 Z 0 , Z n + 1 Z n + τ α n = 0 J 1 p n + 1 α δ ρ β Z 0 , Z n + 1 Z n + τ α n = 0 J 1 k = 0 n p k α δ ρ β ( Z n + 1 k + Z n k ) , Z n + 1 Z n + 1 2 n = 0 J 1 S n + 1 2 , Z n + 1 Z n .
Without loss of generality, according to Lemmas 4 and 5, we can infer that the second term on the right side of Equation (35) is negative, and by applying Young’s inequality, we obtain (33). □

6. Examples

In this section, we present an analytical and numerical solution of some examples for the one-dimensional STFADE (1). For particular values of α , β , s ( ρ , t ) , ϕ ( ρ ) , we demonstrate the computational performance and theoretical findings of our suggested methods.
Example 1.
Consider
α z ( ρ , t ) t α = 0 R D ρ β z ( ρ , t ) + ( ρ + ρ 2 β ) E α , 1 ( t α ) , 0 < α < 1 , 1 < β < 2 z ( 0 , t ) = z ( 1 , t ) = 0 , 0 < t < T , z ( ρ , 0 ) = ρ ( 1 ρ ) , 0 < ρ < 1 .
Especially with (1), by setting s ( ρ , t ) = ( ρ + ρ 2 β ) E α , 1 ( t α ) and ϕ ( ρ ) = ρ ( 1 ρ ) , we obtain (36), so the analytical solution of Equation (36) can be obtained using Theorem 1.
From the relation (4) and (5), one obtains
P m ( ρ ) , Υ m ( ρ ) = 0 1 ρ β 1 E β , β ( λ m ρ β ) ( 1 ρ ) β 1 E β , β ( λ m ( 1 ρ ) β ) d ρ = 1 β E β , 2 β 1 ( λ m ) ( 1 β 1 ) E β , 2 β ( λ m ) ,
ϕ ( ρ ) , Υ m ( ρ ) = 0 1 ( ρ ρ 2 ) ( 1 ρ ) β 1 E β , β ( λ m ( 1 ρ ) β ) d ρ = E β , β + 2 ( λ m ) 2 E β , β + 3 ( λ m ) ,
s ( ρ , t ) , Υ m ( ρ ) = E α , 1 ( t α ) 0 1 ( ρ + ρ 2 β ) ( 1 ρ ) β 1 E β , β ( λ m ( 1 ρ ) β ) d ρ = E α , 1 ( t α ) ( E β , β + 2 ( λ m ) + Γ ( 3 β ) E β , 3 ( λ m ) ) ,
thus, from (12) and (15) can get ϕ m and s m ( t ) , therefor,
0 t τ α 1 E α , α ( λ m τ α ) s m ( t τ ) d τ = β ( E β , β + 2 ( λ m ) + Γ ( 3 β ) E β , 3 ( λ m ) ) E β , 2 β 1 ( λ m ) ( 1 β ) E β , 2 β ( λ m ) 0 t τ α 1 E α , α ( λ m τ α ) E α , 1 ( ( t τ ) α ) d τ = β ( E β , β + 2 ( λ m ) + Γ ( 3 β ) E β , 3 ( λ m ) ) ( E α , 1 ( t α ) E α , 1 ( λ m t α ) ) ( E β , 2 β 1 ( λ m ) ( 1 β ) E β , 2 β ( λ m ) ) ( λ m 1 ) .
Therefore,
z 1 ( ρ , t ) = m = 1 β ( E β , β + 2 ( λ m ) 2 E β , β + 3 ( λ m ) ) E β , 2 β 1 ( λ m ) ( 1 β ) E β , 2 β ( λ m ) E α , 1 ( λ m t α ) ρ β 1 E β , β ( λ m ρ β )
and
z 2 ( ρ , t ) = m = 1 β ( E β , β + 2 ( λ m ) + Γ ( 3 β ) E β , 3 ( λ m ) ) ( E α , 1 ( t α ) E α , 1 ( λ m t α ) ) ( E β , 2 β 1 ( λ m ) ( 1 β ) E β , 2 β ( λ m ) ) ( λ m 1 ) ρ β 1 E β , β ( λ m ρ β )
from (7), we get the analytical solution z ( ρ , t ) of (36).
The compression of the proposed numerical scheme (28) for solving (36) with the numerical solution Z i n and its analytical solution z ( ρ , t ) for β = 1.9 and different α values is shown in Figure 2 at t = T = 1 . Figure 3 shows the graphical analytical and numerical solution of (36). The error is analyzed as the following L 2 -norm
Λ ( τ , ) = i = 1 M 1 | z ( ρ i , t ) Z i | 2 ,
and the order of the convergence can be calculated as
R t log 2 [ Λ ( τ , ) Λ ( 2 τ , ) ] ,
Furthermore, in Table 2, we fix the space step at = 0.02 , T = 1 , and analyse how the error of the proposed numerical scheme for solving (36) ( Λ ( τ , ) ) changes with the maximum less than O ( τ , 2 ) .

7. Conclusions

In this paper, we intend to consider the general form of STFADE with time and fractional derivatives defined in the Caputo and Riemann–Liouville sense by fractional order derivatives 0 < α < 1   a n d   1 < β < 2 , respectively. The analytical solution of one dimension STFADE is carried out by Fourier methods. The justification of the basis for the system of eigenfunctions and their corresponding eigenvalues for the basic Equation (10) is based on the asymptotic properties of zeros of the two-parameter Mittag–Leffler function. The Crank–Nicolson finite difference technique is applied to construct the unconditional stable proposed scheme of the equivalent partial integro-differential Equation (20). The proposed numerical scheme is convergent with first-order accuracy in the temporal direction and second-order accuracy in the spatial direction. All the numerical experiments support our theoretical results. Model (1) may be considered in future work with different definitions of fractional derivatives and different methods for numerical and analytical solutions.

Author Contributions

Conceptualization, E.I.M.; Data curation, E.I.M.; Formal analysis, E.I.M.; Investigation, E.I.M.; Methodology, E.I.M. and T.S.A.; Project administration, T.S.A.; Resources, E.I.M. and T.S.A.; Software, E.I.M.; Supervision, T.S.A.; Validation, T.S.A.; Writing—original draft, E.I.M.; Writing—review & editing, E.I.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  2. Sabatier, J.; Agrawal, O.P.; Machado, A.T. Advances in Fractional Calculus. In Theoretical Developments and Applications in Physics and Engineering; Springer: Dordrecht, The Netherlands, 2007. [Google Scholar]
  3. Aleroev, M.; Aleroev, H.; Aleroev, T. Proof of the completeness of the system of eigenfunctions for one boundary-value problem for the fractional differential equation. AIMS Math. 2019, 4, 714–720. [Google Scholar] [CrossRef]
  4. Gorenflo, R.; Mainardi, F. Random Walk Models for Space Fractional Diffusion Processes. Fract. Calc. Appl. Anal. 1998, 1, 167–191. [Google Scholar]
  5. Mainardi, F.; Luchko, Y.; Pagnini, G. The fundanental solution of the space-time fractional diffusion equation. Fract. Calc. Appl. Anal. 2001, 4, 153–192. [Google Scholar]
  6. Shymanskyi, V.; Sokolovskyy, Y. Variational Method for Solving the Viscoelastic Deformation Problem in Biomaterials with Fractal Structure. In Proceedings of the Information Technology and Implementation (IT&I-2021), CEUR Workshop Proceedings, online, 1–3 December 2021; pp. 125–134. [Google Scholar]
  7. Shymanskyi, V.; Sokolovskyy, Y.; Boretska, I.; Sokolovskyy, I.; Markelov, O.; Storozhuk, O.M. Application of FEM with Piecewise Mittag-Leffler Functions Basis for the Linear Elasticity Problem in Materials with Fractal Structure. In Proceedings of the 2021 IEEE XVIIth International Conference on the Perspective Technologies and Methods in MEMS Design (MEMSTECH), Lviv, Ukraine, 12–16 May 2021; pp. 16–19. [Google Scholar]
  8. Alidousti, J.; Ghafari, E. Dynamic behavior of a fractional order prey-predator model with group defense. Chaos Solitons Fractals 2020, 134, 109688. [Google Scholar] [CrossRef]
  9. Rihan, F.A.; Rajivganthi, C. Dynamics of fractional-order delay differential model of prey-predator system with Holling-type III and infection among predators. Chaos Solitons Fractals 2020, 141, 110365. [Google Scholar] [CrossRef]
  10. Alidousti, J.; Ghahfarokhi, M.M. Stability and bifurcation for time delay fractional predator prey system by incorporating the dispersal of prey. Appl. Math. Model. 2019, 72, 385–402. [Google Scholar] [CrossRef]
  11. Xu, C.; Zhang, W.; Aouiti, C.; Liu, Z.; Yao, L. Further analysis on dynamical properties of fractional-order bi-directional associative memory neural networks involving double delays. Math. Methods Appl. Sci. 2022, 1–19. [Google Scholar]
  12. Alidousti, J. Stability and bifurcation analysis for a fractional prey–predator scavenger model. Appl. Math. Model. 2020, 81, 342–355. [Google Scholar] [CrossRef]
  13. Aleroev, T. Solving the Boundary Value Problems for Differential Equations with Fractional Derivatives by the Method of Separation of Variables. Mathematics 2020, 8, 1877. [Google Scholar] [CrossRef]
  14. Mahmoud, E.; Orlov, V.N. Numerical Solution of Two Dimensional Time-Space Fractional Fokker Planck Equation with Variable Coefficients. Mathematics 2021, 9, 1260. [Google Scholar] [CrossRef]
  15. Sandev, T.; Zivorad, T. The general time fractional wave equation for a vibrating string. J. Phys. A Math. Theor. 2010, 43, 055204. [Google Scholar] [CrossRef]
  16. Aleroev, T.S.; Elsayed, A.M. Analytical and Approximate Solution for Solving the Vibration String Equation with a Fractional Derivative. Mathematics 2020, 8, 1154. [Google Scholar] [CrossRef]
  17. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions Related Topics and Applications; Springer: New York, NY, USA, 2014. [Google Scholar]
  18. Aleroev, T.S.; Elsayed, A.M.; Mahmoud, E.I. Solving one dimensional time-space fractional vibration string equation. IOP Mater. Sci. Eng. Conf. Ser. 2021, 1129, 012030. [Google Scholar] [CrossRef]
  19. Aleroev, T.S.; Kirane, M.; Tang, Y. The boundary-value problem for a differential operator of fractional order. J. Math. Sci. 2013, 194, 499–512. [Google Scholar] [CrossRef]
  20. Ali, M.; Aziz, S.; Malik, S. Inverse source problems for a space-time fractional differential equation. Inverse Probl. Sci. Eng. 2020, 28, 47–68. [Google Scholar] [CrossRef]
  21. Al-Refai, M.; Luchko, Y. Maximum Principle for the Multi-Term Time-Fractional Diffusion Equations with the Riemann-Liouville Fractional Derivatives. Appl. Math. Comput. 2015, 257, 40–51. [Google Scholar] [CrossRef]
  22. Wang, P.D.; Huang, C.M. An energy conservative difference scheme for the nonlinear fractional Schrödinger equations. J. Comput. Phys. 2015, 293, 238–251. [Google Scholar] [CrossRef]
  23. Sousa, E.; Li, C. A weighted finite difference method for the fractional diffusion equation based on the riemann–liouville derivative. Appl. Numer. Math. 2015, 90, 22–379. [Google Scholar] [CrossRef]
  24. Lubich, C. Discretized fractional calculus. SIAM J. Math. Anal. 1986, 17, 704–719. [Google Scholar] [CrossRef]
  25. Tang, T. A finite difference scheme for partial integro-differential equations with a weakly singular kernel. Appl. Numer. Math. 1993, 11, 309–319. [Google Scholar] [CrossRef]
Figure 1. The imaginary and real parts of egenfunctions { P m ( ρ ) } m = 1 5 for Equation (10) at β = 1.1 .
Figure 1. The imaginary and real parts of egenfunctions { P m ( ρ ) } m = 1 5 for Equation (10) at β = 1.1 .
Mathematics 10 03160 g001
Figure 2. The analytical and numerical solution of (36) for α = 0.5 , β = 1.9 and T = 1 .
Figure 2. The analytical and numerical solution of (36) for α = 0.5 , β = 1.9 and T = 1 .
Mathematics 10 03160 g002
Figure 3. The analytical and numerical solution of (36) for different α = 0.2 , 0.5 , 0.9 , β = 1.9 and t = T = 1 .
Figure 3. The analytical and numerical solution of (36) for different α = 0.2 , 0.5 , 0.9 , β = 1.9 and t = T = 1 .
Mathematics 10 03160 g003
Table 1. The first five real and complex eigen values of Equation (10) at different β .
Table 1. The first five real and complex eigen values of Equation (10) at different β .
β = 1.1 1.4 1.6 1.95
λ m 3.86752 4.68735 5.60788 9.00579
  10.7157 ± 7.36402 I 16.2855 19.3929 35.6256
  14.8325 ± 15.3326 I 25.8219 37.7632 78.8307
  18.8242 ± 23.3452 I 44.8121 ± 10.6661 I 61.3084 138.351
  23.0029 ± 31.3436 I 68.9794 ± 25.9039 I 200.772 ± 7.49232 I 213.872
Table 2. The error and convergence rate when = 0.02 and β = 1.9 , T = 1 .
Table 2. The error and convergence rate when = 0.02 and β = 1.9 , T = 1 .
α = 0.2 α = 0.5 α = 0.7 α = 0.9
τ Λ ( τ , ) R t Λ ( τ , ) R t Λ ( τ , ) R t Λ ( τ , ) R t
1/20 5.8547 × 10 4 5.2961 × 10 4 4.9107 × 10 4 4.4985 × 10 4
1/40 2.9077 × 10 4 1.0097 2.5737 × 10 4 1.0410 2.4105 × 10 4 1.0265 2.1869 × 10 4 1.04053
1/80 1.5035 × 10 4 0.9515 1.3682 × 10 4 0.9115 1.1948 × 10 4 1.0125 1.0524 × 10 4 1.0552
1/160 7.9782 × 10 5 0.9142 6.6193 × 10 5 1.0475 5.8716 × 10 5 1.0250 4.9792 × 10 5 1.0797
1/320 3.8468 × 10 5 1.0526 3.1851 × 10 5 1.0553 2.9153 × 10 5 1.0101 2.4790 × 10 5 1.0060
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mahmoud, E.I.; Aleroev, T.S. Boundary Value Problem of Space-Time Fractional Advection Diffusion Equation. Mathematics 2022, 10, 3160. https://doi.org/10.3390/math10173160

AMA Style

Mahmoud EI, Aleroev TS. Boundary Value Problem of Space-Time Fractional Advection Diffusion Equation. Mathematics. 2022; 10(17):3160. https://doi.org/10.3390/math10173160

Chicago/Turabian Style

Mahmoud, Elsayed I., and Temirkhan S. Aleroev. 2022. "Boundary Value Problem of Space-Time Fractional Advection Diffusion Equation" Mathematics 10, no. 17: 3160. https://doi.org/10.3390/math10173160

APA Style

Mahmoud, E. I., & Aleroev, T. S. (2022). Boundary Value Problem of Space-Time Fractional Advection Diffusion Equation. Mathematics, 10(17), 3160. https://doi.org/10.3390/math10173160

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop