O presente estudo apresenta a análise numérica do escoamento do fluido de trabalho em uma turbina Tesla com rotor estacionário. O estudo de independência de malha prevê o uso de aproximadamente 2 milhões de volumes, com desempenho semelhante ao das malhas mais refinadas, apresentando uma economia significativa de esforço computacional. Três modelos de turbulência da abordagem RANS são aplicados com o objetivo de estabelecer uma metodologia para estudos futuros com dados experimentais. Os diferentes modelos de turbulência fornecem resultados para a predição do torque na turbina, com uma variação abaixo de 1 % entre si. Ar é usado como fluido de trabalho a pressão manométrica de 2,5 bar, alcançando velocidades no entorno de 310 m.s-1 na região da garganta do bocal e na região de jato livre, logo após a descarga do bocal. Essa condição permite afirmar que a turbina funciona na sua condição de máxima vazão, com o número de Mach próximo ao valor unitário, com escoamento compressível A velocidade na região interna entre discos chega a um valor máximo de 100 m.s-1. Na ausência de dados experimentais e de literatura, um estudo paramétrico com diferentes condições de operação da turbina é realizado a fim de verificar a qualidade dos resultados simulados. A vazão mássica é estimada com base na temperatura e pressão do fluido de trabalho, modelado como gás ideal. Os resultados preditos pelo modelo numérico para o torque no rotor é de 2,09 N.m com pressão manométrica de 1,5 bar e vazão mássica de 33,58 g/s, 2,22 N.m com pressão manométrica de 2,0 bar e vazão mássica de 40,29 g/s, e 2,38 N.m com pressão manométrica de 2,5 bar e vazão mássica de 53,73 g/s. A temperatura foi de 300 K mantida constante para as três análises. Para os casos analisados, o número de Mach na garganta do bocal convergente apresentou uma tendência ao valor unitário, variando entre 0,7 a 1, o que sugere que o bocal está trabalhando na sua condição máxima de vazão do fluido de trabalho. / The present study presents the numerical analysis of the working fluid flow in a Tesla turbine with stationary rotor. The mesh independence study predicts the use of approximately 2 million volumes, with similar performance to those most refined meshes, presenting a significant saving of computational effort. Several turbulence models of the RANS approach are applied with the aim of establishing a methodology for future studies with experimental data. The different turbulence models provide very close results for turbine torque prediction, with a variation below 1% between them. Air is used as working fluid at a pressure of 2.5 bar gauge, reaching velocities around 310 m.s-1 in the throat region of the nozzle and in the free jet region, just after the discharge of the nozzle. This condition allows to state that the turbine works in its maximum flow condition, with the Mach number close to unitary value, with a compressible flow The velocity in the inner region between disks reaches 100 m.s-1. In the absence of experimental data and literature, a parametric study with different operating conditions of the turbine is performed in order to verify the quality of the simulated results. The mass flow rate is estimated based on the temperature and pressure of the working fluid, modeled as the ideal gas. The results predicted by the numerical model for the torque in the rotor is 2.09 N.m with gauge pressure of 1.5 bar and a mass flow rate of 33.58 g/s, 2.22 N.m with gauge pressure of 2.0 bar and a mass flow rate of 40.29 g/s, and 2.38 N.m with gauge pressure of 2.5 bar and flow mass of 53.73 g/s. The temperature was 300 K kept constant for all three cases. For the analyzed cases, the Mach number in the throat of the convergent nozzle showed a tendency to the unit value, ranging from 0,7 to 1,0 which suggests that the nozzle is working in its maximum flow condition of the working fluid.

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.

