Small-diameter tubes are utilized widely as expansion devices in refrigeration systems. They are employed in either kinds of short-tube orifices or long capillary tubes. Performance of these tubes is reliant upon critical flashing of the two-phase flow that controls the mass flow rate of the refrigeration system resulting in a steep reduction in pressure and temperature. The critical flow condition is approached whenever the mass flow rate increases to an amount whereby the choked-flow phenomenon occurs at the outlet of the tube. Due to their very small tube diameter, the evaporating two-phase flow, and the choked-flow condition, numerical analysis of flow through short-tube orifices is challenging. Accordingly, all available numerical analyses of such flows are performed as one-dimensional and in the majority of them, auxiliary correlations are applied to simplify the solution procedure. Typical approaches include homogeneous flow models and separated flow models, both of which consider the two-phase region in thermal equilibrium. The most comprehensive method for analyzing such flows is the two-fluid model in which there is no assumption of equilibrium between phases. Because of the complicated nature of this model, it has been used in a very limited number of previous investigations. Furthermore, two-phase flow calculations at the entrance and vena contracta region were eliminated. In the current investigation, additional steps utilized to improve the accuracy of computations include the following: (1) applying the most comprehensive two-fluid model including the effect of various two-phase flow patterns and the metastability of liquid phase, and (2) performing a two-phase analysis of the evaporating flow through the entrance and vena contracta regions which involves simulating the region as a converging diverging tube and performing a quasi-one-dimensional solution of governing equations through this region. Results showed more compatibility with experimental data in comparison with those of previous investigations for predicting the critical flow condition of common refrigerants HFC-134a and HFC-410a through short-tube orifices and long capillary tubes.