An experimental investigation of the unsteady aerodynamics of passing highway vehicles. Final report.

National Highway Traffic Safety Administration, Washington, D.C.

An experimental investigation of the unsteady aerodynamics of passing highway vehicles. Final summary.

National Highway Traffic Safety Administration, Washington, D.C.

Panel Method Formulation for Oscillating Airfoils in Sonic Flow

The mathematical model for two-dimensional unsteady sonic flow, based on the classical diffusion equation with imaginary coefficient, is presented and discussed. The main purpose is to develop a rigorous formulation in order to bring into light the correspondence between the sonic, supersonic and subsonic panel method theory. Source and doublet integrals are obtained and Laplace transformation demonstrates that, in fact, the source integral is the solution of the doublet integral equation. It is shown that the doublet-only formulation reduces to a Volterra integral equation of the first kind and a numerical method is proposed in order to solve it. To the authors' knowledge this is the first reported solution to the unsteady sonic thin airfoil problem through the use of doublet singularities. Comparisons with the source-only formulation are shown for the problem of a flat plate in combined harmonic heaving and pitching motion.

Aeroelastic behavior of a wing with a magnetorheological hysteretic model to the damping

This paper shows the application of a hysteretic model for the Magnetorheological Damper (MRD) placed in the plunge degree-of-freedom of aeroelastic model of a wing. This hysteretic MRD model was developed by the researchers of the French Aerospace Lab. (ONERA) and describe, with a very good precision, the hysteretic behavior of the MRD. The aeroelastic model used in this paper do not have structural nonlinearities, the only nonlinearities showed in the model, are in the unsteady flow equations and are the same proposed by Theodorsen and Wagner in their unsteady aerodynamics theory; and the nonlinearity introduced by the hysteretic model used. The main objective of this paper is show the mathematical modeling of the problem and the equations that describes the aeroelastic response of our problem; and the gain obtained with the introduction of this hysteretic model in the equations with respect to other models that do not show the this behavior, through of pictures that represents the time response and Phase diagrams. These pictures are obtained using flow velocities before and after the flutter velocity. Finally, an open-loop control was made to show the effect of the MRD in the aeroelastic behavior.

Asymptotic Description of Flutter Amplitude Saturation by Nonlinear Friction Forces

The computation of the non-linear vibration dynamics of an aerodynamically unstable bladed-disk is a formidable numerical task, even for the simplified case of aerodynamic forces assumed to be linear. The nonlinear friction forces effectively couple dif- ferent travelling waves modes and, in order to properly elucidate the dynamics of the system, large time simulations are typically required to reach a final, saturated state. Despite of all the above complications, the output of the system (in the friction microslip regime) is basically a superposition of the linear aeroelastic un- stable travelling waves, which exhibit a slow time modulation that is much longer than the elastic oscillation period. This slow time modulation is due to both, the small aerodynamic effects and the small nonlinear friction forces, and it is crucial to deter- mine the final amplitude of the flutter vibration. In this presenta- tion we apply asymptotic techniques to obtain a new simplified model that captures the slow time dynamics of the amplitudes of the travelling waves. The resulting asymptotic model is very re- duced and extremely cheap to simulate, and it has the advantage that it gives precise information about the characteristics of the nonlinear friction models that actually play a role in the satura- tion of the vibration amplitude.

Asymptotic Description of Flutter Amplitude Saturation by Nonlinear Friction Forces

The computation of the non-linear vibration dynamics of an aerodynamically unstable bladed-disk is a formidable numerical task, even for the simplified case of aerodynamic forces assumed to be linear. The nonlinear friction forces effectively couple dif- ferent travelling waves modes and, in order to properly elucidate the dynamics of the system, large time simulations are typically required to reach a final, saturated state. Despite of all the above complications, the output of the system (in the friction microslip regime) is basically a superposition of the linear aeroelastic un- stable travelling waves, which exhibit a slow time modulation that is much longer than the elastic oscillation period. This slow time modulation is due to both, the small aerodynamic effects and the small nonlinear friction forces, and it is crucial to deter- mine the final amplitude of the flutter vibration. In this presenta- tion we apply asymptotic techniques to obtain a new simplified model that captures the slow time dynamics of the amplitudes of the travelling waves. The resulting asymptotic model is very re- duced and extremely cheap to simulate, and it has the advantage that it gives precise information about the characteristics of the nonlinear friction models that actually play a role in the satura- tion of the vibration amplitude.

