A Fast Algorithm for Solving the Space-Time Fractional Diffusion Equation
In this paper, we propose a fast algorithm for efficient and accurate solution of the space-time fractional diffusion equations defined in a rectangular domain. The spatial discretization is done by using the central finite difference scheme and matrix transfer technique. Due to its nonlocality, numerical discretization of the spectral fractional Laplacian (−Δ)sα/2 results in a large dense matrix. This causes considerable challenges not only for storing the matrix but also for computing matrix-vector products in practice. By utilizing the compact structure of the discrete system and the discrete sine transform, our algorithm avoids to store the large matrix from discretizing the nonlocal operator and also significantly reduces the computational costs. We then use the Laplace transform method for time integration of the semi-discretized system and a weighted trapezoidal method to numerically compute the convolutions needed in the resulting scheme. Various experiments are presented to demonstrate the efficiency and accuracy of our method.
S. Duo et al., "A Fast Algorithm for Solving the Space-Time Fractional Diffusion Equation," Computers and Mathematics with Applications, vol. 75, no. 6, pp. 1929-1941, Elsevier Ltd, Mar 2018.
The definitive version is available at https://doi.org/10.1016/j.camwa.2017.04.008
Mathematics and Statistics
Center for High Performance Computing Research
Keywords and Phrases
Caputo fractional derivative; Discrete sine transform; Fractional diffusion equation; Matrix transfer method; Matrix-vector product; Spectral fractional Laplacian
International Standard Serial Number (ISSN)
Article - Journal
© 2018 Elsevier Ltd, All rights reserved.
01 Mar 2018