O presente trabalho tem por objetivo a avaliação das características aeroelásticas referentes a flutter de uma aeronave militar de treinamento em sua fase inicial de projeto. Para a execução da análise ée elaborado um modelo aeroelástico da aeronave, composto de dois modelos distintos: um em elementos finitos representativo das características dinâmicas da aeronave e um modelo aerodinâmico não-estacionário baseado na teoria Doublet Lattice. Correções para ajuste do comportamento aerodinâmico estacionário do modelo são feitas com base em resultados obtidos junto ao grupo de aerodinâmica do Projeto Treinador Avançado. A ocorrência do fenômeno de flutter colocará restrições ao projeto. Os detalhes do modelamento dinâmico e aerodinâmico, bem como ajuste aerodinâmico adotado, são também apresentados. O software utilizado na solução dos problemas dinâmico e aeroelástico ée o MSC.Nastran V.70.7. Um estudo paramétrico ée realizado para analisar o comportamento do modelo aeroelástico após alteração de parâmetros relacionados com as superfícies de controle (canard, aileron, leme e empenagem horizontal). Os resultados são apresentados em forma gráfica através dos gráficos V-g-f que mostram a tendência do comportamento do amortecimento e freqüência de cada modo com a variação da velocidade do escoamento.

O escopo deste trabalho ée avaliar o efeito da cablagem na modelagem da distribuição de temperatura em equipamentos eletrônicos em racks de aeronaves através do uso de CFD (Computacional Fluid Dynamic). Para alcançar este objetivo, o estudo se concentra na simulação de quatro casos. Os casos padrão e sem ventilação que não consideram o efeito da cablagem. Os casos de velocidade constante e de pressão constante na entrada empregam o efeito do meio poroso. A cablagem ée modelada como um meio poroso, ou seja, um termo fonte ée adicionado nas equações da continuidade, de quantidade de movimento e de energia. O modelo de turbulência escolhido para a simulação ée o k- padrão por ser utilizado freqüentemente em escoamentos turbulentos semelhantes ao caso em estudo, além de reduzir esforços computacionais se comparado ao modelo das tensões de Reynolds ou o LES (Large Eddy Simulation). Uma malha de 93503 elementos ée gerada para a simulação. A partir das temperaturas experimentais, os fluxos de calor de cada equipamento são calculados e empregados como parâmetros de entrada para cada caso. O caso padrão e sem cablagem mostra o escoamento e a distribuição de temperaturas preliminares. O caso sem ventilação e sem cablagem, isto é, em uma condição sem entrada de ar condicionado indica que o efeito da convecção natural ée importante na distribuição de temperatura dos equipamentos. Dentro do caso de velocidade constante, o efeito do meio poroso ée verificado através do cálculo de diversos termos fonte, sendo que o aumento da pressão não afeta a distribuição de temperatura no interior do rack. No último caso, ou seja, de pressão constante com cablagem, diversos termos fontes são empregados e verifica-se que o aumento de velocidade na entrada modifica a distribuição de temperatura nos componentes. A existência de uma velocidade ótima foi notada, ou seja, um acréscimo na velocidade do escoamento na entrada do rack não altera mais a distribuição de temperatura dos equipamentos eletrônicos.

Simulações numéricas foram realizadas utilizando-se um código computacional desenvolvido para resolver o sistema de equações de Navier-Stokes com média de Reynolds que modela o escoamento compressível turbulento em torno de um aerofólio NACA 0012. Foram utilizadas malhas estruturadas tipo O geradas algebricamente e diversos refinamentos puderam ser feitos. O método de volumes finitos foi empregado para discretizar o sistema de equações diferenciais parciais e os esquemas explícitos de MacCormack e Jameson foram implementados. Termos de viscosidade artificial foram adicionados explicitamente através de um modelo não-linear. O modelo de turbulência de uma equação de Spalart e Allmaras foi implementado para resolver o problema de fechamento da turbulência. Inicialmente, a formulação de Euler foi usada e resultados para a distribuição de pressão e coeficientes aerodinâmicos foram obtidos para quatro casos de escoamentos transônicos não-viscosos sobre o aerofólio. As soluções foram comparadas com os resultados de outros métodos numéricos disponíveis na literatura. Em seguida, um dos casos foi utilizado para avaliar a influência dos parâmetros numéricos como a viscosidade artificial e o refinamento da malha. Outro caso foi utilizado para comparar os esquemas explícitos de MacCormack e Jameson. Por último, o modelo de turbulência de uma equação de Spalart e Allmaras foi utilizado para a formulação de Navier-Stokes e as soluções foram comparadas com os dados experimentais de Harris e outros resultados numéricos obtidos com o modelo de turbulência algébrico de Baldwin e Lomax.

