Nas últimas décadas, as melhorias na tecnologia de extit{hardware} principalmente no que diz respeito a unidade de processamento central, memória e capacidade de armazenamento, diversas ferramentas numéricas aplicadas na engenharia foram desenvolvidas com o intuito de se obter respostas rápidas e mais baratas do que ensaios, experimentos e testes laboratoriais. A Dinâmica dos Fluidos Computacional (DFC) faz uma aproximação numérica como um meio de simular escoamentos. Um dos grandes desafios em soluções numéricas na DFC é a garantia da consistência, estabilidade e convergência garantindo com isso, a confiabilidade do código numérico. As turbomáquinas operam em uma ampla faixa de velocidades no campo do escoamento, podendo ter regiões com baixos números de Mach. Nesse ponto, surge um problema, pois códigos da DFC para o cálculo de escoamento compressível geralmente são programados seguindo formulações matemáticas baseadas no vetor de variáveis conservadas, podendo ter problemas numéricos e apresentar soluções inconsistentes devido a pequena variação da densidade do fluido. Chorin através de alguns estudos constatou que os códigos de DFC possuem uma diminuição de sua convergência e estabilidade quando o número de Mach do escoamento torna-se menor do que o limite de compressibilidade. Com o objetivo de solucionar este problema, Chorin criou a técnica de compressibilidade artificial. Choi e Merkle; Turkel e outros pesquisadores modificaram a técnica de compressibilidade artificial e desenvolveram a técnica de pré-condicionamento para diversas aplicações para garantir a convergência, estabilidade e melhorias na soluções numéricas dos códigos DFC para uma ampla faixa de velocidades do escoamento. No presente trabalho, foi realizada a implementação do método de pré-condicionamento de Turkel, Vatsa e Radaspiel em um programa de DFC desenvolvido por Tomita para calcular escoamentos compressíveis, internos e externos. A fim de melhorar a robustez do código numérico o pré-condicionador foi modificado seguindo o proposto no trabalho de Damorfal e Siu. O pré-condicionamento foi implementado primeiramente para as Equações de Euler e depois para escoamentos viscosos e turbulentos. Diversas experiências numéricas foram realizadas e resultados interessantes utilizando o código de DFC compressível e o código de DFC compressível com o pré-condicionador original e modificado foram obtidos. Todos estes casos estão apresentados e comentados no presente trabalho. A implementação do pré-condicionador em um código de DFC desenvolvido para o cálculo de escoamento tipicamente compressível resultou em uma ferramenta numérica com maior potencial para o cálculo do escoamento numa ampla faixa de número de Mach.

Este trabalho tem o objetivo de avaliar o desempenho aerodinâmico de uma turbina eólica utilizando as técnicas de dinâmica dos fluidos computacional. Inicialmente, a pá ótima de um gerador eólico de 250kW foi desenvolvido com base na aplicação da Teoria do Elemento de Pá de Glauert, levando em consideração os efeitos de perdas de ponta de pá e rotação da esteira. As características da pá ótima foram definidas através da velocidade do vento e razão de velocidades na ponta da pá. O código utilizado também considera o efeito do ângulo de passo (pitch) e torção geométrica. A avaliação numérica de dinâmica dos fluidos computacional (CFD) foi desenvolvida através do software comercial ANSYS Fluent v.14.5. Os métodos utilizados para a discretização do domínio computacional e as técnicas de processamento, tais como condições de contorno, modelos de turbulência e critérios de convergência, estão descritos nesta dissertação. A análise do desempenho do rotor para várias condições de ângulo de passo, rotação e velocidade do vento médio serão devidamente apresentadas. Vários resultados são apresentados e analisados, incluindo a potência líquida extraída pelo rotor, eficiência aerodinâmica da pá, desenvolvimento da esteira do rotor, camada limite, gradiente de pressão, entre outros.

