Turbinas eólicas usualmente são agrupadas em grandes parques, reduzindo o custo de instalação, transmissão da energia e manutenção periódica. A superposição das esteiras sobre turbinas adjacentes normalmente reduz consideravelmente a capacidade total, objeto de estudo de Micrositing. Porém, por vezes o “efeito Venturi” ocasionado pelas turbinas à montante induz maior velocidade às turbinas adjacentes aumentando o potencial eólico disponível nas linhas consecutivas. De forma inovadora empregar o Design Constructal de Bejan, o modelo do disco atuador genérico e a Dinâmica dos Fluidos Computacional (CFD) para obter a melhor disposição geométrica das turbinas em uma área plana e não rugosa, com foco à maior potência extraída por área de turbinas instaladas. Para tal, modelar e predizer o comportamento da esteira é fundamental, assim como conhecer os modelos de esteira e a aplicabilidade dos métodos empregados. O Design Constructal é a fonte dos parâmetros geométricos base das simulações: o espaçamento entre as turbinas e as razões de diâmetros. Após 64 simulações semi-iterativas e mais de 60 iterativas verifica-se que o maior ganho em potência disponível por área é de 7,37% para a configuração V = 7m/s, S = 3D, d/D = 0.5, L = 3D e 8,48% para a configuração V = 11m/s, S = 3D; d/D = 0.25 & 0.5, L= 0.75D, valor relativo à execução de somente um diâmetro de 100 metros. / Usually wind turbines are grouped in large parks, reducing the cost of installation, energy transmission and periodic maintenance. But the overlapping of the aerodynamical wakes on adjacent turbines reduces the total capacity, Micrositing study. However, the "Venturi effect" caused by the turbines upstream sometimes increases the speed to the adjacent turbines increasing the wind potential available in straight lines. Innovatively employing the Design Constructal Bejan, the model of the actuator disc and Computational Fluid Dynamics (CFD) to search the best geometrical layout of the turbines on a roughless and flat area, focus on higher power extracted by area. To do this, model and predict the wake of behavior is fundamental, as well as know the aerodynamical wakes models and the applicability of the methods employed. The Design Constructal is the source of the simulation’s parameters: spacing between the turbines and the diameter’s ratio. After concluded 64 semi-iterative and iterative simulations, and more than 60 verifies, the best gain in available power per area is 7.37% for the configuration V = 7 m/s; S = 3d; d/D = 0.5; L = 3D. And the gain of 8.48% for the configuration V = 11m/s, s = 3D; d/D = 0.25 & 0.50; L = 0.75D, comparing to the implementation of just 100 meters diameter.