Os sistemas de proteção contra gelo são comumente utilizados na aviação. O piccolo tubo, que é um dispositivo com vários furos por onde passam os jatos de ar quente que aquecem as superfícies a serem protegidas, é o mais utilizado. Dependendo do processo de manufatura, os furos podem variar desde um formato cilíndrico até um formato cônico. O objetivo deste trabalho foi analisar a influência da conicidade destes furos na transferência de calor com o lábio da entrada de ar. Os resultados mostraram que o impacto é mínimo. Para tal análise, foram utilizados tanto programas de geração de malha (ICEM), quanto de mecânica de fluidos computacional - (FLUENT). O modelo térmico foi validado com a vasta literatura disponível, que inclui tanto dados experimentais, quanto simulações que também utilizaram mecânica dos fluidos computacional.

In this current work it is presented the two-dimensional axisymmetrical air flow simulations under laminar regime in a conical diffuser using Computer Fluid Dynamics (CFD). The use of diffusers in laminar flow can be seen in micro-pumps applications, especially for micro-electronics cooling. The objective is to analyze the static pressure recovery coefficient (Cp) for Reynolds 64, varying the diffuser expansion angle and the diffuser exit/entrance area ratio (A2/A1 = 1.5 and 2.0) for the cases with and without tail pipe. The diffuser geometric configurations and the Cp formulation are based in the current ESDU 73024 (Engineering Sciences Data Unit) publication. The partial differential equations system (Continuity and Navier-Stoke) was solved using a computer program based in the numerical Finite Element method. For diffusers without tail pipe (A2/A1 = 1.5 and 2.0), the Cp values under turbulent flow are higher than Cp values under laminar flow for the same expansion angles. For diffuser with tail pipe (A2/A1 = 1.5 and 2.0), the Cp values under turbulent flow are higher than Cp values under laminar flow up to 18 diffuser expansion angle. Above 18, the Cp for turbulent and laminar flow follow a similar trend. For diffuser (A2/A1 = 1.5 and 2.0) under turbulent flow with and without tail pipe, the Cp results were similar up to 10 diffuser expansion angle. Above 10, the diffusers with tail pipe presented Cp results higher than diffusers without tail pipe. The same occurs for diffuser (A2/A1 = 1.5 and 2.0) under laminar flow with and without tail pipe. Therefore, the finite element method showed a good agreement to solve this kind of problem and the results are important once static pressure recovery coefficient data for laminar flow is scarce in the literature.

Simulações são realizadas utilizando a dinâmica de fluidos computacional (DFC) para prever o desempenho e características de escoamento em uma bomba ejetora anular. A simulação numérica do escoamento no interior da bomba ejetora anular é realizada utilizando o código FLUENT 6.3.26. No presente estudo, o método de volumes finitos é utilizado para resolver um escoamento tridimensional, em regime permanente e incompressível modelado pelas seguintes equações matemáticas: continuidade, momentum e modelo de turbulência . Uma boa concordância foi obtida na comparação entre os resultados numéricos e experimentais disponíveis. Então, os resultados numéricos são utilizados para analisar as seguintes dependências: razão de altura manométrica versus razão de vazão volumétrica e a eficiência global versus razão de vazão volumétrica. Os perfis de velocidade e pressão estática também são usados para estudar o processo de mistura entre as duas correntes na câmara de mistura da bomba ejetora. A boa concordância com os resultados experimentais e a compreensão do escoamento de água obtidos nesta análise numérica mostram que uma abordagem adequada do DFC pode ser utilizada para melhorar o desempenho da bomba ejetora.

Na presente investigação, o movimento transiente de um vaso cilíndrico, interagindo com um fluido viscoso e um núcleo cilíndrico sólido central, é simulado usando-se um elemento de casca axisimétrico com dois pontos nodais para a parede do cilindro, e elemento isoparamétrico de oito pontos nodais para o líquido. Aplica-se o método dos pesos residuais de Galerkin para as equações de Navier-Stokes expressas em coordenadas polares. Para a casca cilíndrica, emprega-se a teoria clássica linear de Novozhilov. O problema do acoplamento dinâmico transiente é resolvido com o auxílio do procedimento de Newmark da aceleração média para a casca, e no caso do elemento de fluido, emprega-se o esquema "Diferença Retrógrada" de Euler. Para a investigação teórica, são realizados cálculos para determinar o deslocamento, velocidade e aceleração de pontos da parede externa do cilindro, que é engastada na base, aberta e livre na extremidade, sujeita a uma força inicial, que varia linearmente a partir da base. É observada uma resposta com decaimento monótono harmônico, típico de um sistema amortecido, e dessa forma, estima-se as freqüências naturais de acoplamento. Conclui-se que a acurácia dos resultados são obtidas se o refinamento é confinado à região exatamente abaixo da superfície livre. Um estudo adicional está presente, considerando o "Fluxo Circular em Regime Permanente" de um fluido preenchendo outro vaso cilíndrico rígido com um núcleo central. Nesse estudo as hipóteses adotadas na primeira investigação são mantidas, mas nesse caso é suposto que o núcleo interno impõe um movimento rotativo ao fluido. O comportamento das pressões e velocidades são descritos para várias velocidades angulares do núcleo em duas diferentes malhas de elementos finitos.

