Experiments on high-pressure vessel decompression have shown that vaporization occurs in 'boiling shocks' moving with a velocity of ∼10 m s −1 . To explain this phenomenon, a model accounting for bubble breakup was suggested (Ivashnyov, Ivashneva & Smirnov, J. Fluid. Mech., vol. 413, 2000, pp. 149-180). It was shown that the explosive boiling was caused by chain bubble fragmentation, which led to a sharp increase in the interface area and instantaneous transformation of the mixture into an equilibrium state. In the present study, this model is used to simulate nozzle flows with no change in the free parameters chosen earlier for modelling a tube decompression. It is shown that an advanced model ensures the best correspondence to experiments for flashing flows in comparison with an equilibrium model and with a model of boiling at a constant number of centres. It is also shown that the formation of a boiling shock in a critical nozzle flow leads to autovibrations.