Este trabalho objetiva prover subsídios para análise de estabilidade de sistemas eixo-mancal de deslizamento carregados dinamicamente, sob lubrificação hidrodinâmica por filme fluido. A equação da pressão na qual se baseia todo o desenvolvimento apresentado neste trabalho foi concebida por Sommerfeld. Sendo esta fundamentada na teoria de mancais infinitamente longos, assume-se que o mancal de deslizamento considerado insere-se neste contexto. As expressões para as componentes de força radial e tangencial do sistema eixo-mancal em questão podem ser obtidas para as conhecidas soluções de Half Sommerfeld e Full Sommerfeld, integrando a equação da pressão de Sommerfeld. A resolução dessa integração foi obtida através da utilização de fórmulas recursivas tabeladas por Booker, para integrais que apresentam uma mesma forma característica, denominadas Integrais de Booker, e nas quais resulta o desmembramento da integral principal. Linearizando as forças por expansão de Taylor de primeira ordem, geram-se os coeficientes das matrizes características do sistema. A implementação numérica foi feita sobre as forças e matrizes obtidas para Half Sommerfeld utilizando o Método de Newmark, e abrangendo dois casos, ordem zero e primeira ordem. Ambos são analisados e confrontados no que concerne os resultados de convergência e divergência numéricas, e efeito half-speed whirl. Conclusões e considerações finais são tecidas com respeito às características, vantagens e desvantagens, e relevância para aplicação prática. Sobre estes propósitos, observa-se que a formulação em primeira ordem mostra-se mais vantajosa por permitir utilização de intervalos de tempo maiores para descrever a trajetória assumida pelo centro do eixo. Outra contribuição deste trabalho, e não menos importante, é a obtenção de uma formulação geral para as expressões das componentes de força, bem como das matrizes de rigidez e amortecimento. Além de se poder obter as mesmas expressões resultado geradas pelas soluções de Half e Full Sommerfeld, os limites de integração são livres para assumir quaisquer valores de posicionamento angular, e não restringem o intervalo de integração a ? ou a 2?.

Um estudo experimental da interação de dois jatos de ar planos e paralelos, oriundos de bocais retangulares, foi efetuado com um sistema de Velocimetria Doppler, utilizando um laser de He-Ne. Os resultados obtidos apresentam uma boa concordância com trabalhos experimentais anteriores realizados por outros autores, isto é, os perfis de velocidade média exibiram as três regiões distintas conhecidas como região de recirculação, região de jatos duplos e região de jato único. Entretanto, tornou-se claro que alguns melhoramentos podem ser feitos no gerador de jatos, no sistema usado para o mapeamento do campo de fluxo e no Sistema de Velocimetria Doppler, cujo volume de prova deve ser feito tão pequeno quanto possível.

The present work is primarily concerned with studying the influence of an upwind spatial discretization on the capability of representing turbulent flows on aerospace applications, in the context of a flow simulation code that is fairly close to a production code. Therefore, the work addresses the issues of implementing and validating an advanced turbulence model for high Reynolds number aerospace applications in the context of an existing flux-vector splitting simulation tool, which incorporates several advances in current CFD practice, including parallel processing. The flow simulation tool used in the present work was originally developed for high speed, high altitude, hypersonic applications. Hence, the code did not include any provisions for turbulence modeling, since most flows at these conditions can be adequately treated as laminar flows. Moreover, due to the presence of strong shock waves, which are typical of hypersonic applications, a very dissipative spatial discretization scheme, based on the upwind flux vector splitting concept, was employed in the construction of the inviscid numerical fluxes. Therefore, the use of such a tool for the simulation of turbulent aerospace flows requires the implementation of a turbulence closure, as well as an adequate treatment of the excessive artificial dissipation automatically generated by the original spatial discretization scheme. In the present case, the flows of interest are simulated using the three-dimensional Reynolds-averaged Navier-Stokes equations. The turbulence closure considered is the one-equation, eddy viscosity, Spalart-Allmaras model. The work discusses in detail the theoretical and numerical formulation of the selected model, as well as the validation studies. The work also demonstrates how the spatial discretization scheme is selectively modified such that the flow simulation tool remains robust for high speed applications at the same time that it can accurately compute turbulent boundary layers. Furthermore, the work also addresses the parallelization and other high performance computational issues, demonstrating that the resultant flow simulation code can achieve adequate performance on current multi-CPU, multi-core computational clusters. Finally, the work discusses issues that could be considered for the continuation of the research effort here undertaken.

