Efficient numerical solution of the time fractional diffusion equation by mapping from its Brownian counterpart

Abstract

The solution of a Caputo time fractional diffusion equation of order $0<\alpha <1$ is expressed in terms of the solution of a corresponding integer order diffusion equation. We demonstrate a linear time mapping between these solutions that allows for accelerated computation of the solution of the fractional order problem. In the context of an $N$-point finite difference time discretisation, the mapping allows for an improvement in time computational complexity from $O\left(N^2\right)$ to $O\left(N^{\alpha }\right)$, given a precomputation of $O\left(N^{1+\alpha }\ln N\right)$. The mapping is applied successfully to the least-squares fitting of a fractional advection diffusion model for the current in a time-of-flight experiment, resulting in a computational speed up in the range of one to three orders of magnitude for realistic problem sizes.

9 pages, 5 figures; added references for section 2

References in corpus (7)

• Fractional calculus and continuous-time finance
• From Diffusion to Anomalous Diffusion: A Century after Einstein's Brownian Motion
• Observation of anomalous diffusion and fractional self-similarity in one dimension
• Short note on the emergence of fractional kinetics
• A Fractional Lie Group Method For Anomalous Diffusion Equations
• Analytic solution of the fractional advection diffusion equation for the time-of-flight experiment in a finite geometry
• Anomalous diffusion in nonhomogeneous media: Time-subordinated Langevin equation approach

Cited by in corpus (3)

• Bernstein dual-Petrov-Galerkin method: application to 2D time fractional diffusion equation
• Solution of a generalised Boltzmann's equation for non-equilibrium charged particle transport via localised and delocalised states
• Banded operational matrices for Bernstein polynomials and application to the fractional advection-dispersion equation