Journal
JOURNAL OF COMPUTATIONAL PHYSICS
Volume 259, Issue -, Pages 33-50Publisher
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2013.11.017
Keywords
Caputo fractional derivative; Fractional numerical differentiation formula; L1-2 formula; L1 formula; Sub-diffusion
Funding
- National Natural Science Foundation of China [11271068, 11326225, 61304169]
- Natural Science Youth Foundation of Jiangsu Province, China [BK20130860]
- Scientific Research Foundation of Nanjing University of Posts and Telecommunications, China [NY213051]
Ask authors/readers for more resources
In the present work, first, a new fractional numerical differentiation formula (called the L1-2 formula) to approximate the Caputo fractional derivative of order alpha (0 < alpha < 1) is developed. It is established by means of the quadratic interpolation approximation using three points (t(j-2), f(t(j-2))), (t(j-1), f(t(j-1))) and (t(j), f(t(j))) for the integrand f(t) on each small interval [t(j-1), t(j)] (j >= 2), while the linear interpolation approximation is applied on the first small interval [t(0), t(1)]. As a result, the new formula can be formally viewed as a modification of the classical L1 formula, which is obtained by the piecewise linear approximation for f (t). Both the computational efficiency and numerical accuracy of the new formula are superior to that of the L1 formula. The coefficients and truncation errors of this formula are discussed in detail. Two test examples show the numerical accuracy of L1-2 formula. Second, by the new formula, two improved finite difference schemes with high order accuracy in time for solving the time-fractional sub-diffusion equations on a bounded spatial domain and on an unbounded spatial domain are constructed, respectively. In addition, the application of the new formula into solving fractional ordinary differential equations is also presented. Several numerical examples are computed. The comparison with the corresponding results of finite difference methods by the L1 formula demonstrates that the new L1-2 formula is much more effective and more accurate than the L1 formula when solving time-fractional differential equations numerically. (C) 2013 Elsevier Inc. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available