Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
MHD Bioconvection Flow and Heat Transfer of Nanofluid through an Exponentially Stretchable Sheet
Previous Article in Journal
Asymmetry of Motor Behavior and Sensory Perception: Which Comes First?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fourth-Order Difference Approximation for Time-Fractional Modified Sub-Diffusion Equation

1
Department of Mathematics, Al-Fajar University, Mari Indus 42350, Mianwali, Pakistan
2
Department of Applied Mathematics and Statistics, Institute of Space Technology, Islamabad 44000, Pakistan
3
BIC-ESAT, College of Engineering, Peking University, Beijing 100871, China
4
School of Mathematical Sciences, Universiti Sains Malaysia, Pulau Pinang 11800, Malaysia
5
Faculty of Mathematics and Statistics, Ton Duc Thang University, Ho Chi Minh City 72915, Vietnam
6
Department of Mathematics, College of Arts and Sciences, Prince Sattam bin Abdulaziz University, Wadi Al-Dawaser 11991, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(5), 691; https://doi.org/10.3390/sym12050691
Submission received: 21 February 2020 / Revised: 16 March 2020 / Accepted: 19 March 2020 / Published: 1 May 2020

Abstract

:
Fractional differential equations describe nature adequately because of the symmetry properties which describe physical and biological processes. In this article, a fourth-order new implicit difference scheme is formulated and applied to solve the two-dimensional time-fractional modified sub-diffusion equation involving two times Riemann–Liouville fractional derivatives. The stability of the fourth-order implicit difference scheme is investigated using the von Neumann technique. The proposed scheme is shown to be unconditionally stable. Numerical examples are given to illustrate the feasibility of the proposed scheme.

1. Introduction

Consider the two-dimensional time-fractional modified sub-diffusion Equation (2D-TFMSDE) [1]:
u ( x , y , t ) t = A 1 α t 1 α + B 1 β t 1 β 2 u ( x , y , t ) x 2 + 2 u ( x , y , t ) y 2 + f ( x , y , t ) ,
subject to the conditions
u ( x , y , 0 ) = φ ( x , y ) ,
and
u ( 0 , y , t ) = φ 1 ( y , t ) , u ( L , y , t ) = φ 2 ( y , t ) , u ( x , 0 , t ) = φ 3 ( x , t ) , u ( x , L , t ) = φ 4 ( x , t ) , 0 x , y L , 0 t T .
Here, φ , φ 1 , φ 2 , φ 3 , and φ 4 are defined functions; A and B are constants; and 1 α t 1 α and 1 β t 1 β are the fractional order Riemann–Liouville derivatives of order 1 α and 1 β , respectively. The 2D-TFMSDE describes processes that become less anomalous with time. This decrease in anomalous behavior is due to the inclusion of a secondary fractional time derivative which acts on a diffusion operator [2]. This equation arises, for example, in econophysics where there is increasing interest in modeling using continuous time random walks [3]. It should be noted that the crossover between more and less anomalous behavior has been noted in the volatility of some share prices [4].
In the last two decades, there has been an increasing interest in fractional calculus because of applications in various fields such as physics, chemistry, biological science, viscoelastic, and fluid mechanics phenomena [5,6,7,8]. Fractional differential equations have attracted a number of numerical researchers to developed numerical methods. Many researchers have studied the fractional differential equation described in Equation (1) by various numerical techniques. Liu et al. [2] developed an implicit difference technique for the solution of modified sub-diffusion equation within a bounded domain. The stability and convergence was investigated by energy analysis. Li and Wang [9] developed an improved efficient difference technique for a time-fractional modified sub-diffusion equation. They replaced the Riemann–Liouville fractional order derivative, second-order space derivatives and nonlinear inhomogenous part with weighted and shifted Grünwald–Letnikov formula, compact difference approximation, and interpolation formula, respectively. Dehghan et al. [10] applied the difference method for Riemann–Liouville derivative of fractional order. They converted the semi-discrete scheme into full discretized form by integration and then used Legendre spectral element technique for the equation. Cao et al. [11] introduced midpoint implicit technique to solve a modified sub-diffusion equation of fractional order. They used weighted and shifted Grünwald–Letnikov fractional formula and the compact difference formula for Riemann–Liouville fractional derivative and second-order space derivative, respectively. The stability analysis and convergence of the technique was analyzed. Ding and Li [12] developed Riemann–Liouville second-order derivative of fractional order and constructed new types of numerical difference approximations. The stability and convergence were studied by von Neumann analysis. Other researchers have also constructed compact and high order numerical schemes by various techniques for time-fractional modified sub-diffusion equation and analyzed the stability andconvergence [13,14,15,16,17,18]. The advantage of the compact methods is that better accuracy can be obtained with reduced space and time step size as opposed to the standard approach. However compact methods require increased computational costs. Other related studies on fractional calculus and its applications can be found in [19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36].
In this paper, the fourth-order implicit difference method is formulated to solve the 2D-TFMSDE. The construction of the fourth-order difference method improves the rate of convergence to fourth order in space from the usual second order in space [37]. The fourth-order numerical method requires a more complex procedure than the second-order numerical method and generates a smaller size of linear systems. The new discretized form of Riemann–Liouville integral operator is used. Previous studies used Grünwald–Letnikov, weighted and shifted Grünwald–Letnikov formula, compact Riemann–Liouville fractional formula, and discretized Riemann–Liouville integral for integro-differential equation.
There are two novelties of this article: first, the fourth-order difference approximation and the scheme possesses temporal convergence of first order and spatial convergence of fourth order, which reduces the computational cost greatly, and thus improves the computational efficiency; and, second, the new formulation of the Riemann–Liouville fractional derivative based on Jumarie properties. The stability of the proposed fourth-order implicit difference scheme is investigated by von Neumann analysis and a numerical experiment is performed.
The Riemann–Liouville fractional derivative can be defined as [33]:
0 D t 1 α u ( x , y , t n ) = 1 Γ ( α ) t 0 t n u ( x , y , ξ ) ( t n ξ ) 1 α d ξ ,
and
0 D t 1 α u ( x , y , t n ) = t 1 Γ ( α ) 0 t n u ( x , y , ξ ) ( t n ξ ) 1 α d ξ , = t 1 Γ ( α ) 0 t n ( t n ξ ) α 1 u ( x , y , ξ ) d ξ .
Based on Jumarie properties [30], Equation (5) can be written as follows:
= t 1 α Γ ( α ) 0 t n u ( x , y , ξ ) ( d ξ ) α , = t 1 Γ ( α + 1 ) k = 0 n 1 t k t k + 1 u ( x , y , ξ ) ( d ξ ) α ,   = t 1 Γ ( α + 1 ) k = 0 n 1 u ( x , y , t n k ) t k t k + 1 ξ 0 ( d ξ ) α ,
with the property 0 x ξ a ( d ξ ) b = Γ ( a + 1 ) Γ ( b + 1 ) Γ ( a + b + 1 ) x a + b ,
= t τ α Γ ( α + 1 ) k = 0 n 1 u ( x , y , t n k ) ( ( k + 1 ) α ( k ) α ) , 0 D t 1 α u ( x , y , t n ) = τ α 1 Γ ( α + 1 ) k = 0 n 1 b k ( α ) ( u ( x , y , t n k ) u ( x , y , t n k 1 ) ) ,
and b k ( α ) = ( k + 1 ) α ( k ) α , k = 0 , 1 , 2 , , n 1 .
Lemma 1.
The α ( 0 < α < 1 ) order Riemann–Liouville fractional derivative of the function u ( x , y , t ) on [ 0 , T ] can be defined in discretized form as,
0 D t 1 α u ( x , y , t n ) = τ α 1 Γ ( α + 1 ) k = 0 n 1 b j ( α ) ( u ( x , y , t n k ) u ( x , y , t n k 1 ) ) .
Lemma 2.
The coefficients b n ( α ) ( n = 0 , 1 , 2 , 3 ) satisfy the properties as follows [33]:
(i) 
b 0 ( α ) = 1 , b n ( α ) > 0 , n = 0 , 1 , 2 , . ;
(ii) 
b n 1 ( α ) > b n ( α ) , n = 1 , 2 , . ;
(iii) 
There exists a positive constant C > 0 . such that τ C b n ( α ) τ α , n = 1 , 2 , . ; and
(iv) 
k = 0 n b k ( α ) τ α = ( n + 1 ) α T α .
Lemma 1 and 2 are the same for β as well (like α Lemma 1 and 2 are also same for β ).

