Um algoritmo implícito eficiente de fatorização aproximada, concebido para resolver a equações conservativa do potencial completo, foi empregado para cálculo do escoamento transônico em torno de configurações do tipo asa-fuselagem. O algoritmo usa ponderadamente as características da densidade do escoamento à montante ("upwind bias"), também conhecido como esquema de densidade artificial, para estabilidade em regiões de escoamento supersônico. Um gerador de malha algébrico foi desenvolvido para fornecer malhas do tipo O-H para o código aerodinâmico. O gerador permite controle de suavização dos pontos interiores da malha e ajuste do posicionamento da esteira. Uma transformação numérica é usada para estabelecer uma malha de diferenças-finitas que se conforma ao corpo. Resultados de escoamentos para asas montadas entre paredes e com uma das fronteiras laterais em escoamento uniforme são comparados com outros resultados independentes.

No presente trabalho, são apresentados os campos de velocidade na esteira próxima de cunhas. Em uma primeira série de experimentos, utiliza-se um único modelo de cunha e verifica-se o efeito da variação do número de Reynolds (Re) entre 0,6x103 e 3,0x103. Na segunda série de experimentos, o número de Reynolds é fixado e o ângulo de separação (semi-ápice) da cunha é variado 15, 30e 45. Na terceira série, para o mesmo Re foram ensaiados dois modelos de cunha com dimensões diferentes da base (semi-ápice de 45) para obter a influência da razão de bloqueio (9,8% e 19,7%). Um aspecto importante deste trabalho foi o desenvolvimento da metodologia experimental para a obtenção de resultados confiáveis, usando um anemômetro laser. Além disto, foi possível realizar medidas em valores baixos de número de Reynolds, em uma faixa onde resultados experimentais não são disponíveis. A freqüência de liberação de vórtices foi determinada experimentalmente (cálculo do número de Strouhal) com o auxilio do LDA. A partir dos campos de velocidade média (u e v) e flutuante (uRMS e vRMS) é possível determinar a posição do núcleo do vórtice médio e a posição onde ocorrem os máximos de uRMS e vRMS. Estas posições, permitem a ampliação da discussão sobre a posição do final da região de formação. Com os campos de velocidade foi determinada a velocidade no ponto de separação e, indiretamente, o coeficiente de pressão na base.

Este trabalho investiga a eficiência do método multigrid quando usado na obtenção de campos de velocidade e temperatura bidimensionais laminares em domínios retangulares. A análise numérica é baseada no esquema de discretização em volumes finitos aplicado a malhas estruturadas ortogonais regulares. O desempenho do algoritmo multigrid correction storage (CS) é analisado para diferentes números de Reynolds na entrada (Rein) e para um número distinto de malhas. Até quatro malhas computacionais foram usadas para ambos os ciclos V e W. Soluções para os campos de temperatura e velocidade obtidos simultaneamente e não-simultaneamente foram investigados. As vantagens do uso de mais de uma malha computacional são discutidas. Para as soluções simultâneas os resultados indicam um aumento do esforço computacional com o incremento do número de Reynold Rein. O número ótimo de relaxações intermadiárias para os ciclos V e W é discutido.

Neste trabalho é apresentada uma análise numérica para predição de escoamentos turbulentos incompressíveis em geometrias complexas. Modelos de turbulência k - e não lineares que representam um avanço sobre modelos de viscosidade turbulenta clássicos (linear), foram empregados em um sistema de coordenadas generalizadas na solução de diversos tipos de escoamentos bidimensionais e com recirculação. O método numérico utilizado para a discretização das equações foi o método de volumes finitos com arranjo co-localizado. O algoritmo SIMPLE foi usado para relaxação do sistema de equações algébricas obtido. A resolução do sistema linear de equações resultante foi feita através do método SIP. A função de parede clássica e um modelo de Baixo Reynolds foram utilizados para descrever o escoamento próximo à parede. Aqui foi proposto um tratamento numérico inovador para os termos difusivos não lineares das equações de momentum a fim de se evitar dificuldades de convergência em alguns casos específicos. O tratamento consiste em dividir os fluxos difusivos não lineares em um termo implícito, na matriz dos coeficientes e outro explícito agrupando ao termo fonte. Na literatura, não há informação específica sobre o tratamento numérico destes termos não lineares. Resultados para escoamentos turbulentos bidimensionais foram obtidos usando um total de quatro diferentes modelos. As formas linear e não linear do modelo k - e usando as abordagens de Alto e Baixo Reynolds foram empregadas. Assim, as formulações de Alto Reynolds foram designadas respectivamente por L_HRN e NL_HRN, e os mesmos modelos aplicando a formulação Baixo Reynolds foram chamados de L_LRN e NL_LRN. De uma forma geral, o modelo não linear usando a abordagem Baixo Reynolds, NL_LRN, apresentou melhor desempenho que os demais. Também foi analisado o desempenho dos modelos não lineares (NL_HRN e NL_LRN) na predição de escoamento secundário que ocorre em dutos de seção não circular. Nestes casos, modelos isotrópicos são incapazes de descrever o escoamento secundário em tais canais.