Métodos integrais de contorno aplicados à solução de escoamentos aerodinâmicos permanentes e transientes

Abstract

Nessa dissertação é apresentado o método integral de contorno para escoamentos potenciais bidimensionais e tridimensionais incompressíveis. Aplicando o teorema de Green para as regiões internas e externas a um corpo arbitrário e combinando os resultados obtém-se a equação integral de contorno do método em questão. A superfície do corpo é representada por elementos de malha triangulares e quadrilaterais planos com distribuição constante de singularidades (dipolos e sorvedouros) em cada elemento. A condição de Kutta é aplicada de forma simples para o caso de escoamento permanente: as intensidades dos filamentos de vórtice da esteira são determinadas pelas diferenças de intensidade de dipolos nos elementos de malha adjacentes ao bordo de fuga. O método é de baixo custo computacional, mesmo quando usado para simular escoamentos em torno de corpos tridimensionais com geometrias complexas. Os resultados numéricos do método integral de contorno são comparados com resultados analíticos, resultados experimentais, resultados da teoria da linha de sustentação e do método vortex-lattice. Métodos integrais de camada limite delgada bidimensional foram utilizados na determinação dos principais parâmetros de camada limite. Foram simulados a camada limite laminar, um critério de transição e a camada limite turbulenta. Por fim, foi implementada a versão bidimensional do método integral de contorno para escoamentos potenciais transientes. O método em questão considera o aerofólio movendo-se no fluido e trata o problema em termos de um sistema de referência não-inercial fixo no corpo. A condição de Kutta é implementada considerando-se que a pressão é igual nas superfícies superior e inferior do bordo de fuga. A formulação matemática do escoamento potencial transiente explorado neste trabalho é baseada no princípio de Benoulli da conservação da energia mecânica de um escoamento na sua versão não-estacionária. O contorno do aerofólio é substituído por elementos de vórtices e de fontes planos. Em particular, um elemento adicional é incorporado ao bordo de fuga com o objetivo de se representar o início da formação da esteira atrás do corpo. O restante da esteira de vorticidade é representada por pontos de vórtice. O sistema de equações resultantes é não-linear e é resolvido por um método de eliminação de Gauss iterativo. _________________________________________________________________________________ ABSTRACT

In this dissertation we describe a boundary integral method for calculating the incom- pressible potential °ow around arbitrary, lifting, two-dimensional and three-dimensional bodies. By using Green theorems to the inner and outer regions of the body and com- bining the resulting expressions we obtain a general integral representation of the °ow. The body surface is divided into small quadrilateral and triangular elements and each element has a constant singularities distribution of sinks and dipoles. The applica- tion of Kutta's condition for the steady °ow is quite simple; no extra equation or trailing-edge velocity point extrapolation are required. The method is robust with a low computational cost even when it is extended to solve complex three-dimensional body geometries. The boundary integral code developed here is veri¯ed by comparing the numerical predictions with experimental measurements, analytical solutions and results of the lifting-line theory and vortex-lattice method. A two-dimensional integral boundary layer method, was incorporated into the two-dimensional boundary integral method and is intended to give a preliminary estimate of boundary layer properties over a con¯guration. The method comprises a laminar boundary layer analysis, a transition and a turbulent boundary analysis. Finally, an unsteady two-dimensional boundary integral method has been developed. The present approach treats the body as moving through the air and put the mathematical problem in terms of a moving coordinate system attached to it. The Kutta condition is applied by imposing the same pressure on the upper and lower surfaces at the trailing edge. For this end, the unsteady Bernoulli equation is used. The airfoil contour is replaced by straight-line elements and a small straight-line wake element is attached to the trailing edge. The circulation on this element is imposed in the way that the global circulation around the airfoil and the wake is constant. A downstream wake of concentrated vortices is formed from the vorticity shed at earlier times and is convected according to the local velocities. The basic set of equations is nonlinear and an iterative procedure is adopted for its solution.