In the present work, five different subgrid-scale (SGS) models and implicit large-eddy simulation (LES) are evaluated and compared against the DNS for the simulation of the planar turbulent wall-bounded jet with heat transfer. The SGS models tested are the classical constant coefficient Smagorinsky model and its dynamic version, the wall-adaptive local eddy-viscosity (WALE) model, the turbulent kinetic energy one-equation model, and its dynamic version. The effects of using variable turbulent Prandtl number and the near-wall damping function are also studied in these models. The mean, second-order flow and heat-transfer statistics with the evolution of Nusselt number along the jet downstream are used to assess the different SGS models. The quality of resolution of the present LES are evaluated using the activity parameter and the index of resolution quality. Among the models tested, the constant coefficient Smagorinsky model together with Van-Driest damping predicts the solution accurately in the near-wall region as well as in estimating the thermal parameters. However, the dynamic models performed better in evaluating the Reynolds stress profiles away from the wall in the outer region. Capabilities of the models to predict the turbulent kinetic energy budgets, pressure-velocity gradient correlations and triple velocity correlations are also studied. The implemented variable Prandtl number algorithm is noted to have minimal influence on the evolution of the solution.