The object of this work is to apply CFD simulation in the description of the spray burning. As a case study, a pressure swirl injector, characterized and tested by NIST, has been chosen, which atomize liquid kerosene in an atmosphere of gaseous oxygen. The chamber dimensions allow a complete evaporation, avoiding the impact of drops on the circular wall. Swirl-axisymmetric domain and steady state permit to include combustion, a complex process, without requiring of high computational resources. Continuous phase is treated with an Eulerian reference, while fuel drops are tracked following the Lagrangian formulation. Chemical kinetics is reduced to the concept of mixture fraction. This assumption avoids the solution of too many transport equations for all involved species. In the first simulation, the inlet boundary of the continuous phase is obtained from the numerical solution of a fully developed flow transporting the oxidant gas. Then, four cases are proposed and solved, changing the turbulence intensity and swirl velocity on the inlet boundary, each parameter with two different values. Finally, results for the axial velocity, streamlines, drops trajectories, temperature, distribution and total production of selected species are analyzed and compared with other related studies.

A previsão acurada de fenômenos aerodinâmicos não estacionários tem cada vez mais dependido do emprego de técnicas de Mecânica dos Fluidos Computacional (CFD). A aplicação de tais técnicas em problemas não estacionários depende crucialmente de tecnologias de geração de malha que, usualmente, são processos bastante lentos. Para contornar esta limitação, várias metodologias têm sido desenvolvidas a fim de calcular a dinâmica geométrica do problema físico sem gerar novamente uma malha computacional por completo. Neste sentido, este trabalho apresenta um estudo de técnicas eficientes de geração e movimentação de malhas computacionais para simulação não estacionária em CFD, com aplicações em Aeroelasticidade. Para este propósito, é implementada uma abordagem de movimentação de malhas baseada em Funções de Base Radial (RBF). Nesta metodologia, os deslocamentos dos pontos de malha na superfície da estrutura (móvel) são propagados para os pontos volumétricos (interiores), sem alterar a topologia da malha. Este trabalho conduz uma análise comparativa do emprego de diversos tipos de RBF na movimentação de malhas de níveis de refinamento diferentes, bidimensionais e tridimensionais, no que diz respeito à qualidade da malha deformada resultante. Além disso, aborda estratégias de movimentação e sua influência na qualidade de malha gerada. Como aplicação, é conduzida simulação aerodinâmica computacional com código BREXP3D. Este código utiliza formulação de mecânica dos fluidos não estacionária compressível e não viscosa (equações de Euler), implementada numericamente via método de diferenças finitas, o que exige malhas hexaédricas estruturadas. Este trabalho implementa condições de contorno para simulação de asa tridimensional em malhas monobloco e conduz simulação aerodinâmica estacionária em regime transônico, comparando com resultados experimentais. Problemas de malha e convergência não permitiram executar com sucesso a simulação aerodinâmica computacional não estacionária.

A presente dissertação trata da questão importante sobre dispersão no cálculo do coeficiente de arrasto associado a malhas computacionais em aplicações aeroespaciais. Este tema tem sido amplamente discutido na comunidade de Dinâmica dos Fluidos Computacional (CFD) na última década. Determinar com precisão o valor do coeficiente de arrasto de uma configuração sempre foi um dos maiores problemas em CFD na indústria aeroespacial. Tão importante quanto conhecer a grandeza que se deseja calcular é também conhecer quão preciso somos em determinar essa quantidade. Na indústria aeroespacial, programas computacionais baseados nas equações de Navier-Stokes com média de Reynolds (RANS) e malhas computacionais não estruturadas do tipo hexaédricas e do tipo tetraédricas com camadas de prisma para a camada limite aparecem como uma abordagem eficiente. Neste contexto, o presente trabalho conduz estudos de maneira a identificar e a entender as dispersões no cálculo do coeficiente de arrasto causadas por modificações na topologia da malha computacional, utilizando um único código de CFD. A fim de atingir tal objetivo, a asa ONERA-M6 foi escolhida como caso teste deste estudo por representar uma geometria simples e, assim, tornando mais clara a identificação do problema em questão. Com o objetivo de simplificar ainda mais o problema, as simulações não possuem arrasto induzido, que neste caso significa que o ângulo de ataque é zero para todos os casos. Dois números de Mach foram escolhidos para serem investigados neste estudo, Mach 0,3 e 0,8, por representarem regimes subsônico e transônico, respectivamente. Os estudos consideram um número de Reynolds igual a 3 milhões, que é tipicamente usado na indústria aeroespacial durante os estudos experimentais em túnel de vento. As dispersões no cálculo de arrasto encontradas neste estudo foram discutidas e investigadas com o objetivo de identificar as regiões da asa de maior contribuição da dispersão total e, também, apontar as diferenças na malha superficial e volumétrica que poderiam contribuir para tais dispersões. Baseado nos resultados obtidos neste trabalho, foi possível identificar a precisão no cálculo do coeficiente de arrasto de uma configuração simplificada, devido às modificações na topologia da malha que podem ser associadas a variações geométricas em um estudo de otimização aerodinâmica.