O coletor solar de tubos evacuados une uma alta absortividade de radiação solar a um ótimo grau de isolamento térmico. Estas características, aliadas a um custo relativamente baixo, fazem deste tipo de coletor o mais utilizado no mundo. Por isso, diversos tipos de abordagens estão sendo adotadas para descrever seu comportamento térmico. Nesta linha, este trabalho foi desenvolvido através de abordagem experimental e teórica, sendo a última subdividida em numérica por CFD e analítica. A abordagem experimental contou com a construção de uma bancada para medida de temperatura e radiação, em um coletor de 24 tubos evacuados acoplados em um reservatório de 178 L. A abordagem por CFD utilizou um modelo tridimensional transiente. Com o modelo numérico validado, utilizando resultados experimentais, simularam-se diferentes condições de operação, em termos de ângulo de inclinação, fluxo de calor sobre o coletor, tamanho do reservatório e temperatura de entrada da água. Em relação à abordagem analítica, esta é dividida em: modelo de irradiância, modelo do tubo e modelo do reservatório. O modelo de irradiância determina a distribuição da radiação solar ao longo da circunferência do tubo. Parte desta radiação é absorvida pelo coletor e transferida para água. Esta quantidade é determinada com o modelo do tubo, que é baseado no método de resistências térmicas. O modelo do reservatório descreve o comportamento térmico da água em seu interior, tanto em aquecimento quanto em resfriamento, analisando a interação energética com o coletor e com o meio externo. O desenvolvimento do modelo do tubo passa pela avaliação da vazão mássica entre o tubo e o reservatório, além disso, o coeficiente de transferência de calor por convecção no interior do tubo deve ser determinado. Tais variáveis são determinadas a partir de uma correlação para o número de Reynolds, a qual foi obtida com resultados da avaliação por CFD e é função dos números de Nusselt, Grashof e Prandtl. A correlação proposta apresentou bom ajuste com os resultados numéricos. Com a bancada de ensaio experimental foram feitas medidas de temperatura da água no reservatório ao longo de alguns dias. Para as mesmas condições do experimento, a temperatura média da água no reservatório foi estimada com resultados da integração dos modelos de irradiância, do tubo e do reservatório. A diferença entre os resultados experimental e teórico foi de 4,8% para a energia acumulada. / The evacuated tube solar collector combines high solar radiation absorptivity to a great thermal insulation degree. These characteristics, combined with a relatively low cost, make this type of collector the most used in the world. Therefore, various types of approaches are being adopted to describe its thermal behavior. In this way, this work was developed through experimental and theoretical approaches, the latter being subdivided into numeric, by CFD, and analytical analysis. For the experimental approach a test bench was built. The tests was carried on a solar collector with 24 evacuated tubes coupled to a 178 L tank, measuring temperature and solar radiation. The CFD approach used a transient three-dimensional model. After the numerical model validation using experimental data, simulations was carried over different operating conditions in terms of angle, heat flux on the collector, tank size and water inlet temperature. The analytical approach is divided into: irradiance model, tube model and tank model. The irradiance model determines the irradiance distribution of solar radiation along the circumference of the tube. Part of this radiation is absorbed by the collector and transferred to water, this amount is determined with the tube model, using the thermal resistance method. The tank model describes the thermal behavior of inside water, both in heating and in cooling, analyzing energy interaction with the collector and the external environment. The development of the tube model involves the assessment of the mass flow rate between the tube and the tank, furthermore the convection heat transfer coefficient inside the tube must be determined. These variables are determined from a correlation for the Reynolds number, which was obtained with evaluation results by CFD. Proposed Reynolds number is a function of the Nusselt, Prandtl and Grashof numbers. The correlation presented a good agreement with the numerical results. Using the experimental test bench the water temperature was measured into the tank over a few days. For the same experiment conditions, the average temperature of the water into the tank was estimated by results of integration of irradiance, tube and tank models. The difference between the experimental and theoretical results was 4.8% for the accumulated energy.

A simulação numérica de escoamento tem se desenvolvido em diferentes aspectos para melhorar sua exatidão e se aproximar cada vez mais da representação de fenômenos reais. O desenvolvimento científico sobre escoamentos turbulentos multifásicos dependem da simulação numérica, onde o custo computacional destas ainda é um fator determinante, principalmente quando se visa resultados mais precisos. Uma das grandes dificuldades é que códigos numéricos de alta ordem, utilizados para simular escoamentos bifásicos, produzem erros de dispersão numérica presentes principalmente na interface, já que nesta posição se apresenta uma descontinuidade nas propriedades físicas dos fluidos envolvidos. O presente trabalho tem o objetivo de propor uma formulação de alta precisão para escoamentos bifásicos. Para isso, realiza-se modificações em um código já existente, denominado Incompact3d que possui a possibilidade de paralelização para até 100:000 núcleos computacionais, utilizando uma decomposição 2D. Entre estas modificações destacam-se a implementação de uma nova formulação do método Level Set, além de pequenas modificações nas equações de Navier-Stokes, como a adição do termo de gravidade, a segregação do termo de difusividade e a adição de um termo referente a força de tensão superficial. Ainda propõem-se uma nova consideração do termo de pressão, que é separado em pressão não hidrostática e pressão hidrostática. Diversos testes de verificação e validação foram realizados, apresentando a capacidade de se realizar simulações numéricas de alta ordem para escoamentos bifásicos. Com esta metodologia, foi identificar características de fenômenos turbulentos bifásicos. / Flow numerical simulation has been developed in di erent ways to improve its results accuracy and have a better approximation to real phenomenon. To scientific development turbulent two-phase flow depends of numerical simulation, where computational cost is an important factor, mainly when accurate results are required. High order two-phase flow codes issue is that they produce numerical dispersion errors that are amplified with the presence of interface, since in its position exist physical proprieties discontinuity. This research main objective is to propose a high precision formulation to two-phase flow. For that, we modified Incompact3d code that has a 2D decomposition making parallelized calculation up to 100:000 computational cores. Some modifications are the implementation of a new Level Set Method formulation; some minor Navier-Stokes equation modifications as the addition of gravity term; the segregation of di usive term; and the addition of surface tension forces. We also proposed a new pressure term consideration where we separate it in hydrostatic and non-hydrostatic pressure. Many verification and validation tests where developed presenting the ability to perform high order numerical simulations for two-phase flows. Also it was possible to identify multiphase turbulent structures in the phenomena simulated.

