Avaliação de esquemas discretizantes na solução das equações de Navier-Stokes e aplicação ao problema do escoamento na região de entrada de canais planos

Neste trabalho foi desenvolvido um algoritmo para a solução das equações de NavierStokes bidimensionais em variáveis primitivas, empregando o método de resolução semi-implícito descrito por HIRT et al.(1975). Foi utilizada malha deslocada (tipo MAC) regular, em coordenadas cartesianas. As equações de Navier-Stokes foram escritas nas formas divergente econvectiva, cujos desempenhos foram comparados no problema da cavidade quadrada com parede superior em movimento, para números de Reynolds entre 100 e 3200. Para a forma divergente foi utilizado o método dos volumes finitos e quatro esquemas de discretÍzação: central, linear a montante, QUICK, e exponencial. Para a forma convectiva foi utilizado o método das diferenças finitas e o esquema central. Os resultados foram comparados com os de GmA et alo (1982), indicando a superioridade do esquema central divergente sobre o convectivo, bem como a superioridade dos esquemas QUICK e central sobre os esquemas exponencial e linear a montante. Foi estudado o problema da região de entrada em canais planos, com contração abrupta, considerando dois casos: a região de entrada de um canal plano a partir de um reservatório de grandes dimensões, que pode ser considerado um caso de contração plana com razão de contração infinita e o caso de razão de contração finita 2: 1. Para o caso da razão de contração finita, resultados obtidos para números de Reynolds entre 10 e 426, foram comparados com os dados experimentais de DURST et aI. (1985). Os resultados obtidos mostraram que entre os quatro esquemas de discretização utilizados os que se mostraram mais confiáveis e acurados para este problema foram os esquemas QUICK e exponencial. Encontramos diferenças significativas entre os dados experimentais e os dados obtidos pelo algoritmo em tomo do plano de contração, que foram atribuídos aos erros experimentais que apontaram diferenças de até 10% nos fluxos de massa. Para o caso da razão de contração infinita foi utilizada a solução do escoamento de Jeffery-Hame1, como condição de contorno não trivial para a região de entrada. Resultados obtidos para números de Reynolds 100, 200 e 400, mostraram-se qualitativamente compatíveis com a descrição de PANTON (1984)An algorithm for the two-dimensional solution of the Navier-Stokes equations was developed in primitive variables employing the partially implicit method described by HIRT et al (1975). A regular staggered MAC mesh was used in cartesians coordinates. The Navier-Stokes equations were written in both divergent and convective forms, whose performance were compared in the problem of the moving slid square cavity for Reynolds numbers between 100 and 3200. The finite volume method was adopted for the divergent form using four discretization schemes: central differencing, upwind, QUICK and exponential. The finite difference method was adopted for the convective form using the central differencing. The results were compared with those due to GHIA et al (1982), indicating the superior performance of the divergent form central differencing over its convective counterpart, as well as the superior performance of the QUICK and the central differencing over the exponential and the upwind schemes. The problem of the entrance region in plane channels with abrupt contraction was studied considering two cases: the entrance region of the channel after a huge reservoir, which can be considered as a plane contraction of infinite ratio, and the case of a contraction with the finite ratio 2:1. For the case of contraction with results obtained with Reynolds number between 10 and 426 were compared with the experimental data due to DURST et al (1985). The results obtained indicated that the schemes QU1CK and exponential the most accurate and reliable for this problem. Significant differences with the experimental data were obtained around the plane of contraction, which were attributed to the experimental errors which have shown differences up to 10% in the mass flow rates. For the contraction with infinite ratio the Jeffery-Hamel solution was imposed as non-trivial boundary condition in the entrance. Results obtained for Reynolds numbers 100, 200 and 400 were qualitatively compatible with the description of PANTON (1984