O presente trabalho tem como objetivo realizar um estudo comparativo entre os três métodos para obtenção das características aerodinâmicas de tabuleiros de pontes: analítico, experimental e numérico. O tema é motivado pelo crescimento das necessidades de conhecimento científico e tecnológico na área de desempenho aerodinâmico de seções típicas de pontes nas etapas mais iniciais de projeto. A precisão de uma ferramenta de previsão certamente acelera a convergência ao modelo final; porém, qualquer metodologia numérica deve ser extensivamente testada antes de ser utilizada nos escritórios de cálculo. Inicialmente mostra-se a lógica da análise de instabilidades provocadas pelo vento, resume-se a modelagem de tabuleiros de pontes para ensaios em túnel de vento e descreve-se o método numérico implementado no programa comercial de análise de pontes RM. Comparam-se dados disponíveis de testes de modelos reduzidos no túnel de vento com as estimativas de aproximação do Método da Vorticidade em Partículas (MVP). Confrontam-se resultados para casos clássicos, um estudo paramétrico, dados dos relatórios de três pontes Brasileiras já ensaiadas no Túnel de Vento Joaquim Blessmann do LAC – PPGEC/UFRGS: Paulicéia, Guamá e Roberto Marinho. Finalmente, considera-se o método experimental como preciso e julga-se o método numérico comparativamente, analisando a confiabilidade dos resultados e procurando definir a sua faixa de aplicabilidade. Examinando a variabilidade dos resultados, bastante baixa, definiu-se sua faixa de aplicação considerando os resultados obtidos como satisfatórios para projetos básicos e executivos de estruturas; necessitando, porém, para os casos especiais, de comprovação experimental em túnel de vento antes da sua execução. Exemplifica-se então o uso do método para alguns problemas de análises de seções. Os resultados e comparações corroboram verificações anteriores do método dos vórtices discretos para verdadeiras seções de pontes e mostram porque o seu uso vem se tornando cada vez mais aceitável. / The work has the objective to realize a comparative study between three methods to obtain bridge aerodynamic coefficients: analytic, experimental and numerical. The theme is motivated by the increasing need of knowledge about the aerodynamic behaviour of bridge decks in the early stages of the design process. The accuracy of a predictive tool will certainly accelerate the convergence to the final design. However, any numerical methodology must be extensively tested ascertain and hence validated before it becomes of current use by designers. Collaboration is made comparing wind tunnel data of reduced models and numerical estimations by a simplified approach. After a brief overview of wind loading and wind tunnel modeling, the bridge wind resistance design in the commercially available bridge design software package RM is shown. The Discrete Vortex Method and the implementation are shortly described. Initially, classical fluid problems are indicated together with the comparison of parametric cross sectional bridge shapes. After, three Brazilian large cablestayed bridge projects are presented as instance of analysis: Paulicéia, Guamá and Roberto Marinho. Bridges aerodynamic coefficients were experimentally obtained at the Boundary Layer Wind Tunnel Joaquim Blessmann at Universidade Federal do Rio Grande do Sul (UFRGS). Finally, experimental results are considered precise and those results evaluated with CFD are compared analyzing the liability and applicability. According to the results variability, quite low, we considered them appropriated for basic and final structural designs; needing, although, for special cases, of experimental probation in wind tunnel tests before the construction phase. Some CFD analysis problems of cross sections are then presented as application examples. The results and similitude corroborate earlier verifications of the method and this implementation, demonstrating why its use is becoming increasingly acceptable.

