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.

This works presents numerical predictions for the flow and temperature fields formed by two coaxial streams confined in a duct having a gradually varying cross section segment. The well-known marching-forward method of Patankar-Spalding was used for sweepingthe computational domain. The standard k-e model was applied for handling turbulence. Computations are first performed for laminar cases in order to assess the reliability of the numerical code. Both isothermal and heated flow, in a constant cross section duct, were considered. Results are then presented for developed turbulent pipe flow showing good agreement with experimental values. Computations for coaxial confined turbulent jets are then presented and compared with available experimental data. Non-isothermal results follow showing interesting dissimilarity between turbulent kinetic energy and Nusselt number.

The main aim is to understand the gas dynamic process, which occurs inside a two-stage turbine and to propose an algorithm for the simulation of a two-stage turbine at design and off-design points. The explanation of each chapter is constructed to make further works easy, which should be the main of any research work that needs to be continued. Moreover, the structure of the chapters is organized in order to expose the theory in a progressive order, avoiding confusing references to preceding topics. According to this method, the content of each chapter is showed in following topics. The Introduction explains the main goals of the Brazilian Space Agency (AEB), focusing the launching vehicle. Chapter II describes a two-stage impulse turbine designed for a liquid rocket engine and some trubopump basic design criteria. A detailed description of the main losses in two-stage axial turbines of liquid propellant rocket engine is depicted in Chapter III. The explanation are always associated with Russian experimental results showed in graphics or equations relating physical variables neasurable in experiments, and dimensionless factors of losses and efficiency. Experimental evidences have been used to produce experimental equations for the determination of the magnitude of the losses and flow angles in blade rows. Several design solution proposed in Chapter IV for improving the efficiency are based on empirical data. Then, some concepts are introduced in designing gas turbines as: power Nu' specific work Lu'load coefficient Lu.ratio and efficiency hu as function of the ratio of the tangential speed u to the adiabatic speed cad for single-stage an two-stage turbines and turbine specific work LT, turbine load coefficient LT.ratio' overall turbine efficiency hT and turbine power NT as function of the ratio of the tangential speed u to the adiabatic speed cad for two-stage turbines. The flow through nozzles and blade rows is examined in detail in Chapter V, where the process of gas expansion in cascades and features of the flow over blade surfaces are described for subsonic, transonic and supersonic conditions. In Chapter V, it is also discussed the structure of the boundary layer in order to reveal the general laws governing the flow around of blades and the methods used for designing blade profiles to lower shock waves. In the first topic of Chapter VI, the initial parameters of turbine calculations are slected and its values are defined according to initial data analysis. In the second topic, the dimensionless criteria for similarity of gas turbine are described. The basic parameters and the dimensionless criteria for similarity of gas turbine allow the definition of the mass flow and energy charactiristics of the turbines. Special attention is devoted to particular characteristics of high and low-pressure turbines. The relations explained in the preceding chapters are used for implementation of a computer program, using software "Mathcad" for simulation of the turbine working at design and off-design conditions. The algoritms for one and two axial impulse turbines are saved in the files I.mcd and II.mcd respectively in the CD-ROM. The appendix only shows the output data for two-stage axial impulse turbines. In Chapter 7, the validation of the methodology and the computer program is done using test data from Russian bipropellant (kerosene and liquid oxygen) engine RD - 0109 for the design point and using general empirical relations, obtained from linear regression of test data of actual one and two-stage axial impulse turbines for off-design points.

O presente trabalho tem dois objetivos. O primeiro é desenvolver um método numérico para a simulação de escoamentos incompressíveis, e o segundo é aplicar o novo esquema à geometria da cavidade rasa. O método numérico desenvolvido é em diferenças finitas. As equações, escritas em forma de lei de conservação, são discretizadas utilizando diferenças centradas em uma malha co-localizada. Termos de dissipação artificial são introduzidos externamente para controlar o problema do desacoplamento par-ímpar e de eventuais instabilidades não lineares. Na aplicação do método à cavidade rasa bidimensional, 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 é identificar prováveis tendências na troca de calor no piso da cavidade e descobrir qual a condição mais favorável no sentido de se ter o menor efeito convectivo do escoamento externo ao longo do piso da cavidade. Os resultados obtidos demonstram que a troca de calor no interior da cavidade é diretamente afetada pela variação dos parâmetros acima considerados e que existe um regime bem definido no qual a troca de calor por convecção forçada é minimizada. Tal regime se estabelece quando ocorre recirculação sobre todo o piso da cavidade. Esta condição pode ser obtida quando uma única região de recirculação se forma no interior da cavidade ou, também, na situação onde têm-se duas bolhas 'encapsuladas'.

