Show simple item record

dc.contributor.authorXing, Yanyuan*
dc.contributor.authorYan, Yubin*
dc.date.accessioned2018-01-24T11:48:39Z
dc.date.available2018-01-24T11:48:39Z
dc.date.issued2018-01-02
dc.identifier.citationXing, 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.035en
dc.identifier.issn0021-9991
dc.identifier.doi10.1016/j.jcp.2017.12.035
dc.identifier.urihttp://hdl.handle.net/10034/620810
dc.description.abstractGao 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.
dc.language.isoenen
dc.publisherElsevieren
dc.relation.urlhttps://www.sciencedirect.com/science/article/pii/S0021999117309294en
dc.subjectTime fractional partial differential equationen
dc.subjectError estimatesen
dc.subjectLaplace transformen
dc.subjectCaputo fractional derivativeen
dc.titleA higher order numerical method for time fractional partial differential equations with nonsmooth dataen
dc.typeArticleen
dc.contributor.departmentLvliang University; University of Chesteren
dc.identifier.journalJournal of Computational Physics
dc.date.accepted2017-12-24
or.grant.openaccessYesen
rioxxterms.funderUnfundeden
rioxxterms.identifier.projectUnfundeden
rioxxterms.versionAMen
rioxxterms.licenseref.startdate2020-01-02
refterms.dateFCD2019-07-15T09:55:36Z
refterms.versionFCDAM
html.description.abstractGao 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.


Files in this item

Thumbnail
Name:
Publisher version
Thumbnail
Name:
revision_2017_12_01.pdf
Embargo:
2020-01-02
Size:
363.7Kb
Format:
PDF
Description:
Main article

This item appears in the following Collection(s)

Show simple item record