---

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.

---

Analysis of air bubbles entrainment in liquid slug body is one of the most important and complicated phenomena during slug flow regime. In the present attempt, a method is proposed for slug modeling to consider the air bubble entrainment into slug liquid body. The effect of consequencies and their impact on slug behaviour to predict more accurate correlations for slug parameters are estimated and calculated. This method considers a two-fluid single pressure model, combined with population balance model for equal bubble diameter series and solves using volume of fluid. In this regard based on slug and hydraulic jump similarity, a correlation for air bubble entrainment rate and its mechanism selected. This correlation developed in the form of a user defined function code and coupled with other models in FLUENT solver to calculate slug flow. Finally the result of this numerical modeling is validated with the result of other numerical and experimental results where exist in the related literature. The result is consist of the slug flow profile, entrained air bubble profile and their diameter distribution, vortex at slug front, pressure distribution during slugging, slug mixture velocity, turbullent model parameters and etc.

---

Slug flow is one the most complicated flow regime in industrial process that is seen for wide range of fluid flow. However, there are always some differences between experimental and numerical results obtained for slug flow. Proceeding the perivious attempt on the selection of the best turbulent model for numerical simulation, the slug flow is simulated numerically in two dimension by applying implicit VOF method and k-ε RNG turbulent model. To extract the slug flow parameters accurately FLUENT solver is used. The differences of the obtained results with and without turbulent model during slugging is also presented. To overcome this complicated flow behavior, a new user defined function code is developed. This UDF computes and predict slug parameters from FLUENT solver result without increasing the computational cost. The important slug parametrs are presented which are: liquid slug body velocity, liquid film velocity, slug front and tail velocity, slug center position and length, slug front and tail positions, pressure difference across slug, wall shear stress, slug mixture velocity, slug initiation time and position from the duct inlet. These parameters are analyzed and discused in details after that they had been validated