Esta dissertação apresenta a simulação numérica do escoamento tridimensional num canal contendo um feixe de varetas, sustentado por grades de apoio e com desbalanceamentos de 33%/67% e 0%/100% nas vazões na entrada. Os resultados computacionais foram obtidos com o código PHOENICS 2.0. A contribuição deste trabalho foi propor um modelo matemático para resolver as equações de conservação usando o método das porosidades. Com isso foi possível tratar um problema que envolveria a geração de uma malha computacional complexa de uma forma mais simples: usando coordenadas retangulares com menor número de pontos e viabilizando o uso dos microcomputadores disponíveis. Para os desbalanceamentos de vazões 33%/67% na entrada do canal os resultados obtidos foram concordantes com os dados experimentais [4]. Entretanto, para o desbalanceamento de vazões com uma secção de entrada totalmente bloqueada, (0%/100%) verificou-se que todos os valores calculados para as velocidades axiais diferem dos dados experimentais de no máximo 10%, exceto no ponto de coordenada x/a=0.86 no qual esta diferença é no máximo 30%, sugerindo com isso que outras correlações empíricas de perdas de carga em varetas devem ser testadas.

The COBRA series computer codes have been widely used in simulating the flow in sub channels in typical pressurized water reactor. The results of such simulation have been validated against a large number of experimental results. However, in the case of unbalanced flow in parallel, short length, sub channels the agrément between calculated and experimental [4] values were not very good. The presente work shows numerical simulation of the three dimensional flow in a typical rod bundle with spacer grids and inlet sub channel flow unbalance of 33% / 67% and / 100%. The PHOENICS 2.0 computer code was used in the flow simulations. The contribution of this work was to propose a mathematical model to solve the conservation equations using the porosities method which allows to treat the complex grid computational problem using a simple rectangular coordinates system, in order to run the simulations with a available micro computer. For a unbalancead flows with 33% / 67% the calculations provided a very good agrément between calculated and experimental results [4]. However for more severe unbalanced flows with a entrance section entirely blocked (0% / 100%), it is verified that all the calculated values to axial velocities are different from experimental measures in a maximum 10%, except in point X / A = 0.86 which this difference is a maximum 30%; it is meaning that other empirical correlation of pressure drop in rod should be tested.