The model of unsteady nonisothermal droplet vaporization of atomized liquid is formulated, taking into account the temperature gradient and liquid recirculation inside the droplet. An efficient and absolutely stable algorithm is developed for further use in numerical simulations of turbulent evaporating sprays.