Correntes de densidade são fluxos gravitacionais gerados pela diferença de densidade entre dois fluidos. Correntes de turbidez fazem parte de uma sub-classificação das correntes de densidade, na qual o fluido mais denso tem, na sua composição, partículas em suspensão. Muitos trabalhos numéricos já estudaram a dinâmica das correntes de turbidez, mas, nenhum dos encontrados aplicou mudanças de relevo concomitantes com a simulação, causadas pela sedimentação das próprias partículas da corrente e nem alterou o relevo após a passagem de cada evento em um domínio tridimensional. O presente trabalho pretende analisar a alteração no relevo de fundo causada por uma corrente de turbidez. No código Incompact3d, as equações de Navier-Stokes, Continuidade e Transporte e Difusão são resolvidas em uma malha cartesiana tridimensional. A condição inicial adotada é a de Lock-Exchange. As simulações realizadas utilizaram Simulação Numérica Direta (DNS). O código utiliza um esquema compacto centrado de sexta ordem, em diferenças finitas, para o esquema espacial, e Adams-Bashfort de terceira ordem para o esquema temporal. A validação do código foi realizada comparando-se com trabalhos experimental e numéricos. A análise das diferentes proporções granulométricas mostrou que quanto maior é a quantidade de material grosso na condição inicial, maior será seu depósito para um dado tempo. Em consequência, mais relevante se torna a consideração da alteração do relevo de fundo. Além disso, quanto maior o fator de compactação do sedimento, maior será o erro de não considerar a atualização de fundo. Os resultados também apontaram que os erros médios ao não considerar a atualização do fundo são da ordem de 4% da massa de depósito em 20 tempos adimensionais, para os parâmetros utilizados. Ao se propagar uma corrente de turbidez sobre o depósito de outra, os erros se mostram menores. / Gravity currents are gravitational fluxes triggered by density di erence between two fluids. A sub-classification of those are turbidity currents, in which the denser fluid is composed by the lighter fluid plus suspended particles. Many papers had shown turbidity currents dynamics, although none of the papers found had applied changes in the simulated topography due to deposit during the own simulation, neither they had altered a 3D domain topography after each flux, applying the changes caused by the previous current. The present dissertation aims to analyse the turbidity current dynamics alteration caused by the influence of its own deposit, altering the topography during the very simulation. The analysis is conducted in a polidispersed turbidity current. The Incompact3d code solves Navier-Stokes, continuity and transport-di usion equation, in a tridimensional cartesian mesh. Lock-exchange was chosen to be the initial condition. Direct Numerical Simulations (DNS) are performed. Sixth order compact finite-di erence schemes are used on the spatial domain, while third order Adams-Bashfort is applied for the temporal evaluation. Comparisons with numerical and experimental papers were performed for code verification. Results showed the coarser the particles on the starting lock-exchange, the higher its deposit is, and the more the terrain will be altered. Nevertheless, the bigger the compacting factor, the bigger the error of not considering bathymetry alteration. Results also point that the average errors of not considering the update are in order of 4% on the mass deposit, after 20 dimensionless times, for the used parameters. When a current propagates over the deposit of a previous one, these errors are smaller.

