A novel Douglas alternating direction implicit (ADI) method is proposed in this work to solve a two-dimensional (2D) heat equation with interfaces. The ADI scheme is a powerful finite difference method for solving parabolic equations, due to its unconditional stability and high efficiency. However, it suffers from a serious accuracy reduction in space for interface problems with different mater...