A difference scheme for a conjugate-operator model of the heat conduction problem in the polar coordinates

In the polar coordinates, a discrete analog of the conjugate-operator model of the heat conduction problem preserves the structure of the original model. The difference scheme converges with the second order of accuracy for the cases of discontinuous parameters of the medium in the Fourier law and irregular grids. An efficient algorithm for solving the discrete conjugate-operator model in the case when the thermal conductivity tensor is a single operator.