Diversos processos práticos apresentam uma seqüência de diferentes configurações de escoamentos bifásicos, com diferentes padrões de escoamento. Dependendo das razões entre as velocidades superficiais do gás e do líquido, diferentes padrões de escoamento bifásico podem ser encontrados. O padrão de escoamento slug pode ser formado a partir do padrão estratificado, devido ao crescimento das instabilidades hidrodinâmicas ou devido às ondulações no terreno onde se instala a tubulação. O presente trabalho consiste na análise do escoamento bifásico no padrão slug, ao longo de tubulações horizontais, mediante a aplicação do modelo de dois fluidos em sua forma transiente e unidimensional. Através de uma análise de estabilidade de Kelvin-Helmholtz para escoamento estratificado, estima-se a possível faixa de operação para a obtenção de um modelo matemático bem posto. Diferentes técnicas de discretização foram implementadas e testadas. Comparações com dados disponíveis na literatura foram realizadas.