Numerical Analysis has been done for the supersonic off-design jet flow due to the pressure difference between the jet and the ambient fluid. The difference of pressure generates an oblique shock or an expansion wave at the nozzle exit. The waves reflect repeatedly on the center axis and the sonic surface in the shear layer. The pressure difference is resolved across these reflected waves. In this paper, the axi-symmetric Navier-Stokes equation has been used with the κ-ε turbulence model. The second order TVD scheme with flux limiters, based on the flux vector split with the smooth eigenvalue split, has been used to capture internal shocks and other discontinuities. Numerical calculations have been done to analyze the off-design jet flow due to the pressure difference. The variation of pressure along the flow axis is compared with an experimental result and other numerical result. The characteristics of the interaction between the shock cell and the turbulence mixing layer have been analyzed.