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 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.

CPL (Capillary Pumped Loop) é um sistema de transporte de calor na forma de um circuito fechado por onde circula um fluido, bombeado de pressão capilar. O calor é transportado utilizando-se os processos de evaporação e condensação que o fluido sofre ao longo do circuito. Neste trabalho, uma nova configuração de CPL é estudada experimental e teoricamente. Com o objetivo de criar uma interface fixa entre o líquido e o vapor dentro do circuito, o condensador convencional tipo serpentina é substituída por um contendo material poroso no seu interior. A idéia é ter um sistema simples, leve e confiável, sem os problemas operacionais do CPL convencional. Um protótipo é fabricado e testado em diferentes condições de operação, demonstrando a funcionalidade e indicando o desempenho deste tipo de CPL. Um modelo matemático, que simula o comportamento térmico e hidrodinâmico do dispositivo em função das condições de operação, é desenvolvido e validado com dados experimentais. Utilizando o modelo, estudam-se as características opercaionais deste sitema, incluindo levantamento das condições necessárias para operação e inicialização, identificação e caracterização do modo de operação, e análise de sensibilidade em relação aos principais parâmetros de projeto.

Dutos retangulares são muito utilizados em dispositivos trocadores de calor. Este trablaho visa investigar, usando um modelo computacional, os efeitos da geração de entropia provocadas pelas irreversibilidades decorrentes do atrito viscoco do escoamento e das resultantes da transferência de calor em escoamentos laminares com perfis térmicos e hidrodinamicamente desenvolvidos. Esta investigação permite avaliar a geração de entropia, a qual é influenciada pelo campo de temperatura no contorno, pelo tipo de escoamento, pelas propriedades do fluido e por parâmetros geométricos. Primeiro, o modelo calcula o campo de velocidade. Em seguida é resolvida a equação da energia, obtendo-se assim, o campo de temperatura do fluido para uma dada distribuição imposta no contorno, não necessariamente uniforme. Uma vez conhecidos o campo de velocidade e de temperatura é calculado o campo de geração de entropia por unidade de volume. Com uma dupla integração é calculada a geração de entropia global por unidade de comprimento na seção analisada. A consistência do modelo computacional é verificada pelo refino da malha e pela diminuição da tolerância. São realizadas validações dos métodos utilizados para os cálculos dos campos de velocidade e de temperatura do fluido. Várias razões de aspecto (altura por largura do duto) são estudadas para vários tipos de fluidos. A influência devida aos efeitos viscosos e da troca de calor sobre a geração de entropia é estudada para vários tipos de fluidos, variando-se o número de Prandtl.

Esta dissertação apresenta um conjunto de equações que modelam a dinâmica dos seguintes elementos básicos de um motor-foguete a propelente líquido (MFPL): câmara de empuxo (cabeçote de injeção e câmara de combustão); linha de alimentação (duto de seção transversal circular ou geometricamente equivalente); constritor (orifício calibrado ou válvula); regulador (regulador de pressão ou regulador de empuxo) e turbobomba (bomba centrífuga e turbina de impulso com um estágio). A abordagem usada neste trabalho é adequada para estudar variações temporais durante a fase de operação normal do motor, isto é, excluindo-se as fases de partida e de corte do mesmo. Adicionalmente, podem ser investigados problemas relacionados a instabilidades de baixa freqüência no sistema hidropneumático, as quais podem degenerar a eficiência do motor ou causar danos ao mesmo. As equações finais, relativas a cada um dos elementos básicos, são implementadas como diagramas de blocos encapsulados usando-se o software Matlab-Simulink e, com eles, é criada uma biblioteca de blocos para aplicações específicas em MFPL. Tal recurso dá a possiblidade de se fazer simulações computacionais de diferentes configurações de motor através de simples modificações na interação entre os blocos, além de permitir se verificar os efeitos causados pela modificação de quaisquer parâmetros do motor. Para ilustrar algumas possíveis aplicações, são estudadas duas configurações distintas: um MFPL pressurizado por gás e um por bomba. Em ambos os casos, respostas a degrau e respostas em freqüência são determinadas e usadas para analisar os efeitos do tempo de preparação da mistura (necessário para injetar, pulverizar, misturar e evaporar os propelentes) sobre a pressão da câmara de combustão.

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'.

Univariate Statistical Process Control (SPC) has been extensively employed in industry for fault detection and product quality control. However, for complex systems, better results have been obtained by the use of more elaborate Multivariate Statistical Process Control (MSPC) techniques. This work proposes the use of MSPC techniques (Hotelling T2 and Runger U2 statistics) for health monitoring of the Embraer 190 horizontal stabilizer system. For illustration, these techniques are applied in a simulated case study representing a degradation scenario with progressive increase in friction. For comparison, Univariate SPC is also employed. The results indicate that MSPC provides earlier indications of degradation as compared to Univariate SPC. In addition, a sensitivity analysis reveals that the MSPC methods are more robust with respect to changes in the time profile of the degradation, as well as alterations in the influence of temperature on the friction coefficient. In all cases, the best results are obtained by using the Runger U2 statistic.

