This study used Computational Fluid Dynamics (CFD) to simulate and analyze the working fluid in magnetic centrifugal refrigerant compressors using R-134a to mixed refrigerant: R-513A and HFO (Hydrofluoroolefins) Hydrofluoroolefin refrigerant: R-1234yf, and the impact on integrated part-load performance, Integrated Part Load Value (IPLV) and internal flow field. This study used a single-stage 280 USRT maglev centrifugal refrigerant compressor as a simulation model. Three different refrigerants were used: R-134a, R-513A, and R-1234yf, as presented in the National Institute of Standards and Technology (NIST) real gas database. The refrigerant was used to set the IPLV working conditions and change the compressor speed and mass flow rate to simulate the compressor’s characteristic curve after replacing the refrigerant. The compressor working conditions were the fixed refrigeration cycle condensation and evaporation following the same capacity standards. This study used the CFD software by Ansys software company to simulate the flow field. The k-omega turbulence software was used to model the turbulence. The results show that the maglev centrifugal refrigerant compressor efficiency dropped significantly when the refrigerant was directly replaced. Based on R-134a, the full load efficiency of R-1234yf dropped 13.21%, the full load efficiency of R-513A dropped 9.97%, and the partial load efficiency was similar to R-134a.

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.