Este trabalho apresenta uma comparação numérica de formulações estabilizadas Streamline-Upwind/Petrov-Galerkin (SUPG) com operador de captura de descontinuidades aplicadas a problemas de escoamentos bi-fásicos de fluidos imiscíveis ligeiramente compressíveis em meios porosos, onde os efeitos da gravidade e pressão capilar são considerados.
O conjunto acoplado de equações diferenciais parciais resultante é aproximado no espaço pelo método dos elementos finitos. A integração no tempo do sistema semi-discreto é feita através de um algoritmo bloco-iterativo preditor-multicorretor, onde os sistemas lineares resultantes são resolvidos através de técnicas iterativas elemento-por-elemento. As diferentes formulações são avaliadas do ponto de vista de precisão e desempenho computacional.