Next Article in Journal
Asymmetric Lévy Flights Are More Efficient in Random Search
Previous Article in Journal
Novel Approaches for Differentiable Convex Functions via the Proportional Caputo-Hybrid Operators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Parallelized Computation Method of HASC-N Difference Scheme for Inhomogeneous Time Fractional Fisher Equation

Institute of Information and Computation, School of Mathematics and Physics, North China Electric Power University, Beijing 102206, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(5), 259; https://doi.org/10.3390/fractalfract6050259
Submission received: 30 March 2022 / Revised: 1 May 2022 / Accepted: 2 May 2022 / Published: 7 May 2022

Abstract

:
The fractional Fisher equation has a wide range of applications in many engineering fields. The rapid numerical methods for fractional Fisher equation have momentous scientific meaning and engineering applied value. A parallelized computation method for inhomogeneous time-fractional Fisher equation (TFFE) is proposed. The main idea is to construct the hybrid alternating segment Crank-Nicolson (HASC-N) difference scheme based on alternating segment difference technology, using the classical explicit scheme and classical implicit scheme combined with Crank-Nicolson (C-N) scheme. The unique existence, unconditional stability and convergence are proved theoretically. Numerical tests show that the HASC-N difference scheme is unconditionally stable. The HASC-N difference scheme converges to O ( τ 2 α + h 2 ) under strong regularity and O ( τ α + h 2 ) under weak regularity of fractional derivative discontinuity. The HASC-N difference scheme has high precision and distinct parallel computing characteristics, which is efficient for solving inhomogeneous TFFE.

1. Introduction

The time fractional Fisher equation (TFFE) is a nonlinear physical model with linear diffusion and nonlinear growth. Derived from population dynamics, chemical dynamics and other fields, it describes phenomena such as mutant gene reproduction, nonlinear evolution of population and autocatalytic chemical reaction [1,2]. Exact solutions of TFFE are difficult to be given explicitly and most of them contain special functions, such as the multivariable Mittag-Leffler function [3,4,5]. In the past two decades, with the deepening of the application of TFFE, the rapid numerical solution for TFFE has become an urgent research work [6,7].
At present, the finite difference method is still the more widely used and mature numerical method for solving TFFE. The finite difference method can achieve the precision and stability of simulation requirements well [8]. Zhang et al. (2014) [9] constructed a fully discrete scheme of TFFE by combining the finite difference method and locally discontinuous Galerkin finite element method, and discussed the stability and error estimation of the method. Alquran et al. (2015) [10] numerically solved the TFFE based on the self-collocation method and finite difference method, and analyzed the analytical and numerical solutions of the equation. Mejía and Piedrahita (2019) [11] proposed an implicit finite difference scheme for approximating TFFE with variable coefficients, and the numerical results verified the correctness of the theoretical analysis. There are also many research results on other numerical solutions of TFFE [12,13,14,15], but in most of the above numerical methods, computational efficiency has not been paid enough attention.
Due to the improvement of cluster technology and the increasing number of CPU cores, the parallelized numerical method has become one of the main methods for fast computing [16,17]. For the past few years, parallel computing has been widely used in the field of rapid numerical solutions for fractional partial differential equations (FPDEs). At present, there are three kinds of parallel algorithms using fractional differential equations: algebraic parallel algorithm, Parareal algorithm and parallel difference scheme.
Based on the algebraic parallel algorithm, Gong et al. (2013) [18] came up with a parallelized calculation method of explicit difference schemes for fractional reaction-diffusion equations, which was mainly used for parallel calculation of matrix and vector in algebraic matrix equation. Sweilam et al. (2014) [19] proposed an algebraic parallel algorithm for the time-fractional parabolic equation. This method solved the algebraic equation matrix after discrete-time in parallel. Biala and Khaliq (2018) [20] developed a C-N scheme similar to integer-order parabolic equations for nonlinear spatio-temporal fractional parabolic equations, and used the precursor-correction method respectively in MPI, OpenMP and Hybrid Version.
Using Parareal algorithm, Fu and Wang (2019) [21] constructed a Parareal algorithm to solve the space-time FPDE that models an anomalous diffusion process in a one-dimensional tube. The numerical advantages of the traditional Parareal algorithm were well preserved in this method. Yue et al. (2019) [22] proposed a multi-grid time reduction (MGRIT) algorithm based on time-varying time-grid propagators for two-dimensional fractional diffusion equations, and presented the two-level convergence theory of the algorithm. Liu et al. (2020) [23] proposed the finite volume method for time-varying fractional parabolic equations, and parallelized it with the parallel-In-time method to improve the computational efficiency of the finite volume method. Based on the Parareal method, Lorin (2020) [24] constructed the Parareal-Gorenflo algorithm for space-time FPDEs, and the spatial parallelization of this method relied on the parallelization of Riesz derivative and fast Fourier transform.
For the study of parallel difference schemes, Wang et al. (2016) [25] parallelized the implicit difference scheme of the Caputo fractional reaction-diffusion equation, and changed the serial algorithm to parallel as far as possible without changing the original serial difference scheme, to reasonably allocate computing tasks. Yang and Wu (2020) [26] proposed a parallel nature difference method for a multi-term time-fractional diffusion equation and proved that the method was unconditionally stable and convergent through theoretical analysis. Numerical experiments showed that the scheme proposed by Yang and Wu is an efficient scheme for the multi-term time-fractional diffusion equation.
To solve the problem of large computation of fractional Fisher parabolic equation, we explore the parallelization of the difference scheme for the inhomogeneous TFFE. A new parallelized computation method is proposed by using an alternate technique appropriately, which ensures the unconditional stability and spatial convergence order O ( h 2 ) of the new algorithm, and is easy to be used in many types of parallel machines.

2. HASC-N Difference Scheme for Inhomogeneous TFFE

2.1. Inhomogeneous Time Fractional Fisher Equation

Consider the inhomogeneous TFFE as follows [27,28,29]:
α u ( x , t ) t α = 2 u ( x , t ) x 2 + u ( x , t ) ( 1 u ( x , t ) ) + g ( x , t ) , ( x , t ) 0 , L × 0 , T , u ( x , 0 ) = ϕ ( x ) , x 0 , L , u ( 0 , t ) = φ 1 ( t ) , u ( L , t ) = φ 2 ( t ) , t 0 , T ,
where ϕ x , φ 1 t , φ 2 t are the given functions with suitable smoothness. The nonlinear source term u ( 1 u ) is a nonlinear function, 0 < α 1 .
For brevity, let f u x , t , x , t = u 1 u be Lipschitz continuous with respect to u, and there exists a Lipschitz coefficient l such that f ( u 1 ) f ( u 2 ) l u 1 u 2 . The inhomogeneous term g x , t is a known function independent of u. α u ( x , t ) t α is the fractional derivative in the sense of Caputo:
α u ( x , t ) t α = 1 Γ ( 1 α ) 0 t u ( x , s ) τ d s ( t s ) α , 0 < α < 1 ,
where Γ ( ) is the Gamma function. When α = 1 , equation is
u x , t t = 2 u x , t x 2 + u x , t ( 1 u x , t ) .
The above Equation (3) is called the classical Fisher equation in general. As α tends to 1, according to the conclusion in reference [30,31], solution u x , t of TFFE tends to u ˜ x , t ( u ˜ x , t is the solution of the classical Fisher equation).

2.2. Construction of HASC-N Difference Scheme for Inhomogeneous TFFE

