In this paper, the results of numerical analysis of the thermal consolidation of a two phase medium, under the assumption of independent heat transfer in fluid and the solid phase of the medium, are presented. Three cases of pore fluid were considered: liquid, represented by water, and gas, represented by air and carbon dioxide. The mathematical model was derived from irreversible thermodynamics, with the assumption of a constant heat transfer between the phases. In the case of the accepted geometry of the classical dimensions of the soil sample and boundary conditions, the process leads to equalization of temperatures of the skeleton on the pore fluid. Heat transfer is associated with the fluid flow in the pores of the medium. In the case of gas as the pore fluid, a non-linear mathematical model of gas filtration through the pores of the medium was accepted. For the computing process, relationships between viscosity or density and temperature proposed by other authors were taken into account. Despite accepting mechanical constants of the solid phase that do not depend on temperature, the obtained model is nonlinear and develops the classical Biot–Darcy model.