Fast, Accurate and Robust Adaptive Finite Difference Methods for Fractional Diffusion Equations: The Size of the Timesteps does Matter

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Yuste, Santos B., Quintana-Murillo, Joaquin
Format: Preprint
Publié: 2014
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866910502918029312
author Yuste, Santos B.
Quintana-Murillo, Joaquin
author_facet Yuste, Santos B.
Quintana-Murillo, Joaquin
contents The computation time required by standard finite difference methods with fixed timesteps for solving fractional diffusion equations is usually very large because the number of operations required to find the solution scales as the square of the number of timesteps. Besides, the solutions of these problems usually involve markedly different time scales, which leads to quite inhomogeneous numerical errors. A natural way to address these difficulties is by resorting to adaptive numerical methods where the size of the timesteps is chosen according to the behaviour of the solution. A key feature of these methods is then the efficiency of the adaptive algorithm employed to dynamically set the size of every timestep. Here we discuss two adaptive methods based on the step-doubling technique. These methods are, in many cases, immensely faster than the corresponding standard method with fixed timesteps and they allow a tolerance level to be set for the numerical errors that turns out to be a good indicator of the actual errors.
format Preprint
id arxiv_https___arxiv_org_abs_1403_4434
institution arXiv
publishDate 2014
record_format arxiv
spellingShingle Fast, Accurate and Robust Adaptive Finite Difference Methods for Fractional Diffusion Equations: The Size of the Timesteps does Matter
Yuste, Santos B.
Quintana-Murillo, Joaquin
Numerical Analysis
The computation time required by standard finite difference methods with fixed timesteps for solving fractional diffusion equations is usually very large because the number of operations required to find the solution scales as the square of the number of timesteps. Besides, the solutions of these problems usually involve markedly different time scales, which leads to quite inhomogeneous numerical errors. A natural way to address these difficulties is by resorting to adaptive numerical methods where the size of the timesteps is chosen according to the behaviour of the solution. A key feature of these methods is then the efficiency of the adaptive algorithm employed to dynamically set the size of every timestep. Here we discuss two adaptive methods based on the step-doubling technique. These methods are, in many cases, immensely faster than the corresponding standard method with fixed timesteps and they allow a tolerance level to be set for the numerical errors that turns out to be a good indicator of the actual errors.
title Fast, Accurate and Robust Adaptive Finite Difference Methods for Fractional Diffusion Equations: The Size of the Timesteps does Matter
topic Numerical Analysis
url https://arxiv.org/abs/1403.4434