×

Approximate inversion method for time-fractional subdiffusion equations. (English) Zbl 1499.65108

Summary: The finite-difference method applied to the time-fractional subdiffusion equation usually leads to a large-scale linear system with a block lower triangular Toeplitz coefficient matrix. The approximate inversion method is employed to solve this system. A sufficient condition is proved to guarantee the high accuracy of the approximate inversion method for solving the block lower triangular Toeplitz systems, which are easy to verify in practice and have a wide range of applications. The applications of this sufficient condition to several existing finite-difference schemes are investigated. Numerical experiments are presented to verify the validity of theoretical results.

MSC:

65F08 Preconditioners for iterative methods
65F05 Direct numerical methods for linear systems and matrix inversion
35R11 Fractional partial differential equations

Software:

Matlab
Full Text: DOI

References:

[1] ZhaoX, XuQ. Efficient numerical schemes for fractional sub‐diffusion equation with the spatially variable coefficient. Appl Math Model. 2014;38:3848-3859. · Zbl 1429.65210
[2] PodlubnyI. Fractional Differential Equations. New York: Academic Press; 1999. · Zbl 0924.34008
[3] LanglandsT, HenryB. The accuracy and stability of an implicit solution method for the fractional diffusion equation. J Comput Phys. 2005;205:719-736. · Zbl 1072.65123
[4] ChenC, LiuF, TurnerI, AnhV. A fourier method for the fractional diffusion equation describing sub‐diffusion. J Comput Phys. 2007;227:886-897. · Zbl 1165.65053
[5] ZhuangP, LiuF, AnhV, TurnerI. New solution and analytical techniques of the implicit numerical method for the anomalous subdiffusion equation. SIAM J Numer Anal. 2008;46:1079-1095. · Zbl 1173.26006
[6] CuiM. Compact finite difference method for the fractional diffusion equation. J Comput Phys. 2009;228:7792-7804. · Zbl 1179.65107
[7] GaoG, SunZ. A compact finite difference scheme for the fractional sub‐diffusion equations. J Comput Phys. 2011;230:586-595. · Zbl 1211.65112
[8] JiangY, MaJ. High‐order finite element methods for time‐fractional partial differential equations. J Comput Appl Math. 2011;235:3285-3290. · Zbl 1216.65130
[9] LinY, XuC. Finite difference/spectral approximations for the time‐fractional diffusion equation. J Comput Phys. 2007;225:1533-1552. · Zbl 1126.65121
[10] MurioD. Implicit finite difference approximation for time fractional diffusion equations. Comput Math Appl. 2008;56:1138-1145. · Zbl 1155.65372
[11] MustaphaK. An implicit finite‐difference time‐stepping method for a sub‐diffusion equation, with spatial discretization by finite elements. IMA J Numer Anal. 2011;31:719-739. · Zbl 1219.65091
[12] YusteS, AcedoL. An explicit finite difference method and a new von Neumann‐type stability analysis for fractional diffusion equations. SIAM J Numer Anal. 2005;42:1862-1874. · Zbl 1119.65379
[13] ZhangY, SunZ, WuH. Error estimates of Crank‐Nicolson‐type difference schemes for the sub‐diffusion equation. SIAM J Numer Anal. 2011;49:2302-2322. · Zbl 1251.65132
[14] ZengF, LiC, LiuF, TurnerI. The use of finite difference/element approaches for solving the time‐fractional subdiffusion equation. SIAM J Sci Comput. 2013;35:2976-3000. · Zbl 1292.65096
[15] DiethelmK. The analysis of Fractional Differential Equations. Berlin: Springer; 2010. · Zbl 1215.34001
[16] KeR, NgMK, SunH. A fast direct method for block triangular Toeplitz‐like with tri‐diagonal block systems from time‐fractional partial differential equations. J Comput Phys. 2015;303:203-211. · Zbl 1349.65404
[17] LuX, PangH, SunH. Fast approximate inversion of a block triangular Toeplitz matrix with applications to fractional sub‐diffusion equations. Numer Linear Algebra Appl. 2015;22:866-882. · Zbl 1349.65104
[18] GolubG, LoanC. Matrix Computations. 4th ed. Baltimore: Johns Hopkins University Press; 2013. · Zbl 1268.65037
[19] BiniD, LatoucheG, MeiniB. Numerical Methods for Structured Markov Chains. New York: Oxford University Press Inc; 2005. · Zbl 1076.60002
[20] McLeanW. Fast summation by interval clustering for an evolution equation with memory. SIAM J Sci Comput. 2012;34:A3039-A3056. · Zbl 1276.65059
[21] HackbushchW, NowakZP. On the fast matrix multiplication in the boundary element method by panel clustering. Numer Math. 1989;54:463-491. · Zbl 0641.65038
[22] HornR, JohnsonC. Topics in Matrix Analysis. Cambridge: Cambridge University Press; 1991. · Zbl 0729.15001
[23] BiniD, DendievelS, LatoucheG, MeiniB. Computing the exponential of large block‐triangular block‐Toeplitz matrices encountered in fluid queues. Linear Algebra Appl. 2016;502:387-419. · Zbl 1386.65129
[24] LinX, LuX, NgMK, SunH. A fast accurate approximation method with multigrid solver for two‐dimensional fractional sub‐diffusion equation. J Comput Phys. 2016;323:204-218. · Zbl 1415.65196
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.