摘要

In this paper, a new kind of alternating direction implicit (ADI) Crank-Nicolsontype orthogonal spline collocation (OSC) method is formulated for the two-dimensional fractional evolution equation with a weakly singular kernel arising in the theory of linear viscoelasticity. The novel OSC method is used for the spatial discretization, and ADI Crank-Nicolsontype method combined with the second order fractional quadrature rule are considered for the temporal component. The stability of proposed scheme is rigourously established, and nearly optimal order error estimate is also derived. Numerical experiments are conducted to support the predicted convergence rates and also exhibit expected super-convergence phenomena.