Bombas de alta velocidade, i.e., bombas com velocidade angular entre 300 a 6000 rad/s, são largamente empregadas na aviação, mísseis, área naval, indústria petroquímica, plantas de geração de energia como turbinas a gás ou vapor. Por causa das altas velocidades angulares, elas podem ser aplicadas sem a aplicação de redutores em motores foguetes a propelente líquido. O uso de bombas de elevada rotação se deve ao fato das necessidades de obter elevadas pressões (50 a 60 MPa) com pequenas dimensões, reduzida massa e poucos estágios. O objetivo deste trabalho é apresentar uma metodologia de dimensionamento dos principais parâmetros e características de uma bomba utilizada em Motores Foguetes a Propelente Líquido. A partir das dimensões e características energéticas da bomba de oxigênio líquido, é determinado os parâmetros de desempenho, através da dinâmica dos fluidos computacional. A partir dos resultados numéricos gerados pela simulação são propostos melhorias no projeto a fim de melhorar e otimizar a eficiência da bomba através da análise do comportamento do escoamento no interior da mesma. Para geração da malha do volume foi utilizado o programa Gambit. Para simular o escoamento na bomba centrífuga, as equações de Navier-Stokes com condições de contorno adequadas são resolvidas para a geometria tridimensional através do programa comercial de volumes finitos Fluent, assumindo o escoamento permanente, incompressível e sem a presença de gases dissolvidos no líquido. O modelo ?-? foi utilizado para simular a turbulência no escoamento. Baseado nos resultados obtidos como campo de velocidade, distribuição de pressão para diferentes fluxos mássicos foram analisados.

O presente trabalho tem por objetivo desenvolver uma formulação numérica baseada em Análise Isogeométrica para o estudo de escoamentos incompressíveis isotérmicos de fluidos newtonianos. Com o emprego desta metodologia, os procedimentos de pré-processamento e análise são unificados, melhorando as condições de continuidade das funções de base empregadas tanto na discretização espacial do problema como na aproximação das variáveis do sistema de equações. O sistema de equações fundamentais do escoamento é formado pelas equações de Navier-Stokes e pela equação de conservação de massa, descrita segundo a hipótese de pseudo-compressibilidade, além de uma equação constitutiva para fluidos viscosos de acordo com a hipótese de Stokes. Para problemas com escoamentos turbulentos emprega-se a Simulação de Grandes Escalas - LES (Large Eddy Simulation), na qual o modelo clássico de Smagorinsky é utilizado para a representação das escalas inferiores à resolução da malha. O esquema explícito de dois passos de Taylor-Galerkin é aplicado no contexto da Análise Isogeométrica para a discretização das equações governantes, sendo que a discretização espacial é realizada empregando-se funções NURBS (Non Uniform Rational Basis B-Splines). Essas funções base apresentam vantagens em relação às tradicionais funções utilizadas no MEF (Método dos Elementos Finitos), principalmente no que diz respeito à facilidade de obtenção de continuidade superior a C0 entre os elementos e representação precisa das geometrias. Propõe-se também o desenvolvimento de ferramentas de pré e pós-processamento baseadas na estrutura de dados da Análise Isogeométrica para a geração de malhas e visualização de resultados. Alguns problemas clássicos da Dinâmica dos Fluidos Computacional são analisados para a validação da metodologia apresentada. Os resultados apresentados demonstram boa aproximação da formulação em relação a dados de referência, além de maior versatilidade quanto à discretização espacial dos problemas em comparação com as tradicionais formulações baseadas em elementos finitos. / This work aims to develop a numerical formulation based on Isogeometric Analysis for the study of incompressible flows of Newtonian fluids under isothermal conditions. By using this methodology, pre-processing and analysis procedures are unified, improving the conditions of continuity of the basis functions utilized in the approximations of the equation variables and spatial discretization of the problem. The system of fundamental equations of the fluid flow is constituted by the Navier-Stokes equations and the mass conservation equation, which is described according to the pseudo-compressibility hypothesis. In addition, a constitutive equation for viscous fluids according to Stokes' hypothesis is also provided. Turbulent flows are analyzed using LES (Large Eddy Simulation), where the Smagorinsky’s model is adopted for sub-grid scales. The explicit two-step Taylor-Galerkin method is applied into the context of Isogeometric Analysis for the discretization of the flow equations, where spatial discretization is carried out taking into account Non Uniform Rational Basis B-Splines (NURBS) basis functions. These basis functions have advantages over traditional functions employed in the FEM (Finite Element Method). Particularly, it is easier to obtain continuity order higher than C0 between adjacent elements and geometry representation is more accurate. Pre and post-processing tools for mesh generation and results visualization are also proposed considering the data structure inherent to Isogeometric Analysis. Some classic problems of Computational Fluid Dynamics are analyzed in order to validate the proposed methodology. Results obtained here show that the present formulation has good approximation when compared with predictions obtained by reference authors. Moreover, Isogeometric Analysis is more versatile than traditional finite element formulations when spatial discretization procedures are considered.

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.

