A wind tunnel based investigation into the effects of unsteady yaw angles on the aerodynamics of a simple automotive body has been carried out to increase the understanding of the effects of unsteady onset conditions similar to those experienced in normal driving conditions. Detailed flow field measurements have been made using surface pressure tappings and PIV around a simple automotive model in steady state conditions and these have been compared to measurements made whilst the model was oscillating in the yaw plane. The oscillating motion was created by a motored crank which was used to produce consistent and repeated motion which produced a reduced frequency that indicated that a quasi-static response should be expected. The PIV data are used to compare the wake flow structures and the surface pressures are used to infer aerodynamic loads and investigate the development of the flow structures across the surfaces of the model. This includes a comprehensive comparison of the surface pressures on the sides of the model during a transient and quasi-static yaw angel oscillation. These results show differences between the two test conditions with the oscillating model results containing hysteresis and the greatest differences in the flow field occurring on the leeside of the model. Two configurations of the same model with different rear pillar geometries were used to isolate model specific effects. Square rear pillars create strong and stable trailing vortices which are less affected by the model motion whereas radiused rear pillars created weaker and less steady vortices that mixed with the quasi-2D wake behind the model base and were affected to a greater extent by the model motion. The unsteadiness in the trailing vortex separation feeds upstream into the A-pillar vortex demonstrating that small geometry changes at the rear can affect the entire flow field around the model.

Um modelo numérico para simulação da resposta aeroelástica de asas fixas é proposto. A estratégia adotada no trabalho é a de tratar a aerodinâmica e a dinâmica estrutural separadamente e então acoplá-las na equação de movimento. A caracterização dinâmica de uma asa protótipo é feita pelo método dos elementos finitos e a equação de movimento é escrita em função das coordenadas modais. O carregamento aerodinâmico não-estacionário é determinado pelo método de malha de vórtices. A troca de informações entre as malhas estrutural e aerodinâmica é feita através do método de interpolação por splines de superfície e a equação de movimento é resolvida iterativamente no domínio do tempo, utilizando-se um método preditor-corretor. As teorias de aerodinâmica, dinâmica estrutural e do acoplamento entre elas são apresentadas separadamente, juntamente com os respectivos resultados obtidos. A resposta aeroelástica da asa protótipo é representada por curvas de deslocamentos modais em função do tempo para várias velocidades de vôo e a ocorrência de flutter é verificada quando estas curvas divergem (i.e. as amplitudes aumentam progressivamente). Transformadas de Fourier destas curvas mostram o acoplamento de freqüências característico do fenômeno de flutter. / A numerical model for the simulation of the aeroelastic response of fixed wings is proposed. The methodology used in the work is to treat the aerodynamic and the structural dynamics separately and then couple them in the equation of motion. The dynamic characterization of a prototype wing is done by the finite element method and the equation of motion is written in modal coordinates. The unsteady aerodynamic loads are predicted using the vortex lattice method. The exchange of information between the aerodynamic and structural meshes is done by the surface splines interpolation scheme, and the equation of motion is solved interactively in the time domain, employing a predictor-corrector method. The aerodynamic and structural dynamics theories, and the methodology to couple them, are described separately, together with the corresponding obtained results. The aeroelastic response of the prototype wing is represented by time histories of the modal coordinates for different airspeeds, and the flutter occurrence is verified when the time histories diverge (i.e. the amplitudes keep growing). Fast Fourier Transforms of these time histories show the coupling of frequencies, typical of the flutter phenomenon.

This research presents a detailed study of non-synchronous vibration (NSV) and the development of an efficient design method for NSV. NSV occurs as a result of the complex interaction of an aerodynamic instability with blade vibrations. Two NSV design methods are considered and applied to three test cases: 2-D circular cylinder, 2-D airfoil cascade tip section of a modern compressor, and 3-D high pressure compressor cascade that encountered NSV in rig testing. The current industry analysis method is to search directly for the frequency of the instability using CFD analysis and then compare it with a fundamental blade mode frequency computed from a structural analysis code. The main disadvantage of this method is that the blades' motion is not considered and therefore, the maximum response is assumed to be when the blade natural frequency and fluid frequency are coincident. An alternate approach, the enforced motion method, is also presented. In this case, enforced blade motion is used to promote lock-in of the blade frequency to the fluid natural frequency at a specified critical amplitude for a range of interblade phase angles (IBPAs). For the IBPAs that are locked-on, the unsteady modal forces are determined. This mode is acceptable if the equivalent damping is greater than zero for all IBPAs. A method for blade re-design is also proposed to determine the maximum blade response by finding the limit cycle oscillation (LCO) amplitude. It is assumed that outside of the lock-in region is an off-resonant, low amplitude condition. A significant result of this research is that for all cases studied herein, the maximum blade response is not at the natural fluid frequency as is assumed by the direct frequency search approach. This has significant implications for NSV design analysis because it demonstrates the requirement to include blade motion. Hence, an enforced motion design method is recommended for industry and the current approach is of little value. / Dissertation