O objetivo deste trabalho é estudar algoritmos paralelos para a solução de problemas de Mecânica dos Sólidos, Mecânica dos Fluídos e Interação Fluido-Estrutura empregando o Método dos Elementos Finitos para uso em configurações de memória distribuída e compartilhada. Dois processos para o particionamento da estrutura de dados entre os processadores e divisão de tarefas foram desenvolvidos baseados na aplicação do método de particionamento em faixas e do método da bissecção coordenada recursiva não sobre a geometria da malha mas sim diretamente sobre o sistema de equações, através de reordenações nodais para minimização da largura da banda. Para ordenar a comunicação entre os processadores, foi desenvolvido um algoritmo simples e genérico baseado em uma ordenação circular e alternada que permite a organização eficiente dos processos mesmo em cenários nos quais cada processador precisa trocar dados com todos os demais. Os algoritmos selecionados foram todos do tipo iterativo, por sua adequabilidade ao paralelismo de memória distribuída. Foram desenvolvidos códigos paralelos para o Método dos Gradientes Conjugados utilizado em problemas de Mecânica dos Sólidos, para o esquema explícito de Taylor-Galerkin com um passo e iterações utilizado na simulação de escoamentos compressíveis em regime transônico e supersônico, para o esquema explícito de Taylor- Galerkin com 2 passos para simulação de escoamentos incompressíveis em regime subsônico e para interação fluído-estrutura usando o esquema explícito de dois passos para o fluído e o método implícito de Newmark no contexto do método de estabilização α-Generalizado para a estrutura, com acoplamento particionado. Numerosas configurações foram testadas com problemas tridimensionais utilizando elementos tetraédricos e hexaédricos em clusters temporários e permanentes, homogêneos e heterogêneos, com diferentes tamanhos de problemas, diferentes números de computadores e diferentes velocidades de rede. / Analysis and development of distributed memory parallel algorithms for the solution of Solid Mechanics, Fluid Mechanics and Fluid-Structure Interaction problems using the Finite Element Method is the main goal of this work. Two process for mesh partitioning and task division were developed, based in the Stripwise Partitioning and the Recursive Coordinate Bisection Methods, but applied not over the mesh geometry but over the resultant system of equations through a nodal ordering algorithm for system bandwidth minimization. To schedule the communication tasks in scenarios where each processor must exchange data with all others in the cluster, a simple and generic algorithm based in a circular an alternate ordering was developed. The algorithms selected to be parallelized were of iterative types due to their suitability for distributed memory parallelism. Parallel codes were developed for the Conjugate Gradient Method ( for Solid Mechanics analysis), for the explicit one-step scheme of Taylor-Galerkin method (for transonic and supersonic compressible flow analysis), for the two-step explicit scheme of Taylor-Galerkin method (for subsonic incompressible flow analysis) and for a Fluid-Structure Interaction algorithm using a coupling model based on a partitioned scheme. Explicit two-step scheme of Taylor-Galerkin were employed for the fluid and the implicit Newmark algorithm for the structure. Several configurations were tested for three-dimensional problems using tetrahedral and hexahedral elements in uniform and nonuniform clusters and grids, with several sizes of meshes, numbers of computers and network speeds.

BARBOSA, Charles Welton Ferreira. Uma análise comparativa de métodos aproximativos baseados em smoothed particle hydrodynamics para animação de fluidos. 2013. 59 f. : Dissertação (mestrado) - Universidade Federal do Ceará, Centro de Ciências, Departamento de Computação, Fortaleza-CE, 2013. / Submitted by guaracy araujo (guaraa3355@gmail.com) on 2016-06-14T18:05:54Z No. of bitstreams: 1 2013_dis_cwfbarbosa.pdf: 4026755 bytes, checksum: 3808122d6233e5855d6891854317c373 (MD5) / Approved for entry into archive by guaracy araujo (guaraa3355@gmail.com) on 2016-06-14T18:07:17Z (GMT) No. of bitstreams: 1 2013_dis_cwfbarbosa.pdf: 4026755 bytes, checksum: 3808122d6233e5855d6891854317c373 (MD5) / Made available in DSpace on 2016-06-14T18:07:17Z (GMT). No. of bitstreams: 1 2013_dis_cwfbarbosa.pdf: 4026755 bytes, checksum: 3808122d6233e5855d6891854317c373 (MD5) Previous issue date: 2013 / Animations of fluids, such as water or smoke, are used to add details to virtual games or films. Various methods exist for simulating fluids using particle systems, including the Smoothed Particle Hydrodynamics (SPH) method. The SPH method is a good choice for simulating fluids for its easy description and implementation. This work describes some methods based on SPH for general fluid simulation, and for the interaction between fluids and solids. These methods are analysed and compared on their realism, and a few weak points of each method are presented. This analysis may be used for the creation of better methods based on SPH, and for the realistic animation of fluids / Animações de fluidos, como a água ou fumaça, são utilizadas para introduzir detalhes em jogos virtuais ou filmes. Vários métodos existem para a simulação de fluidos utilizando sistemas de partículas e, em especial, utilizando o método Smoothed Particle Hydrodynamics, ou SPH. O método SPH é uma boa escolha para a simulação de fluidos por sua fácil descrição e implementação. Este trabalho apresenta alguns métodos baseados no SPH para a simulação de fluidos gerais, e para a interação entre fluidos e sólidos. Esses métodos são analisados e comparados relativos ao seu realismo, apresentando alguns pontos fracos de cada método. Essa análise pode ser utilizada para guiar a implementação de melhores métodos baseados no SPH, e para a animação realista de fluidos.

