Nesta dissertação, estuda-se a transição de esteira no escoamento ao redor de um aerofólio NACA 0012 com ângulos de ataque de zero dez e vinte graus. Dois casos são considerados: fólio isolado e fólio alinhado com um cilindro. Nas duas configurações, analisa-se a estabilidade linear em relação a perturbações tridimensionais. Tais perturbações foram estudadas usando a teoria de Floquet para um conjunto de números de Reynolds e ângulos de ataques. O escoamento base é calculado usando o método de elementos finitos espectrais para a discretização espacial. Dos resultados de estabilidade no caso do aerofólio isolado, pode-se observar dois picos de instabilidade com diferentes comprimentos de onda na envergadura. A simetria dos modos instáveis é também apresentada. Um dos modos instáveis presente na esteira do aerofólio isolado foi também observado no caso do cilindro alinhado com o fólio, enquanto o outro modo foi suprimido em tal geometria / Study of hydrodynamic stability of the flow around a cylinder aligned with on airfoil

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.

This paper will analyze the effects of using hybrid optimization methods for optimizing objective functions that are determined by computational fluid dynamics solvers for compressible viscous flow for optimal design of airfoils. Previous studies on this topic by the authors had examined the application of deterministic optimization methods and stochastic optimization methods such as Simulated Annealing and Simultaneous Perturbation Stochastic Analysis (SPSA). The studies indicated that SPSA method has a greater or equal efficiency as compared with SA method in reaching optimal airfoil designs for the design problem in question. However, in some situations SPSA method has a tendency to demonstrate an oscillatory behavior in the vicinity of a local optima. To overcome this tendency, a hybrid method designed to take full advantage of SPSA’s high rate of reduction of the objective function at the inception of the design process to drive the design cycles towards the optimal zone at first, and then combining with other methods to perform the final stages of the convergence towards the optimal solutions is considered. SPSA method has been combined with the gradient-based Broydon-Fletcher-Goldfarb-Shanno (BFGS) method as well as Simulated Annealing method for the transonic inverse airfoil design problem that is concerned with the specification of a target airfoil surface pressure distribution and starting from an initial guess of an airfoil shape, the target airfoil shape is reached by way of minimization of a quantity that depends on the difference between the target and current airfoil surface pressure distribution. For a typical transonic flow test case, the effects of using hybrid optimization techniques such as SPSA+BFGS and SPSA+SA as opposed to using SPSA alone can be seen in Figure 1. After 800 design cycles using SPSA, the hybrid SPSA+SA method took 2521 function evaluations of SA while the SPSA+BFGS method took 271 function evaluations to reach similar values which are much better than that reached by using SPSA alone in the entire minimization process. Results indicate that both of the two hybrid methods have capability to find a global optimum more efficiently than the SPSA method. The paper will address issues related to hybridization and its impact on the optimal airfoil shape designs in various contexts. / Singapore-MIT Alliance (SMA)

In this thesis we describe the development of Reynolds- averaged Navier Stokes code for the flow past two- dimensional configuration. Particularly, emphasis has been laid on the study of low Reynolds number airfoil aerodynamics. The thesis consists of five chapters covering the back ground history, problem formulation, method of solution and discussion of the results and conclusion. Chapter I deals with a detailed background history of low Reynolds number aerodynamics, problem associated with it, state of the art, its importance in practical applications in aircraft industries. Chapter II describes the mathematical model of the flow physics and various levels of approximations. Also it gives an account of complexity of the equations at low Reynolds number regarding flow separation, transition and reattachment. Chapter III describes method of solution, numerical algorithm developed, description of various upwind schemes, grid system, finite volume discrieti-zation of the governing equations described in Chapter II. Chapter IV describes the application of the newly developed Navier Stokes code for the test cases from GAMM Workshop proceedings. Also it describes validation of the code for Euler solutions, Blasius solution for the flow past flat plate and compressible Navier Stokes solution for the flow past NACA 0012 Airfoil at low Reynolds number. Chapter V describes the application of the Navier Stokes code for the more test cases of current practical interest . In this chapter laminar separation bubble characteristics are investigated in detail regarding formation, growth and shedding in an unsteady environment. Finally the conclusion is drawn regarding the robustness of the newly developed code in predicting the airfoil aerodynamic characteristics at low Reynolds number both in steady and unsteady environment. Lastly, suggestion for future work has been highlighted.

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.

