O constante crescimento na demanda de combustíveis fósseis e as últimas descobertas de petróleo no Brasil têm tornado a modelagem de escoamentos multifásicos em tubulações e equipamentos na indústria do Petróleo numa tarefa crucial. Dentro as configurações possíveis que o escoamento multifásico pode adotar, o regime das golfadas possui uma enorme relevância devido às suas características intermitentes e seus fortes gradientes de pressão produzidos nas tubulações por onde escoa. Neste trabalho, estuda-se numericamente, com o método de volumes finitos, o regime das golfadas, determinando-se o campo de velocidade e pressão, assim como a forma do nariz da bolha de Taylor e a velocidade de translação de uma unidade básica de golfada tridimensional em uma tubulação horizontal. Considerou-se o escoamento isotérmico, e utilizouse o referencial coincidente com a bolha de Taylor. Para modelar o escoamento turbulento utilizou-se o modelo K−E RNG. A interface foi determinada com o modelo VOF (Volume of Fluid). Os resultados das simulações foram comparados com dados experimentais disponíveis, apresentando uma boa concordância.