摘要

In this article, we consider the finite element methods (FEM) for Grwunwald-Letnikov time-fractional diffusion equation, which is obtained from the standard two-dimensional diffusion equation by replacing the first-order time derivative with a fractional derivative (of order , with 0<<1). The proposed method is based on high-order FEM for spacial discretization and finite difference method for time discretization. We prove that the method is unconditionally stable, and the numerical solution converges to the exact one with order O(h(r+1)+(2-)), where h, and r are the space step size, time step size and polynomial degree, respectively. A numerical example is presented to verify the order of convergence.

全文