To construct the HASC-N difference scheme of inhomogeneous TFFE (1), the solution region Ω = ( x , t ) | 0 x L , 0 t T is meshed. Take the space step h = L M and time step τ = T N , where M and N are positive integers. Thus, x j = j h ( j = 1 , 2 , , M ) , M h = L ,   t k = k τ ( k = 1 , 2 , , N ) ,   N τ = T and the grid node is ( x j , t k ) . Define u j k = u ( x j , t k ) , f j k = f ( u ( x j , t k ) , x j , t k ) , g j k = g ( x j , t k ) .
Lemma 1
([7,32]). Suppose 0 < α < 1 , let y C 2 [ 0 , t n ] . Then
1 Γ ( 1 α ) 0 t n y ( ξ ) d ξ ( t ξ ) α τ α Γ ( 2 α ) [ y n k = 1 n 1 ( l n k 1 l n k ) y k l n 1 y 0 ] 1 Γ ( 2 α ) [ 1 α 12 + 2 2 α 2 α ( 1 + 2 α ) ] max 0 t t n y ( t ) τ 2 α ,
where l i = ( i + 1 ) ( 1 α ) i ( 1 α ) , i = 0 , 1 , 2 , , N .
The discrete formula is defined by Lemma 1:
D t α u ( x j , t k + 1 ) = τ α Γ ( 2 α ) l 0 u ( x j , t k + 1 ) i = 1 k ( l i 1 l i ) u ( x j , t k i + 1 ) l k u ( x j , t 0 ) .
The method of processing nonlinear source term f ( u ) is derived from references [33,34]:
f j k = 2 f j k 1 f j k 2 + O ( τ 2 ) .
Define the space second derivative discrete formula:
δ x 2 u j k : = 1 h 2 u j 1 k 2 u j k + u j + 1 k ,
δ x 2 u j k + 1 : = 1 h 2 u j 1 k + 1 2 u j k + 1 + u j + 1 k + 1 ,
D u j k : = 1 2 h 2 u j 1 k + 1 2 u j k + 1 + u j + 1 k + 1 + u j 1 k 2 u j k + u j + 1 k .
Three difference schemes are obtained:
Classical explicit scheme,
D t α u ( x j , t k + 1 ) = δ x 2 u j k + f j k + g j k .
Classical implicit scheme,
D t α u ( x j , t k + 1 ) = δ x 2 u j k + 1 + f j k + 1 + g j k + 1 .
Classical C-N scheme,
D t α u ( x j , t k + 1 ) = D u j k + f j k + 1 2 + g j k + 1 2 .
Further collate the above three difference schemes, we get
u j k + 1 = a u j 1 k + ( b 1 2 a ) u j k + a u j + 1 k + i = 2 k b i u j k i + 1 + l k u j 0 + c f j k + c g j k ,
a u j 1 k + 1 + ( 1 + 2 a ) u j k + 1 a u j + 1 k + 1 = b 1 u j k + i = 2 k b i u j k i + 1 + l k u j 0 + c f j k + 1 + c g j k + 1 ,
a 2 u j 1 k + 1 + ( 1 + a ) u j k + 1 a 2 u j + 1 k + 1 = a 2 u j 1 k + ( b 1 a ) u j k + a 2 u j + 1 k + i = 2 k b i u j k i + 1 + l k u j 0 + c f j k + 1 2 + c g j k + 1 2 ,
where b j = l j 1 l j , c = τ α Γ ( 2 α ) , a = c h 2 .
According to the thought in references [35,36], the HASC-N difference scheme for inhomogeneous TFFE (1) is constructed:
M + 1 points are taken at each time layer, except for the first point and the M + 1 point on the boundary, the remaining M 1 points to be calculated at the same layer are divided into B segments (B is odd without losing generality). If there are n points in each segment, n B = M 1 (n and B are positive integers and n 3 , B 3 ). The classical explicit scheme and classical implicit scheme are used alternately at the boundary points of two adjacent time layers. At the inner boundary points of two adjacent time layers, the classical explicit scheme and the classical implicit scheme are used alternately. The C-N scheme is used at the remaining points of two adjacent time layers. ● is the classical explicit scheme, ◯ is the classical implicit scheme, and ■ is the classical C-N scheme. HASC-N difference scheme construction principle is shown in Figure 1:
The HASC-N difference scheme for inhomogeneous TFFE (1) can be as follows:
I + A 1 G U k + 1 = b 1 I A 2 G U k + c k + i = 2 k b i U k i + 1 + l k U 0 + c A 1 F k + 1 + c A 2 F k , I + A 2 G U k + 2 = b 1 I A 1 G U k + 1 + c k + 1 + i = 2 k + 1 b i U k i + 2 + l k + 1 U 0 + c A 2 F k + 2 + c A 1 F k + 1 , k = 0 , 2 , 4
where
G = 2 a a a 2 a a a 2 a a a 2 a a a 2 a M 1 × M 1
A 1 = θ 1 θ 2 θ 3 θ M 2 θ M 1 M 1 × M 1
θ j = 0 , j = n , 2 n , , ( B 1 ) n , 1 , j = n + 1 , 2 n + 1 , , ( B 1 ) n + 1 , 1 2 , e l s e w h e r e . , U k = u 1 k , u 2 k , , u M 1 k T , c k = a u 0 k , 0 , , 0 , a u M k T , f k = ( f 1 k , f 2 k , , f M 1 k ) T , g k = ( g 1 k , g 2 k , , g M 1 k ) T , F k = f k + g k , A 2 = I A 1 , I is identity matrix.

3. Existence and Uniqueness of Solution to HASC-N Difference Scheme for Inhomogeneous TFFE

Lemma 2.
The matrices A 1 G and A 2 G in HASC-N difference scheme (16) are non-negative definite matrices.
Proof. 
According to
A 1 G = 2 θ 1 a θ 1 a θ 2 a 2 θ 2 a θ 2 a θ 3 a 2 θ 3 a θ 3 a θ M 2 a 2 θ M 2 a θ M 2 a θ M 1 a 2 θ M 1 a ( M 1 ) × ( M 1 ) ,
we get that A 1 G is the diagonally dominant tridiagonal matrix, and the main diagonal elements are non-negative real numbers, So A 1 G is a non-negative definite matrix. Similarly, A 2 G is a non-negative definite matrix. Lemma 2 is proved. □
Theorem 1.
The solution of HASC-N difference scheme (16) for inhomogeneous TFFE (1) is existing and unique.
Proof. 
According to Lemma 2, the inverse matrices I + A 1 G 1 and I + A 2 G 1 of I + A 1 G and I + A 2 G exist, the HASC-N difference scheme (16) has a unique solution. Therefore, the Theorem 1 is proved. □

4. Stability of HASC-N Difference Scheme for Inhomogeneous TFFE

