-A Bell type nozzle is most commonly used shape for rocket nozzles. This type of nozzle not only offers significant advantages in terms of size and performance over the conical nozzle but also reduces complexity compared to annular nozzles. The nozzle uses the stagnation temperature (To) and stagnation pressure (Po) generated in the combustion chamber to create thrust by accelerating the combustion gases to a high supersonic velocity. The nozzle expansion ratio was governed by the exit velocity. During flight, the jet flow is ideally expanded and adapted to the surrounding flow only during a short period. The rest of the time, the rocket engine operates in off-design conditions. The present work incorporates 2D axisymmetric flow analysis within the bell type nozzle, at design and off-design conditions, by using computational fluid dynamic software GAMBIT 2.4.6 and FLUENT 6.3.26. A computer code, with the use of the method of characteristics and stream function, is developed to define the higher efficiency nozzle contours for analysis. Simulation has been carried out separately for two different flow conditions i.e. cold and hot. Shear Stress Transport k-ω turbulence model has been chosen for flow analysis. The converged solutions captured asymmetric lambda shock in the nozzles at higher nozzle pressure ratios (NPR) for viscous flows. It also predicted aftershock and flow separation depending upon NPR. The strength of the normal shock, at Mach stem in viscous prediction, generally increases with an increase in NPR. Good agreement is observed between predicted simulation and analytical results in terms of shock structure, shock location, the size of normal shock, aftershock, and asymmetric lambda shocks.