Simulações numéricas do escoamento atmosférico em microescala constituem o foco principal deste estudo. Estas simulações são abordadas tendo em vista aplicações para o setor eólico, em especial para avaliações de produção de energia em parques eólicos. Existem diversas categorias de incertezas associadas às estimativas de produção de energia para um projeto eólico, mas na maioria dos casos, a incerteza associada ao modelo de escoamento é a maior e mais relevante de todas. Dentro do setor eólico, o termo “modelo de escoamento” refere-se à ferramenta numérica utilizada para extrapolar o recurso eólico medido na posição das torres anemométricas (e sensores remotos) até as posições projetadas para os aerogeradores. Diversos autores sugerem através de estudos comparativos que os modelos tipo “CFD RANS k-ε” atualmente representam o “estado da arte” para aplicações em parques eólicos e são os mais utilizados comercialmente no setor. Contudo, o escoamento atmosférico livre é intrinsicamente turbulento, e a dinâmica dos escoamentos turbulentos é um campo científico que ainda não foi totalmente dominado pelo conhecimento humano. O presente estudo demonstra que a maioria dos “modelos de escoamento” atualmente disponíveis possuem pontos fracos, em especial quando aplicados em simulações do escoamento atmosférico livre sobre áreas com topografia e rugosidade complexas Uma das fraquezas presentes na maioria dos modelos de microescala para escoamento atmosférico é a “incapacidade” de simular com precisão o escoamento que ocorre durante períodos de “estabilidade atmosférica”. Diversos locais com elevado potencial eólico apresentam ciclos durante os quais as características do escoamento são afetadas pela ocorrência de estratificação térmica dentro da Camada Limite Atmosférica. Tendo como objetivo principal melhorar as simulações do escoamento nestas condições, propõe-se através deste estudo algumas modificações na modelagem “CFD RANS k-ε” tradicionalmente empregada. Dentre estas, destacam-se a inclusão de um perfil estratificado de temperatura potencial como condição de contorno, a inclusão dos efeitos das forças de empuxo no equacionamento “k-ε” e a solução simultânea das equações para balanço de energia e para o fluxo de temperatura potencial. Este modelo foi chamado de “RANS estável”. Para validação deste modelo foram utilizadas cinco torres anemométricas instaladas em um local com topografia complexa. Estas torres foram montadas e instrumentadas conforme as melhores práticas internacionais Os dados anemométricos registrados por essas torres demonstram a presença de ciclos diários de estabilidade atmosférica. Os erros de previsão cruzada foram calculados comparando-se as previsões de cada modelo com as medições reais registradas na posição das torres. O erro global médio de previsão cruzada entre torres anemométricas obtido com a composição dos modelos RANS “estável + neutro” foi de 3,8% enquanto o erro obtido apenas com o modelo RANS k-ε tradicional foi de 5,2%. Para o modelo linear WAsP, amplamente utilizado no setor eólico, o erro foi de 7,1%. Além dos erros de previsão cruzada entre torres, os perfis verticais de velocidade e os fatores de aceleração direcionais obtidos com a composição dos modelos RANS “estável + neutro” também sugerem que esta é uma alternativa versátil e promissora para capturar os ciclos de estabilidade atmosférica utilizando simulações numéricas em regime permanente. / Microscale numerical simulations of the atmospheric wind flow are the central focus of this study. These simulations are analysed from the wind energy perspective. Special attention is given to the usage and application of these simulations in energy production assessments for proposed wind farms. There are multiple uncertainty categories associated with energy production forecasts for future wind farms. However, in most cases the uncertainty factors related with wind flow modelling are the largest and most relevant of them all. The wording “flow model” refers to the numerical simulations (or “models”) that are used to extrapolate the anemometric data recorded at meteorological masts positions to the proposed wind turbine positions. Several authors have demonstrated through comparative studies that the “CFD RANS k-ε” models currently represent the “state of the art” when it comes to microscale wind flow simulations targeted at wind farms. Nonetheless, the atmospheric wind flow is turbulent by nature, and the dynamics of turbulent flows represent one of the scientific fields that have not yet been fully dominated by the human knowledge. The present study demonstrates that the majority of flow models currently available to mankind still lack in precision, even more so when it comes to modelling free atmospheric wind flow over complex terrain. One of the major weak spots of most microscale wind flow models is their inability to precisely simulate the wind flow that occurs during periods of atmospheric stability Numerous locations with large potential for wind energy production present cyclic periods of thermal stratification inside the atmospheric boundary layer. These cycles alter the dynamics and characteristics of the wind stream. With the purpose of improving wind flow simulations under stable atmospheric conditions, some modifications to the standard “RANS k-ε” model implementation are proposed. The most significant of these modifications are the usage of a potential temperature profile among the boundary conditions, the inclusion of the buoyancy forces in the “k-ε” equations and the simultaneous solution of the equations for energy balance and for potential temperature transport. This “modified” model was named “stable RANS”. It was validated using five well mounted meteorological masts installed in a location with complex topography. The anemometric data measured by these site masts suggest the existence of strong daily cycles of atmospheric stability. Cross prediction errors were calculated by comparing the forecasts (outputs) from each flow model against real wind data measured at each mast position The global average cross prediction error yielded by the RANS “stable +neutral” model was around 3,8%, whereas the error yielded by the traditional “RANS k-ε” implementation was near 5,2%. For the linear model WAsP the error was calculated to be 7,1%. In addition to cross prediction errors, the vertical wind speed profiles and speed-up factors calculated with the RANS “stable +neutral” model composition also suggest that it is a promising and versatile alternative for capturing the effects from atmospheric stability on wind flow using steady state numerical simulations.