O formato do perfil de velocidades do vento varia de acordo com as características locais da superfície terrestre e de rugosidade do terreno, parâmetros que definem o perfil da Camada-Limite Atmosférica (CLA). As características do escoamento do ar atmosférico sobre e ao redor de acidentes geográficos, tais como morros e colinas, são de grande interesse para aplicações relacionadas à Engenharia de Turbinas e Parques Eólicos. No topo de morros, ocorre a aceleração do vento, fenômeno que pode representar um fator decisivo para a instalação de aerogeradores. Este trabalho dedica-se ao estudo do comportamento da CLA como função da inclinação e rugosidade superficial da elevação, fazendo uso da Dinâmica de Fluidos Computacional (CFD) para construir perfis de velocidade do vento e de intensidade de turbulência. O problema de fechamento das Equações Médias de Reynolds (RANS) é contornado com o uso do modelo de turbulência k-ω SST; os resultados numéricos obtidos são comparados com dados experimentais medidos em túnel de vento sobre modelos em escala dos morros. São testados oito modelos de morros com declives que variam de 25° a 64° para dois tipos de categorias de terreno, em 2D e 3D, e são aplicados dois códigos analíticos para representar o perfil de velocidades de entrada. Resultados numéricos para os perfis de velocidade apresentam diferença inferior a 4% em relação aos respectivos dados obtidos experimentalmente. Os perfis de intensidade de turbulência apresentam diferença máxima na casa dos 7% em comparação aos dados experimentais, o que é explicado pelo fato de que não é possível inserir o perfil de entrada de intensidade de turbulência nas simulações numéricas. Em alternativa, foi usado um valor constante resultado da média dos valores dos perfis usados no túnel de vento. Os modelos de morro em 3D apresentam maior concordância nos resultados de velocidade que os modelos em 2D e que ademais quanto maior é a inclinação do morro maior é a concordância com as medições experimentais. / The shape of the wind velocity profile changes according to local features of terrain shape and roughness, which are parameters responsible for defining the Atmospheric Boundary Layer (ABL) profile. Air flow characteristics over and around landforms, such as hills, are of considerable importance for applications related to Wind Farm and Turbine Engineering. The air flow is accelerated on top of hills, which can represent a decisive factor for Wind Turbine placement choices. The present work focuses on the study of ABL behavior as a function of slope and surface roughness of hill-shaped landforms, using the Computational Fluid Dynamics (CFD) to build wind velocity and turbulent intensity profiles. Reynolds-Averaged Navier-Stokes (RANS) equations are closed using the SST k-ω turbulence model; numerical results are compared to experimental data measured in wind tunnel over scale models of the hills under consideration. Eight hill models with slopes varying from 25° to 64° were tested for two types of terrain categories in 2D and 3D, and two analytical codes are used to represent the inlet velocity profiles. Numerical results for the velocity profiles show differences under 4% when compared to their respective experimental data. Turbulent intensity profiles show maximum differences around 7% when compared to experimental data, this can be explained by not being possible to insert inlet turbulent intensity profiles in the simulations. Alternatively, constant values based on the averages of the turbulent intensity at the wind tunnel inlet were used. The 3D models present greater concordance in the speed results than the 2D models and that in addition the greater the slope of the hill, the greater the agreement with the experimental measurements.

Neste trabalho, é apresentada a simulação de escoamentos compressíveis turbulentos no entorno de corpos móveis rígidos ou deformáveis empregando-se técnicas adaptativas. As simulações numéricas são conduzidas utilizando-se o método dos elementos finitos. A discretização espaço-temporal é desenvolvida através do método das linhas ou direções características (Characteristic-Based Split - CBS) e a modelagem da turbulência é feita através de um modelo de Simulação de Grandes Escalas (SGE, ou na terminologia em inglês, Large Eddy Simulation – LES) com o coeficiente de Smagorinsky variável no tempo e espaço (SGE ou LES dinâmico). A análise estrutural de corpos deformáveis imersos no fluido é realizada através de um modelo de elementos finitos triangulares para análise de placas e cascas com não linearidade geométrica, usando materiais elásticos com comportamento linear. Conjuntamente, um método de adaptação anisotrópica transiente de malhas é empregado para obter resultados com boa resolução a baixos custos computacionais. A consideração do movimento relativo de corpos imersos no escoamento é feita através de um método híbrido de movimento da malha que emprega interpolação com funções de base radial. Exemplos bidimensionais e tridimensionais são apresentados de forma a validar cada uma das metodologias desenvolvidas. Por fim, exemplos de simulações complexas são investigados, comparando-se os resultados obtidos com resultados experimentais e numéricos presentes na literatura. / In this work, the simulation of compressible turbulent flows around rigid and flexible moving bodies is presented using adaptative techniques. The numerical simulations are solved employing the finite element method. The space-time discretization is performed using the Characteristic-Based Split scheme (CBS) and turbulence is modelled with Large Eddy Simulation (LES) and a dynamic Smagorinsky sub-grid model. The structural analysis of deformable bodies immersed on the flow is performed using a triangular finite element model for the analysis of geometrically non-linear elastic plates and shells. An anisotropic mesh adaptation algorithm for transient simulations is coupled with the solver to achieve results with good resolution and low computational costs. The consideration of the relative movement of immersed bodies on the flow is performed employing an hybrid method of mesh movement based on radial basis function interpolation. Twodimensional and three-dimensional examples are presented in order to validate the proposed methodologies. Finally, complex simulations are investigated, where results are compared with experimental and numerical data available in the literature.

