Nos últimos anos, a utilização da energia eólica vem apresentando uma tendência de aumento. 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 com terrenos complexos 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 uma superfície de topografia complexa. Resultados numéricos com diferentes alternativas de modelagem do problema são comparados com dados de um experimento em túnel de vento, visando determinar a metodologia adequada para avaliação do problema proposto. As simulações numéricas do escoamento de ar sobre o terreno 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). O estudo é dividido em três casos. No primeiro caso, a rugosidade superficial é negligenciada e o problema de fechamento é contornado com a utilização do modelo de turbulência k ω SST. No segundo caso, a rugosidade superficial é estipulada de acordo com modelo utilizado no ensaio experimental e o problema de fechamento é contornado com a utilização do modelo de turbulência k - ε. No terceiro caso, a rugosidade superficial também é estipulada de acordo com modelo utilizado no ensaio experimental e o problema de fechamento é contornado com a utilização do modelo de turbulência k ω SST. Os resultados das simulações são apresentados de forma que se possa observar o perfil de velocidades adimensional sobre a superfície da geometria para cada caso, para que seja possível verificar o campo de velocidades sobre a superfície em estudo. O resultado das simulações são comparados com dados experimentais obtidos em túnel de vento: verifica-se um comportamento similar nos perfis de velocidade alcançados. Através da análise do campo de velocidades sobre a superfície em estudo, pode-se obter a localização que apresenta o melhor potencial eólico de uma região. Este processo é conhecido como Micrositing. / In recent years the use of wind energy has shown an increasing. A key aspect to determine the technical and economic viability for the wind power plant is the accurate assessment of the distribution of wind speeds in the area of utilization. The installation of wind turbines in areas with complex terrain has determined the necessity of improve the methodology for the prediction of wind velocity field in order to better determine the distribution of equipment and utilization of existing potential. 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 a model of a complex surface topography. Numerical results with different alternatives for modeling the problem are compared with data from an experiment in wind tunnel to determine the appropriate methodology for evaluation of the problem. The numerical simulations of the air flow over the terrain are performed using the ANSYS Fluent 13.0 which uses the finite volumes method for solving the Reynolds Averaged Navier Stokes (RANS) equation. The study is divided in three cases. In the first one, the surface roughness is neglected and the closure problem is solved by k ω SST turbulence model. In the second case, the surface roughness is stipulated according to the model used in the experimental test and the closure problem is solved by k ε turbulence model. In the third case, the surface roughness is stipulated according to the model used in the experimental test and the closure problem is solved by k ω SST turbulence model. The simulation results are presented so that one can observe the dimensionless velocity profile on the surface in each case in order to check the velocity field on the surface under investigation. These results are compared with experimental data obtained in wind tunnel which shows consistency with them. By analyzing the velocity field on the surface, it can be verified that the exact location where there is the best wind potential of a region. This process is called Micrositing.