Este trabalho aborda a utilização de uma geometria representativa de um motor mono cilíndrico, com injeção direta de etanol hidratado e ignição assistida por centelha. A partir da utilização de um código comercial de Dinâmica dos Fluidos Computacional, a malha móvel é criada e as equações governantes são resolvidas. O escopo do trabalho inclui a caracterização do escoamento de ar durante a compressão e a formação de mistura estratificada. É analisado o efeito da geometria da câmara na formação da mesma e sua distribuição ao longo da câmara de combustão. É analisada a influência de parâmetros como: pressão do combustível injetado, início da injeção e massa de combustível injetada. O objetivo final é a obtenção da concentração de combustível ao longo do cilindro, durante o desenvolvimento do jato e o movimento do pistão, além da formação de filme nas paredes da câmara. Tais resultados são de extrema importância para a avaliação de estratégias de injeção que tornem possível a partida a frio de veículos em baixas temperaturas sem a aplicação de sistemas auxiliares. Os resultados, como já eram previstos, mostraram que a pressão e o início da injeção de combustível são parâmetros que influenciam diretamente a evaporação de combustível, assim como a quantidade de filme formado, principalmente no pistão.

Submitted by William Justo Figueiro (williamjf) on 2015-07-08T23:40:19Z No. of bitstreams: 1 19b.pdf: 2923099 bytes, checksum: f04244aa991c1a65ebc9273c29c4aa38 (MD5) / Made available in DSpace on 2015-07-08T23:40:19Z (GMT). No. of bitstreams: 1 19b.pdf: 2923099 bytes, checksum: f04244aa991c1a65ebc9273c29c4aa38 (MD5) Previous issue date: 2012-11 / Nenhuma / Nas últimas décadas, a miniaturização de dispositivos de natureza eletroeletrônica em diversas áreas de aplicação, como biomédica, química e de tecnologia de computadores, tem proporcionado alta eficiência de espaço em equipamentos. Paralelamente, essa redução em espaço físico é contrabalançada pelo alto desempenho exigido com relação aos sistemas de refrigeração nestes equipamentos. Por isso, o controle térmico é uma das áreas mais críticas para o desenvolvimento da microeletrônica moderna. Diversos estudos experimentais e numéricos foram realizados por vários pesquisadores, nas últimas décadas, com vistas a investigar a hidrodinâmica e a transferência de calor em microescala. Porém, os resultados obtidos revelam divergências entre si. Em geral, estas podem ser visualizadas pela análise dos números de Poiseuille e de Nusselt, quando estes são comparados aos resultados previstos pela teoria convencional. Normalmente, as divergências relatadas com relação a medições em microescala estão associadas a fatores geométricos do microcanal, como a sua razão de aspecto, o seu diâmetro hidráulico e condições superficiais de rugosidade. Alguns desvios encontrados com relação ao fator de fricção foram atribuídos às variações da seção transversal dos microcanais. O objetivo deste trabalho foi analisar numericamente como as características hidrodinâmicas e de transferência de calor podem ser influenciadas com relação a variações na seção transversal dos microcanais. Os resultados obtidos para o escoamento laminar monofásico de água em microcanais com deformidades na seção transversal foram comparados aos obtidos para microcanais geometricamente perfeitos, por meio dos números de Poiseuille e de Nusselt local. Desvios para os números de Poiseuille e de Nusselt local, por meio de imperfeições e variações da seção transversal dos microcanais, foram verificados. Alguns desvios para o número de Poiseuille se mostraram dependentes do número de Reynolds. No entanto, alguns resultados obtidos para o número de Nusselt local mostraram que o mesmo é mais sensível à forma da seção transversal do microcanal do que o número de Poiseuille. Esses resultados foram determinados através das equações de conservação da massa, Navier-Stokes e energia, por dinâmica dos fluidos computacional (CFD). / In recent years, the reduction of electronic devices in several application fields, such as biomedicine, chemistry and computer technology has been providing high efficiency related to the space in equipment. At the same time, this reduction in physical space is counterweighted by the high performance required at the refrigeration systems in such equipment. Therefore, the thermal control is one of the most critical areas for the development of modern microelectronic devices. A lot of experimental and numerical studies have been done by several researchers, in the last decades, seeking to investigate the hydrodynamic and heat transfer in microscale. However, the results show divergences among them. In general, these related diversions can be viewed through the Poiseuille and Nusselt numbers, when compared to the predicted results through conventional theory. Usually, the reported divergences related to the measurements in microscale are associated to the microchannel geometric factors, such as channel aspect ratio, channel hydraulic diameter and surface roughness. Some deviations related to the friction factor were attributed to cross-section variations of the microchannels. The aim of this work was verify numerically how the hydrodynamic and heat transfer characteristics can be influenced by cross-section variations of the microchannels. The results obtained for the single-phase laminar flow of water in microchannels with deformities at the cross-section were compared to the perfect ones, through Poiseuille and local Nusselt numbers. Deviations at Poiseuille and local Nusselt numbers, through imperfections and variations of the cross-section of microchannels, were verified. Some deviations at Poiseuille number were dependent on Reynolds number. However, some results obtained for the local Nusselt number showed that it is more sensitive to the shape of the cross-section of the microchannel than Poiseuille number. Such results were obtained through mass conservation, Navier-Stokes and energy equations, by computational fluid dynamics (CFD).

