Uma análise numérica de escoamento anular vertical ascendente é realizada utilizando o Modelo de Dois Fluidos unidimensional. Para escoamentos verticais, na ausência de mecanismos estabilizadores, o sistema de equações resultante é incondicionalmente mal posto. Dessa forma, visando tornar o sistema de equações bem posto, adicionou-se às equações de quantidade de movimento um termo de pressão dinâmica. Dois modelos disponíveis na literatura são investigados. O primeiro só considera a pressão dinâmica na equação de quantidade de movimento do líquido, a qual é baseada em uma simples expressão para a velocidade da onda na interface. O segundo modelo inclui uma contribuição da pressão dinâmica em ambas equações de quantidade de movimento. No presente trabalho, um terceiro modelo é proposto, o qual é baseado no primeiro modelo, com uma avaliação mais realista da velocidade da onda na interface. O efeito estabilizante da pressão dinâmica é demonstrado através de um rigoroso teste de convergência de malha. As equações de conservação são discretizadas com o método de volumes finitos, com uma integração temporal de primeira ordem e uma discretização espacial TVD de segunda ordem. Tanto o segundo quanto o terceiro modelo considerados apresentaram solução independente da malha. Parâmetros do escoamento como gradiente de pressão, espessura do filme, e variáveis características da onda obtidos numericamente são comparados com dados experimentais disponíveis na literatura, apresentando boa concordância.