Next Article in Journal
Blind Color Image Watermarking Using Fan Beam Transform and QR Decomposition
Next Article in Special Issue
Noether Symmetries of a Generalized Coupled Lane-Emden-Klein-Gordon-Fock System with Central Symmetry
Previous Article in Journal
An EDAS Method for Multiple Attribute Group Decision-Making under Intuitionistic Fuzzy Environment and Its Application for Evaluating Green Building Energy-Saving Design Projects
Previous Article in Special Issue
A Sharp Oscillation Criterion for a Linear Differential Equation with Variable Delay
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Finite Difference Approximation Method for a Space Fractional Convection–Diffusion Equation with Variable Coefficients

by
Eyaya Fekadie Anley
1,2 and
Zhoushun Zheng
1,*
1
School of Mathematics and Statistics, Central South University, Changsha 410083, China
2
College of Natural and Computational Science, Department of Mathematics, Arba-Minch University, Arba-Minch 21, Ethiopia
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(3), 485; https://doi.org/10.3390/sym12030485
Submission received: 17 January 2020 / Revised: 22 February 2020 / Accepted: 1 March 2020 / Published: 23 March 2020

Abstract

:
Space non-integer order convection–diffusion descriptions are generalized form of integer order convection–diffusion problems expressing super diffusive and convective transport processes. In this article, we propose finite difference approximation for space fractional convection–diffusion model having space variable coefficients on the given bounded domain over time and space. It is shown that the Crank–Nicolson difference scheme based on the right shifted Grünwald–Letnikov difference formula is unconditionally stable and it is also of second order consistency both in temporal and spatial terms with extrapolation to the limit approach. Numerical experiments are tested to verify the efficiency of our theoretical analysis and confirm order of convergence.

1. Introduction

Fractional differential equations (FDE) have attracted the attention of many researchers and scientists due to their importance in different fields of study such as viscoelasticity, fluid mechanics, physics, biology, engineering, and flows in porous media (see [1,2,3,4,5,6] and the references cited therein). As different experiments and implementations have shown, non-integer space derivatives have been used to develop anomalous diffusion to which a particle spreads at a rate inconsistent with the integer Brownian motion problem in the direction of both time and space. When non-integer order is replaced by the second order derivative in a diffusion equation, it acts to enhance the process which we call super-diffusion [7,8,9,10,11,12]. Laboratory experiments and field-scale tracer dispersion breakthrough curves (BTCs) are suitable for exhibiting early time arrivals that are not captured by the integer order derivatives and these non-Fickian phenomena can be controlled by non-classical order convection–diffusion and dispersion equations (FCDE) as it was explained in [13]. To increase the number of applications, there should be significant interest in constructing numerical schemes to solve a well known space fractional convection–diffusion model that has space variable coefficients. In most cases, non-integer order differential problems have no exact solution, so various iterative and numerical approximations [3,9,14] must be pointed out in advance. In general, these kinds of approaches have become important in finding the approximate solutions of fractional differential equations, so extensive numerical methods have been developed for space fractional convection–diffusion equations such as the spectral method [15], finite volume method [16,17], finite difference method [2,9,14,18,19,20,21,22,23,24,25,26], finite element method [27,28,29,30] and collocation method [31,32].
When the discretization of domain over the region (which belongs to the geometry) is not complex, finite difference approximations are easier and faster than other methods (see [16,33] for further details) to get numerical solutions. In [34], the author used an unconditional stable difference method for time−space fractional convection–diffusion problems with space variable coefficients with first order convergence both in time and space. The Crank–Nicolson finite difference method for one-sided space fractional diffusion equations using an extrapolation method to get second order convergence was studied in [23]. In [9], the explicit and implicit finite difference methods are discussed for a one-sided space fractional convection–diffusion equation with first order convergence in both time and space. A first-order implicit finite difference discretization method for a two-sided space fractional diffusion equation (SFDE) is also applied in [10]. Recently, an unconditionally stable second order accurate difference method for a two-sided time−space fractional convection–diffusion equation was constructed in [35] using the weighted and Shifted Grünwald–Letnikov difference approximation. It is not suitable to apply the weighted combined with shifted Grünwald–Letnikov difference approximation for one-sided Riemann–Liouville fractional derivative to have second order accurate in space. To deal with such issues, it is important to develop a numerical scheme that leads to evaluate a one-sided space fractional convection–diffusion problem. Thus, the main focus of our study is to have temporal and spatial second order convergence estimates for one-sided space fractional convection–diffusion equations based on a stable finite difference method and using spatial extrapolation to the limit approach. The scheme has been treated using the Crank–Nicholson method with the novel Shifted Grünwald–Letnikov difference approximation and the algorithm has been examined both theoretically and experimentally.
Let us consider space-fractional convection–diffusion equation with variable coefficients:
u ( x , t ) t + c ( x ) u ( x , t ) x = d ( x ) α u ( x , t ) x α + p ( x , t ) , x ( L , R ) , t ( 0 , T ] , α ( 1 , 2 ] ;
with the given initial condition:
u ( x , 0 ) = g ( x ) , L x R ,
and homogeneous Dirichlet boundary conditions:
u ( L , t ) = 0 , u ( R , t ) = 0 , 0 t T ,
where c ( x ) , d ( x ) and g ( x ) are continuous functions on [ L , R ] and p ( x , t ) is continuous function on [ L , R ] × [ 0 , T ] . Here u ( x , t ) is the concentration, d ( x ) > 0 is the variable diffusion coefficient, c ( x ) > 0 is the fluid variable velocity which means the system is evolving in space due to a velocity field and p ( x , t ) is sink term so that the fluid transport is from left to right. For the case of integer order ( α = 2 ), Equation (1) gives to the classical convection–diffusion equation(CDE). In this study, we have only considered the fractional derivative case which describes a physical meaning in [36] and it involves only a left-sided fractional order derivative. We have assumed that this one-dimensional space fractional convection–diffusion problem has sufficiently smooth and unique enough solutions.
The structure of this paper is arranged as follows. In Section 2, we introduce some preliminary remarks, lemmas and definitions and we show the formulation of the new Crank–Nicolson with right Shifted Grünwald–Letnikov difference scheme in Section 3. In Section 4, we describe the unconditional stability using Gerschgorin Theorem and convergence order analysis of the scheme. In Section 5, numerical tests are implemented to show the relevance of our theoretical study and the conclusions are put in Section 6.

2. Preliminary Remarks

