Compact difference scheme for the fractional sub-diffusion equation with Neumann boundary conditions

  • Authors:
  • Jincheng Ren;Zhi-Zhong Sun;Xuan Zhao

  • Affiliations:
  • Department of Mathematics, Southeast University, Nanjing 210096, PR China and Department of Mathematics, Shangqiu Normal University, Henan Shangqiu 476000, PR China;Department of Mathematics, Southeast University, Nanjing 210096, PR China;Department of Mathematics, Southeast University, Nanjing 210096, PR China

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

Quantified Score

Hi-index 31.46

Visualization

Abstract

An effective finite difference scheme is considered for solving the time fractional sub-diffusion equation with Neumann boundary conditions. A difference scheme combining the compact difference approach the spatial discretization and L"1 approximation for the Caputo fractional derivative is proposed and analyzed. Although the spatial approximation order at the Neumann boundary is one order lower than that for interior mesh points, the unconditional stability and the global convergence order O(@t^2^-^@a+h^4) in discrete L"2 norm of the compact difference scheme are proved rigorously, where @t is the temporal grid size and h is the spatial grid size. Numerical experiments are included to support the theoretical results, and comparison with the related works are presented to show the effectiveness of our method.