No presente trabalho são empregadas ferramentas de fluidodinâmica computacional (CFD) como base à criação de uma rede de reatores equivalente para o estudo de um reator de leito fluidizado borbulhante em escala piloto pertencente à CIENTEC-RS, empregado para a gaseificação de carvão das minas de Candiota-RS. Inicialmente, foram conduzidas simulações fluidodinâmicas em meio reacional gasoso (utilizando metano e ar) para análise do comportamento do escoamento no gaseificador. A partir de medições do equipamento real e dados de projeto do mesmo foi implementada no software CFX 13.0 a geometria do reator e foram realizados testes de convergência de malhas computacionais, de forma a escolher a mais apropriada para a realização das simulações em CFD. A avaliação de parâmetros como perfis de velocidade e concentrações ao longo do reator serviram para a escolha das malhas para a continuidade do estudo. Diferentes variações na geometria foram propostas, de modo a simplificá-la, reduzindo o custo computacional das simulações, porém mantendo a qualidade dos resultados. Foram, ainda, avaliados os padrões de escoamento típicos do reator, operando em fase gasosa e na presença de um leito fluidizado. De posse dos padrões de escoamentos obtidos foi proposta uma rede de reatores equivalentes – ERN (equivalent reactor network) a fim de modelar o reator da CIENTEC, empregando um modelo cinético detalhado para determinação das taxas de reação na fase gasosa. Para fins de avaliação da metodologia, conduziram-se simulações envolvendo a combustão de metano, analogamente ao estudo em CFD, e envolvendo o processo de gaseificação de carvão, com as características do obtido nas minas de Candiota. O objetivo do estudo foi avaliar os pontos cruciais no desenvolvimento de um modelo ERN e prever as características dos produtos obtidos para diferentes condições operacionais. Espera-se que este modelo seja aprimorado e usado para o melhoramento e otimização do processo, com o propósito de obter um gás de potencial energético e com características adequadas à utilização como insumo para processos carboquímicos. / In this dissertation, computational fluid dynamics (CFD) simulations were applied to develop equivalent reactor networks (ERN) for the study of a bubbling fluidized bed reactor from CIENTEC-RS, used for coal gasification. Initially, the computational geometry of the gasifier was created using the software CFX 13.0, based on field measurements and design data. Mesh convergence tests were carried out in order to choose the most appropriate for obtaining accurate solutions, using gas-phase reactions. Different variations in geometry were proposed in order to reduce the computational cost of the simulations, while maintaining the quality of the results. The simulations also evaluated the typical flow patterns of the reactor operating in the gaseous phase and in the presence of a fluidized bed. The flow patterns were used as a basis for proposing equivalent reactor networks for the CIENTEC reactor, using a detailed kinetic model to determine the reaction rates in the gas phase. For methodology evaluation purposes, simulations involving the combustion of methane were conducted, similar to the study in CFD, as well as the coal gasification process. The aim of the study was to evaluate the crucial points in the development of an ERN model and to predict the characteristics of the gaseous product streams obtained for different operating conditions. It is expected that this model can be used for the improvement and optimization of the gasification process, in order to obtain a gas with energetic potential and suitable for using it as a feedstock for carbochemical processes.

