Currently, most micrositing techniques aim to maximize annual energy production (AEP) or minimize cost of energy (COE) with no direct regard to revenue. This research study developed a method that utilizes the seasonal electricity price and wind data to microsite wind farms in terms of profitability. To accomplish this, six candidate wind farms with differing layouts and spacing were selected at a given location. They were then simulated using a wake modeling software to produce expected power outputs at different wind speeds, wind directions, and turbulence intensities. By interpolating the power output tables with wind data, a power time-series was created for each wind farm over a multi-year period. Electrical price was then incorporated with the power time-series to produce a revenue time-series of the revenue produced at each hour over the same time-period. Each relative wind farm was then rotated in increments to evaluate new candidate wind farms and revenue totals. This method is site specific and results may differ dependent on location and seasonal correlation between wind and electrical data. Overall, the method looks to exploit a different approach to the micrositing problem.

This dissertation discusses the potential for using a conventional main lateral-force resisting system, combined with the reserve strength in the gravity framing, and or auxiliary collapse-inhibiting mechanisms deployed throughout the building, or enhanced shear tab connections, to provide adequate serviceability performance and collapse safety for seismic and wind hazards in the central and eastern United States. While the proposed concept is likely applicable to building structures of all materials, the focus of this study is on structural steel-frame buildings using either non-ductile moment frames with fully-restrained flange welded connections not specifically detailed for seismic resistance, or ductile moment frames with reduced beam section connections designed for moderate seismic demands. The research shows that collapse prevention systems were effective at reducing the conditional probability of seismic collapse during Maximum Considered Earthquake (MCE) level ground motions, and at lowering the seismic and wind collapse risk of a building with moment frames not specifically detailed for seismic resistance. Reserve lateral strength in gravity framing, including the shear tab connections was a significant factor. The pattern of collapse prevention component failure depended on the type of loading, archetype building, and type of collapse prevention system, but most story collapse mechanisms formed in the lower stories of the building. Collapse prevention devices usually did not change the story failure mechanism of the building. Collapse prevention systems with energy dissipation devices contributed to a significant reduction in both repair cost and downtime. Resilience contour plots showed that reserve lateral strength in the gravity framing was effective at reducing recovery time, but less effective at reducing the associated economic losses. A conventional lateral force resisting system or a collapse prevention system with a highly ductile moment frame would be required for regions of higher seismicity or exposed to high hurricane wind speeds, but buildings with collapse prevention systems were adequate for many regions in the central and eastern United States. / Ph. D.

In this thesis, the tonal noise produced by flow over perforated plates at oblique angles of incidence is studied experimentally. A two-dimensional model of a perforated plate is used, where the circular holes of a typical perforated plate are replaced by a series of long rectangular Aluminum slats with an adjustable gap width between them. The slats are 3.175 mm thick and the gap width between them is set to 3.175 mm, 6.35 mm, and 12.7 mm. This simplified model is mounted at the exit of an open-loop wind tunnel and tested at angles of incidence of 0° to 40° and flow velocities of 0 to 30 m/s. An angle of 0° is defined as flow parallel to the plate. The acoustic response is studied using microphone measurements, and flow visualization is done using particle image velocimetry. The effect of the angle of incidence, flow velocity, gap width, and streamwise position are investigated. The flow visualization reveals that tonal noise is produced by the periodic shedding and impingement of vortices at the trailing edge of the gaps. Vortices form in the unstable free shear layer originating at the leading edge of the gap and impinge on the downstream side of the gap. At the downstream corner, these vortices separate into vortex pairs, consisting of one positively rotating and one negatively rotating vortex. These vortices are shed periodically, leading to the production of tonal noise at the shedding frequency. The effect of the angle of incidence is investigated by changing the angle of the plate with respect to the flow. For a given gap width, tones are produced only for a specific range of angles. Depending on the plate geometry, this range of angles is typically around 5° to 30°. Within this range of angles, the free shear layer impinges on the downstream side of the gap. For angles which are too small or too large, the free shear layer misses this downstream side and tones are not produced. For a larger gap width, tones are produced at smaller angles of incidence. Similarly, for a given plate geometry, there is a preferred range of flow velocities at which tonal noise is produced. The velocity at which the free shear layer is the most unstable at the tone frequency produces the strongest vortices and the loudest tones. The optimal velocity is lower for larger gap widths. Finally, it is found that the magnitude of the produced tones increases in the streamwise direction over repeated gaps along the length of the plate. This is due to the local flow conditions changing in the streamwise direction, only reaching the optimal conditions after a certain length of the plate. / Thesis / Master of Applied Science (MASc)

