Skip to main content
Article
A Radial Basis Functions method for fractional diffusion equations
Journal of Computational Physics (2013)
  • Cecile M Piret, Michigan Technological University
  • Emmanuel Hanert, Universit√© catholique de Louvain
Abstract
One of the ongoing issues with fractional diffusion models is the design of an efficient high order numerical discretization. This is one of the reasons why fractional diffusion models are not yet more widely used to describe complex systems. In this paper, we derive a radial basis functions (RBF) discretization of the one-dimensional space-fractional diffusion equation. In order to remove the ill-conditioning that often impairs the convergence rate of standard RBF methods, we use the RBF-QR method [1,33]. By using this algorithm, we can analytically remove the ill-conditioning that appears when the number of nodes increases or when basis functions are made increasingly flat. The resulting RBF-QR-based method exhibits an exponential rate of convergence for infinitely smooth solutions that is comparable to the one achieved with pseudo-spectral methods. We illustrate the flexibility of the algorithm by comparing the standard RBF and RBF-QR methods for two numerical examples. Our results suggest that the global character of the RBFs makes them well-suited to fractional diffusion equations. They naturally take the global behavior of the solution into account and thus do not result in an extra computational cost when moving from a second-order to a fractional-order diffusion model. As such, they should be considered as one of the methods of choice to discretize fractional diffusion models of complex systems.
Keywords
  • Anomalous Diffusion,
  • Radial Basis Functions,
  • Fractional Derivative,
  • Fractional Diffusion
Publication Date
April 1, 2013
DOI
10.1016/j.jcp.2012.10.041
Publisher Statement
© 2013 Journal of Computational Physics. Deposited in compliance with publisher policies. Publisher's version of record: http://www.sciencedirect.com/science/article/pii/S0021999112006444
Citation Information
Cecile M Piret and Emmanuel Hanert. "A Radial Basis Functions method for fractional diffusion equations" Journal of Computational Physics Vol. 238 (2013) p. 71 - 81
Available at: http://works.bepress.com/cecile_piret/2/