Made available in DSpace on 2014-12-17T14:58:20Z (GMT). No. of bitstreams: 1 BrunoCD_DISSERT.pdf: 4492853 bytes, checksum: 8fad50750c043cd4dab6a1e27f414d81 (MD5) Previous issue date: 2013-02-15 / One of the current major concerns in engineering is the development of aircrafts that have low power consumption and high performance. So, airfoils that have a high value of Lift Coefficient and a low value for the Drag Coefficient, generating a High-Efficiency airfoil are studied and designed. When the value of the Efficiency increases, the aircraft s fuel consumption decreases, thus improving its performance. Therefore, this work aims to develop a tool for designing of airfoils from desired characteristics, as Lift and Drag coefficients and the maximum Efficiency, using an algorithm based on an Artificial Neural Network (ANN). For this, it was initially collected an aerodynamic characteristics database, with a total of 300 airfoils, from the software XFoil. Then, through the software MATLAB, several network architectures were trained, between modular and hierarchical, using the Back-propagation algorithm and the Momentum rule. For data analysis, was used the technique of cross- validation, evaluating the network that has the lowest value of Root Mean Square (RMS). In this case, the best result was obtained for a hierarchical architecture with two modules and one layer of hidden neurons. The airfoils developed for that network, in the regions of lower RMS, were compared with the same airfoils imported into the software XFoil / Uma das maiores preocupa??es atuais na Engenharia ? o desenvolvimento de aeronaves que possuam baixo consumo e alto desempenho. Para isso, s?o estudados e projetados perfis aerodin?micos que tenham um valor elevado de coeficiente de sustenta??o e um valor baixo para o coeficiente de arrasto, gerando um perfil de alta efici?ncia. Quanto maior o valor da efici?ncia, menor ser? o consumo de combust?vel da aeronave, melhorando assim, o seu desempenho. Neste sentido, este trabalho objetiva desenvolver uma ferramenta para cria??o de perfis aerodin?micos a partir de caracter?sticas desejadas, como coeficiente de sustenta??o e de arrasto e efici?ncia m?xima, utilizando-se um algoritmo baseado em uma Rede Neural Artificial (RNA). Para isso, inicialmente foram coletados uma base de dados de caracter?sticas aerodin?micas de um total de 300 perfis, a partir do software XFoil. Ent?o, atrav?s de uma rotina implementada no software MATLAB, foram treinadas diversas arquiteturas de redes, entre modulares e hier?rquicas, utilizando-se o algoritmo de Retropropaga??o e a regra do Momento. Para an?lise dos resultados, foi utilizada a t?cnica de valida??o cruzada, avaliando a rede que possuiu o menor valor de Erro M?dio Quadr?tico (EMQ). Neste caso, o melhor resultado obtido foi para uma arquitetura hier?rquica com dois m?dulos e uma camada de neur?nios ocultos. Os perfis aerodin?micos desenvolvidos por essa rede, nas regi?es de menor EMQ, foram comparados aos mesmos perfis importados ao software XFoil. O presente trabalho oferece como contribui??o, em rela??o a outros trabalhos que envolvem RNA aplicada ? mec?nica dos fluidos, o desenvolvimento de perfis aerodin?micos a partir de suas caracter?sticas aerodin?micas

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.