Uma das mais recentes abordagens propostas na literatura para tratar o problema do alto número de Weissenberg (We) é a Formulação Log-Conformação (FLC). Nesta formulação, a equação constitutiva viscoelástica utilizada é reescrita em termos de uma variável Ψ, que é o logaritmo do tensor conformação. Apesar do potencial de aplicação da FLC, pouca atenção tem sido dirigida para análise da acurácia da solução obtida para o campo de tensões quando se utiliza esta formulação. Assim, o objetivo do presente trabalho foi estudar a acurácia da solução obtida pela FLC na análise de escoamentos de fluidos viscoelásticos usando duas geometrias padrão de estudo: placas paralelas e cavidade quadrada com tampa móvel. Primeiramente, a FLC foi implementada no pacote de CFD OpenFOAM. Em seguida foram verificados os limites do número de Weissenberg na formulação numérica padrão (Welim,P), onde para a geometria de placas paralelas foi encontrado Welim,P = 0,3 e para a geometria da cavidade quadrada com tampa móvel foi encontrado Welim,P = 0,8. Depois o código implementado foi aplicado em ambas as geometrias, comparando-se a solução obtida pela FLC com aquela da formulação padrão na faixa de We < Welim,P. Os resultados obtidos na geometria de placas paralelas apresentaram boa concordância com a solução padrão e solução analítica. Para a geometria da cavidade quadrada com tampa móvel, que não possui solução analítica, boa concordância dos resultados também foi observada em comparação com a solução padrão. Posteriormente foram comparados os resultados obtidos pela FLC na faixa de We > Welim,P. Na geometria de placas paralelas, além da boa concordância com a solução analítica, obteve-se convergência em todos os casos estudados neste trabalho, com o maior número de Weissenberg utilizado sendo igual a 8 Os resultados da geometria da cavidade quadrada com tampa móvel também apresentaram boa concordância em comparação com dados da literatura, porém a convergência foi obtida até para We = 2. Com respeito à comparação das formulações numéricas com a solução analítica, feita apenas na geometria de placas paralelas, foi observado um erro máximo de 7,57% na solução padrão e de 12,33% na FLC. Em relação à análise da qualidade das tensões usando os resíduos da equação constitutiva viscoelástica como critério de acurácia, foi verificado nas duas geometrias que os valores de tensão obtidos usando a FLC são menos acurados que aqueles obtidos pela formulação explícita no tensor das tensões nos casos em que esta última converge. Também foi observado que a acurácia diminui com o aumento do We. Esse efeito pôde ser melhor notado na geometria de placas paralelas. Uma razão para a perda de acurácia da tensão provavelmente ocorre devido à natureza matemática da transformação algébrica inversa de Ψxx para τxx. O novo solver implementado neste trabalho apresentou convergência e soluções corretas para as duas geometrias, logo foi implementado corretamente. Ele também potencializa o solver de partida viscoelastiFluidFoam ao estender simulações para uma faixa maior do número de Weissenberg. / A recent approach proposed in the literature to deal with the High Weissenberg Number Problem is the Log-Conformation formulation (LCF). In this formulation the viscoelastic constitutive equation is rewritten in terms of the logarithm of the conformation tensor Ψ. Despite the great potential application of the LCF, little attention has been given in the literature to the accuracy of the obtained stress fields. The purpose of this work was to study the solution obtained by LCF in the analysis of viscoelastic flows using two benchmark geometries: parallel plates and lid driven cavity. Firstly, the LCF was implemented in the OpenFOAM CFD package. Then, the limits of Weissenberg number for the standard numerical formulation (Welim,P) were verified, obtaining Welim,P = 0.3 for the parallel plates and Welim,P = 0.8 for the lid driven cavity. When comparing the solution obtained by the LCF with that of the standard formulation in a range of We < Welim,P, the results obtained for the parallel plates geometry showed good agreement with the standard solution and the analytical solution. For the lid driven cavity geometry, for which there is not analytical solution, good agreement with the standard solution was also observed. For We > Welim,P in the parallel plates geometry, in addition to the good agreement with the analytical solution, it was possible to obtain convergence in all the cases studied in this work, with the largest number of Weissenberg used being equal to 8 The results of the lid driven cavity geometry also presented good agreement in comparison with literature data, but convergence was obtained up to We = 2. With respect to the comparison of the numerical formulations with the analytical solution for the parallel plates geometry, a maximum error of 7.57% was observed in the standard solution and of 12.33% in the LCF. When using the residues of the viscoelastic constitutive equation as a criterion of accuracy, it was verified that for the two geometries the stress values obtained using the LCF were less accurate than those obtained by the explicit formulation in the stress tensor. It has also been observed that accuracy decreases with increasing of We. One reason for the loss of stress accuracy probably occurs because of the mathematical nature of the inverse algebraic transformation from Ψxx to τxx. The new solver implemented in this work presented convergence and correct solutions for the two geometries, so it was implemented correctly. It also potentiates the viscoelastiFluidFoam starting solver by extending simulations to a larger range of Weissenberg number.