2. Fourth-Order Implicit Difference Scheme

A fourth-order implicit difference scheme is constructed to solve the modified fractional sub-diffusion equation (Equation (1)). The discretized form of fractional derivatives and fourth-order difference approximation is used for second-order derivatives. The grid point ( x i , y j , t n ) with x i = i Δ x , in the x-direction with i = 1 , 2 , , M x 1 , Δ x = L M x . In the y-direction y j = j Δ y , with j = 1 , 2 , , M y 1 , Δ y = L M y . The time step is t n = n τ , n = 1 , 2 , , N where τ = T N .
To obtain the compact difference scheme for the space derivatives, u i + 1 , j n and u i 1 , j n are expanded about u ( x i , y j , t n ) by Taylor series. After adding and simplifying the calculation, as in [14], we have
2 u x 2 | i n = δ x 2 / h 2 I + 1 12 δ x 2 u i , j n + O ( h x 4 ) ,
similarly
2 u y 2 | j n = δ y 2 / h 2 I + 1 12 δ y 2 u i , j n + O ( h y 4 ) ,
For the discretization of Equation (1), using Lemma 1 and fourth-order difference approximation in Equations (7) and (8), we have
I + 1 12 δ y 2 I + 1 12 δ x 2 u i , j n u i , j n 1 = S 1 k = 0 n 1 b k ( α ) I + δ y 2 12 δ x 2 u i , j n k u i , j n k 1 + S 2 k = 0 n 1 b k ( α ) I + δ x 2 12 δ y 2 u i , j n k u i , j n k 1 + S 3 k = 0 n 1 b k ( β ) I + δ y 2 12 δ x 2 u i , j n k u i , j n k 1 + S 4 k = 0 n 1 b k ( β ) I + δ x 2 12 δ y 2 u i , j n k u i , j n k 1 + τ I + 1 12 δ y 2 I + 1 12 δ x 2 f i , j n .
Here,
S 1 = A τ α Γ ( α + 1 ) Δ x 2 , S 2 = A τ α Γ ( α + 1 ) Δ y 2 , S 3 = B τ β Γ ( β + 1 ) Δ x 2 , S 4 = B τ β Γ ( β + 1 ) Δ y 2 ,
and
δ x 2 u i , j n = u i + 1 , j n 2 u i , j n + u i 1 , j n , δ y 2 u i , j n = u i , j + 1 n 2 u i , j n + u i , j 1 n .
After simplification of Equation (9), the implicit difference scheme for Equations (1)–(3) with the associated conditions is obtained, as follows:
1 + 1 12 δ x 2 + 1 12 δ y 2 + 1 144 δ x 2 δ y 2 u i , j n S 1 I + 1 12 δ y 2 δ x 2 u i , j n S 2 I + 1 12 δ x 2 δ y 2 u i , j n S 3 I + 1 12 δ y 2 δ x 2 u i , j n S 4 I + 1 12 δ x 2 δ y 2 u i , j n = 1 + 1 12 δ x 2 + 1 12 δ y 2 + 1 144 δ x 2 δ y 2 u i , j n 1 S 1 b n 1 α I + 1 12 δ y 2 δ x 2 u i , j 0 S 2 b n 1 α I + 1 12 δ x 2 δ y 2 u i , j 0 S 3 b n 1 β I + 1 12 δ y 2 δ x 2 u i , j 0 S 4 b n 1 β I + 1 12 δ x 2 δ y 2 u i , j 0 k = 1 n 1 ( b k 1 ( α ) b k ( α ) ) S 1 I + 1 12 δ y 2 δ x 2 u i , j n k + S 2 I + 1 12 δ x 2 δ y 2 u i , j n k k = 1 n 1 ( b k 1 ( β ) b k ( β ) ) S 3 I + 1 12 δ y 2 δ x 2 u i , j n k + S 4 I + 1 12 δ x 2 δ y 2 u i , j n k + τ ( 1 + 1 12 δ x 2 + 1 12 δ y 2 + 1 144 δ x 2 δ y 2 ) f i , j n ,
where i = 1 , 2 , , M x 1 , j = 1 , 2 , , M y 1 and n = 1 , 2 , , N .
with
u i , j 0 = φ ( x i , y j ) ,
u 0 , j n = φ 1 ( y j , t n ) , u M x , j n = φ 2 ( y j , t n ) , u i , 0 n = φ 3 ( x i , t n ) , u i , M y n = φ 4 ( x i , t n ) , i = 1 , 2 , M x 1 , j = 1 , 2 , M y 1 a n d n = 1 , 2 , N .

