Alternating direction implicit-Euler method for the two-dimensional fractional evolution equation

  • Authors:
  • Limei Li;Da Xu

  • Affiliations:
  • -;-

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

Quantified Score

Hi-index 31.45

Visualization

Abstract

In this paper, we present a fast and efficient numerical method to solve a two-dimensional fractional evolution equation on a finite domain. This numerical method combines the alternating direction implicit (ADI) approach with the second-order difference quotient in space, the backward Euler in time and order one convolution quadrature approximating the integral term. By using the discrete energy method, we prove that the ADI scheme is unconditionally stable and the numerical solution converges to the exact one with order O(k+h"x^2+h"y^2), where k is the temporal grid size and h"x,h"y are spatial grid sizes in the x and y directions, respectively. Two numerical examples with known exact solution are also presented, and the behavior of the error is analyzed to verify the order of convergence of the ADI-Euler method.