Loading...
A higher order numerical method for time fractional partial differential equations with nonsmooth data
Xing, Yanyuan Yan, Yubin
Xing, Yanyuan
Yan, Yubin
Advisors
Editors
Other Contributors
Affiliation
EPub Date
Publication Date
2018-01-02
Submitted Date
Collections
Files
Loading...
Main article
Adobe PDF, 363.72 KB
Other Titles
Abstract
Gao et al. (2014) introduced a numerical scheme to approximate the Caputo fractional derivative with the convergence rate $O(k^{3-\alpha}), 0< \alpha <1$ by directly approximating the integer-order derivative with some finite difference quotients in the definition of the Caputo fractional derivative, see also Lv and Xu (2016), where $k$ is the time step size. Under the assumption that the solution of the time fractional partial differential equation is sufficiently smooth, Lv and Xu (2016) proved by using energy method that the corresponding numerical method for solving time fractional partial differential equation has the convergence rate $O(k^{3-\alpha}), 0< \alpha <1$ uniformly with respect to the time variable $t$. However, in general the solution of the time fractional partial differential equation has low regularity and in this case the numerical method fails to have the convergence rate $O(k^{3- \alpha}), 0 < \alpha <1$ uniformly with respect to the time variable $t$. In this paper, we first obtain a similar approximation scheme to the Riemann-Liouville fractional derivative with the convergence rate $O(k^{3- \alpha}), 0 < \alpha <1$ as in Gao \et \cite{gaosunzha} (2014) by approximating the Hadamard finite-part integral with the piecewise quadratic interpolation polynomials. Based on this scheme, we introduce a time discretization scheme to approximate the time fractional partial differential equation and show by using Laplace transform methods that the time discretization scheme has the convergence rate $O(k^{3- \alpha}), 0 < \alpha <1$ for any fixed $t_{n}>0$ for smooth and nonsmooth data in both homogeneous and inhomogeneous cases. Numerical examples are given to show that the theoretical results are consistent with the numerical results.
Citation
Xing, Y., & Yan, Y. (2018). A higher order numerical method for time fractional partial differential equations with nonsmooth data. Journal of Computational Physics, 357, 305-323. https://doi.org/10.1016/j.jcp.2017.12.035
Publisher
Elsevier
Journal
Journal of Computational Physics
Research Unit
DOI
10.1016/j.jcp.2017.12.035
PubMed ID
PubMed Central ID
Type
Article
Language
en
Description
Series/Report no.
ISSN
0021-9991