Stability

We next investigate the stability of the proposed scheme by von Neumann analysis method following the approach in [1,18]. Let U i , j n be the exact solution for Equation (12). We have
1 + 1 12 δ x 2 + 1 12 δ y 2 + 1 144 δ x 2 δ y 2 U i , j n S 1 I + 1 12 δ y 2 δ x 2 U i , j n S 2 I + 1 12 δ x 2 δ y 2 U i , j n S 3 I + 1 12 δ y 2 δ x 2 U i , j n S 4 I + 1 12 δ x 2 δ y 2 U i , j n = 1 + 1 12 δ x 2 + 1 12 δ y 2 + 1 144 δ x 2 δ y 2 U i , j n 1 S 1 b n 1 α I + 1 12 δ y 2 δ x 2 U i , j 0 S 2 b n 1 α I + 1 12 δ x 2 δ y 2 U i , j 0 S 3 b n 1 β I + 1 12 δ y 2 δ x 2 U i , j 0 S 4 b n 1 β I + 1 12 δ x 2 δ y 2 U i , j 0 k = 1 n 1 ( b k 1 ( α ) b k ( α ) ) S 1 I + 1 12 δ y 2 δ x 2 U i , j n k + S 2 I + 1 12 δ x 2 δ y 2 U i , j n k k = 1 n 1 ( b k 1 ( β ) b k ( β ) ) S 3 I + 1 12 δ y 2 δ x 2 U i , j n k + S 4 I + 1 12 δ x 2 δ y 2 U i , j n k + τ ( 1 + 1 12 δ x 2 + 1 12 δ y 2 + 1 144 δ x 2 δ y 2 ) f i , j n ,
where i = 1 , 2 , , M x 1 , j = 1 , 2 , , M y 1 and n = 1 , 2 , , N 1 .
Next, the error is defined as
e i , j n = u i , j n U i , j n ,
where e i , j n satisfies Equation (15) and
1 + 1 12 δ x 2 + 1 12 δ y 2 + 1 144 δ x 2 δ y 2 e i , j n S 1 I + 1 12 δ y 2 δ x 2 e i , j n S 2 I + 1 12 δ x 2 δ y 2 e i , j n S 3 I + 1 12 δ y 2 δ x 2 e i , j n S 4 I + 1 12 δ x 2 δ y 2 e i , j n = 1 + 1 12 δ x 2 + 1 12 δ y 2 + 1 144 δ x 2 δ y 2 e i , j n 1 S 1 b n 1 α I + 1 12 δ y 2 δ x 2 e i , j 0 S 2 b n 1 α I + 1 12 δ x 2 δ y 2 e i , j 0 S 3 b n 1 β I + 1 12 δ y 2 δ x 2 e i , j 0 S 4 b n 1 β I + 1 12 δ x 2 δ y 2 e i , j 0 k = 1 n 1 ( b k 1 ( α ) b k ( α ) ) S 1 I + 1 12 δ y 2 δ x 2 e i , j n k + S 2 I + 1 12 δ x 2 δ y 2 e i , j n k k = 1 n 1 ( b k 1 ( β ) b k ( β ) ) S 3 I + 1 12 δ y 2 δ x 2 e i , j n k + S 4 I + 1 12 δ x 2 δ y 2 e i , j n k .
The error conditions can be defined as
e 0 , j n = e i , 0 n = e i , j 0 = 0 , e i , M y n = e M x , j n = 0 .
The grid functions for n = 1 , 2 , , N are
e n ( x , y ) = e i , j n , w h e n x i Δ x 2 < x x i + Δ x 2 , y j Δ y 2 < y y j + Δ y 2 , 0 , w h e n 0 x Δ x 2 o r L Δ x 2 x L , 0 , w h e n 0 y Δ y 2 o r L Δ y 2 y L .
Now, e n ( x , y ) is defined in a Fourier series as
e n ( x , y ) = l 1 , l 2 = λ n ( l 1 , l 2 ) e 2 1 π ( l 1 x / L + l 2 y / L ) ,
where
λ n ( l 1 , l 2 ) = 1 L 0 L 0 L e n ( x , y ) e 2 1 π ( l 1 x / L + l 2 y / L ) d x d y .
From the definition of l 2 norm and Parseval equality, we have
e n 2 = i = 1 M x 1 j = 1 M y 1 Δ x Δ y | e i , j n | 2 = l 1 , l 2 = | λ n ( l 1 , l 2 ) | 2 .
Supposing that
e i , j n = λ n e 1 ( σ 1 i Δ x + σ 2 j Δ y ) ,
where σ 1 = 2 π l 1 / L , σ 2 = 2 π l 2 / L and substituting Equation (23) into Equation (17), we obtain
λ n = 1 ( 144 ( 1 + μ 2 + μ 3 ) 12 ( μ 1 + μ 4 + μ 5 ) + μ 0 ) ( λ n 1 ( 144 12 μ 1 + μ 0 ) + 144 ( μ 2 b n 1 ( α ) + μ 3 b n 1 ( β ) ) 12 ( μ 4 b n 1 ( α ) + μ 5 b n 1 ( β ) ) λ 0 + ( 144 μ 2 12 μ 4 ) k = 1 n 1 ( b k 1 ( α ) b k ( α ) ) λ n k + ( 144 μ 3 12 μ 5 ) k = 1 n 1 ( b k 1 ( β ) b k ( β ) ) λ n k ) .
where
μ 0 = 16 s i n 2 ( σ 1 Δ x 2 ) s i n 2 ( σ 2 Δ y 2 ) 0 ,
μ 1 = 4 s i n 2 ( σ 1 Δ x 2 ) + s i n 2 ( σ 2 Δ y 2 ) 0 ,
μ 2 = 4 S 1 s i n 2 ( σ 1 Δ x 2 ) + S 2 s i n 2 ( σ 2 Δ y 2 ) 0 ,
μ 3 = 4 S 3 s i n 2 ( σ 1 Δ x 2 ) + S 4 s i n 2 ( σ 2 Δ y 2 ) 0 .
μ 4 = 16 ( S 1 + S 2 ) s i n 2 ( σ 1 Δ x 2 ) s i n 2 ( σ 2 Δ y 2 ) 0 ,
μ 5 = 16 ( S 3 + S 4 ) s i n 2 ( σ 1 Δ x 2 ) s i n 2 ( σ 2 Δ y 2 ) 0 .
Proposition 1.
If λ n ( n = 1 , 2 , , N ) satisfy Equation (24), then | λ n | | λ 0 | .
Proof. 
Mathematical induction is used. Let us take n = 1 in Equation (24)
λ 1 = ( 144 ( 1 + μ 2 + μ 3 ) 12 ( μ 1 + μ 4 + μ 5 ) + μ 0 ) λ 0 ( 144 ( 1 + μ 2 + μ 3 ) 12 ( μ 1 + μ 4 + μ 5 ) + μ 0 ) ,
and as μ 0 , μ 1 , μ 2 , μ 3 , μ 4 , μ 5 0 and, in Lemma 2, b 0 ( α ) = b 0 ( β ) = 1 , then
| λ 1 | | λ 0 | .
Now, assuming that
| λ s | | λ 0 | ; s = 1 , 2 , , n 1 ,
and as 0 < α , β < 1 , from Equation (24) and Lemma 2, we obtain
| λ n | 1 ( 144 ( 1 + μ 2 + μ 3 ) 12 ( μ 1 + μ 4 + μ 5 ) + μ 0 ) [ | λ n 1 | ( 144 12 μ 1 + μ 0 ) + 144 ( μ 2 b n 1 ( α ) + μ 3 b n 1 ( β ) ) 12 ( μ 4 b n 1 ( α ) + μ 5 b n 1 ( β ) ) | λ 0 | + ( 144 μ 2 12 μ 4 ) k = 1 n 1 ( b k 1 ( α ) b k ( α ) ) | λ n k | + ( 144 μ 3 12 μ 5 ) k = 1 n 1 ( b k 1 ( β ) b k ( β ) ) | λ n k | ] , 1 ( 144 ( 1 + μ 2 + μ 3 ) 12 ( μ 1 + μ 4 + μ 5 ) + μ 0 ) [ 144 12 μ 1 + μ 0 12 ( μ 4 b n 1 ( α ) + μ 5 b n 1 ( β ) ) + 144 ( μ 2 b n 1 ( α ) + μ 3 b n 1 ( β ) ) + ( 144 μ 2 12 μ 4 ) k = 1 n 1 ( b k 1 ( α ) b k ( α ) ) + ( 144 μ 3 12 μ 5 ) k = 1 n 1 ( b k 1 ( β ) b k ( β ) ) ] | λ 0 | , = 1 ( 144 ( 1 + μ 2 + μ 3 ) 12 ( μ 1 + μ 4 + μ 5 ) + μ 0 ) [ 144 12 μ 1 + μ 0 12 ( μ 4 b n 1 ( α ) + μ 5 b n 1 ( β ) ) + 144 ( μ 2 b n 1 ( α ) + μ 3 b n 1 ( β ) ) + ( 144 μ 2 12 μ 4 ) ( 1 b n 1 ( α ) ) + ( 144 μ 3 12 μ 5 ) ( 1 b n 1 ( β ) ) ] | λ 0 | , = 144 ( 1 + μ 2 + μ 3 ) 12 ( μ 1 + μ 4 + μ 5 ) + μ 0 144 ( 1 + μ 2 + μ 3 ) 12 ( μ 1 + μ 4 + μ 5 ) + μ 0 | λ 0 | , | λ n | | λ 0 | .
Hence, Equation (12) satisfies
λ n 2 λ 0 2 .
This shows that the proposed scheme in Equation (12) is unconditionally stable. □

