This paper considers both separating and impinging turbulent flows with heat transfer. The performance of the k-epsilon and k-omega turbulence models is assessed, especially how the low Reynolds number regions are resolved. In the case of the k-epsilon model, three wall functions and six low Reynolds number near-wall closures are evaluated. The turbulence models are solved in conjunction with the Reynolds-averaged momentum and energy equations using a control volume method and QUICK differencing scheme. The results indicate that the k-omega model demonstrates superior performance for prediction of convection heat transfer in complex turbulent flows and numerically is easy to implement.