The present paper addresses a comparative analysis of twelve different versions of low Reynolds number k-epsilon turbulence models. The predictive capability of the models have been tested on the basis of the flow patterns and energy balance. Numerical simulations were performed at the Reynolds numbers of 7400, 22 000 and 500 000. The predicted mean axial velocity and turbulent kinetic energy were compared with the experimental data of Durst et al. (1995) and Schildknecht et al. (1979) for the Reynolds number of 7400 and 22 000 respectively. The overall energy balance was established at three Reynolds numbers of 7400, 22 000 and 500 000. A comparison of all the models has been predicted.