In this study, the performance of a pump-jet propulsor (PJP) with pre-swirl stator in open water is numerically investigated. Both full-scale and model-scale configurations are considered. The Reynolds-averaged Navier–Stokes equations and shear stress transport k−ω turbulence model are used in the numerical calculation. The computational domain is discretized using structured grids, and a rotating grid is affixed to the rotor to deal with the relative motion between the rotor and stationary components. The mesh quality is determined based on a grid uncertainty analysis. The numerical method is validated using model-scale experimental data. The simulation results reveal the influences of the scale size on the hydrodynamic performance and the distributions of the velocity, pressure and vorticity under three advance coefficients. With the increase in the advance coefficients, the scale influences on the efficiency become more obvious, and the efficiency of the full-scale PJP is always higher than that of the model-scale PJP. The full-scale configuration is found with a more significant instability in the gap vortex development, because it presents larger interaction between tip leakage vortex (TLV) and the inner wall of the duct. As the main velocity increases, the TLV shedding is delayed. Finally, the development process of gap vortices is analyzed for the difference operation conditions.