Show simple item record

dc.contributor.authorLiu, Yanzhi*
dc.contributor.authorRoberts, Jason A.*
dc.contributor.authorYan, Yubin*
dc.date.accessioned2017-08-10T10:47:06Z
dc.date.available2017-08-10T10:47:06Z
dc.date.issued2017-09-21
dc.identifier.citationYanzhi, L., Roberts, J., & Yan, Y. (2018). Detailed error analysis for a fractional Adams method with graded meshes. Numerical Algorithms, 78(4), 1195-1216. https://doi.org/10.1007/s11075-017-0419-5en
dc.identifier.issn1572-9265
dc.identifier.doi10.1007/s11075-017-0419-5
dc.identifier.urihttp://hdl.handle.net/10034/620594
dc.descriptionThe final publication is available at Springer via http://dx.doi.org/10.1007/s11075-017-0419-5
dc.description.abstractWe consider a fractional Adams method for solving the nonlinear fractional differential equation $\, ^{C}_{0}D^{\alpha}_{t} y(t) = f(t, y(t)), \, \alpha >0$, equipped with the initial conditions $y^{(k)} (0) = y_{0}^{(k)}, k=0, 1, \dots, \lceil \alpha \rceil -1$. Here $\alpha$ may be an arbitrary positive number and $ \lceil \alpha \rceil$ denotes the smallest integer no less than $\alpha$ and the differential operator is the Caputo derivative. Under the assumption $\, ^{C}_{0}D^{\alpha}_{t} y \in C^{2}[0, T]$, Diethelm et al. \cite[Theorem 3.2]{dieforfre} introduced a fractional Adams method with the uniform meshes $t_{n}= T (n/N), n=0, 1, 2, \dots, N$ and proved that this method has the optimal convergence order uniformly in $t_{n}$, that is $O(N^{-2})$ if $\alpha > 1$ and $O(N^{-1-\alpha})$ if $\alpha \leq 1$. They also showed that if $\, ^{C}_{0}D^{\alpha}_{t} y(t) \notin C^{2}[0, T]$, the optimal convergence order of this method cannot be obtained with the uniform meshes. However, it is well known that for $y \in C^{m} [0, T]$ for some $m \in \mathbb{N}$ and $ 0 < \alpha 1$, we show that the optimal convergence order of this method can be recovered uniformly in $t_{n}$ even if $\, ^{C}_{0}D^{\alpha}_{t} y$ behaves as $t^{\sigma}, 0< \sigma <1$. Numerical examples are given to show that the numerical results are consistent with the theoretical results.
dc.language.isoenen
dc.publisherSpringeren
dc.relation.urlhttps://link.springer.com/article/10.1007/s11075-017-0419-5en
dc.rights.urihttp://creativecommons.org/licenses/by/4.0/en
dc.subjectFractional differential equationsen
dc.subjectCaputo derivativeen
dc.subjectAdams methoden
dc.titleDetailed error analysis for a fractional adams method with graded meshesen
dc.typeArticleen
dc.contributor.departmentLvliang University; University of Chester;en
dc.identifier.journalNumerical Algorithms
dc.date.accepted2017-07-21
or.grant.openaccessYesen
rioxxterms.funderUnfundeden
rioxxterms.identifier.projectUnfundeden
rioxxterms.versionAMen
rioxxterms.licenseref.startdate2018-09-21
html.description.abstractWe consider a fractional Adams method for solving the nonlinear fractional differential equation $\, ^{C}_{0}D^{\alpha}_{t} y(t) = f(t, y(t)), \, \alpha >0$, equipped with the initial conditions $y^{(k)} (0) = y_{0}^{(k)}, k=0, 1, \dots, \lceil \alpha \rceil -1$. Here $\alpha$ may be an arbitrary positive number and $ \lceil \alpha \rceil$ denotes the smallest integer no less than $\alpha$ and the differential operator is the Caputo derivative. Under the assumption $\, ^{C}_{0}D^{\alpha}_{t} y \in C^{2}[0, T]$, Diethelm et al. \cite[Theorem 3.2]{dieforfre} introduced a fractional Adams method with the uniform meshes $t_{n}= T (n/N), n=0, 1, 2, \dots, N$ and proved that this method has the optimal convergence order uniformly in $t_{n}$, that is $O(N^{-2})$ if $\alpha > 1$ and $O(N^{-1-\alpha})$ if $\alpha \leq 1$. They also showed that if $\, ^{C}_{0}D^{\alpha}_{t} y(t) \notin C^{2}[0, T]$, the optimal convergence order of this method cannot be obtained with the uniform meshes. However, it is well known that for $y \in C^{m} [0, T]$ for some $m \in \mathbb{N}$ and $ 0 < \alpha <m$, the Caputo fractional derivative $\, ^{C}_{0}D^{\alpha}_{t} y(t) $ takes the form \lq \lq $\, ^{C}_{0}D^{\alpha}_{t} y(t) = c t^{\lceil \alpha \rceil -\alpha} + \mbox{\emph{smoother terms}}$\rq\rq \cite[Theorem 2.2]{dieforfre}, which implies that $\, ^{C}_{0}D^{\alpha}_{t} y $ behaves as $t^{\lceil \alpha \rceil -\alpha}$ which is not in $ C^{2}[0, T]$. By using the graded meshes $t_{n}= T (n/N)^{r}, n=0, 1, 2, \dots, N$ with some suitable $r > 1$, we show that the optimal convergence order of this method can be recovered uniformly in $t_{n}$ even if $\, ^{C}_{0}D^{\alpha}_{t} y$ behaves as $t^{\sigma}, 0< \sigma <1$. Numerical examples are given to show that the numerical results are consistent with the theoretical results.
rioxxterms.publicationdate2017-09-21


Files in this item

Thumbnail
Name:
liurobertsyan_2017_06_21.pdf
Size:
383.1Kb
Format:
PDF
Request:
PDF of Article

This item appears in the following Collection(s)

Show simple item record

http://creativecommons.org/licenses/by/4.0/
Except where otherwise noted, this item's license is described as http://creativecommons.org/licenses/by/4.0/