Assessment and validation of RANS turbulence models are conducted for the optimal analysis of supersonic converging-diverging nozzle through the comparison between computational results and experimental data. One/two equation turbulence closures such as Spalart-Allmaras, RNG k-${varepsilon}$, and k-${omega}$ SST are employed to simulate the two-dimensional nozzle flow. Computational results with the turbulence models mentioned fairly well predict shock structure of the nozzle-inside and pressure distribution along the wall. Especially, SST model among the employed ones shows the best agreement to experimental results.