Esta dissertação investiga numericamente a influência dos efeitos viscosos e elásticos em escoamentos de materiais viscoplásticos no interior de uma cavidade dirigida. O modelo mecânico empregado é constituído pelas equações de conservação de massa e pelo princípio da quantidade de movimento linear, para fluidos incompressíveis, acoplado à equação constitutiva. Esta equação modifica o modelo viscoelástico de Oldroyd-B de modo a acomodar que os tempos de relaxação e retardo do material, bem como sua viscosidade viscoplástica, dependam das mudanças de sua microestrutura. A aproximação numérica do modelo emprega o método multi-campos de Galerkin mínimos-quadrados em termos do tensor de tensão extra, do vetor velocidade e do campo de pressão. Os resultados objetivam a determinção do tamanho e localização das regiões aparentemente não-escoadas do material, bem como sua deformação elástica, intensidade de tensão, e a sua vorticidade no interior da cavidade. Os resultados claramente indicam que o padrão do escoamento é fortemente influenciado pela variação dos efeitos elásticos (variação do tempo de relaxação adimensional, θ0 * ), viscosos (variação do índice de power-law, n) e cinemáticos (variação da velocidade adimensional, U* , do escoamento) no interior da cavidade. / This dissertation investigated numerically the influence of viscous and elastic effects on flows of viscoplastic materials within a lid-driven cavity. The mechanical model used is made up of mass and momentum balance equations, coupled with the constitutive equation. This equation modifies the viscoelastic Oldroyd-B model to accommodate both relaxation and retardation times, and viscosity function, dependent on the microstructure changes. Numerical approximations of the model make use a three-field Galerkin least squares method in terms of the extra stress tensor, velocity vector and pressure field. Computations focus on the determination of the size and position of apparently unyielded regions as well as the elastic deformation, stress intensity, and the vorticity within of the cavity. Results clearly indicate that the flow pattern is strongly influenced by the elastic (variation of the dimensionless relaxation time, θ0 * ), viscous (variation of the power-law index, n) and kinematic (variation of the dimensionless flow velocity, U* ) effects within the cavity.

