O sangue é um fluido complexo cujas propriedades são significativamente afetadas pelo arranjo, orientação e deformabilidade dos glóbulos vermelhos. O estudo do escoamento sangüíneo é extremamente importante, visto que o desenvolvimento de muitas patologias cardiovasculares estão relacionadas com a variação de suas características frente à circulação sanguínea. Nesse sentido a hemodinâmica têm sido freqüentemente utilizada para obter melhor compreensão dos aspectos funcionais, diagnóstico e terapêuticos. Devido à complexidade da reologia do escoamento sanguíneo, não existe um modelo universal para representar a propriedade viscosa do sangue, entretanto encontramos diferentes classificações que descrevem o comportamento não-Newtoniano do sangue com base em diversas variáveis e graus de precisão. No presente estudo, o método de volumes finitos é utilizado para analisar o escoamento sanguíneo em uma bifurcação da carótida humana. Foi adotado um modelo tridimensional de fluido incompressível em regime permanente e transiente. A equação da quantidade de movimento foi resolvida empregando três funções de viscosidade não-newtonianas: Power-Law, Carreau e Casson, bem como o modelo newtoniano. Foram obtidos resultados para a distribuição da tensão de cisalhamento na parede e perfis de velocidade na direção axial e radial, utilizando o modelo newtoniano e os modelos de viscosidade não-Newtoniana. Posteriormente, o estudo analisou o escoamento sanguíneo sob regime pulsátil utilizando o modelo de Casson. A solução numérica foi obtida utilizando uma condição de velocidade dependente do tempo na entrada da artéria carótida com objetivo de incluir o efeito pulsátil ao escoamento. Os resultados demonstraram que é fundamental a utilização do modelo não-Newtoniano para analisar o escoamento sangüíneo em artérias com variações geométricas por apresentar regiões com baixos valores de velocidade e tensão de cisalhamento.

O trabalho apresenta um estudo do princípio da superposição através da aplicação da Integral de Duhamel para a obtenção da resposta a uma excitação harmônica do ângulo de ataque. Como entrada para a Integral de Duhamel foram utilizadas respostas indiciais obtidas com a utilização de uma ferramenta de CFD desenvolvida pelo IAE/CTA. As respostas indiciais foram obtidas a partir de excitações em degrau de baixa amplitude no ângulo de ataque de um perfil NACA0012 em escoamento com número de Mach igual a 0.3. O estudo mostrou que mesmo que a ferramenta de CFD seja capaz de gerar uma resposta indicial é importante a validação dessa resposta indicial como uma resposta de boa qualidade para a utilização de sua derivada na Integral de Duhamel. Os resultados obtidos com a Integral de Duhamel utilizando uma boa resposta indicial como entrada foram comparados com uma simulação direta de um movimento harmônico realizada com a mesma ferramenta de CFD. A comparação demonstrou que os resultados foram bastante satisfatórios.

Os ejetores possuem diferentes aplicações como na propulsão, refrigeração e aeroespacial. Este dispositivo não possui partes móveis e por isso é um equipamento simples, confiável, com potencial de vida longa e pode ser fabricado a um custo relativamente baixo, além disso, podem ser utilizados materiais baratos para sua construção. Por estas razões, pesquisas vêm sendo realizadas com a finalidade de melhorar o desempenho deste dispositivo. O objetivo do presente trabalho é analisar ejetores supersônicos a gás utilizando a ferramenta de CFD (Computational Fluid Dynamic - Dinâmica dos Fluidos Computacional) para estudar os fenômenos físicos, os efeitos da geometria e das condições de operação no escoamento variando os seguintes parâmetros: vazão na entrada secundária, o comprimento e o tipo da câmara de mistura. As equações do modelo matemático (continuidade, quantidade de movimento, energia e modelos de turbulência k-?-realizável e k-?-SST) são resolvidas numericamente empregando o método de volumes finitos. Um estudo com várias malhas foi realizado para garantir independência dos resultados em relação a estas e para validação do procedimento numérico foram utilizados dados experimentais disponíveis na literatura. Os resultados da validação mostraram que o modelo de turbulência k- ?- SST e o esquema de discretização de 1 ordem apresentaram melhor desempenho para reproduzir os dados experimentais disponíveis para este ejetor. Os resultados para o comprimento da câmara de mistura de pressão constante indicam que para este modelo do ejetor, o comprimento ideal encontrado de acordo com uma equação, pois este é influenciado pela vazão na entrada secundária. Outros resultados também mostraram que com o aumento da vazão na entrada secundária, a diferença de pressão entre a saída e a entrada secundária diminui.

