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