This work is focused on the aerodynamics for a proprietary laminar flow airfoil for Unmanned Aerial Vehicle (UAV) applications. The two main focuses are (1) aerodynamic performance at Reynolds number on the order of 10,000, (2) the effect of a conventional hot-wire probe on laminar separation bubbles. For aerodynamic performance, pressure and wake velocity distributions were measured at Re = 40,000 and 60,000 for a range of angles of attack. The airfoil performed poorly for these Reynolds numbers due to laminar boundary layer separation. 2-D boundary layer trips significantly improved the lift-to-drag ratio. For probe effects, three Reynolds numbers were investigated (Re = 100,000, 150,000, and 200,000), with three angles of attack for each. Pressure and surface shear distributions were measured. Flow upstream of the probe tip was not affected. Transition was promoted downstream due to the additional disturbances in the separated shear layer.

This work is focused on the aerodynamics for a proprietary laminar flow airfoil for Unmanned Aerial Vehicle (UAV) applications. The two main focuses are (1) aerodynamic performance at Reynolds number on the order of 10,000, (2) the effect of a conventional hot-wire probe on laminar separation bubbles. For aerodynamic performance, pressure and wake velocity distributions were measured at Re = 40,000 and 60,000 for a range of angles of attack. The airfoil performed poorly for these Reynolds numbers due to laminar boundary layer separation. 2-D boundary layer trips significantly improved the lift-to-drag ratio. For probe effects, three Reynolds numbers were investigated (Re = 100,000, 150,000, and 200,000), with three angles of attack for each. Pressure and surface shear distributions were measured. Flow upstream of the probe tip was not affected. Transition was promoted downstream due to the additional disturbances in the separated shear layer.

Micro Air Vehicles (MAV’s) have small size and extreme maneuverability which makes them ideal for surveillance. Propulsion mechanisms include propellers, rotors, and flapping airfoils. Flapping motions, along with biologically-inspired wing profiles, are of interest due to their use of natural physics. Corrugated airfoil structures appears to have poor aerodynamic performance at higher Reynolds numbers, but serve well at Re<10,000. Understanding flow structures around corrugated profiles and comparing them to a standard airfoil will aid in understanding how these corrugated profiles perform well and have been adopted by some of nature’s most acrobatic flyers. Motivation for this investigation is to compare static flow visualizations of corrugated profiles to a standard National Advisory Committee for Aeronautics (NACA) airfoil from low to high angles of attack and further observe flow structure development of a pitching and plunging flat plate at a Re<10,000 and a Strouhal number relevant to natural fliers. The static visualization was conducted at Re=1,000 with a NACA 0012 airfoil and two corrugated models. The Pitch and Plunge Fixture (PPF) developed was constructed by simplifying flapping wings as a two degree of freedom motion in plunge (translation) and pitch (rotation). Results obtained from the PPF were compared with a numerical simulation.

Use of an inflatable/rigidizable wing is explored for Mars airplane designs. The BIG BLUE (Baseline Inflatable-wing Glider Balloon Launched Unmanned airplane Experiment) project was developed at the University of Kentucky, with an objective to demonstrate feasibility of this technology with a flight-test of an high-altitude glider with inflatable/rigidizable wings. The focus of this thesis research was to design and analyze the wing for this project. The wings are stowed in the fuselage, inflate during ascent, and rigidize with exposure to UV light. The design of wings was evaluated by using aerodynamic and finite element software and wind tunnel testing. The profile is chosen based upon aerodynamic results and consideration of manufacturability of the inflatable wing structures. Flow over prototypes of inflatable/rigidizable and ideal shaped wings were also examined in the wind tunnel. Flow visualization, lift and drag measurements, and wake survey testing methods were performed. Results from the wind tunnel testing are presented along with suggestions in improving the inflatable/rigidizable wings aerodynamic efficiency and use on a low Reynolds number platform. In addition, high altitude wing deployment tests and low altitude flight tests of the inflatable/rigidizable wing were conducted.