Uma análise matemática da interação entre um fluido incompressível e viscoso e o elemento estrutural presente em sua interface é procedida. Para isso, o trabalho fundamentou-se nas equações clássicas de teoria de cascas para a análise estrutural e nas equações de continuidade e de Navier-Stokes para o estudo do fluido. As variáveis no domínio decorrem tanto quando se considera o escoamento em regime permanente como a resposta transiente. No escoamento em regime os termos inerciais estão associados aos coeficientes convectivos, enquanto que na resposta transiente os termos convectivos não aparecem. Uma breve abordagem do problema em regime é descrita de modo a justificar a dificultade quando esta situação é encontrada. Na condição transiente um acoplamento das respectivas componentes do sistema casca-fluido é realizado. O elemento estrutural considerado é a casca axisimétrica de curvatura meridional variável enquanto que o domínio fluido é discretizado através de elementos quadrilaterais de oito pontos nodais. O problema é descrito matematicamente por um sistema de equações diferenciais lineares, onde, por conveniência, está representado através de coordenadas cilíndricas. Um procedimento numérico é utilizado, onde o método de elementos finitos é aplicado para se encontrar a solução. Uma boa concordância é encontrada em termos quantitativos para o escoamento em regime permanente sujeito a um baixo gradiente de pressão. Poucas iterações conduzem aos resultados analíticos. Com relação à resposta transiente, em termos de freqüência, pode também ser verificado boa concordância quando se procede o refinamento da malha de elementos finitos. É considerada a influência da viscocidade no comportamento dinâmico da casca. Em todas variáveis pertencentes ao regime transiente é observada a presença do amortecimento ao longo do tempo.

Um medoto de diferencas finitas para simulacao deescoamentos independente do regime de velocidade considerada ehapresentado. O algoritmo proposto utiliza um esquema implicito defatoracao aproximada em forma delta tipicamente desenvolvido paraescoamentos compressiveis. A habilidade de tratar escoamentos emqualquer numero de Mach e obtida atraves de uma escolha apropriadadas variaveis dependentes, chamadas variaveis primitivas, a saber,pressao, componentes cartesianas de velocidade e temperatura. As equacoes governantes sao escritas em forma de lei de conservacaopara um sistema de coordenadas curvilineas generalizadas esimultaneamente resolvidas por um processo implicito de marcha notempo com linearizacao local por via de matrizes jacobianas. Operadores de diferencas centradas sao empregados para as derivadasespaciais e estabilidade numerica e assegurada pela escolha de ummodelo de dissipacao artificial. O esquema numerico resultante ehsimilar ao algoritmo tradicional de Beam & Warming, com excessao deuma nova matriz jacobiana relacionando variaveis conservadas eprimitivas. O metodo e aplicado a uma grande variedade de testes,incluindo escoamentos internos e externos, e velocidades que variamdo regime incompressivel ao supersonico. Todos os casos analisadosmostram concordancia muito boa com os dados disponiveis daliteratura. A eficiencia computacional do algoritmo proposto etambem avaliada e revela-se tao boa quanto a do de Beam & Warming. Apesar de terem sido considerados apenas escoamentos nao viscosos,ou seja, escoamentos governados pelas equacoes de Euler, nao haqualquer dificuldade conceitual na extensao do presente algoritmopara as equacoes de Navier-Stokes.

A transferência de calor utilizando jatos impingentes tem se estabelecido como uma técnica efetiva para aquecimento ou resfriamento de uma superfície e é empregada em muitas aplicações de engenharia. Na indústria aeronáutica é empregado principalmente em sistemas de proteção, evitando a formação de gelo em superfícies expostas. A compreensão deste fenômeno é muito importante para se prever o desempenho do sistema de proteção. Este trabalho apresenta um estudo numérico de um jato impingindo numa superfície em regime permanente. Foram analisadas duas configurações do jato: um circular impingindo numa superfície plana e um jato expandindo através de uma fenda (slot) impingindo numa superfície côncava. As equações de conservação (massa, quantidade de movimento, energia) e turbulência são resolvidas pela Técnica de Volumes Finitos. Uma análise do caso bidimensional, laminar e turbulento, foi realizada. Há uma boa concordância dos resultados obtidos com os da literatura para o caso laminar. Já para o turbulento, necessitou-se de uma técnica de refinamento sucessivo de malha para se obter convergência e assim concordância com os dados da literatura. Um estudo tridimensional, laminar e turbulento, para o caso da placa plana também foi realizado. Adotou-se um procedimento para geração de malha totalmente hexahédrica para os casos tridimensionais. Os resultados obtidos apresentam boa concordância tanto para o caso laminar quanto para o turbulento. Para o caso da superfície côncava os estudos se concentraram no caso bidimensional, laminar e turbulento. Os resultados mostram que devido a curvatura, o número de Nusselt máximo ocorre no ponto de estagnação tanto para o laminar quanto para o turbulento. Verificou-se que o problema do jato impingente é geometricamente simples mas devido a grande curvatura das linhas de correntes e dos fortes gradientes das variáveis de fluxo, o problema se torna num grande desafio para os modelos de turbulência atuais. Mesmo os modelos de turbulência que foram implementados para o caso de jato impingente precisam ser usados com cautela.