English

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

Computational Physics 2015-04-28 v3 Statistical Mechanics

Abstract

The solution of a Caputo time fractional diffusion equation of order 0<α<10<\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 NN-point finite difference time discretisation, the mapping allows for an improvement in time computational complexity from O(N2)O\left(N^2\right) to O(Nα)O\left(N^\alpha\right), given a precomputation of O(N1+αlnN)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.

Keywords

Cite

@article{arxiv.1408.3246,
  title  = {Efficient numerical solution of the time fractional diffusion equation by mapping from its Brownian counterpart},
  author = {Peter W. Stokes and Bronson Philippa and Wayne Read and Ronald D. White},
  journal= {arXiv preprint arXiv:1408.3246},
  year   = {2015}
}

Comments

9 pages, 5 figures; added references for section 2