Introduction. Calculation of dynamics of the anisothermal moisture transfer processes in axisymmetric formulation is essential in the study of wet soils condition around, for example, vertical drains, wells, piles, etc. In this paper, we formulate the initial boundary value problem for the system of moisture and heat transfer nonstationary equations. The problem is considered for isotropic medium in cylindrical coordinate system under the inhomogeneous mixed boundary conditions. The obtained results are important for future research in cylindrical coordinates of problems that model the migration of moisture during the seasonal freezing of the soil, taking into account phase transitions from unfrozen water to ice in the entire volume of the soil mass without highlighting the crystallization front. In this case moisture exchange and heat transfer characteristics appear as functions of the total humidity. Consequently, the equation of moisture transfer is written relative to the "fictitious" moisture content. Because of the main direction of moisture migration relative to the freezing/melting front, the convective heat transfer along the vertical coordinate axis is considered to be essential that leads to sufficient coincidence with the experimental data. The purpose of the paper is to formulate the appropriate generalized problem in the Galorkin form for the axisymmetric initial-boundary value problem. The important goal is to investigate the accuracy of the continuous in time and completely discrete approximate generalized solutions based on the finite elements method. Results. The algorithm for constructing of approximate generalized solution of the axisymmetric initial-boundary value problem for the system of filtration and heat transfer equations is proposed. The estimates of the convergence rate for the continuous in time and discrete approximate solutions based on the finite elements method are obtained.