This report concerns the research, design, implementation and application of an airfoil trailing-edge noise prediction tool in the development of new, quieter airfoil for large-size wind turbine application. The tool is aimed at enabling comparative acoustic performance assessment of airfoils during the early development cycle of new blades and rotors for wind turbine applications. The ultimate goal is to enable the development of quieter wind turbines by the Wind Energy Industry. The task was accomplished by developing software that is simultaneously suitable for comparative design, computationally efficient and user-friendly. The tool was integrated into a state-of-the-art wind turbine design and analysis code that may be downloaded from the web, in compiled or source code form, under general public licensing, at no charge. During the development, an extensive review of the existing airfoil trailing-edge noise prediction models was accomplished, and the semi-empirical BPM model was selected and modified to cope with generic airfoil geometry. The intrinsic accuracy of the original noise prediction model was evaluated as well as its sensitivity to the turbulence length scale parameter, with restrictions imposed accordingly. The criterion allowed comparison of performance of both CFD-RANS and a hybrid solver (XFLR5) on the calculation of the turbulent boundary layer data, with the eventual adjustment and selection of the latter. After all the elements for assembling the method had been selected and the code specified, a collaboration project was made effective between Poli-USP and TU-Berlin, which allowed the seamless coupling of the new airfoil TE noise module, \"PNoise\", to the popular wind turbine design/analysis integrated environment, \"QBlade\". After implementation, the code calculation routines were thoroughly verified and then used in the development of a family of \"silent profiles\" with good relative acoustic and aerodynamic performance. The sample airfoil development study closed the initial design cycle of the new tool and illustrated its ability to fulfill the originally intended purpose of enabling the design of new, quieter blades and rotors for the advancement of the Wind Energy Industry with limited environmental footprint. / Este trabalho descreve a pesquisa de elementos iniciais, o projeto, a implantação e a aplicação de uma ferramenta de predição de ruído de bordo de fuga, no desenvolvimento de aerofólios mais silenciosos para turbinas eólicas de grande porte. O objetivo imediato da ferramenta é permitir a comparação de desempenho acústico relativo entre aerofólios no início do ciclo de projeto de novas pás e rotores de turbinas eólicas. O objetivo mais amplo é possibilitar o projeto de turbinas eólicas mais silenciosas, mas de desempenho aerodinâmico preservado, pela indústria da Energia Eólica. A consecução desses objetivos demandou o desenvolvimento de uma ferramenta que reunisse, simultaneamente, resolução comparativa, eficiência computacional e interface amigável, devido à natureza iterativa do projeto preliminar de um novo rotor. A ferramenta foi integrada a um ambiente avançado de projeto e análise de turbinas eólicas, de código aberto, que pode ser livremente baixado na Web. Durante a pesquisa foi realizada uma ampla revisão dos modelos existentes para predição de ruído de bordo de fuga, com a seleção do modelo semi-empírico BPM, que foi modificado para lidar com geometrias genéricas. A precisão intrínseca do modelo original foi avaliada, assim como sua sensibilidade ao parâmetro de escala de turbulência transversal, com restrições sendo impostas a esse parâmetro em decorrência da análise. Esse critério permitiu a comparação de resultados de cálculo provenientes de método CFD-RANS e de método híbrido (XFLR5) de solução da camada limite turbulenta, com a escolha do último. Após a seleção de todos os elementos do método e especificação do código, uma parceria foi estabelecida entre a Poli-USP e a TU-Berlin, que permitiu a adição de um novo módulo de ruído de bordo de fuga, denominado \"PNoise\", ao ambiente de projeto e análise integrado de turbinas eólicas \"QBlade\". Após a adição, as rotinas de cálculo foram criteriosamente verificadas e, em seguida, aplicadas ao desenvolvimento de aerofólios mais silenciosos, com bons resultados acústicos e aerodinâmicos relativos a uma geometria de referência. Esse desenvolvimento ilustrou a capacidade da ferramenta de cumprir a missão para a qual foi inicialmente projetada, qual seja, permitir à Indústria desenvolver pás mais silenciosas que irão colaborar com o avanço da energia eólica através da limitação do seu impacto ambiental.

