A fast second-order accurate method for a two-sided space-fractional diffusion equation with variable coefficients. (15th March 2017)

• Record Type: Journal Article
• Authors: Feng, L.B.; Zhuang, P.; Liu, F.; Turner, I.; Anh, V.; Li, J.
• Is Part Of: Computers & mathematics with applications. Volume 73:issue 6(2017)
• Page Start: 1155
• Page End: 1171
• Publication Date: 2017-03-15
• Subjects: Finite difference method -- Riemann–Liouville fractional derivative -- Fractional diffusion equation -- Crank–Nicolson scheme -- Variable coefficients -- Fast Bi-CGSTAB algorithm; Electronic data processing -- Periodicals; Mathematics -- Data processing -- Periodicals; 510.28541
• DOI: 10.1016/j.camwa.2016.06.007
• Languages: English
• ISSNs: 0898-1221

Abstract

Abstract: In this paper, we consider a type of fractional diffusion equation (FDE) with variable coefficients on a finite domain. Firstly, we utilize a second-order scheme to approximate the Riemann–Liouville fractional derivative and present the finite difference scheme. Specifically, we discuss the Crank–Nicolson scheme and solve it in matrix form. Secondly, we prove the stability and convergence of the scheme and conclude that the scheme is unconditionally stable and convergent with the second-order accuracy of O ( τ 2 + h 2 ) . Furthermore, we develop a fast accurate iterative method for the Crank–Nicolson scheme, which only requires storage of O ( m ) and computational cost of O ( m log m ) while retaining the same accuracy and approximation property as Gauss elimination, where m = 1 / h is the partition number in space direction. Finally, several numerical examples are given to show the effectiveness of the numerical method, and the results are in excellent agreement with the theoretical analysis.