Improved Finite Difference Results for the Caputo Time-Fractional Diffusion Equation
Abstract
We begin with a treatment of the Caputo time-fractional diffusion equation, by using the Laplace transform, to obtain a Volterra intego-differential equation where we may examine the weakly singular nature of this convolution kernel.\iffalse The order of fractional derivative, , is tied to finite difference methods and is of great interest.\fi We examine this new equation and utilize a numerical scheme that is derived in parallel to the L1-method for the time variable and a usual fourth order approximation in the spatial variable. The main method derived in this paper has a rate of convergence of for , which improves previous estimates by a factor of . We also present a novel alternative method for a first order approximation in time, which allows us to relax our regularity assumption to , while exhibiting order of convergence slightly less than in time. This allows for a much wider class of functions to be analyzed which was previously not possible under the L1-method. We present numerical examples demonstrating these results and discuss future improvements and implications by using these techniques.
Cite
@article{arxiv.1811.12910,
title = {Improved Finite Difference Results for the Caputo Time-Fractional Diffusion Equation},
author = {Wesley Davis and Richard Noren and Ke Shi},
journal= {arXiv preprint arXiv:1811.12910},
year = {2020}
}
Comments
22 pages, 3 tables, 1 appendix