A Finite Difference Method with Non-uniform Timesteps for Fractional Diffusion Equations

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Yuste, Santos B., Quintana-Murillo, Joaquín
Format: Preprint
Published: 2011
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909231935913984
author Yuste, Santos B.
Quintana-Murillo, Joaquín
author_facet Yuste, Santos B.
Quintana-Murillo, Joaquín
contents An implicit finite difference method with non-uniform timesteps for solving the fractional diffusion equation in the Caputo form is proposed. The method allows one to build adaptive methods where the size of the timesteps is adjusted to the behaviour of the solution in order to keep the numerical errors small without the penalty of a huge computational cost. The method is unconditionally stable and convergent. In fact, it is shown that consistency and stability implies convergence for a rather general class of fractional finite difference methods to which the present method belongs. The huge computational advantage of adaptive methods against fixed step methods for fractional diffusion equations is illustrated by solving the problem of the dispersion of a flux of subdiffusive particles stemming from a point source.
format Preprint
id arxiv_https___arxiv_org_abs_1109_6622
institution arXiv
publishDate 2011
record_format arxiv
spellingShingle A Finite Difference Method with Non-uniform Timesteps for Fractional Diffusion Equations
Yuste, Santos B.
Quintana-Murillo, Joaquín
Numerical Analysis
Statistical Mechanics
Computational Physics
An implicit finite difference method with non-uniform timesteps for solving the fractional diffusion equation in the Caputo form is proposed. The method allows one to build adaptive methods where the size of the timesteps is adjusted to the behaviour of the solution in order to keep the numerical errors small without the penalty of a huge computational cost. The method is unconditionally stable and convergent. In fact, it is shown that consistency and stability implies convergence for a rather general class of fractional finite difference methods to which the present method belongs. The huge computational advantage of adaptive methods against fixed step methods for fractional diffusion equations is illustrated by solving the problem of the dispersion of a flux of subdiffusive particles stemming from a point source.
title A Finite Difference Method with Non-uniform Timesteps for Fractional Diffusion Equations
topic Numerical Analysis
Statistical Mechanics
Computational Physics
url https://arxiv.org/abs/1109.6622