The diffuser of a reactor coolant pump was optimized using an orthogonal approach with numerical simulation to improve the pump hydraulic performance. Steady simulation was conducted by solving Reynolds-averaged Naiver-Stokes equations with the SST k-omega turbulence model using CFX code. The influence of the diffuser geometric parameters, namely, S, phi, alpha (4), b (4), delta (2), R (t) and R (4), on the pump performance were determined. L-18 (3(7)) orthogonal table was chosen for the optimization process. Best indicators were determined, and range analysis of energy losses, head, and efficiency at the rated condition was performed. Optimal parameters of the diffuser were S = 490 mm, phi = 36A degrees, alpha (4) = 30A degrees, b (4) = 200 mm, delta (2) = 20 mm, R (t) = 5 mm and R (4) = 565 mm. The final design was experimentally tested. Simulation results showed more remarkable performance than the experimental result. However, the numerical predictions and experimental results were consistent, validating the design procedure. Loading of the impeller and diffuser blades was analyzed to investigate the direct impact on the hydrodynamic flow field. The head was 14.74 m, efficiency was 79.6 %, and efficiency of the prototype pump was 83.3 % when the model pump functioned at the rated conditions. Optimization results showed that efficiency and head were improved at the design condition.