The primary objective of the present study is evaluation of the k-ε-vv-f turbulence model for prediction of natural convection in a rectangular cavity. As a comparative study, the two-layer k-ε model is also considered. Both models, with and without algebraic heat flux model, are applied to the analysis of natural convection in a rectangular cavity. The performances of turbulence models are investigated through comparison with available experimental data. The predicted results of vertical velocity component, turbulent heat fluxes, turbulent shear stress, local Nusselt number and wall shear stress are compared with experimental data. It is shown that, among the turbulence models considered in the present study, the k-ε-vv-f model with an algebraic heat flux model predicts best the vertical mean velocity and velocity fluctuation, and the inclusion of algebraic heat flux model slightly improves the accuracy of results.