As a preliminary design study to achieve target aerodynamic performance, this work was conducted on an original nozzle with 9 flutes in order to design a fluted nozzle with 12 flutes. The thrust and rolling moment of the nozzle with 12 flutes were analyzed using a CFD code according to the depth and rotation angle of the flutes. Based on this, a fluted nozzle with 12 flutes was optimized to yield the same thrust as that of the original nozzle with 9 flutes. The response surface method was applied for shape optimization of the fluted nozzle and design variables were selected to determine the depth angle and rotation angle of the flutes. An optimized shape that led to a thrust as strong as that of the original nozzle was obtained.