Micrositing é o campo do posicionamento estratégico dos aerogeradores na área do parque eólico visando a con guração mais promissora em termos econômicos ou de produção. Dado que em terrenos complexos as características do vento variam na área do parque eólico de forma não-linear, emprega-se a modelagem numérica do vento por CFD para extrapolar os dados medidos para toda a região. O presente trabalho consiste no desenvolvimento de uma metodologia de micrositing em terrenos complexos capaz de auxiliar no projeto do layout e seleção do tipo e altura de eixo do aerogerador que maximiza o fator de capacidade (FC) utilizando softwares livres de código aberto. A metodologia consiste na simulação do vento para várias direções de incidência, assimilação das medições, convers ão de velocidade em densidade de potência, ponderação por frequência de ocorrência de cada direção, sobreposição, seleção das coordenadas dos aerogeradores e cálculo do FC para diversas con gurações de tipo e altura de eixo dos aerogeradores. Veri cação, validação e seleção das constantes do modelo de turbulência é realizada anteriormente às simulações Veri cou-se que o modelo k- produziu um escoamento horizontalmente homog êneo e que o melhor desempenho na validação foi obtido com a escolha de constantes para escoamentos atmosféricos. A metodologia foi demonstrada em uma região de terreno complexo em que o FC do parque eólico proposto apresentou caráter convergente com o re no progressivo da malha, porém oscilatório em termos do número de direções simuladas. Por m, obteve-se FC brutos superiores a 40% para as cinco melhores con gurações e de aproximadamente 52% no melhor caso, indicando bom potencial eólico. A metodologia foi capaz de preencher uma lacuna na literatura cientí ca de micrositing ao possibilitar o planejamento do layout, tipo de aerogerador e altura de eixo, bem como a estimativa da produção e FC brutos da usina em terrenos complexos. Além disso, a estrutura de trabalho com o uso de recursos computacionais livres e de código aberto reforça o caráter de desenvolvimento contínuo, compartilhamento e transparência da metodologia. / Micrositing is the eld concerned with the strategic positioning of wind turbines in the wind farm area aimed at the most promissing con guration economically- or yield-wise. Given the wind characteristics vary non-linearly across the wind farm area in complex terrain, numerical wind modeling with CFD is employed to extrapolate the measured data to the whole site. The present work consists of the development of a micrositing methodology in complex terrain capable of assisting the layout project and selection of wind turbine type and hub height that maximizes the capacity factor (CF) using free and open-source software. The methodology consists of simulating the wind for a number of incoming directions, assimilation of measurements, conversion of wind speed into power density, weighing by frequency of occurrence of each direction, overlapping, selection of wind turbine coordinates and CF calculation for a number of wind turbine types and hub heights. Veri cation, validation and selection of turbulence model constants is performed previous to the simulations It was veri ed that the k- model is able to sustain horizontally-homogeneous ow and that the classic atmospheric ow constants performed best in the validation step. The methodology was demonstrated in a complex terrain region for which the proposed wind farm CF showed converging behavior with progressive mesh re nement, however oscillating with the number of wind directions simulated. Ultimately, CF greater than 40% were obtained with the ve best performing con gurations and approximately 52% in the best case scenario, suggesting good wind potential. The methodology was capable of lling a major gap in the scienti c literature of micrositing for allowing the layout planning, selection of wind turbine type and hub height, as well as gross production estimates and CF for the wind farm in complex terrain. Additionally, the free and open-source-based framework strengthens the continuous development, sharing and transparency of the methodology.

Orientador: Helcio Francisco Villa Nova / Dissertação (mestrado) - Universidade Federal do ABC. Programa de Pós-Graduação em Engenharia Mecânica, 2013