This study employed Computational Fluid Dynamics, CFD solution to configure the entire shape of the nozzle. The problems arising from the initial geometrical shape obtained using one-dimensional approach were highlighted. Then with a two-dimensional approach and coded FORTRAN programs, the convergent portion was handled with time matching CFD while the divergent portion was handled with space matching CFD. A post-processing results analysis shows that the results are in agreement with that obtained with the analytical approach. This research has demonstrated significant and tangible benefits of the use of CFD numerical experimentation to optimize the shape of the nozzle. These benefits are not solely limited to performance enhancements, but solution reliability and algorithm development.