Polímeros sintéticos ocupam uma posição de grande importância no estilo de vida moderno, servindo como matérias-primas para a construção de uma variedade de utensílios. Apesar do grande número de operações de processamento e produtos disponíveis, o planejamento de produtos e a otimização dos processos de produção raramente constituem-se de tarefas triviais. Isso deve-se ao fato da maioria das operações aplicadas na indústria de processamento de polímeros envolverem geometrias e padrões de escoamento complexos, além da dificuldade intrínseca relacionada ao comportamento reológico complexo de polímeros fundidos ou soluções poliméricas. Devido a estes fatores, o desenvolvimento de técnicas de dinâmica de fluido computacional (computational fluid dynamics – CFD) para a simulação de escoamentos de fluidos poliméricos e etapas de operações de processamento tem sido assunto de numerosos estudos durante as últimas décadas. Sob esta perspectiva, o solver viscoelasticFluidFoam, merece destaque. Ele é capaz de resolver simulações de escoamentos de fluidos viscoelásticos utilizando diferentes equações constitutivas. Contudo, apesar de resultados existentes na literatura apresentarem um bom potencial de aplicação, uma análise extensiva de seu desempenho numérico ainda não foi realizada. Neste contexto, a proposta do presente trabalho é a análise da influência de parâmetros de malha, numéricos e constitutivos no comportamento do solver. As bases para os testes compreendem uma geometria simples – escoamento laminar entre duas placas paralelas – o modelo constitutivo de Oldroyd-B e respectivas soluções analíticas para os campos de velocidade e tensão. Mesmo os testes demonstrando a inegável versatilidade do solver, eles revelam limitações em lidar com algumas configurações de malha e parâmetros constitutivos, principalmente com relação ao refinamento na direção perpendicular ao escoamento, diminuição do número de Reynolds e aumento do número de Weisenberg. Estas limitações podem ser parcialmente contornadas com escolha adequada de parâmetros de relaxação das variáveis e da razão de aspecto dos volumes de controle. Tais dificuldades não estão presentes em simulações de escoamentos de fluidos newtonianos em condições semelhantes, sugerindo que trabalhos futuros devem focar em implementações mais robustas do solver viscoelasticFluidFoam. / Synthetic polymers hold a position of great importance in modern lifestyle, serving as raw materials for the construction of a wide variety of appliances. Despite the large number of processing operations and products available, product planning and optimization of production processes rarely constitute a trivial task. This is due to the fact of operations applied in polymer processing industry involve complex geometries and flow patterns, plus the intrinsic difficulty related to the molten polymers or polymer solutions complex rheological behavior. Because of these factors, the development of techniques of computational fluid dynamics (CFD) for the simulation of flows of polymeric fluids and stages of processing operations has been the subject of numerous studies during the last decades. From this perspective, the viscoelasticFluidFoam solver deserves mention. The solver is capable of resolving simulations of viscoelastic fluid flows using different constitutive equations. However, despite the existing results in the literature present a great potential for application, an extensive analysis of their numerical performance has not been performed yet. The purpose of this paper is to examine the influence of mesh, numerical and constitutive parameters in the behavior of the solver. Bases for the tests comprise a simple geometry – laminar flow between two parallel plates – the constitutive model of Oldroyd-B and its analytical solutions for the velocity and stress fields. Although the tests show the undeniable versatility of the solver, they also reveal limitations in dealing with some mesh settings and constitutive parameters, particularly with respect to refinement in the direction perpendicular to the flow, decreasing in the Reynolds number and increasing in the Weisenberg number. This limitation can be partially circumvented with proper choice of variables relaxation parameters and aspect ratio of the control volumes. Such difficulties are not present in simulations of Newtonian fluids flows under similar conditions, suggesting that future works should focus on more robust implementations of the viscoelasticFluidFoam solver.

