Estudio y validación experimental del efecto de aire interno en un perfil aerodinámico hueco

La presente investigación se llevó a cabo en la Universidad Politécnica de Madrid (España) conjuntamente con la Universidad Nacional Experimental del Táchira (Venezuela). El estudio consistió en diseñar una cavidad interna dentro del perfil aerodinámico 2415-3s, el cual fue desarrollado en la Universidad Técnica Checa (Praga, República Checa). Se realizó un estudio computacional, mediante la técnica del CFD, de diferentes modelos de cavidades internas en este perfil, para seleccionar el diseño más adecuado, fabricando un prototipo en 3D; logrando de esta manera validar la simulación computacional con los datos experimentales obtenidos con los ensayos en el túnel de viento AF6109 de la Universidad Nacional Experimental del Táchira. También se aplicaron técnicas de visualización en el túnel de viento, como líneas de corriente de humo y películas de aceite sobre el perfil aerodinámico. Dicho procedimiento permitió corroborar la validación de la simulación computacional. El perfil aerodinámico seleccionado se denominó 2415-3s-TC, cuya característica principal consiste en tres canales independientes entre sí, alojados dentro de la cavidad interna, permitiendo que el flujo de aire forzado a través de la cavidad, cambiara de dirección, para desembocar lo más tangencialmente, así como, lo más perpendicularmente posible al escalón del perfil aerodinámico 2415-3s. Esta configuración de diseñó permitió elevar el coeficiente de sustentación para ángulos de ataque mayores a 8º, así como para ángulos cercanos al ángulo crítico. ABSTRACT This research was conducted at the Polytechnic University of Madrid (Spain) together with the National Experimental University of Táchira (Venezuela). The study was to design an internal cavity within the airfoil 2415-3s, which was developed in the Czech Technical University (Prague, Czech Republic). A computational study was performed using CFD technique, different models of internal cavities in the profile to select the most appropriate design, manufacturing a prototype 3D; thus achieving validate the computer simulation with experimental data obtained from the tests in the wind tunnel AF6109 of the National Experimental University of Táchira. Visualization techniques were also applied in the wind tunnel, as streamlines smoke and oil films on the airfoil. This procedure corroborated validation of computational simulation. The airfoil selected denominated 2415-3s-TC, whose main characteristic consists of three independent channels each other, housed within the inner cavity, allowing the forced air flow through the cavity, change direction, to lead as more tangentially and, as perpendicular as possible to the step 2415-3s aerofoil. This configuration designed allowed increasing the lift coefficient for higher angles of attack to 8º, and for angles near the critical angle.

On the circulation and the position of the forward stagnation point on airfoils

Lift and velocity circulation around airfoils are two aspects of the same phenomenon when airfoils are not stalled and the Kutta—Joukowski theorem applies. This theorem establishes a linear dependence between lift and circulation, which breaks when stalling occurs. As the angle of attack increases beyond this point, the circulation vanishes. Since the circulation determines to a great extent the position of the forward stagnation point on an airfoil, the measurement of this position is an easy and simple way to determine the circulation, which is of help in understanding the role of the latter in the generation of aerodynamic forces on airfoils.

Three-component force measurements on a large scramjet in a shock tunnel

A stress-wave force balance for measurement of thrust, lift, and pitching moment on a large scramjet model (40 kg in mass, 1.165 in in length) in a reflected shock tunnel has been designed, calibrated, and tested. Transient finite element analysis was used to model the performance of the balance. This modeling indicates that good decoupling of signals and low sensitivity of the balance to the distribution of. the load can be achieved with a three-bar balance. The balance was constructed and calibrated by applying a series of point loads to the model. A good comparison between finite element analysis and experimental results was obtained with finite element analysis aiding in the interpretation of some experimental results. Force measurements were made in a shock tunnel both with and without fuel injection, and measurements were compared with predictions using simple models of the scramjet and combustion. Results indicate that the balance is capable of resolving lift, thrust, and pitching moments with and without combustion. However vibrations associated with tunnel operation interfered with the signals indicating the importance of vibration isolation for accurate measurements.

Measurement of the aerodynamic performance of a scramjet-powered vehicle in a shock tunnel

