Neste trabalho são apresentados os detalhes de projeto e construção de um equipamento de tocha de plasma gerada por uma descarga em arco de baixa potência, 300 - 2000W, para sintetização de filmes de diamantes "CVD", assim como são comentados os resultados da caracterização do arco elétrico, da tocha de plasma e dos filmes produzidos. O sistema foi projetado e construido de acordo com a dinâmica de gases unidimensional, teve como diretriz a necessidade de se utilizar componentes e equipamentos disponíveis no Laboratório Associado de Plasma do INPE. A descarga foi caracterizada eletricamente, ou seja, foram determinadas as curvas características da descarga, tensão em função da corrente, para diversas condições de operação, as quais foram comparadas com as curvas características obtidas por um modelo teórico simples desenvolvido para este sistema apresentando boa concordância. Verificou-se que o sistema de descarga apresenta dois modos de operação, que foram denominados modo luminescente e modo arco. No modo luminescente foram determinados os parâmetros da tocha de plasma: densidade e temperatura eletrônica usando sonda eletrostática dupla e no modo arco foi determinada a freqüência angular de rotação do arco utilizando um anodo dividido em quatro setores independentes. Foi realizado também uma análise da superfície dos eletrodos. Tanto no modo luminescente como no modo arco foram crescidos filmes de diamante que foram caracterizados pelas técnicas de microscopia eletrônica de varredura, dispersão em energia de raio-X, difração de raio-X e espectroscopia Raman. Os filmes produzidos no modo luminescente apresentam uma fase cristalina misturada com uma fase amorfa e pouca uniformidade ao passo que os filmes produzidos no modo "arco" são bem uniformes porém apresentam somente a fase amorfa.

Este trabalho apresenta métodos para simulação de escoamentos a altos números de Mach utilizando diversos esquemas de discretização espacial. Dentre estes esquemas, considerou-se um esquema centrado e os esquemas de separação de vetores de fluxo de va Leer e de Liou. Os esquemas "upwind"foram implementados tanto em suas versões de primeira ordem quanto de segunda ordem de precisão espacial. As versões de segunda ordem de precisão estão baseadas numa abordagem "MUSCL", onde os estados à direita e à esquerda são reconstruídos linearmente e afetados por um limitador. O trabalho também discute uma técnica de geração de malhas não estruturadas bidimensionais, baseada no método de avanço de frente de geração. A malha é gerada através da introdução automática de pontos de domínio computacional. Utilizam-se termos fonte locais, que auxiliam na concentração de pontos em regiões específicas de domínio. Foram implementadas algumas técnicas de melhoramento de malha, tais como troca de diagonais e suavização. A metodologia de simulação de escoamento implementada utiliza ainda um procedimento de refinamento adaptativo, sendo que o processo de refinamento de triângulos usa um sensor baseado em propriedades físicas do escoamento. Os resultados foram obtidos através da simulação do escoamento em uma configuração típica de entrada de ar supersônica/hipersônica. Os diferentes esquemas de discretização espacial utilizados neste trabalho foram capazes de reproduzir adequadamente as características do escoamento.

