Three commonly used two equation turbulence models were adopted. The computed velocity field was compared with the data from corresponding laser Doppler anemometer measurement to analyze the accuracy of these models. By adjusting the model parameters, the simulation data can fit into the experimental result. The optimal viscosity coefficients Cμin these used eddy viscosity models were presented. It is shown that the RNG k-ε model is more sensitive to Cμ, while the results from Standard k-ε model vary almost linearly with the change of the Cμ value. Study on the simulated flow field gives that if the peak of turbulence kinetic energy appears within the radial span, it will affect the velocity distribution along the jet axis; otherwise, it will modify the velocity magnitude. The proposed numerical scheme reduces the physical complexities involved in ultra high pressure injection process so that it can be applied to archive fast production from relevant CAE workflow.