É considerado o problema de encontrar descrições de Poisson (formulações Hamiltonianas generalizadas) associadas a modelos físicos. Aspectos básicos e aplicações dos sistemas de Poisson são explanados utilizando a linguagem da geometria diferencial. Sobre geometria diferencial, consta um capítulo com noções fundamentais. São consideradas as Mecânicas de Nambu e Birkho:ff e suas relações com a Mecânica Hamiltoniana generalizada. A questão da estábilidade é discutida do ponto de vista das formulações de Poisson. Os métodos existentes atualmente para derivação de estruturas Hamiltonianas generalizadas são expostos. Em particular, o processo de redução é estudado. Propõe-se uma abordagem dedutiva e inédita para construção de formulações de Poisson. O novo método é capaz de resolver (localmente) a questão de como encontrar descrições Hamiltonianas de sistemas dinâmicos com no máximo três dimensões. Nos casos tridimensionais nos quais é conhecida uma superfície à qual as trajetórias são sempre tangentes, a nova estratégia reduz esta questão à solução de uma equação diferencial parcial de primeira ordem linear. Deste modo demonstra-se a existência (local) genérica de estruturas de Poisson para sistemas tridimensionais. O caso tridimensional é analizado com detalhe, par ticularmente no concernente à invari ância conforme da identid ade de Jacobi nesta dimensionalidade. A abordagem tratada nesta dissertação é aplicada a vários sistemas tridimensionais de interesse. / The problem of finding Poisson descriptions (generalized Hamiltonian formulations) assoei ateei with physical models is considered. The basic features anel aplications of Poisson systems are explained in the language of differential geometry. One chapter is included with the fundamental notions on differential geometry. The Nambu anel Birkhoff's Mechanics anel their relationship with the generalized Hamiltonian Mechanics are considered. The question of stability is discussed from the point of view of the Poisson formulations. The currently existing methods for derivation of generalized Hamiltonian structures are reviewed. Particularly, the reduction process is analized. A deductive approach is proposed for the construction of Poisson formu lations. The new method can solve (locally) the question of how to finei Hamiltonian descriptions of dynamical systems in, at most, three dimensions. When a surface to wich the motion is always tangent is known , in three dimensions the new approach reduces the problem to the solution of a linear partia! differential equation of first order. This demonstrates the general existence (local) of Poisson structures for tridimensional systems. The tridimensional case is analized in detail, particularly in what concerns the conformai invariance of the Jacobi identity in this dimensionality. The approach proposed in this dissertation is applied to various tridimensional systems of interest.

Esta dissertação tem por objetivo a implementação de um método adaptativo anisotrópico de malhas e seu acoplamento a um programa de simulação de escoamentos de fluidos compressíveis numa faixa ampla de número de Mach utilizando o Método dos Elementos Finitos (Bubnov-Galerkin). Mais especificamente, implementa-se um código para simulação de escoamentos compressíveis utilizando o método das linhas ou direções características (CBS), o qual possui boa performance para uma amplo espectro de números de Mach. A discretização espacial é feita utilizando-se elementos finitos triangulares lineares e tetraédricos lineares para problemas bidimensionais e tridimensionais, respectivamente. A discretização temporal é feita por uma expansão em séries de Taylor. Implementa-se uma adaptação anisotrópica com capacidade de refinamento e desrefinamento do domínio computacional para se obter soluções precisas a baixos custos computacionais, onde a estimativa de erro é conduzida por uma métrica Riemanniana, a qual permite o tratamento contínuo do erro através de um tensor de erro anisotrópico. O código é modificado para possibilitar o uso de arquiteturas de memória compartilhada a fim de aumentar a velocidade de processamento da simulação. Por fim, diversos estudos, tanto para problemas bidimensionais quanto tridimensionais, são realizados, validando a metodologia e realizando estudos adicionais em problemas mais complexos. / The objectives of this dissertation are the implementation of an anisotropic mesh adaptation procedure and its coupling with a compressible flow simulation code with large Mach range using the Finite Element Method (Bubnov-Galerkin). Specifically, a computer code using the Characteristic Based Split method called CBS, which has a good performance for a wide range of Mach numbers, is implemented. The spatial discretization is performed using linear triangular and linear tetrahedral finite elements for two and three dimensional cases, respectively. The temporal discretization is done with Taylor series. An anisotropic mesh adaption technique, based on a Riemann metric to evaluate the errors (which allows continuous error analysis through an anisotropic error tensor evaluation), including mesh refinement and mesh coarsening procedures is also implemented looking for accurate and computationally cheaper solutions. The code is able to use shared memory architecture in order to reduce computer process time. Finally, two-dimensional and three-dimensional studies are performed in order to validate the proposed methodology and additional complex tests are also analysed.

