A model of two-component two-phase critical flow is presented. The modelling approach is based on one-dimensional homogeneous gas-liquid two-phase isentropic flow of mixture. The homogeneous model is modified by taking into account the void fraction and two-phase mixture density dependence on velocity slip. The velocity slip is calculated using Chisholm correlation that depends on the gas phase quality and Zivi correlation for the prediction of the maximum velocity slip values. At the location of the critical flow the two-phase mixture velocity equals sonic velocity and it is calculated with the so-called ''frozen sonic velocity' model. The model is validated against data measured in air water flow at the PUMA experimental facility. Obtained results are presented together with the predictions by the well-known Fauske model. It is shown that Fauske model overpredicts measured critical mass fluxes, while the present model shows acceptable agreement with the measured data.