O mercado de biorreatores com membranas (BRMs) têm exibido alto crescimento. Contudo, o fouling diminui o desempenho desses sistemas drasticamente. A aeração promove a mitigação do fouling, mas possui alto custo operacional (de até 70% do total da demanda energética) e é utilizada de forma otimizada apenas 10% das vezes, gerando a necessidade de estudos sobre a hidrodinâmica em BRMs. Ferramentas de dinâmica de fluidos computacional (CFD) são úteis para esse tipo de análise. Um dos primeiros passos para encontrar uma solução apropriada em simulações numéricas é propor uma modelagem correta. Dentre os fenômenos a serem modelados, os efeitos da turbulência são provavelmente um dos mais importantes; porém, nenhum artigo que comparasse predições com base na escolha de modelo de turbulência para simulações de sistemas BRM foi encontrado. Dessa forma, o objetivo desse trabalho foi verificar a importância da escolha do modelo de turbulência para simulações de biorreatores com membranas através de CFD. Para isso, as predições obtidas de velocidade local próxima às superfícies das membranas e de tensão cisalhante nessas superfícies para duas taxas de aeração, 5 e 15 m³/h, empregando dois modelos de turbulência, k-ϵ com funções de parede para alto (aR) e para baixo número de Reynolds (bR) e k-ω SST (Shear Stress Transport) com funções de parede automáticas, na análise de um sistema BRM foram comparadas entre si e com dados experimentais e simulados disponíveis na literatura. Os perfis temporais da velocidade e da tensão cisalhante exibiram alta variabilidade no período das oscilações em certos pontos, exigindo um longo tempo de escoamento para a convergência das variáveis locais. Identificou-se a necessidade de outorgar maior importância à definição do intervalo de tempo de coleta de dados experimentais, de modo a adquirir médias representativas do perfil dinâmico das variáveis e destes perfis serem também analisados para comparações mais definitivas entre resultados de simulações e medições experimentais. As diferenças, entre as medições experimentais da literatura e predições, obtidas pelas simulações deste trabalho foram, no geral, de ordem similar ou menores que as obtidas pelas simulações na literatura. Além disso, maior atenção deve ser dada à escolha da estratégia de modelagem de turbulência, visto que houve alta sensibilidade das predições, que variaram em até 21,6% dependendo dessa escolha. / Membrane bioreactors (MBR) market has been showing high growth rates over recent years. However, membrane fouling drastically reduces MBR overall performance. Aeration promotes fouling mitigation, but at a high operational cost (up to 70% of the MBR energy demand) and it is optimally employed only in 10% of the cases. This created the need of studies focused on MBR hydrodynamic. Computational fluid dynamics (CFD) is a useful tool for hydrodynamic analysis. One of the first steps in finding a proper solution for numerical simulation is proposing a correct modelling. Among the phenomena to be modelled for MBR simulations, turbulence effects are probably one of the most important; nevertheless, no paper comparing the predictions based on the turbulence model choice for MBR simulations was found. In sight of that, this work aimed to verify the relevance of the choice of turbulence model for MBR simulations through CFD. Predictions of local velocities near membranes surfaces and of local shear stress on those surfaces, for two aeration rates (5 and 15 m³/h), employing k-ϵ with wall functions for high (aR) and low (bR) Reynolds number and k-ω SST with automatic wall functions, on the analysis of a MRB system, were compared between each other and with experimental and simulated data available in the literature. The velocity and shear stress temporal profiles showed oscillations with highly variable periods in some points, which required a long process real time to verify the local variables convergence. It was identified the need to give more importance to the definition of the time interval for experimental data collection in order to acquire reliable temporal means; also, one must properly analyze the temporal profiles for more definitive comparisons between predictions and experimental measurements. The differences, between experimental data and predictions, obtained through this work simulations were, in general, of similar order or smaller than the ones reported in the literature. Besides, more attention must be given to the turbulence modelling choices, since the predictions obtained here were highly sensitive to those choices, showing differences up to 21,6% among them.

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