O presente trabalho realiza um estudo do desempenho aerodinâmico da turbina eólica de duas pás e 10 m de diâmetro da UAE Phase VI, testada pelo NREL no túnel de vento NASA Ames, através do uso de dinâmica dos fluidos computacional (CFD). Para tanto são apresentados conceitos fundamentais para análise da performance do aerogerador, a metodologia aplicada para a análise numérica e os resultados obtidos. O modelo geométrico da turbina foi inserido em ambiente computacional num domínio com as mesmas dimensões da seção de testes do túnel de vento. A esse domínio foi empregada uma malha de elementos tetraédricos. A modelagem numérica e as simulações foram efetuadas através do código comercial ANSYS FLUENT 13.0 e utilizaram as equações médias de Reynolds e modelos de turbulência k  SST e Spalart-Allmaras. A turbina simulada possui ângulo de giro (yaw) e de cone de 0°. Os casos simulados foram de turbina com ângulos de passo das pás de 0° e 3°, para velocidades de ar na entrada do domínio de 5 m/s, 7 m/s e 9 m/s. As simulações foram feitas em regime transiente utilizando o método de malha móvel. Dados provenientes da simulação numérica são comparados aos dados experimentais de empuxo, torque e potência, divulgados pelo laboratório. Os modelos de turbulência testados apresentaram boa concordância com os resultados experimentais de empuxo. O torque e a potência foram bem previstos para velocidades de entrada do domínio baixas, mas foram subdimensionados para velocidades mais altas. A separação da camada limite foi prevista por ambos os modelos. Na comparação com resultados obtidos com outros autores para o mesmo caso, o início de descolamento da camada limite se deu antes do previsto. A metodologia adotada nas simulações numéricas realizadas mostrou-se adequada e representativa nas análises desse estudo. / This study makes a CFD analysis of the aerodynamic performance of the NREL UAE Phase VI two-bladed wind turbine, which have been tested in the NASA Ames wind tunnel. Fundamental concepts on the evaluation of turbine performance, the applied methodology for the numerical analysis and the results are presented. A numerical model has been inserted in a computational environment that has the same dimensions as the real wind tunnel section, and a tetrahedral mesh has been created to fill this domain. Numerical modeling and the simulations have been performed using the ANSYS FLUENT 13.0 commercial code, making use of the RANS equations and the k  SST and Spalart-Allmaras turbulence models. The simulated turbine has 0 degrees of cone and yaw angle. Simulations were performed in unsteady state using the moving mesh technique. Results are compared to experimental data regarding to thrust, torque and power. The k  SST and Spalart-Allmaras turbulence models have shown good agreement with experimental values of thrust. Torque and power have been adequately predicted to lower free flow velocities, but have been underpredicted to higher velocities. Separation of the boundary layer has been adequately predicted by both models, but the beginning of the separation occurred earlier than expected in comparison with results of other authors to the same case. The methodology used in the numerical simulations proved to be adequate and representative in this study.

Escoamentos sem inércia de fluido elasto-viscoplástico, dentro de uma cavidade, são numericamente analisados. As soluções visam compreender a influência dos efeitos viscosos e elásticos na topologia de superfícies de escoamento. Assumindo-se que o colapso da microestrutura do material é instantâneo, o modelo mecânico é constituído pelas equações governantes de massa e momentum para fluidos incompressíveis, associado a uma equação hiperbólica para o tensor tensão extra, baseado na equação do modelo Oldroyd-B (Nassar et al, 2011). A principal característica do modelo é considerar a viscosidade e o tempo de relaxação como função da taxa de deformação, permitindo a pseudoplasticidade de viscosidade e restringindo os efeitos elásticos para as regiões não deformadas do material. As simulações numéricas são realizadas através do método de Galerkin mínimos quadrados a três campos: tensor tensão extra, pressão e velocidade. Os resultados mostram que as superfícies de escoamento do material são fortemente influenciadas pela ação combinada entre os efeitos elásticos e viscosos, estando em conformidade com a recente visualização experimental dos fluxos elasto-viscoplásticos. / Elasto viscoplastic uid ows without inertia, within a cavity, are numerically analyzed. The solutions aim to understand the in uence of viscous and elastic e ects on the topology of yield surfaces. Assuming that the collapse of the material microstructure is instantaneous, the mechanical model consists of the governing equations of mass and momentum for incompressible uids, associated with a hyperbolic equation for the extra stress tensor, based on the equation of the Oldroyd-B model (Nassar et al, 2011). The main feature of model is to consider the viscosity and the relaxation time as a function of shear rate, allowing the shear-thinning of viscosity and restricting elastic e ects for regions not deformed material. Numerical simulations are performed by the method of Galerkin Least Squares to three elds: extra stress tensor, pressure and velocity. The results show that the yield surfaces of material are strongly in uenced by the combined action between the elastic and viscous e ects, complying with the recent experimental visualization of elasto-viscoplastic ows.

