The prediction of distillation zone is very important in steam turbine blades and steam nozzles. In identification of distillery with equilibrium method, as the steam flow contacts the two-phase dome, the second phase formes and flow properties will pass the distillery without any jumping, therefor after crossing the saturation curve, the droplet formation transpires, but in non-equilibrium method by a sudden increase in pressure, called “condensation shock” a discontinuity in the flow characteristics is seen and after crossing the saturation curve, the formation of droplets starts. In this paper, numerical analysis of a vapor-liquid two-phase transonic flow in a convergent-divergent nozzle with and without shock is investigated. Effects of stagnation temperature at nozzle inlet, viscosity and geometry is studied using thermodynamic equilibrium and non-equilibrium methods and results compared with experimental datas. Roe numerical method is used for vapor-liquid two-phase flow numerical solution. The main properties of the flow at the boundary of elements is extrapolated by MUSCL third order acuracy and time discretization is performed using Lax-Wendroff explicit two-step method of second order accuracy. It is observed that the results of non-equilibrium solution, has more correspondence to experimental results and Condensation starts earlier in the nozzle with further expansion rate. By increasing the temperature at nozzle inlet, the place at which condensation starts goes forward. Also in comparision with non-viscous flow, the shock location in viscous flow comes closed to the throat.