O presente trabalho apresenta um estudo numérico sobre o desempenho de trocadores de calor com microcanais. Primeiramente, para a verificação da ferramenta numérica, uma comparação com um estudo experimental conhecido é realizado. Em seguida, aperfeiçoamentos disponíveis na literatura são testados e comparados. As simulações numéricas são realizadas com o software Fluent utilizando o Método dos Volumes Finitos (MVF). O resultado numérico da resistência térmica do dissipador do estudo experimental em que a primeira parte desse trabalho foi baseada ficou em 0,097 °C/W (com uma temperatura máxima de 373 K), enquanto o resultado experimental ficou em 0,090 °C/W, o que representa uma diferença de 7,2 %. São realizados aperfeiçoamentos integrando o trocador de calor a um sistema de refrigeração. Isso reduz a temperatura de entrada e, consequentemente, a temperatura máxima encontrada no sistema (350 K para uma mesma potência dissipada) não havendo mudança na resistência térmica, porém trazendo um ganho em confiabilidade. São apresentadas mudanças na geometria do dissipador diminuindo a área de entrada, com uma sensível piora nos resultados da resistência térmica, fazendo com que a resistência aumentasse de 0,097 °C/W para o canal original para 0,272 °C/W para o canal com menor altura. Apresenta-se um trocador de calor cujos microcanais em suas dimensões longitudinais possuem ondulações, porém os efeitos esperados de melhor mistura e maior área de contato para o fluxo de calor não foram observados para as condições de contorno do experimento original; a resistência térmica calculada foi de 0,102 °C/W. A última análise apresenta uma idealização onde o fluxo de calor é distribuído uniformemente no volume de silício, encontrando-se uma resistência térmica de 0,084 °C/W. / This paper presents a numerical study on the performance of heatsinks with microchannels. For the validation of the numerical tool, a comparison with a known experimental study is carried out, then improvements available in the literature on the subject are tested and compared. The numerical simulations are performed with the Fluent software using the Finite Volumes Method (MVF). The numerical result of the thermal resistance of the heatsink from the experimental study which the first part of this paper is based was 0.097 °C/W, whereas the experimental result was 0.090 °C/W (with a maximum temperature of 373 K), representing a difference of 7.2%. An improvement is performed by integrating the heatsink to a HVAC (Heating, ventilation and air conditioning) system, which lowers the inlet temperature and consequently lowers the highest temperature found in the system (350 K for the same dissipated power) without changing the thermal resistance, but providing gains in reliability. Changes in the heatsink geometry are presented by decreasing the inlet area, with a significant worsening in the thermal resistance results, causing the resistance to vary from 0.097 °C/W for the original channel to 0.272 °C/W for the channel that has the lowest height. A heatsink whose microchannels have waves in their longitudinal dimensions is shown, but the expected effects of better mixing and greater contact area for the heat flux were not observed for the boundary conditions of the original experiment, the calculated thermal resistance was of 0.102 °C/W. The last analysis presents an idealization where the heat flux is evenly distributed in the silicon volume, with a thermal resistance of 0.084 °C/W.

A simulação numérica de escoamentos reativos, como a combustão, tem um caráter altamente não-linear devido a presença de diversas reações químicas que acontecem entre as espécies que descrevem o processo de oxidação do combustível. Além disso, tais processos ocorrem a nível molecular, tornando o sistema de equações governantes rígido, o que implica na necessidade de esquemas numéricos de alta ordem bem como malhas finas e passo de tempo pequeno, aumentando consideravelmente o custo computacional. Neste sentido, o uso de mecanismos de oxidação detalhados na simulação numérica é proibitivo, e técnicas de redução química são necessárias de modo a desenvolver modelos reduzidos com menos variáveis e rigidez moderado, mantendo a precisão e abrangência do modelo detalhado. O objetivo do presente trabalho é obter uma comparação dos resultados obtidos para duas técnicas de redução química diferentes, Directed Relation Graph - DRG, baseada no desenvolvimento de mecanismos esqueletos, e a Reaction Diffusion Manifolds - REDIM, baseada na separação das escalas de tempo. Como validação dos modelos propostos, simulações numéricas 1D de chamas pré-misturadas e não pré-misturadas, bem como de reatores homogêneos, são desenvolvidas. Além disso, uma estratégia que une as duas técnicas de redução é apresentada, com o objetivo de ser aplicada em mecanismos cinéticos grandes. / Numerical simulation of reactive flows, such as combustion, has a highly non-linear character due to the presence of several chemical reactions that occur among the chemical species that describe the process of fuel’s oxidation. Besides, such processes occur at a molecular level, making the system of governing equations stiff, which implies in the need of high order numerical schemes as well as fine meshes and small time step, enhancing considerably the computational cost. In this sense, the use of detailed oxidation mechanisms in the numerical simulation is prohibitive, and chemical reduction techniques are needed in order to develop reduced models with less variables and moderate stiffness, while maintaining the accuracy and comprehensiveness of the detailed model. The objective of the present works if to obtain a comparison between two chemical reduction techniques, the Directed Relation Graph - DRG, based on the skeletal mechanisms generation, and the Reaction Diffusion Manifolds - REDIM, based on the separation of time scales. As validation of the proposed models, one-dimensional numerical simulations of premixed and non-premixed flames, as well as homogeneous reactors, are carry out. Besides, a coupled methodology between DRG and REDIM is presented, that will provide a useful tool for simulation of fuels with very large detailed kinetic mechanisms.