É notório que a degradação ambiental vem ao longo do tempo posicionando-se como um dos principais problemas do mundo moderno. Dentre as várias questões do interesse ambiental podemos destacar a ascensão da bolha de metano, em reservatórios hidrelétricos, desde o sedimento anóxico no fundo do reservatório até a interface água atmosfera. Neste contexto, a presente tese vêm propor uma nova modelagem matemática para a ascensão da bolha axissimétrica em fluidos newtonianos/nãonewtonianos e mostrar resultados numéricos simulados. Desta forma, o estado da arte estaria elevado a posição de permitir, via Matemática e Simulação Numérica-Computacional, a análise do transporte de metano em reservatórios de hidrelétricas através da bolha / It is well-known that environmental degradation has come along the time positioning as one of the main problems from modern world. Among several questions of the environmental interest may emphasize the methane bubble rise in hydroelectric reservoirs from the anoxic sediment in the bottom of reservoir until water interface atmosphere. In this context, the current thesis has come to propose a new mathematic modeling to the rise of the axisymmetric bubble in the newtonian/non-newtonian fluids and display numerical results simulated. Therefore, the state of art would be ascended to a position to permit, via Mathematics and Computing-Numeric Simulation, the analysis of transport of methane in hydroelectric reservoirs through of bubble

As bombas de fluxo, além das aplicações clássicas em Engenharia, são instrumentos importantes em áreas como a biologia, farmácia e medicina. Um novo princípio para o bombeamento de fluidos está sendo estudado dentro do Departamento deEngenharia Mecatrônica e de Sistemas Mecânicos da Escola Politécnica da Universidade de São Paulo. O presente projeto trata deste princípio: a utilização de um atuador piezelétrico bilaminar inserido num meio fluido para geração de vazão. O objetivo deste projeto é verificar a viabilidade desta bomba piezelétrica através de simulações computacionais, estudando a sensibilidade do sistema a diversos parâmetros e realizando a otimização do mesmo de forma a maximizar seudesempenho. O software ANSYS é utilizado para a simulação computacional do problema de dinâmica de fluidos e para modelar o atuador piezelétrico bilaminar, e o software Altair Hyperstudy na etapa de otimização. O texto apresenta as metodologias empregadas e discute os resultados obtidos, de forma a analisar os fenômenos físicos em questão e validar este novo princípio proposto. / Flow pumps, besides their classical applications in Engineering, are important devices in areas such as Biology, Pharmacy and Medicine. A novel principle for fluid pumping has been studied at the Department of Mechatronic Engineering and Mechanical Systems of the Engineering School of the University of Sao Paulo. The current project deals with this principle: the use of a bimorph piezoelectric actuator in a fluid environment for flow generation. The objective of this project is to verify the viability of this pump through computational simulations, by performing a sensitivity analysis for various parameters and an optimization to maximize its performance. The ANSYS software is used for the computational fluid dynamics simulations and for modeling the bimorph piezoelectric actuator, and the Altair Hyperstudy software for the optimization. The text presents the employed methodologies and discusses the obtained data in order to analyze the physical phenomena involved and to validate this novel principle.

