O escoamento anular é caracterizado por um núcleo gasoso fluindo a alta velocidade com um filme líquido no seu entorno, molhando a parede do duto. A presença de gotículas líquidas no núcleo gasoso resulta em impacto relevante em características do escoamento anular, como gradiente de pressão e propriedades das ondas presentes no filme líquido. A formação de gotículas se dá usualmente na crista das ondas de perturbação presentes na interface líquido-gás. No presente trabalho, é realizado um estudo do regime anular com presença de gotículas em tubulações verticais utilizando o Modelo de Dois Fluidos unidimensional. Um modelo de transferência de massa das gotículas é desenvolvido e acoplado ao modelo de Dois Fluidos. O modelo resultante permite capturar a evolução automática da interface gás-líquido e a formação de ondas de filme líquido e sua influência no desprendimento e deposição de gotículas. Analisa-se o desempenho de três modelos de entranhamento de gotículas disponíveis na literatura, além de um modelo de deposição de gotículas. Considerando que gotículas são criadas por cisalhamento nas cristas das ondas de perturbação, modificações dos modelos são propostas com a finalidade de melhor capturar a influência das ondas do filme líquido nos processos de entranhamento e deposição de gotículas. Parâmetros do escoamento como gradiente de pressão, espessura do filme do líquido e variáveis relacionadas com as ondas interfaciais são avaliados, mostrando boa concordância com dados experimentais disponíveis na literatura.