Na matriz elétrica brasileira, as termelétricas vêm sendo utilizadas como fontes de energia complementar. Para minimizar os impactos ambientais na geração de energia a partir do carvão mineral, diversos estudos apontam para o uso da tecnologia de leito fluidizado como uma alternativa flexível, eficiente e limpa. Neste tipo de reator a fluidodinâmica determina a distribuição das partículas e do gás, o que influencia na mistura entre as fases, e consequentemente, na eficiência do processo, porém, não é totalmente compreendida. Baseando-se nisto, esta tese tem como objetivo caracterizar a fluidodinâmica do leito fluidizado borbulhante, a partir de modelos matemáticos e análise experimental. Busca-se desenvolver uma maior compreensão da influência fluidodinâmica no desempenho da planta em processos de combustão do carvão mineral brasileiro, o qual possui alto teor de cinzas. Para isso, os modelos matemáticos foram validados na planta de bancada a temperatura ambiente, e posteriormente, aplicados para a planta piloto a 850 °C. O desempenho de cada operação foi avaliado a partir do rendimento de conversão e teor de carbono residual nas cinzas para diferentes alturas de leito e velocidade superficial do gás. Os resultados obtidos mostram que para a obtenção de um maior tempo de residência, a velocidade superficial do gás mais adequada depende da altura do leito fixo. O teor de carbono residual nas cinzas foi menor para o leito mais alto, o que é explicado pelas características do combustível utilizado e pelo ponto de alimentação do carvão. Não foi observada relação direta entre os coeficientes de troca gasosa para a transferência de massa no leito, com o teor de carbono residual nas cinzas, ou com o rendimento de conversão. Para uma análise mais completa da influência fluidodinâmica no desempenho da planta, a região de freeboard deve ser incluída na modelagem. A análise em planta de bancada a temperatura ambiente foi fundamental para compreender como ocorre a movimentação das partículas e das bolhas. As variações na altura do leito e velocidade do gás geraram diferentes variações nos coeficientes de troca gasosa para a transferência de massa, no tempo de residência e no teor de carbono nas cinzas. / In the Brazilian electric matrix, thermoelectric power plants are used as complementary energy sources. In order to minimize the environmental impacts of coal power generation, several studies point to the use of fluidized bed technology as a flexible, efficient and clean alternative. In this type of reactor, the fluid dynamics determines the distribution of particles and the gas, which influences in the mixing between the phases, and consequently, in the efficiency of the process, however, is not fully understood. Based on this information, this thesis aims to characterize the fluid dynamics of the bubbling fluidized bed, from mathematical models and experimental analysis. Aim to develop a better understanding of the fluid dynamics influence on the plant performance in the combustion processes of the Brazilian coal, which has a high ash content. For this purpose, the mathematical models were validated in the bench plant at room temperature, and later, applied to the pilot plant at 850 °C. The performance of each operation was evaluated from the combustion efficiency and residual carbon content in the ashes for different bed heights and gassurface velocity. The obtained results show that to attain a longer residence time the most suitable gas-surface velocity depends on the height of the fixed bed. The residual carbon content in the ashes was lower in the upper bed, which is explained by the characteristics of the fuel used and the point of coal feeding. Besides, was not observed a direct relation between the gas exchange coefficients for mass transfer in the bed with the residual carbon content in the ashes or with the combustion efficiency. Toward a more complete analysis of fluid dynamics influence on the plant performance, the freeboard region should be included in the modeling. The analysis in bench plant at room temperature was fundamental to understand how the movement of particles and bubbles occurs. The variations in bed height and gas velocity generated different variations in the gas exchange coefficients for mass transfer, residence time and carbon content in the ashes.