Definition 1.
The Riemann fractional derivative operator D α with order α is written as:
D α u ( x ) = 1 Γ ( r α ) d r d x r L x u ( t ) ( x t ) α r + 1 d t , α > 0
where r 1 < α < r , r N , t > 0 .
Definition 2.
The left hand side and the right hand side fractional order derivatives, respectively, in Equation (1) are the Riemann–Liouville fractional derivatives with order α which are given by:
( D + α u ) ( x ) = 1 Γ ( r α ) d r d x r L x ( x s ) r α 1 u ( s ) d s ( D α u ) ( x ) = ( 1 ) r Γ ( r α ) d r d x r x R ( s x ) r α 1 u ( s ) d s
for r 1 < α < r , x .
Definition 3
([3]). Let u be given on ℜ. The standard Grünwald–Letnikov estimate for 1 < α 2 with positive order α is defined by the formula,
D α u ( x , t ) 1 h α k = 0 N x ω k ( α ) u ( x k h , t ) ,
we also define the Grünwald–Letnikov difference operator as:
h α ( Δ h α u ) ( x , t ) k = 0 N x ω k ( α ) u ( x k h , t ) , h > 0 , x ,
where
ω k ( α ) = α ( α 1 ) ( α k + 1 ) k ! ,
is called Grünwald–Letnikov coefficient which is the Taylor series expansion ω ( z ) = ( 1 z ) α which is the generating function. We can expressed the coefficients by the following recursive relations.
ω 0 ( α ) = 1 , ω k ( α ) = ( 1 α + 1 k ) ω k 1 ( α ) , k = 1 , 2 , .
Lemma 1
([37]). Assume that 1 < α 2 , then Grünwald–Letnikov coefficients ω k ( α ) satisfy:
ω 0 ( α ) = 1 , ω 1 ( α ) = α < 0 , ω 2 ( α ) = α ( α 1 ) 2 > 0 1 ω 2 ( α ) ω 3 ( α ) 0 , k = 0 ω k ( α ) = 0 , k = 0 N x ω k ( α ) < 0 , N x 1 .
The Shifted Grünwald–Letnikov difference operator expression is suitable for our purpose because, it allows us to estimate D α u ( x ) , which is defined in Equation (2), numerically in an accurate way. According to [14], right shifted Grünwald–Letnikov difference operator with p shifts for α t h order Left R-L fractional derivative of u ( x , t ) , x [ L , R ] at x = x m can be expressed as:
D α u ( x , t ) 1 h α k = 0 x m L h + p ω k ( α ) u ( x ( k p ) h , t )
where
x m = L + m h , h = R L N x , m = 0 , 1 , 2 , N x .
Lemma 2
([38,39]). Let u C 2 n ( ) that has a finite degree of smoothness with ( D + α u ) ( x ) which is approximated by h α Δ h α u ( x ) possesses an asymptotic expansion in integer powers of the step-length h, then an expansion in even powers of h for the Shifted operator can be written in the form:
Δ h , p α u ( x ) = j = 0 1 j j α u x + α h 2 j h , h > 0 .
Lemma 3
([39]). Let u C n + 3 ( ) all derivative of u up to the order n + 4 belong to L 1 ( ) . Then the Fourier transform of the Grünwald–Letnikov difference operator defined in Equation (5), is
ϕ ^ ( x ) = ϕ ( t ) e i x t d t .
Theorem 1.
Let u C 2 n + 3 ( ) with all derivatives of u up to order 2 n + 3 belong to L 1 ( ) . For p 0 define the shifted Grünwald–Letnikov operator:
( Δ h , p α ) u ( x ) = k = 0 ω k ( α ) u x k p h ,
with ω k ( α ) = ( 1 ) 2 k α a 2 k = α a 2 k . Then, if L = in Equation (2), for any computable coefficient a 2 k , which is independent of h , u and x, we have
h α Δ h , p α u ( x ) = D + α u ( x ) + k = 1 n 1 b 2 k D + α + 2 k u ( x ) h 2 k + O ( h 2 n )
uniformly in x .
Proof of Theorem 1.
We closely follow the result described in [9,10] for the unshifted Grünwald–Letnikov formula and also in [23] for the shifted Grünwald–Letnikov formula. We can see that with the Riemann-Lebesgue lemma, the assumptions on u indicates for real positive constant C 1 and from the condition which is imposed on u, we have
u ˜ ( t ) C 1 ( 1 + | t | ) 2 n 3 .
From Lemma 3 for all t the Fourier transform for u ( x ) of the Grünwald–Letnikov approximation is
u ˜ ( t ) = u ( x ) e i x t d x .
From the definition of Fourier transform, we have observed that for a constant a , we have:
F ( [ u ( x a ) ] ) ( t ) = e i a t u ˜ ( t ) .
The function
1 e z z α e z p = ω α , p ( z ) ,
have the Taylor expansion
ω α , p ( z ) = k = 0 a 2 k z 2 k ,
where a 2 k = ( 1 ) 2 k α a 2 k = α a 2 k , converges absolutely for z 1 since the function ω α , p ( z ) is bounded on ℜ. The shifted Grünwald difference approximation ( Δ h , p ) u ( x ) L 1 ( ) .
Thus, we have
F ( h α Δ h , p α u ) ( t ) = h α e i t p h k = 0 2 k α e i k t h u ˜ ( t ) = h α e i t p h 1 e i t p h α u ˜ ( t ) = ( i t ) α 1 e i t h i t h α e i t p h u ˜ ( t ) = ( i t ) α ω α , p ( i t h ) u ˜ ( t )
since ω α , p ( i t h ) is analytic around the origin, we express it as an even power expansions
ω α , p ( z ) = k = 0 a 2 k z 2 k
which absolutely convergent for all | z | R . For this a bounded function ω α , p ( z ) on , there exist a real positive constant C 2 which satisfy:
1 e i x i x α k = 0 n 1 a 2 k i x 2 k C 2 x 2 n
is bounded uniformly in x . For any value | x | R , we have
( ω α , p ( i x ) k = 0 n 1 a 2 k ( i x ) 2 k = k = n a 2 k ( i x ) 2 k x 2 n k = n a 2 k α x 2 ( k n ) C 3 x 2 n
which is bounded on . For the other case | x | > R also, we have
ω α , p i x = 1 e i x i x α e i p x 2 α R α < C 4 x 2 n
where C 4 = 2 α R α + 2 n < and also
k = 0 n 1 a 2 k ( i x ) 2 k x 2 n k = 0 n 1 | ( a 2 k α ) | x 2 ( k n ) C 5 x 2 n
with C 5 = k = 0 n 1 a 2 k α R 2 k 2 n < . Now, we set that
C 2 = max k = 0 a 2 k R 2 k 2 n , 2 α R α + 2 n + k = 0 n 1 a 2 k R 2 k 2 n
since
k = 0 a 2 k R 2 k 2 n = k = 0 n 1 a 2 k R 2 k 2 n + k = n a 2 k R 2 k 2 n
C 2 = 2 α R α + 2 n + k = 0 n 1 a 2 k R 2 k 2 n
Then, this implies that Equation (15) holds for all x . From Equation (17), we can write
F ( h α Δ h , p α u ) ( t ) = k = 0 n 1 a 2 k i t α + 2 k h 2 k u ˜ ( t ) + φ ˜ ( t , h )
where
φ ˜ ( k , h ) = i t α ω α , p i t h k = 0 n 1 a 2 k i t h 2 k u ˜ ( t )
since
i t α + 2 k u ˜ ( t ) = D + α + 2 k u ˜ ( t ) .
Therefore, we have
( i t ) α + 2 k u ˜ ( t ) L 1 ( ) .
Moreover, we see that
φ ˜ ( t , h ) L 1 ( ) ,
and with the conditions imposed on u, we can say that ( 1 + x 2 n + 3 u ˜ ( t ) is bounded on .
Thus, t 2 α 3 | u ˜ ( t ) | L 1 ( R ) . This implies that,
φ ˜ ( t , h ) C h 2 n 1 + t 2 α 3
for k with C = C 1 C 2 . Therefore using the Fourier inversion transform, we have
h α Δ h , p α u ( x ) = D + α u ( x ) + k = 1 n 1 a 2 k D + α + 2 k u ( x ) h 2 k + φ ( x , h ) ,
where
φ ( x , h ) = C R e i t x φ ˜ ( t , h ) d t C R φ ˜ ( t , h ) d t C h 2 n .
At last, we have
h α Δ h , p α u ( x ) = D + α u ( x ) + k = 1 n 1 a 2 k ( D + α + 2 k u ) ( x ) h 2 k + O ( h 2 n ) .
 □
Remark 1.
From Equation (10), it can be seen that for p = α / 2 , the error takes its minimum value and a second order convergence is achieved. We need the grid points x m ( k p ) h to find an optimal positive integer p that makes p α / 2 is minimum. It is numerically proved in [3] that for the value 0 < α 1 , p = 0 is acceptable; while for 1 < α 2 , p = 1 is optimal.
Remark 2.
Theorem 1 is the base of Extrapolation to the limit. Therefore one can apply it the Shifted Grünwald–Letnikov difference operator to obtain the convergence rate with arbitrary high order h k , k = 1 , 2 , 3 , , n such that
h α ( q α Δ q h , p α u ) ( x ) q ( Δ h , p α u ) ( x ) 1 q , 0 < q < 1
(q is fixed) converges to ( D + α u ) ( x ) + O ( h 2 ).

3. Problem Formulation of the Scheme

Consider the following one-dimensional space fractional convection–diffusion problem:
u ( x , t ) t = c ( x ) u ( x , t ) x + d ( x ) α u ( x , t ) x α + p ( x , t ) , ( x , t ) ( L , R ) × ( 0 , T ] u ( x , 0 ) = g ( x ) , x [ L , R ] u ( L , t ) = 0 , u ( R , t ) = 0 , t [ 0 , T ]
which is based on shifted Grünwald–Letnikov difference method with 1 < α 2 on a finite domain L < x < R .

Crank–Nicolson Scheme for Time and Shifted Grünwald Difference Scheme for Space Discretization

We partition the finite interval [ L , R ] with a uniform mesh in the space size step h = ( R L ) / N x and the time step τ = T / N t , in which N x , N t are non-negative integers and the set of grid size points is symbolized by x m = m h and t n = n τ for 0 m N x , 0 n N t . Set t n + 1 / 2 = ( t n + 1 + t n ) / 2 with 0 n N t 1 .
We use the following notations:
u m n = u ( x m , t n ) , p m n + 1 / 2 = p ( x m , t n + 1 / 2 ) , δ t u m n = u m n + 1 u m n τ , c m = c ( x m ) , d m = d ( x m ) .
Applying the C-N technique for the time discretization of Equation (20) gives to
δ t u m n = c m 4 h u m + 1 n + 1 u m 1 n + 1 + u m + 1 n u m 1 n + d m 2 h α z = 0 1 k = 0 N x 1 ω k ( α ) u m k + 1 n + z = p m n + 1 / 2 + O ( τ 2 ) .
In space discretization we have used the central finite difference method for the convection term and the Shifted Grünwald–Letnikov operator for the space fractional derivative with the approach of spatial Extrapolation to the limit, respectively.
See the full discretization of the scheme:
u m n + 1 u m n τ = c m u m + 1 n u m 1 n + u m + 1 n + 1 u m 1 n + 1 4 h + d m 2 h α z = 0 1 k = 0 m + 1 ω k ( α ) u m k + 1 n + z + p m n + p m n + 1 2 .
Multiplying Equation (22) by τ the discretization equation, we have
u m n + 1 u m n = c m τ 4 h ( u m + 1 n u m 1 n + u m + 1 n + 1 u m 1 n + 1 ) + d m τ 2 h α z = 0 1 k = 0 m + 1 ω k ( α ) u m k + 1 n + z + τ p m n + 1 / 2
The above equation is used to predict the values of u ( x , t ) at time n + 1 , so all the values of u at time n are assumed to be known. For simplification
μ m = c m τ h , η m = d m τ h α , then we have
1 η m 2 ω 1 α u m n + 1 + μ m 2 η m 2 ω 2 α u m 1 n + 1 + μ m 2 η m 2 ω 0 α u m + 1 n + 1 η m 2 k = 3 m + 1 ω k α u m k + 1 n + 1 = 1 + η m 2 ω 1 α u m n + μ m 2 + η m 2 ω 2 α u m 1 n + η m 2 ω 0 α + μ m 2 u m + 1 n + η m 2 k = 3 m + 1 ω k ( α ) u m k + 1 n + τ p m n + 1 2 .
Both the convection and diffusion variable coefficients are ( N x 1 ) × ( N x 1 ) diagonal matrices which are defined by
μ m = τ 2 h d i a g C 1 , C 2 , C 3 , C N x 1 , η m = τ h α d i a g d 1 , d 2 , d 3 , d N x 1 .
These discretization together with Dirichlet boundary conditions which results in a linear system of equations for which the coefficient matrix is the sum of lower triangular and upper-diagonal matrices. The above discretization can be re-arranged to yield:
1 η m 2 ω 1 α u m n + 1 + ( μ m 2 η m 2 ω 2 α ) u m 1 n + 1 + μ m 2 η m 2 ω 0 α ) u m + 1 n + 1 η m 2 ( k = 3 m + 1 ω k α u m k + 1 n + 1 = ( 1 + η m 2 ω 1 α ) u m n + μ m 2 + η m 2 ω 2 α u m 1 n + ( η m 2 ω 0 α + μ m 2 ) u m + 1 n + η m 2 ( k = 3 m + 1 ω k α u m k + 1 n ) + τ ( P m n + 1 2 ) .
Denoting U m n as the numerical approximation of u m n , we can construct the C-N scheme for Equation (20)
1 η m 2 ω 1 α U m n + 1 + μ m 2 η m 2 ω 2 α U m 1 n + 1 + μ m 2 η m 2 ω 0 α U m + 1 n + 1 η m 2 k = 3 m + 1 ω k α U m k + 1 n + 1 = 1 + η m 2 ω 1 α U m n + μ m 2 + η m 2 ω 2 α U m 1 n + η m 2 ω 0 α + μ m 2 U m + 1 n + η m 2 k = 3 m + 1 ω k α U m k + 1 n + τ ( P m n + 1 2 ) .
I is the N x 1 × N t 1 identity matrix with A m , n as the matrix coefficients. These coefficients, for m = 1 , 2 , 3 , , N x 1 , n = 1 , 2 , , N t 1 are given by:
A m , n = 0 , n m + 2 μ m 2 η m 2 ω 0 ( α ) , n = m + 1 ( 1 η m 2 ω 1 ( α ) ) , n = m ( η m 2 ω 2 ( α ) μ m 2 ) , n = m 1 η m 2 ω m n + 1 ( α ) n m 1 .
The finite difference scheme (24) and (26) defines a linear system of equations as
( I + A ) U n + 1 = ( I A ) U n + τ ( p m n + 1 2 ) U n + 1 = [ u 1 n + 1 , u 2 n + 1 , , u N x 1 n + 1 ] U n + τ P m n + 1 2 = [ 0 , τ p 1 n + 1 2 , τ p 2 n + 1 2 , , τ p N x 1 n + 1 2 + ( η N x 1 2 + μ N x 1 2 ) , 0 ] .
Theorem 2.
Suppose that 1 < α 2 , the coefficient matrix defined in Equations (24)–(27), then the diagonal matrix and the coefficient matrix satisfy:
A m , m > n = 0 , m 1 N x 1 | A m , n | , m = 1 , 2 , 3 , , N x 1 .
Proof of Theorem 2.
As we have seen from the coefficient matrix defined in Equation (27),
A m , m + 1 = μ m 2 η m 2 ω 0 ( α ) = μ m 2 η m 2 < 0
A m , m 1 = η m 2 ω 2 ( α ) μ m 2 = η m 2 ( α 2 α 2 ) , but from Lemma 1, α 2 α 2 > 0 for 1 < α 2 mean that η m 2 ( α 2 α 2 ) < 0 .
When n < m 1 , we have, η m 2 ω m n + 1 ( α ) < 0 and when n = m , A m , m = 1 η m 2 ω 1 ( α ) = 1 + η m 2 α > 0 . This implies that n = 0 , m 1 N x 1 | A m , n | < A m , m .
Therefore, the diagonal matrix is strictly dominant. □

4. Theoretical Analysis of Finite Difference Scheme

In general for analyzing convergence and stability, we consider the following description.
Let χ h = ν : ν = ν m : x m = m h m = 0 N x , ν 0 = ν N x = 0 be the grid function.
For any ν = ν m χ h , we define our point-wise maximum norm as
| | ν | | = m a x 1 m N x | ν m | ,
and the discrete L 2 -norm
ν = h m = 1 N x 1 ν m 2 .

4.1. Boundedness of the Fractional Scheme

The Classical Crank–Nicolson scheme combines the stability of an implicit finite difference method with its accuracy which produce second order convergence in both space and time.
Theorem 3.
Crank–Nicolson scheme for solving space fractional convection–diffusion equations given by the following problem:
u ( x , t ) t + c ( x ) u ( x , t ) x = d ( x ) α u ( x , t ) x α + p ( x , t ) .
which is based on shifted Grünwald–Letnikov difference approximation scheme is bounded for 1 < α 2 .
Proof of Theorem 3.
Consider C-N scheme for the space-fractional convection–diffusion problem for 1 < α 2
u m n + 1 u m n τ = c m u m n u m 1 n + u m n + 1 u m 1 n + 1 2 h + d m 2 h α j = 0 1 k = 0 m + 1 ω k ( α ) u m k + 1 n + j + p m n + p m n + 1 2 .
Here, we have shown the convergence and boundedness of the scheme by taking the smaller time-step in terms of Lax–Richtmyer stability analysis that uses a weaker bound (see [40]). Our matrix A has an eigenvalues of λ that have positive real parts, and, we also have found a strictly dominant matrix. These eigenvalues which are centered in the disks at each diagonal entries as:
A m , m = ( 1 η m 2 ω 1 α ) = 1 + α η m 2 .
with μ m = c m τ h , η m = τ d m h α . From the Gerschgorin Theorem in [41], the radius of the matrix can be expressed as
n = 0 , m 1 N x A m , n 2 2 = η m 2 μ m 2 n = 0 m + 1 ω m n + 1 ( α ) 2 η m 2 μ m 2 2 n = 0 m + 1 ω m n + 1 ( α ) 2 .
Since from the Grünwald coefficients we have ω m n + 1 ( α ) ω 1 ( α ) and ω 1 ( α ) = α , we have that:
n = 0 , m 1 N x ( A m , n ) 2 2 n = 0 , m 1 N x ( A m , m ) 2 2 A m , m 2 2 ( η m 2 μ m 2 ) 2 2 ω 1 ( α ) 2 2 1 + η m 2 α 2 2 .
For a bounded ratio of time-step τ and space-step h with n τ T , we have
A m , m n 2 1 + η m 2 α n / 2 .
From the relation of Parseval’s Theorem, [40]
A m , m 2 1 + η m 2 α n / 2 e α T / 2 .
which shows that the scheme is bounded. □

4.2. Stability Analysis

Theorem 4.
Let U m n be the numerical approximation of the exact solution u m n , then the C-N finite difference scheme (28) is unconditionally stable.
Proof of Theorem 4.
Consider the matrix coefficient of the difference approximation for the problem (20) can be written as described above
( I + A ) U n + 1 = ( I A ) U n + τ p m n + 1 / 2 .
Let e n = e 1 n , e 2 n , e 3 n , , e N x 1 n , and take the relation between the error e n + 1 in U n + 1 and the error e n in U n which is given by the linear system
e n + 1 = ( I + A ) 1 ( I A ) e n .
First of all, we must show that the (non-real valued) eigenvalues of the coefficient matrices A have positive real parts. For ω 1 ( α ) = α with fractional order 1 < α < 2 and k 1 ; we have ω k ( α ) > 0 . In addition to this, ω 1 α = α k = 0 , k 1 N ω k α for the value N > 1 . As stated in Gerschgorin Theorem ([41], pp. 136–139), the eigenvalues of the given matrix A are inside the disks centered at each diagonal entry.
A m , m = ( 1 η m 2 ω 1 ( α ) ) = 1 + η m 2 α > 0 ,
with radius
r m = n = 0 , m 1 N x | A m , n | = η m 2 n = 0 m + 1 ω m n + 1 ( α ) < ( 1 + η m 2 ) .
These Gerschgorin disks are belong to the right half of the complex plane. Thus, the eigenvalue of the coefficient matrix A has positive real part which implies that A has an eigenvalue λ if and only if ( I A ) has an eigenvalue ( 1 λ ) if and only if ( I + A ) 1 ( I A ) has an eigenvalue 1 λ 1 + λ . From the first part of this sentence, we have seen that all the eigenvalues of the matrix given by ( I + A ) have a radius larger than unity which implies the matrix is invertible. Now we can see from the above description the real part of λ is non-negative which we can conclude that ( 1 λ ) ( 1 + λ ) < 1 .
Thus, the spectral radius of the system matrix ( I + A ) 1 ( I A ) is strictly less than unity which implies that the difference scheme is unconditionally stable. □

4.3. Convergence Analysis

First of all we have given the Truncation error of the C-N scheme. It is obvious to conclude that:
u ( x m , t n + 1 ) u ( x m , t n ) τ = u ( x , t ) t n + 1 / 2 + O ( τ 2 ) . c ( x ) u ( x , t ) x + d ( x ) α u ( x , t ) x α m n + 1 / 2 = 1 2 c m u ( x m , t n + 1 ) x + d m α u ( x m , t n + 1 x α + 1 2 c m u ( x m , t n ) x + d m α u ( x m , t n x α + O ( τ 2 ) .
c ( x m ) u ( x , t ) x u ( x m + 1 , t n + 1 ) u ( x m 1 , t n + 1 ) 2 h + O ( h 2 ) .
From the above Extrapolation to the limit Theorem for n = 1 , we got
α u ( x , t ) x α k = 0 m + 1 g k ( α ) u m k + 1 + O ( h 2 ) .
Therefore the local truncation error of (20) is given by T m n + 1 = O ( τ 2 + τ h )
Theorem 5.
Let u m n be the exact solution of problem (20), and U m n be the solution of the finite difference scheme (26), then for all 1 n N t , we have the estimate
u m n U m n c ( τ 2 + h )
where u m n U m n = m a x 1 m N x u m n U m n , c is a non-negative constant independent of h and τ with | | . | | stands for the discrete L 2 -norm.
Proof of Theorem 5.
Denote e n = u m n U m n where e n = ( e 1 n , e 2 n , , e N x 1 n ) . We have e 0 = 0 , we have from Equations (26) and (27) if n = 0 ,
R m 1 = μ m 2 η m 2 ω 0 ( α ) e m 1 1 + 1 + η m 2 α e 1 m + μ m 2 η m 2 ω 2 ( α ) e m + 1 1 η m 2 k = 3 N x ω k ( α ) e m n + 1 1 .
if n > 0 ,
R m n + 1 = μ m 2 η m 2 ω 0 ( α ) e m 1 n + 1 + 1 + η m 2 α e n + 1 m + μ m 2 η m 2 ω 2 ( α ) e m + 1 n + 1 η m 2 k = 3 N x ω k ( α ) e m n + 1 n + 1 .
where R m n + 1 c ( τ 2 + h ) , m = 1 , 2 , , N x 1 , n = 1 , 2 , 3 , , N t 1 , c is non-negative constant independent of h and τ .
We can use the mathematical induction to prove the Theorem. Let n = 1 and assume | e j | = m a x 1 m N x 1 | e m 1 | , we have the following expression.
| | e 1 | | = e j 1 μ j 2 η j 2 ω 0 ( α ) e j 1 1 + 1 + η j 2 α e 1 j + μ j 2 η j 2 ω 2 ( α ) e j + 1 1 η j 2 k = 3 N x ω k ( α ) e j n + 1 1 μ j 2 η j 2 ω 0 ( α ) e j 1 1 + 1 + η j 2 α e 1 j + μ j 2 η j 2 ω 2 ( α ) e j + 1 1 η j 2 k = 3 N x ω k ( α ) e j n + 1 1 = R j 1 c ( τ 2 + h ) .
Suppose that if n r , | | e r | | c ( τ 2 + h 2 ) hold and assume n = r + 1 , let e j r + 1 = m a x 1 m N x 1 e m r + 1 , notice that from Lemma 1, we have k = 0 N x ω k ( α ) < 0 , m = 1 , 2 , , N x . Therefore,
e r + 1 = e j r + 1 μ j 2 η j 2 ω 0 ( α ) e j 1 r + 1 + 1 + η j 2 α e r + 1 j + μ j 2 η j 2 ω 2 ( α ) e j + 1 r + 1 η j 2 k = 3 N x ω k ( α ) e j n + 1 r + 1 μ j 2 η j 2 ω 0 ( α ) e j 1 r + 1 + 1 + η j 2 α e r + 1 j + μ j 2 η j 2 ω 2 ( α ) e j + 1 r + 1 η j 2 k = 3 N x ω k ( α ) e j n + 1 r + 1 = R j r + 1 c ( τ 2 + h )
which completes the proof. □
Remark 3.
The Crank–Nicolson scheme, for classical convection–diffusion equation, provides stable C-N finite difference method that is second order convergence in time and space. Also a study based on C-N finite difference method with the spatial extrapolation to the limit method, see Theorem 1, is used to get temporal and spatial second order for one-sided SFCDEs with space variable coefficients.

5. Numerical Tests

Problem test 1
1. Consider the space-fractional diffusion type of problem:
u ( x , t ) t = d ( x ) α u ( x , t ) x α + p ( x , t )
with initial condition
u ( x , 0 ) = ( x 2 x 3 ) ; 0 x 1
homogeneous Dirichlet boundary condition
u ( 0 , t ) = 0 = u ( 1 , t )
with variable diffusion coefficient,
d ( x ) = Γ ( 1.2 ) x α ,
and source term
p ( x , t ) = ( 6 x 3 3 x 2 ) e t
The exact solution is
u ( x , t ) = ( x 2 x 3 ) e t
All numerical experiments are implemented using Theorem 1 and C-N scheme with the space domain, 0 < x < 1 and time domain, 0 < t < T . Figure 1 shows the maximum error produced by C-N scheme for large enough time domain and numerical solution is close enough to the exact solution using C-N scheme with α = 1.5 in Figure 2. The maximum error and second order convergence for the fractional diffusion and fractional convection–diffusion equation with variable coefficients are given in Table 1, Table 2 and Table 3.
Problem test 2
2. Consider the space-fractional convection–diffusion type of equation with variable coefficients:
u ( x , t ) t + c ( x ) u ( x , t ) x = d ( x ) α u ( x , t ) x α + p ( x , t )
with initial condition
u ( x , 0 ) = ( x α x ) ; 0 x 1
homogeneous Dirichlet boundary condition
u ( 0 , t ) = 0 = u ( 1 , t )
with variable convection–diffusion coefficients respectively,
c ( x ) = x 1 5 , d ( x ) = x 1 100 ,
and source term
p ( x , t ) = e 2 t ( 2 ( x x α ) Γ ( α ) + Γ ( α + 1 ) Γ ( α ) x α 1 1 )
The exact solution is
u ( x , t ) = e 2 t ( x α x )
Figure 3 and Figure 4 show the numerical and exact solutions for fractional diffusion and fractional convection–diffusion problems with large enough time domain in example 1 and 2, respectively.The exact and numerical solution of fractional convection–diffusion equation by C-N scheme is also given in Figure 5. In Table 4, the maximum error and first order convergence in space is obtained using C-N scheme without extrapolation to the limit approach by fixing the time step.
Problem Test 3
3. Consider the space-fractional convection–diffusion type of equation with variable coefficients:
u ( x , t ) t + c ( x ) u ( x , t ) x = d ( x ) α u ( x , t ) x α + p ( x , t )
with initial condition
u ( x , 0 ) = x 2 ( 1 x )
homogeneous Dirichlet boundary condition
u ( 0 , t ) = 0 = u ( 1 , t )
with variable convection–diffusion coefficients respectively,
c ( x ) = x 0.6 , d ( x ) = Γ ( 2.8 ) x 3 / 4
and the forcing function
p ( x , t ) = 2 x 2 ( 1 x ) t 1.3 / Γ ( 2.3 ) + 0.3 x 1.8 e t
The exact solution is
u ( x , t ) = x 2 ( 1 x ) e t
Problem test 4
4. Consider the space-fractional convection–diffusion equation with variable coefficients:
u ( x , t ) t + c ( x ) u ( x , t ) x = d ( x ) α u ( x , t ) x α + p ( x , t )
with initial condition
u ( x , 0 ) = x α ( 1 x )
homogeneous Dirichlet boundary condition
u ( 0 , t ) = 0 = u ( 1 , t )
with variable convection–diffusion coefficients respectively,
c ( x ) = x 3 / 5 , d ( x ) = x 3 / 4
and the forcing function
p ( x , t ) = 2 x α ( 1 x ) t 1.3 / Γ ( 2.3 ) + 0.3 x 1.8 e t
The exact solution is
u ( x , t ) = x α ( 1 x ) e t
Problem test 4 is experimented with the grid size reduction extrapolation approach stated in [23]. We have smooth enough numerical and exact solutions by using C-N scheme in Figure 6, and Table 5 shows the maximum error with the error rate is given for space fractional convection–diffusion problem with a grid size reduction extrapolation method.

6. Conclusions

The one dimension space fractional diffusion and fractional convection–diffusion problem with space variable coefficients is solved by the fractional C-N scheme based on the Extrapolation to the limit approach of right shifted Grünwald–Letnikov approximation. The fractional C-N method, for the fractional diffusion problem and fractional convection–diffusion equation with space variable coefficients, is consistent and unconditionally stable with second order convergence. Numerical examples confirmed that the C-N method is suitable for the space fractional convection–diffusion problem even for a large value of time domain.

Author Contributions

The authors contributed equally to the writing and approved the final manuscript of this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the National Key Research (Grant No. 2017YFB0305601) and Development Program of China (Grant No. 2017YFB0701700).

Acknowledgments

The authors would like to thank the editor and the anonymous reviewers for their helpful comments for revising the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Diethelm, K. The Analysis of Fractional Differential Equations; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  2. Fomin, S.; Chugunov, V.; Hashida, T. Application of fractional differential equations for modeling the anomalous diffusion of contaminant from fracture into porous rock matrix with bordering alteration zone. Trans. Porous Med. 2010, 81, 187–205. [Google Scholar] [CrossRef]
  3. Guo, B.; Pu, X.; Huang, F. Fractional Partial Differential Equations and Their Numerical Solutions; World Scientific: Singapore, 2015. [Google Scholar]
  4. Podlubny, I. Fractional Differential Equations Academic; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  5. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: New York, NY, USA, 1998; Volume 198. [Google Scholar]
  6. Salman, W.; Gavriilidis, A.; Angeli, P. A model for predicting axial mixing during gas–liquid Taylor flow in microchannels at low Bodenstein number. Chem. Eng. J. 2004, 101, 391–396. [Google Scholar] [CrossRef]
  7. Berkowitz, B.; Cortis, A.; Dror, I.; Scher, H. Laboratory experiments on dispersive transport across interfaces. Water Resour. Res. 2009, 45. [Google Scholar] [CrossRef] [Green Version]
  8. Cortis, A.; Berkowitz, B. Computing “anomalous” contaminant transport in porous media. Ground Water 2005, 43, 947–950. [Google Scholar] [CrossRef]
  9. Meerschaert, M.M.; Tadjeran, C. Finite difference approximations for fractional advection–dispersion flow equations. J. Comput. Appl. Math. 2004, 172, 65–77. [Google Scholar] [CrossRef] [Green Version]
  10. Meerschaert, M.M.; Tadjeran, C. Finite difference approximations for two-sided space fractional partial differential equations. Appl. Numer. Math. 2006, 56, 80–90. [Google Scholar] [CrossRef]
  11. Ren, J.; Sun, Z.Z.; Zhao, X. Compact difference scheme for the fractional sub-diffusion equation with Neumann boundary conditions. J. Comput. Phys. 2013, 232, 456–467. [Google Scholar] [CrossRef]
  12. Wang, K.; Wang, H. A fast characteristic finite difference method for fractional advection–diffusion equations. Adv. Water. Resour. 2011, 34, 810–816. [Google Scholar] [CrossRef]
  13. Singh, J.; Swroop, R.; Kumar, D. A computational approach for fractional convection–diffusion equation via integral transforms. Ain Shams Eng. J. 2016, 9, 1019–1028. [Google Scholar] [CrossRef] [Green Version]
  14. Liu, F.; Zhuang, P.; Burrage, K. Numerical methods and analysis for a class of fractional advection–dispersion models. Comput. Math. Appl. 2012, 64, 2990–3007. [Google Scholar] [CrossRef] [Green Version]
  15. Tian, W.; Deng, W.; Wu, Y. Polynomial spectral collocation method for space fractional advection–diffusion equation. Numer. Methods Partial Differ. Equ. 2014, 30, 514–535. [Google Scholar] [CrossRef] [Green Version]
  16. Hejazi, H.; Moroney, T.; Liu, F. Comparison of finite difference and finite volume methods for solving the space fractional advection–dispersion equation with variable coefficients. ANZIAM J. 2013, 54, 557–573. [Google Scholar] [CrossRef] [Green Version]
  17. Liu, F.; Zhuang, P.; Turner, I.; Burrage, K.; Anh, V. A new fractional finite volume method for solving the fractional diffusion equation. Appl. Math. Model. 2014, 38, 15–16. [Google Scholar] [CrossRef]
  18. Li, C.; Zeng, F. Numerical Methods for Fractional Calculus; Chapman and Hall/CRC: Boca Raton, FL, USA, 2015. [Google Scholar]
  19. Li, C.; Zeng, F. Finite difference methods for fractional differential equations. Int. J. Bifurcat. Chaos 2012, 22, 1230014. [Google Scholar] [CrossRef]
  20. Luchko, Y. Multi-dimensional fractional wave equation and some properties of its fundamental solution. Ind. Math. 2014, 6. [Google Scholar] [CrossRef]
  21. Lynch, V.E.; Carreras, B.A.; del-Castillo-Negrete, D.; Ferreira-Mejias, K.M.; Hicks, H.R. Numerical methods for the solution of partial differential equations of fractional order. J. Comput. Phys. 2003, 192, 406–421. [Google Scholar] [CrossRef]
  22. Sweilam, N.H.; Khader, M.M.; Mahdy, A.M.S. Crank–Nicolson finite difference method for solving time-fractional diffusion equation. J. Fract. Calc. Appl. 2012, 2, 1–9. [Google Scholar]
  23. Tadjeran, C.; Meerschaert, M.M.; Scheffler, H.P. A second-order accurate numerical approximation for the fractional diffusion equation. J. Comput. Phys. 2006, 201, 205–213. [Google Scholar] [CrossRef]
  24. Wang, Y.M. A compact finite difference method for a class of time fractional convection–diffusion-wave equations with variable coefficients. Num. Algor. 2015, 70, 625–651. [Google Scholar] [CrossRef]
  25. Wang, Y.M.; Wang, T. A compact Alternative difference implicit method and its extrapolation for time fractional sub-diffusion equations with nonhomogeneous Neumann boundary conditions. Comput. Math. Appl. 2018, 75, 721–739. [Google Scholar] [CrossRef]
  26. Zhou, F.; Xu, X. The third kind Chebyshev wavelets collocation method for solving the time-fractional convection–diffusion equations with variable coefficients. Appl. Math. Comput. 2016, 280, 11–29. [Google Scholar] [CrossRef]
  27. Jin, B.; Lazarov, R.; Zhou, Z. A Petrov–Galerkin finite element method for fractional convection–diffusion equations. SIAM J. Numer. Anal. 2016, 54, 481–503. [Google Scholar] [CrossRef] [Green Version]
  28. Gao, F.; Yuan, Y.; Du, N. An upwind finite volume element method for nonlinear convection–diffusion problem. AJCM 2011, 1, 264. [Google Scholar] [CrossRef] [Green Version]
  29. Aboelenen, T. A direct discontinuous Galerkin method for fractional convection–diffusion and Schrödinger-type equations. Eur. Phys. J. Plus 2018, 133, 316. [Google Scholar] [CrossRef]
  30. Xu, Q.; Hesthaven, J.S. Discontinuous Galerkin method for fractional convection–diffusion equations. SIAM J. Numer. Anal. 2014, 52, 405–423. [Google Scholar] [CrossRef] [Green Version]
  31. Bhrawy, A.H.; Baleanu, D. A spectral Legendre−Gauss–Lobatto collocation method for a space-fractional advection–diffusion equations with variable coefficients. Rep. Math. Phys. 2013, 72, 219–233. [Google Scholar] [CrossRef]
  32. Saadatmandi, A.; Dehghan, M.; Azizi, M.R. The Sinc-Legendre collocation method for a class of fractional convection–diffusion equations with variable coefficients. Commun. Nonlinear. Sci. Numer. Simulat. 2012, 17, 4125–4136. [Google Scholar] [CrossRef]
  33. Pang, G.; Chen, W.; Sze, K.Y. A comparative study of finite element and finite difference methods for two-dimensional space-fractional advection–dispersion equation. Adv. Appl. Math. Mech. 2016, 8, 166–186. [Google Scholar] [CrossRef]
  34. Zhang, Y. A finite difference method for fractional partial differential equation. Appl. Math. Comput. 2009, 2015, 524–529. [Google Scholar] [CrossRef]
  35. Gu, X.M.; Huang, T.Z.; Ji, C.C. Carpentieri B, Alikhanov AA, Fast iterative method with a second-order implicit difference scheme for time-space fractional convection–diffusion equation. J. Sci. Comput. 2017, 72, 957–985. [Google Scholar] [CrossRef]
  36. Benson, D.A.; Wheatcraft, S.; Meerschaert, M.M. Application of a fractional advection–dispersion equation. Water Resour. Res. 2000, 36, 1403–1412. [Google Scholar] [CrossRef] [Green Version]
  37. Tian, W.; Zhou, H.; Deng, W. A class of second order difference approximations for solving space fractional diffusion equations. Math. Comp. 2015, 84, 1703–1727. [Google Scholar] [CrossRef] [Green Version]
  38. Henrici, P. Elements of Numerical Analysis; John Wiley and Sons: New York, NY, USA, 1984. [Google Scholar]
  39. Tuan, V.K.; Gorenflo, R. Extrapolation to the Limit for Numerical Fractional Differentiation. Z. Angew. Math. Mech. 1995, 75, 646–648. [Google Scholar] [CrossRef]
  40. LeVeque, R.J. Finite Difference Methods for Ordinary and Partial Differential Equations; SIAM: Philadelphia, PA, USA, 2007; Volume 98. [Google Scholar]
  41. Isaacson, E.; Keller, H.B. Analysis of Numerical Methods; Courier Corporation: Chelmsford, MA, USA, 2012. [Google Scholar]
Figure 1. The Maximum error by C-N scheme at ( T = 10 , M a x E r r o r = 6.5276 e 07 ) , ( T = 20 , M a x E r r o r = 1.7244 e 08 ) , α = 1.5 left to right, respectively, for Example 1.
Figure 1. The Maximum error by C-N scheme at ( T = 10 , M a x E r r o r = 6.5276 e 07 ) , ( T = 20 , M a x E r r o r = 1.7244 e 08 ) , α = 1.5 left to right, respectively, for Example 1.
Symmetry 12 00485 g001
Figure 2. The exact (left) and numerical (right) solution by C-N scheme at T = 1 , α = 1.5 , τ = 0.01 = h for Example 1.
Figure 2. The exact (left) and numerical (right) solution by C-N scheme at T = 1 , α = 1.5 , τ = 0.01 = h for Example 1.
Symmetry 12 00485 g002
Figure 3. Numerical and exact solution by C-N scheme at α = 1.5 , τ = h = 0.01 , with ( T = 10 , T = 30 , T = 40 ) left to right-down respectively, for Example 1.
Figure 3. Numerical and exact solution by C-N scheme at α = 1.5 , τ = h = 0.01 , with ( T = 10 , T = 30 , T = 40 ) left to right-down respectively, for Example 1.
Symmetry 12 00485 g003
Figure 4. The exact (left) and numerical (right) solution by C-N scheme for the FCDE at ( h = τ = 0.005 , α = 1.5 , ( t = 5 , max e r r o r = 4.0657 e 05 ) for Example 2.
Figure 4. The exact (left) and numerical (right) solution by C-N scheme for the FCDE at ( h = τ = 0.005 , α = 1.5 , ( t = 5 , max e r r o r = 4.0657 e 05 ) for Example 2.
Symmetry 12 00485 g004
Figure 5. The exact (left) and numerical (right) solution by C-N scheme for the FCDE at ( h = τ = 0.01 , ( t = 2 , max E r r o r = 4.2158 e 04 ) , α = 1.75 ) for Example 2.
Figure 5. The exact (left) and numerical (right) solution by C-N scheme for the FCDE at ( h = τ = 0.01 , ( t = 2 , max E r r o r = 4.2158 e 04 ) , α = 1.75 ) for Example 2.
Symmetry 12 00485 g005
Figure 6. The exact (left) and numerical (right) solution by C-N scheme at ( h = τ = 0.0025 , ( t = 0.1 , max E r r o r = 1.4 e 03 , α = 1.1 ) for Example 4.
Figure 6. The exact (left) and numerical (right) solution by C-N scheme at ( h = τ = 0.0025 , ( t = 0.1 , max E r r o r = 1.4 e 03 , α = 1.1 ) for Example 4.
Symmetry 12 00485 g006
Table 1. The maximum error and convergence order of the C-N scheme for FDE in example 1.
Table 1. The maximum error and convergence order of the C-N scheme for FDE in example 1.
α = 1.25 α = 1.5 α = 1.8
Δ t Δ x Max-ErrorOrderMax-ErrorOrderMax-ErrorOrder
1/501/504.9807e−044.0046e−041.4048e−04
1/1001/1001.0660e−042.22418.8946e−052.17073.6848e−051.9307
1/2001/2002.4413e−052.12652.0643e−052.10739.4393e−061.9648
1/4001/4005.8239e−062.06764.9592e−062.05752.3887e−061.9825
1/8001/8001.4211e−062.03501.2146e−062.02966.0078e−071.9913
Table 2. The maximum error and convergence order for FCDE in example 2.
Table 2. The maximum error and convergence order for FCDE in example 2.
T = 1 T = 5
Δ t Δ x Max-ErrorOrderMax-ErrorOrder
1/501/501.4048e−042.5297e−05
1/1001/1003.6848e−051.93077.4748e−061.7589
1/2001/2009.4393e−061.96482.0122e−061.8933
1/4001/4002.3887e−061.98254.9017e−072.0374
1/8001/8006.0078e−071.99131.0620e−072.2065
Table 3. The maximum error and convergence order by C-N for SFCDE in example 2 at T = 1 , α = 1.55 .
Table 3. The maximum error and convergence order by C-N for SFCDE in example 2 at T = 1 , α = 1.55 .
Δ t Δ x Max-ErrorOrder
1/501/502.6e−03
1/1001/1007.695e−041.7563
1/1501/1502.144e−041.8436
1/2001/2005.688e−051.9143
Table 4. The Maximum error and convergence order produced by C-N scheme for example 3 at T = 1 , N t = 100 .
Table 4. The Maximum error and convergence order produced by C-N scheme for example 3 at T = 1 , N t = 100 .
α = 1.35 α = 1.5 α = 1.75
Δ x Max-ErrorOrderMax-ErrorOrderMax-ErrorOrder
1/504.5e−032.8e−031.7e−03
1/1002.7e−030.73701.6e−030.80748.9641–040.97224
1/2001.6e−030.75498.6405e−040.88894.6491e−040.8981
1/4009.5896e−040.73854.7955e−040.84942.4086e−040.9488
1/8005.7034e−040.74962.6609e−040.84981.2473e−040.9494
Table 5. The Maximum error and error-rate produced by C-N scheme for example 4 at t = 0.1 .
Table 5. The Maximum error and error-rate produced by C-N scheme for example 4 at t = 0.1 .
α = 1.25 α = 1.55
Δ t Δ x Max-ErrorError-RateMax-ErrorError-Rate
1/501/501.91e−021.52e−02
1/1001/1009.9e−031.937.9e−031.9
1/2001/2005.2e−031.904.3e−031.84
1/4001/4002.8e−031.862.4e−031.79
1/8001/8001.6e−031.751.4e−031.7

Share and Cite

MDPI and ACS Style

Anley, E.F.; Zheng, Z. Finite Difference Approximation Method for a Space Fractional Convection–Diffusion Equation with Variable Coefficients. Symmetry 2020, 12, 485. https://doi.org/10.3390/sym12030485

AMA Style

Anley EF, Zheng Z. Finite Difference Approximation Method for a Space Fractional Convection–Diffusion Equation with Variable Coefficients. Symmetry. 2020; 12(3):485. https://doi.org/10.3390/sym12030485

Chicago/Turabian Style

Anley, Eyaya Fekadie, and Zhoushun Zheng. 2020. "Finite Difference Approximation Method for a Space Fractional Convection–Diffusion Equation with Variable Coefficients" Symmetry 12, no. 3: 485. https://doi.org/10.3390/sym12030485

APA Style

Anley, E. F., & Zheng, Z. (2020). Finite Difference Approximation Method for a Space Fractional Convection–Diffusion Equation with Variable Coefficients. Symmetry, 12(3), 485. https://doi.org/10.3390/sym12030485

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