Neste trabalho são apresentadas soluções numéricas para o escoamento laminar e turbulento em um canal aletado com material poroso, para um escoamento de fluido monofásico em regime permanente, no qual a região porosa é modelada como sendo homogênea e isotrópica. A condição de periodicidade espacial ao longo do canal é empregada. As equações do movimento e continuidade de massa são integradas em um volume elementar representativo acarretando em um único conjunto de equações governantes válido para todo o domínio computacional. Estas equações são discretizadas pelo método de volume de controle e o sistema de equações resultante é resolvido pelo método SIP [Strongly Implicity Procedur

In this study the influence of physical properties on heat transfer between the fluid and solid phases for laminar and turbulent flow in a square cavity filled with porous material, heated by one side and cooled by the opposite one is investigated. In order to simulate flow and heat transfer between the phases models, one and two-energy-equation models are used. The transport equations are discretized using the control volume method and the system of algebraic equations is relaxed via the SIMPLE algorithm. Validations were perfomed first for the case of clean cavity, using Laminar and High Reynolds Turbulent model, and the results are in agreement with the existing literature. In the case of porous cavities, simulation runs were performed using the 1EEM and 2EEM. Overall, this study showed that as porosity increases, the value of the average Nusselt at hot wall, for the same Ram, decreases. It is also shown that the average Nusselt decreases, from a critical Rayleigh , as the ratio between the thermal conductivities of the phases increases. Comparisons were made between both energy models and 2EEM shown smaller results for average Nusselt as ks/kf increases. Also found that the critical Rayleigh depends directly on the ratio ks/kf, being larger as it also grows. Still object of this work, it has been discovered that as the porosity decreases and the ratio between thermal conductivities increases, the turbulence within cavity increases as well. In general, the results and analysis obtained by this work can be used in real engineering situations where the porous cavities may be used.

Sabe-se que um bom projeto de coletor de admissão é um fator importante para um bom desempenho de um motor de combustão interna. Modelos computacionais têm sido usados para simular escoamentos não-permanentes e compressíveis fornecendo assim suporte para melhor projeto do coletor. Neste trabalho, o estudo de escoamento compressível, intermitente e turbulento de fluido newtoniano no interior de cilindros de seção circular tem por objetivo entender os fenômenos que ocorrem dentro de um sistema complexo de dutos que conduzem o ar até os cilindros de um motor. Para tal é necessário primeiramente conhecer quais são estes fenômenos e identificar os parâmetros cruciais que influenciam no processo de solução do modelo matemático em uma geometria simplificada e assim possibilitar a avaliação do sistema em operação em uma geometria complexa. O modelo matemático usado para o escoamento compressível pulsátil dentro do duto circular reto é o URANS (Unsteady Reynolds Averaged Navier-Stokes). Para resolver as equações governantes do modelo matemático (continuidade, quantidade de movimento, energia e modelo de turbulência) é usado o software ANSYS Fluent CFD. Esta ferramenta é baseada na discretização de volumes finitos. Este trabalho é composto por três etapas. Na primeira etapa a abordagem baseada na pressão é usada resolvendo as equações governantes de forma segregada. Após um estudo de malha, os resultados obtidos mostraram que um duto excitado por um único pulso tem uma resposta periódica com frequência constante, e certo grau de amortecimento ao longo do comprimento do duto. Foi mostrado que o comprimento do duto é o parâmetro de maior influencia no comportamento dinâmico do sistema. Na segunda etapa a mesma geometria foi empregada e o mesmo tipo de abordagem foi usada. Foi observada a resposta dinâmica de uma condição de contorno de pressão oscilante na saída do duto. É possível verificar a importância da compressibilidade do ar que atenua as amplitudes das oscilações de pressão e vazão em massa ao longo do duto. E finalmente, a terceira etapa consistiu em uma geometria complexa, formada por um duto principal, um plenum e dois runners. Como condições de contorno na saída dos dois dutos secundários foram impostas perturbações de pressão de mesma amplitude e fases variáveis, e assim verificada a vazão em massa para cada defasagem pré-estabelecida.

A aerodinâmica não-estacionária foi muito estudada, principalmente devido aos problemas em Aeroelasticidade. No início do Século XX muitos cientistas como Theodorsen, Sears, Küssner e Wagner se dedicaram a estudar este fenômeno. Tradicionalmente os métodos desenvolvidos para determinar as cargas aerodinâmicas não-estacionárias para escoamentos subsônicos e supersônicos são baseados em formulações lineares. Porém estudos e ensaios mostram que esta premissa deixa de ser válida principalmente no regime transônico, com o resultado da formação de ondas de choque e/ou separações da camada limite causado pela onda de choque do regime transônico, o método linearizado vai ficando progressivamente pior à medida que se aproxima do alto transônico. A mecânica dos fluidos computacional (CFD-computional fluid dynamics) tornou-se uma área de grande importância da Aerodinâmica atual, sendo possível trabalhar com as equações não lineares sem restrições físicas ou geométricas. Este trabalho apresenta o modelamento e análise da resposta indicial de um perfil aerodinâmico exposto a uma perturbação do tipo degrau em ângulo de ataque utilizando ferramentas de CFD. A análise do comportamento quanto a variações de espessura, numero de Mach e amplitude além de, uma comparação com a teoria da aerodinâmica linearizada. Desse modo é possível obter uma sensibilidade quanto à influência desses parâmetros na resposta transiente dos perfis estudados, além de um estudo quanto às capacidades, benefícios e limitações de se utilizar ferramentas de CFD para modelar e estudar fenômenos transientes.

Neste trabalho apresenta-se um estudo cinético do transporte de pequenas partículas, isto é, em escala micrométrica ou nanométrica, em gases diluídos e no regime molecular livre. Com base na teoria cinética de Boltzmann, foram obtidas as forças de arrasto e de termoforese em dois modelos colisionais distintos, denominados especular e difuso. Mostrou-se como adaptar a teoria cinética de Boltzmann, que é válida para colisões moleculares, para descrever as colisões partícula/moléculas. Em especial, analisa-se o fenômeno de termoforese que atua nas partículas quando suspensas no meio gasoso com uma distribuição não uniforme de temperatura. Neste estudo, considera-se explicitamente a influência do potencial de interação entre as partículas e as moléculas do gás, e mostra-se que ele influencia consideravelmente as forças envolvidas no transporte quando se compara aos resultados obtidos com o potencial da esfera rígida. À partir da condição de equilíbrio das forças obtém-se também a velocidade de transporte, denominada velocidade de termoforese, e encontra-se que ela é muito dependente da energia potencial de interação entre as partículas e as moléculas do gás.

The present work is concerned with studying techniques which would allow the identification of a multiple degree of freedom aeroelastic system from a single computational fluid dynamics (CFD) unsteady simulation. This data is, then, used to generate the root locus for aeroelastic stability analysis of the dynamic system. The system being considered in the present work is a NACA 0012 airfoil-based typical section in the transonic regime. The CFD calculations are based on the Euler equations and the code uses a finite volume formulation for general unstructured grids. A centered spatial discretization with added artificial dissipation is used, and an explicit Runge-Kutta time marching method is employed. Unsteady calculations are performed for several types of excitation on the plunge and pitch degrees of freedom of the dynamic system. These inputs are mostly based on step and orthogonal Walsh functions. System identification techniques are used to allow the splitting of the aerodynamic coeficient time histories into the contributions of each individual mode to the corresponding aerodynamic transfer functions. Such transfer functions are, then, represented by rational polynomials and used in an aeroelastic stability analysis in the frequency domain. The work compares the results provided for each case and attempts to contribute with guidelines for such analyses.

A proposta desta pesquisa é realizar um estudo numérico da utilização do etanol como combustível alternativo em turbinas a gás, seguindo a tendência mundial na busca por combustíveis menos poluentes. Esta pesquisa tem como principal objetivo, identificar o impacto da mudança do combustível na metodologia corrente de projeto de câmaras de combustão de turbinas a gás. Como ferramenta principal, é utilizado pacote computacional de CFD - "Computational Fluid Dynamics" - para simulação dos escoamentos reativos no interior dos combustores analisados. O presente trabalho é parte de um esforço para o desenvolvimento de metodologia para projeto de câmaras de combustão de turbinas a gás voltadas a utilização de etanol como combustível, de forma eficiente. Em grande parte, isto depende do posicionamento otimizado dos furos de admissão de ar para as zonas primária e secundária, no que tange o grau de homogeneização da temperatura dos gases à saída do combustor, ou "temperature traverse quality" (T.T.Q.), ou TQ simplesmente. Para tanto, utiliza-se critérios de projeto unidimensionais para projeto preliminar de combustores ';genéricos';, equivalentes a combustores existentes, tomados como referência. Posteriormente, procede-se à verificação da qualidade do escoamento no interior da câmara de combustão através de simulações computacionais (CFD) para cálculo do escoamento 3D, bifásico (spray), viscoso, compressível, turbulento, radiante e reativo. Os resultados demonstram a viabilidade técnica de utilização do etanol como combustível. As dimensões do combustor a etanol são compatíveis com aquelas obtidas no projeto de combustores a querosene de aviação. É possível a identificação do posicionamento ótimo das fileiras de furos das zonas primária e secundária ao longo do eixo longitudinal do combustor. Os valores de eficiência térmica, TQ e de emissões para alguns poluentes são compatíveis com aqueles encontrados na literatura. Resultados dos cálculos em CFD e parâmetros utilizados no dimensionamento preliminar para mínimo TQ são correlacionados ao final.

Este trabalho apresenta um estudo comparativo entre três ferramentas computacionais diferentes na determinação do coeficiente de pressão em um caso típico de aerofólio supercrítico, em particular o Boeing A4, com sua conseqüente comparação com os dados experimentais obtidos em túnel de vento para este aerofólio. Esta comparação contribui para o processo de validação de uma ferramenta de mecânica dos fluidos computacional (CFD) que vem sendo desenvolvida no presente grupo de trabalho. Os códigos de CFD empregados resolvem uma formulação de Navier-Stokes com média de Reynolds (RANS) ou uma formulação de Euler acoplada a soluções de camada limite. Neste último caso, tem-se a opção de fixar a transição sobre a superfície do aerofólio ou deixá-la livre. Todos os códigos utilizam modelos de turbulência do tipo viscosidade de vórtice para levar em consideração os efeitos de transporte de turbulência. Em particular, o modelo de uma equação de Spalart-Allmaras e o modelo de duas equações k-e RT são empregados nos códigos com formulação RANS. O código de camada limite utiliza o modelo algébrico de Cebeci-Smith. O trabalho apresenta os aspectos da integração no tempo e a integração espacial, bem como os procedimentos de aceleração de convergência usados no código que está sendo desenvolvido, o que possibilita numa análise posterior estabelecer relações de tendências de resultados de acordo com os modelamentos e equações que estão sendo utilizados. Foram utilizados onze casos teste diferentes, todos utilizando a configuração do aerofólio Boeing A4. Os resultados obtidos permitiram analisar a qualidade e a confiabilidade dos resultados dos códigos no que se diz respeito à correta identificação de possíveis ondas de choque, valores de coeficientes de pressão, posicionamento da transição de escoamento laminar para turbulento e tendências que os códigos possam apresentar de oscilações numéricas na solução. Estes resultados e análises realizadas dão subsídios importantes no processo de validação do código em desenvolvimento.