A fast second-order implicit difference method for time-space fractional advection-diffusion equation

DOI10.1080/01630563.2019.1627369arXiv1704.06733WikidataQ127635178 ScholiaQ127635178MaRDI QIDQ6285809FDOQ6285809

Authors: Yong-Liang Zhao, Ting-Zhu Huang, Xian-Ming Gu, Wei-Hua Luo 

Publication date: 21 April 2017

Abstract: In this paper, we consider a fast and second-order implicit difference method for approximation of a class of time-space fractional variable coefficients advection-diffusion equation. To begin with, we construct an implicit difference scheme, based on $L2-1_{sigma}$ formula [A. A. Alikhanov, A new difference scheme for the time fractional diffusion equation, emph{J. Comput. Phys.}, 280 (2015)] for the temporal discretization and weighted and shifted Gr"{u}nwald method for the spatial discretization. Then, unconditional stability of the implicit difference scheme is proved, and we theoretically and numerically show that it converges in the $L_2$-norm with the optimal order $mathcalO(a{u}^2+h^2)$ with time step $au$ and mesh size $h$. Secondly, three fast Krylov subspace solvers with suitable circulant preconditioners are designed to solve the discretized linear systems with the Toeplitz matrix. In each iterative step, these methods reduce the memory requirement of the resulting linear equations from $mathcalO(N^2)$ to $mathcalO(N)$ and the computational complexity from $mathcalO(N^3)$ to $mathcalO(NlogN)$, where $N$ is the number of grid nodes. Finally, numerical experiments are carried out to demonstrate that these methods are more practical than the traditional direct solvers of the implicit difference methods, in terms of memory requirement and computational cost.

Mathematics Subject Classification ID

Numerical analysis (65-XX) Fluid mechanics (76-XX) Mechanics of particles and systems (70-XX)

This page was built for publication: A fast second-order implicit difference method for time-space fractional advection-diffusion equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6285809)