The later stages of separated shear layer transition within separation bubbles developing over a NACA0018 airfoil operating at a chord Reynolds number of 105 and at angles of attack of 0, 5, 8, and 10 degrees were investigated experimentally in a wind tunnel. Several experimental tools, including a rake of six boundary-layer hot-wire anemometers, were used to perform measurements over the model. Novel high-speed flow visualization performed with a smoke-wire placed within the separated shear layer showed that roll-up vortices are shed within separation bubbles forming on the suction side of the airfoil. The structures were found to convect downstream and eventually break down during laminar-to-turbulent transition. Top view visualizations revealed that, at angles of attack of 0, 5, and 8 degrees, roll-up vortices form coherently across the span and undergo significant spanwise deformations prior to breaking down. At angles of attack of 5 and 8 degrees, rows of streamwise-oriented structures were observed to form during vortex breakdown. Statistics regarding the formation and development of shear layer roll-up vortices were extracted from high-speed flow visualization sequences and compared to the results of boundary layer measurements. It was found that, on the average, roll-up vortices form following the initial exponential growth of unstable disturbances within the separated shear layer and initiate the later stages of transition. The onset of these nonlinear stages was found to occur when the amplitude of velocity disturbances reached approximately 10% of the free-stream velocity. The rate of vortex shedding was found to fall within the frequency band of the unstable disturbances and lie near the central frequency of this band. The formation of vortices has been linked to the generation of harmonics of these unstable disturbances in velocity signals acquired ahead of mean transition. Once shed, vortices were found to drift at speeds between 33% and 44% of the edge velocity. Vortex merging at an angle of attack of 5?? was investigated. It was found that the majority of roll-up vortices proceed to merge with either one or two other vortices. Vortex merging between two and three vortices was found to occur periodically in a process similar to vortex merging in plane mixing layers undergoing subharmonic forcing of the most amplified disturbance. The flapping motion of the separated shear layer was investigated by performing a cross-correlation analysis on the high-speed flow visualization sequences to extract vertical displacement signals of the smoke within the shear layer. The frequency of flapping was found to correspond to the unstable disturbance band. At an angle of attack of 5??, it was found that the separated shear layer has a low-frequency component of flapping that matches a strong peak in velocity and surface pressure spectra that lies outside the unstable disturbance frequency band. The spanwise development of disturbances was assessed in the aft portion of the separation bubbles by performing a cross-correlation analysis on signals acquired simultaneously across the span with the rake of hot-wires. The spanwise correlations between signals was found to be well-correlated ahead of shear layer roll-up, after which disturbances became rapidly uncorrelated ahead of mean reattachment. These results were found to be linked to the coherent roll-up and subsequent breakdown of roll-up vortices.

The accuracy of existing methods for estimating the drag based on experimental flow field measurements were assessed for two-dimensional bodies. The effects of control volume boundary placement and inherent simplifying assumptions were also investigated. Wind tunnel experiments were performed on a circular cylinder operating at a Reynolds number of 8,000 and 20,000, and on a NACA 0018 airfoil operating at a chord Reynolds number of 100,000 for three angles of attack (α), specifically, 5◦, 10◦, and 15◦. The circular cylinder experiments fall within the the shear layer transition flow regime. Airfoil investigations span both types of flow development common to low Reynolds number airfoil operation. For α = 5◦ and 10◦, a separation bubble forms on the upper surface of the airfoil, while, for α = 15◦, the flow separates without reattaching, resulting in a stalled airfoil. Wake velocity and pressure measurements were performed at several downstream locations to investigate the impact of control volume boundary placement. Wake profiles were measured between 3 and 40 diameters downstream from the circular cylinder axis and between 1 and 4.5 chord lengths from the trailing edge of the airfoil. In addition to wake profiles, the outer flow velocity variation was quantified to investigate the appropriate location to measure freestream flow characteristics in a test section with streamwise-varying outer flow conditions. The results show that drag estimates are strongly dependent on the streamwise position of the measured wake profile for all methods investigated. Drag estimates improved, and streamwise variation decreased, with increasing streamwise position of the flow measurements. For the pressure based method examined, wake measurements should be taken at least 10 times the projected model height downstream of the model. In the case of the circular cylinder, this is equivalent to 10 diameters and, for the airfoil investigated, it is approximately 1 chord length from the trailing edge. For the methods relying on velocity measurements, acceptable estimates of drag were possible when based on measurements taken at least 30 projected heights downstream, i.e., 30 diameters for the circular cylinder and 3 chord lengths for the airfoil model investigated. The findings highlight the importance of providing a detailed description of the methodology and experimental implementation for drag estimates based on flow field measurements. Finally the study offers guidelines for implementing momentum integral based drag calculations in future investigations.