Theorem 2.
The HASC-N difference scheme (16) for inhomogeneous TFFE (1) is unconditionally stable.
Proof. 
Assume that u j k is the HASC-N difference scheme solution for inhomogeneous TFFE (1), U ¯ j k is the approximate solution of HASC-N difference scheme for inhomogeneous TFFE (1). Error ε j k is defined as ε j k = U ¯ j k u j k , let ε 0 k = ε M k = 0 , E k = ε 1 k , ε 2 k , , ε M 1 k , k = 0 , 1 , 2 , , N . Substitute the approximate solution U ¯ j k of HASC-N difference scheme and the HASC-N difference scheme solution u j k into scheme (16), respectively, to get two equations, and make the difference between the two equations, we get
I + A 1 G E k + 1 = b 1 I A 2 G E k + i = 2 k b i E k i + 1 + l k E 0 + c A 1 ( f ¯ k + 1 f k + 1 ) + c A 2 ( f ¯ k f k ) , I + A 2 G E k + 2 = b 1 I A 1 G E k + 1 + i = 2 k + 1 b i E k i + 2 + l k + 1 E 0 + c A 2 ( f ¯ k + 2 f k + 2 ) + c A 1 ( f ¯ k + 1 f k + 1 ) , k = 0 , 2 , 4
where f k = ( f 1 k , f 2 k , , f M k ) , f ¯ k = ( f ¯ 1 k , f ¯ 2 k , , f ¯ M k ) , f ¯ j k = f ( u ¯ j k , x j , t k ) .
Since f ( u , x , t ) satisfifies the Lipschitz condition l, there is the Lipschitz conditional constant, we have
f ¯ k f k l U ¯ k u k l E k .
Substitute Equation (19) into Equation (18) and we get
I l c A 1 + A 1 G E k + 1 b 1 I + l c A 2 A 2 G E k + i = 2 k b i E k i + 1 + l k E 0 , I l c A 2 + A 2 G E k + 2 b 1 I + l c A 1 A 1 G E k + 1 + i = 2 k + 1 b i E k i + 2 + l k + 1 E 0 . k = 0 , 2 , 4
So for simplicity, let α θ j = l c θ j , β θ j = l c 1 θ j , we have
( 1 α θ j ) I + A 1 G E k + 1 ( b 1 + β θ j ) I A 2 G E k + i = 2 k b i E k i + 1 + l k E 0 , ( 1 β θ j ) I + A 2 G E k + 2 ( b 1 + α θ j ) I A 1 G E k + 1 + i = 2 k + 1 b i E k i + 2 + l k + 1 E 0 . k = 0 , 2 , 4
Define norm U k = U k = max 1 j ( M 1 ) u j k . Known by the definition of matrices A 1 , A 2 , G , A 1 G and A 2 G 5 are non-negative definite matrices, and they have different non-negative characteristic values. Let the characteristic value of A 1 G be λ j and the characteristic value of A 2 G be γ j , λ j H 1 , γ j H 2 , H 1 and H 2 are constants, γ j = λ j + K j , K j is constant, j = 1 , 2 , , M 1 .
According to reference [14], there is an unequal relationship between the time process T of the TFFE and Lipschitz coefficient l. Assume that in the unequal relation between time process T and Lipschitz coefficient l, the following inequality holds: 2 l 1 1 α θ j min 1 , λ j , 2 l 1 1 + K j β θ j min 1 + K j , γ j , where γ j = λ j + K j . This is bound to affect the length of time process T, however, in order to ensure the stability of HASC-N scheme, the following proofs are carried out under the premise that the above assumption is true. The results of numerical tests also confirm the feasibility of this assumption.
We will prove E k E 0 by mathematical induction.
When k = 0 , namely ( 1 α θ j ) I + A 1 G E 1 ( b 1 + β θ j ) I A 2 G E 0 , ( 1 β θ j ) I + A 2 G E 2 ( b 1 + α θ j ) I A 1 G E 1 + l 1 E 0 .
Firstly, We discuss ( 1 α θ j ) I + A 1 G E 1 ( b 1 + β θ j ) I A 2 G E 0 . Solve for E 1 and take the norm of both sides, we get
E 1 ( 1 α θ j ) I + A 1 G 1 ( b 1 + β θ j ) I A 2 G E 0 max ( b 1 + β θ j ) γ j ( 1 α θ j ) + λ j E 0 .
Case 1, b 1 + β θ j > γ j ,
max ( b 1 + β θ j ) γ j ( 1 α θ j ) + λ j ( b 1 + β θ j ) γ j ( 1 α θ j ) + λ j 1 ( γ j β θ j ) 1 + ( λ j α θ j ) 1 .
Case 2, b 1 + β θ j γ j ,
max ( b 1 + β θ j ) γ j ( 1 α θ j ) + λ j γ j ( b 1 + β θ j ) λ j + ( 1 α θ j ) 1 .
According to (22) and (23), we have max ( b 1 + β θ j ) λ j ( 1 α θ j ) + λ j 1 , E 1 E 0 .
Secondly, we discuss ( 1 β θ j ) I + A 2 G E 2 ( b 1 + α θ j ) I A 1 G E 1 + l 1 E 0 . Solve for E 2 and take the norm of both sides, we get
E 2 ( 1 β θ j ) I + A 2 G 1 [ ( ( b 1 + α θ j ) I A 1 G ) E 1 + l 1 E 0 ] max ( b 1 + α θ j ) λ j + l 1 ( 1 β θ j ) + γ j E 0 .
Case 1, b 1 + α θ j > λ j ,
max ( b 1 + α θ j ) λ j + l 1 ( 1 β θ j ) + γ j ( b 1 + α θ j ) λ j + l 1 ( 1 β θ j ) + γ j = 1 ( λ j α θ j ) 1 + ( γ j β θ j ) 1 .
Case 2, b 1 + α θ j λ j ,
max ( b 1 + α θ j ) λ j + l 1 ( 1 β θ j ) + γ j λ j ( b 1 + α θ j ) + l 1 ( 1 β θ j ) + γ j = λ j ( 1 2 l 1 ) + α θ j λ j + ( 1 + K j β θ j ) 1 .
According to (24) and (25), we have max ( b 1 + α θ j ) λ j + l 1 ( 1 β θ j ) + γ j 1 , E 2 E 0 .
Finally, assuming that the previous layers are all true, namely E k E 0 .
When the time layer is layer k + 1 and layer k + 2 , ( 1 α θ j ) I + A 1 G E k + 1 ( b 1 + β θ j ) I A 2 G E k + i = 2 k b i E k i + 1 + l k E 0 , ( 1 β θ j ) I + A 2 G E k + 2 ( b 1 + α θ j ) I A 1 G E k + 1 + i = 2 k + 1 b i E k i + 2 + l k + 1 E 0 .
Solve for E k + 1 and E k + 2 and take the norm of both sides, we get
E k + 1 ( 1 α θ j ) I + A 1 G 1 [ ( b 1 + β θ j ) I A 2 G E k + i = 2 k b i E k i + 1 + l k E 0 ] , E k + 2 ( 1 β θ j ) I + A 2 G 1 [ ( b 1 + α θ j ) I A 1 G E k + 1 + i = 2 k + 1 b i E k i + 2 + l k + 1 E 0 ] .
According to E k E 0 and b j = l j 1 l j , we get
E k + 1 = ( 1 α θ j ) I + A 1 G 1 [ ( b 1 + β θ j ) I A 2 G E k + i = 2 k b i E k i + 1 + l k E 0 ] = ( 1 α θ j ) I + A 1 G 1 [ ( b 1 + β θ j ) I A 2 G E k + b 2 E k 1 + + b k E 1 + l k E 0 ] ( 1 α θ j ) I + A 1 G 1 [ ( b 1 + β θ j ) I A 2 G E 0 + b 2 E 0 + + b k E 0 + l k E 0 ] = ( 1 α θ j ) I + A 1 G 1 [ ( b 1 + β θ j ) I A 2 G E 0 + l 1 E 0 ] ( 1 α θ j ) I + A 1 G 1 [ ( b 1 + β θ j ) I A 2 G + l 1 ] E 0 max ( b 1 + β θ j ) γ j + l 1 ( 1 α θ j ) + λ j E 0 .
Case 1, b 1 + α θ j > λ j ,
max ( b 1 + β θ j ) γ j + l 1 ( 1 α θ j ) + λ j ( b 1 + β θ j ) γ j + l 1 ( 1 α θ j ) + λ j = 1 ( γ j β θ j ) 1 + ( λ j α θ j ) 1 .
Case 2, b 1 + α θ j λ j ,
max ( b 1 + β θ j ) γ j + l 1 ( 1 α θ j ) + λ j γ j ( b 1 + β θ j ) + l 1 ( 1 α θ j ) + λ j = λ j ( 1 2 l 1 ) + β θ j K j λ j + ( 1 α θ j ) 1 .
According to (26) and (27), we have max ( b 1 + β θ j ) γ j + l 1 ( 1 α θ j ) + λ j 1 , E k + 1 E 0 .
Similarly, according to E k + 1 E 0 and b j = l j 1 l j , we get
E k + 2 = ( 1 β θ j ) I + A 2 G 1 [ ( b 1 + α θ j ) I A 1 G E k + 1 + i = 2 k + 1 b i E k i + 2 + l k + 1 E 0 ] = ( 1 β θ j ) I + A 2 G 1 [ ( b 1 + α θ j ) I A 1 G E k + 1 + b 2 E k + + b k + 1 E 1 + l k + 1 E 0 ] ( 1 β θ j ) I + A 2 G 1 [ ( b 1 + α θ j ) I A 1 G E 0 + b 2 E 0 + + b k + 1 E 0 + l k + 1 E 0 ] = ( 1 β θ j ) I + A 2 G 1 [ ( b 1 + α θ j ) I A 1 G E 0 + l 1 E 0 ] ( 1 β θ j ) I + A 2 G 1 [ ( b 1 + α θ j ) I A 1 G + l 1 ] E 0 max ( b 1 + α θ j ) λ j + l 1 ( 1 β θ j ) + γ j E 0 E 0 .
In summary, stability is proved. □

5. Convergence of HASC-N Difference Scheme for Inhomogeneous TFFE

