Abstract: A numerical method for calculating a fully-coupled model of one-dimensional two-phase nozzle flow is described in brief. Governing equations are derived. The Runge-Kutta-Gill method is adopted to solve the equations numerically. The solution for uniform acceleration of gas phase is utilized as an approximate solution for the dimensionless distance x*(=x/L) = 0 to 0.02,i.e.ug = ax*. The initial value of a is determined by a trial and error approach and then checked at the throat. The solution is continued to the throat where the conditions must be satisfied. To eliminate the singularity of the equation, a superseding variable is introduced. The computations are performed for velocity-lag as well as temperature-lag of five parsticle sizes and two ratios of particle to gas flow rates.Temperature, pressure and velocity distributions of both gas and particle along the length of the nozzle are also obtained.In addition, various modes of one-dimensional two-phase nozzle flow have also been briefly reviewed.