Year: 2024
Author: Boyang Yu, Yonghai Li, Jiangguo Liu
International Journal of Numerical Analysis and Modeling, Vol. 21 (2024), Iss. 3 : pp. 431–458
Abstract
This paper develops novel fast numerical solvers for subdiffusion problems with spatial interfaces. These problems are modeled by partial differential equations that contain both fractional order and conventional first order time derivatives. The former is non-local and approximated by L1 and L2 discretizations along with fast evaluation algorithms based on approximation by sums of exponentials. This results in an effective treatment of the “long-tail” kernel of subdiffusion. The latter is local and hence conventional implicit Euler or Crank-Nicolson discretizations can be used. Finite volumes are utilized for spatial discretization based on consideration of local mass conservation. Interface conditions for mass and fractional fluxes are incorporated into these fast solvers. Computational complexity and implementation procedures are briefly discussed. Numerical experiments demonstrate accuracy and efficiency of these new fast solvers.
You do not have full access to this article.
Already a Subscriber? Sign in as an individual or via your institution
Journal Article Details
Publisher Name: Global Science Press
Language: English
DOI: https://doi.org/10.4208/ijnam2024-1017
International Journal of Numerical Analysis and Modeling, Vol. 21 (2024), Iss. 3 : pp. 431–458
Published online: 2024-01
AMS Subject Headings: Global Science Press
Copyright: COPYRIGHT: © Global Science Press
Pages: 28
Keywords: Caputo and Riemann-Liouville derivatives fast numerical solvers fractional order fluxes interface problems subdiffusion sum of exponentials (SOE).