New Developments in Fluid Mechanics and Its Engineering ApplicationsView this Special Issue
Research Article | Open Access
Nonlinear Characteristics of Helicopter Rotor Blade Airfoils: An Analytical Evaluation
Some results are presented about the study of airloads of the helicopter rotor blades, the aerodynamic characteristics of airfoil sections, the physical features, and the techniques for modeling the unsteady effects found on airfoil operating under nominally attached flow conditions away from stall. The unsteady problem was approached on the basis of Theodorsen's theory, where the aerodynamic response (lift and pitching moment) is considered as a sum of noncirculatory and circulatory parts. The noncirculatory or apparent mass accounts for the pressure forces required to accelerate the fluid in the vicinity of the airfoil. The apparent mass contributions to the forces and pitching moments, which are proportional to the instantaneous motion, are included as part of the quasi-steady result.
The most important component of the helicopter is the main rotor for which there is a great deal of activity in developing new and improved mathematical models that predict the flow physics. A high tip speed gives the rotor a high level of stored rotational kinetic energy and reduces the rotor torque required for a given power, but there are two important factors that work against the use of a high tip speed: compressibility effects and noise.
The additional effects of compressibility on the overall rotor profile power requirements, when the tip of the advancing blade approaches and exceeds the drag divergence Mach number, were estimated using blade element theory combined with the airfoil section characteristics . A more detailed analysis of compressibility effects on the rotor must represent the actual nonlinear airfoil characteristics as functions of Mach number through stall at each blade element followed by numerical integration.
The classical unsteady aerodynamic theories describing the observed behavior have formed the basis for many types of rotor analysis. The tools for the analysis of 2D, incompressible, and unsteady aerodynamic problems were extended to compressible flows, being a basis for developing linearized unsteady aerodynamic models applicable to compressible flows . But, while the classical theories assume linearity in the airloads, the assumption of linearity can probably be justified for many of the problems encountered on the rotor, in practice. The advent of nonlinear methods based on CFD solutions to the Euler and Navier-Stokes equations has provided new results that justify and define the limits of the linear models and may give guidance in developing improved and more practical unsteady aerodynamic models for future use in helicopter rotor blade airloads prediction, aeroelastic analysis, and rotor design . At the blade element level, the various sources of unsteady effects can be decomposed into perturbations to the local angle of attack and velocity field. At low angle of attack with fully attached flow, the various sources of unsteady effects manifest as moderate amplitude and phase variations relative to the quasi-steady airloads. At higher angles of attack when time-dependent flow separation from the airfoil may be involved, a phenomenon characterized by large overshoots in the values of the lift, drag, and pitching moments relative to the quasi-steady stall values may occur .
The helicopter rotor airfoil must assure a high maximum lift coefficient, a high drag divergence Mach number, a good lift-to-drag ratio over a wide range of Mach number, and a low pitching moment. At high angles of attack, the adverse pressure gradients produced on the upper surface of the airfoil result in a progressive increase in the thickness of the boundary layer and cause some deviation from the linear lift versus angle of attack behavior. On many airfoils, the onset of flow separation and stall occurs gradually with increasing angle of attack, but on some airfoils (those with sharp leading edges) the flow separation may occur suddenly.
2. Helicopter Rotor Blade Aerodynamics
The region of the rotor disk affected by compressibility effects is shown in Figure 1 and is defined on the surface where the incident Mach number of the flow that is normal to the leading edge of the blade exceeds the drag divergence Mach number, . If is the hover tip Mach number, then the region of the disk affected by compressibility effects is defined by
The angular or rotational speed of the rotor is denoted by , the rotor radius by , the advanced ratio by and the radial distance from the rotational axis by . In the above equation the nondimensional quantities are and , where is the angle between the forward velocity and the plane of the rotor [1, 2]. The azimuth angle for the onset drag divergence, , can be obtained by setting , so that and .
The increment in profile power associated with this region on the disk is where is the extra drag on the blade section when it exceeds the drag divergence Mach number, , and is rotor solidity coefficient which represents the ratio of the blades area to the rotor disk area. For the NACA 0012 airfoil, Prouty  suggests that this can be approximated by The rotor limits may be determined by two conditions, one condition given by advancing blade compressibility effects and the other condition given by retreating blade stall. In either case the advancing blade operates at low angle of attack (AoA) but at high subsonic or transonic conditions, whereas the retreating blade operates at low Mach numbers and high lift coefficients.
The aerodynamic behavior of airfoils in the high AoA regime is important for predicting the adverse effects produced in the reverse flow regime on the rotor. In the reverse flow region, the direction of the relative flow vector changes from the trailing edge toward the leading edge of the airfoil. While the fundamental process of the blade wake and tip vortex formation is similar to that found with a fixed wing, one difference with helicopter tip vortices is that they are curved, and so they experience a self-induced effect. Another complication with helicopter rotors is that the wakes and tip vortices from other blades can lie close to each other and to the plane of blade rotation, and so they have large induced effects on the blade lift distribution.
If the wake is assumed to be undistorted in the tip path plane and no wake contraction occurs in the radial direction (Figure 2), then the tip vortex trajectories are described by the equations where is the position of the blade when the vortex was formed and is the position of the vortex element relative to the blade.
One important parameter used in the description of unsteady aerodynamics and unsteady airfoil behavior is the reduced frequency, , defined as , where is the angular frequency, is the chord of the airfoil, and is the flow velocity . According to the dimensional analysis, the resultant force, , on the airfoil chord can be written in functional form as . For , the flow is steady, and for , the flow can be considered quasi-steady; that is, unsteady effects are generally small. Flows with characteristic reduced frequencies above 0.05 are considered unsteady [5, 6]. For a helicopter rotor in forward flight (Figure 3), the local sectional velocity, which appears in the denominator of the reduced frequency expression, is constantly changing.
The approach to modeling of unsteady aerodynamic effects through an extension of steady, 2D thin airfoil theory gives a good level of analysis of the problem and provides considerable insight into the physics responsible for the underlying unsteady behavior. The Laplace’s equation for incompressible flow is eliptic; therefore, the unsteady aerodynamic theories cannot be obtained in a corresponding analytical form.
3. The Apparent Mass Tensor
The rate of change of the impulse vector, in general, is not in the direction of the acceleration of the body. The external force applied to the body to translate it through the fluid has to be applied in a direction different from that of the acceleration of the body through the fluid . Physical conditions that should be satisfied on given boundaries of the fluid (boundary conditions) depend on the assumptions made with regard to the nature of the fluid, more specifically on the nature of the differential equations that are assumed to govern the motion of the fluid. For a solid-fluid boundary, at each point of the solid-fluid surface, at every instant, the component normal to the surface of the relative velocity between the fluid and the solid must be vanish, , where represents the relative velocity and the normal to the surface (Figure 4).
If the surface is represented by a scalar function of position and time, , then the total time rate of change is zero: on .
The fluid force acting on a rigid body of arbitrary shape translating with a velocity is given by where denotes the surface of the body and is the pressure on the surface of the body. In general, the body may be translating, rotating, and deforming; consequently, the velocity is a function of position on the surface and time. If the body is rigid and is in translatory motion, then is a function of time but uniform over the surface of the body. The mathematical problem is to determine the external force applied to the body to translate it through the fluid. According to Newton’s second law, we have where is the mass of the body. The above equation may be rewritten as or where is the impulse applied on the fluid and The fluid force acting on the body is The integral is related to the circulation around the body, and is the velocity potential.
The unit vector is normal to the cutting planes, the unit vector is tangent to the curve of intersection between the blade element surface and the cutting plane, and the unit vector is tangent to the blade element surface.
The component of the vector in the direction is and the vector corresponding to the surface element may be written as , where , , and .
On the other hand, Since the unit vectors and are normal, it follows that and The scalar product , is the normal distance between the cutting planes of the solid body (Figure 5). If we denote , it follows that Here, is the distance measured along the fixed direction , and is the circulation around the curve of intersection between the body surface and the cutting plane. The limits and denote the extremities of the body measured along the direction . It follows that for motions without circulation the force on the body is given by where the velocity potential is the solution of the system Since the equation and boundary condition for are linear, the solution could have the form where each of the functions , , and is a solution of the equations where , 2, or 3.
In the Cartesian coordinate system, the vectors and have the expressions Because time enters through , it is convenient to set , so the system (21) takes the form With these considerations the impulse becomes The components of the impulse are The surface integral in above equations may be written as follows: According to Green’s theorem, if and are two harmonic functions, then Introducing the symbol , with ; the components of the impulse are therefore given by and the force applied externally to the body is or The coefficients form a set of nine numbers which may be displayed as an array and may be referred to as a virtual mass tensor or virtual masses that need to be added to the mass of the body in order to find the force that must be applied to translate it through the fluid. Introducing the symbol defined by if and for , (30) may be rewritten as For any body, there are three perpendicular directions such that for , so with respect to such axes, (33) becomes The sum represents the apparent mass for translation in the -direction, and the corresponding is the additional apparent mass.
4. The Airloads on an Oscillating Airfoil
The oscillatory motion of the airfoil can be decomposed into contributions associated with angle of attack which is equivalent to a pure plunging motion (Figure 6) and contributions associated with pitching (Figure 7).
A plunge velocity produces a uniform velocity perturbation that is normal to the chord, , and the pitch-rate term produces a linear variation in normal perturbation velocity.
For a pitch rate imposed about an axis at “” semichords from the midchord, then , so that the induced chamber is a parabolic arc.
The problem of finding the airloads on an oscillating airfoil was solved by Theodorsen, who gave a solution to the unsteady airloads on a 2D harmonically oscillated airfoil in inviscid, incompressible flow, with the assumption of small disturbances . Both the airfoil and its shed wake were represented by a vortex sheet with the shed wake extending as a planar surface from the trailing edge downstream to infinity. The assumption of planar wake is justified if the angle of attack disturbances remain relatively small. As with the standard quasi-steady thin airfoil theory, the bound vorticity, , can sustain a pressure difference and, therefore, a lift force. The wake vorticity, , must be force-free with zero net pressure jump over the sheet. According to the Theodorsen’s theory, the solution for the loading on the airfoil surface under harmonic forcing conditions is obtained from integral equation where is the downwash on the airfoil surface. At the trailing edge, , and the airfoil circulation is given by So long as the circulation about the airfoil is changing with respect to time, the circulation is continuously shed into the wake and will continuously affect the aerodynamic loads on the airfoil.
For a general motion, where an airfoil of chord is undergoing a combination of pitching and plunging motion in a flow of steady velocity , Theodorsen’s solution [8–10] for the lift coefficient and pitching moment coefficient corresponding to midchord, , is where is the pitch axis location relative to the midchord of the airfoil, measured in terms of semichord and is the complex transfer function (known as Theodorsen’s function) which accounts for the effects of the shed wake on the unsteady airloads: with , , , and being Bessel functions of the first and second kinds, respectively (Figure 8).
The Hankel functions in above expression are The real and imaginary parts of function are plotted in Figure 9.
It could be appreciated that function serves to introduce an amplitude reduction and phase lag effect on the circulatory part of the lift response compared to the result obtained under quasi-steady conditions .
This effect can be seen if a pure oscillatory variation in angle of attack is considered, that is, , so the circulatory part of the airfoil lift coefficient is given by For , the steady-state lift behavior is obtained; that is, is linearly proportional to . As is increased, the lift plots develop into hysteresis loops, and these loops rotate such that the amplitude of the lift response (half of the peak-to-peak value) decreases with increasing reduced frequency.
These loops are circumvented in a counterclockwise direction such that the lift is lower than the steady value, when is decreasing with time (i.e., there is a phase lag). For infinite reduced frequency, the circulatory part of the lift amplitude is half that at , and there is no phase lag angle.
5. Pure Angle of Attack Oscillations
For a harmonic variation in , that is, , the lift is or, in terms of the lift coefficient, the result is The term inside the square brackets can be considered the lift transfer function, which accounts for the difference between the unsteady and quasi-steady airloads.
The first term inside the brackets is the circulatory term, and the second term is the apparent mass contribution, which is proportional to the reduced frequency and leads the forcing by a phase angle of . The noncirculatory or apparent mass terms arise from the nonsteady term of the velocity potential and account for the pressure forces required to accelerate the fluid in the vicinity of the airfoil.
The normalized lift amplitude is The normalized lift amplitude and phase of lift for pure angle of attack oscillations are presented in Figures 10 and 11, where the significance of the apparent mass contribution to both the amplitude and phase can be appreciated.
At lower values of reduced frequency, the circulatory terms dominate the solution. At higher values of reduced frequency, the apparent mass forces dominate.
6. Pure Plunging Oscillations
For a harmonic plunging motion such as that contributed by blade flapping, the forcing is so that and . Substituting into the expression for the lift and solving for the lift coefficient, it gives The complete term inside the square brackets can be considered as the lift transfer function . The circulatory part of the lift response leads the forcing displacement by a phase angle of . Also, the apparent mass force leads the circulatory part of the response by a phase angle of or the forcing by a phase angle of . The corresponding pitching moment about midchord for this case is
7. Pitching Oscillations
For harmonic pitch oscillations, additional terms involving pitch rate appear in the equations for the aerodynamic response. The forcing is given by and the pitch rate by . In this case, the lift coefficient is The lift amplitude initially decreases with increasing because of the effects of the shed wake; and then, for , it begins to increase, as the apparent mass forces begin to dominate the airloads. This is also shown by the phase angle, which exhibits an increasing lead for .
Von Karman and Sear analyzed the problem of a thin airfoil moving through a sinusoidal vertical gust field, where the gust can be considered as an upwash velocity that is uniformly convected by the free stream. The forcing function in this case is where is the gust frequency. If the gust is referenced to the airfoil leading edge, then and becomes , and if the gust is referenced to the midchord, then and Therefore, , which is equivalent to a phase shift. In this case, the lift coefficient can be written as where is known as Sears function and the gust encounter frequency, , is given by and is the wavelength of the gust.
In terms of Bessel functions, Sears’s function is given by The terms of real and imaginary parts are If the gust is referenced to the leading edge of the airfoil, the function will be called and can be written as The two results are plotted in Figure 12. The peculiar spiral shape of the transfer function arises only when the gust front is referenced to the midchord of the airfoil.
The kinematics of the pitching and plunging airfoil of a typical blade element is the resultant of a combination of forcing from collective and cyclic blade pitch, twist angle, elastic torsion, blade flapping velocity, and elastic bending. At low angles of attack with fully attached flow, the various sources of unsteady effects manifest primarily as moderate amplitude and phase variations relative to the quasi-steady airloads. At higher angles of attack when time-dependent flow separation from the airfoil may be involved, the dynamic stall may occur. The amplitude and phase effects produced by the stalled airloads can lead to various aeroelastic problems on the helicopter rotor that may limit its performance. The need to control the aerodynamic forces on the rotor requires that the pitch of each blade be changed individually as the blades rotate about the shaft.
The first flap frequency of a helicopter rigid blade is about , then the reduced frequency at the 75% radius location is . For a blade aspect ratio , the reduced frequency is in excess of 0.07, which is in the unsteady range. In the case of the first elastic torsion mode, which is typically about at the blade tip, the reduced frequency is in excess of 0.2. At these reduced frequencies, there is a significant amplitude and phasing introduced into the airloads by the effects of the unsteady aerodynamics.
When a wing’s angle of attack is increased rapidly, it can momentarily generate a higher maximum lift coefficient than it could if the angle of attack was increased slowly. This overshoot can be related to the change in angle of attack during the time required for the air to travel one chord length. The dynamic overshoot is attributed to two effects (for the airfoils that stall first at the leading edge): the delay in the separation of the boundary layer and the momentary existence of a vortex shed at the leading edge after the boundary layer does separate. The delay in separation corresponds to the finite time required for the aft edge of the separation bubble to move forward to its bursting position. On the other hands, the airfoil can generate high lift as a result of a vortex that is shed at the leading edge at the instant of stall. The vortex travels back over the top of the airfoil carrying with it a low pressure wave that accounts for the very large lift coefficient. Airfoils that stall first at the trailing edge also exhibit a dynamic overshoot but considerably less than those airfoils that have leading edge stall.
- J. G. Leishman, Principles of Helicopter Aerodynamics, Cambridge University Press, 2007.
- J. Seddon and S. Newman, Basic Helicopter Aerodynamics, AIAA Education Series, 2001.
- M. A. Ashraf, J. Young, and I. C. S. Lai, “Oscillation frequency and amplitude effects on plunging airfoil propulsion and flow periodicity,” The American Institute of Aeronautics and Astronautics Journal, pp. 2308–2324, 2012.
- R. Prouty, Helicopter Performance, Stability and Control, Krieger, Gainesville, Fla, USA, 2002.
- P. D. Gareth, Helicopter Flight Dynamics, AIAA Education Series, 2007.
- B. J. Lee and M. S. Liou, “Unsteady adjoint approach for design optimization of flapping airfoils,” The American Institute of Aeronautics and Astronautics Journal, pp. 2460–2491, 2012.
- K. Karamcheti, Principles of Ideal-Fluid Aerodynamics, Krieger, Gainesville, Fla, USA, 1966.
- T. Theodorsen, “General theory of aerodynamic instability and the mechanism of flutter,” Tech. Rep. 496, NACA, 1934.
- I. E. Garrick, “On some reciprocal relations in the theory of nonstationary flows,” NACA Report 629, 1938.
- J. Katz and A. Plotkin, Low Speed Aerodynamics, Cambridge University Press, 2010.
- C. Rotaru, I. Circiu, and M. Boscoianu, “Computational methods for the aerodynamic design,” Review of the Air Force Academy, vol. 2, no. 17, pp. 43–48, 2010.
- C. Rotaru and V. Jaulin, “Numerical simulations of Helicopter rotor blade airloads,” MTA Review, vol. 22, no. 2, pp. 81–90, 2012.
Copyright © 2013 Constantin Rotaru. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.