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.

Um desafio constante no projeto aerodinâmico de uma superfície é obter a forma geométrica que permite, baseado em uma determinada medida de mérito, o melhor desempenho possível. No contexto de projeto de aeronaves de transporte, o desempenho ótimo em cruzeiro é a principal meta do projetista. Nesse cenário, o uso da Dinâmica do Fluidos Computacional como não só uma ferramenta de análise mas também de síntese torna-se uma forma atrativa para melhorar o projeto de aeronaves que é uma atividade dispendiosa em termos de tempo e recursos financeiros. O método adotado para projeto aerodinâmico é baseado na teoria de controle ótimo. Essa abordagem para o problema de otimização aerodinâmica foi inicialmente proposta por Jameson (1997) e é chamada de método adjunto. Esse método apresenta uma grande diminuição de custo computacional se comparado com a abordagem de diferenças finitas para a otimização baseada em gradiente. Essa dissertação apresenta o método adjunto contínuo aplicado às equações de Euler. Tal método está inserido no contexto de um ciclo de projeto inverso aerodinâmico. Nesse ciclo, tanto o código computacional de solução das equações do escoamento quanto o código de solução das equações adjuntas foram desenvolvidos ao longo desse trabalho. Além disso, foi adotada uma metodologia de redução do gradiente da função de mérito em relação às variáveis de projeto. O algorítmo utilizado para a busca do mínimo da função de mérito é o steepest descent. Os binômios de Bernstein foram escolhidos para representar a geometria do aerofólio de acordo com a parametrização proposta por Kulfan e Bussoletti (2006). Apresenta-se um estudo dessa parametrização mostrando suas características relevantes para a otimização aerodinâmica. Os resultados apresentados estão divididos em dois grupos: validação do ciclo de projeto inverso e aplicações práticas. O primeiro grupo consiste em exercícios de projeto inverso nos quais são estabelecidas distribuições de pressão desejadas obtidas a partir de geometrias conhecidas, desta forma garante-se que tais distribuições são realizáveis. No segundo grupo, porém, as distribuições desejadas são propostas pelo projetista baseado em sua experiência e, portanto, não sendo garantida a realizabilidade dessas distribuições. Em ambos os grupos, incluem-se resultados nos regimes de escoamento transônico e subsônico incompressível. / A constant endeavor in aerodynamic design is to find the shape that yields optimum performance, according to some context-dependent measure of merit. In particular for transport aircrafts, an optimum cruise performance is usually the designers main goal. In this scenario the use of the Computational Fluid Dynamics (CFD) technique as not only an analysis tool but as a design tool becomes an attractive aid to the time and financial resource consuming activity that is aircraft design. The method adopted for aerodynamic design is based on optimal control theory. This approach to the design problem was first proposed by Jameson (1997) and it is called adjoint method. It shows a great computational cost advantage over the finite difference approach to gradient-based optimization. This dissertation presents an Euler adjoint method implemented in context of an inverse aerodynamic design loop. In this loop both the flow solver and the adjoint solver were developed during the course of this work and their formulation are presented. Further on, a gradient reduction methodology is used to obtain the gradient of the cost function with respect to the design variables. The method chosen to drive the cost function to its minimum is the steepest descent. Bernstein binomials were chosen to represent the airfoil geometry as proposed by Kulfan and Bussoletti (2006). A study of such geometric representation method is carried on showing its relevant properties for aerodynamic optimization. Results are presented in two groups: inverse design loop validation and practical application. The first group consists of inverse design exercises in which the target pressure distribution is from a known geometry, this way such distribution is guaranteed to be realizable. On the second group however, the target distribution is proposed based on the designers knowledge and its not necessarily realizable. In both groups the results include transonic and subsonic incompressible conditions.

