An O(N log2N) alternating-direction finite difference method for two-dimensional fractional diffusion equations

  • Authors:

  • Hong Wang;Kaixin Wang

  • Affiliations:

  • School of Mathematics, Shandong University, Jinan, Shandong 250100, China and Department of Mathematics, University of South Carolina, Columbia, SC 29208, USA;School of Mathematics, Shandong University, Jinan, Shandong 250100, China

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2011

Quantified Score

Hi-index 31.48

Visualization

Abstract

Fractional diffusion equations model phenomena exhibiting anomalous diffusion that cannot be modeled accurately by the second-order diffusion equations. Because of the nonlocal property of fractional differential operators, the numerical methods for fractional diffusion equations often generate dense or even full coefficient matrices. Consequently, the numerical solution of these methods often require computational work of O(N^3) per time step and memory of O(N^2) for where N is the number of grid points. In this paper we develop a fast alternating-direction implicit finite difference method for space-fractional diffusion equations in two space dimensions. The method only requires computational work of O(N log^2N) per time step and memory of O(N), while retaining the same accuracy and approximation property as the regular finite difference method with Gaussian elimination. Our preliminary numerical example runs for two dimensional model problem of intermediate size seem to indicate the observations: To achieve the same accuracy, the new method has a significant reduction of the CPU time from more than 2 months and 1 week consumed by a traditional finite difference method to 1.5h, using less than one thousandth of memory the standard method does. This demonstrates the utility of the method.