Made available in DSpace on 2016-12-23T14:08:08Z (GMT). No. of bitstreams: 1 Martins RS Part 1.pdf: 212782 bytes, checksum: 40ca14a4b97179f5b3af71f964306b97 (MD5) Previous issue date: 2012-05-23 / Oil and gas industry requires accurate flow measurements since they are stated by law. Nevertheless, curves and other obstacles are commonly found in such industry field, which may affect the quality of flow measurement due to flow disturbances, such as swirl and velocity profile asymmetries. Single-path ultrasonic flow meters are often used in flare gas installations, despite being sensitive to such disturbances. The present work use commercial CFD codes to obtain disturbed flow fields downstream from single and double elbow pipe installations, aiming to investigate both magnitude and behaviour of such effects on ultrasonic flow measurement. Numerical integration is applied for several acoustic path arrangements, simulating single-path ultrasonic flow meters in different situations in order to evaluate its correction factor deviation under disturbed conditions. Reynolds numbers from 1X10 4 to 2 X10 6 are considered. Transducers mounting angles from 0° to 180° are tested and axial positions up to 80D downstream from the curve are evaluated. Results indicate that single-path ultrasonic flow meters are sensitive to installation effects. Correction factor deviations usually showed to be significantly higher than 2% for axial positions shorter than 20D, as recommended by several manufacturers or regulations. Nevertheless, deviations may reach 0.01% in some specific configurations, which suggests that ultrasonic flow measurement might be improved by rearranging flow meter device in favourable angular position and mainly by implementation of specific functions for correction factors under disturbed conditions / A indústria de petróleo e gás requer medições de vazão de baixa incerteza, uma vez que são estabelecidas por lei. Contudo, curvas e outros obstáculos são comumente encontrados nesse cenário e podem afetar a qualidade da medição de vazão em função de perturbações no escoamento, tais como swirl e assimetrias no perfil de velocidades. Medidores de vazão por ultrassom de um canal são frequentemente utilizados em instalações de gás de queimadores, apesar de serem sensíveis a tais perturbações. O presente trabalho usa códigos comerciais de CFD para obter o escoamento à jusante de instalações com uma curva e duas curvas, visando a investigar a magnitude e o comportamento de tais efeitos na medição de vazão. Integração numérica é utilizada para diversos arranjos de caminho acústico, simulando medidores de vazão por ultrassom de um canal em várias condições para avaliar o desvio do fator de correção em escoamentos perturbados. Números de Reynolds de 1×10 4 a 2×10 6 são considerados. Ângulos de montagem dos transdutores de 0° a 180° são testados e posições axiais até 80D à jusante do obstáculo são avaliadas. Os resultados indicam que medidores de vazão por ultrassom são sensíveis aos efeitos de acidente de linha. O desvio do fator de correção mostra-se, em geral, consideravelmente maior que 2% em distâncias menores que 20D, conforme recomendado por alguns fabricantes e por leis. Não obstante, tais desvios podem atingir 0,01% em algumas configurações específicas, o que sugere que a medição de vazão por ultrassom pode ser melhorada pelo rearranjo do aparato em posição angular favorável e, principalmente, pela implementação de funções específicas para fatores de correção em condições perturbadas

RESUMO Simulações de aeroacústica computacional demandam uma quantidade considerável de tempo, o que torna complicada a realização de estudos paramétricos. O presente trabalho propõe uma metodologia viável para otimização aeroacústica. Através da análise numérica utilizando dinâmica dos fluidos computacional, foi estudada a aplicação de uma placa separadora desacoplada como método de controle passivo da esteira turbulenta de um cilindro e avaliou-se a irradiação de ruído causado pela interação do escoamento com ambos os corpos, empregando ferramentas de aeroacústica computacional baseadas no método de Ffowcs-Williams e Hawkings. Algumas abordagens distintas de metodologias de otimização de projeto foram aplicadas neste problema, com o objetivo de chegar a uma configuração otimizada que permita a redução do nível sonoro ao longe. Assim, utilizando uma ferramenta de otimização multidisciplinar, pode-se avaliar a capacidade de modelos heurísticos e a grande vantagem do emprego de algoritmos baseados em método de superfície de resposta quando aplicados em um problema não linear, pois requerem a avaliação de um menor número de alternativas para se obter um ponto ótimo. Além disso, foi possível identificar e agrupar os resultados em 5 clusters baseados em seus parâmetros geométricos, nível de pressão sonora global e o valor quadrático médio do coeficiente de arrasto, confirmando a eficiência da aplicação de placas separadoras longas desacopladas posicionadas próximas ao cilindro na estabilização da esteira turbulenta, enquanto que o posicionamento de placas acima de um espaçamento crítico aumentou o nível de pressão acústica irradiado devido à formação de vórtices no espaço entre o cilindro e a placa separadora. / Computational aeroacoustics simulations require a considerable amount of time, which makes the comparison of a large number of different geometric designs a difficult task. The goal of the present study is to provide a suitable methodology for aeroacustic optimization. By means of numerical analyses using computational fluid dynamics tools, the application of a detached splitter plate as a passive control method for the turbulent wake of a circular cylinder was investigated. The irradiation of noise caused by the interaction between the flow and both bodies was evaluated using computational aeroacoustics tools based on the Ffowcs-Williams and Hawkings method. Various design optimization methodologies were applied to this flow in order to achieve a possible optimal configuration, i.e., one which is capable of reducing the far field noise level without increasing the aerodynamic forces. Using a multidisciplinary optimization tool, it was possible to evaluate the behavior of heuristic optimization algorithms and the major advantage of algorithms based on response surface methods when applied to a nonlinear aeroacoustics problem, since they require a smaller number of calculated designs to reach the optimal configuration. In addition, it was possible to identify and group the outcomes into 5 clusters based on their geometric parameters, overall sound pressure level and drag coefficient, confirming the efficiency of the application of long detached splitter plates placed next to the cylinder in stabilizing the turbulent wake, whereas the positioning of splitter plates at a distance larger than a critical gap increased the overall sound pressure level radiated due to the formation of vortices in the gap.

