Pipeline and equipment multiphase flow modeling in oil industry has become a crucial task due to the constant increasing fossil fuel demand and latest oil discovers in Brazil. Among all possible flow configurations that multiphase flow can take, the slug flow regime plays an important role due to its intermittent characteristics and its strong pressure gradients in pipelines, where this flow regime can be found. In the present work, a numerical study of a slug flow is performed, with the finite volume method. The flow and pressure field are determined, as well as, the Taylor bubble nose shape and the translational velocity of a three dimensional single unit slug in a horizontal pipeline. The flow field was considered isothermal and it was employed a coordinate system coincident with the Taylor bubble. To model the turbulent flow, the RNG K−E model was selected. The interface was determined with the Volume of Fluid models (VOF). The results from the simulations were compared with available experimental data, presenting good agreement.