Esta dissertação tem por objetivo simular o escoamento pulsátil de sangue em artéria humana variando-se o grau de estenose (obstrução) e a freqüência do pulso de entrada. O escoamento será modelado resolvendo-se o sistema de equações de continuidade e Navier-Stokes, em regime laminar, com um sistema de coordenadas axisimétrico. O método de volumes finitos é utilizado para a discretização do domínio, com uma formulação acoplada para o cálculo dos campos de pressão e velocidade. São desenvolvidas subrotinas em linguagem C para incorporar a condição de contorno de velocidade oscilatória na entrada da artéria, e compiladas no programa Fluent como UDFs (User Defined Function). Para a simulação com o modelo de viscosidade de Casson, será implementada uma UDF, enquanto o modelo Newtoniano utiliza um valor de viscosidade constante. Serão testados três números de Womersley: baixo (Wo = 4), moderado (Wo = 8) e alto (Wo = 12), o que implica em redução dos períodos, com rápidas acelerações e desacelerações à medida que o número de Womersley cresce. Além disso, varia-se o grau de obstrução da artéria (25%, 50% e 75%), com o objetivo de verificar em que condições os modelos de viscosidade Newtoniano e Casson produzem respostas similares. São obtidos resultados para (i) as linhas de corrente no interior do domínio, (ii) contornos e distribuição de velocidade axial bem como (iii) perfis da tensão de cisalhamento ao longo da parede da artéria. Os resultados permitem concluir que a aproximação entre os modelos de viscosidade Newtoniano e Casson ocorre para alto valor de número de Womersley (Wo = 12) e obstrução de 75%.

O objetivo deste trabalho é o projeto preliminar de uma aeronave comercial focado no pré-dimensionamento e definições aerodinâmicas iniciais. A metodologia utilizada compreende a abordagem clássica combinada com modernas análises de dinâmica dos fluidos computacional para determinar dimensões, peso, tração e perfis aerodinâmicos da asa.

The Computational Fluid Dynamics is extensively employed in aerospace industry, which has become feasible due to the evolution of the computers (increased capacity) and to the relatively low computational cost. A great challenge is the inclusion of viscosity effects in the flow simulations, also using a mesh fine enough to catch all the velocity fluctuations due to turbulence. The demand for computational power makes it unfeasible for aerospace applications. Due to the coming of the turbulence models and Reynolds-averaged Navier-Stokes equations, the use of CFD has become common in aerospace industry, allowing optimization and reducing cycles of wind tunnel testing. This work presents the results of two-dimensional flow simulations made using three different turbulence models, as well as two distinct airfoils - a conventional and a supercritical. The simulated cases were chosen taking into account the availability of data for validation. A discussion concerning the results is presented, as well as the conclusion.

O conhecimento fisico e a possibilidade de calculo do escoamento transonico turbulento em torno de aerofólios é de grandeimportancia para projetos aeronauticos. O trabalho aqui representadoe uma proposta de ferramenta numerica para analise desse tipo de escoamento. A modelagem matematica do problema fisico corresponde a versao tipo camada fina das equacoes de Navier-Stokes com media de Reynolds. A integracao numerica e implementada por meio do algoritmode Beam e Warming, o qual incorpora um esquema centrado de diferencas finitas, com aplicacao de fatoracao aproximada. A marcha no tempo e efetuada segundo o metodo de Euler implicito, com passo no tempo variavel para acelerar a convergencia ate o estado estacionario. Cuidados especiais sao dedicados a implementacao numerica das condicoes de contorno. A malha de calculo e obtida atraves de geradores algebricos. Os efeitos de turbulencia sao simulados atraves da aplicacao de dois modelos tipo viscosidade de vortice: um modelo algebrico devido a Baldwin-Lomax e outro de transporte devido a Johnson-King. Esses dois modelos sao os mais utilizados atualmente em Dinamica dos Fluidos Computacional. Inicialmente, o codigo desenvolvido e validado, utilizando-se para isso escoamentos tipicos de referencia. Posteriormente sao apresentados resultados para escoamentos sobre perfis, comparando-seos valores obtidos entre si e com resultados experimentais e numericos conhecidos. Os desempenhos dos dois modelos sao analisdos e é estudada a influencia de fatores numericos e fisicos sobre as solucoes obtidas.