O aumento do uso de códigos numéricos baseados em dinâmica de fluidos computacional (CFD) atrelado a problemas de engenharia naval e a necessidade de se estabelecer a proximidade dos resultados obtidos nas simulações com os experimentos físicos serviu de motivação para execução deste trabalho. O tema do presente trabalho é pertinente uma vez que o número de trabalhos publicados que simulam com precisão ensaios de embarcações em condição de manobra e obtêm as derivadas hidrodinâmicas que caracterizam o movimento é muito pequeno, atendendo assim a uma demanda existente por novas confrontações numérico/experimentais. Com isso, um conjunto de ensaios de PMM (Planar Motion Mechanism), mais especificamente os ensaios de sway e yaw puros e yaw com drift, e arrasto oblíquo foram executados na embarcação KVLCC2, nas condições de manobra, similares aos ensaios realizados pelo laboratório italiano ISEAN que tinha por objetivo fornecer uma base de comparação para trabalhos em CFD. Os resultados das forças e momentos das simulações foram comparados com os obtidos nos ensaios, bem como as derivadas hidrodinâmicas que caracterizam a embarcação estudada. De maneira geral, os resultados obtidos nas simulações mostraram-se muito próximos aos obtidos a partir do ensaio, fornecendo assim um bom resultado nas derivadas hidrodinâmica, principalmente com relação às derivadas em função das velocidades de translação e rotação. Na expectativa de ampliar a aplicação do método de solução proposto, um novo conjunto de simulações foi realizado na escala real da embarcação. As dificuldades encontradas devido às incompatibilidades nos números de Reynolds mostraram-se presentes, de maneira similar às encontradas nos ensaios físicos evidenciando inclusive que o fator de forma está diretamente atrelado ao número de Reynolds, contrariando a hipótese clássica e alinhando-se com trabalhos mais recentes publicados sobre o tema. Os resultados encontrados a partir da aplicação da metodologia mostram coerência em relação às respostas encontradas nos ensaios, fornecendo uma alternativa ao método clássico de obtenção das características hidrodinâmicas de uma embarcação de superfície em condição de manobra. Ainda fornece uma opção de cálculos, que ainda necessita de melhorias e validação mais elaborada, nas condições da escala real, fato impossível sem o uso de simulações. / The growing use of numerical codes based on computational fluid dynamics (CFD) in problems related with marine engineering and the need to establish the reliability of the results from the simulations with those derived from physical tests served as the motivation for performing the work. The subject of this work is relevant since the number of published works that accurately simulate and test the conditions of vessels in maneuvering and obtains the hydrodynamic derivatives that characterize the movement is very small, thus meeting an existing demand for the new numeric confrontations / experimental. Therefore, a set of PMM tests (Planar Motion Mechanism), more specifically: pure sway, pure yaw and yaw whit drift, and static drift, were performed on the vessel KVLCC2 under conditions of maneuver, similar to tests conducted by the Italian laboratory ISEAN which was intended to provide a basis of comparison for work with CFD. The results of the forces and moments of simulations were compared with those from the tests, as well as the hydrodynamic derivatives which characterize the vessel studied. In general, the results from simulations were very close to those from the experiment, thus providing a good result in hydrodynamic, derived primarily related to derivatives, depending on the speed of translation and rotation. Expecting to expand the application of the proposed method, a new solution set of simulations was carried out in a full-scale vessel. The difficulties encountered due to incompatibilities in the Reynolds numbers were similarly found in the physical testing. This included demonstrating the form factor directly linked to the Reynolds number, contrary to the classical hypothesis and aligning with the latest publications on the subject. The results from the application of the methodology show consistency with the solutions found in the tests, providing an alternative to the classical method of obtaining the hydrodynamic characteristics of a surface vessel in a condition of maneuver. We also provide an option of calculations, which still needs improvement and a more elaborate validation under full-scale conditions, which would be impossible without the use of simulations.

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.

A evolução do uso da energia eólica nas últimas décadas está diretamente relacionada ao desenvolvimento da tecnologia empregada na conversão e projeto das instalações. A viabilização de uma instalação eólica de grande porte depende da avaliação correta do potencial eólico. A fim de avaliar a capacidade de conversão de energia cinética do vento em torque, é fundamental a modelagem da esteira aerodinâmica das turbinas eólicas. Este trabalho apresenta um estudo sobre a esteira aerodinâmica formada por uma turbina eólica dimensionada conforme a teoria de otimização de Betz. Para tanto, realiza-se inicialmente uma pesquisa sobre a evolução da transformação da energia contida no vento em energia mecânica e métodos de análise adotados. Para modelagem da esteira, realiza-se a simulação numérica do escoamento sobre uma turbina de eixo horizontal empregando o Método de Volumes Finitos. Através do uso da Dinâmica dos Fluidos Computacional são resolvidas as Equações de Navier-Stokes com Médias de Reynolds (RANS) e a utilização dos modelos de turbulência k-, k- RNG, k- e k- SST. Para a solução das equações é utilizado o programa ANSYS-CFX. Define-se o perfil NACA 4412 como perfil aerodinâmico para projeto das pás da turbina, e modela-se a turbina através da teoria de dimensionamento ótimo de Betz. O domínio, discretizado em um número finito de volumes de controle, possui duas regiões distintas, uma estática e outra rotacional, representado o rotor da turbina. Inicialmente são apresentadas simulações com os diferentes modelos de turbulência e definido o modelo que apresenta os melhores resultados, o k- SST. Para o modelo escolhido são realizadas simulações incluindo estudos com a turbina inclinada em relação à direção de incidência do vento, para verificar a alteração no perfil da esteira gerada e na capacidade da turbina em transformar a energia do vento em Torque. Os resultados obtidos para os campos de velocidade e pressão são comparados com os de outros autores e mostraram-se coerentes, indicando que a simulação feita é capaz de representar o fenômeno analisado. / The evolution of the use of the wind energy in recent decades is directly related to the technology in the facilities conversion and design. The feasibility of a large wind farm depends on the correct assessment of wind potential. To evaluate the capacity of converting wind kinetic energy in torque, it is essential to model the wake aerodynamics of wind turbines. This paper presents a study on the aerodynamic wake formed by a wind turbine sized according to the Betz optimization theory. For this, initially it is performed a research on the evolution of the energy contained in wind into mechanical energy and the adopted analysis methods. For wake modeling, a horizontal axis wind turbine flow numerical simulation is done using the Finite Volume Method. Using the Computational Fluid Dynamics, the Reynolds Averaged Navier-Stokes equations are solved with the use of the k-, k- RNG, k- e k- SST turbulence models and using the software ANSYS-CFX. It is defined the NACA 4412 profile as aerodynamic profile for turbine blades and the turbine is modeled using the Betz optimization theory. The domain, discretized into a finite number of control volumes, has two distinct regions, one static and one rotational, represented the turbine rotor. Initially the simulations with different turbulence models are presented and the model k- SST is defined as the model that gives best results. For the chosen model, simulations are performed, including studies on the turbine yawed with respect to the wind incidence direction, to verify the profile change in the generated wake and turbine capacity to convert wind energy into torque. The obtained results for the velocity and pressure fields are compared to other authors are very consistent, indicating that the proposed simulation is capable of representing the analyzed phenomenon.