A crescente demanda de energia dos dias de hoje, aliada aos graves problemas ambientais advindos da sua produção tem estimulado a procura de formas alternativas de energia. Dentre as diversas fontes de energia, as renováveis ganham expressão particularmente significativa. Neste contexto, a energia eólica revela-se uma aposta correta, a nível econômico, energético e ambiental. Por outro lado, a microprodução é apontada como uma das soluções para os problemas supracitados. Assim, a combinação de ambas as áreas tem sido alvo de especial interesse. Este trabalho visa caracterizar o vento em ambiente urbano, com particular interesse na possibilidade de seu aproveitamento energético. Um dos principais aspectos para determinar a viabilidade técnica e econômica de uma instalação eólica é a avaliação precisa da distribuição das velocidades de vento na área de aproveitamento. A instalação de turbinas eólicas em áreas urbanas tem determinado a necessidade de aprimorar a metodologia de previsão do campo de velocidades do vento visando à melhor determinação da distribuição dos equipamentos e aproveitamento do potencial existente. Neste contexto, esta dissertação apresenta um estudo sobre a aplicação da Dinâmica dos Fluidos Computacional - CFD para avaliação do potencial eólico e o comportamento do vento sobre um modelo de área urbana. Resultados numéricos com diferentes alternativas de modelagem do problema são comparados com dados de um experimento em túnel de vento e dados obtidos com CFD publicados na literatura, visando determinar uma metodologia adequada para avaliação do problema proposto. As simulações numéricas do escoamento de ar em área urbana são realizadas com o uso do programa ANSYS Fluent 13.0, que utiliza o método de volumes finitos para a solução das equações de Navier-Stokes com médias de Reynolds (RANS) e Simulação de Grandes Escalas (LES). O estudo é dividido em dois casos. No primeiro caso, a simulação da área em questão é realizada com as médias de Reynolds de forma a validar o modelo criado. No segundo caso, a simulação é realizada com o método de grandes escalas, cujos resultados serão explorados e comparados com o método RANS. Os resultados obtidos mostram essencialmente duas realidades, a primeira diz respeito à intensidade do recurso eólico e a segunda a sua particular complexidade. / The growing demand for today’s energy, combined with the serious environmental problems of its production, has stimulated the demand for alternative forms of energy. Among the various forms of energy, those that are based on renewable sources have been gaining a particularly significant space for discussion. In this context, the wind energy has proved to be a reliable source, energetically, economically and environmentally speaking. Furthermore, the micro production has been appointed as one of the solutions for the aforementioned problems. Thus, the combination of both these areas has particular interest. This work aims to characterize the wind in the urban environment, with particular interest for the possibility to carry out its energy. In this context this work presents a study on the application of computational fluid dynamics to evaluate the wind potential and the behavior of the wind on urban area. Numerical results with different alternatives for modeling the problem are compared with data from an experiment in wind tunnel and CFD simulation to determine one appropriate methodology for evaluation of the problem. The numerical simulations of the air flow over urban area are performed using the ANSYS Fluent 13.0 which uses the finite volumes method for solving the Reynolds Averaged Navier-Stokes (RANS) and Large-Eddy Simulation (LES). The study is divided in two cases. In the first one, the simulation of the area is performed with the Reynolds Averaged Navier-Stokes to validate the model created. In the second case, the simulation is performed with the Large-Eddy Simulation, which will be explored and compared with the results with RANS method. The results obtained show essentially two realities, the first is related to the intensity of wind resource and the second due to its high complexity. It was noted that in the urban areas turbulence intensity is significantly higher than in the rural environment. The non-horizontal flow components showed a high dependence on geometry of the building. The brief comparison of results obtained by simulation software shows an acceptable approach to the wind velocity field. In terms of turbulence intensity, it’s over estimated by the software. It is concluded from the results obtained that the energy use in urban environment provides a delicate task, due to the low wind intensity and complexity of the resource.

