A comparative study on the treatment of the turbulent heat flux with the elliptic blending second-moment closure for a natural convection flow is performed. Three cases of different treating the turbulent heat flux are considered. Those are the generalized gradient diffusion hypothesis (GGDH), the algebraic flux model (AFM) and the differential flux model (DFM). The constants in the models are adjusted with a primary emphasis placed on the accuracy of predicting the local Nusselt number. These models are implemented in a computer code specially designed for evaluation of turbulent models. Calculations are performed for a turbulent natural convection in the 1:5 rectangular cavity and the calculated results are compared with the available experimental data. The results show that the three models produce nearly the same accuracy of solutions. These results show that the GGDH, AFM and DFM models for treating the turbulent heat flux are sufficient for this simple shear flow where the shear production is dominant. It is observed that, in the weakly stratified region at the center zone of the cavity, the vertical velocity fluctuation is nearly zero in the GGDH solutions, which shows that the GGDH model may not be suitable for the strongly stratified flow. Thus, further study on the strongly stratified flow should be followed.