Neste trabalho foi desenvolvido um programa para simulação numérica do fluxo bifásico em meios porosos com transferência de massa entre as duas fases e transporte dos solutos em cada fase. O programa calcula o aumento da mobilidade da fase não aquosa e da solubilidade de seus componentes na água como resultado da adição de um co- solvente. O objetivo foi simular problemas de contaminação no subsolo com líquidos imiscíveis com a água e a remediação através da técnica de injeção de co-solventes, assim como a contaminação a partir de vazamentos de gasolina acrescida de álcool em postos de abastecimento. Podem-se ainda resolver problemas de fluxo bifásico envolvendo uma fase gasosa. Foi desenvolvido e implementado um modelo para a transferência de massa durante o fluxo das duas fases. O cálculo das propriedades das duas fases ao longo do tempo é feito utilizando o diagrama ternário de fases. Para resolver as equações do fluxo e do transporte de solutos foi empregado o método de elementos finitos.