Numerical algorithm for solving diffusion equations on the basis of multigrid methods

A new effective algorithm based on multigrid methods is proposed for solving parabolic equations. The algorithm preserves implicit-scheme advantages (such as stability, accuracy, and conservativeness) while it involves a considerably reduced amount of arithmetic operations at every time level. The absolute stability, conservativeness, and convergence of the algorithm is proved theoretically using one- and two-dimensional initial-boundary value model problems for the heat equation. The error of the solution is estimated. The good accuracy of the method is demonstrated using two-dimensional model problems, including ones with discontinuous coefficients.