O presente trabalho apresenta a modelagem do escoamento de ar sobre um terreno complexo empregando, e comparando os resultados, três programas comerciais de previsão de potencial eólico, WaSP, Meteodyn WT e WindSim. Por se tratar de um modelo simplificado, o programa WaSP nem sempre é adequado para previsão de energia em terrenos muito complexos pois não é capaz de prever turbulência, separação de escoamento e efeitos de estabilidade existentes neste tipo de terreno. Para que se consiga modelar a turbulência existente de maneira mais coerente, utilizam-se programas de mecânica dos fluidos computacional como o Meteodyn WT e o WindSim que modelam o escoamento através das equações médias de Reynolds. A finalidade do estudo é identificar a capacidade de previsão das velocidades de vento por meio de cada um dos programas em um parque eólico localizado na Chapada da Diamantina, tendo em vista que a correta previsão da produção de energia depende de um entendimento correto e detalhado dos recursos disponíveis. Avaliando apenas condição atmosférica neutra e comparando os resultados obtidos, conclui-se que o programa WindSim obteve o menor erro na estimativa da velocidade do vento, -3,7% quando comparado aos outros programas WaSP e Meteodyn, -5,1% e -6,5% respectivamente. O maior erro de aproximação foi obtido pelo programa WaSP, +11% quando comparado aos outros programas Meteodyn e WindSim, +10,6% e +7,9% respectivamente. Para o programa Meteodyn foram avaliados casos com diferentes condições de estabilidade, o erro máximo e mínimo foi de +10,6% e -6,5%, respectivamente em atmosfera estável, mantendo-se elevado. O programa WindSim permite a alteração de variáveis que são fixas no Meteodyn e também a simulação dos efeitos de estabilidade através da inclusão da equação de energia inicializada através do comprimento de Monin-Obukhov o que, implicou em uma melhora significativa nos resultados obtidos e em um maior tempo computacional. O erro, menor e maior, apresentado pelo programa WindSim foi reduzido a, -1,8% e +4,4% respectivamente. Destaca-se que as conclusões obtidas se aplicam ao caso estudado. / This work studies the air flow modeling over a complex terrain using and comparing the results of three commercial programs that estimates the wind field in a site: WaSP, Meteodyn WT and WindSim. Being WAsP a simplified model, it is not always suitable to estimate the flow in complex terrain because it is not able to predict the turbulence, flow separations and stability effects present in these kind of terrains. To model the existing turbulence two computational fluid mechanics programs i.e. Meteodyn WT and WindSim, which models the flow through the Reynolds Averaged Navier-Stokes equations have been used. This study aims to identify the wind speed predictive capacity of each program in a windfarm located on Chapada da Diamantina, as a correct production prediction depends on a correct estimation of the wind field. Evaluating only neutral atmospheric setup and comparing it with the measures, WindSim had the lowest wind speed estimation error, -3.7% compared with WaSP and Meteodyn programs, that provided -5.1% and -6.5% respectively. The greatest error was obtained by WaSP program, +11%, followed by Meteodyn, +10.6% and WindSim, +7.9%. WindSim enables changes in some variables that are fixed in Meteodyn and enables also the simulation of stability effects through the introduction of the energy equation with a Monin-Obukhov length initialization. For Meteodyn program different stabilities were evaluated, but the results obtained were not satisfactory. The maximum and minimum error were reduced to +10.6% and -6.5%, respectively, in stable atmosphere. Due to the specific site studied conditions the simulation of stability effects resulted in a significant improvement in the results even though a larger computational time. Enabling the energy equations setting up stable atmospheric stability, both WindSim errors are reduced to -1.8% and +4.4% respectively. It is noteworthy that the conclusions apply to the case studied.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Turbomáquinas são máquinas operacionais que transferem energia mecânica entre um rotor e um fluido. Estas máquinas têm muitas aplicações industriais. Um dos componentes de uma turbomáquina responsável pela transferência da energia, ou receber a rotação do eixo e transformar em energia de fluido em caso de bomba ou transferir a energia do fluido para o eixo em caso de uma turbina, é o impelidor ou rotor. O fenómeno da cavitação envolve escoamento bifásico: o líquido a ser bombeado e as bolhas de vapor que são formadas durante o processo de bombeamento. O processo de formação dessas bolhas é complexo, mas ocorre principalmente devido a presença de regiões de pressões muito baixas. O colapso dessas bolhas pode muitas vezes levar a deterioração do material, dependendo da intensidade ou da velocidade de colapso das bolhas. O principal objetivo deste trabalho foi estudar o comportamento hidrodinâmico do escoamento nos canais do impelidor de uma turbomáquina do tipo radial usando recursos de fluidodinâmica computacional (CFD). Uma abordagem Euler-Lagrange acoplada com o modelo da equação de Langevin foi empregada para estimar a trajetória das bolhas. Resultados das simulações mostram as particularidades de um escoamento líquido-bolha de vapor passando em um canal de geometria curva, fornecendo assim informações que podem nos ajudar na prevenção da cavitação nessas máquinas. / Turbomachines are operational machines that transfer mechanical energy between a rotor and a fluid. This type of machinery is employed in many industries. One of the main components of a turbomachine responsible for the energy transference, either receiving the rotation of the shaft and transforming it into fluid energy in the case of a pump or transferring energy from the fluid to the shaft in the case of a turbine, is the impeller or rotor. The cavitation phenomenon involves two-phase flow: the liquid to be pumped and the vapor bubbles which are formed during pumping. The formation process of these bubbles is complex, but occurs mainly due to the presence of regions of very low pressure. The collapse of the bubbles can often lead to a deterioration of the material, depending on the intensity or speed of bubbles collapse. The main objectives of this work was to study the hydrodynamic behavior of the flow in the impeller channels of a turbomachine (radial flow turbopump) using computational fluid dynamics (CFD resources). An Euler-Lagrange approach coupled with the Langevin equation model, was employed to estimate the bubbles tracking trajectory. Results of the simulations show the details of liquid-vapor bubble flow in a curved channel, providing insights that help us in the cavitation prevent of this machines.