This work presents a method for building a three-dimensional mesh from two- dimensional topologically identical layers, for use in aeroelastic simulations. The method allows modeling of large deformations of the wing in both the span direction and deformations in the cord of the wing. In addition, the method allows for the modeling of wings attached to fuselages. The mesh created is a hybrid mesh, which allows cell clustering in the viscous region. The generated mesh is high quality and allows capturing of nonlinear uid structure interactions in the form of limit cycle oscillation.

The aerodynamics of two-dimensional and three-dimensional flapping motion in hover is analyzed in incompressible, laminar flow at low Reynolds number regime. The aim of this study is to understand the physics and the underlying mechanisms of the flapping motion using both numerical tools (Direct Numerical Simulation) and experimental tools (Particle Image Velocimetry PIV technique). Numerical analyses cover both two-dimensional and three-dimensional configurations for different parameters using two different flow solvers. The obtained results are then analyzed in terms of aerodynamic force coefficients and vortex dynamics. Both symmetric and cambered airfoil sections are investigated at different starting angle of attacks. Both numerical and experimental simulations are carried out at Reynolds number 1000. The experimental analysis is carried out using Particle Image Velocimetry (PIV) technique in parallel with the numerical tools. Experimental measurements are taken for both two-dimensional and three-dimensional wing configurations using stereoscopic PIV technique.

ABSTRACT EXPERIMENTAL AND NUMERICAL INVESTIGATION OF FLOW FIELD AROUND FLAPPI G AIRFOILS MAKING FIGURE-OF-EIGHT IN HOVER BASKAN, &Ouml / zge M.Sc., Department of Aerospace Engineering Supervisor: Prof. Dr. H. Nafiz Alemdaroglu September 2009, 94 pages The aim of this study is to investigate the flow field around a flapping airfoil making figure-of-eight motion in hover and to compare these results with those of linear flapping motion. Aerodynamic characteristics of these two-dimensional flapping motions are analyzed in incompressible, laminar flow at very low Reynolds numbers regime using both the numerical (Computational Fluid Dynamics, CFD) and the experimental (Particle Image Velocimetry, PIV) tools. Numerical analyses are performed to investigate the effect of different parameters such as the amplitude of motion in y-direction, angle of attack, Reynolds number and camber on the aerodynamic force coefficients and vortex formation mechanisms. Both symmetric and cambered airfoil sections are investigated at three different starting angles of attack for five different amplitudes of motion in y-direction including linear flapping motion. Experimental simulations are performed in order to verify the numerical results only for linear motion at Reynolds number of 1000 for symmetric and cambered airfoils at three different angles of attack. Computed vortical structures are then compared to vorticity contours obtained from the experiments and advantages of figure-of&ndash / eight motion over linear motion are discussed.

The scope of the thesis is to numerically investigate the aerodynamics of flapping airfoils in hover and forward flight. The flowfields around flapping airfoils are computed by solving the governing equations on moving and/or deforming grids. The effects of Reynolds number, reduced frequency and airfoil geometry on unsteady aerodynamics of flapping airfoils undergoing pure plunge and combined pitch-plunge motions in forward flight are investigated. It is observed that dynamic stall of the airfoil is the main mechanism of lift augmentation for both motions at all Reynolds numbers ranging from 10000 to 60000. However, the strength and duration of the leading edge vortex vary with airfoil geometry and reduced frequency. It is also observed that more favorable force characteristics are achieved at higher reduced frequencies and low plunging amplitudes while keeping the Strouhal number constant. The computed flowfields are compared with the wide range of experimental studies and high fidelity simulations thus it is concluded that the present approach is applicable for investigating the flapping wing aerodynamics in forward flight. The effects of vertical translation amplitude and Reynolds number on flapping airfoils in hover are also studied. As the vertical translation amplitude increases, the vortices become stronger and the formation of leading edge vortex is pushed towards the midstroke of the motion. The instantaneous aerodynamic forces for a given figure-of-eight motion do not alter significantly for Reynolds numbers ranging from 500 to 5500.