Este trabalho apresenta a simulação numérica do escoamento turbulento em torno de uma turbina eólica de eixo vertical de pequeno porte, Savonius tipo helicoidal com torção de 180° nas pás. Com o intuito de avaliar a metodologia computacional empregada os resultados numéricos obtidos são comparados com os resultados experimental e numérico contidos no estado da arte. Também, compara-se o coeficiente de toque da turbina Savonius helicoidal com a turbina Savonius convencional. As simulações numéricas são baseadas no Método de Volumes Finitos, e para tal emprega-se o programa Fluent /Ansys versão 13.0 que resolve as equações da continuidade e as equações de Navier-Stokes com médias de Reynolds, juntamente com o modelo de turbulência . As simulações são desenvolvidas empregando diferentes malhas computacionais em estudos transientes, tridimensionais, com a turbina estacionária. A avaliação da qualidade da malha é realizada através do método de Índice de Convergência de Malha (GCI) o qual analisa o quão longe os resultados estão da solução assintótica para a malha utilizada. Após a análise da qualidade de malha, realizam-se simulações com a turbina em rotação as quais fazem uso da malha contendo uma região móvel possibilitando a imposição de uma velocidade angular ao rotor. O coeficiente de torque é obtido nas simulações e a partir dele calcula-se o coeficiente de potência. Além da análise do desempenho do rotor realiza-se uma análise qualitativa das características do escoamento sobre a turbina. A turbina Savonius helicoidal apresenta um valor de coeficiente de potência de 0,175 para a razão de velocidade de ponta de 0,58 considerando correção do efeito de bloqueio. Os resultados obtidos apresentam boa concordância com os resultados publicados por outros autores. / This dissertation presents the numerical simulation of the turbulent flow around of a small sized vertical axis wind turbine, consisting in a helical Savonius type with a 180° degree of blade twist. In order to evaluate the used methodology the obtained results are compared with the state of the art numerical and experimental data. It will be also presented the comparison between the torque coefficient of the conventional Savonius turbine and the helical Savonius turbine. The numerical simulations are based on the Finite Volume Method (FVM), using the commercial code Fluent/ANSYS version 13.0, which solves the continuity and Navier-Stokes through the Reynolds time-averaged methodology, including the turbulence model. The simulations are developed using different computational meshes for transient and three-dimensional studies with the stationary turbine. The evaluating the quality of the mesh is performed by of Grid Convergence Index (GCI) method which analyzes how far the results are the asymptotic solution to the mesh used. After the evaluation of the mesh quality, it was simulated a case considering the rotor motion using the moving mesh configuration, allowing the imposition of an angular velocity to the turbine. In the post-processing stage, it is possible to obtain the torque coefficient on the rotor shaft, allowing the calculation of the power coefficient for the turbine. In addition to the performance analysis, it is also made a qualitative analysis of the flow characteristics over the turbine rotor and in both cases presenting a good correspondence with the results in the literature. The helical Savonius turbine presents a value of power coefficient of 0.175 to a tip speed ratio of 0.58 whereas blocking effect correction.

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.