The ice slurry heat transfer process in flow can be described as forced convection heat transfer of a non-Newtonian fluid subject to phase change (melting of ice crystals). This paper looks at the application of the enthalpy-porosity method to the modelling of the melting process occurring during homogenous tube flow of ice slurry. It will discuss experimentally verified results of thermal calculations, including the distribution of temperature values, heat transfer coefficients and changes to the content of solid particles. The results of the calculations apply to the laminar and turbulent flow areas, various initial mass fractions of ice and various tube cross-section shapes.