Os dispositivos oscilantes por translação de ondas correspondem a um dos sistemas de maior potencial de geração de energia elétrica a partir das ondas do oceano, razão pela qual diversos engenheiros e pesquisadores dedicam-se ao desenvolvimento e aperfeiçoamento desta tecnologia. O princípio de funcionamento destes conversores, cujos principais exemplos são o Oyster, o Waveroller e o Langlee, baseia-se no movimento oscilatório de um sistema placapistão, que promove a geração de eletricidade por meio de uma turbina. Os trabalhos existentes acerca destes dispositivos ainda não reúnem informações suficientes que possibilitem a estimativa do comportamento hidrodinâmico e da captação de energia dos mesmos em situações reais de funcionamento. Desta forma, o presente trabalho busca preencher estas lacunas por meio de modelagem numérica bidimensional, baseada em médias de Reynolds (RANS), de ondas regulares atuando sobre os respectivos dispositivos. As simulações são realizadas pelo código computacional livre e aberto OpenFOAM v 4.1 e sua extensão OLAFOAM, ambos fundamentados em volumes finitos e que utilizam a metodologia VOF para o tratamento da superfície livre, solucionando as equações de Navier-Stokes e da continuidade. A representação da dinâmica de corpo rígido dos dispositivos é tratada por uma metodologia que permite a deformação dos elementos de malha e o movimento do fundo do domínio, que oscila juntamente com as placas dos conversores no intuito de assegurar a convergência numérica. Tal metodologia é validada com base na comparação dos resultados numéricos com valores experimentais presentes na literatura. As simulações realizadas permitiram concluir que o método de fundo deformável se mostra adequado no estudo da hidrodinâmica e captação de energia destes conversores, sendo uma alternativa a outras metodologias existentes. Um aumento na largura das placas destes conversores relaciona-se a uma intensificação das velocidades angulares e numa maior captação de energia, sendo a largura ideal de placa aquela cuja dimensão é próxima a duas vezes a sua altura. Um aumento na espessura das placas ocasiona uma diminuição na oscilação e na captação energética dos dispositivos, sendo a espessura ideal com dimensão entre 10 % e 20 % da altura da placa. A captação de energia tende a diminuir com o aumento da profundidade, sendo que, para lâminas de água superiores a duas vezes a altura da placa do conversor, a influência da profundidade torna-se pouco significativa. O aumento do período de onda é responsável por intensificar a ação do campo de velocidades nas regiões próximas ao conversor, sendo o período de máxima captação de energia em torno de 11 s em escala real. Maiores alturas de onda são responsáveis por intensificar o movimento e a captação de energia dos conversores. / The oscillating wave surge converters represents one of the most effective systems to transform the energy of ocean waves into electric power, reason why many researchers and engineers are dedicated to the development and improvement of this technology. The main examples of those kind of converters are Oyster, Waveroller and Langlee, which energy production is related to a oscillating plate-piston mechanism, that drives a secondary hydraulic system, moving a turbine and generating electricity. Recent papers about this technology does not bring enough information to predict the hydrodynamic and the power capture in real operation cases. In order to fulfill those gaps, the present study is focused on bidimensional numerical modelling of regular ocean waves interaction over the converters, based on Reynolds Average Navier Stokes equations (RANS). The numerical simulations are performed by the opensource code OpenFOAM v 4.1 and its extension OLAFOAM, both based on finite volume discretization and VOF method for free surface representation, which solves the Navier-Stokes and continuity equations. The rigid body dynamics is represented by a methodology that uses the mesh morphing method combined with an oscillating bottom, which moves according to the plate movement, in order to ensure numerical convergence. The present method is validated by the comparison of the numerical model results with previous experimental studies. The numerical results led to the conclusion that the oscillating bottom method seems to be adequate to perform a study based on the hydrodynamic and power capture of the oscillating wave surge converters, representing an alternative method to modelling the dynamic of those devices. An increase in plate width are related to an increase on the device oscillation and its power capture. The ideal width is approximately twice the plate height. Thicker plates causes lower angular velocities, as a result, lower electric generation. The more efficient thickness is around 10 % to 20 % of plate height. The power capture decreases according to an increase on water depth, and the depths with values around twice of the plate height does not seem to modify the device hydrodynamics. Higher wave periods causes an intensification of the velocity field near the converter. The pitching period, which is related to the maximum power captured, is near to 11 s in real scale. The highest wave heights are responsible to increase the converter oscillation and electricity generation.

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.

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.

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.