O presente trabalho tem por objetivo estudar as Interações Turbulência-Radiação em um escoamento não reativo para meios participantes. Estas interações caracterizam-se por um complexo fenômeno transiente, devido à combinação de dois fenômenos, unindo as características das flutuações da turbulência e da elevada não linearidade do fenômeno da radiação térmica. O estudo consiste em análise numérica do problema por dinâmica de fluidos computacional, através da utilização do Fire Dynamics Simulator (FDS), um software Open-Source, na qual a modelagem da turbulência é feita através da Simulação de Grandes Escalas. Como se trata de um software novo, bem como sendo introduzido no grupo de pesquisa, primeiramente é realizada a simulação de um caso benchmark para verificação e avaliação da formulação numérica. A análise do TRI é realizada em um problema proposto baseado em trocadores de calor reais utilizados em máquinas térmicas, como por exemplo, geradores de vapor ou coletores de escapamento de motores, envolvendo transferência combinada de convecção forçada e radiação térmica. A metodologia de avaliação consiste em comparar o fluxo radiante médio nas fronteiras obtido através da simulação transiente e compará-lo com o fluxo obtido por meio do campo médio temporal de temperaturas. São avaliadas a influência da intensidade de turbulência na entrada do escoamento, assim como a da espessura óptica, ambos relevantes para os efeitos do TRI. Conforme descrito pela literatura, neste tipo de problema as interações podem ser negligenciadas, confirmando os resultados obtidos, da ordem de 2% para o fluxo radiante. / This dissertation has the objective of analyzing the Turbulence-Radiation Interaction for a non-reactive flow with a participating media. These interactions are characterized by complex transient effects, due to the combination of two phenomena, coupling the scalar fluctuations of the turbulence and the highly non-linearity of thermal radiation. The study consists in a numerical analysis through Computational Fluid Dynamics, using the Fire Dynamics Simulator (FDS), an Open-Source software, which employs the Large Eddy Simulation method. Because the software is under development and new in the research group, it will be performed the simulation of a benchmark case for verification and evaluation of the numerical methodology. The TRI analysis will be performed in a proposed problem, based on real heat exchangers, as an example, steam generators or exhaust manifold of combustion engines, involving combined heat transfer between forced convection and radiative heat transfer. The methodology consists in evaluating the radiative mean heat flux obtained by the transient simulation and compare it with the flux obtained with the time-averaged temperature field. It will be evaluated the influence of the turbulence intensity at the inlet and the optical thickness, both very important for the TRI effects. According to the literature, in this case the TRI effects could be neglected, confirming the obtained results, around 2% for the radiative heat flux.

