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.

Neste trabalho, são inicialmente discutidas as dificuldades referentes à obtenção de resultados numéricos para a operação de uma turbina Savonius independentes do grau de discretização, tamanho de domínio de cálculo e de máximo tempo físico simulado. Também são relatadas as divergências entre as metodologias numéricas e experimentais adotadas por diversos autores, que dificultam análises e comparações dos resultados obtidos por meio dessas metodologias com os resultados próprios obtidos. Devido a esses fatos, no presente trabalho, uma série de procedimentos experimentais e numéricos são realizados para efetuar análises do escoamento sobre uma turbina eólica Savonius. Nos experimentos em canal aerodinâmico, perfis de velocidade e parâmetros da turbulência são obtidos pela técnica de anemometria de fio quente. Medições com o uso de tubos de Pitot e manômetros eletrônicos são efetuadas para avaliar a variação da pressão e os perfis de velocidade média em posições selecionadas. Além de dados para análise, informações úteis para uso como condições de contorno nas simulações numéricas também são obtidas. Os fenômenos são reproduzidos através de simulações numéricas pelo Método de Volumes Finitos, que solucionam as equações da continuidade, de Navier-Stokes com médias de Reynolds e do modelo de turbulência k-ω SST. Análises experimentais e numéricas considerando o escoamento sobre um cilindro, que mantém semelhanças com o escoamento sobre o rotor, também são realizadas. Simulações numéricas do escoamento sobre o cilindro são efetuadas, fornecendo resultados representativos do escoamento real, quando geometrias tridimensionais são aplicadas na modelagem numérica. Nas simulações do escoamento sobre o rotor Savonius em condição estática, resultados representativos do escoamento real são obtidos com o uso de uma modelagem que leva em consideração a rugosidade das pás do rotor, estacionado na posição angular de 90°. Para posições angulares menores, não se obteve uma boa concordância entre os resultados experimentais e numéricos. A realização deste trabalho fornece informações úteis para a análise do fenômeno e tem potencial para contribuir com futuros trabalhos desse tema. / This research work initially presents a discussion about the difficulties related to obtaining numerical results for the operation of a turbine Savonius independent of the degree of discretization, calculation domain size and maximum physical time of the simulation. The differences between the numerical and experimental methodologies adopted by various authors difficult the analysis and comparisons of the results obtained through these methods with the results obtained by the methodology. Due to these facts, in this research work, a series of experimental and numerical procedures are performed to conduct analyzes of flow over a Savonius wind turbine. In the experiments on aerodynamic channel, velocity and turbulence profiles parameters are obtained by the technique of hot wire anemometry. Measurements using Pitot tubes and electronic manometers are made to evaluate the variation of pressure and mean velocity profiles at selected positions. In addition to data analysis, useful information for use as boundary conditions in the numerical simulations are also obtained. The phenomena are reproduced through numerical simulations by Finite Volume Method, that solve the equations of continuity, Reynolds-averaged Navier–Stokes equations and the equation of the turbulence model k-ω SST. Experimental and numerical analyzes considering the flow over a cylinder, which holds similarities with the flow over the rotor, are also performed. Numerical simulations of the flow over the cylinder are made, providing results representative of the actual flow when three-dimensional geometries are applied in numerical modeling. In flow simulations over the Savonius rotor in static condition, at 90°, representative results of the actual flow are obtained using a model that takes into account the roughness of the rotor blades. For smaller angular positions a good agreement between experimental and numerical results was not obtained. This work provides useful information for the analysis of the phenomenon and has the potential to contribute to future work on this theme.

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.

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 ê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.

O presente trabalho trata da modelagem de combustão em um queimador industrial de gás natural utilizando cadeias de reatores químicos idealizados. A partir de simulações de CFD e dados experimentais do ensaio denominado BERL 300 kW, sucessivas cadeias de reatores são propostas, baseando-se em premissas sobre as condições de temperatura e fluxo de massa entre as diversas regiões do escoamento reativo. O objetivo é obter uma cadeia de reatores que reproduza as emissões de NOx do ensaio. O erro relativo das emissões é de 98,23% na primeira cadeia, de 82,30% na segunda cadeia e de 1,54% na terceira cadeia. A partir da terceira cadeia de reatores proposta, considerada representativa do ensaio, realiza-se uma análise de sensibilidade e um estudo paramétrico, para avaliar a sensibilidade do modelo aos ajustes de segmentação de correntes e a capacidade de reproduzir os dados experimentais de emissão em condições de pré-aquecimento do ar de combustão e em operação fora do ponto de projeto. O modelo produz resultados satisfatórios para os casos de pré-aquecimento do ar até 110°C e para potências térmicas até 33% superiores ao ponto de projeto, não oferecendo bons resultados para as demais condições de operação. / This work shows the development and application of chemical reactor networks (CRN) to the modeling of a natural gas burner. The CRN development is based on experimental and CFD simulation results of the BERL 300 kW test. The CRN element arrangement, reactor volumes and flow splits are adjusted based on the best agreement with characteristic temperatures of the reactive flow, aiming to reproduce the experimental NOx emissions data. The relative error in the emissions predictions are 98,23% for the first CRN proposed, 82,30% for the second CRN and 1,54% for the third one. A sensitivity analysis and a parametric study of the third CRN are carried out to evaluate its sensitivity to the flow splits adjustments and its ability to predict emission with air preheat and turndown operation. The NOx emissions predicted by CRN are in good agreement until 110ºC preheating and for a turndown ratio of 1,33.