Experimental aeroacoustic research was conducted on a wind turbine specific airfoil at low Reynolds numbers. The goal of this thesis was to study trailing edge noise generation from the airfoil and investigate correlations between the noise and the flow field. Before experiments were performed the current wind tunnel had to be modified in order to make it more suitable for aeroacoustic tests. Sound absorbing foam was added to the inside of the tunnel to lower the background noise levels and turbulence reduction screens were added which lowered the turbulence. An S822 airfoil was chosen because it is designed for low Reynolds flows attainable in the wind tunnel which are on the order of 104. Smoke wire flow visualization was used to gain insight into the airfoil wake development and oil film flow visualization was used to qualitatively assess the boundary layer development. Laser Doppler anemometry (LDA) was used to measure two components of velocity at high data rates in the airfoil wake. Wake profiles were measured in addition to single point measurements to determine the velocity spectrum. A microphone was mounted inside the test section in order to measure the trailing edge noise. Initial plans included measuring the trailing edge noise with a microphone array capable of quantifying and locating noise sources. Although an array was built and beamforming code was written it was only used in preliminary monopole source tests. Oil film results showed the behaviour of the boundary layer to be consistent with previous low Reynolds number experiments. LDA results revealed sharp peaks in the velocity spectra at 1100 Hz from U0 = 15–24 m/s, and 3100 and 3800 Hz, from U0 = 25–35 m/s, which were inconsistent with vortex shedding results of previous researchers. Also present were a series of broad peaks in the spectra that increase from 1200–1700 Hz in the U0 = 25–35 m/s range. The shedding frequency from the smoke wire flow visualization was calculated to be 1250 Hz at U0 = 26 m/s. These sharp peaks were also present in the acoustic spectrum. It was reasoned that these peaks are due to wind tunnel resonance which is a common occurrence in hard wall wind tunnels. In particular the tone at 1100 Hz is due to a standing wave with a wavelength equal to half the tunnel width. The shedding frequency from the smoke wire flow visualization was calculated to be 1100 Hz at U0 = 20 m/s. These tones exhibited a “ladder-like” relationship with freestream velocity, another aspect indicative of wind tunnel resonance. It was reasoned that the wind tunnel resonance was forcing the shedding frequency of the airfoil in the U0 = 15–24 m/s range, and in the U0 = 25–35 m/s range, the shedding frequency corresponded to the broad peaks in the LDA spectra.

Designing horizontal-axis wind turbine (HAWT) blades to achieve satisfactory levels of performance starts with knowledge of the aerodynamic forces acting on the blades. In this thesis, HAWT blade design is studied from the aspect of aerodynamic view and the basic principles of the aerodynamic behaviors of HAWTs are investigated. Blade-element momentum theory (BEM) known as also strip theory, which is the current mainstay of aerodynamic design and analysis of HAWT blades, is used for HAWT blade design in this thesis. Firstly, blade design procedure for an optimum rotor according to BEM theory is performed. Then designed blade shape is modified such that modified blade will be lightly loaded regarding the highly loaded of the designed blade and power prediction of modified blade is analyzed. When the designed blade shape is modified, it is seen that the power extracted from the wind is reduced about 10% and the length of modified blade is increased about 5% for the same required power. BLADESIGN which is a user-interface computer program for HAWT blade design is written. It gives blade geometry parameters (chord-length and twist distributions) and design conditions (design tip-speed ratio, design power coefficient and rotor diameter) for the following inputs / power required from a turbine, number of blades, design wind velocity and blade profile type (airfoil type). The program can be used by anyone who may not be intimately concerned with the concepts of blade design procedure and the results taken from the program can be used for further studies.

