An accurate LOD scheme for two-dimensional parabolic problems

We propose a high order locally one-dimensional scheme for solving parabolic problems. The method is fourth-order in space and second-order in time, and provides a computationally efficient implicit scheme. It is shown through a discrete Fourier analysis that the method is unconditionally stable. Numerical experiments are conducted to test its high accuracy and to compare it with other schemes.