In order to study the flow behavior of multiple jets, numerical prediction of the three-dimensional domain of round jets from the nozzle edge up to the turbulent region is essential. The previous numerical studies on the round jet are limited to either two-dimensional investigation with Reynolds-averaged Navier–Stokes (RANS) models or three-dimensional prediction with higher turbulence models such as large eddy simulation (LES) or direct numerical simulation (DNS). The present study tries to evaluate different RANS turbulence models in the three-dimensional simulation of the whole domain of an isothermal, low Re (Re = 2125, 3461, and 4555), free, turbulent round jet. For this evaluation the simulation results from two two-equation (low Re and low Re shear stress transport (SST) ), a transition three-equation (), and a transition four-equation (SST) eddy-viscosity turbulence models are compared with hot-wire anemometry measurements. Due to the importance of providing correct inlet boundary conditions, the inlet velocity profile, the turbulent kinetic energy (), and its specific dissipation rate () at the nozzle exit have been employed from an earlier verified numerical simulation. Two-equation RANS models with low Reynolds correction can predict the whole domain (initial, transition, and fully developed regions) of the round jet with prescribed inlet boundary conditions. The transition models could only reach to a good agreement with the measured mean axial velocities and its rms in the initial region. It worth mentioning that the round jet anomaly is still present in the turbulent region of the round jet predicted by the low Re . By comparing the and the predicted by different turbulence models, the blending functions in the cross-diffusion term is found one of the reasons behind the more consistent prediction by the low Re SST .