O presente trabalho apresenta o estudo da esteira de turbinas eólicas de eixo horizontal, utilizando um modelo de turbina eólica testada em túnel de vento. Trata-se de uma turbina eólica de duas pás com 10 m de diâmetro, conhecida como UAE Phase VI, testada pelo NREL no túnel de vento NASA Ames. Desta turbina foi obtida a geometria e dados experimentais, como o coeficiente de empuxo, além das características do túnel de vento. Utilizando este experimento como modelo, são realizadas simulações com equações médias de Reynolds e modelo de turbulência k-ω SST, com o uso de malha móvel e abordagens permanente e transiente, com o software de volumes finitos Star-CCM+®. São empregadas malhas poliédricas em um domínio de grande extensão na esteira, para apenas uma velocidade do vento. São apresentados três modelos de esteira, dois analíticos e outro considerado numérico, para fins de comparação com os resultados simulados. Estes modelos de esteira tem o objetivo de descrever o campo de velocidade a jusante da turbina. Os dados de entrada dos modelos são os mesmos, retirados da simulação. Além dos modelos para o déficit de velocidade, são apresentados equacionamentos para se obter a intensidade de turbulência na esteira, com bons resultados. No modelo analítico de esteira PARK, a constante de decaimento é baixa, pois o caso estudado possui escoamento sem grandes perturbações. Verifica-se que esta constante deve variar ao longo da esteira para se ajustar ao campo médio de velocidades simulado. Outro modelo analítico proposto por Werle foi estudado, não apresentando boa concordância com a simulação, pois apresenta um rápido incremento da velocidade na linha central na esteira distante. Este modelo se ajustou bem com o campo não perturbado a montante da turbina. O modelo de viscosidade turbulenta é um modelo de esteira com equacionamento original de segunda ordem para obter um campo de velocidade. Contudo, é apresentado e calculado com uma simplificação para primeira ordem, obtendo-se somente a velocidade na linha de centro. Este modelo tem o melhor ajuste com a simulação. A metodologia adotada para simulação numérica do escoamento apresenta boa concordância com os resultados experimentais de CT, já os dados anemométricos medidos na esteira próxima não foram reproduzidos com a mesma qualidade, contudo aproximam-se de forma qualitativa. A comparação com os modelos de esteira da literatura confirma que o modelo de viscosidade turbulenta é o que mais se aproxima dos resultados da simulação. / This study presents wakes of horizontal axis wind turbines, using as model a wind turbine tested in wind tunnel. It is a two bladed wind turbine with 10 m diameter, known as UAE Phase VI, tested by NREL in the NASA Ames wind tunnel. From this experiment was obtained the geometry and measured data, as the thrust coefficient, and the characteristics of the wind tunnel. Using this experiment as model, simulations are performed with Reynolds Averaged Navier-Stokes equations and k-ω SST turbulence model, using moving mesh and permanent and transient approaches, with help of the finite volume software Star-CCM+®. Polyhedral meshes are employed in a domain with large extent downwind, for just one wind speed. There are presented three wake models, two analytical and other considered numerical, for purposes of comparison with simulated results. These wake models aims to describe the velocity field downstream of the turbine. The input data of the models are the outputs of the simulations. Besides the models for the velocity deficit, are given equations to derive the turbulence intensity on the wake, with good results. In the analytical PARK wake model, the decay constant is low, as the study case has low perturbed flow. It is noted that this constant should vary along the wake to fit the average speeds field simulated. Another analytical model proposed by Werle is presented, not showing good results when comparing with the simulation, since it shows a rapid increase in the centerline velocity on the far wake. This model fits well with the undisturbed field upwind of the turbine. The turbulent viscosity is a wake model with the original second order equation to obtain the velocity field. However, it is presented and calculated in a simplified way to first order, obtaining only the centerline velocity. This model has the best fit with the simulation. The methodology adopted for numerical simulation of the flow shows good agreement with the experimental results of CT, since the wind data measured on the near wake are not reproduced with the same quality, yet close to a qualitative manner. The comparison with the wake models from the literature confirms that the eddy viscosity wake model is closest to the simulation results.

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.