Lemma 3
([6]). Suppose 0 < α < 1 , let y C 2 [ 0 , t n + 1 ] . Then we have
α + 1 y ( t n + 1 ) t α + 1 = 1 Γ ( 1 α ) 0 t 2 y ( ξ ) ξ 2 d ξ t n + 1 ξ α = 1 Γ ( 1 α ) j = 0 n j τ j + 1 τ 2 y ( ξ ) ξ 2 d ξ t n + 1 ξ α 1 Γ ( 1 α ) max 0 t t n + 1 2 y ( t ) t 2 j = 0 n j τ j + 1 τ d ξ t n + 1 ξ α = C y Γ ( 1 α ) j = 0 n j τ j + 1 τ d ξ t n + 1 ξ α . = C y τ 1 α Γ ( 2 α ) j = 0 n j τ j + 1 τ n + 1 j 1 α n j 1 α n + 1 1 α C y Γ ( 2 α ) τ 1 α , where C y = max 0 t t n + 1 2 y t t 2 .
The solution of the inhomogeneous TFFE (1) satisfies the strong regularity condition as follows,
γ u t γ C 0 , L × 0 , T , δ u x δ C 0 , L × 0 , T ,
where γ 1 , 2 and δ 0 , 1 , 2 , 3 , 4 .
Consider the explicit scheme on the time layer k + 1 ,
D t α u ( x j , t k + 1 ) = 1 h 2 ( u j 1 k 2 u j k + u j + 1 k ) + F j k ,
and the implicit scheme on the time layer k + 2 ,
D t α u ( x j , t k + 2 ) = 1 h 2 ( u j 1 k + 2 2 u j k + 2 + u j + 1 k + 2 ) + F j k + 2 .
Taylor expansion is performed at u j k + 1 for scheme (29) and scheme (30) to obtain truncation error,
R 1 ( τ , h ) = α u ( x j , t k + 1 ) t α u x x + τ u x x t 1 12 h 2 u x x x x τ F u t + O ( τ 2 α + h 2 ) ,
R 2 ( τ , h ) = α u ( x j , t k + 1 ) t α + τ α + 1 u ( x j , t k + 1 ) t α + 1 u x x τ u x x t 1 12 h 2 u x x x x + τ F u t + O ( τ 2 α + h 2 ) .
Consider the C-N scheme on the time layer k + 1 ,
D t α u ( x j , t k + 1 ) = 1 2 h 2 ( u j 1 k + 1 2 u j k + 1 + u j + 1 k + 1 + u j 1 k 2 u j k + u j + 1 k ) + 1 2 F j k + 1 + F j k ,
and the C-N scheme on the time layer k + 2 ,
D t α u ( x j , t k + 2 ) = 1 2 h 2 ( u j 1 k + 2 2 u j k + 2 + u j + 1 k + 2 + u j 1 k + 1 2 u j k + 1 + u j + 1 k + 1 ) + 1 2 F j k + 2 + F j k + 1 .
Taylor expansion is performed at u j k + 1 for scheme (33) and scheme (34) to obtain truncation error,
R 3 ( τ , h ) = α u ( x j , t k + 1 ) t α u x x + τ 2 u x x t 1 12 h 2 u x x x x τ 2 F u t + O ( τ 2 α + h 2 ) ,
R 4 ( τ , h ) = α u ( x j , t k + 1 ) t α + τ α + 1 u ( x j , t k + 1 ) t α + 1 u x x τ 2 u x x t 1 12 h 2 u x x x x + τ 2 F u t + O ( τ 2 α + h 2 ) .
According to Lemma 1, the calculation precision of α u ( x j , t k + 1 ) t α is O ( τ 2 α ) , the calculation precision of τ α + 1 u ( x j , t k + 1 ) t α + 1 is also O ( τ 2 α ) according to Lemma 3. By using explicit and implicit schemes alternately at the inner boundary points of adjacent time layers, two basic error components with opposite signs are generated, and the two partially cancel each other, so as to obtain ideal calculation precision.
Add (31) and (32) to get
R 1 ( τ , h ) + R 2 ( τ , h ) = 2 α u ( x j , t k + 1 ) t α + τ α + 1 u ( x j , t k + 1 ) t α + 1 2 u x x 1 6 h 2 u x x x x + O ( τ 2 α + h 2 ) .
At the inner boundary points, the calculation precision is O ( τ 2 α + h 2 ) . Similarly, C-N scheme is alternately used at interior points of adjacent time layers, (35) and (36) are added to obtain R 3 ( τ , h ) + R 4 ( τ , h ) = 2 α u ( x j , t k + 1 ) t α + τ α + 1 u ( x j , t k + 1 ) t α + 1 2 u x x 1 6 h 2 u x x x x + O ( τ 2 α + h 2 ) .
So the precision at the interior points is also O ( τ 2 α + h 2 ) .
Theorem 3.
Assuming that the solution of Equation (1) satisfies the strong regularity condition (28), the HASC-N difference scheme (16) for inhomogeneous TFFE (1) is convergent, and e n C ( τ 2 α + h 2 ) , n = 1 , 2 , , N , C > 0 .
Proof. 
Let U j k = U ( x j , t k ) be the exact solution of inhomogeneous TFFE (1) at t = t k , x = x j under strong regularity. Define e j k = U j k u j k ,   1 j M 1 ,   e 0 k = e M k = 0 , e k = ( e 1 k , , e M 1 k ) ,   e 0 = 0 .
Substitute the exact solution U j k and the HASC-N difference scheme solution u j k into scheme (16), respectively, to get two equations, and make the difference between the two equations, we get
I + A 1 G e k + 1 = b 1 I A 2 G e k + i = 2 k b i e k i + 1 + c A 1 ( f ¯ k + 1 f k + 1 ) + c A 2 ( f ¯ k f k ) + R ˜ k + 1 , I + A 2 G e k + 2 = b 1 I A 1 G e k + 1 + i = 2 k + 1 b i e k i + 2 + c A 2 ( f ¯ k + 2 f k + 2 ) + c A 1 ( f ¯ k + 1 f k + 1 ) + R ˜ k + 2 , k = 0 , 2 , 4
where R ˜ k + 1 = τ α O ( τ 2 α + h 2 ) , R ˜ k + 1 C 1 τ α ( τ 2 α + h 2 ) = C 1 ( τ 2 + τ α h 2 ) , C 1 is a real constant.
Similar to stability analysis, convergence is studied:
When k = 0 , ( 1 α θ j ) I + A 1 G e 1 ( b 1 + β θ j ) I A 2 G e 0 + R ˜ 1 , ( 1 β θ j ) I + A 2 G e 2 ( b 1 + α θ j ) I A 1 G e 1 + R ˜ 2 .
Solve for e 1 and e 2 and take the norm, we get
e 1 ( 1 α θ j ) I + A 1 G 1 [ ( b 1 + β θ j ) I A 2 G e 0 + R ˜ 1 ] , e 2 ( 1 β θ j ) I + A 2 G 1 [ ( b 1 + α θ j ) I A 1 G e 1 + R ˜ 2 ] .
Firstly, e 1 ( 1 α θ j ) I + A 1 G 1 [ ( b 1 + β θ j ) I A 2 G e 0 + R ˜ 1 ] , and, e 0 = 0 , we get e 1 ( 1 α θ j ) I + A 1 G 1 R ˜ 1 max 1 1 + λ j α θ j R ˜ 1 l 0 1 C 1 ( τ 2 + τ α h 2 ) .
Secondly, e 2 ( 1 β θ j ) I + A 2 G 1 [ ( b 1 + α θ j ) I A 1 G e 1 + R ˜ 2 ] l 1 1 ( 1 β θ j ) I + A 2 G 1 [ ( b 1 + α θ j ) I A 1 G + l 1 ] || R ˜ 2 || l 1 1 max b 1 + α θ j λ j + l 1 1 β θ j + γ j C 1 ( τ 2 + τ α h 2 ) l 1 1 C 1 ( τ 2 + τ α h 2 ) .
Assuming that e k l k 1 1 C 1 ( τ 2 + τ α h 2 ) is true for all the previouslayers. When the time layer are k + 1 and k + 2 ,
( 1 α θ j ) I + A 1 G e k + 1 ( b 1 + β θ j ) I A 2 G e k + i = 2 k b i e k i + 1 + R ˜ k + 1 , ( 1 β θ j ) I + A 2 G e k + 2 ( b 1 + α θ j ) I A 1 G e k + 1 + i = 2 k + 1 b i e k i + 2 + R ˜ k + 2 .
Solve for e k + 1 and e k + 2 and take the norm, we get
e k + 1 ( 1 α θ j ) I + A 1 G 1 [ ( b 1 + β θ j ) I A 2 G e k + i = 2 k b i e k i + 1 + R ˜ k + 1 ] , e k + 2 ( 1 β θ j ) I + A 2 G 1 [ ( b 1 + α θ j ) I A 1 G e k + 1 + i = 2 k + 1 b i e k i + 2 + R ˜ k + 2 ] .
Then e k + 1 ( 1 α θ j ) I + A 1 G 1 [ ( b 1 + β θ j ) I A 2 G e k + i = 2 k b i e k i + 1 + R ˜ k + 1 ] l k 1 ( 1 α θ j ) I + A 1 G 1 [ ( b 1 + β θ j ) I A 2 G + b 2 + + b k + l k ] || R ˜ k + 1 || l k 1 max ( b 1 + β θ j ) γ j + l 1 ( 1 α θ j ) + λ j C 1 ( τ 2 + τ α h 2 ) l k 1 C 1 ( τ 2 + τ α h 2 ) ,
e k + 2 ( 1 β θ j ) I + A 2 G 1 [ ( b 1 + α θ j ) I A 1 G e k + 1 + i = 2 k + 1 b i e k i + 2 + R ˜ k + 2 ] l k + 1 1 ( 1 β θ j ) I + A 2 G 1 [ ( b 1 + α θ j ) I A 1 G + b 2 + + b k + 1 + l k + 1 ] || R ˜ k + 2 || l k + 1 1 max ( b 1 + α θ j ) λ j + l 1 ( 1 β θ j ) + γ j C 1 ( τ 2 + τ α h 2 ) l k + 1 1 C 1 ( τ 2 + τ α h 2 ) .
In conclusion, we prove that e n l n 1 1 C 1 ( τ 2 + τ α h 2 ) ,   n = 1 , 2 , , N . From lim n l n 1 n α = lim n n α n ( 1 α ) ( n 1 ) ( 1 α ) = lim n n 1 1 ( 1 1 n ) ( 1 α ) = 1 1 α , there exists C = ( n τ ) α 1 α C 1 such that e n l n 1 1 C 1 ( τ 2 + τ α h 2 ) C ( τ 2 α + h 2 ) ,   n = 1 , 2 , , N . Therefore, proof is completed, and the convergence order is O ( τ 2 α + h 2 ) .  □
Remark 1.
The exact solution u ( x , t ) of the inhomogeneous TFFE (1) satisfies the strong regularity condition (28). In this case, the convergence order of HASC-N difference scheme (16) is O ( τ 2 α + h 2 ) . In general, the condition of strong regularity is too harsh. The exact solution of inhomogeneous TFFE (1) cannot meet this requirement under some conditions, such as solving the solution of inhomogeneous TFFE with initial singularities, and the corresponding theoretical analysis has obvious limitations. Nevertheless, the conclusion is significant because Theorem 3 at least rigorously proves the theoretical correctness of HASC-N difference scheme (16) in a certain range.
Remark 2.
The time fractional derivative of the exact solution u ( x , t ) for the inhomogeneous TFFE (1) is a discontinuous function at initial time, namely γ u t γ does not exist in some regions of [ 0 , L ] × [ 0 , T ] , where γ { 0 , 1 , 2 } . In this case, the strong regularity condition (28) cannot be satisfied, resulting in the initial singularity of inhomogeneous TFFE (1). Consider the following two cases:
(1) The partial derivative of the solution u ( x , t ) in the spatial direction satisfies δ u x δ C ( [ 0 , L ] × [ 0 , T ] ) ,   δ { 0 , 1 , 2 , 3 , 4 } . The HASC-N difference scheme (16) converges to O ( τ α ) in the temporal direction (consistent with the conclusions of references [37,38,39]), and it converges to O ( h 2 ) in the spatial direction.
(2) The partial derivative of the solution u ( x , t ) in the spatial direction is a discontinuous function, that is, δ u x δ does not exist in some regions of [ 0 , L ] × [ 0 , T ] , where δ { 0 , 1 , 2 , 3 , 4 } . In this case, the local truncation error of HASC-N difference scheme (16) lacks clear overall control. Even if the loose discrete L 2 norm is used as a measure, the order of local truncation error is not clear [7,40]. Therefore, the analytic path of spatial and temporal convergence order based on strong regularity condition is no longer effective.