Flight loads analysis of a maneuvering transport aircraft

The paper provides a method applicable for the determination of flight loads for maneuvering aircraft, in which aerodynamic loads are calculated based on doublet lattice method, which contains three primary steps. Firstly, non-dimensional stability and control derivative coefficients are obtained through solving unsteady aerodynamics in subsonic flow based on a doublet lattice technical. These stability and control derivative coefficients are used in second step. Secondly, the simulation of aircraft dynamic maneuvers is completed utilizing fourth order Runge-Kutta method to solve motion equations in different maneuvers to gain response parameters of aircraft due to the motion of control surfaces. Finally, the response results calculated in the second step are introduced to the calculation of aerodynamic loads. Thus, total loads and loads distribution on different components of aircraft are obtained. According to the above method, abrupt pitching maneuvers, rolling maneuvers and yawing maneuvers are investigated respectively.

Physics of Vibrating Airfoils at Low Reduced Frequency

The unsteady aerodynamics of low pressure turbine vibrating airfoils in flap mode is studied in detail using a frequency domain linearized Navier-Stokes solver. Both the travelling-wave and influence coefficient formulations of the problem are used to highlight key aspects of the physics and understand different trends such as the effect of reduced frequency and Mach number. The study is focused in the low-reduced frequency regime which is of paramount relevance for the design of aeronautical low-pressure turbines and compressors. It is concluded that the effect of the Mach number on the unsteady pressure phase can be neglected in first approximation and that the unsteadiness of the vibrating and adjacent airfoils is driven by vortex shedding mechanisms. Finally a simple model to estimate the work-per-cycle as a function of the reduced frequency and Mach Number is provided. The edge-wise and torsion modes are presented in less detail but it is shown that acoustic waves are essential to explain its behaviour. The non-dimensional work-per-cycle of the edge-wise mode shows a large dependence with the Mach number while in the torsion mode a large number of airfoils is needed to reconstruct the work-per-cycle departing from the influence coefficients.

Design of an Experimental Rig to Characterize the Gust Response of Small Wind Turbines and Autorotating Seeds

This thesis investigates the rotational behavior of abstracted small-wind-turbine rotors exposed to a sudden increase in oncoming flow velocity, i.e. a gust. These rotors consisted of blades with aspect ratios characteristic of samara seeds, which are known for their ability to maintain autorotation in unsteady wind. The models were tested in a towing tank using a custom-built experimental rig. The setup was designed and constructed to allow for the measurement of instantaneous angular velocity of a rotor model towed at a prescribed kinematic profile along the tank. The conclusions presented in this thesis are based on the observed trends in effective angle-of-attack distribution, tip speed ratio, angular velocity, and time delay in the rotational response for each of rotors over prescribed gust cases. It was found that the blades with the higher aspect ratio had higher tip speed ratios and responded faster than the blades with a lower aspect ratio. The decrease in instantaneous tip speed ratio during the onset of a prescribed gust correlated with the time delay in each rotor model's rotational response. The time delays were found to increase nonlinearly with decreasing durations over which the simulated gusts occurred.

Non-Linear Unsteady Aerodynamic Response Approximation Using Multi-Layer Functionals

