O objetivo deste trabalho é ser uma contribuição ao entendimento dos mecanismos envolvidos na migração por gravidade dos compostos orgânicos chamados de DNAPLs, quando eles são liberados em meios porosos e em meios porosos fraturados, para aportar ao desenvolvimento de tecnologias efetivas orientadas principalmente à localização e à remediação do sistema subterrâneo contaminado. Primeiramente são apresentados os conceitos elementares envolvidos nos modelos matemáticos que descrevem o fluxo bifásico em meios porosos, o processo de modelagem de um problema geral da natureza, os modelos conceituais, os matemáticos e os numéricos, e a aplicabilidade dos modelos conceituais conforme a considerações de escala. Em segundo lugar, são desenvolvidas as equações matemáticas que governam os fenômenos em estudo e são apresentadas as soluções às equações governantes a partir de técnicas computacionais e esquemas de integração numérica. As equações do fluxo bifásico são resolvidas mediante técnicas de elementos finitos mistos hibridizados (EFHM) e elementos finitos descontínuos (GD), e as equações do transporte de contaminantes são resolvidas mediante a técnica dos elementos finitos convencionais (EF). Seguidamente são avaliados numericamente problemas de transporte de contaminantes em 1D e 2D, problemas de transporte de contaminantes com transferência de massa, problemas de fluxo bifásico em 2D, e problemas acoplados envolvendo tanto fluxo bifásico como transporte de contaminantes com transferência de massa. Finalmente, são apresentadas as conclusões do trabalho desenvolvido bem como sugestões para trabalhos futuros.