6. Numerical Tests

The numerical tests are based on Intel Core I5-5200 CPU @2.20 GHz, dual-core processor, and carried out in MatlabR2018b environment. Numerical tests verify the correctness of the above theoretical analysis.
Example 1
([41]). Consider the inhomogeneous TFFE with a smooth solution:
α u ( x , t ) t α = 2 u ( x , t ) x 2 + u ( x , t ) ( 1 u ( x , t ) ) + g 1 x , t , ( x , t ) 0 , 1 × 0 , 1 , u ( x , 0 ) = 0 , x 0 , 1 , u ( 0 , t ) = u ( 1 , t ) = 0 , t 0 , 1 .
where g 1 x , t = 24 t ( 4 α ) s i n ( 2 π x ) / Γ ( 5 α ) + 4 π 2 t 4 s i n ( 2 π x ) t 4 s i n ( 2 π x ) ( 1 t 4 s i n ( 2 π x ) ) ,   0 < α 1 . Exact solution of the inhomogeneous TFFE (38): u ( x , t ) = t 4 s i n ( 2 π x ) .
When α = 0.7 ,   N = 100 ,   M = 71 , the exact solution surface, C-N scheme solution surface and HASC-N scheme solution surface are as follows:
According to Figure 2, Figure 3 and Figure 4, the surfaces of the two schemes are consistent with those of the exact solution and the surface of the HASC-N difference scheme is smooth. It is shown below that when α is of different values, the HASC-N scheme solution is compared with the exact solution at t = 0.5 . The HASC-N scheme solution approximates the exact solution well, and the calculation results are shown in Table 1:
Let U j k be the exact solution, and U ˜ j k be the solution of C-N scheme and HASC-N difference scheme. Define the Sum of Relative Error for every Time layer, S R E T k = j = 1 M U ˜ j k U j k U j k . The purpose is to test the stability of HASC-N difference scheme.
The SRET values of the two schemes are shown in Figure 5. When α = 0.7 ,   N = 100 ,   M = 71 , the SRET values of the two schemes tend to 0 with the increase of the time grid numbers. Therefore, the C-N scheme and HASC-N difference scheme of inhomogeneous TFFE (38) are stable, and the results in Figure 5 verify the correctness of Theorem 2.
The spatial convergence order and the temporal convergence order of HASC-N difference scheme are compared. The error E m , τ , the error E h , n , the spatial convergence order O r d e r 1 and temporal convergence order O r d e r 2 are defined as follows [42,43]:
E m , τ = max 0 k N U ˜ m k U m k , 0 m M , E h , n = max 0 j M U ˜ j n U j n , 0 n N .
O r d e r 1 = l n ( E h 1 , n E h 1 , n E h 2 , n E h 2 , n ) l n ( h 1 h 1 h 2 h 2 ) , O r d e r 2 = l n ( E ( m , τ 1 ) E ( m , τ 1 ) E ( m , τ 2 ) E ( m , τ 2 ) ) l n ( τ 1 τ 1 τ 2 τ 2 ) .
To verify the spatial convergence order of HASC-N difference scheme, take M = 21 , 41 , 81 , 161 and τ = h 2 / 4 .  Table 2 shows that the spatial convergence order of HASC-N difference scheme is O ( h 2 ) , and its error decreases gradually with the increase of space step. The theoretical analysis is validated by numerical test data.
Calculate the temporal convergence order of HASC-N difference scheme. Fixed space step h = 1/101, namely, M = 101 and let N = 16, 32, 64, 128. As can be seen from Table 3, the temporal convergence order of HASC-N scheme reaches O ( τ 2 α ) , and the error of HASC-N difference scheme decreases gradually with the increase of time step.
The numerical experimental data in Table 2 and Table 3 can correspond to the conclusion of spatial convergence order O ( h 2 ) and temporal convergence order O ( τ 2 α ) in Theorem 3.
Speed-up ratio S p = T / T p (T is the CPU time in C-N scheme, T p is the CPU time in HASC-N difference scheme) and efficiency E p = S p / P (p is the number of processor cores) [17]. Take α = 0.7 ,   N = 100 , space grid points M = 201 , 401, 601, 801, 1001, 1201. Table 4 shows the CPU time of C-N scheme solution and HASC-N scheme solution, speed-up ratio ( S p ) and efficiency ( E p ) of HASC-N scheme solution.
According to the comparative analysis in Table 4, the computational efficiency of the HASC-N difference scheme of inhomogeneous TFFE (38) is obviously better than that of the C-N scheme. With the encryption of the spatial grid, the computational time advantage of the HASC-N difference scheme is more and more prominent than that of the C-N scheme. The speed-up ratio of the HASC-N difference scheme and C-N scheme is above 3, and the efficiency is about 1.5 . The results show that the HASC-N difference scheme has obvious parallel computing characteristics.
Example 2.
Consider the inhomogeneous TFFE for the discontinuity of the time fractional derivative at the initial time:
α u ( x , t ) t α = 2 u ( x , t ) x 2 + u ( x , t ) ( 1 u ( x , t ) ) + g 2 x , t , ( x , t ) 0 , 1 × 0 , 1 , u ( x , 0 ) = 0 , x 0 , 1 , u ( 0 , t ) = u ( 1 , t ) = 0 , t 0 , 1 .
where g 2 x , t = Γ ( 1 + α ) x x 1 2 t α x x 1 t α ( 1 x x 1 t α ) ,   0 < α 1 . Exact solution of the inhomogeneous TFFE (39): u ( x , t ) = x x 1 t α .
When α = 0.5 ,   N = 100 ,   M = 71 , the exact solution surface and HASC-N scheme solution surface are as follows:
When α = 0.5 , the solution of inhomogeneous TFFE (39) has initial singularity near t = 0 , and the solution is smooth away from t = 0 , as shown in Figure 6 and Figure 7.
In Table 5, the number of space grids M = 21 , 41 , 81 , 161 and let τ = h 2 / 4 . In Table 6, the number of time grids N = 8 , 16 , 32 , 64 , and the fixed spatial step h = 1 / 101 , namely, M = 101 . When α is set to different values, it can be seen from Table 5 that the spatial convergence order of HASC-N difference scheme is O ( h 2 ) , and the temporal convergence order given in Table 6 can reach O ( τ α ) . Therefore, the fractional derivative α u ( x , t ) t α of inhomogeneous TFFE (39) is a discontinuous function (satisfying the weak regularity conditions), the partial derivative of spatial direction satisfies δ u x δ C 0 , L × 0 , T ( δ 0 , 1 , 2 , 3 , 4 ) , and the solution of HASC-N difference scheme converges to O ( τ α + h 2 ) , which verifies the first statement in Remark 2.
Example 3.
Consider the inhomogeneous TFFE for the discontinuity of the temporal fractional derivative at the initial time and the discontinuity of the spatial derivative:
α u ( x , t ) t α = 2 u ( x , t ) x 2 + u ( x , t ) ( 1 u ( x , t ) ) + g 3 x , t , ( x , t ) 0 , 1 × 0 , 1 , u ( x , 0 ) = 0 , x 0 , 1 , u ( 0 , t ) = u ( 1 , t ) = 0 , t 0 , 1 .
where 0 < α 1 ,
g 3 x , t = Γ ( 2 ) Γ ( 2 α ) sin ( π x ) + Γ ( 3 ) Γ ( 3 α ) [ x 1 x ] 0.5 t 2 α + π 2 t α sin ( π x ) t 2 2 [ 1 2 [ x 1 x ] 1.5 ( 1 2 x ) 2 2 [ x 1 x ] 0.5 ] ( t α sin ( π x ) + t 2 [ x 1 x ] 0.5 ) ( 1 ( t α sin ( π x ) + t 2 [ x 1 x ] 0.5 ) ) .
Exact solution of the inhomogeneous TFFE (40):
u ( x , t ) = t α sin ( π x ) + t 2 [ x 1 x ] 0.5 .
When α is set to different values, N = 100 ,   M = 71 , the HASC-N scheme solution is compared with the exact solution at t = 0.5 , and the calculation results are shown in Table 7:
According to the inhomogeneous TFFE (40) and its exact solution u ( x , t ) (41), the equation has the initial singularity. In addition, it has singularity near the boundary of u ( 0 , t ) and u ( 1 , t ) . The Equation (40) meets the weak regularity condition and produces a certain disturbance to the HASC-N difference scheme. According to the analysis of Table 7, although there is some error between the HASC-N scheme solution and the exact solution, the approximation effect is still satisfactory.
The solution (41) of the inhomogeneous TFFE (40) has an initial singularity (satisfying the weak regularity condition), and the partial derivative of the spatial direction δ u x δ ( δ 0 , 1 , 2 , 3 , 4 ) is a discontinuous function on 0 , L × 0 , T . In order to explore whether the truncation error of HASC-N scheme solution for inhomogeneous TFFE (40) has a clear global control, The loose L 2 norm is used as a measure, and the L 2 norm is defined as follows [44]:
E 2 m , τ = k = 1 N U ˜ m k U m k 2 τ 1 2 , 0 m M ,
E 2 h , n = j = 1 M U ˜ j n U j n 2 h 1 2 , 0 n N .
So the spatial convergence order O r d e r 3 and temporal convergence order O r d e r 4 are defined as:
O r d e r 3 = l n ( E 2 h 1 , n E 2 h 1 , n E 2 h 2 , n E 2 h 2 , n ) l n ( h 1 h 1 h 2 h 2 ) , O r d e r 4 = l n ( E 2 ( m , τ 1 ) E 2 ( m , τ 1 ) E 2 ( m , τ 2 ) E 2 ( m , τ 2 ) ) l n ( τ 1 τ 1 τ 2 τ 2 ) .
As shown in Table 8 and Table 9, the local truncation error of the HASC-N difference scheme lacks a clear overall control. Even if the loose discrete L 2 norm is used as the measure, the local truncation error has no definite order. Therefore, the analysis path of convergence order based on strong regularity condition (28) is no longer effective. This verifies the second statement in Remark 2.