This report concerns the research, design, implementation and application of an airfoil trailing-edge noise prediction tool in the development of new, quieter airfoil for large-size wind turbine application. The tool is aimed at enabling comparative acoustic performance assessment of airfoils during the early development cycle of new blades and rotors for wind turbine applications. The ultimate goal is to enable the development of quieter wind turbines by the Wind Energy Industry. The task was accomplished by developing software that is simultaneously suitable for comparative design, computationally efficient and user-friendly. The tool was integrated into a state-of-the-art wind turbine design and analysis code that may be downloaded from the web, in compiled or source code form, under general public licensing, at no charge. During the development, an extensive review of the existing airfoil trailing-edge noise prediction models was accomplished, and the semi-empirical BPM model was selected and modified to cope with generic airfoil geometry. The intrinsic accuracy of the original noise prediction model was evaluated as well as its sensitivity to the turbulence length scale parameter, with restrictions imposed accordingly. The criterion allowed comparison of performance of both CFD-RANS and a hybrid solver (XFLR5) on the calculation of the turbulent boundary layer data, with the eventual adjustment and selection of the latter. After all the elements for assembling the method had been selected and the code specified, a collaboration project was made effective between Poli-USP and TU-Berlin, which allowed the seamless coupling of the new airfoil TE noise module, \"PNoise\", to the popular wind turbine design/analysis integrated environment, \"QBlade\". After implementation, the code calculation routines were thoroughly verified and then used in the development of a family of \"silent profiles\" with good relative acoustic and aerodynamic performance. The sample airfoil development study closed the initial design cycle of the new tool and illustrated its ability to fulfill the originally intended purpose of enabling the design of new, quieter blades and rotors for the advancement of the Wind Energy Industry with limited environmental footprint. / Este trabalho descreve a pesquisa de elementos iniciais, o projeto, a implantação e a aplicação de uma ferramenta de predição de ruído de bordo de fuga, no desenvolvimento de aerofólios mais silenciosos para turbinas eólicas de grande porte. O objetivo imediato da ferramenta é permitir a comparação de desempenho acústico relativo entre aerofólios no início do ciclo de projeto de novas pás e rotores de turbinas eólicas. O objetivo mais amplo é possibilitar o projeto de turbinas eólicas mais silenciosas, mas de desempenho aerodinâmico preservado, pela indústria da Energia Eólica. A consecução desses objetivos demandou o desenvolvimento de uma ferramenta que reunisse, simultaneamente, resolução comparativa, eficiência computacional e interface amigável, devido à natureza iterativa do projeto preliminar de um novo rotor. A ferramenta foi integrada a um ambiente avançado de projeto e análise de turbinas eólicas, de código aberto, que pode ser livremente baixado na Web. Durante a pesquisa foi realizada uma ampla revisão dos modelos existentes para predição de ruído de bordo de fuga, com a seleção do modelo semi-empírico BPM, que foi modificado para lidar com geometrias genéricas. A precisão intrínseca do modelo original foi avaliada, assim como sua sensibilidade ao parâmetro de escala de turbulência transversal, com restrições sendo impostas a esse parâmetro em decorrência da análise. Esse critério permitiu a comparação de resultados de cálculo provenientes de método CFD-RANS e de método híbrido (XFLR5) de solução da camada limite turbulenta, com a escolha do último. Após a seleção de todos os elementos do método e especificação do código, uma parceria foi estabelecida entre a Poli-USP e a TU-Berlin, que permitiu a adição de um novo módulo de ruído de bordo de fuga, denominado \"PNoise\", ao ambiente de projeto e análise integrado de turbinas eólicas \"QBlade\". Após a adição, as rotinas de cálculo foram criteriosamente verificadas e, em seguida, aplicadas ao desenvolvimento de aerofólios mais silenciosos, com bons resultados acústicos e aerodinâmicos relativos a uma geometria de referência. Esse desenvolvimento ilustrou a capacidade da ferramenta de cumprir a missão para a qual foi inicialmente projetada, qual seja, permitir à Indústria desenvolver pás mais silenciosas que irão colaborar com o avanço da energia eólica através da limitação do seu impacto ambiental.

