On the high-frequency response of unsteady lift and circulation: A dynamical systems perspective

In this paper, we consider the unsteady aerodynamics of a two-dimensional airfoil as a dynamical system whose input is the angle of attack (or airfoil motion) and output is the lift force. Based on this view, we discuss the evolution of lift and circulation from a purely dynamical perspective through step response, frequency response, transfer function, etc. In particular, we point to the relation between the high-frequency gain of the transfer function and the physics of the development of lift and circulation. Based on this view, we show that the circulatory lift dynamics is different from the circulation dynamics. That is, we show that the circulatory lift is not lift due to circulation. In fact, we show that the circulatory–non-circulatory classification is arbitrary. By comparing the steady and unsteady thin airfoil theory, we show that the circulatory lift possesses some acceleration (added-mass) effects. Finally, we perform simulations of Navier–Stokes equations to show that a non-circulatory maneuver in the absence of a free stream induces viscous circulation over the airfoil.     

Asymptotic and spectral properties of operator‐valued functions generated by aircraft wing model

The present paper is devoted to the asymptotic and spectral analysis of an aircraft wing model in a subsonic air flow. The model is governed by a system of two coupled integro‐differential equations and a two parameter family of boundary conditions modelling the action of the self‐straining actuators. The differential parts of the above equations form a coupled linear hyperbolic system; the integral parts are of the convolution type. The system of equations of motion is equivalent to a single operator evolution–convolution equation in the energy space. The Laplace transform of the solution of this equation can be represented in terms of the so‐called generalized resolvent operator, which is an operator‐valued function of the spectral parameter. More precisely, the generalized resolvent is a finite‐meromorphic function on the complex plane having a branch‐cut along the negative real semi‐axis. Its poles are precisely the aeroelastic modes and the residues at these poles are the projectors on the generalized eigenspaces. The dynamics generator of the differential part of the system has been systematically studied in a series of works by the second author. This generator is a non‐selfadjoint operator in the energy space with a purely discrete spectrum. In the aforementioned series of papers, it has been shown that the set of aeroelastic modes is asymptotically close to the spectrum of the dynamics generator, that this spectrum consists of two branches, and a precise spectral asymptotics with respect to the eigenvalue number has been derived. The asymptotical approximations for the mode shapes have also been obtained. It has also been proven that the set of the generalized eigenvectors of the dynamics generator forms a Riesz basis in the energy space. In the present paper, we consider the entire integro‐differential system which governs the model. Namely, we investigate the properties of the integral convolution‐type part of the original system. We show, in particular, that the set of poles of the adjoint generalized resolvent is asymptotically close to the discrete spectrum of the operator that is adjoint to the dynamics generator corresponding to the differential part. The results of this paper will be important for the reconstruction of the solution of the original initial boundary‐value problem from its Laplace transform and for the analysis of the flutter phenomenon in the forthcoming work. Copyright © 2004 John Wiley & Sons, Ltd.     

Unsteady response of airfoils due to small-scale pitching motion with considerations for foil thickness and wake motion

Unsteady pressures, forces, and pitching moments generated by foils experiencing vibratory motion in an incompressible, attached flow configuration are studied within this work. Specifically, two-dimensional, unsteady potential flow and unsteady Reynolds-Averaged Navier–Stokes calculations are performed on various Joukowski foils undergoing sinusoidal, variable amplitude, small-scale pitching motion at a chord-based Reynolds number of 10${}^6$ over a range of reduced frequencies between 0.01–100. These calculated results from both approaches are compared directly to predictions from implementing the Theodorsen model, which treats foils as infinitely thin, flat plates that shed a planar sheet of vorticity. The effects of relaxing these seemingly strict conditions are explored, and the particular terms which control the unsteady responses are identified and discussed. For increasing pitch amplitudes and reduced frequencies the shed wake is seen to become quite non-planar and to form coherent vortex structures. Despite this wake behavior, the normalized airfoil responses at the disturbance reduced frequency are seen to be largely unaffected. However, non-negligible responses are generated across a wide range of other frequencies. Potential flow calculations for symmetric Joukowski foils show that there is marginal effect of foil thickness at reduced frequencies less than one. For higher reduced frequency conditions however, the unsteady lift response is seen to experience both an amplification of level and a phase shift relative to the Theodorsen model. A specific augmenting expression is developed for this behavior through analysis within the potential flow framework.     

Flutter suppression of a plate-like wing via parametric excitation

The present work is motivated by the well known stabilizing effect of parametric excitation of some dynamical systems such as the inverted pendulum. The possibility of suppressing wing flutter via parametric excitation along the plane of highest rigidity in the neighborhood of combination resonance is explored. The nonlinear equations of motion in the presence of incompressible fluid flow are derived using Hamilton's principle and Theodorsen's theory for modeling aerodynamic forces. In the presence of air flow, the bending and torsion modes possess nearly the same frequency. Under parametric excitation and in the absence of air flow, each mode oscillates at its own natural frequency. In the neighborhood of combination resonance, the nonlinear response is determined using the multiple scales method at the critical flutter speed and at slightly higher airflow speed. The domains of attraction and bifurcation diagrams are obtained to reveal the conditions under which the parametric excitation can provide stabilizing effect. The basins of attraction for different values of excitation amplitude reveal the stabilizing effect that takes place above a critical excitation level. Below that level, the response experiences limit cycle oscillations, cascade of period doubling, and chaos. For flow speed slightly higher than the critical flutter speed, the response experiences a train of spikes, known as ‘firing,’ a term that is borrowed from neuroscience, followed by ‘refractory’ or recovery effect, up to an excitation level above which the wing is stabilized. The results of the multiple scales method are verified using numerical simulation of the original nonlinear differential equations.