3. Numerical Experiments

We consider a 2D TFMSDE to test the numerical scheme. The E errors between the numerical and exact solutions are compared with the previous literature. The E error is defined as
E = max 0 i M x 1 , 0 j M y 1 , 0 k N | u ( x i , y j , t k ) u i , j k | .
the formula for computational orders of the method presented in time variables is [18]:
C 1 -order= l o g 2 E ( 2 τ , h ) E ( τ , h ) ,
and in space variables
C 2 -order= l o g 2 E ( 16 τ , 2 h ) E ( τ , h ) .
Example 1.
The following 2D-TFMSDE is considered [1].
u ( x , y , t ) t = 1 α t 1 α + 1 β t 1 β 2 u ( x , y , t ) x 2 + 2 u ( x , y , t ) y 2 + f ( x , y , t ) , 0 t 1 ,
where
f ( x , y , t ) = s i n ( x + y ) ( 1 + α + β ) t α + β + 2 Γ ( 2 + α + β ) Γ ( 1 + 2 α + β ) t 2 α + β + 2 Γ ( 2 + α + β ) Γ ( 1 + α + 2 β ) t α + 2 β ,
subject to the conditions
u ( x , y , 0 ) = 0 , 0 x , y 1 ,
u ( 0 , y , t ) = t 1 + α + β s i n ( y ) , u ( 1 , y , t ) = t 1 + α + β s i n ( 1 + y ) , u ( x , 0 , t ) = t 1 + α + β s i n ( x ) , u ( x , 1 , t ) = t 1 + α + β s i n ( 1 + x ) , 0 x , y 1 , 0 t T .
The exact solution is
u ( x , y , t ) = t 1 + α + β s i n ( x + y ) .
The proposed fourth-order new implicit difference scheme is applied to a 2D-TFMSDE. Table 1 compares the numerical results of the proposed scheme with the previous studies. In [17], the problem was solved using a compact difference scheme and, in [1], by a finite difference scheme. The scheme proposed in this paper shows better accuracy as the space and time steps are reduced. The order of the proposed method is O ( ( Δ x ) 4 + ( Δ y ) 4 ) in space variables and O ( τ ) in time variable [37]. The computational C 1 -order and C 2 -order are close to theoretical order, as can be seen in Table 1 and Table 2, respectively. For more confirmation, the proposed method is also extended to 2D RSP-HGSGF and the numerical results compared with the exact solution. The obtained results in Table 3 show very good accuracy for different values of γ , space, and time steps for Example 2. The 3D plots in Figure 1 and Figure 2 show the comparison of approximate and exact solution, and the 2D plot in Figure 3 shows the difference between the approximate and exact solution at α = 0.25 , β = 0.45 , T = 1 , M x = M y = 10 , and N = 160 for Example 1 with good agreement and high accuracy.
For a wide range of τ , the results remain stable, indicating that the developed scheme is unconditionally stable (as shown in the theoretical analysis).
Example 2.
The following 2D Rayleigh–Stokes problem for heated generalized second-grade fluid (RSP-HGSGF) [22] is considered
u ( x , y , t ) t = 1 γ t 1 γ 2 u ( x , y , t ) x 2 + 2 u ( x , y , t ) y 2 + 2 u ( x , y , t ) x 2 + 2 u ( x , y , t ) y 2 + f ( x , y , t ) ,
where
f ( x , y , t ) = ( 1 + γ ) t 1 + γ 2 Γ ( 2 + γ ) Γ ( 1 + 2 γ ) t 2 γ 2 t 1 + 2 γ e x + y ,
subject to the initial condition
u ( x , y , 0 ) = 0 , 0 x , y 1 ,
u ( 0 , y , t ) = e y t 1 + γ , u ( 1 , y , t ) = e 1 + y t 1 + γ , u ( x , 0 , t ) = e x t 1 + γ , u ( x , 1 , t ) = e 1 + x t 1 + γ ,   0 x , y 1 , 0 t T .
The exact solution is
u ( x , y , t ) = e x + y t 1 + γ .

