A simulação numérica de escoamentos de fluidos em uma grande variedade de aplicações requer a utilização de técnicas numéricas de alta eficência e recursos computacionais de alto desempenho. O objetivo deste trabalho é iniciar uma investigação de escoamentos de fluido durante o enchimento de compartimentos de reservatórios. Uma abordagem inicial foi tratar problemas de escoamento em um canal, rebuscando a geometria do domínio para contemplar problemas mais complexos. Este trabalho apresenta o desenvolvimento e os resultados obtidos de um método numérico para simulação de escoamento de fluido incompressível em um domínio tridimensional, onde as equações de Navier-Stokes são desenvolvidas em uma formulação euleriana e discretizadas pelo método de elementos finitos. Os termos convectivos destas equações foram tratados pelo método semi-lagrangeano e o método de Galerkin foi utilizado para discretização espacial, um método baseado em decomposição LU foi utilizado para desacoplar as componentes de velocidade e pressão, sendo esta última calculada utilizando-se uma aproximação hidrostática. O domínio tridimensional foi representado por uma malha manipulada por uma estrutura de dados topológica, formada por células que definem elementos prismáticos lineares. Foram realizados experimentos sob várias alterações na geometria do domínio e também sob diferentes condições iniciais. Os resultados mostraram uma boa aproximação do método, quando analisado comparativamente a uma solução analítica

The numerical simulation of fluid flow over many applications require the use of numerical techniques of high efficiency and demand high computational power. This work aims at initiating an investigation about fluid flows while filling reservoirs. The initial approach was to deal with fluid flows in a retangular duct, as increasing the complexity of its geometry in order to model more complex cases. This document describes the development of a numerical method for the simulation of incompressible fluid flow over a threedimentional domain, where the Navier-Stokes equations were written under an Eulerian formulation and discretized by the Finite Elements Method. The semi-Lagrangean method was used to discretize the convective terms and the components of velocity and pressure were decoupled through the use of a method based on LU decomposition, where the final pressure was determined by using a hidrostatic aproximation. The threedimentional domain was represented by a mesh, manipulated by a topologic data structure, formed with cells that define linear prismatic elements. Many experiments were performed under different geometries of the domain and also under different initial conditions. The result showed a good approximation of the described method, when compared with an analitical solution