Non-linear functional representation of the aerodynamic response provides a convenient mathematical model for motion-induced unsteady transonic aerodynamic loads response, that accounts for both complex non-linearities and time-history effects. A recent development, based on functional approximation theory, has established a novel functional form; namely, the multi-layer functional. For a large class of non-linear dynamic systems, such multi-layer functional representations can be realised via finite impulse response (FIR) neural networks. Identification of an appropriate FIR neural network model is facilitated by means of a supervised training process in which a limited sample of system input-output data sets is presented to the temporal neural network. The present work describes a procedure for the systematic identification of parameterised neural network models of motion-induced unsteady transonic aerodynamic loads response. The training process is based on a conventional genetic algorithm to optimise the network architecture, combined with a simplified random search algorithm to update weight and bias values. Application of the scheme to representative transonic aerodynamic loads response data for a bidimensional airfoil executing finite-amplitude motion in transonic flow is used to demonstrate the feasibility of the approach. The approach is shown to furnish a satisfactory generalisation property to different motion histories over a range of Mach numbers in the transonic regime.

FastAero – A Precorrected FFT – Fast Multipole Tree Steady and Unsteady Potential Flow Solver

In this paper a precorrected FFT-Fast Multipole Tree (pFFT-FMT) method for solving the potential flow around arbitrary three dimensional bodies is presented. The method takes advantage of the efficiency of the pFFT and FMT algorithms to facilitate more demanding computations such as automatic wake generation and hands-off steady and unsteady aerodynamic simulations. The velocity potential on the body surfaces and in the domain is determined using a pFFT Boundary Element Method (BEM) approach based on the Green’s Theorem Boundary Integral Equation. The vorticity trailing all lifting surfaces in the domain is represented using a Fast Multipole Tree, time advected, vortex participle method. Some simple steady state flow solutions are performed to demonstrate the basic capabilities of the solver. Although this paper focuses primarily on steady state solutions, it should be noted that this approach is designed to be a robust and efficient unsteady potential flow simulation tool, useful for rapid computational prototyping.

Unsteady flow simulation of a vertical axis augmented wind turbine: a two-dimensional study

As the integration of vertical axis wind turbines in the built environment is a promising alternative to horizontal axis wind turbines, a 2D computational investigation of an augmented wind turbine is proposed and analysed. In the initial CFD analysis, three parameters are carefully investigated: mesh resolution; turbulence model; and time step size. It appears that the mesh resolution and the turbulence model affect result accuracy; while the time step size examined, for the unsteady nature of the flow, has small impact on the numerical results. In the CFD validation of the open rotor with secondary data, the numerical results are in good agreement in terms of shape. It is, however, observed a discrepancy factor of 2 between numerical and experimental data. Successively, the introduction of an omnidirectional stator around the wind turbine increases the power and torque coefficients by around 30–35% when compared to the open case; but attention needs to be given to the orientation of the stator blades for optimum performance. It is found that the power and torque coefficients of the augmented wind turbine are independent of the incident wind speed considered.

Coupling time-stepping numerical methods and standard aerodynamics codes for instability analysis of flows in complex geometries

The development of a global instability analysis code coupling a time-stepping approach, as applied to the solution of BiGlobal and TriGlobal instability analysis 1, 2 and finite-volume-based spatial discretization, as used in standard aerodynamics codes is presented. The key advantage of the time-stepping method over matrix-formulation approaches is that the former provides a solution to the computer-storage issues associated with the latter methodology. To-date both approaches are successfully in use to analyze instability in complex geometries, although their relative advantages have never been quantified. The ultimate goal of the present work is to address this issue in the context of spatial discretization schemes typically used in industry. The time-stepping approach of Chiba 3 has been implemented in conjunction with two direct numerical simulation algorithms, one based on the typically-used in this context high-order method and another based on low-order methods representative of those in common use in industry. The two codes have been validated with solutions of the BiGlobal EVP and it has been showed that small errors in the base flow do not have affect significantly the results. As a result, a three-dimensional compressible unsteady second-order code for global linear stability has been successfully developed based on finite-volume spatial discretization and time-stepping method with the ability to study complex geometries by means of unstructured and hybrid meshes

Application of the pressure sensitive paint technique to steady and unsteady flow /

"June 1996."

The unsteady flow field about a right circular cone in unsteady flight.

"These studies were conducted by the General Electric Company, Reentry Systems Department, for the Stability and Control Section of the Flight Dynamics Laboratory of the Air Force Research and Technology Division."