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

O presente trabalho 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.

O presente trabalho apresenta uma metodologia de simulação numérica de perfis aerodinâmicos bidimensionais com foco na utilização para o projeto e otimização de pás e rotores de pequenas turbinas eólicas de eixo horizontal, bem como o emprego desses métodos em simulações nas quais efeitos de alta solidez do rotor e baixos números de Reynolds são avaliados. Essa metodologia inclui geração de malhas, seleção de métodos numéricos e validação, tendo as escolhas sido guiadas pelas práticas mais bem sucedidas na simulação de perfis aerodinâmicos, e foi aplicada na simulação dos aerofólios NACA 0012, S809 e SD7062. O código comercial ANSYS Fluent foi utilizado em todas as simulações. Na simulação de aerofólios isolados a altos números de Reynolds dos perfis NACA 0012 e S809, o modelo Transition SST (γ-Reθ) apresentou resultados mais próximos a dados experimentais do que aqueles apresentados pelo modelo k-ω SST para CL e CD, além de produzir resultados para CP que mostraram boa precisão quando comparados aos mesmos dados experimentais. Resultados de CL, CD, CF e CP são apresentados para 20 diferentes condições de operação às quais o perfil SD7062 foi submetido, com números de Reynolds variando entre 25.000 e 125.000. As distribuições dos dois últimos coeficientes sobre os dorsos do aerofólio evidenciam com clareza a presença e magnitude da bolha de separação laminar. Os coeficientes de sustentação e arrasto mostram o impacto negativo da presença da bolha nessa faixa de números de Reynolds. Além disso, nos casos simulados, o arrasto aumenta em função da diminuição do Re. Um design de pá produzido com o auxílio do código de otimização SWRDC, baseado em algoritmos genéticos, é apresentado. Três seções ao longo da envergadura dessa pá foram simuladas em uma bateria de 45 simulações, sob diversas condições de operação em função de solidez, ângulo de ataque e razão de velocidade de ponta de pá. Esses resultados mostram que a bolha de separação laminar se move na direção do bordo de ataque com o aumento da solidez, do ângulo de ataque e da TSR. Além disso, distribuições do CP mostram aumento de pressão em ambos os dorsos do perfil quando submetido aos efeitos da solidez, embora esses efeitos tenham sido responsáveis por um aumento na relação CL/CD nos casos estudados. / This thesis presents a methodology of two-dimensional airfoil simulation focusing on its application on the design and optimization of blades and rotors of small horizontal axis wind turbines, and its application in a set of numerical simulations involving high rotor solidity and low-Re effects. This methodology includes grid generation, selection of numerical methods and validation, reflecting the most successful practices in airfoil simulation, and was applied in the simulation of the NACA 0012, S809 and SD7062 airfoils. The ANSYS Fluent commercial code was used in all simulations. Results for the isolated NACA 0012 and S809 airfoils at high Reynolds numbers show that the Transition SST (γ-Reθ) turbulence model produces results closer to experimental data than those yielded by the SST k-ω model for CL and CD, having also produced CP plots that show good agreement to the same experimental data. Plots of CL, CD, CF and CP for the SD7062 airfoil are presented, for simulations at 20 different operating conditions. The CF and CP distributions evidence the negative impact of the laminar separation bubble in the range of Reynolds numbers evaluated. Results show that, for Re between 25,000 and 125,000, drag increases with decreasing Re. A blade design generated using the SWRDC optimization code, based on genetic algorithms, is presented. Three sections of the resulting blade shape were selected and were tested in a set of 45 simulations, under an array of operating conditions defined by solidity, angle of attack and TSR. Results show that the laminar separation bubble moves towards the leading edge with increasing solidity, angle of attack and TSR. Furthermore, CP plots show an increase in pressure on both surfaces when the airfoil is subject to solidity effects, although these effects show an increase in the lift-to-drag ratio at the conditions evaluated.

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.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / This dissertation addresses the analysis and control of permanent magnet synchronous generators (PMSG) found in wind energy conversion system. The operation modes of the wind turbine the current and voltage constraints of a PWM rectifier are defined and a methodology to determine the generator current is presented. The generator currents obtained guarantee the maximization of wind energy extraction as well as the minimization of the PMSG and rectifier losses. In addition, a space-vector modulation technique in overmodulation region is reviewed in detail and applied aiming: (i) to maximize the generator power; (ii) to utilize the entire PWM rectifier DC link voltage resource, and (iii) to obtain a smooth transition to six-step operation. Furthermore, a single current controller is designed for all operating modes being capable to avoid the controller states overload as well as the undesired effects due to the operation in the overmodulation region. The stability analysis of the closed loop system is investigated by means of linear matrix inequality conditions to estimate a region of operation with guaranteed stability. Simulation results are presented to demonstrate the performance of the proposed system. / Esta dissertação trata da análise e do controle de geradores sincronos a íımã permanente usados em sistema de conversão de energia eólica. Os modos de operação da turbina eólica sob restrição de tensão e corrente são definidos e uma metodologia para determinar as correntes do gerador ´e apresentada. As correntes obtidas do gerador garantem a maximização da potência extraída do vento bem como a minimização das perdas do conjunto gerador e retificador PWM. Além disso, a técnica de modulação vetorial na região de sobremodulação é revisada em detalhes e aplicada com objetivo de: (i) maximizar a potência do gerador, (ii) utilizar todo o recurso do barramento CC disponível pelo retificador PWM e (iii) obter uma transição suave para operação do retificador para o modo de seis-pulsos. Uma técnica de controle vetorial é apresentada e comparada. Um único controlador de corrente é projetado para todos os modos operacionais, sendo capaz de evitar tanto sobrecarga dos estados do controlador quanto os efeitos indesejáveis provenientes da operação na região de sobremodulação. A análise da estabilidade do sistema em malha fechada é investigada por meio de desigualdades matriciais lineares para determinar uma região de estabilidade garantida. Por fim, resultados de simulação são apresentados para demonstrar o bom desempenho do sistema proposto.

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