O processo de condensação de vapor no interior de turbinas danifica sua superfície interna e as palhetas dos rotores, além de diminuir a eficiência termodinâmica nessa região. Para estudar o processo de condensação de vapor em estágios de turbinas, foi proposto neste trabalho uma simplificação da geometria, estabelecendo o escoamento de vapor úmido em uma grade de bocais baseado no modelo de rotor de uma turbina a vapor real. Um estudo preliminar do bocal supersônico de Laval 2D (convergente-divergente) foi realizado devido à possibilidade de validar o procedimento de uso do código CFD com base em dados experimentais adequados. O objetivo desta tese é avaliar o escoamento e investigar como a variação dos parâmetros térmicos, como temperatura e pressão, alteram o ponto de início da condensação, analisando a região de nucleação e crescimento das gotículas entre as palhetas da grade. O fenômeno da condensação foi modelado com base na teoria da nucleação clássica. Nesse trabalho, o modelo matemático deste escoamento compressível multifásico foi resolvido numericamente usando o método dos volumes finitos com uma abordagem baseada na densidade com os modelos de turbulência realizável e SST. Os resultados das simulações mostraram que, tanto para o bocal de Laval quanto para a grade de bocais, o fenômeno da mudança de fase é responsável por perdas termodinâmicas significativas no escoamento e que a fração de massa líquida é influenciada diretamente pela temperatura de admissão do vapor e pela diferença de pressão estabelecida entre a entrada e a saída do domínio, atingindo valores elevados em regiões próximas às superfícies das palhetas.

Por anos o voo dos insetos intrigou os seres humanos. Em alguns aspectos não encontra par entre os voos de outros animais, especialmente nas questões de desempenho e manobrabilidade. Devido a estas características e às deficiências de configurações convencionais de aeronaves em regime de baixo número de Reynolds, têm sido desenvolvidas configurações de micro veículos aéreos (MAVs) baseadas em insetos. Entretanto, ainda há muitos desafios que envolvem o estudo dos mecanismos relacionados com o voo dos insetos. Um entendimento maior destes mecanismos possibilitaria o desenvolvimento de MAVs mais eficientes e eficazes. Neste contexto, a eficiência do movimento periódico inspirado no voo dos insetos foi investigada utilizando duas metodologias: o modelo quase-estacionário e simulação numérica bidimensional. Verificou-se a influência do desvio em relação ao plano de batimento (stroke deviation) horizontal em uma situação de voo pairado, modelando a asa da mosca Drosophila melanogaster por um perfil elíptico bidimensional executando um movimento formando a "figura-de-oito". Novas equações cinemáticas foram propostas para representar a trajetória da asa de modo a permitir uma variação maior de parâmetros do movimento no eixo vertical, o que permitiu obter, por meio de otimização utilizando o modelo quase-estacionário, os limites teóricos da eficiência aerodinâmica. Com base nestas otimizações e em variação dos parâmetros cinemáticos, resultados obtidos por estudos anteriores foram explicados e verificou-se a possibilidade de redução de gasto energético no voo pairado por meio da implementação do desvio do plano horizontal. Simulações numéricas por dinâmica dos fluidos computacional foram realizadas para investigar esta possibilidade e notou-se que a potência do momento aerodinâmico adiciona uma parcela significativa à potência total, o que não é contabilizado no modelo quase-estacionário devido a simplificações. Posteriormente, concluiu-se que, nos casos estudados, as diferenças na eficiência não foram significativas e que, em razão de efeitos não estacionários, em especial a interação entre asa e esteira, a força vertical média se comporta de maneira contrária à prevista pelo modelo quase-estacionário. Finalmente, o mecanismo de captura de esteira foi analisado e foram identificadas mudanças significativas em seus efeitos, induzidas pelo desvio do plano horizontal, as quais puderam ser explicadas devido às alterações na cinemática do perfil.

Um dos componentes principais de uma turbina a gás é a turbina, pois produzem potência para o acionamento do compressor assim como de outras cargas externas. A turbina trabalha em um ambiente hostil, recebendo gases quentes de alta pressão da câmara de combustão, tornando-a um componente crítico em usinas de geração de energia. A turbina afeta o desempenho e o custo total do motor, sendo necessário realizar no início do programa de desenvolvimento do motor a predição de desempenho para garantir que o motor possa satisfazer todos os requisitos de operação de forma eficiente. Com a dinâmica do fluido computacional é possível obter informações do desempenho da turbina axial através da simulação do escoamento entre pás. Para avaliar o desempenho de um estágio de turbina em toda faixa de operação são utilizados mapas de desempenho que serão construídos através das simulações numéricas. A simulação do escoamento viscoso em regime permanente foi realizada no primeiro estágio de uma turbina axial. Adotou-se abordagens bidimensional e tridimensional para resolver as equações de Navier-Stokes com média de Reynolds que utiliza o modelo de turbulência de Spalart-Allmaras. A simulação em vários pontos de operação possibilitou a construção do mapa de desempenho do primeiro estágio da turbina. O mapa foi comparado com uma metodologia de predição de desempenho que contabiliza as perdas por meio do modelo de perdas de Denton.