Numerical approximation of Caputo-type advection-diffusion equations via Sylvester equations
Fuente:
arXiv
Saved in:
| Main Authors: | , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866910786652209152 |
|---|---|
| author | de la Hoz, Francisco Muniain, Peru |
| author_facet | de la Hoz, Francisco Muniain, Peru |
| contents | In this paper, we approximate numerically the solution of Caputo-type advection-diffusion equations of the form $D_t^α u(t,x) = a_1(x)u_{xx}(t,x) + a_2(x)u_x(t,x) + a_3u(t,x) + a_4(t,x)$, where $D_t^α u$ denotes the Caputo fractional derivative of order $α\in(0,1)$ of $u$ with respect to $t$, $t\in[0, t_f]$ and the spatial domain can be the whole real line or a closed interval. First, we propose a method of order $3 - α$ to approximate Caputo fractional derivatives, explain how to implement an FFT-based fast convolution to reduce the computational cost, and express the numerical approximation in terms of an operational matrix. Then, we transform a given Caputo-type advection-diffusion equation into a Sylvester equation of the form $\mathbf A\mathbf U + \mathbf U \mathbf B = \mathbf C$, and special care is given to the treatment of the boundary conditions, when the spatial domain is a closed interval. Finally, we perform several numerical experiments that illustrate the adequacy of our approach. The implementation has been done in Matlab, and we share and explain in detail the majority of the actual codes that we have used. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2501_09180 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Numerical approximation of Caputo-type advection-diffusion equations via Sylvester equations de la Hoz, Francisco Muniain, Peru Numerical Analysis 15A24, 26A33, 65T50, 35K57 In this paper, we approximate numerically the solution of Caputo-type advection-diffusion equations of the form $D_t^α u(t,x) = a_1(x)u_{xx}(t,x) + a_2(x)u_x(t,x) + a_3u(t,x) + a_4(t,x)$, where $D_t^α u$ denotes the Caputo fractional derivative of order $α\in(0,1)$ of $u$ with respect to $t$, $t\in[0, t_f]$ and the spatial domain can be the whole real line or a closed interval. First, we propose a method of order $3 - α$ to approximate Caputo fractional derivatives, explain how to implement an FFT-based fast convolution to reduce the computational cost, and express the numerical approximation in terms of an operational matrix. Then, we transform a given Caputo-type advection-diffusion equation into a Sylvester equation of the form $\mathbf A\mathbf U + \mathbf U \mathbf B = \mathbf C$, and special care is given to the treatment of the boundary conditions, when the spatial domain is a closed interval. Finally, we perform several numerical experiments that illustrate the adequacy of our approach. The implementation has been done in Matlab, and we share and explain in detail the majority of the actual codes that we have used. |
| title | Numerical approximation of Caputo-type advection-diffusion equations via Sylvester equations |
| topic | Numerical Analysis 15A24, 26A33, 65T50, 35K57 |
| url | https://arxiv.org/abs/2501.09180 |