A force balance system for measuring lift, thrust and pitching moment has been used to measure the performance of fueled scramjet-powered vehicle in the T4 Shock Tunnel at The University of Queensland. Detailed measurements have been made of the effects of different fuel flow rates corresponding to equivalence ratios between 0.0 and 1.5. For proposed scramjet-powered vehicles, the fore-body of the vehicle acts as part of the inlet to the engine and the aft-body acts as the thrust surface for the engine. This type of engine-integrated design leads to a strong coupling between the performance of the engine and the lift and trim characteristics of the vehicle. The measurements show that the lift force increased by approximately 50% and centre-of-pressure changed by approximately 10% of the chord of the vehicle when the equivalence ratio varied from 0.0 to 1.0. The results demonstrate the importance of engine performance to the overall aerodynamic characteristics of engine-integrated scramjet vehicles and that such characteristics can be measured in a shock tunnel.

Physical Insights, Steady Aerodynamic Effects, and a Design Tool for Low-Pressure Turbine Flutter

The successful, efficient, and safe turbine design requires a thorough understanding of the underlying physical phenomena. This research investigates the physical understanding and parameters highly correlated to flutter, an aeroelastic instability prevalent among low pressure turbine (LPT) blades in both aircraft engines and power turbines. The modern way of determining whether a certain cascade of LPT blades is susceptible to flutter is through time-expensive computational fluid dynamics (CFD) codes. These codes converge to solution satisfying the Eulerian conservation equations subject to the boundary conditions of a nodal domain consisting fluid and solid wall particles. Most detailed CFD codes are accompanied by cryptic turbulence models, meticulous grid constructions, and elegant boundary condition enforcements all with one goal in mind: determine the sign (and therefore stability) of the aerodynamic damping. The main question being asked by the aeroelastician, ``is it positive or negative?'' This type of thought-process eventually gives rise to a black-box effect, leaving physical understanding behind. Therefore, the first part of this research aims to understand and reveal the physics behind LPT flutter in addition to several related topics including acoustic resonance effects. A percentage of this initial numerical investigation is completed using an influence coefficient approach to study the variation the work-per-cycle contributions of neighboring cascade blades to a reference airfoil. The second part of this research introduces new discoveries regarding the relationship between steady aerodynamic loading and negative aerodynamic damping. Using validated CFD codes as computational wind tunnels, a multitude of low-pressure turbine flutter parameters, such as reduced frequency, mode shape, and interblade phase angle, will be scrutinized across various airfoil geometries and steady operating conditions to reach new design guidelines regarding the influence of steady aerodynamic loading and LPT flutter. Many pressing topics influencing LPT flutter including shocks, their nonlinearity, and three-dimensionality are also addressed along the way. The work is concluded by introducing a useful preliminary design tool that can estimate within seconds the entire aerodynamic damping versus nodal diameter curve for a given three-dimensional cascade.

Model reduction for nonlinear aeroelasticity using the Discrete Empirical Interpolation Method

A novel surrogate model is proposed in lieu of computational fluid dynamic (CFD) code for fast nonlinear aerodynamic modeling. First, a nonlinear function is identified on selected interpolation points defined by discrete empirical interpolation method (DEIM). The flow field is then reconstructed by a least square approximation of flow modes extracted by proper orthogonal decomposition (POD). The proposed model is applied in the prediction of limit cycle oscillation for a plunge/pitch airfoil and a delta wing with linear structural model, results are validate against a time accurate CFD-FEM code. The results show the model is able to replicate the aerodynamic forces and flow fields with sufficient accuracy while requiring a fraction of CFD cost.

Coupled CFD Shape Optimization for aerodynamic profiles

The present document deals with the optimization of shape of aerodynamic profiles -- The objective is to reduce the drag coefficient on a given profile without penalising the lift coefficient -- A set of control points defining the geometry are passed and parameterized as a B-Spline curve -- These points are modified automatically by means of CFD analysis -- A given shape is defined by an user and a valid volumetric CFD domain is constructed from this planar data and a set of user-defined parameters -- The construction process involves the usage of 2D and 3D meshing algorithms that were coupled into own- code -- The volume of air surrounding the airfoil and mesh quality are also parametrically defined -- Some standard NACA profiles were used by obtaining first its control points in order to test the algorithm -- Navier-Stokes equations were solved for turbulent, steady-state ow of compressible uids using the k-epsilon model and SIMPLE algorithm -- In order to obtain data for the optimization process an utility to extract drag and lift data from the CFD simulation was added -- After a simulation is run drag and lift data are passed to the optimization process -- A gradient-based method using the steepest descent was implemented in order to define the magnitude and direction of the displacement of each control point -- The control points and other parameters defined as the design variables are iteratively modified in order to achieve an optimum -- Preliminary results on conceptual examples show a decrease in drag and a change in geometry that obeys to aerodynamic behavior principles