Title: Solution of inverse problem for a flow around an airfoil Author: Mgr. Jan Šimák Department: Department of Numerical Mathematics Supervisor: prof. RNDr. Miloslav Feistauer, DrSc., dr. h. c., Department of Numerical Mathematics Abstract: The method described in this thesis deals with a solution of an inverse problem for a flow around an airfoil. It can be used to design an airfoil shape according to a specified velocity or pressure distribution along the chord line. The method is based on searching for a fixed point of an operator, which combines an approximate inverse and direct operator. The approximate inverse operator, derived on the basis of the thin airfoil theory, assigns a corresponding shape to the specified distribution. The resulting shape is then constructed using the mean camber line and thickness function. The direct operator determines the pressure or velocity distribution on the airfoil surface. We can apply a fast, simplified model of potential flow solved using the Fredholm integral equation, or a slower but more accurate model of RANS equations with a k-omega turbulence model. The method is intended for a subsonic flow.

Este trabalho consiste na aplicação de métodos de otimização e de dinâmica dos fluidos computacional a turbinas eólicas. O grande crescimento no mercado de energias renováveis exige que turbinas cada vez mais potentes sejam criadas e que o projeto e análise destas seja cada vez mais preciso. A presente dissertação tem como objetivos a otimização um aerofólio para turbinas eólicas, a simulação de um aerofólio de uma turbina eólica com alto ângulo de ataque e a simulação de uma turbina tridimensional. A otimização de aerofólios foi feita com simulações bidimensionais permanentes, utilizando as equações médias de Reynolds e o modelo de turbulência de Spalart-Allmaras, com algoritmos genéticos acoplados a redes neurais artificiais. O cálculo de um aerofólio com alto ângulo de ataque foi feito utilizando simulações de grandes escalas com o modelo dinâmico de Smagorinsky. As simulações de uma turbina tridimensional foram feitas empregando as equações médias de Reynolds em forma permanente, com um termo adicional representando as forças de Coriolis, também com o modelo de turbulência de Spalart-Allmaras. Da primeira etapa pode-se concluir que as simulações bidimensionais permanentes são muito precisas para o aerofólio de referência, com boa concordância nos coeficientes de arrasto, sustentação e pressão. Os algoritmos genéticos geraram bons resultados, com cerca de 8% de aumento da razão sustentação/arrasto e com aproximadamente 50% de economia no tempo computacional ao se utilizar redes neurais artificiais. Na segunda etapa, o cálculo de um aerofólio com alto ângulo de ataque demonstrou necessidade de simulações tridimensionais transientes, pela alta variação dos coeficientes aerodinâmicos ao longo do tempo e alta tridimensionalidade da esteira. Na última etapa, a simulação de uma turbina tridimensional mostrou resultados muito próximos dos experimentais. Muita atenção foi dada na discretização deste caso, chegando a uma malha com 700 mil elementos, enquanto outros autores utilizaram de 3 a 38 milhões de elementos para o mesmo caso. / The present work consists in the application of optimization methods and computational fluid dynamics to wind turbines. The massive growth in renewable energies demands more powerful turbines and more accuracy in their design and analysis. This work has three objectives: optimization of an airfoil for wind turbines, simulation of a wind turbine airfoil in deep stall, and simulation of a three-dimensional wind turbine. The airfoil optimization is accomplished by means of two-dimensional steady-state Reynolds averaged Navier-Stokes simulations with the Spalart-Allmaras turbulence model, with genetic algorithms coupled with artificial neural networks. The airfoil in deep stall is calculated with unsteady three-dimensional Large Eddy Simulations with the dynamic Smagorinsky model. The simulation of a wind turbine is also done by means of the Reynolds averaged Navier-Stokes equations, with an additional term to take the Coriolis forces into account, and the Spalart-Allmaras turbulence model. In the first application, it can be confirmed that the two-dimensional steady state simulations are very accurate for the reference airfoil, with good agreement for drag, lift, and pressure coefficients. Genetic algorithms improved the lift-to-drag ratio about 8%, with a 50% decrease in computational time when using artificial neural networks. For the second application, the airfoil with a high angle of attack showed that transient three-dimensional simulations were indeed required, with a high variation of aerodynamic coefficient as a function of time and the highly three-dimensional wake. In the final part, the three-dimensional wind turbine showed very good agreement with experimental results. A great deal of attention was devoted to the creation of the grid and a mesh with only 700 thousand elements was achieved, while other authors used from 3 to 38 million elements for the same case.