This thesis describes a series of 1:75 scale wind tunnel experiments investigating the wind flow over, and through, three different forest models and the resultant wind loading on individual model trees. The experiments were designed to lead to a quantitative assessment of the wind stability of the particular forest arrangements and also to permit a study of the coherent gust structures in the flow. Forest canopy flow is dominated by a plane mixing layer flow regime with a shear layer close to the canopy top. It has been confirmed that data can be correlated usefully in terms of a shear length, Ls, related to the form of this shear layer. Frequency analysis has confirmed that the flow structures have the same frequency as the swaying of the tallest trees in each forest. A mechanism is proposed whereby upstream turbulence induces swaying of trees at and near the upwind edge region of the forest, which in turn perturbs the air in the unstable shear layer. This leads to a roll-up of the shear layer and the creation of coherent flow structures. Conditional sampling of the gust structures, using wavelet analysis, has also supported the theory of a plane mixing layer type flow. An eddy-pair structure was revealed, the arrangement of which accounts for the intermittent strong downward sweeps of air into the canopy that have been reported by many observers. The large downward sweep of air was also shown to be responsible for the highest bending moments experienced by individual trees. Assessments of the different forest formations showed that in a forest consisting of a 50/50 mix of 200 mm and 100 mm model trees, gusts did not penetrate the lower forest. This arrangement should improve the protection of younger trees and may be worth investigating in field trials.

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 objetivo deste trabalho é estudar e desenvolver procedimentos numéricos adequados para a análise de problemas da Engenharia do Vento Computacional (EVC). O escoamento é analisado a partir das equações de Navier-Stokes para um fluido Newtoniano e de uma equação de conservação de massa considerando a hipótese de pseudo-compressibilidade, ambas em um processo isotérmico. Na presença de escoamentos turbulentos emprega-se a Simulação de Grandes Escalas (“LES”) com os modelos clássico e dinâmico de Smagorinsky para as escalas inferiores à resolução da malha. Dois modelos numéricos de Taylor-Galerkin para a análise do escoamento são estudados: o esquema explícito de dois passos e o esquema explícito-iterativo. O Método dos Elementos Finitos (MEF) é empregado para a discretização do domínio espacial utilizando o elemento hexaédrico trilinear isoparamétrico com integração reduzida das matrizes em nível de elemento. Em problemas envolvendo efeitos de interação fluido-estrutura emprega-se um esquema de acoplamento particionado com características superiores de conservação, permitindo, inclusive, o uso de subciclos entre as análises do fluido e da estrutura e de malhas não compatíveis na interface. A estrutura é considerada como um corpo deformável constituído de um material elástico linear com a presença de nãolinearidade geométrica. O MEF é também usado para a discretização da estrutura, empregando-se para tanto o elemento hexaédrico trilinear isoparamétrico com integração reduzida e controle de modos espúrios. A equação de equilíbrio dinâmico é integrada no tempo utilizando o método implícito de Newmark no contexto do método de estabilização α- Generalizado. Na presença de estruturas deformáveis, o escoamento é descrito através de uma formulação arbitrária Lagrangeana-Euleriana (ALE). Ao final, comparações com exemplos numéricos e experimentais são apresentadas para demonstrar a viabilidade dos algoritmos desenvolvidos, seguindo-se com as conclusões do trabalho e as sugestões para trabalhos futuros. / Analysis and development of numerical tools to simulate Computational Wind Engineering (CWE) problems is the main goal of the present work. The isothermal flow is analyzed using the Navier-Stokes equations for viscous fluids and a mass conservation equation obtained according to the pseudo-compressibility assumption. Turbulent flows are simulated employing Large Eddy Simulation (LES) with the classical and dynamic Smagorinsky’s models for subgrid scales. Two Taylor-Galerkin models for the flow analysis are investigated: the explicit two-step scheme and the explicit-iterative scheme. The Finite Element Method (MEF) is employed for spatial discretizations using the eight-node hexahedrical isoparametric element with one-point quadrature. Fluid-structure interaction problems are analyzed with a coupling model based on a conservative partitioned scheme. The Finite Element Method (MEF) is employed for spatial discretizations using the eight-node hexahedrical isoparametric element with one-point quadrature. Fluid-structure interaction problems are analyzed with a coupling model based on a conservative partitioned scheme. Subcycling and nonmatching meshes for independent discretizations of the fluid and structure domains are also available. The structure is considered as a deformable body constituted by a linear elastic material with geometrically nonlinear effects. The FEM is used for the spatial discretization of the structure as well. Eight-node hexahedrical isoparametric elements with one-point quadrature and hourglass control are adopted in this process. The implicit Newmark algorithm within the framework of the α-Generalized method is employed for the numerical integration of the dynamic equilibrium equation. An arbitrary Lagrangean-Eulerian (ALE) description is adopted for the kinematic description of the flow when deformable structures are analyzed. Numerical and experimental examples are simulated in order to demonstrate the accuracy of the developed algorithms. Concluding remarks and suggestions for future works are pointed out in the last chapter of the present work.

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