The elliptic conceptual second moment models for turbulent heat fluxes, which are proposed on the basis of elliptic-blending and elliptic-relaxation equations, are applied to calculate the combined forced and natural turbulent convection in a vertical plane channel. The models satisfy the near-wall balance between viscous diffusion, viscous dissipation and temperature-pressure gradient correlation, and also have the characteristics of approaching its respective conventional high Reynolds number model far away from the wall. Also the models are closely linked to the elliptic blending model which is used for the prediction of Reynolds stress. In order to calibrate the heat flux models, firstly, the distributions of mean temperature and scala flux in fully developed channel flow with constant wall difference temperature are solved by the present models. The buoyancy effect on the turbulent characteristics including the mean velocity and temperature, the Reynolds stress tensor, and the turbulent heat flux vector are examined. In the opposing flow, the turbulent transport is greatly enhanced with both the Reynolds stresses and the turbulent heat fluxes being remarkably increased; whereas, in the aiding flow, the opposite change is observed. The results of prediction are directly compared to the DNS to assess the performance of the model predictions and show that the behaviors of the turbulent heat transfer in the whole flow region are well captured by the present models.