Probabilistic design optimization of horizontal axis wind turbine rotors

Considerable interest in renewable energy has increased in recent years due to the concerns raised over the environmental impact of conventional energy sources and their price volatility. In particular, wind power has enjoyed a dramatic global growth in installed capacity over the past few decades. Nowadays, the advancement of wind turbine industry represents a challenge for several engineering areas, including materials science, computer science, aerodynamics, analytical design and analysis methods, testing and monitoring, and power electronics. In particular, the technological improvement of wind turbines is currently tied to the use of advanced design methodologies, allowing the designers to develop new and more efficient design concepts. Integrating mathematical optimization techniques into the multidisciplinary design of wind turbines constitutes a promising way to enhance the profitability of these devices. In the literature, wind turbine design optimization is typically performed deterministically. Deterministic optimizations do not consider any degree of randomness affecting the inputs of the system under consideration, and result, therefore, in an unique set of outputs. However, given the stochastic nature of the wind and the uncertainties associated, for instance, with wind turbine operating conditions or geometric tolerances, deterministically optimized designs may be inefficient. Therefore, one of the ways to further improve the design of modern wind turbines is to take into account the aforementioned sources of uncertainty in the optimization process, achieving robust configurations with minimal performance sensitivity to factors causing variability. The research work presented in this thesis deals with the development of a novel integrated multidisciplinary design framework for the robust aeroservoelastic design optimization of multi-megawatt horizontal axis wind turbine (HAWT) rotors, accounting for the stochastic variability related to the input variables. The design system is based on a multidisciplinary analysis module integrating several simulations tools needed to characterize the aeroservoelastic behavior of wind turbines, and determine their economical performance by means of the levelized cost of energy (LCOE). The reported design framework is portable and modular in that any of its analysis modules can be replaced with counterparts of user-selected fidelity. The presented technology is applied to the design of a 5-MW HAWT rotor to be used at sites of wind power density class from 3 to 7, where the mean wind speed at 50 m above the ground ranges from 6.4 to 11.9 m/s. Assuming the mean wind speed to vary stochastically in such range, the rotor design is optimized by minimizing the mean and standard deviation of the LCOE. Airfoil shapes, spanwise distributions of blade chord and twist, internal structural layup and rotor speed are optimized concurrently, subject to an extensive set of structural and aeroelastic constraints. The effectiveness of the multidisciplinary and robust design framework is demonstrated by showing that the probabilistically designed turbine achieves more favorable probabilistic performance than those of the initial baseline turbine and a turbine designed deterministically.

Effect of Geometric Uncertainties on the Aerodynamic Characteristic of Offshore Wind Turbine Blades

Offshore wind turbines operate in a complex unsteady flow environment which causes unsteady aerodynamic loads. The unsteady flow environment is characterized by a high degree of uncertainty. In addition, geometry variations and material imperfections also cause uncertainties in the design process. Probabilistic design methods consider these uncertainties in order to reach acceptable reliability and safety levels for offshore wind turbines. Variations of the rotor blade geometry influence the aerodynamic loads which also affect the reliability of other wind turbine components. Therefore, the present paper is dealing with geometric uncertainties of the rotor blades. These can arise from manufacturing tolerances and operational wear of the blades. First, the effect of geometry variations of wind turbine airfoils on the lift and drag coefficients are investigated using a Latin hypercube sampling. Then, the resulting effects on the performance and the blade loads of an offshore wind turbine are analyzed. The variations of the airfoil geometry lead to a significant scatter of the lift and drag coefficients which also affects the damage-equivalent flapwise bending moments. In contrast to that, the effects on the power and the annual energy production are almost negligible with regard to the assumptions made.