期刊
JOURNAL OF COMPUTATIONAL PHYSICS
卷 259, 期 -, 页码 33-50出版社
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2013.11.017
关键词
Caputo fractional derivative; Fractional numerical differentiation formula; L1-2 formula; L1 formula; Sub-diffusion
资金
- 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]
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.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据