Many models have been applied to slug flow using laminar flow condition. The results obtained from these models are not consistent with the physical behavior of slug flow. Furthermore, discussion on the turbulent models is very rare or not related to the such flow regime. The slug regime is a complicated regime with shear flow and high strain in addition to some vorticity at some sections of the flow. In the present attempt at first stage, the turbulent models differences, the initial assumptions to drive, privilages and shorcomings have been considered with details. Then, its consistency with the physics of slug flow was analysed with high accuracy. In the second stage, simulations using different turbulent models were conducted. The obtained results were compared to each other and with the experimental results of other investigators. Finally, the most consistent model with the physics of the slug flow was selected. The turbulent model of RNG k-ε showed more reliable in compare to other turbulent models. Thus, it was selected and used to obtain slug flow behavior with higher accuracy. The parameters as pressure distribution during slugging, slug mixture velocity, slug initiation time and position from the duct inlet with RNG model conducted and presented with detailed explanations.