Este trabalho tem como objetivo a comparação de metodologias normalmente empregadas pelas EMBRAER (Empresa Brasileira de Aeronáutica) na determinação do arrasto de onda de asas operando em regime transônico. Os métodos avaliados neste trabalho são métodos que, devido à sua boa precisão e rápida velocidade de processamentos, são bastante adequados para serem utilizados no desenvolvimento de uma aeronave, principalmente durante seu projeto preliminar. A parcela do arrasto de onda é bastante significativa para aeronaves voando com velocidade acima de Mach 0.6 e sua importância é cada vez maior à medida que sua velocidade aumenta. No entanto, devido ao grande número de configurações possíveis para a aeronave em fase inicial de concepção, torna-se bastante difícil e custoso a determinação do arrasto de onda através de cálculos computacionais complexos ou ensaios em túnel. Técnicas que permitam uma estimativa rápida, porém confiável, da importância deste tipo de arrasto devem então ser empregadas. A metodologia utilizada deve ser de fácil aplicação, de modo a se avaliar, em um curto espaço de tempo, diversas configurações possíveis e, deste modo, auxiliar na escolha da melhor configuração. Tendo isso em mente, neste trabalho foi analisada a metodologia descrita pelo ESDU, cuja formulação remete a resultados já consagrados na indústria aeronáutica e/ou baseados em procedimentos experimentais, e uma metodologia numérica, BLWF, que utiliza-se de malhas computacionais e do cálculo do escoamento potencial. Inicialmente será apresentada uma breve descrição dos tipos de arrasto que se fazem presentes durante o vôo, bem como a importância de cada uma dessas parcelas. Em seguida, com finalidade ilustrativa, serão apresentados os cálculos e a distribuição do arrasto total em uma aeronave teórica voando na faixa transônica. As descrições dos métodos empregados são apresentada nos capítulos 5 e 6, que é seguido por uma confrontação e análise dos resultados obtidos utilizando-se as duas metodologias apresentadas. Os resultados obtidos pelo ESDU, embora extremamente rápidos, apresentam boa precisão para cálculos preliminares. No entanto, para asas que utilizam perfis supercríticos, a metodologia ESDU tende a sub-avaliar a contribuição do arrasto de onda. O BLWF, por exigir um maior número de parâmetros de entrada e tempo de simulação, é recomendado para uma fase já mais avançada do projeto. Deve-se levam em consideração que, apesar de importante, o arrasto da aeronave não é o único fator relevante na escolha de uma determinada configuração. Parâmetros como qualidade de vôo e até mesmo custos de fabricação jamais devem ser negligenciados.

O acúmulo de gelo nas superfícies é uma preocupação, pois pode causar deterioração considerável na performance aerodinâmica e no controle da aeronave, comprometendo a segurança da mesma. A forma e localização do gelo dependem das condições meteorológicas e de fatores específicos do avião. A operação segura de aeronaves em condição de gelo requer o uso de sistemas de proteção contra gelo que são tipicamente instalados nas superfícies aerodinâmicas críticas como asas, empenagem e entrada do motor. Para prevenir a formação do gelo, geralmente, os sistemas de degelo usam meios mecânicos ou térmicos enquanto que sistemas de antigelo utilizam fontes térmicas. O estudo para o conhecimento do fenômeno físico da formação de gelo e possíveis soluções para o mesmo, pode ser feito através de técnicas experimentais e também através de modelos computacionais. Em se tratando deste último, a resolução de problemas de escoamentos pode ser feito através da dinâmica dos fluidos computacional, do inglês CFD Computational Fluid Dynamics. Os benefícios do método estão na possibilidade de se modelar os fenômenos físicos fluidos que não podem ser facilmente simulados ou medidos com um experimento e também a possibiliade de se investigar o problema de forma mais rápida do que com os procedimentos experimentais, reduzindo os custos de projeto. Neste contexto o objetivo do presente trabalho é estudar o fenômeno físico de um sistema de proteção contra gelo com ar quente para assim obter a temperatura na superfície de um perfil com um tubo piccolo via mecânica dos fluidos computacional e estudar a viabilidade de adaptação de malha sem considerar a presença de camadas de prisma na superfície. Os resultados da simulação com a malha refinada apresentam boa concordância com os resultados experimentais, o que demonstra a confiabilidade dos modelos computacionais. O critério de adaptação de malha empregado neste trabalho, e especificamente para as condições aqui simuladas, não se mostrou favorável para a resolução do problema.