Neste trabalho, são apresentados a discussão de conceitos fundamentais, a metodologia e os resultados de simulações numéricas baseadas no Método de Volumes Finitos do escoamento de ar sobre algumas opções de configurações de turbinas eólicas do tipo Savonius, com e sem estatores, em operação e, também, em condições estáticas, como as encontradas nas partidas das mesmas. Comparam-se os resultados para diferentes domínios computacionais, bem como alternativas de discretização espacial e temporal, visando apresentar a influência desses sobre os valores obtidos e estabelecer os parâmetros computacionais adequados para a análise das turbinas em estudo. Nas simulações numéricas, desenvolvidas empregando o programa comercial Star-CCM+, a equação da continuidade e as equações de Navier-Stokes com médias de Reynolds são resolvidas, juntamente com as equações de um modelo de turbulência adequado, que é escolhido, para a obtenção dos campos de pressão e de velocidade do escoamento. Emprega-se um domínio contendo uma região com malha móvel, na qual o rotor é inserido. A cada simulação, a velocidade angular da região de malha móvel é especificada de maneira a variar a razão de velocidade de ponta do rotor. Através da integração das forças ocasionadas devido aos gradientes de pressão e das forças originadas pelo atrito viscoso sobre as pás do rotor eólico, obtém-se o coeficiente de torque em cada simulação. O torque e as forças atuantes no rotor também são obtidos de forma semelhante. Com esses dados, outros parâmetros como a potência e o coeficiente de potência são obtidos. Análises dos principais parâmetros de desempenho do rotor Savonius são realizadas e indicam uma boa concordância com resultados experimentais e de simulações numéricas realizadas por outros autores. Os resultados obtidos nas simulações apresentaram-se bastante representativos do fenômeno analisado. / This research work presents a discussion of basic concepts, the methodology and the results of numerical simulations based on Finite Volume Method for the air flow through some configuration options of the Savonius wind turbines, with and without stators, in operation, and also under static conditions, such as those found in the self starting. Results for different computational domains, as well as alternative spatial and temporal discretization are compared, in order to present the influence of these on the obtained values from the computational analysis of the turbines in study. In the numerical simulations, performed using the commercial software Star-CCM+, the equation of continuity and the Reynolds Averaged Navier-Stokes Equations were solved, together with the equations of a turbulence model appropriate, which is chosen, so that the fields of pressure and velocity could be found. It was used, in the calculations, a domain containing a region with a moving mesh, in which the rotor was inserted. In each simulation, the rotational rate of the moving mesh region was specified so as to vary the tip speed ratio of rotor. Through the integration of the forces arising due to the pressure gradients and the forces originated from the viscous friction on the wind rotor blades, the moment coefficient could be obtained in each simulation. The moment and forces acting on the rotor were also obtained similarly. With these data, other parameters such as the power and the power coefficient of the wind rotor could be obtained. Analysis of the principals performance parameters of the Savonius wind rotor were performed and indicated a good agreement with experimental results and numerical simulations performed by other authors. The simulations results are quite representative of the phenomenon analyzed.