O método WENO (Esquema Ponderada Essencialmente Não-Oscilatório) é o produto evolutivo de uma sequência de estudos adicionais, que começou com o trabalho seminal de Godunov nos últimos anos da década de 50 do século passado. A partir da reconstrução clássica de Riemann de Godunov os esquemas de alta resolução passaram pelos esquemas MUSCL ("Esquema Monotonicamente Centrado para Leis de Conservação"), o ENO ("Esquema Essencialmente Não-Oscilatório"), chegando finalmente ao estágio do esquema WENO. Esquemas WENO são esquemas de ordem mais alta precisão projetados para problemas com soluções contínuas por partes contendo descontinuidades. A idéia-chave reside no nível de aproximação, quando um procedimento não-linear adaptativo é utilizado para escolher automaticamente o estêncil mais suave localmente, evitando assim que descontinuidades surjam no procedimento de interpolação como tanto quanto possível. Seguindo de perto este caminho evolutiva, uma outra característica muito importante está relacionada com a aproximação tão acurada quanto possível, dos fluxos de Riemann nas junções entre as fronteiras das células. Muitas idéias para essa aproximação apareceram ao longo do caminho. Os esquemas de Lax-Friedrichs, Roe, Steger-Warming, AUSM, AUSMPW, HLL, HLLE, estão todos em uso hoje. Um código computacional baseado em algoritmos de alta ordem e alta resolução, baseado nos algoritmos WENO e MUSCL, chamado HYNE2D, para a simulação de escoamentos viscosos em altas velocidades e com eventualmente soluções constantes por partes, foi desenvolvido pelos autores. Muitos dos esquemas mais eficientes para o cálculo do fluxo foram implementadas, e, neste trabalho, uma comparação completa entre Lax-Friedrichs, Roe, AUSM+ e AUSMPW será explorado estritamente para uma formulação em volumes finitos baseada em arestas em malhas não-estruturadas. Os problemas de validação são todos os casos teste de aferição normalmente encontrados na literatura.

Este trabalho trata da simulação numérica do escoamento incompressível, laminar e turbulento, ao longo de cavidades rasas bidimensionais. O estudo é focado, principalmente, na análise da influência dos seguintes parâmetros: (i) alongamento da cavidade; (ii) número de Reynolds e (iii) nível de turbulência do escoamento não-perturbado. A meta principal é descobrir qual a condição mais favorável no sentido de se ter o menor efeito convectivo do escoamento ao longo do piso da cavidade. As equações de transporte que modelam o problema físico são discretizadas utilizando a técnica de volumes finitos. O algoritmo SIMPLER, desenvolvido por Patankar, é utilizado para solucionar as equações discretizadas em uma malha deslocada. Os resultados obtidos demonstram que a topologia do escoamento na cavidade é diretamente afetada pela varição dos parâmetros acima considerados e que existe um regime bem definido de baixa convecção na cavidade.