4. Conclusions

In this paper, a fourth-order implicit difference scheme for time-fractional modified sub-diffusion equation is developed, analyzed, and applied. The von Neumann analysis is used to show the scheme is unconditionally stable. The results of an application to a particular examples are discussed. The numerical results seem to confirm the theoretical results and indicate the effectiveness and feasibility of the new formulated scheme. The numerical method in this paper can be applied to other types of fractional differential equations and higher-dimensional problems.

Author Contributions

Conceptualization, U.A., M.U., I.K. and F.A.A.; methodology, U.A., M.S., K.S.N., and M.U.; software, M.U., and U.A.; validation, I.K., and F.A.A.; formal analysis, U.A., I.K., M.U., and K.S.N.; investigation, U.A., and F.A.A.; writing–original draft preparation, U.A., M.S., K.S.N., and M.S.; writing–review and editing, U.A., M.S., M.U., F.A.A., I.K., and K.S.N.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors are thankful to the reviewers for the useful suggestions to improve the quality of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ali, U.; Abdullah, F.A.; Mohyud-Din, S.T. Modified implicit fractional difference scheme for 2D modified anomalous fractional sub-diffusion equation. Adv. Diff. Equ. 2017, 1–14. [Google Scholar] [CrossRef]
  2. Liu, F.; Yang, C.; Burrage, K. Numerical method and analytical technique of the modified anomalous subdiffusion equation with a nonlinear source term. J. Comput. Appl. Math. 2009, 231, 160–176. [Google Scholar] [CrossRef] [Green Version]
  3. Langlands, T.A.M. Solution of a modified fractional diffusion equation. Physica A 2006, 367, 136–144. [Google Scholar] [CrossRef] [Green Version]
  4. Jiang, W.; Chen, Z. A collocation method based on reproducing kernel for a modified anomalous subdiffusion equation. Numer. Meth. Part. Differ. Equ. 2014, 30, 289–300. [Google Scholar] [CrossRef]
  5. Miller, K.S.; Ross, B. An Introductional the Fractional Calculus and Fractional Differential Equations; Academic Press: New York, NY, USA; London, UK, 1974. [Google Scholar]
  6. Oldham, K.B.; Spanier, J. The Fractional Calculus; Academic Press: New York, NY, USA; London, UK, 1974. [Google Scholar]
  7. Liangliang, M.A.; Dongbing, L. An implicit difference approximation for fractional cable equation in high-dimensional case. J. Liao. Tech. Univ. 2014, 4, 24. [Google Scholar]
  8. Zhuang, P.; Liu, F.; Anh, V.; Turne, I. New solution and analytical techniques of the implicit numerical methods for the anomalous sub-diffusion equation. SIAM J. Numer. Anal. 2008, 46, 1079–1095. [Google Scholar] [CrossRef] [Green Version]
  9. Li, Y.; Wang, D. Improved efficient difference method for the modified anomalous sub-diffusion equation with a nonlinear source term. Int. J. Comput. Math. 2016, 94, 821–840. [Google Scholar] [CrossRef]
  10. Dehghan, M.; Abbaszadeh, M.; Mohebbi, A. Legendre spectral element method for solving time fractional modified anomalous sub-diffusion equation. Appl. Math. Model 2015, 40, 3635–3654. [Google Scholar] [CrossRef]
  11. Cao, X.; Xianxian, C.; Wen, L. The implicit midpoint method for the modified anomalous sub-diffusion equation with a nonlinear source term. J. Comput. Appl. Math. 2016, 1–12. [Google Scholar] [CrossRef]
  12. Ding, H.; Li, C. High-order compact difference schemes for the modified anomalous subdiffusion equation. Numer. Meth. Part. Differ. Equ. 2015, 213–242. [Google Scholar] [CrossRef] [Green Version]
  13. Wang, Z.; Vong, S. Compact difference schemes for the modified anomalous fractional sub-diffusion equation and the fractional diffusion-wave equation. J. Comput. Phys. 2014, 277, 1–15. [Google Scholar] [CrossRef] [Green Version]
  14. Al-Shibani, F.S.; Ismail, A.I.M. Compact Crank-Nicolson and Du Fort-Frankel method for the solution of the anomalous time fractional diffusion equation. Int. J. Comput. Methods 2013, 12, 1550041. [Google Scholar] [CrossRef]
  15. Adell, J.A.; Lekuona, A. Note on two extensions of the classical formula for sums of powers on arithmetic progressions. Adv. Differ. Equ. 2017, 1–5. [Google Scholar] [CrossRef] [Green Version]
  16. Ali, U.; Abdullah, F.A.; Ismail, A.I. Crank-Nicolson finite difference method for two-dimensional fractional sub-diffusion equation. J. Interpolat. Approx. Sci. Comput. 2017, 2, 18–29. [Google Scholar] [CrossRef]
  17. Abbaszadeh, M.; Mohebbi, A. A fourth-order compact solution of the two-dimensional modified anomalous fractional sub-diffusion equation with a nonlinear source term. Comput. Math. Appl. 2013, 66, 1345–1359. [Google Scholar] [CrossRef]
  18. Mohebbi, A.; Abbaszadeh, M.; Dehghan, M. A high-order and unconditionally stable scheme for the modified anomalous fractional sub-diffusion equation with a nonlinear source term. J. Comput. Phys. 2013, 240, 36–48. [Google Scholar] [CrossRef]
  19. Gao, W.; Yel, G.; Baskonus, H.M.; Cattani, C. Complex solitons in the conformable (2+ 1)-dimensional Ablowitz-Kaup-Newell-Segur equation. AIMS Math. 2019, 5, 507–521. [Google Scholar] [CrossRef]
  20. Gao, W.; Veeresha, P.; Prakasha, D.G.; Baskonus, H.M.; Yel, G. A powerful approach for fractional Drinfeld–Sokolov–Wilson equation with Mittag-Leffler law. Alex. Eng. J. 2019, 58, 1301–1311. [Google Scholar] [CrossRef]
  21. Gao, W.; Ghanbari, B.; Baskonus, H.M. New numerical simulations for some real world problems with Atangana–Baleanu fractional derivative. Chaos Solitons Fract. 2019, 128, 34–43. [Google Scholar] [CrossRef]
  22. Chen, C.M.; Liu, F.; Anh, V. Numerical analysis of the Rayleigh–Stokes problem for a heated generalized second grade fluid with fractional derivatives. Appl. Math. Comput. 2008, 204, 340–351. [Google Scholar] [CrossRef] [Green Version]
  23. Ali, U.; Abdullah, F.A. Modified implicit difference method for one-dimensional fractional wave equation. In Proceedings of the AIP Conference Proceedings, Penang, Malaysia, 10–12 December 2018; Volume 2184, p. 060021. [Google Scholar]
  24. Sokolov, I.M.; Klafter, J. From diffusion to anomalous diffusion: A century after Einstein’s Brownian motion. Chaos 2005, 15, 026103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Zhuang, P.; Liu, F.; Turner, L.; Anh, V. Galerkin finite element method and error analysis for the fractional cable equation. Numer. Algorithms 2016, 72, 447–466. [Google Scholar] [CrossRef] [Green Version]
  26. Oldham, K.B.; Spanier, J. The Fractional Calculus: Theory and Application of Differentiation and Integration to Arbitrary Order; Academic Press: Cambridge, MA, USA, 1974. [Google Scholar]
  27. Metzler, R.; Klafter, J. The restaurant at the end of the random walk: Recent developments in the description of anomalous transport by fractional dynamics. J. Phys. A 2004, 37, R161–R208. [Google Scholar] [CrossRef]
  28. Hu, Z.; Zhang, L. Implicit compact difference schemes the fractional cable equation. Appl. Math. Model 2012, 36, 4027–4043. [Google Scholar] [CrossRef]
  29. Gang, G.; Kun, L.; Yuhui, W. Exact solutions of a modified fractional diffusion equation in the finite and semi-infinite domains. Physica A 2015, 417, 193–201. [Google Scholar]
  30. Jumarie, G. Modified Riemann-Liouville derivative and fractional Taylor series of non-differentiable functions further results. Comput. Math. Appl. 2006, 51, 1367–1376. [Google Scholar]
  31. Podulbny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  32. Ali, U.; Abdullah, F.A. Explicit Saul’yev finite difference approximation for two-dimensional fractional sub-diffusion equation. In Proceedings of the AIP Conference Proceedings, Pahang, Malaysia, 27–29 August 2018; Volume 1, p. 020111. [Google Scholar]
  33. Liu, F.; Yang, Q.; Turner, I. Two new implicit numerical methods for the fractional cable equation. J. Comput. Nonlinear Dyn. 2011, 6, 1–7. [Google Scholar] [CrossRef]
  34. Chen, C.M.; Liu, F.; Turner, I.; Anh, V. Numerical methods with fourth-order spatial accuracy for variable-order nonlinear Stocks’ first problem for a heated generalized second grade fluid. Comput. Math. Appl. 2011, 62, 971–986. [Google Scholar] [CrossRef] [Green Version]
  35. Liu, Q.; Liu, F.; Turner, I.; Anh, V. Finite element approximation for a modified anomalous subdiffusion equation. Appl. Math. Model 2011, 35, 4103–4116. [Google Scholar] [CrossRef]
  36. Oliveira, E.C.; Machado, J.A.T. A Review of Definitions for Fractional Derivatives and Integral. Math. Probl. Eng. 2014, 2014, 238459. [Google Scholar] [CrossRef] [Green Version]
  37. Ali, U. Numerical Solutions for two Dimensional Time-Fractional Differential Sub-Diffusion Equation. Ph.D. Thesis, University Sains Malaysia, Penang, Malaysia, 2019; pp. 1–200. [Google Scholar]
