This paper focuses on CFD simulations of the internal flow field of jet pumps, which are passive pumping devices.A numerical method based on the RANS equations is proposed, and its abilities are analysed by comparing the obtained results with the available experimental data.The turbulent flow is simulated using Boussinesq hypothesis which relates Reynolds stresses to the mean velocity gradients via turbulent viscosity.The k-ε model (along with the Standard wall function) and the k-ω SST model are used to calculate the turbulent viscosity for solving the flow in near-wall regions.Before checking the validity of the numerical results, Grid Convergence Index (GCI) is estimated to evaluate the mesh.The analysis showed that the maximum GCI is at the cell near the wall with the value of 0.08%.Comparison of the numerical results and experimental data shows that the model captures the jet pump efficiency range with average relative errors of 7.3% and 8.47% for the k-ε and k-ω SST turbulence models, respectively.The numerical results also confirmed that increasing the flow ratio causes the mixing location between the primary and secondary flows to move toward the outlet of the jet pump, which was observed in the previous experimental investigations.The static pressure coefficient, which indicates the pumping effect of the jet pump, is calculated via both of the turbulence models and the comparison with the experimental data showed that the average relative error is 10.68% and 3.21% for the k-ε model and the k-ω SST model, respectively.The study shows that accurate modelling of jet pump parameters is feasible using RANS approach.