Dissertação (mestrado)—Universidade de Brasília, Programa de Pós-Graduação em Arquitetura e Urbanismo, 2012. / Submitted by Alaíde Gonçalves dos Santos (alaide@unb.br) on 2012-05-31T12:31:51Z No. of bitstreams: 1 2012_GustavodeLunaSales.pdf: 12763978 bytes, checksum: 2726d5a3fd10ac2247cf90d06e13dacb (MD5) / Approved for entry into archive by Elzi Bittencourt(elzi@bce.unb.br) on 2012-06-02T13:32:32Z (GMT) No. of bitstreams: 1 2012_GustavodeLunaSales.pdf: 12763978 bytes, checksum: 2726d5a3fd10ac2247cf90d06e13dacb (MD5) / Made available in DSpace on 2012-06-02T13:32:32Z (GMT). No. of bitstreams: 1 2012_GustavodeLunaSales.pdf: 12763978 bytes, checksum: 2726d5a3fd10ac2247cf90d06e13dacb (MD5) / Esta dissertação aborda um estudo de ventilação natural e a identificação da influência do tamanho de aberturas sobre a velocidade e distribuição do fluxo de ar no edifício da Secretaria de Cultura da cidade de Bragança, no estado do Pará, escolhida como estudo de caso. Os métodos utilizados para alcançar os objetivos estabelecidos foram: revisão bibliográfica sobre o tema; coleta de dados e medições no local dos aspectos relativos à ventilação natural; informações coletadas por meio de questionários aplicados aos funcionários da edificação e anotações do pesquisador; simulações computacionais no programa de dinâmica dos fluidos CFX do pacote ANSYS/Worbench 2012. Os resultados apontaram aspectos positivos da escala urbana e da edificação no sentido de favorecer o aproveitamento da ventilação natural no estudo de caso. Com relação aos modelos e cenários propostos, verificou-se que o aumento de 20% no tamanho das aberturas pode possibilitar uma melhor distribuição e níveis mais elevados de velocidade do ar no interior de edificação. _______________________________________________________________________________ ABSTRACT / This dissertation presents a study on natural ventilation and to identify the influence of size openings on the speed and air distribution in the Secretariat of Culture building`s in the city of Bragança in Brazil. The following methods were used: literature review; data collection; on-site measurements; questionnaires for users; researched notes and computer simulations in fluid dynamics software ANSYS/CFX. The results indicated positive aspects in urban and building scale favoring the use of natural ventilation in the case study. The models and scenarios indicated that the 20% increase in the openings improves the distribution and increases the air speed in the building.

Tese (doutorado)—Universidade de Brasília, Faculdade de Tecnologia, Programa de Pós-Graduação em Ciências Mecânicas, 2016. / Submitted by Albânia Cézar de Melo (albania@bce.unb.br) on 2016-04-06T12:53:32Z No. of bitstreams: 1 2016_BraulioGutierrezPimenta.pdf: 58058787 bytes, checksum: a0cdeeb6fca3ecffb891becf62a09da5 (MD5) / Approved for entry into archive by Raquel Viana(raquelviana@bce.unb.br) on 2016-04-12T20:30:53Z (GMT) No. of bitstreams: 1 2016_BraulioGutierrezPimenta.pdf: 58058787 bytes, checksum: a0cdeeb6fca3ecffb891becf62a09da5 (MD5) / Made available in DSpace on 2016-04-12T20:30:54Z (GMT). No. of bitstreams: 1 2016_BraulioGutierrezPimenta.pdf: 58058787 bytes, checksum: a0cdeeb6fca3ecffb891becf62a09da5 (MD5) / O presente assunto desta tese se trata de simulações numéricas e novas considerações na modelagem física do ruído gerado por interação do rotor e estator em turbofans aeronáuticos. É feito um estudo teórico preliminar de propagação acústica em dutos com escoamento uniforme, onde é modelada a fonte acústica a partir de preceitos cíclicos da interação rotor e estator. Então são avaliada as características de propagação desses ditos modos de interação, juntamente com a aproximação de dutos anulares finos a ser utilizada nas simulações numéricas, simplificando o grau de liberdade radial do problema de geração e propagação acústica. É feito um estudo teórico da hipótese de escoamento fora do equilíbrio termodinâmico, quando são envolvidas altas frequências características no escoamento e elevadas amplitudes de oscilação. A partir de um modelo pré-existente de viscosidade expansional, é feito um estudo teórico unidimensional no domínio da frequência para baixas amplitudes, em que todos os grupos adimensionais identificados são avaliados na sensibilidade de seus parâmetros. Para o estudo numérico preliminar da viscosidade expansional em condições de propagação unidimensional e no domínio do tempo, é proposto um esquema numérico baseado em diferenças finitas compactas e passo temporal do tipo Runge-Kutta, ambos de alta ordem de precisão numérica. A validação desse esquema numérico é feita com a modificação das equações governantes para se adequar a equação de Burgers viscosa, onde a solução do caso de um salto de descontinuidade é comparada com os valores numéricos obtidos. Modificações também são propostas para o já existente código de simulação numérica VAT (Virtual Aeroacoustic Tunnel), onde um novo esquema de interpolação das faces dos volumes de controle é proposto com argumentos espectrais, visando melhorar as características de propagação do código. Também são propostas modificações no esquema de viscosidade artificial para melhorar sua estabilidade numérica e é feita uma nova implementação do código em CUDA Fortran para utilização de GPUs (Graphical Processing Units) para o cálculo numérico. No estudo unidimensional numérico da viscosidade expansional, é feita a análise de sensibilidade dos parâmetros sobre a dissipação da amplitude da onda estudada. Também é feita uma comparação do decaimento de amplitude de onda para vários valores iniciais de amplitude, em que são comparados os valores numéricos e teóricos e identificada a barreira de linearidade do regime de propagação. São feitas simulações bidimensionais de interação rotor e estator com a fronteira imersa móvel que equivalem a aproximação anular e seus resultados são comparados com a teoria desenvolvida nos regimes de escoamento subsônico, transônico e supersônico. Todos os regimes apresentaram excelente concordância teórica no quesito geração e propagação modal e ainda na propagação de ondas de choque em variadas condições geométricas de geração, incluindo a condição realista de ruído de serra elétrica. Uma modificação do modelo de viscosidade expansional é proposta para a sua aplicação em métodos numéricos que efetuem a marcha temporal. Uma primeira aplicação é feita e seus parâmetros testados para o caso supersônico, de elevada amplitude de onda. Os resultados indicaram que o modelo atuou somente no conteúdo relacionado ao divergente do campo de velocidade do escoamento, onde foi mantida a característica dos outros fatores do escoamento, como a vorticidade e a definição dos corpos no escoamento com a metodologia de fronteira imersa móvel. _______________________________________________________________________________________________ ABSTRACT / The subject of this thesis is about numerical simulation and new physical considerations on the noise generated by rotor and stator interaction in aeronautic turbofans. A preliminary theoretical study on noise propagation in ducts with uniform flow was carried, where the acoustic noise source is characterized after cyclic conditions in time and space of the rotor and stator interaction. Then the propagation characteristics of these modes from the noise sources were assessed with the thin duct approximation, simplifications that are used in the numerical simulations, resulting in a reduction of the radial degree of freedom from the noise generation and propagation. A theoretical study of the sound propagation out of the thermodynamic equilibrium was carried out for cases with high frequency and large amplitude waves. A one dimensional theoretical study in frequency domain for waves with infinitesimal amplitude was carried from an already existent expansional viscosity model, where the nondimensional parameters ware identified and their influence on the flow was assessed. For the numerical study of the expansional viscosity in one dimensional flows, a numerical scheme based on compact finite difference for space derivatives and Runge-Kutta time stepping was proposed, both with high order of numerical precision. The validation of this numerical scheme was done with some simplifications on the governing equations to result in the viscous Burgers’ equation. An exact solution of traveling discontinuity on velocity is compared with the obtained numerical solution. Scheme modifications were proposed for the already existent numerical code VAT (Virtual Aeroacoustic Tunnel), where a new control volume face interpolation based on spectral characteristics to improve the code propagating capabilities was used. Modifications in the artificial viscosity were also proposed to achieve a better numerical stability condition. A code implementation in CUDA Fortran was done to use GPUs (Graphical Processing Unit) in the numerical calculations. On the one dimensional numerical study of the expansional viscosity, a sensitivity analysis of the physical parameters is carried to assess its influence on the wave amplitude decay. Also the decay of several initial wave amplitudes of the numerical results were compared with the theoretical values, and the barrier of the linear propagation regime of initial wave amplitude was identified. Two dimensional cases of rotor and stator interaction were simulated with the moving immersed boundary methodology, where its results were compared with the theoretical ones with the thin duct approximation for the subsonic, transonic and supersonic flow regimes. All the flow regimes were in excellent agreement with the generation and propagation theory of the interaction modes, including shock-wave generation and propagation in several rotor geometric conditions, resulting in the buzz-saw noise. A modification of the expansional viscosity for numerical codes in time domain was proposed. Simulations were carried out and its physical parameters were tested for the supersonic rotor case. Preliminary results indicated that the expansional viscosity model acted only on the velocity divergence related content, where other flow features, such as vorticity and the immersed boundary were kept intact.