Turbulent thermal convection between two plates, bottom of which is at higher temperature $T_h$ and the upper plate is at lower temperature $T_1$, is numerically investigated. Model equations are abridged Reynolds stress equations; full Reynolds stress equations are simplified to yield algebraic relations.
Boundary conditions for K and $\bar{θ^2}$ at the plate surface are set to be zero and those of ε and $ε_θ$ are assumed at first grid point nearest to the boundary surface, whose values are approximated by inviscid estimates.
Results show that temperature profiles are in good agreement with experimental data except, in a so called, transition region, in which temperature is overpredicted. Such discrepancy is larger for lower Rayleigh number.
Nusselt numbers, which are calculated from the temperature gradients at the boundary surface, are in good agreement with experimental data.
