Made available in DSpace on 2014-06-11T19:24:47Z (GMT). No. of bitstreams: 0 Previous issue date: 2005-02-21Bitstream added on 2014-06-13T19:31:51Z : No. of bitstreams: 1 pereira_vd_me_ilha.pdf: 1970071 bytes, checksum: 74646e1883439e38fa62ff7f34d06488 (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Neste trabalho foram feitas simulações de escoamentos incompressiveis por um método de elementos finitos de mínimos quadrados (LSFEM – Least Squares Finite Element Method), usando as formulações velocidade-pressão-vorticidade e velocidade-pressão-tensão, denominadas na literatura de formulações u − p −ω e u = p −τ respectivamente. Estas formulações são preferidas por resultarem em sistemas de equações diferenciais de primeira ordem, o que é mais conveniente para implementação pelo LSFEM. O objetivo principal deste trabalho é a simulação computacional de escoamentos laminares, transicionais e turbulentos através da aplicação da metodologia de simulação de grandes escalas (LES – Large Eddy Simulation) com o modelo de viscosidade turbulenta de Smagorinky para modelar as tensões submalha. Alguns problemas padrões foram resolvidos para validar um código computacional desenvolvido e os resultados são apresentados e comparados com resultados disponíveis na literatura. / In this work simulations of incompressible fluid flows have been done by a Least Squares Finite Element Method (LSFEM) using the velocity-pressure-vorticity and velocity-pressurestress formulations, named, in the literature, u − p −ω and u = p −τ formulations respectively. These formulations are preferred because the resulting equations are partial differential equations of first order, which is more convenient for implementation by LSFEM. The main purpose of this work are the numerical computations of laminar, transitional and turbulent fluid flows through the application of large eddy simulation (LES) methodology using the LSFEM. The Navier- Stokes equations in u − p −ω and u = p −τ formulations are filtered and the eddy viscosity model of Smagorinsky is used for modeling the sub-grid-scale stresses. Some benchmark problems are solved for validate a developed numerical code and the preliminary results are presented and compared with available results from the literature.

Flow induced self-oscillations cause acoustic pressure oscillations of large amplitude in pipe flows. If Reynolds number is treated as a parameter, these floinduced oscillations occur only at discrete and critical values of Reynolds number. However, for a small range of Reynolds numbers around such a critical value, such periodic oscillations may appear intermittently. If intermittency, which is a precursor to these self-oscillations, can be detected, prediction of an impending instability may be possible. In experiments done by Vineeth and Sujith (Int. J. Aeroacoustics, 2016) on flow in a duct orifice arrangement, where flow enters through the duct inlet, and leaves into the atmosphere through the orifice exit, “whistling” was observed at a Reynolds number of 4200 (based on the orifice thickness and flow speed within the orifice), where large amplitude pressure oscillations were observed. At slightly lower Reynolds numbers, bursts of relatively smaller amplitudes of pressure oscillations were observed to appear intermittently. For a similar configuration, Large Eddy Simulations (LES) have been carried out with explicit filtering as a sub­ grid scale model here. Both whistling and intermittency are observed in the simulations. As air flows from the duct into the orifice, it turns sharply around the corner at the duct­ orifice interface. Due to this sharp turn, flow separation occurs, and hence, a shear layer is formed at the mouth of the orifice. The mechanism of whistling is found to be this shear layer within the orifice flapping about and hitting the trailing edge of the orifice periodically, thus causing the shear layer to break and roll up into a vortex. At Reynolds numbers where intermittency is observed, the shear layer is found to very mildly come in contact with the edges of the orifice walls, thus disturbing it. In the simulations, time series data of pressure are recorded at various probe locations. In a given time series, if scale invariance behaviour exists, it can be quantified by measuring the Hurst exponent. The numerical value of the Hurst exponent is an index of “long range or short range dependence” in a time series. Hurst exponent is measured in the time series data obtained. It is found to drop to zero as the flow approaches the state of a self-sustained oscillation, since the growth rates of all the long term as well as short term trends in the time series vanish. A loss of multifractality in the time series is also observed as the flow approaches whistling. As a part of the this thesis, new, split high resolution schemes of high order are designed following the Hixon Turmel Proposal.

Orientador: João Batista Campos Silva / Banca: João Batista Aparecido / Banca: Luiz Felipe Mendes de Moura / Resumo: Neste trabalho foram feitas simulações de escoamentos incompressiveis por um método de elementos finitos de mínimos quadrados (LSFEM - Least Squares Finite Element Method), usando as formulações velocidade-pressão-vorticidade e velocidade-pressão-tensão, denominadas na literatura de formulações u − p −ω e u = p −τ respectivamente. Estas formulações são preferidas por resultarem em sistemas de equações diferenciais de primeira ordem, o que é mais conveniente para implementação pelo LSFEM. O objetivo principal deste trabalho é a simulação computacional de escoamentos laminares, transicionais e turbulentos através da aplicação da metodologia de simulação de grandes escalas (LES - Large Eddy Simulation) com o modelo de viscosidade turbulenta de Smagorinky para modelar as tensões submalha. Alguns problemas padrões foram resolvidos para validar um código computacional desenvolvido e os resultados são apresentados e comparados com resultados disponíveis na literatura. / Abstract: In this work simulations of incompressible fluid flows have been done by a Least Squares Finite Element Method (LSFEM) using the velocity-pressure-vorticity and velocity-pressurestress formulations, named, in the literature, u − p −ω and u = p −τ formulations respectively. These formulations are preferred because the resulting equations are partial differential equations of first order, which is more convenient for implementation by LSFEM. The main purpose of this work are the numerical computations of laminar, transitional and turbulent fluid flows through the application of large eddy simulation (LES) methodology using the LSFEM. The Navier- Stokes equations in u − p −ω and u = p −τ formulations are filtered and the eddy viscosity model of Smagorinsky is used for modeling the sub-grid-scale stresses. Some benchmark problems are solved for validate a developed numerical code and the preliminary results are presented and compared with available results from the literature. / Mestre

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / The turbulent flow simulation through the Boussinesq s hypothesis is represented, currently, by two distinct methodologies, the Large-Eddy Simulation (LES) and the Reynolds Averaged Navier-Stokes Equations (RANS). New Hybrid RANS/LES methods are in development, taking off advantage of LES and RANS potentialities through a one only model. The present work deals with the evaluation of these three methodologies, LES, RANS and Hybrid RANS/LES through the turbulent backward-facing step flow simulation. This classical flow is a benchmark for new turbulence models due to the fact that, despite its simple geometry, it presents a very complex generation of three-dimensional structures, influencing the transition phenomenon and properties such as characteristics frequencies of vortex emission and reattachment length. Parallel to this, an inlet turbulent boundary condition influence study showed that the statistical and topological content of the inlet boundary layer profile can modify substantially results like reattachment length and pressure coefficient. A recycling method for generating three-dimensional, time-dependent turbulent boundary layer inflow data for Large-Eddy and Hybrid RANS/LES simulation is employed. Results for the three methodologies disclose that Large-Eddy Simulation and Hybrid RANS/LES methods present very similar descriptions for the turbulent backward-facing step flow, differing from the RANS s results, where the second order statistical moments are totally suppressed, with absence of three-dimensional and transient structures. / A simulação numérica de escoamentos turbulentos através da hipótese de Boussinesq é representada, atualmente, por duas grandes metodologias distintas, a Simulação de Grandes Escalas (LES Large-Eddy Simulation) e as Equações Médias de Reynolds (RANS Reynolds Averaged Navier-Stokes). Uma nova metodologia, chamada de Híbrida RANS/LES, está em desenvolvimento, tirando proveito das potencialidades das metodologias tradicionais LES e RANS através de um único modelo. O presente trabalho trata da avaliação das três metodologias, LES, RANS e Híbrida RANS/LES de modelagem da turbulência através da simulação numérica do escoamento turbulento sobre um degrau. Os modelos são avaliados através deste escoamento, que apesar de simples geometricamente, é capaz de gerar um escoamento complexo, com regiões de escoamento parietal e cisalhante livre. Juntamente com a modelagem da turbulência, um estudo de imposição de condições de contorno turbulentas na entrada do domínio utilizado revelou que tão importante quanto o modelo de turbulência, as condições de contorno empregadas modificam substancialmente os resultados obtidos. Foi implementado um modelo de geração de contorno baseado no escalonamento de informações internas do escoamento de forma a satisfazer estatística e topologicamente o caráter turbulento da condição de contorno na entrada. Resultados para as três metodologias revelam que a Simulação de Grandes Escalas e métodos Híbridos RANS/LES apresentam descrições muito semelhantes para o escoamento turbulento sobre o degrau, diferindo dos resultados da metodologia RANS, onde momentos estatísticos de segunda ordem são suprimidos, com ausência de estruturas tridimensionais e transientes. / Mestre em Engenharia Mecânica

O experimento Unsteady Aerodynamics Experiment Phase VI, realizado no ano de 2000 pelo Laboratório Nacional norte-americano para as Energias Renováveis (NREL) no túnel de vento Ames da NASA, foi reproduzido numericamente neste trabalho. O objetivo é o estudo das características da esteira aerodinâmica produzida pela turbina eólica de duas pás e 10 metros de diâmetro, operando à velocidade de rotação constante de 72 RPM, sujeita a uma velocidade de corrente livre do vento uniforme de 9 m/s, em um túnel de vento cuja seção de testes mede 36,6 m de largura por 24,4 m de altura e o comprimento mede 170 m. Para isso, foi utilizado o programa comercial ANSYS FLUENT versão 13.0, baseado no Método dos Volumes Finitos para a solução numérica das Equações de Navier-Stokes em regime transiente em conjunto com a Simulação de Grandes Escalas (SGE) para resolver a turbulência. As geometrias de todos os componentes da máquina foram criadas em software CAD. Um domínio móvel em forma de disco, contendo as pás do rotor e o hub da máquina, foi criado separadamente, e posteriormente inserido no domínio principal, estático, usando a ferramenta Moving Mesh disponível no software FLUENT. Ambos os domínios foram preenchidos por malhas compostas por tetraedros. Dados provenientes das simulações numéricas foram comparados aos dados experimentais de velocidade fornecidos por dois anemômetros sônicos instalados 5,8 m à jusante do rotor, ao que foi verificada boa concordância, com diferenças da ordem de 1% para o anemômetro 1 e 6% para o anemômetro 2. Resultados de velocidade na linha de centro do túnel e perfis de velocidade à jusante foram comparados com recente estudo numérico, e revelam diferenças importantes entre dados obtidos pela SGE, principalmente no que se refere à detecção de picos e flutuações relacionados às escalas turbulentas, e dados obtidos através da modelagem clássica da turbulência, RANS. As perturbações ultrapassaram a marca dos 10 diâmetros à jusante e atingiram o final do domínio localizado a 15 diâmetros. A esteira não apresentou simetria axial, e o ponto de maior redução na velocidade do escoamento foi detectado fora da linha de centro do rotor. / The Unsteady Aerodynamics Experiment Phase VI, which has been carried out in 2000 by the US National Renewable Energy Laboratory (NREL) at the NASA Ames wind tunnel, has been numerically reproduced. The purpose of this work is to study the characteristics of the wind wake produced by the 10 meter two bladed wind turbine, operating at a constant rotational speed of 72 RPM, subject to a free stream wind velocity of 9 m/s, inside a wind tunnel in which dimensions are 36.6 m in width, 24.4 m in height and length of 170 m. To achieve that, the ANSYS FLUENT version 13.0 commercial code, based in the Finite Volume Method to numerically solve the Navier-Stokes equations in transient state, has been used, together with the Large Eddy Simulation (LES) to characterize the turbulence. Geometries of all the machine components have been created in CAD software. A disc shaped moving domain, containing the blades and hub, has been created separately, and later inserted into the main, static domain, using the Moving Mesh tool available in the software. Both domains have been filled with meshes composed by tetrahedra. Data collected at the numerical simulations have been compared to experimental wind speed data provided by two sonic anemometers installed 5.8 m downstream from the rotor, for which a good agreement has been found, with differences of approximately 1% to the anemometer 1 and 6% to the anemometer 2. Results of wind velocity at the tunnel centerline and velocity profiles downstream have been compared with recent numerical study, and show important differences between data obtained by LES, especially with regard to the detection of peaks and fluctuations related to the turbulent scales, and data obtained by the classic turbulence modeling, RANS. Disturbances have passed the 10 diameter mark and reached the end at the domain located at 15 diameters. The wake did not show axial symmetry and the point of maximum reduction in the flow speed was detected outside the rotor centerline.

O presente trabalho tem como objetivo estudar escoamentos de tornados e sua ação sobre corpos imersos empregando ferramentas numéricas da Engenharia do Vento Computacional (EVC). Os tornados constituem-se atualmente em uma das causas de desastres naturais no Brasil, especialmente nas regiões sul e sudeste do país, como também em alguns países vizinhos. Os efeitos gerados são geralmente localizados e de curta duração, podendo ser devastadores dependendo da escala do tornado. Tais características dificultam a realização de estudos detalhados a partir de eventos reais, o que levou ao desenvolvimento de modelos experimentais e numéricos. A abordagem numérica é utilizada neste trabalho para a simulação de tornados, a qual se baseia nas equações de Navier-Stokes e na equação de conservação de massa, considerando a hipótese de pseudo-compressibilidade e condições isotérmicas. Para escoamentos com turbulência utiliza-se a Simulação Direta de Grandes Escalas com o modelo clássico de Smagorinsky para as escalas inferiores à resolução da malha (Large Eddy Simulation ou LES em inglês). A discretização das equações fundamentais do escoamento se realiza com um esquema explícito de dois passos de Taylor-Galerkin, onde o Método dos Elementos Finitos é empregado na discretização espacial utilizando-se o elemento hexaédrico trilinear isoparamétrico com um ponto de integração e controle de modos espúrios Na presença de corpos imersos que se movem para simular os deslocamentos dos tornados, o escoamento é descrito cinematicamente através de uma formulação Arbitrária Lagrangeana-Euleriana (ALE) que inclui um esquema de movimento de malha. Tornados são reproduzidos através da simulação numérica de dispositivos experimentais e do Modelo de Vórtice Combinado de Rankine (RCVM). Exemplos clássicos da Dinâmica dos Fluidos Computacional são apresentados inicialmente para a verificação das ferramentas numéricas implementadas. Finalmente, problemas envolvendo tornados móveis e estacionários são analisados, incluindo sua ação sobre corpos imersos. Nos modelos baseados em experimentos, a variação da relação de redemoinho determinou os diferentes padrões de escoamento observados no laboratório. Nos exemplos de modelo de vórtice, quando o tornado impactou o corpo imerso gerou picos de forças em todas as direções e, após a passar pelo mesmo, produziu uma alteração significativa na estrutura do vórtice. / Analyses of tornado flows and its action on immersed bodies using numerical tools of Computational Wind Engineering (CWE) are the main aims of the present work. Tornadoes are currently one of the causes of natural disasters in Brazil, occurring more frequently in the southern and southeastern regions of the country, as well as in some neighboring countries. Effects are usually localized, presenting a short time interval, which can be devastating depending on the scale of the tornado. These characteristics difficult to carry out detailed studies based on real events, leading to the development of experimental and numerical models. The numerical approach is used in this work for the simulation of tornadoes, which is based on the Navier-Stokes equations and the mass conservation equation, considering the hypothesis of pseudo-compressibility and isothermal conditions. For turbulent flows, Large Eddy Simulation (LES) is used with the classical Smagorinsky model for sub-grid scales Discretization is performed the explicit two-step Taylor-Galerkin scheme, where the Finite Element Method is used in spatial discretization using isoparametric trilinear hexahedral elements with one-point quadrature and hourglass control. In the presence of immersed bodies that are moving in order to simulate translating tornadoes, the flow is kinematically described through a Lagrangian-Eulerian Arbitrary (ALE) formulation, which includes a mesh motion scheme. Tornadoes are reproduced using numerical simulation of experimental devices and the Rankine Combined Vortex Model (RCVM). Classical examples of Computational Fluid Dynamics are presented initially for the verification of the numerical tools implemented here. Finally, problems involving moving and stationary tornadoes are analyzed, including their actions on immersed bodies. For models based on experiments, the variation of the swirl ratio determined the different flow patterns observed in the laboratory. In the vortex model examples, when the tornado impacted on the immersed body, peaks of forces were generated in all directions and, after passing over it, a significant change in the structure of the vortex was produced.

A interação turbulência-radiação (TRI, do inglês Turbulence-Radiation Interaction) resulta do acoplamento altamente não linear entre flutuações da intensidade de radiação e flutuações da temperatura e da composição química do meio, e tem-se demonstrado experimentalmente, teoricamente e numericamente que este é um fenômeno relevante em diversas aplicações envolvendo altas temperaturas, especialmente em problemas reativos. Neste trabalho, o TRI é analisado em um escoamento interno não reativo de um gás participante que se desenvolve em um duto de seção transversal quadrada, para diferentes intensidades de turbulência do escoamento e considerando duas espécies distintas para a composição do fluido de trabalho (dióxido de carbono e vapor de água). O objetivo central é avaliar como a inclusão ou não da variação espectral das propriedades radiativas do meio no cálculo influencia a magnitude do TRI. Isso é feito através de simulações numéricas no código de dinâmica dos fluidos computacional Fire Dynamics Simulator (FDS), que resolve, através do método dos volumes finitos, as equações fundamentais que regem o problema – isto é, os balanços de massa, de quantidade de movimento e de energia e a equação de estado – em uma formulação adequada para baixos números de Mach, utilizando um algoritmo de solução explícito e de segunda ordem no tempo e no espaço. A turbulência é modelada através da simulação de grandes escalas (LES, do inglês Large Eddy Simulation), empregando-se o modelo de Smagorinsky dinâmico para o fechamento dos termos submalha; para a radiação térmica, o método dos volumes finitos é utilizado na discretização da equação da transferência radiativa e os modelos do gás cinza e da soma-ponderada-de-gases-cinza (WSGG, do inglês Weighted-Sum-of-Gray-Gases) são implementados como forma de desconsiderar e de incluir a dependência espectral das propriedades radiativas, respectivamente. A magnitude do TRI sobre o problema é avaliada através de diferenças entre as médias temporais dos fluxos de calor superficiais e do termo fonte radiativo obtidas em cálculos que consideram os efeitos do fenômeno e cálculos que os negligenciam. Em geral, a interação turbulência-radiação mostrou ser pouco importante em todos os casos considerados, o que concorda com resultados de outros estudos sobre o tema em escoamento não reativos. Com o modelo WSGG, as contribuições do fenômeno foram maiores do que com a hipótese do gás cinza, evidenciando que a inclusão da variação espectral na solução do problema radiativo tem um impacto sobre a magnitude dos efeitos do TRI. Além disso, é feita uma discussão, em parte inédita no contexto do TRI, sobre diferentes metodologias para a análise do fenômeno. Finalmente, é proposto um fator de correção para o termo fonte radiativo médio no modelo WSGG, que é validado através de sua implementação nos casos simulados. Em estudos futuros, uma análise de sensibilidade sobre os termos constituintes desse fator de correção pode levar a um melhor entendimento de como as flutuações de temperatura se correlacionam com o fenômeno da interação turbulência-radiação. / Turbulence-radiation interaction (TRI) results from the highly non-linear coupling between fluctuations of radiation intensity and fluctuations of temperature and chemical composition of the medium, and its relevance in a number of high-temperature problems, especially when chemical reactions are included, has been demonstrated experimentally, theoretically, and numerically. In the present study, the TRI is analyzed in a channel flow of a non-reactive participating gas for different turbulence intensities of the flow at the inlet and considering two distinct species for the medium composition (carbon dioxide and water vapor). The central objective is to evaluate how the inclusion or not of the spectral variation of the radiative properties of a participating gas in the radiative transfer calculations affects the turbulence-radiation interaction. With this purpose, numerical simulations are performed using the computational fluid dynamics Fortranbased code Fire Dynamics Simulator, that employs the finite volume method to solve a form of the fundamental equations – i.e., the mass, momentum and energy balances and the state equation – appropriate for low Mach number flows, through an explicit second-order (both in time and in space) core algorithm. Turbulence is modeled by the large eddy simulation approach (LES), using the dynamic Smagorinsky model to close the subgrid-scale terms; for the thermal radiation part of the problem, the finite volume method is used for the discretization of the radiative transfer equation and the gray gas and weighted-sum-of-gray-gases (WSGG) models are implemented as a way to omit and consider the spectral dependence of the radiative properties, respectively. The TRI magnitude in the problem is evaluated by differences between values for the time-averaged heat fluxes at the wall (convective and radiative) and for the time-averaged radiative heat source calculated accounting for and neglecting the turbulence-radiation interaction effects. In general, TRI had little importance over all the considered cases, a conclusion that agrees with results of previous studies. When using the WSGG model, the contributions of the phenomenon were greater that with the gray gas hypothesis, demonstrating that the inclusion of the spectral variance in the solution of the radiative problem has an impact in the TRI effects. Furthermore, this paper presents a discussion, partly unprecedented in the context of the turbulence-radiation interaction, about the different methodologies that can be used for the TRI analysis. Finally, a correction factor is proposed for the time-averaged radiative heat source in the WSGG model, which is then validated by its implementation in the simulated cases. In future studies, a sensibility analysis on the terms that compose this factor can lead to a better understanding of how fluctuations of temperature correlate with the turbulence-radiation interaction phenomenon.

O objetivo deste trabalho é estudar, dentro da Engenharia do Vento Computacional (EVC), a dispersão de poluentes em zonas urbanas, empregando-se um modelo numérico baseado em técnicas da Dinâmica dos Fluidos Computacional para escoamentos incompressíveis, não isotérmicos e com transporte de massa. Um esquema explícito de dois passos é usado para a discretização temporal das equações governantes, considerando expansões em séries de Taylor de segunda ordem para as derivadas no tempo. O processo de discretização espacial é realizado através da aplicação do Método dos Elementos Finitos (MEF), onde hexaedros de oito nós com um ponto de integração são utilizados. A turbulência é tratada numericamente através da Simulação de Grandes Escalas (LES) e os modelos clássico e dinâmico de Smagorinsky são empregados na modelagem das escalas inferiores à resolução da malha. Efeitos de temperatura sobre o escoamento são considerados na forma de forças de flutuação presentes na equação de balanço de momentum, as quais são calculadas a partir da aproximação de Boussinesq. Técnicas de paralelização em memória compartilhada (OpenMP) são também usadas a fim de melhorar a eficiência computacional do presente modelo para problemas com grande número de elementos. Exemplos clássicos de Dinâmica de Fluidos e Fenômenos de Transporte são inicialmente analisados para teste das ferramentas numéricas implementadas. Problemas de dispersão de poluentes com e sem a inclusão dos efeitos de temperatura são abordados para configurações geométricas bi e tridimensionais de street canyons, representando a unidade geométrica básica encontrada em centros urbanos de grandes cidades. / The main goal of the present work is to study the pollutant dispersion in urban areas using a numerical model based on techniques developed by Computational Fluid Dynamics, where applications of Computational Wind Engineering (CWE) are analyzed considering incompressible flows with heat and mass transport. A two-step explicit scheme is adopted for the time discretization of the governing equations considering second order Taylor series expansions of the time derivative terms. Spatial discretization is performed by applying the Finite Element Method (FEM), where eight-node hexahedral elements with one-point quadrature are utilized. Turbulence is numerically analyzed by using Large Eddy Simulation (LES) with the classical and dynamic Smagorinsky’s models for subgrid scale modeling. Thermal effects on the flow field are taken into account through buoyancy forces acting on the momentum balance equation, which are calculated considering the Boussinesq approximation. Shared memory parallelization techniques (OpenMP) are also employed in order to improve computational efficiency for problems with large number of elements. Classic examples of Fluid Dynamics and Transport Phenomena are first analyzed to verify the numerical tools implemented. Problems involving pollutant dispersion with and without the inclusion of thermal effects are investigated for two and three-dimensional geometric configurations of street canyons, which represent the basic geometric unit observed in urban centers of large cities.

O objetivo deste trabalho é estudar, dentro da Engenharia do Vento Computacional (EVC), a dispersão de poluentes em zonas urbanas, empregando-se um modelo numérico baseado em técnicas da Dinâmica dos Fluidos Computacional para escoamentos incompressíveis, não isotérmicos e com transporte de massa. Um esquema explícito de dois passos é usado para a discretização temporal das equações governantes, considerando expansões em séries de Taylor de segunda ordem para as derivadas no tempo. O processo de discretização espacial é realizado através da aplicação do Método dos Elementos Finitos (MEF), onde hexaedros de oito nós com um ponto de integração são utilizados. A turbulência é tratada numericamente através da Simulação de Grandes Escalas (LES) e os modelos clássico e dinâmico de Smagorinsky são empregados na modelagem das escalas inferiores à resolução da malha. Efeitos de temperatura sobre o escoamento são considerados na forma de forças de flutuação presentes na equação de balanço de momentum, as quais são calculadas a partir da aproximação de Boussinesq. Técnicas de paralelização em memória compartilhada (OpenMP) são também usadas a fim de melhorar a eficiência computacional do presente modelo para problemas com grande número de elementos. Exemplos clássicos de Dinâmica de Fluidos e Fenômenos de Transporte são inicialmente analisados para teste das ferramentas numéricas implementadas. Problemas de dispersão de poluentes com e sem a inclusão dos efeitos de temperatura são abordados para configurações geométricas bi e tridimensionais de street canyons, representando a unidade geométrica básica encontrada em centros urbanos de grandes cidades. / The main goal of the present work is to study the pollutant dispersion in urban areas using a numerical model based on techniques developed by Computational Fluid Dynamics, where applications of Computational Wind Engineering (CWE) are analyzed considering incompressible flows with heat and mass transport. A two-step explicit scheme is adopted for the time discretization of the governing equations considering second order Taylor series expansions of the time derivative terms. Spatial discretization is performed by applying the Finite Element Method (FEM), where eight-node hexahedral elements with one-point quadrature are utilized. Turbulence is numerically analyzed by using Large Eddy Simulation (LES) with the classical and dynamic Smagorinsky’s models for subgrid scale modeling. Thermal effects on the flow field are taken into account through buoyancy forces acting on the momentum balance equation, which are calculated considering the Boussinesq approximation. Shared memory parallelization techniques (OpenMP) are also employed in order to improve computational efficiency for problems with large number of elements. Classic examples of Fluid Dynamics and Transport Phenomena are first analyzed to verify the numerical tools implemented. Problems involving pollutant dispersion with and without the inclusion of thermal effects are investigated for two and three-dimensional geometric configurations of street canyons, which represent the basic geometric unit observed in urban centers of large cities.

O experimento Unsteady Aerodynamics Experiment Phase VI, realizado no ano de 2000 pelo Laboratório Nacional norte-americano para as Energias Renováveis (NREL) no túnel de vento Ames da NASA, foi reproduzido numericamente neste trabalho. O objetivo é o estudo das características da esteira aerodinâmica produzida pela turbina eólica de duas pás e 10 metros de diâmetro, operando à velocidade de rotação constante de 72 RPM, sujeita a uma velocidade de corrente livre do vento uniforme de 9 m/s, em um túnel de vento cuja seção de testes mede 36,6 m de largura por 24,4 m de altura e o comprimento mede 170 m. Para isso, foi utilizado o programa comercial ANSYS FLUENT versão 13.0, baseado no Método dos Volumes Finitos para a solução numérica das Equações de Navier-Stokes em regime transiente em conjunto com a Simulação de Grandes Escalas (SGE) para resolver a turbulência. As geometrias de todos os componentes da máquina foram criadas em software CAD. Um domínio móvel em forma de disco, contendo as pás do rotor e o hub da máquina, foi criado separadamente, e posteriormente inserido no domínio principal, estático, usando a ferramenta Moving Mesh disponível no software FLUENT. Ambos os domínios foram preenchidos por malhas compostas por tetraedros. Dados provenientes das simulações numéricas foram comparados aos dados experimentais de velocidade fornecidos por dois anemômetros sônicos instalados 5,8 m à jusante do rotor, ao que foi verificada boa concordância, com diferenças da ordem de 1% para o anemômetro 1 e 6% para o anemômetro 2. Resultados de velocidade na linha de centro do túnel e perfis de velocidade à jusante foram comparados com recente estudo numérico, e revelam diferenças importantes entre dados obtidos pela SGE, principalmente no que se refere à detecção de picos e flutuações relacionados às escalas turbulentas, e dados obtidos através da modelagem clássica da turbulência, RANS. As perturbações ultrapassaram a marca dos 10 diâmetros à jusante e atingiram o final do domínio localizado a 15 diâmetros. A esteira não apresentou simetria axial, e o ponto de maior redução na velocidade do escoamento foi detectado fora da linha de centro do rotor. / The Unsteady Aerodynamics Experiment Phase VI, which has been carried out in 2000 by the US National Renewable Energy Laboratory (NREL) at the NASA Ames wind tunnel, has been numerically reproduced. The purpose of this work is to study the characteristics of the wind wake produced by the 10 meter two bladed wind turbine, operating at a constant rotational speed of 72 RPM, subject to a free stream wind velocity of 9 m/s, inside a wind tunnel in which dimensions are 36.6 m in width, 24.4 m in height and length of 170 m. To achieve that, the ANSYS FLUENT version 13.0 commercial code, based in the Finite Volume Method to numerically solve the Navier-Stokes equations in transient state, has been used, together with the Large Eddy Simulation (LES) to characterize the turbulence. Geometries of all the machine components have been created in CAD software. A disc shaped moving domain, containing the blades and hub, has been created separately, and later inserted into the main, static domain, using the Moving Mesh tool available in the software. Both domains have been filled with meshes composed by tetrahedra. Data collected at the numerical simulations have been compared to experimental wind speed data provided by two sonic anemometers installed 5.8 m downstream from the rotor, for which a good agreement has been found, with differences of approximately 1% to the anemometer 1 and 6% to the anemometer 2. Results of wind velocity at the tunnel centerline and velocity profiles downstream have been compared with recent numerical study, and show important differences between data obtained by LES, especially with regard to the detection of peaks and fluctuations related to the turbulent scales, and data obtained by the classic turbulence modeling, RANS. Disturbances have passed the 10 diameter mark and reached the end at the domain located at 15 diameters. The wake did not show axial symmetry and the point of maximum reduction in the flow speed was detected outside the rotor centerline.