7. Conclusions

Most schemes with parallelism are not unconditionally stable for a long time, or the stability meets the requirements but the space has only precision O ( h ) [45,46]. The HASC-N difference scheme for inhomogeneous TFFE is constructed in this paper, which is unconditionally stable. The convergence order of HASC-N difference scheme is O ( τ 2 α + h 2 ) under the strong regularity condition, and O ( τ α + h 2 ) under the weak regularity condition that the time-fractional derivative is discontinuous at the initial time and the space derivative is continuous. Under the weak regularity condition that the time-fractional derivative is discontinuous at the initial time and the spatial derivative is discontinuous, the error of the HASC-N difference scheme lacks a clear global control and does not specify the convergence order. Therefore, the analysis path of convergence order based on strong regularity conditions is no longer effective.
The HASC-N difference scheme has obvious parallel computing properties. The localization characteristics of the HASC-N difference scheme in computing and communication will become more and more remarkable with the continuous encryption of space grid points, which is suitable for parallelized computing systems with distributed storage. The numerical tests verify the theoretical analysis and show that the HASC-N difference scheme in this paper is high-efficient in solving inhomogeneous TFFE.

Author Contributions

Conceptualization, R.L., X.Y. and P.L.; methodology, R.L., X.Y. and P.L.; software, R.L., X.Y. and P.L.; validation, R.L., X.Y. and P.L.; formal analysis, R.L., X.Y. and P.L.; investigation, R.L., X.Y. and P.L.; resources, R.L., X.Y. and P.L.; data curation, R.L., X.Y. and P.L.; writing—original draft preparation, R.L., X.Y. and P.L.; writing—review and editing, R.L., X.Y. and P.L.; visualization, R.L., X.Y. and P.L.; supervision, X.Y.; project administration, X.Y.; funding acquisition, X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No.11371135).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this research. Data sharing does not apply to this research.

Acknowledgments