Este trabalho propõe o desenvolvimento e estudo de um modelo matemático do sistema de alimentação do banco de teste de motores foguete a propelente líquido para corpo de prova de 15 kN de empuxo do laboratório de propulsão líquida da APE/IAE (Divisão de Propulsão Espacial/Instituto de Aeronáutica e Espaço) com o objetivo de analisar a influência de efeitos fluido dinâmicos nos principais componentes do motor, levantando resultados da dinâmica de pressão e vazão para ensaios a frio e a quente. Para tal utilizou-se a modelagem em grafos de ligação, importante ferramenta para a simulação dinâmica de sistemas independentemente do seu domínio físico. Esta metodologia é fundamentada em elementos resistivos, capacitivos e de inertância, os quais interagem entre si através de ligações de energia, transmitindo potência e energia, sendo adequada para representação de sistemas hidráulicos. Neste trabalho foi utilizado uma amostra de pressão de ensaio a frio para o corpo de prova do propulsor de 15 kN de empuxo; esta será útil para aproximação da dinâmica de pressão do modelo com a real. Com essa aproximação foi possível simular ensaios a quente para dois pares-bipropelentes, Etanol/LOX e Querosene/LOX. O desenvolvimento deste trabalho e seus resultados contribuem para o melhor conhecimento da dinâmica fluida do sistema de alimentação da bancada de teste, o que torna a ferramenta apta para análise da influência dos componentes hidráulicos e alternativa para aplicação de futuros aperfeiçoamentos do sistema real.

Bombas de alta velocidade, i.e., bombas com velocidade angular entre 300 a 6000 rad/s, são largamente empregadas na aviação, mísseis, área naval, indústria petroquímica, plantas de geração de energia como turbinas a gás ou vapor. Por causa das altas velocidades angulares, elas podem ser aplicadas sem a aplicação de redutores em motores foguetes a propelente líquido. O uso de bombas de elevada rotação se deve ao fato das necessidades de obter elevadas pressões (50 a 60 MPa) com pequenas dimensões, reduzida massa e poucos estágios. O objetivo deste trabalho é apresentar uma metodologia de dimensionamento dos principais parâmetros e características de uma bomba utilizada em Motores Foguetes a Propelente Líquido. A partir das dimensões e características energéticas da bomba de oxigênio líquido, é determinado os parâmetros de desempenho, através da dinâmica dos fluidos computacional. A partir dos resultados numéricos gerados pela simulação são propostos melhorias no projeto a fim de melhorar e otimizar a eficiência da bomba através da análise do comportamento do escoamento no interior da mesma. Para geração da malha do volume foi utilizado o programa Gambit. Para simular o escoamento na bomba centrífuga, as equações de Navier-Stokes com condições de contorno adequadas são resolvidas para a geometria tridimensional através do programa comercial de volumes finitos Fluent, assumindo o escoamento permanente, incompressível e sem a presença de gases dissolvidos no líquido. O modelo ?-? foi utilizado para simular a turbulência no escoamento. Baseado nos resultados obtidos como campo de velocidade, distribuição de pressão para diferentes fluxos mássicos foram analisados.

Difusores de ar é um tema de particular interesse na indústria dos sistemas de ar condicionado e climatização. O difusor swirl (ou tangencial) é um tipo de difusor já utilizado em alguns ambientes climatizados. O seu comportamento é mais conhecido em aplicações no campo da combustão, mas, em aplicações de sistemas de ar condicionado o Número de Reynolds é bem menor, não há a combustão e as condições de contorno são diferentes. Além disso, têm-se poucos estudos voltados para estes difusores num domínio 3D. Com esta motivação, o presente trabalho apresenta as simulações de um difusor tangencial em vazões típicas de aplicações de conforto térmico personalizado, utilizando a modelo de turbulência Reynolds Stress em um domínio 3D. Algumas simulações em um domínio 2D são realizadas a fim de se obter algumas características essenciais do escoamento, como abertura e comprimento do jato. Porém, comparados ao experimento, os resultados 2D precisam ser melhorados. Esquemas de discretização de maior ordem são utilizados para se avaliar o desempenho. Nas simulações no domínio 3D, verifica-se que um resultado melhor é alcançado quando se refina a malha na região central do jato, logo abaixo do difusor. Dois métodos de especificação da condição de contorno de entrada são estudados: o primeiro consiste em utilizar os dados experimentais obtidos na saída do difusor para simular o escoamento sem a geometria do difusor e o segundo simula o difusor completo, aplicando a magnitude da velocidade perpendicularmente à superfície de entrada com base na vazão calculada pelos dados do experimento do PIV (Particle Image Velocimetry). Os resultados numéricos são comparados com o experimento. Verifica-se que o método de simulação sem o difusor apresenta resultados mais precisos com relação ao experimento e apresenta maiores vantagens na simulação numérica. / Air diffusers are a topic of particular interest in the industry of acclimatization and air conditioning systems. The swirl (or tangential) diffuser is a type of device already used in some air conditioned environments. Their behavior is best known in combustion applications, but in air conditioning systems applications, the Reynolds number is much lower, there is no combustion and the boundary conditions are different. In addition, there have been few studies on these diffusers in 3D domain. With this motivation, this work presents simulations of a tangential flow diffuser for applications of thermal comfort. The numerical study uses the Reynolds stress turbulence model in a 3D domain. Some simulations in a 2D domain are performed in order to obtain some essential features of the flow, as the width and length of the jet. However, compared to the experiment, the 2D results need to be improved. Higher order discretization schemes are used to evaluate performance. In 3D domain simulations, it is verified that a better result is achieved when the mesh is refined in the jets central region, just below the diffuser. Two methods of the inlet boundary conditions are studied: the first consists of using the experimental data obtained at the exit of the diffuser to simulate the flow without the geometry of the diffuser and the second method simulates the diffuser completely, applying the velocity magnitude perpendicular to the inlet surface based on the calculated flow rate with experimental data of PIV (Particle Image Velocimetry). The numerical results are compared with experiment. It is noted that a simulation method without the geometry of the diffuser provides more accurate results with the experiment and has major advantages in the numerical simulation.