Figure 1. Exact solution.
Figure 1. Exact solution.
Symmetry 12 00691 g001
Figure 2. Approximate solution.
Figure 2. Approximate solution.
Symmetry 12 00691 g002
Figure 3. Comparison of the compact difference scheme in Equation (35) and the exact solution in Equation (38) at α = 0.25 , β = 0.45 , T = 1 , y = 0.1 , and N = 160 .
Figure 3. Comparison of the compact difference scheme in Equation (35) and the exact solution in Equation (38) at α = 0.25 , β = 0.45 , T = 1 , y = 0.1 , and N = 160 .
Symmetry 12 00691 g003
Table 1. Comparison of fourth-order implicit difference scheme for Equation (35), the error E at T = 1.0 , Δ x = Δ y = 1 10 , α = 0.25 , β = 0.45 .
Table 1. Comparison of fourth-order implicit difference scheme for Equation (35), the error E at T = 1.0 , Δ x = Δ y = 1 10 , α = 0.25 , β = 0.45 .
τ [17][1]Scheme C 1 -Order of Scheme
1 / 10 1.5729 × 10 2 5.1782 × 10 3 4.9806 × 10 3 -
1 / 20 7.8976 × 10 3 2.4390 × 10 3 2.2450 × 10 3 1.1497
1 / 40 3.9560 × 10 3 1.1772 × 10 3 9.8504 × 10 4 1.1885
1 / 80 1.9794 × 10 3 5.9755 × 10 4 4.0617 × 10 4 1.2782
1 / 160 9.8995 × 10 4 3.3033 × 10 4 1.3919 × 10 4 1.5451
Table 2. The order of convergence of fourth-order implicit difference scheme for Equation (35), the error E at T = 1.0 , α = 0.2 , β = 0.8 .
Table 2. The order of convergence of fourth-order implicit difference scheme for Equation (35), the error E at T = 1.0 , α = 0.2 , β = 0.8 .
Scheme C 2 -Order
τ = h x = h y = 1 / 4 1.6406 × 10 2 -
τ = 1 / 8 , h x = h y = 1 / 64 7.3344 × 10 4 4.4838
τ = h x = h y = 1 / 8 8.2233 × 10 3 -
τ = 1 / 16 , h x = h y = 1 / 128 3.8694 × 10 4 4.4099
Table 3. Comparison of fourth-order implicit difference scheme for Equation (39) with the exact solution (43), the error E at T = 1.0 , Δ x = Δ y = 1 10 , α = 0.25 , β = 0.45 .
Table 3. Comparison of fourth-order implicit difference scheme for Equation (39) with the exact solution (43), the error E at T = 1.0 , Δ x = Δ y = 1 10 , α = 0.25 , β = 0.45 .
γ ( Δ x ) 2 = ( Δ y ) 2 = τ = 1 / 16 ( Δ x ) 2 = ( Δ y ) 2 = τ = 1 / 64 ( Δ x ) 2 = ( Δ y ) 2 = τ = 1 / 144
0.5 5.0794 × 10 3 1.4043 × 10 3 6.2458 × 10 4
0.6 5.9521 × 10 3 1.6387 × 10 3 6.8873 × 10 4
0.7 6.7049 × 10 3 1.8575 × 10 3 7.5239 × 10 4
0.8 7.4635 × 10 3 2.1414 × 10 3 8.1996 × 10 4
0.9 8.2610 × 10 3 2.3888 × 10 3 8.9478 × 10 4

Share and Cite

MDPI and ACS Style

Ali, U.; Sohail, M.; Usman, M.; Abdullah, F.A.; Khan, I.; Nisar, K.S. Fourth-Order Difference Approximation for Time-Fractional Modified Sub-Diffusion Equation. Symmetry 2020, 12, 691. https://doi.org/10.3390/sym12050691

AMA Style

Ali U, Sohail M, Usman M, Abdullah FA, Khan I, Nisar KS. Fourth-Order Difference Approximation for Time-Fractional Modified Sub-Diffusion Equation. Symmetry. 2020; 12(5):691. https://doi.org/10.3390/sym12050691

Chicago/Turabian Style

Ali, Umair, Muhammad Sohail, Muhammad Usman, Farah Aini Abdullah, Ilyas Khan, and Kottakkaran Sooppy Nisar. 2020. "Fourth-Order Difference Approximation for Time-Fractional Modified Sub-Diffusion Equation" Symmetry 12, no. 5: 691. https://doi.org/10.3390/sym12050691

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