A new second-moment closure model for turbulent heat fluxes is proposed on the basis of the elliptic equation. The new model satisfies the near-wall balance between viscous diffusion, viscous dissipation and temperature-pressure gradient correlation, and also has the characteristics of approaching its respective conventional high Reynolds number model far away from the wall. The predictions of turbulent heat transfer in a channel flow have been carried out with constant wall heat flux and constant wall temperature difference boundary conditions respectively. The velocity field variables are supplied from the DNS data and the differential equations only fur the mean temperature and the scalar flux are solved by the present calculations. The present model is tested by direct comparisons with the DNS to validate the performance of the model predictions. The prediction results show that the behavior of the turbulent heat fluxes in the whole region is well captured by the present model.