O presente trabalho objetiva o desenvolvimento, implementação e análise de métodos numéricos baseados em técnicas de malhas de blocos múltiplos. As malhas de blocos múltiplos se dividem em malhas de blocos múltiplos justapostos e malhas de blocos múltiplos sobrepostos. Dentre as classes existentes destes métodos, o trabalho vai se concentrar em técnicas de malhas sobrepostos. Dentre as classes existentes destes métodos, o trabalho vai se concentrar em técnicas de malhas sobrepostas, porque estas oferecem uma maior flexibilidade para o procedimento de geração de malhas estruturadas. Estas técnicas são particularmente adequadas para configurações realisticamente complexas. O presente trabalho analisa a técnica Chimera, explicando os processos a serem feitos para a criação de uma malha adequada para aplicação desta técnica. O trabalho compreeende também o desenvolvimento e validação de códigos de simulação baseados em uma formulação bidimensional e tridimensional. O algoritmo numérico para a solução de equações Euler/Navier-Stokes usa um esquema centrado com dissipação artificial no contexto de diferenças finitas. As aplicações realizadas com esta técnica foram feitas sobre geometrias de interesse em aerodinâmica.

O objetivo principal desta proposta é a de analisar o uso do squeeze-films como um apoio flexível de um sistema eixo-mancal, com um enfoque computacional. As forças produzidas pela pressão dinâmica do lubrificante são obtidas pela solução da equação de Reynolds para o filme de fluido. Para tratar das equações que governam o movimento dos eixos rotores apoiados por squeeze-films, o método de Newmark foi empregado. Parâmetros associados à massa e velocidade do eixo e às características físicas e geométricas do mancal, tal como a viscosidade do fluido, a folga radial, comprimento e diâmetro do mancal, foram variados para permitir avaliações do comportamento orbital do eixo. Determinados grupos de parâmetros foram selecionados e, foi estudada a influência da força que produz o desequilíbrio do eixo. Resultados revelaram uma sensibilidade na estabilidade de movimento e o tamanho orbital para todos os parâmetros.

Neste trabalho são apresentados resultados para o campo hidrodinâmico e térmico para a convecção natural turbulenta em cavidade contendo material poroso. As equações microscópicas do escoamento turbulento são integradas em um volume elementar representativo para se obter equações macroscópicas válidas também no domínio poroso. Um único conjunto de equações é então discretizado, e a solução do sistema de equações algébricas obtido seguem o método SIMPLE. A intensidade da corrente convectiva através da matriz porosa é observada com o aumento do número de Rayleigh. A existência de uma fina camada limite próximas às paredes de toda a cavidade é detectada assim como a estratificação do campo de temperaturas para Ra >109.

A formulação do método de elementos de contorno é desenvolvida para aplicação em problemas de escoamento de fluidos, utilizando a técnica dos resíduos ponderados. As equações integrais de contorno e as soluções fundamentais são obtidas para o estudo de escoamentos em regime permanente e não-permanente. A aplicação do método de elementos de contorno a problemas bidimensionais de fluidos viscosos incompressíveis, em regime permanente, é estendida para fluidos viscosos compressíveis. Uma abordagem especial dos termos não-lineares das equações de Navier-Stokes é feita para permitir a aplicação do método a escoamentos com altos números de Reynolds. O método é aplicado para a solução de diversos tipos de escoamentos.

O presente trabalho teve por objetivo estudar a possibilidade de se utilizar esquemas "upwind" para simular escoamentos supersônicos e hipersônicos, sendo adotado, para isto, o esquema de separação de vetores de fluxo de Van Leer. Foram realizadas diversas simulações numéricas utilizando o escoamento interno de uma entrada de ar hipersônica e compararam-se resultados com os obtidos empregando-se o esquema centrado. O escoamento foi modelado pelas equações de Euler em duas dimensões, utilizando-se um método de volumes finitos aplicado a um contexto de malhas não-estruturadas, sendo a marcha no tempo do sistema de equações feita por meio de um método explícito de Runge-Kutta híbrido de cinco passos. A malha computacional foi obtida utilizando-se um esquema de avanço de geração com capacidade limitada de refinamento. Os resultados obtidos possibilitam a avaliação e comparação dos métodos com relação a formação da onda de choque oblíqua, bem como o comportamento em relação à pressão, à densidade e ao número de Mach. Finalmente, foram discutidas as características observadas para cada método.