A ênfase deste trabalho está voltada ao desenvolvimento de um novo sistema de limpeza para a colheitadeira de grãos visando, principalmente, a utilização de uma nova geometria. O design de um sistema com peneiras rotativas vem de encontro a todos os sistemas conhecidos até o momento, este sistema visa eliminar sistemas eletrônicos cada vez mais complexos e caros. Uma descrição sobre os sistemas de colheita e máquinas agrícolas, principalmente colheitadeiras, é apresentada neste trabalho, já que a bibliografia sobre este assunto é bastante restrita. A utilização de métodos computacionais, voltados à análise de fluídos foi utilizada com o intuito de nortear o desenvolvimento de protótipos reais. O método computacional dos fluídos dinâmicos simula diversas hipóteses de maneira mais rápida e mais barata até que sejam alcançados resultados suficientemente bons que compensem a fabricação de um protótipo real. Com o protótipo real os resultados virtuais podem ser validados e comparados com outros sistemas. Os resultados obtidos foram bons, a peneira rotativa tem um ótimo desempenho em testes realizados no plano, com perdas inferiores comparativamente às peneiras planas. Nos testes realizados com inclinação o desempenho não foi bom, devido a alguns problemas no duto de ar. Com relação à qualidade dos grãos separados, não foram detectadas quaisquer avarias que pudessem ter sido resultadas da peneira. / The emphasis of this work is dedicated to the development of a new cleaning system for grain harvester aiming mainly to use a new geometry in this part of the machine. The design of a system with rotary sieves comes against all systems known at this time; this system seeks to eliminate electronic systems increasingly complex and expensive. A description of the harvesting systems and agricultural machinery, harvesters mainly, is been presented in this work, since the literature on this subject is very limited. The use of computational methods, focused on the analysis of fluid was used in order to guide the development of prototypes. The method of computational fluid dynamic simulates scenarios faster and cheaper to achieve results that are good enough to compensate the fabrication of a real prototype. With the real prototype, virtual results can be validated and compared with other systems. The results achieved were good. The rotary sieve cleaning system shows excellent performance when testing on level ground with lower losses than a flat sieve system. In tests on slopes, the performance was not as good due to some problems in the air duct. With regard to grain quality, no mechanical damage that could have been produced by the screens was detected in either test.

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.