An important phenomenon in the gas–liquid two-phase mixtures is a sudden change in the flow that may lead to transition of the flow regime from non-homogeneous slip flow to homogenous flow. This phenomenon is called mixing shock and has been investigated by several researchers. In the present paper, a more comprehensive model is proposed by including the entrainment ratio in the governing equations. Moreover, parametric studies are performed to assess the importance of this parameter in different conditions. The results of the present study indicate that at Euler numbers less than four the effect of gas mass flow is insignificant. However, at higher Euler numbers the differences of the present and previous models are noticeable. Since two possible solutions for the aftershock state exist, four criteria—entropy change across the shock, possibility of an expansion shock, positivity of the Euler number, and choking flow condition—are considered to identify the correct solution. The results indicate that the flow after the shock could only be of a subsonic type, and the mixing shock is compressive. A comparison of the possible realizable zones for the developed model and the previous models indicates that the developed model predicts a larger area for the realizable solution of the mixing shock. A comparison between analytical results and experimental data shows that the developed model predicts reasonable results.