A CCD-ADI method for unsteady convection-diffusion equations

With a combined compact difference scheme for the spatial discretization and the Crank-Nicolson scheme for the temporal discretization, respectively, a high-order alternating direction implicit method (ADI) is proposed for solving unsteady two dimensional convection diffusion equations. The method is sixth-order accurate in space and second-order accurate in time. The resulting matrix at each ADI computation step corresponds to a triple-tridiagonal system which can be effectively solved with a considerable saving in computing time. In practice, Richardson extrapolation is exploited to increase the temporal accuracy. The unconditional stability is proved by means of Fourier analysis for two dimensional convection diffusion problems with periodic boundary conditions. Numerical experiments are conducted to demonstrate the efficiency of the proposed method. Moreover, the present method preserves the higher order accuracy for convection-dominated problems. (C) 2013 Elsevier B.V. All rights reserved.