The aerodynamics phenomena of flapping motion in hover are considered in view of the future Micro Air Vehicle applications. The aim of this work is to characterize the vortex dynamics generated by the wing in motion using direct numerical simulation and experimental analysis then to propose a simplified analytical model for prediction of the forces in order to optimize the parameters of the motion leading to maximum force. A great number of cases are investigated corresponding to different angles of attack, location of start of change of incidence, location of start of change of velocity, axis of rotation, and Re number. The airfoil used is symmetrical. The flow is assumed to be incompressible and laminar with the Reynolds numbers between 500 and 2000. The experimental results obtained by the laser sheet visualization and the Particle Image Velocimetry (PIV) techniques are used in parallel with the direct numerical simulation results for the phenomenological analysis of the flow. The model developed for the aerodynamic forces is an indicial method based on the use of the Duhamel Integral and the results obtained by this model are compared with the ones of the numerical simulations.

The present thesis summarizes work on developing a method for unsteady aerodynamic analysis primarily for aeroelastic simulations. In contrast to widely used prediction tools based on frequency-domain representations, the current approach aims to provide a time-domain simulation capability which can be readily integrated with possibly nonlinear structural and control system models. Further, due to the potential flow model underlying the computational method, and the solution algorithm based on an efficient boundary element formulation, the computational effort for the solution is moderate, allowing time-dependent simulations of complex configurations. The computational method is applied to simulate a number of wind-tunnel experiments involving highly flexible models. Two of the experiments are utilized to verify the method and to ascertain the validity of the unsteady flow model. In the third study, simulations are used for the numerical optimization of a configuration with multiple control surfaces. Here, the flexibility of the model is exploited in order to achieve a reduction of induced drag. Comparison with experimental results shows that the numerical method attains adequate accuracy within the inherent limits of the potential flow model. Finally, rather extensive aeroelastic simulations are performed for the ASK 21 sailplane. Time-domain simulations of a pull-up maneuver and comparisons with flight test data demonstrate that, considering modeling and computational effort, excellent agreement is obtained. Furthermore, a flutter analysis is performed for the same aircraft using identified frequency-domain loads. Results are found to deviate only slightly from critical speed and frequency obtained using an industry-standard aeroelastic analysis code. Nevertheless, erratic results for control surface hinge moments indicate that the accuracy of the present method would benefit from improved control surface modeling and coupled boundary layer analysis. / QC 20100531

Um modelo numérico para simulação da resposta aeroelástica de asas fixas é proposto. A estratégia adotada no trabalho é a de tratar a aerodinâmica e a dinâmica estrutural separadamente e então acoplá-las na equação de movimento. A caracterização dinâmica de uma asa protótipo é feita pelo método dos elementos finitos e a equação de movimento é escrita em função das coordenadas modais. O carregamento aerodinâmico não-estacionário é determinado pelo método de malha de vórtices. A troca de informações entre as malhas estrutural e aerodinâmica é feita através do método de interpolação por splines de superfície e a equação de movimento é resolvida iterativamente no domínio do tempo, utilizando-se um método preditor-corretor. As teorias de aerodinâmica, dinâmica estrutural e do acoplamento entre elas são apresentadas separadamente, juntamente com os respectivos resultados obtidos. A resposta aeroelástica da asa protótipo é representada por curvas de deslocamentos modais em função do tempo para várias velocidades de vôo e a ocorrência de flutter é verificada quando estas curvas divergem (i.e. as amplitudes aumentam progressivamente). Transformadas de Fourier destas curvas mostram o acoplamento de freqüências característico do fenômeno de flutter. / A numerical model for the simulation of the aeroelastic response of fixed wings is proposed. The methodology used in the work is to treat the aerodynamic and the structural dynamics separately and then couple them in the equation of motion. The dynamic characterization of a prototype wing is done by the finite element method and the equation of motion is written in modal coordinates. The unsteady aerodynamic loads are predicted using the vortex lattice method. The exchange of information between the aerodynamic and structural meshes is done by the surface splines interpolation scheme, and the equation of motion is solved interactively in the time domain, employing a predictor-corrector method. The aerodynamic and structural dynamics theories, and the methodology to couple them, are described separately, together with the corresponding obtained results. The aeroelastic response of the prototype wing is represented by time histories of the modal coordinates for different airspeeds, and the flutter occurrence is verified when the time histories diverge (i.e. the amplitudes keep growing). Fast Fourier Transforms of these time histories show the coupling of frequencies, typical of the flutter phenomenon.