O presente trabalho 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.

O presente trabalho apresenta uma metodologia de simulação numérica de perfis aerodinâmicos bidimensionais com foco na utilização para o projeto e otimização de pás e rotores de pequenas turbinas eólicas de eixo horizontal, bem como o emprego desses métodos em simulações nas quais efeitos de alta solidez do rotor e baixos números de Reynolds são avaliados. Essa metodologia inclui geração de malhas, seleção de métodos numéricos e validação, tendo as escolhas sido guiadas pelas práticas mais bem sucedidas na simulação de perfis aerodinâmicos, e foi aplicada na simulação dos aerofólios NACA 0012, S809 e SD7062. O código comercial ANSYS Fluent foi utilizado em todas as simulações. Na simulação de aerofólios isolados a altos números de Reynolds dos perfis NACA 0012 e S809, o modelo Transition SST (γ-Reθ) apresentou resultados mais próximos a dados experimentais do que aqueles apresentados pelo modelo k-ω SST para CL e CD, além de produzir resultados para CP que mostraram boa precisão quando comparados aos mesmos dados experimentais. Resultados de CL, CD, CF e CP são apresentados para 20 diferentes condições de operação às quais o perfil SD7062 foi submetido, com números de Reynolds variando entre 25.000 e 125.000. As distribuições dos dois últimos coeficientes sobre os dorsos do aerofólio evidenciam com clareza a presença e magnitude da bolha de separação laminar. Os coeficientes de sustentação e arrasto mostram o impacto negativo da presença da bolha nessa faixa de números de Reynolds. Além disso, nos casos simulados, o arrasto aumenta em função da diminuição do Re. Um design de pá produzido com o auxílio do código de otimização SWRDC, baseado em algoritmos genéticos, é apresentado. Três seções ao longo da envergadura dessa pá foram simuladas em uma bateria de 45 simulações, sob diversas condições de operação em função de solidez, ângulo de ataque e razão de velocidade de ponta de pá. Esses resultados mostram que a bolha de separação laminar se move na direção do bordo de ataque com o aumento da solidez, do ângulo de ataque e da TSR. Além disso, distribuições do CP mostram aumento de pressão em ambos os dorsos do perfil quando submetido aos efeitos da solidez, embora esses efeitos tenham sido responsáveis por um aumento na relação CL/CD nos casos estudados. / This thesis presents a methodology of two-dimensional airfoil simulation focusing on its application on the design and optimization of blades and rotors of small horizontal axis wind turbines, and its application in a set of numerical simulations involving high rotor solidity and low-Re effects. This methodology includes grid generation, selection of numerical methods and validation, reflecting the most successful practices in airfoil simulation, and was applied in the simulation of the NACA 0012, S809 and SD7062 airfoils. The ANSYS Fluent commercial code was used in all simulations. Results for the isolated NACA 0012 and S809 airfoils at high Reynolds numbers show that the Transition SST (γ-Reθ) turbulence model produces results closer to experimental data than those yielded by the SST k-ω model for CL and CD, having also produced CP plots that show good agreement to the same experimental data. Plots of CL, CD, CF and CP for the SD7062 airfoil are presented, for simulations at 20 different operating conditions. The CF and CP distributions evidence the negative impact of the laminar separation bubble in the range of Reynolds numbers evaluated. Results show that, for Re between 25,000 and 125,000, drag increases with decreasing Re. A blade design generated using the SWRDC optimization code, based on genetic algorithms, is presented. Three sections of the resulting blade shape were selected and were tested in a set of 45 simulations, under an array of operating conditions defined by solidity, angle of attack and TSR. Results show that the laminar separation bubble moves towards the leading edge with increasing solidity, angle of attack and TSR. Furthermore, CP plots show an increase in pressure on both surfaces when the airfoil is subject to solidity effects, although these effects show an increase in the lift-to-drag ratio at the conditions evaluated.

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

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.