In northern seas the gas temperature in the pipeline may be lower than the water-ice phase transition temperature, so the glaciation process must be considered. The coefficient smoothing method for the problem of glaciation of the cylinder immersed in seawater is considered. The model of glaciation process is presented as a problem for the linear heat equation in domain with unknown moving boundary. The method of solution is based on the transition to the Dirichlet problem for the nonlinear two-dimensional heat equation in domain with fixed boundaries. The splitting method is applied to the solution of the problem for two-dimensional heat equation. Two approximations for the Dirac delta function are proposed. Calculation of uniform glaciation process was carried out with the using of nonlinear implicit finite-difference schemes. Time moments of the glaciation for different spatial layers are obtained. Results of calculations are compared with the results obtained by the front-tracking method.