We would like to thank Lifei Wu of North China Electric Power University for many helpful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Uchaikin, V.V. Fractional Derivatives for Physicists and Engineers, Volume II: Applications; Springer: Berlin, Germany, 2013. [Google Scholar]
  2. Rashid, S.; Hammouch, Z.; Aydi, H.; Ahmad, A.G.; Alsharif, A.M. Novel computations of the time-fractional Fisher’s model via generalized fractional integral operators by means of the Elzaki transform. Fractal Fract. 2021, 5, 94. [Google Scholar] [CrossRef]
  3. Li, C.P.; Zeng, F.H. Numerical Methods for Fractional Calculus; CRC Press: New York, NY, USA, 2015. [Google Scholar]
  4. Sabatier, J.; Agrawal, O.P. Advances in Fractional Calculus: Theoretical Developments and Applications in Physics and Engineering; Machado, J.A.T., Ed.; Beijing World Publishing Corporation: Beijing, China, 2014. [Google Scholar]
  5. Diethelm, K. The Analysis of Fraction Differential Equations; Springer: Berlin, Germany, 2010. [Google Scholar]
  6. Liu, F.W.; Zhuang, P.H.; Liu, Q.X. Numerical Methods and its Application of Fractional Partial Differential Equation; Science Press: Beijing, China, 2015. (In Chinese) [Google Scholar]
  7. Sun, Z.Z.; Gao, G.H. Finite Difference Methods for Fractional Differential Equations, 2nd ed.; Science Press: Beijing, China, 2021. (In Chinese) [Google Scholar]
  8. Guo, B.L.; Pu, X.K.; Huang, F.H. Fractional Partial Differential Equations and Their Numerical Solutions; Science Press: Beijing, China, 2015. [Google Scholar]
  9. Zhang, X.D.; He, Y.N.; Wei, L.L.; Tang, B.; Wang, S.L. A fully discrete local discontinuous Galerkin method for one-dimensional time-fractional Fisher’s equation. Int. J. Comput. Math. 2014, 91, 2021–2038. [Google Scholar] [CrossRef]
  10. Alquran, M.; Al-Khaled, K.; Sardar, T.; Chattopadhyay, J. Revisited Fisher’s equation in a new outlook: A fractional derivative approach. Phys. A 2015, 438, 81–93. [Google Scholar] [CrossRef]
  11. Mejía, C.E.; Piedrahita, A. A numerical method for a time-fractional advection–dispersion equation with a nonlinear source term. J. Appl. Math. Comput. 2019, 61, 593–609. [Google Scholar] [CrossRef]
  12. Al Qurashi, M.M.; Korpinar, Z.; Baleanu, D.; Mustafa Inc. A new iterative algorithm on the time-fractional Fisher equation: Residual power series method. Adv. Mech. Eng. 2017, 9, 1–8. [Google Scholar] [CrossRef] [Green Version]
  13. Khader, M.M.; Saad, K.M. A numerical approach for solving the fractional Fisher equation using Chebyshev spectral collocation method. Chaos Solitons Fractals 2018, 110, 169–177. [Google Scholar] [CrossRef]
  14. Li, D.F.; Liao, H.L.; Sun, W.W.; Wang, J.L.; Zhang, J.W. Analysis of L1-Galerkin FEMs for time-fractional nonlinear parabolic problems. Commun. Comput. Phys. 2018, 24, 86–103. [Google Scholar] [CrossRef] [Green Version]
  15. Veeresha, P.; Prakasha, D.G.; Baskonus, H.M. Novel simulations to the time-fractional Fisher’s equation. Math. Sci. 2019, 13, 33–42. [Google Scholar] [CrossRef]
  16. Petter, B.; Mitchell, L. Parallel Solution of Partial Differential Equations; Springer: New York, NY, USA, 2000. [Google Scholar]
  17. Pacheco, P. An Introduction to Parallel Programming; Morgan Kaufmann: Burlington, VT, USA, 2011. [Google Scholar]
  18. Gong, C.Y.; Bao, W.M.; Tang, G.J. A parallel algorithm for the Riesz fractional reaction-diffusion equation with explicit finite difference method. Fract. Calc. Appl. Anal. 2013, 16, 654–669. [Google Scholar] [CrossRef]
  19. Sweilam, N.H.; Moharram, H.; Moniem, N.K.A.; Ahmed, S. A parallel Crank–Nicolson finite difference method for time-fractional parabolic equation. J. Numer. Math. 2014, 22, 363–382. [Google Scholar] [CrossRef]
  20. Biala, T.A.; Khaliq, A.Q.M. Parallel algorithms for nonlinear time-space fractional parabolic PDEs. J. Comput. Phys. 2018, 375, 135–154. [Google Scholar] [CrossRef]
  21. Fu, H.F.; Wang, H. A preconditioned fast parareal finite difference method for space-time fractional partial differential equation. J. Sci. Comput. 2019, 78, 1724–1743. [Google Scholar] [CrossRef]
  22. Yue, X.Q.; Shu, S.; Xu, X.W.; Bu, W.P.; Pan, K.J. Parallel-in-time multigrid for space-time finite element approximations of two-dimensional space-fractional diffusion equations. Comput. Math. Appl. 2019, 78, 3471–3484. [Google Scholar] [CrossRef]
  23. Liu, H.; Cheng, A.J.; Wang, H. A parareal finite volume method for variable-order time-fractional diffusion equations. J. Sci. Comput. 2020, 85, 19. [Google Scholar] [CrossRef]
  24. Lorin, E. A parallel algorithm for space-time-fractional partial differential equations. Adv. Differ. Equ. 2020, 2020, 283. [Google Scholar] [CrossRef]
  25. Wang, Q.L.; Liu, J.; Gong, C.Y.; Tang, X.T.; Fu, G.T.; Xing, Z.C. An efficient parallel algorithm for Caputo fractional reaction-diffusion equation with implicit finite-difference method. Adv. Differ. Equ. 2016, 2016, 207. [Google Scholar] [CrossRef] [Green Version]
  26. Yang, X.; Wu, L. A New Kind of Parallel Natural Difference Method for Multi-Term Time Fractional Diffusion Model. Mathematics 2020, 8, 596. [Google Scholar] [CrossRef] [Green Version]
  27. Henry, B.I.; Langlands, T.A.M.; Wearne, S.L. Anomalous diffusion with linear reaction dynamics: From continuous time random walks to fractional reaction-diffusion equations. Phys. Rev. E 2006, 74, 031116. [Google Scholar] [CrossRef] [Green Version]
  28. Angstmann, C.N.; Henry, B.I. Time Fractional Fisher–KPP and Fitzhugh–Nagumo Equations. Entropy 2020, 22, 1035. [Google Scholar] [CrossRef]
  29. Sandev, T.; Tomovski, Z. Fractional Equations and Models: Theory and Applications; Springer: Berlin, Germany, 2019. [Google Scholar]
  30. Ngoc, T.B.; Tri, V.V.; Hammouch, Z.; Can, N.H. Stability of a class of problems for time-space fractional pseudo-parabolic equation with datum measured at terminal time. Appl. Numer. Math. 2021, 167, 308–329. [Google Scholar] [CrossRef]
  31. Dang, D.T.; Nane, E.; Nguyen, D.M.; Tuan, N.H. Continuity of Solutions of a Class of Fractional Equations. Potential Anal 2018, 49, 423–478. [Google Scholar] [CrossRef] [Green Version]
  32. Zhang, Y.N.; Sun, Z.Z. Alternating direction implicit schemes for the two-dimensional fractional sub-diffusion equation. J. Comput. Phys. 2011, 230, 8713–8728. [Google Scholar] [CrossRef]
  33. Liu, N.; Liu, Y.; Li, H.; Wang, J. Time second-order finite difference/finite element algorithm for nonlinear time-fractional diffusion problem with fourth-order derivative term. Comput. Math. Appl. 2018, 75, 3521–3536. [Google Scholar] [CrossRef]
  34. Wang, Y.; Liu, Y.; Li, H.; Wang, J. Finite element method combined with second-order time discrete scheme for nonlinear fractional cable equation. Eur. Phys. J. Plus 2016, 131, 61. [Google Scholar] [CrossRef]
  35. Zhang, B.L.; Gu, T.X.; Mo, Z.Y. Principles and Methods of Numerical Parallel Computation; National Defence Industry Press: Beijing, China, 1999. (In Chinese) [Google Scholar]
  36. Zhou, Y.L. Difference schemes with intrinsic parallelism for quasi-linear parabolic systems. Sci. China Math. 1997, 40, 270–278. [Google Scholar] [CrossRef]
  37. Stynes, M.; O’Riordan, E.; Gracia, J.L. Error analysis of a finite difference method on graded meshes for a time-fractional diffusion equation. SIAM J. Numer. Anal. 2017, 55, 1057–1079. [Google Scholar] [CrossRef]
  38. Shen, J.Y.; Sun, Z.Z.; Du, R. Fast finite difference schemes for time-fractional diffusion equations with a weak singularity at initial time. East Asian J. Appl. Math. 2018, 8, 834–858. [Google Scholar] [CrossRef]
  39. Yang, X.Z.; Liu, X.L. Numerical analysis of two new finite difference methods for time-fractional telegraph equation. Discrete Contin. Dyn. Syst. Ser. B 2021, 26, 3921–3942. [Google Scholar] [CrossRef]
  40. Evans, G.; Blackledge, J.; Yardley, P. Numerical Methods for Partial Differential Equations; Springer: Berlin, Germany, 2012. [Google Scholar]
  41. Kumar, D.; Chaudhary, S.; Srinivas Kumar, V.V.K. Fractional Crank–Nicolson–Galerkin finite element scheme for the time-fractional nonlinear diffusion equation. Numer. Methods Partial Differ. Equ. 2019, 35, 2056–2075. [Google Scholar] [CrossRef]
  42. Vong, S.; Lyu, P.; Wang, Z.B. A compact difference scheme for fractional sub-diffusion equations with the spatially variable coefficient under neumann boundary conditions. J. Sci. Comput. 2016, 66, 725–739. [Google Scholar] [CrossRef]
  43. Tverdyi, D.; Parovik, R. Investigation of Finite-Difference Schemes for the Numerical Solution of a Fractional Nonlinear Equation. Fractal Fract. 2022, 6, 23. [Google Scholar] [CrossRef]
  44. Long, L.D.; Luc, N.H.; Tatar, S.; Baleanu, D.; Can, N.H. An inverse source problem for pseudo-parabolic equation with Caputo derivative. J. Appl. Math. Comput. 2022, 68, 739–765. [Google Scholar] [CrossRef]
  45. Yuan, G.W.; Sheng, Z.Q.; Hang, X. The unconditional stability of parallel difference schemes with second order convergence for nonlinear parabolic system. J. Partial Differ. Equ. 2007, 20, 45–64. [Google Scholar]
  46. Deng, W.H.; Zhang, Z.J. High Accuracy Algorithms for the Differential Equation Governing Anomalous Diffusion, Algorithm and Models for Anomalous Diffusion; World Scientific: Singapore, 2019. [Google Scholar]
