O tema central desta tese é a proposição de modelos de fechamento para a densidade de força interfacial necessária à solução do modelo de dois fluidos empregado amplamente na descrição de escoamentos multifásicos. os modelos de fechamento empregados até então fazem uso da decomposição dos diversos fenômenos que ocorrem na interação das fases como o arrasto, tração, desenvolvimento de camadas limites hidrodinâmicas, gradientes adversos de pressão e vorticidade, possuindo cada um sua própria descrição matemática e correlações para seus respectivos parâmetros empíricos. estes últimos são dependentes de uma série de números adimensionais, entre eles o número de reynolds, eötvos e morton, sendo estes os principais. as correlações geralmente são dependentes do regime e do grau de interação das fases, bem como da sua dispersão ou não no seio fluido do escoamento. esta descrição particular e as diversas correlações criam novas formas de se resolver problemas até então mal resolvidos e/ou mesmo sem solução; contudo, elas não são genéricas e sua aplicabilidade fora dos limites em que foram calibradas deve ser bem avaliada e nem sempre os resultados após um longo tempo de tentativas e erros são satisfatórios. o presente trabalho trata de forma diferenciada o problema de fechamento da densidade de força interfacial modelando como o tensor tensão atua sobre a interface de uma forma global, sem separar os efeitos categorizando-os nos tipos supracitados, e, ao invés de um somatório de efeitos e o uso de uma série de correlações, deve-se apenas responder à questão de como descrever matematicamente o tensor sobre a interface, e deixar ao cargo do próprio escoamento “decidir” que forças atuam sobre a interface dos fluidos. a metodologia aqui proposta é descrita em seus pormenores físicos e matemáticos e alguns casos de teste, “benchmarks” numéricos, são escolhidos para o seu emprego. o modelo foi implementado no simulador comercial ansys cfx. os resultados obtidos demonstraram a eficácia da abordagem em problemas onde a interface é bem definida, como os de superfície livre, mesmo com mudança local de morfologia, sem, todavia, apresentar o mesmo desempenho em escoamentos dispersos e com baixas frações volumétricas. não obstante, a atuação conjunta de ambas as abordagens (a clássica e a proposta) neste último caso, mostrou-se muito eficiente na captura do comportamento físico de um problema complexo como uma coluna de bolhas.