Analysis of unsteady flows is a very challenging topic of research. A decade ago, potential flow equations were used to predict unsteady pressures on oscillating bodies. Recognising the fact that nonlinear aerodynamics is essential to analyse unsteady flows accurately, particularly in transonic and supersonic flows, different Euler formulations operating on moving grids have emerged recently as important CFD tools for unsteady aerodynamics. Numerical solution of Euler equations on moving grids based on upwind schemes such as the ones due to van Leer and Roe have been developed for the purpose of numerical simulation of unsteady transonic and supersonic flows. In the present work, Euler computations based on yet another recent robust upwind scheme (for steady flows) namely Kinetic Flux Vector Splitting (KFVS) scheme due to Deshpande and Mandal is chosen for further development of a time accurate Euler solver to operate on problems involving moving boundaries. The development of an Euler code based on this scheme is likely to be highly useful to analyse problems of unsteady aerodynamics and computational aeroelasiticity especially when it is noted that KFVS has been found to be an extremely robust scheme for computation of subsonic, transonic, supersonic and hypersonic flows. The KFVS scheme, basically exploits the connection between the linear scalar Boltzmann equation of kinetic theory of gases and the nonlinear vector conservation law, that is, Euler equations of fluid dynamics through moment method strategy. The KFVS scheme has inherent simplicity in splitting the flux even on moving grids due to underlying particle model. The inherent simplicity of KFVS for moving grid problems is due to its relationship with the Boltzmann equation. If a surface is moving with velocity w and a particle has velocity v, then it is quite reasonable to do the splitting based on (v-w)<0 or >0. Only particles having velocity v greater than w will cross the moving surface from left to right and similar arguments hold good for particles moving in opposite direction. It is therefore quite natural to extend KFVS by splitting the Maxwellian velocity distribution at Boltzmann level based on the sign of the normal component of the relative velocity. The relative velocity is the difference between the molecular velocity (v) and the velocity of the moving surface(w). This inherent simplicity of the Kinetic Flux Vector Splitting scheme on Moving Grids (KFMG) method has prompted us to extend the same ideas to 2-D and 3-D problems leading to the present KFMG method. If w is set to zero then KFMG formulation reduces to the one corresponding to KFVS. Thus KFMG formulations axe generalisation of the KFVS formulation. In 2-D and 3-D cases, in addition to the KFMG formulation, the method to move the grids, the appropriate boundary conditions for treating moving surfaces and techniques to improve accuracy in space and time are required to be developed. The 2-D and 3-D formulations based on Kinetic Flux Vector Splitting scheme on Moving Grids method have been developed for computing unsteady flows. Between two successive time steps, the body changes its orientation in case of an oscillation or it deforms when subjected to, aerodynamic loads. In either of these cases the grid corresponding to the first time step has to be moved or regenerated around the displaced or deformed body. There are several approaches available to generate grids around moving bodies. In the present work, the 'spring analogy method' is followed to obtain grid around deflected geometries within the frame work of structured grid. Using this method, the grids are moved from previous time to the current time. This method is capable of tackling any kind of aeroelastic deformation of the body. For oscillating bodies, a suitable boundary condition enforcing the flow tangency on the body needs to be developed. As a first attempt, the body surface has been treated as an 1-D piston undergoing compression and expansion. Then, a more general Kinetic Moving Boundary Condition(KMBC) has been developed. The KMBC uses specular reflection model of kinetic theory of gases. In order to treat fixed outer boundary, Kinetic Outer Boundary Condition(KOBC) has been applied. The KOBC is more general in the sense that, it can treat different type of boundaries (subsonic, supersonic, inflow or out flow boundary). A 2-D cell-centered finite volume KFMG Euler code to operate on structured grid has been developed. The time accuracy is achieved by incorporating a fourth order Runge-Kutta time marching method. The space accuracy has been enhanced by using high resolution scheme as well as second order scheme using the method of reconstruction of fluxes. First, the KFMG Euler code has been applied to standard test cases for computing steady flows around NACA 0012 and NACA 64AQ06 airfoils in transonic flow. For these two airfoils both computational and experimental results are available in literature. It is thus possible to verify (that is, prove the claim that code is indeed solving the partial differential equations + boundary conditions posed to the code) and validate(that is, comparison with experimental results) the 2-D KFMG Euler code. Having verified and validated the 2-D KFMG Euler code for the standard test cases, the code is then applied to predict unsteady flows around sinusoidally oscillating NACA 0012 and NACA 64A006 airfoils in transonic flow. The computational and experimental unsteady results are available in literature for these airfoils for verification and validation of the present results. The unsteady lift and normal force coefficients have been predicted fairly accurately by all the CFD codes. However there is some difficulty about accurate prediction of unsteady pitching moment coefficient. Even Navier-Stokes code could not predict pitching moment accurately. This issue needs further in depth study and probably intensive computation which have not been undertaken in the present study. Next, a two degrees of £reedom(2-DOF) structural dynamics model of an airfoil undergoing pitch and plunge motions has been coupled with the 2-D KFMG Euler code for numerical simulation of aeroelastic problems. This aeroelastic analysis code is applied to NACA 64A006 airfoil undergoing pitch and plunge motions in transonic flow to obtain aeroelastic response characteristics for a set of structural parameters. For this test case also computed results are available in literature for verification. The response characteristics obtained have showed three modes namely stable, neutrally stable and unstable modes of oscillations. It is interesting to compare the value of airfoil-to-air mass ratio (Formula) obtained by us for neutrally stable condition with similar values obtained by others and some differences between them are worth mentioning here. The values of \i for neutral stability are different for different authors. The differences in values of (Formula) predicted by various authors are primarily due to differences which can be due to grid as well as mathematical model used. For example, the Euler calculations, TSP calculations and full potential calculations always show differences in shock location for the same flow problem. Changes in shock location will cause change in pressure distribution on airfoil which in turn will cause changes in values of \L for conditions of neutral stability. The flutter speed parameter(U*) has also been plotted with free stream Mach number for two different values of airfoil - to - air mass ratio. These curves shown a dip when the free stream Mach number is close to 0.855. This is referred as "Transonic Dip Phenomenon". The shock waves play a dominant role in the mechanism of transonic dip phenomenon. Lastly, cell-centered finite volume KFMG 3-D Euler code has been developed to operate on structured grids. The time accuracy is achieved by incorporating a fourth order Runge-Kutta method. The space accuracy has been enhanced by using high resolution scheme. This code has 3-D grid movement module which is based on spring analogy method. The KMBC to treat oscillating 3-D configuration and KOBC for treating 3-D outer boundary have also been formulated and implemented in the code. The 3-D KFMG Euler code has been first verified and validated for 3-D steady flows around standard shapes such as, transonic flow past a hemisphere cylinder and ONERA M6 wing. This code has also been used for predicting hypersonic flow past blunt cone-eylinder-flare configuration for which experimental data are available. Also, for this case, the results are compared with a similar Euler code. Then the KFMG Euler code has been used for predicting steady flow around ogive-cylinder-ogive configuration with elliptical cross section. The aerodynamic coefficients obtained have been compared with those of another Euler code. Thus, the 3-D KFMG Euler code has been verified and validated extensively for steady flow problems. Finally, the 3-D KFMG based Euler code has been applied to an oscillating ogive-cylinder-ogive configuration in transonic flow. This test case has been chosen as it resembles the core body of a flight vehicle configuration of interest to DRDO,India. For this test case, the unsteady lift coefficients are available in literature for verifying the present results. Two grid sizes are used to perform the unsteady calculations using the present KFMG 3-D Euler code. The hysteresis loops of lift and moment coefficients confirmed the unsteady behaviour during the oscillation of the configuration. This has proved that, the 3-D formulations are capable of predicting the unsteady flows satisfactorily. The unsteady results obtained for a grid with size of 45x41x51 which is very close to the grid size chosen in the reference(Nixon et al.) are considered for comparison. It has been mentioned in the reference that, a phase lag of (Formula) was observed in lift coefficients with respect to motion of the configuration for a free stream Mach number of 0.3 with other conditions remaining the same. The unsteady lift coefficients obtained using KFMG code as well as those available in literature are plotted for the same flow conditions. Approximately the same phase lag of (Formula) is present (for (Formula)) between the lift coefficient curves of KFMG and due to Nixon et al. The phase lag corrected plot of lift coefficient obtained by Nixon et al. is compared with the lift coefficient versus time obtained by 3-D KFMG Euler code. The two results compare well except that the peaks are over predicted by KFMG code. It is nut clear at this stage whether our results should at all match with those due to Nixon et al. Further in depth study is obviously required to settle the issue. Thus the Kinetic Flux Vector Splitting on Moving Grids has been found to be a very good and a sound method for splitting fluxes and is a generalisation of earlier KFVS on fixed grids. It has been found to be very successful in numerical simulation of unsteady aerodynamics and computational aeroelasticity.