Figure 1. Construction principle of HASC-N difference scheme.
Figure 1. Construction principle of HASC-N difference scheme.
Fractalfract 06 00259 g001
Figure 2. Exact solution surface for Example 1.
Figure 2. Exact solution surface for Example 1.
Fractalfract 06 00259 g002
Figure 3. C-N scheme solution surface for Example 1. ( α = 0.7 ).
Figure 3. C-N scheme solution surface for Example 1. ( α = 0.7 ).
Fractalfract 06 00259 g003
Figure 4. HASC-N scheme solution surface for Example 1. ( α = 0.7 ).
Figure 4. HASC-N scheme solution surface for Example 1. ( α = 0.7 ).
Fractalfract 06 00259 g004
Figure 5. SRET values in two schemes for Example 1. ( α = 0.7 ).
Figure 5. SRET values in two schemes for Example 1. ( α = 0.7 ).
Fractalfract 06 00259 g005
Figure 6. Exact solution surface for Example 2. ( α = 0.5 ).
Figure 6. Exact solution surface for Example 2. ( α = 0.5 ).
Fractalfract 06 00259 g006
Figure 7. HASC-N scheme solution surface for Example 2. ( α = 0.5 ).
Figure 7. HASC-N scheme solution surface for Example 2. ( α = 0.5 ).
Fractalfract 06 00259 g007
Table 1. Comparison between exact solution and HASC-N scheme solution for Example 1 ( t = 0.5 ).
Table 1. Comparison between exact solution and HASC-N scheme solution for Example 1 ( t = 0.5 ).
α x
0.250.50.75
0.3Exact solution5.752108 × 10 2 2.549964 × 10 3 −5.763390 × 10 2
HASC-N scheme solution5.517192 × 10 2 2.529212 × 10 3 −5.600727 × 10 2
0.5Exact solution5.752108 × 10 2 2.549964 × 10 3 −5.763390 × 10 2
HASC-N scheme solution5.565180 × 10 2 2.490347 × 10 3 −5.608690 × 10 2
0.7Exact solution5.752108 × 10 2 2.549964 × 10 3 −5.763390 × 10 2
HASC-N scheme solution5.605673 × 10 2 2.486836 × 10 3 −5.629021 × 10 2
Table 2. Numerical error and spatial convergence order of HASC-N difference scheme for Example 1.
Table 2. Numerical error and spatial convergence order of HASC-N difference scheme for Example 1.
α MN E h , n Order 1
0.5211005.805049 × 10 3
414001.535960 × 10 3 1.987252
8116003.930984 × 10 4 2.001614
16164009.953251 × 10 5 1.999513
0.7211005.272847 × 10 3
414001.326208 × 10 3 2.062997
8116003.261035 × 10 4 2.060378
16164007.997470 × 10 5 2.045991
0.9211004.857505 × 10 3
414001.220010 × 10 3 2.065118
8116003.029576 × 10 4 2.045922
16164007.518864 × 10 5 2.028652
Table 3. Numerical error and temporal convergence order of HASC-N difference scheme for Example 1.
Table 3. Numerical error and temporal convergence order of HASC-N difference scheme for Example 1.
α MN E m , τ Order 2
0.5101161.035923 × 10 2
323.687553 × 10 3 1.490181
641.303842 × 10 3 1.499895
1284.594546 × 10 4 1.504774
0.7101169.512878 × 10 3
323.879871 × 10 3 1.293873
641.578013 × 10 3 1.297899
1286.392140 × 10 4 1.303738
0.9101168.555828 × 10 3
323.996333 × 10 3 1.098231
641.865814 × 10 3 1.098872
1288.681042 × 10 4 1.103865
Table 4. CPU time, speed-up ratio and efficiency of the schemes for Example 1. ( α = 0.7 ,   N = 100 ).
Table 4. CPU time, speed-up ratio and efficiency of the schemes for Example 1. ( α = 0.7 ,   N = 100 ).
M20140160180110011201
T ( s ) 0.1232510.3251160.4667550.8794791.3573761.891728
T p ( s ) 0.0413640.1067370.1515070.2821100.4217850.582479
S p 2.9796683.0459543.0807493.1175043.2181703.247719
E p 1.4898341.5229771.5403741.5587521.6090851.623860
Table 5. Numerical error and spatial convergence of HASC-N difference scheme for Example 2.
Table 5. Numerical error and spatial convergence of HASC-N difference scheme for Example 2.
α MN E h , n Order 1
0.5211001.820529 × 10 3
414004.481606 × 10 4 2.095107
8116001.114969 × 10 4 2.043181
16164002.783025 × 10 5 2.020325
0.7211002.197871 × 10 3
414005.454330 × 10 4 2.083051
8116001.359905 × 10 4 2.040008
16164003.395738 × 10 5 2.019752
0.9211002.440139 × 10 3
414006.093245 × 10 4 2.073776
8116001.522336 × 10 4 2.036982
16164003.804131 × 10 5 2.018681
Table 6. Numerical error and temporal convergence of HASC-N difference scheme for Example 2.
Table 6. Numerical error and temporal convergence of HASC-N difference scheme for Example 2.
α MN E m , τ Order 2
0.5101161.939613 × 10 2
321.355892 × 10 2 0.516527
649.510714 × 10 3 0.511617
1286.674831 × 10 3 0.510822
0.7101161.137173 × 10 2
326.929028 × 10 3 0.714727
644.230193 × 10 3 0.711930
1282.590579 × 10 3 0.707449
0.9101161.424087 × 10 2
327.566100 × 10 3 0.912416
644.024197 × 10 3 0.910849
1282.154974 × 10 3 0.901031
Table 7. Comparison between exact solution and HASC-N scheme solution for Example 3 ( t = 0.5 ).
Table 7. Comparison between exact solution and HASC-N scheme solution for Example 3 ( t = 0.5 ).
α x
0.250.50.75
0.3Exact solution6.540829 × 10 1 9.271849 × 10 1 6.816085 × 10 1
HASC-N scheme solution6.538316 × 10 1 9.323716 × 10 1 6.823412 × 10 1
0.5Exact solution5.807393 × 10 1 8.198668 × 10 1 6.048696 × 10 1
HASC-N scheme solution5.844902 × 10 1 8.303472 × 10 1 6.094448 × 10 1
0.7Exact solution5.171474 × 10 1 7.268177 × 10 1 5.383339 × 10 1
HASC-N scheme solution5.199174 × 10 1 7.356311 × 10 1 5.417396 × 10 1
Table 8. Numerical error and spatial convergence order of HASC-N difference scheme for Example 3.
Table 8. Numerical error and spatial convergence order of HASC-N difference scheme for Example 3.
α MN E 2 h , n Order 3
0.1211004.456840 × 10 2
414003.011746 × 10 2 0.565421
8116002.088088 × 10 2 0.528418
16164001.448888 × 10 2 0.527236
0.2211002.061921 × 10 2
414008.043848 × 10 2 1.358031
8116003.662202 × 10 2 1.135174
16164001.709010 × 10 3 1.099551
0.3211001.064295 × 10 2
414005.183752 × 10 3 1.037830
8116002.538664 × 10 3 1.029927
16164001.330512 × 10 3 0.932088
0.4211001.206877 × 10 2
414005.867005 × 10 3 1.040582
8116002.594804 × 10 3 1.176999
16164001.241889 × 10 3 1.063090
0.5211001.505094 × 10 2
414006.260401 × 10 3 1.265527
8116002.542224 × 10 3 1.300164
16164001.017199 × 10 3 1.321489
0.6211001.108872 × 10 2
414006.254087 × 10 3 0.826221
8116003.460837 × 10 3 0.853678
16164002.050544 × 10 3 0.755114
0.7211001.210763 × 10 2
414005.822212 × 10 3 1.056277
8116002.310618 × 10 3 1.333288
16164001.009494 × 10 3 1.194646
0.8211001.428816 × 10 2
414005.795162 × 10 3 1.301900
8116002.306303 × 10 3 1.329267
16164001.063531 × 10 3 1.116719
0.9211001.789562 × 10 2
414007.830313 × 10 3 1.192465
8116002.261498 × 10 3 1.791791
16164009.787323 × 10 4 1.208293
Table 9. Numerical error and temporal convergence order of HASC-N difference scheme for Example 3.
Table 9. Numerical error and temporal convergence order of HASC-N difference scheme for Example 3.
α MN E 2 m , τ Order 4
0.1101163.235754 × 10 1
321.322441 × 10 1 1.290898
643.805932 × 10 2 1.796882
1281.137961 × 10 2 1.741799
0.2101163.400465 × 10 1
321.126743 × 10 2 1.593573
643.902410 × 10 2 1.529722
1281.266599 × 10 2 1.623405
0.3101161.441621 × 10 1
324.722499 × 10 2 1.610070
641.306512 × 10 2 1.853830
1284.114175 × 10 3 1.667045
0.4101164.485998 × 10 2
322.509118 × 10 2 0.838249
641.128387 × 10 2 1.152918
1283.700395 × 10 3 1.608511
0.5101168.566342 × 10 2
324.305125 × 10 2 0.992624
641.101626 × 10 2 1.966421
1283.663051 × 10 3 1.588516
0.6101161.175347 × 10 1
324.492863 × 10 2 1.387380
641.211080 × 10 2 1.891341
1283.668064 × 10 3 1.723204
0.7101161.384476 × 10 1
324.957889 × 10 2 1.481543
641.458671 × 10 2 1.765071
1284.035364 × 10 3 1.853884
0.8101161.572810 × 10 1
325.660295 × 10 2 1.474396
641.851979 × 10 2 1.611810
1285.576710 × 10 3 1.731581
0.9101161.769765 × 10 1
326.605010 × 10 2 1.421925
642.405686 × 10 2 1.457112
1288.758845 × 10 3 1.457636
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, R.; Yang, X.; Lyu, P. A New Parallelized Computation Method of HASC-N Difference Scheme for Inhomogeneous Time Fractional Fisher Equation. Fractal Fract. 2022, 6, 259. https://doi.org/10.3390/fractalfract6050259

AMA Style

Liu R, Yang X, Lyu P. A New Parallelized Computation Method of HASC-N Difference Scheme for Inhomogeneous Time Fractional Fisher Equation. Fractal and Fractional. 2022; 6(5):259. https://doi.org/10.3390/fractalfract6050259

Chicago/Turabian Style

Liu, Ren, Xiaozhong Yang, and Peng Lyu. 2022. "A New Parallelized Computation Method of HASC-N Difference Scheme for Inhomogeneous Time Fractional Fisher Equation" Fractal and Fractional 6, no. 5: 259. https://doi.org/10.3390/fractalfract6050259

APA Style

Liu, R., Yang, X., & Lyu, P. (2022). A New Parallelized Computation Method of HASC-N Difference Scheme for Inhomogeneous Time Fractional Fisher Equation. Fractal and Fractional, 6(5), 259. https://doi.org/10.3390/fractalfract6050259

Article Metrics

Back to TopTop