Numerical Algorithm With High Spatial Accuracy for the Fractional Diffusion-Wave Equation With Neumann Boundary Conditions

  • Authors:

  • Jincheng Ren;Zhi-Zhong Sun

  • Affiliations:

  • Department of Mathematics, Southeast University, Nanjing, People's Republic of China 210096 and Department of Mathematics, Shangqiu Normal University, Henan, People's Republic of China 476000;Department of Mathematics, Southeast University, Nanjing, People's Republic of China 210096

  • Venue:
  • Journal of Scientific Computing
  • Year:
  • 2013

Quantified Score

Hi-index 0.00

Visualization

Abstract

A fourth-order compact algorithm is discussed for solving the time fractional diffusion-wave equation with Neumann boundary conditions. The $$L1$$ discretization is applied for the time-fractional derivative and the compact difference approach for the spatial discretization. The unconditional stability and the global convergence of the compact difference scheme are proved rigorously, where a new inner product is introduced for the theoretical analysis. The convergence order is $$\mathcal{O }(\tau ^{3-\alpha }+h^4)$$ in the maximum norm, where $$\tau $$ is the temporal grid size and $$h$$ is the spatial grid size, respectively. In addition, a Crank---Nicolson scheme is presented and the corresponding error estimates are also established. Meanwhile, a compact ADI difference scheme for solving two-dimensional case is derived and the global convergence order of $$\mathcal{O }(\tau ^{3-\alpha }+h_1^4+h_2^4)$$ is given. Then extension to the case with Robin boundary conditions is also discussed. Finally, several numerical experiments are included to support the theoretical results, and some comparisons with the Crank---Nicolson scheme are presented to show the effectiveness of the compact scheme.