A computational study on a strongly stratified natural convection is performed with the elliptic blending second-moment closure. The turbulent heat flux is treated by both the algebraic flux model (AFM) and the differential flux model (DFM). Calculations are performed for a turbulent natural convection in a square cavity with conducting top and bottom walls and the calculated results are compared with the available experimental data. The results show that both the AFM and DFM models produce very accurate solutions with the elliptic-blending second-moment closure without invoking any numerical stability problems. These results show that the AFM and DFM models for treating the turbulent heat flux are sufficient for this strongly stratified flow. However, a slight difference between two models is observed for some variables.