O presente trabalho destina-se ao estudo numérico da geometria de tubos de seção elíptica que facilite a transferência de calor adimensional e diminua a queda de pressão adimensional (Δ̃) sofrida pelo escoamento. O método aplicado é o Design Construtal, que visa determinar a geometria que apresentará a menor resistência ao escoamento, ou seja, busca-se determinar a razão de aspecto da elipse (=⁄) que favorece a transferência de calor e diminui a queda de pressão do escoamento. O fluido empregado neste estudo apresenta características de viscoplasticidade. A relação entre a tensão cisalhante e a taxa de deformação obedece ao modelo de Herschel-Bulkley modificado. Considera-se que o escoamento é incompressível, laminar, bidimensional, externo e ocorre em regime permanente. A solução numérica do problema proposto é realizada com um código comercial baseado no método dos volumes finitos. É investigada a influência do índice de potência, , sobre a seção elíptica que facilita o escoamento e, para isso, este índice é variado de 0,4 a 1. A influência dos números de Reynolds (√), Herschel-Bulkley (√) e Prandtl (√) sobre o comportamento do escoamento também é avaliada. √ é variado de 1 a 40, √ é variado de 1 a 100 e √ é variado de 0,1 a 100 Os resultados mostram que, para um escoamento com √=1, √=1 e √=1, o aumento do índice de potência influencia negativamente na transferência de calor adimensional e a seção elíptica, que maximiza esta transferência de calor adimensional, tende a ser mais alongada na direção do escoamento. Já e influenciam positivamente na transferência de calor adimensional. Para um escoamento com √=1, √=1, =0,4 conclui-se que com o aumento de a razão de aspecto ótima (q,opt), do ponto de vista térmico, diminui. Quando é considerado um escoamento com √=1, √=1, =0,4 conclui-se que q,opt diminui com o aumento de , ou seja, a elipse torna-se mais alongada no sentido do escoamento. A variação de √ em um escoamento com √=1, √=1, =0,4 mostra que o aumento deste parâmetro acarreta em aumento da taxa de transferência adimensional e de Δ̃. / The present work is aimed at the numerical study of the geometry of elliptic section tubes that facilitates the dimensionless heat transfer and decreases the dimensionless pressure drop (Δ̃) suffered by the flow. The applied method is the Construtal Design, which aims to determine the geometry that will present the least resistance to the flow, that is, to determine the aspect ratio of the ellipse (=⁄) that favors heat transfer and decreases the flow pressure drop. The fluid used in this study has viscoplasticity characteristics. The relationship between shear stress and strain rate follows the modified Herschel-Bulkley model. It is considered that the flow is incompressible, laminar, two-dimensional, external and occurs in steady state. The numerical solution of the proposed problem is carried out with a commercial code based on the finite volume method. The influence of the power index, n, on the elliptical section facilitating the flow is investigated, and for this, the index is varied from 0.4 to 1. The influence of the Reynolds number (√), Herschel-Bulkley number (√) and Prandtl number (√) on the flow behavior is also evaluated √ is varied from 1 to 40, √ is varied from 1 to 100 and √ is varied from 0.1 to 100. The results show that for a flow with √=1, √=1 and √=1, the increase of the power index negatively influences the dimensionless heat transfer and the elliptic section, which maximizes this dimensionless heat transfer, tends to be more elongated in the direction of flow. Already √ and √ influence positively the dimensionless heat transfer. For a flow with √=1, √=1, =0.4 it is concluded that with the increase of √ the optimum aspect ratio (q,opt), from the thermal point of view, decreases. When a flow is considered with √=1, √=1, =0.4 it is concluded that q,opt decreases with the increase of √, that is, ellipse becomes more elongated in the flow direction. The variation of √ in a flow with √=1, √=1, =0.4 shows that the increase of this parameter causes an increase of the dimensionless transfer rate and Δ̃.

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.