This study used computational fluid dynamics for low greenhouse effect refrigerant (R-513A) simulation analysis in the two-stage 150 USRT class oil-free centrifugal refrigerant compressor using integrated part load value (IPLV) and internal flow field. The compressor rotor speed and mass flow rate for the IPLV working conditions with various loads were planned using Stage 1 and Stage 2 simulations, respectively. The performance and flow field numerical simulation analyses for the two-stage centrifugal compressor is discussed. This study used Ansys-CFX software for numerical simulation analysis and the conservation form of a 3D steady-state Navier–Stokes equation set with the finite volume discretization method for computation. The computing mode produced better computing efficiency and flexible mesh setting using the k-omega (k-ω) model, which has better computational performance in the near wall boundary layer and low Reynolds number flow field (used as the turbulence model) for simulation. The R-513A refrigerant property setting was calculated using the Soave–Redlich–Kwong gas equation. This study discusses the shaft work, pressure ratio, and isentropic efficiency and also describes the main research findings with the meridional pressure, isentropic efficiency contour maps, and flow field velocity vector diagrams. According to the numerical simulation results, in Stage 1 and Stage 2 simulations, the isentropic efficiency produced the highest numerical results in the 75% load case, which are 88.19% and 89.06%, respectively. The isentropic efficiency decreased to 75.93% and 82.26%, respectively, in the 25% load case. The flow field velocity vector diagram shows that in the 25% load case, there was back-flow field distribution near the impeller shroud. The compressor performance was also analyzed.