Convergence analysis of the finite difference ADI scheme for variable coefficient parabolic problems with nonzero Dirichlet boundary conditions

Since the invention by Peaceman and Rachford, more than 60 years ago, of the well celebrated ADI finite difference scheme for parabolic initial-boundary problems on rectangular regions, many papers have been concerned with prescribing the boundary values for the intermediate approximations at half time levels in the case of nonzero Dirichlet boundary conditions. In the present paper, for variable coefficient parabolic problems and time-stepsize sufficiently small, we prove second order accuracy in the discrete $L^2$ norm of the ADI finite difference scheme in which the intermediate approximations do not involve the so called “perturbation term”. As a byproduct of our stability analysis we also show that, for variable coefficients and time-stepsize sufficiently small, the ADI scheme with the perturbation term converges with order two in the discrete $H^1$ norm. Our convergence results generalize previous results obtained for the heat equation.