Two-phase gas-solid particles nonstationary axisymmetric flows through converging-diverging nozzles and the free overexpanded jets emerging from the nozzles are studied numerically. The solution of the flowfield is carried out until a steady flow is established. The Eulerian approach is used to describe the flowfield and both phases are treated in the homogeneous mixtures approximation. A monotone second order accurate in space and time W-modification of Godunov's scheme is applied for the numerical solution of the governing equations. The evolution of the discontinuities is investigated as well as the effects of the particle size and the loading ratio on the flow pattern. The two-phase flowfield is compared with a similar pure gas flowfield in a two-dimensional nozzle and in the plume. Limiting cases of the two-phase flow, frozen and equilibrium, are considered. Intermediate regimes of the two-phase flow are compared with the frozen and the equilibrium flows.