Exponentially Convergent Numerical Method for Abstract Cauchy Problem with Fractional Derivative of Caputo Type
Fuente:
arXiv
Saved in:
| Main Authors: | , |
|---|---|
| Format: | Preprint |
| Published: |
2023
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866912311862624256 |
|---|---|
| author | Sytnyk, Dmytro Wohlmuth, Barbara |
| author_facet | Sytnyk, Dmytro Wohlmuth, Barbara |
| contents | We present an exponentially convergent numerical method to approximate the solution of the Cauchy problem for the inhomogeneous fractional differential equation with an unbounded operator coefficient and Caputo fractional derivative in time. The numerical method is based on the newly obtained solution formula that consolidates the mild solution representations of sub-parabolic, parabolic and sub-hyperbolic equations with sectorial operator coefficient $A$ and non-zero initial data. The involved integral operators are approximated using the sinc-quadrature formulas that are tailored to the spectral parameters of $A$, fractional order $α$ and the smoothness of the first initial condition, as well as to the properties of the equation's right-hand side $f(t)$. The resulting method possesses exponential convergence for positive sectorial $A$, any finite $t$, including $t = 0$ and the whole range $α\in (0,2)$. It is suitable for a practically important case, when no knowledge of $f(t)$ is available outside the considered interval $t \in [0, T]$. The algorithm of the method is capable of multi-level parallelism. We provide numerical examples that confirm the theoretical error estimates. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2304_13099 |
| institution | arXiv |
| publishDate | 2023 |
| record_format | arxiv |
| spellingShingle | Exponentially Convergent Numerical Method for Abstract Cauchy Problem with Fractional Derivative of Caputo Type Sytnyk, Dmytro Wohlmuth, Barbara Numerical Analysis Mathematical Software Analysis of PDEs Classical Analysis and ODEs 34A08, 35R11, 34G10, 35R20, 65L05, 65J08, 65J1 We present an exponentially convergent numerical method to approximate the solution of the Cauchy problem for the inhomogeneous fractional differential equation with an unbounded operator coefficient and Caputo fractional derivative in time. The numerical method is based on the newly obtained solution formula that consolidates the mild solution representations of sub-parabolic, parabolic and sub-hyperbolic equations with sectorial operator coefficient $A$ and non-zero initial data. The involved integral operators are approximated using the sinc-quadrature formulas that are tailored to the spectral parameters of $A$, fractional order $α$ and the smoothness of the first initial condition, as well as to the properties of the equation's right-hand side $f(t)$. The resulting method possesses exponential convergence for positive sectorial $A$, any finite $t$, including $t = 0$ and the whole range $α\in (0,2)$. It is suitable for a practically important case, when no knowledge of $f(t)$ is available outside the considered interval $t \in [0, T]$. The algorithm of the method is capable of multi-level parallelism. We provide numerical examples that confirm the theoretical error estimates. |
| title | Exponentially Convergent Numerical Method for Abstract Cauchy Problem with Fractional Derivative of Caputo Type |
| topic | Numerical Analysis Mathematical Software Analysis of PDEs Classical Analysis and ODEs 34A08, 35R11, 34G10, 35R20, 65L05, 65J08, 65J1 |
| url | https://arxiv.org/abs/2304.13099 |