Modern Helicopter Aerodynamics
Modern Helicopter Aerodynamics
MODERN HELICOPTER
AERODYNAMICS
A. T. Conlisk
Department of Mechanical Engineering, The Ohio State University, Columbus,
Ohio 43210-1107
KEY WORDS: rotor aerodynamics, vortex wakes, tip-vortex, computational fluid dynamics,
experiments, dynamic stall, blade-vortex interaction
ABSTRACT
Modern helicopter aerodynamics is challenging because the flow field gener-
ated by a helicopter is extremely complicated and difficult to measure, model,
and predict; moreover, experiments are expensive and difficult to conduct. In
this article we discuss the basic principles of modern helicopter aerodynamics.
Many sophisticated experimental and computational techniques have been em-
ployed in an effort to predict performance parameters. Of particular interest is
the structure of the rotor wake, which is highly three-dimensional and unsteady,
and the rotor-blade pressure distribution, which is significantly affected by the
strength and position of the wake. We describe the various modern methods of
computation and experiment which span the range from vortex techniques to full
three-dimensional Navier-Stokes computations, and from classical probe meth-
ods to laser velocimetry techniques. Typical results for the structure of the wake
and the blade pressure distribution in both hover and forward flight are presented.
Despite the complexity of the helicopter flow, significant progress has been made
within the last ten years and the future will likely bring marked advances.
1. INTRODUCTION
For over 40 years the helicopter has played an important role in both military
and civilian air transportation. In this article we discuss the basic principles
of modern helicopter aerodynamics. In the past, the term “helicopter aerody-
namics” has been used synonymously with rotor-blade aerodynamics; for the
most part, in this article we will consider this to be the case. However, as will
515
0066-4189/97/0115-0515$08.00
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
516 CONLISK
2. OVERVIEW
The flow past a helicopter is particularly complicated for several reasons. First,
unlike the case of flow over a fixed wing which can often be analyzed by lin-
ear aerodynamics, the flow past a rotary wing is never what aerodynamicists
consider to be “linear”. This poses significant problems in modeling since
numerical simulations need to be iterative in character and experimental obser-
vations of highly nonlinear phenomena are often difficult to interpret because of
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
518 CONLISK
Advancing Side
ΩR
Blade
R V
ψ
ψ=180 o ψ=0
o
Hinges
Rotor Disk
Retreating Side
(a)
α
Ω
Pitch angle θ T L
Lag angle
θ=θ(t) θ
D
Flap angle
θ=α+φ
(b)
(c)
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
advancing blade side), the velocity near the blade is large and since the lift is
proportional to the velocity squared, the angle of attack need not be large to
achieve sufficient lift. On the other hand, as the blade moves in a direction
opposite to the direction of flight (the retreating blade side) the relative velocity
is smaller and the angle of attack must thus be larger to achieve the same total
lift. Thus, without a moment-balancing mechanism, the helicopter would tend
to roll. To balance the forces and moments, the rotor needs to be trimmed; that
is, the angle of attack of the blades on the advancing and retreating sides must
be adjusted periodically throughout each blade rotation cycle so that there is
a balance of moments. This is called cyclic pitch. The collective pitch of the
blades is a control in which the angle of attack (AOA) of each of the blades is
increased simultaneously to achieve a higher lift; an increase of the collective
pitch, for example, results in climb. In hover, theoretically, trim and flap are not
required to balance forces on an isolated rotor; however, non-uniformities and
the presence of the fuselage make them necessary. In addition, rotor blades are
twisted and often tapered; a twisted blade is one in which the local geometric
pitch angle varies along the span.
To provide trim capability and for aeroelastic stress relief, helicopter rotors
are often hinged in the sense that the rotor blades must be permitted to bend out
of the rotor disk plane as well as pitch to satisfy trim requirements; a sketch of a
simple hinging mechanism is also depicted in Figure 1(c). There are two modes
in which the rotor is hinged; the lead-lag hinge permits motion of the blade
within the rotor-disk plane. The flapping hinge permits the flapping motion
of the blades out of the rotor-disk plane. A rotor having both types of hinges
is said to be fully articulated. When the blades flap, they no longer trace out
a single planar “disk.” In this case we speak of a tip-path plane which is the
plane whose boundary is defined by the trajectory of the blade tips.
Rotor blades have a large span-to-chord ratio and thus severe stresses can be
communicated to the hub if the blades are not permitted to flap. However, if
the blades are aeroelastically soft, then hub stresses can be kept to a minimum
and both types of hinges can be eliminated. In such cases, the rotor is said to be
hingeless. Cyclic pitch changes result in changes in the flapping motion. Blade
aeroelastic effects play a major role in determining helicopter performance;
blade and helicopter aeroelasticity is discussed in Johnson (1980).
The complexity of the flow induced by a helicopter is illustrated by the
presence of so many fundamental fluid dynamic research problems. A sketch
←−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
Figure 1 A single rotor blade in forward flight. (a) Advancing and retreating sides of the rotor
disk. (b) Definition of lift, drag and thrust in hover or vertical climb and lag and flap angles. (c)
Sketch of a typical hinge system of a fully articulated rotor; sketch from Johnson (1980).
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
520 CONLISK
Figure 2 A summary of specific flow problems which occur on a helicopter. From Caradonna
(1992).
Figure 3 Sketch of a helicopter rotor wake for a single blade. From Gray (1956).
the tip-vortices, and the airframe. Note that the sense of circulation of the
inboard sheet is opposite to that of the tip-vortex so that unsteady interaction
between the two will occur. There is also a root vortex (not shown in Figure 2)
which emanates from the inboard edge of the rotor blade; however, because of
the relatively low vorticity near the root, this region is usually not a large factor
in design. In addition, there is a wake shed from the rotor hub; hub drag can
be a significant portion of the overall drag. However, for brevity we will not
discuss the hub flow in this review.
The primary task in rotorcraft aerodynamic design is to determine the lift
and drag coefficients of the rotor blades because these two quantities determine
the thrust and power required for given speed in forward flight or hover. There
are two components to the drag: pressure or form drag, and viscous drag.
In situations where loads are generated by three-dimensional vortex systems,
the pressure drag is usually called induced drag. Lift is comparatively easy to
predict because it is usually found from a surface pressure integration, although
this is not the case when the influence of blade-vortex interaction is strong. On
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
522 CONLISK
the other hand, the power loss due to drag is very hard to predict because
it is a much smaller force and is thus sensitive to small changes in pressure
(Ramachandran et al 1989).
From this discussion it is seen that the flow past a helicopter rotor blade
features a wide range of velocities from low subsonic speeds to the transonic
regime. Moreover, important length scales range from the blade length to the
size of the vortex core and thickness of the inboard sheet; these length scales
can span several orders of magnitude. Thus modeling and experimentation
of helicopter flows are extremely challenging, time-consuming, and costly.
Because of these complexities it is difficult to incorporate the dynamic nature
of the entire rotor flow in the presence of the helicopter airframe in one single
numerical computation or experimental measurement program. For this reason,
rotorcraft research tends to be focused on one or two specific aspects of the rotor
flow field and tends to have both experimental and computational components.
For example, many computational and experimental approaches have focused
on the rotor wake flow for the case of two or four rigid blades rotating at relatively
low tip-speeds. Under these conditions, it is often not difficult to obtain good
results for the blade pressure distribution. On the other hand, at high tip-speeds
under forward flight and descent conditions this is often not possible, and a
much more fundamental understanding of these flight conditions is required.
Momentum Theory
For both hover and climb (or descent), the analysis of the mechanics of the
helicopter began by drawing an analogy with the study of propellers. In the
mid-ninteenth century, theories were developed to meet the steady growth of
the ship propeller industry. Rankine (1865) developed a simple model of a
propeller flow field by applying linear momentum theory derived from the
basic relationships of Newtonian mechanics. Subsequently, this early theory
was applied to rotorcraft.
During hover, which is the simplest helicopter flight regime, the rotor pro-
duces an upward thrust by pushing a column of air downwards through the
rotor-disk. If the flow is assumed to be steady, inviscid, and incompressible,
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
524 CONLISK
For small angles of attack, exceptionally simple formulas for C T can be deduced.
The power coefficient is defined in terms of the torque produced in rotating
the blades
P
CQ = ,
ρ A(R)3
and following a similar procedure to that described above, the result is
Z 1
1
CQ = σ (λC L r 2 + C D r 3 )dr, (6)
2 0
vi
where λ is termed the inflow factor and for the case of hover is given by λ = R .
√
Note that in hover, from the definition of the thrust coefficient and the induced
velocity, λ = C2T . The power coefficient, which is the measure of how much
power is required to produce lift and to rotate the blades, depends crucially on
the drag coefficient.
From this discussion, it is evident that rotor performance depends critically
on sectional flow properties; namely, local blade angle of attack and local
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
526 CONLISK
wake is not truly predictive. The prescribed wake concept has been extended to
numerical prediction of rotor airloads for a range of flight conditions by Egolf
and Landgrebe (1983).
In a free wake calculation, the vortex system motion is calculated directly
from the effects of all the other wake components and the influence of the
blade. In this method, the wake is allowed to develop in time and initial-value,
Lagrangian methods are used to determine its position at each time step. This
is the industry standard today, and as computers have become faster, the free
wake calculation has become more affordable.
There are two means of determining the structure of a given flow field. In a
Lagrangian description, individual fluid particles are followed forward in time.
In an Eulerian description (not to be confused with the Euler equations), the
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
Figure 5 Vortex lattice representation of the rotor wake as described by Caradonna (1992).
solution for the velocity field is calculated at a number of specified field points.
Thus an Eulerian description involves a grid system which is not required in a
Lagrangian description. In what follows, we consider both sets of methods for
calculating the rotor wake.
528 CONLISK
In equations (7) and (8) 0 is the circulation, r is the radial coordinate, and av
is the vortex core radius. Use of these profiles in rotorcraft calculations has
recently been discussed by Leishman et al (1995); their results indicate little
difference in the profiles when compared with experimental data, although the
Rankine profile seems to overpredict the maximum velocity in the vortex when
compared with data (Leishman et al 1995).
Scully (1967) used this approach to numerically compute solutions for the
rotor wake. Using the swirl velocity given by equation (8), he used the straight
line segments of Figure 5 to model the tip-vortex. The inboard vortex sheet
was represented by a single large-core vortex located at about midspan. It was
found that the tip-vortex is the dominant component of the wake velocity field
because its strength is much greater than that of the inboard sheet; this can be
seen from the sketch of the bound circulation depicted in Figure 5.
As noted above, the free wake calculation is now the industry standard. This
means that the wake is allowed to evolve based on the velocity felt at each point
on the tip-vortex and vortex sheet. Once the velocity distribution is determined
at each point, the wake must be advanced forward in time using an initial-value
problem solver. Typical methods used include first-order Euler, Euler predictor-
corrector, and more accurate methods. Because the Biot-Savart Law is only
valid for incompressible flow, rotor codes generally employ the Prandtl-Glauert
compressibility correction for local sectional aerodynamic loading (Johnson
1980, pp. 262–263).
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
The advantage of vortex methods is that the wake structure can usually be
calculated without numerical dissipation because the size of the vortex core is
specified. However, the strength of the wake as measured by the strength of
the tip-vortex and the inboard sheet is not known a priori. A common way to
specify this is to view the rotor blade as a source of lift in which the spanwise
concentration of vorticity shed into the wake is concentrated on a single line,
usually the quarter-chord of the airfoil. This is called lifting-line theory. This
method is simple and adds little computational cost since, for the linearized
case, an analytical solution for the circulation is available. The method for
a fixed wing in incompressible flow is described in aerodynamics text books
(Katz and Plotkin 1991) and the extension to compressible flow for a rotary
wing is straightforward (Johnson 1980). For incompressible flow, the amount
of bound circulation generated by the lifting line can be calculated by satisfying
the Kutta condition (just as for a fixed wing). The result for a small angle of
attack at each two-dimensional slice of the blade inboard of the blade tip is
0 = π cyα (9)
where α is the angle of attack measured from the zero-lift angle of attack. The
sectional circulation corresponds to the fixed-wing result with external flow
speed U = y and is depicted in Figure 5; of course, equation (9) is not valid
close to the blade tip.
While lifting-line theory is computationally simple and fast, it is well known
that it is not valid for large variations in the downwash velocity caused by the
vortex wake. To remedy this, some rotor codes use a more accurate lifting-
surface method in which the rotor blade is viewed as a wing of zero thickness
to calculate the circulation (Johnson 1971). However, if the angle of attack
is not small and the blade is near stall, the lift coefficient must often be ob-
tained from experimental data or from a simplified lifting model (Leishman
& Beddoes 1989). Airfoil lookup tables are also used to provide the lift and
drag coefficients. The static stall angle for many rotor blades is generally in the
range of 12–15◦ although dynamic stall angles may be above 20◦ .
Egolf (1988) has used a modified lifting-line theory along with a vortex box
model to calculate the wake for a forward flight condition. The wake is divided
into a number of square “boxes” whose boundaries are vortex filaments of equal
circulation. Results for a number of quantities such as local section lift, blade
circulation, and BVI-induced blade loading are presented. Bliss et al (1987)
describe the implementation of the curved vortex elements of Figure 5 first
described by Bliss et al (1983). The motivation of this approach is that curved
vortex elements are more natural to describe the wake; they indicate that, using
this approach, far fewer vortex segments are required than for the straight-line
case.
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
530 CONLISK
Torok and Berezin (1993) compare several methods for computing wake-
induced airloads. The results are depicted in Figure 6. Figures 6a,b,c depict
results for a vortex lattice approach using straight line segments (Figure 6a),
a rigid inner wake and a free tip-vortex (Figure 6b), and a curved vortex line
approach based on the constant vorticity contour approach (Figure 6c), respec-
tively. The constant vorticity contour method (Bliss et al 1987) consists of
discretizing the wake with curved vortex elements along lines of constant vor-
ticity. Note the considerable differences in the wake geometry for each of the
three methods. The differences are particularly acute near the blade tips where
the vorticity field is strongest. For Figures 6a,b only the tip-vortex is shown
although the entire wake is calculated; for Figure 6c the constant vorticity con-
tours along the entire blade are shown. Airload results (not shown) indicate
that the constant vorticity contour method compares better with experimental
data than the others, especially at low advance ratio, but all the results for the
blade loads exhibit significant deviation from experimental data.
Lifting-line and lifting-surface methods for the representation of the body
are special cases of a more general class of methods referred to as panel meth-
ods. These methods consist of solving for the velocity potential in terms of
the fundamental solution of Laplace’s equation, which can be a source, vor-
tex, or doublet. The solution is evaluated on the body and the integrals are
discretized to lead to a system of linear equations for the unknown elemen-
tal strengths which may be solved using standard techniques. Panel meth-
ods have been used extensively in helicopter aerodynamics for a number of
years. The primary advantage of these methods is that, theoretically, arbi-
trary body shapes may be considered. Panel methods are usually coupled with
a far-wake model or with a grid-based full potential method (see below, in
Section 4). The use of panel methods in computational fluid dynamics is re-
viewed by Hess (1990); applications to rotary wing aerodynamics are discussed
by Morino and Gennaretti (1991) and by Caradonna (1992). The method
has been used by Maskew (1986) and a new boundary-element formulation
of the rotor blade and wake has been described by Gennaretti and Morino
(1992).
The methods described so far involve the calculation of the wake-induced
airloads in which the wake vortex system is described in terms of the Biot-Savart
integral and/or superposition of related elementary singularities. As such, no
computational grid in the wake is required, and for this reason the methods
described here are normally not viewed as computationally intensive. With the
advances in speed and memory of computational hardware, increasingly larger
problems may be solved with a finite-difference or finite-volume approach in
which the velocities in the wake are calculated directly. These methods are
described next.
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
532 CONLISK
(a)
(b)
(c)
Figure 6 Comparison of various wake models which utilize vortex methods to describe the rotor
wake. The tip-vortex only is shown. (a) Vortex lattice, rigid blade; (b) Scully wake, straight line
segments; (c) constant vorticity contour method. From Torok and Berezin (1994).
534 CONLISK
Figure 7 Computed wake geometry at rotor phase angle 90 and 180 for advance ratio µ = 0.19
using the full potential equation with an adaptive grid. No wake model is employed. From
Ramachandran et al (1993).
are incorporated and the grid dynamically adjusts to the motion of the blade.
Unlike much of the work using a grid-required numerical method, there is no
external wake model required. The wake geometry at two rotor azimuths after
six revolutions is depicted on Figure 7. The strength of the tip-vortex is seen by
the fact that the wake lines bunch near the blade tips; however, the tip-vortex
and the inboard sheet do not appear to be distinctly separate as in Figure 1.
performance variables such as the thrust and power coefficients. Two relatively
focused reviews of these methods have been published. Landgrebe (1994) out-
lines the primary contributions of Navier-Stokes and Euler calculations in the
US. A very recent assessment of Euler and Navier-Stokes methods and pos-
sibilities of developing new methodologies has been given by Srinivasan and
Sankar (1995).
There are several important elements to the calculation of solutions to the
Navier-Stokes equations in rotorcraft applications. First, some sort of body-
fitted grid generation module is required; this is usually done by some prescribed
means, often as the solution to a linear partial differential equation. The gen-
erated grids need to be fine in regions where the velocity varies rapidly; these
regions include the blade, nose, and tail regions at the tip of the blade and in
the viscous boundary layer on the blade. The grid may be structured or un-
structured; an unstructured grid is a grid system in which the nodal points are
specified in an arbitrary manner. This leaves the wake flow to be covered by
a grid system in which the local mesh size is relatively large compared with
the thickness of the inboard vortex sheet and the tip-vortex. For this reason,
in grid-based computations of the rotor wake, the wake is often smeared out
and the location of the wake is sometimes hard to pinpoint. This is in contrast
to vortex methods in which the inboard sheet and the tip-vortex structure are
specified to a large extent.
Second, to discretize the convective terms in the Navier-Stokes equations,
some form of upwinding is required; typically third-order upwinding is used
although fifth-order upwinding is now becoming popular in an effort to preserve
accuracy (Bangalore and Sankar 1996). Finally, some time-advancing scheme
is required and this can be either explicit or implicit. Typically, modern rotor
codes employ a first- or second-order implicit time advance algorithm. There
is generally an option to calculate solutions to the Euler equations rather than
Navier-Stokes; in general, the Navier-Stokes solutions compare better with
experiment. While not explicitly discussed in this review, rotor codes include
a turbulence model; different turbulence models can lead to different results
(Srinivasan et al 1995).
The Navier-Stokes equations for an isolated rotor in hover are usually solved
as an initial-value problem with the rotor started from a given state; the equations
need to be integrated to steady state. This is very difficult, especially in hover
and low-speed forward flight, because of the number of grid points required to
resolve the flow for a relatively large number of time steps. Typically, more
than two turns of the rotor are required to establish steady state, and by this
time, due to numerical diffusion, the strength of the tip-vortex and the inboard
sheet are much weaker than those seen in experiments. This is especially true
at the high disk loadings required in helicopters today.
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
536 CONLISK
Complicating the computation is the need to trim the rotor. Many of the
Navier-Stokes computations described here are, at most, only partially trimmed
using an external comprehensive rotor design code. At the present time, com-
putational time limitations prevent adding a trim module to the codes. Cyclic
and collective pitch settings are calculated from the comprehensive rotor code
and then input to the CFD rotor blade code and the flow field is computed. The
thrust and power coefficients may then be calculated and the values reinserted
into the rotor design code in which the rotor is trimmed again. This process
may be continued until the rotor and trim solutions are compatible.
Wake and Sankar(1989), in a paper first presented at the American Helicopter
Society National Specialists’ meeting on Aerodynamics and Aeroacoustics in
1987, were the first to present solutions to the unsteady Navier-Stokes equations
for a rigid rotor blade. They solve the compressible three-dimensional Navier-
Stokes equations and compare their computational results with the experimental
data of Caradonna and Tung (1981) in the nonlifting, transonic regime for an
ONERA blade; results are also produced for a lifting, NACA-0012 blade. A C-
grid is used to describe the domain near the blade. The numerical procedure to
solve the equations is a fully implicit procedure in space, and is based on a Beam
and Warming scheme (Beam and Warming 1978). The far-wake field is either
extrapolated from the interior of the computational domain or set equal to zero.
Rotor inflow conditions are specified by the transpiration velocity technique
which requires the effective angle of attack of the wake; this parameter is fixed
externally using a rotor design code. The computational results for the blade
pressure in hover are depicted in Figure 8a. Forward flight results are also
presented, and the comparisons with the experimental results for rotor phase
angles in which the flow remains subsonic are good; at higher free-stream Mach
numbers where the flow may be locally supersonic over a good portion of the
blade, the agreement is not as good.
Srinivasan and McCroskey (1988a) solve the unsteady thin-layer Navier-
Stokes equations for the same conditions as Wake and Sankar (1989). The
thin-layer Navier-Stokes equations are a subset of the Navier-Stokes equations
in which the streamwise and blade-spanwise derivatives are neglected in the
viscous terms. A typical result for the pressure distribution is depicted in
Figure 8b; these results are very similar to those of Wake and Sankar (1989)
thus indicating that, at least for the blade pressure distribution, solutions for
thin-layer Navier-Stokes equations do not differ significantly from solutions
for the full Navier-Stokes equations.
As mentioned earlier, all current Navier-Stokes calculations of the rotor wake
suffer from numerical diffusion in the sense that the vortex system is consid-
erably smeared. Thus, the fact that the blade pressure distributions depicted in
Figure 8 seem to agree well with experiment is somewhat surprising. However,
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
(a)
(b)
Figure 8 Surface pressure distribution for a lifting rotor in hover. Mti p = 0.44, collective pitch
8◦ , Re = 106 for the experimental data of Caradonna and Tung (1981). (a) From Wake and
Sankar (1989) using Navier-Stokes. (b) From Srinivasan and McCroskey (1988) for thin-layer
Navier-Stokes; the solid line is the computation.
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
538 CONLISK
note that there are significant relative errors in the pressure especially near the
tip and at the root. Moreover, it must be pointed out that the results of Figure 8
depict pressures at only one section of the rotor blade; these local section errors
often lead to large errors in the integrated lift and normal force coefficients and
pitching moments. Industry design requirements suggest that it is necessary
for computations to agree with experiments to about 1% for blade loads. Sig-
nificant improvements to the computational scheme described by Srinivasan
and McCroskey (1988a) were made and reported by Srinivasan et al (1992).
Wake and Egolf (1990) have formulated the problem for use on a massively
parallel (SIMD) machine using thousands of processors. Massively parallel
architecture is a means for making these computations more affordable.
The influence of the far-field boundary conditions used in a given Navier-
Stokes computation have been discussed by Srinivasan et al (1993). The usual
boundary condition for an isolated rotor in hover is to assume that outside a
suitable computational cylinder, the velocity vanishes. This means that the fluid
merely recirculates within the computational box; clearly, this seems unphysical
in the sense that the rotor continuously draws fluid into the rotor-disk from
outside the box. To remedy this inconsistency, Srinivasan et al (1993) model
this process with a three-dimensional sink to satisfy mass flow requirements.
A sketch of this boundary condition is depicted in Figure 9a; the influence of
this new boundary condition is depicted in Figure 9b where the sectional thrust
distribution is depicted. Note that the results for the new boundary condition
are in better agreement with experiment especially near the blade tip where
errors are normally larger.
Solutions for the forward flight regime have also been produced. For this
condition, blade-vortex interaction may occur and, even in the absence of
significant blade-vortex interaction, forward flight computations are much more
difficult than hover. Srinivasan and Baeder (1993) have produced solutions for a
forward flight condition of µ = 0.2 and a blade tip Mach number of Mti p = 0.8.
In this case, the flow is locally supersonic near the nose of the blade and so a
shock forms. The presence of the shock is indicated by the very large surface
pressure gradient on the upper surface of the blade as shown in Figure 10. Both
Euler and Navier-Stokes solutions are computed using the Baldwin-Lomax tur-
bulence model; the agreement with experiment is similar to that of Figure 8.
The good agreement with the Euler solution just aft of the shock is probably
fortuitous.
Ahmad and Duque (1996) include moving embedded grids in their calcula-
tion of the solution for the rotor system of the AH-1G helicopter. The embedded
grid procedure allows a more efficient and more accurate calculation of the flow
near the blade under both pitching and flapping conditions; the rotor is partially
trimmed externally. The time-accurate calculation is started from freestream
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
(a)
(b)
Figure 9 Effect of far field boundary conditions on the rotor wake calculation for a single rotor.
(a) Schematic of the new set of boundary conditions. (b) Sectional thrust distributions for the
UH-60 rotor; here Mti p = 0.63, θc = 9◦ , Re = 2.75 × 106 . From Srinivasan et al (1993).
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
540 CONLISK
Figure 10 Instantaneous surface pressure distribution at ψ = 120◦ for a nonlifting rotor in forward
flight. Here Mti p = 0.8, µ = 0.2, Re = 2.89 × 106 , at y/R = 0.89. From Srinivasan and Baeder
(1993).
conditions and five complete rotor revolutions are calculated. A typical result
for the rotor wake is shown in Figure 11. The wake streaklines exhibit periodic-
ity in about two rotor revolutions; note that while the streakline patterns clearly
show the trajectory of the tip-vortex, the magnitude of the vorticity in the tip-
vortex may be very small, indicating significant diffusion. The blade pressure
and section normal force show significant differences when compared to exper-
iment and the power is overpredicted by 15%. These comparisons are typical
of forward flight computations. The complete unsteady calculation takes a total
of 45 hours of single-processor CPU time on a Cray C-90 supercomputer and
generates 40 Gb of flowfield data.
The results discussed above for the Navier-Stokes equations include a time-
stepping algorithm; as such computational error grows with time and in general,
solution accuracy degrades substantially after only about one turn of the rotor.
Since numerical diffusion increases with time, the accuracy of blade loads is
substantially degraded and this is exacerbated at full scale by the aeroelastic
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
Figure 11 Streaklines for a rotor wake at µ = 0.19. Flow periodicity is established in about two
rotor revolutions. From Ahmad and Duque (1996).
542 CONLISK
(a)
(b)
Figure 12 Surface pressure coefficient for Mti p = 0.628, µ = 0.3, for the UH-60A baseline rotor
at (a) y/R = 0.775, ψ = 30◦ ; (b) on the retreating side at y/R = 0.4, ψ = 320◦ . From Bangalore
and Sankar (1996).
significant vorticity (i.e. the inboard vortex sheet and the tip-vortex), and n̂ is a
unit vector in the direction of the normal to the boundary of the non-zero vorticity
region. With relevance to helicopter aerodynamics, the method has been applied
to incompressible blade-vortex interactions (Steinhoff and Raviprakash 1995)
and shows promise for compressible wake calculations.
Performance parameters can be calculated directly from the CFD calculations
discussed in this section. Tung and Lee (1994) have compared performance
parameters such as the figure of merit and the section thrust and torque coeffi-
cients for several different methods of calculating the rotor wake of an isolated
rotor in hover. Generally the agreement is adequate with the model-scale data
set, although significant differences of up to 10–20% near the blade tip are
observed. Moreover, moments due to drag effects are more difficult to predict:
Little or no experimental data for drag is available for comparison with the
computed results.
The rotor wake is complicated by two additional physical problems which
are exacerbated by the pitching and torsional motion of the rotor blades: These
are blade-vortex interactions and dynamic stall. These are significant issues in
their own right and the problem of dynamic stall has been discussed in reviews
of the problem on fixed-wing aircraft. Here we discuss the problems briefly
within a rotorcraft perspective.
Blade-Vortex Interactions
One of the most difficult problems with helicopter operation is the occurrence
of rotor-blade tip-vortex interaction. Blade-vortex interaction (BVI) is defined
as the interaction between a tip-vortex shed from a given blade and another
following blade and is most severe when the vortex approaches the blade ap-
proximately aligned with the spanwise axis of the blade. This means that the
interaction between the vortex and the rotor blade is nearly two-dimensional.
A sketch of a direct collision is depicted in Figure 13 from McCroskey (1995).
As the blade rotates, pitches, and flaps, the origin of the tip-vortex at the blade
tip varies in position; thus, combined with its self-induced motion, the shed
tip-vortex may pass very near or collide with a following blade. Blade-vortex
interaction is rare in hover, can occur in forward flight, but may be particularly
severe during maneuvers, in vertical descent, and in landing (forward flight
descent). There is a large body of work in the literature on various aspects of
BVI and the noise produced as a result. Due to space limitations we discuss
only briefly the main characteristics of BVI. Reducing the number and intensity
of blade-vortex interactions is critical for reducing rotor noise.
BVI noise is one of the two major components of impulsive noise associated
with the flow past the rotor blades, the other being high-speed impulsive noise
due to the high tip Mach number on the advancing side of the rotor (McCroskey
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
544 CONLISK
1995, Gallman et al 1995, Heller et al 1994, Gorton et al 1995a). BVI can occur
on both the advancing and retreating blade sides, but from an acoustic point
of view, the interactions on the advancing side are more important because of
the higher Mach number there. As pointed out in the review by McCroskey
(1995) and which can be verified by a simple dimensional analysis, the most
important parameters in the blade-vortex encounter are the strength of the vortex
and its distance from the blade. Prediction of the magnitude of the acoustic
signal in the far acoustic field is limited by the accuracy of the calculation of
the location and strength of the tip-vortex; small changes in miss distance can
result in significant differences in the nature of the acoustic field. From this
short discussion, it is evident that a highly accurate computation of the local
blade loads during BVI is necessary for an accurate calculation of the noise
field.
Generally, for a two-bladed rotor, BVI can commence late in the first quadrant
of the blade motion (ψ ∼ 60◦ ) and is completed near ψ ∼ 180◦ . The precise
extent of the influence of BVI depends on a number of factors including forward
flight speed and tip Mach number. Computational models have been developed
for the prediction of the loads during BVI; a few of the many papers include
Caradonna et al (1988), Srinivasan and McCroskey (1988b), and Srinivasan
and Baeder (1992). In these papers, the vortex structure is specified and fixed
Rotor Tip-
Path Plane
Γ
Ω
Γ Ψ = 90
o
BVI
o
Ψ=0
Figure 13 A sketch of a direct collision between a blade and a vortex. From McCroskey (1995).
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
in time, and the position of the vortex, and hence the blade loads, can therefore
be calculated fairly accurately. For example, Caradonna et al (1988) solve the
three-dimensional full potential equations; good agreement with experimental
data was obtained for both parallel and oblique vortex interactions when the
core size of the vortex is fixed in time. However, when the structure of the vortex
is calculated along with the balance of the solution of the rotor wake, the results
are not nearly as good. Moreover, a number of blade-vortex interactions are
expected to occur over a single blade cycle.
Lorber (1991) has experimentally investigated BVI in the forward flight–
descent regime for a four-bladed rotor. He was able to map regions where
BVI is expected; both retreating and advancing BVI events were found. He
was also able to develop a model using a simple wake distortion model and
the agreement with experiment is encouraging although some BVI events in
the experiments were not predicted in the computations. See also Gorton et al
(1995a) for additional experimental results.
A detailed study of multiple self-generated BVI encounters has been given
by Hassan et al (1992). They solve the Euler equations to produce results for a
number of BVI encounters during a single rotor revolution. The solutions are
found to be extremely sensitive to the value of the vortex core radius chosen. A
typical result at two inboard stations is depicted in Figure 14. In these results,
BVI commences at around ψ ∼ 60◦ . Note that agreement with the experimental
values is good until about ψ = 90◦ ; note the differences between vortices of
different core radius. These results are typical of results at other inboard stations
as well and reflect the fact that the surface pressure is extremely sensitive to the
position and structure of the vortex system relative to each blade.
The strength of the tip-vortex as measured by its circulation is a strong
function of the blade tip shape from which it was shed. Consequently, much
research has been done to determine if BVI noise can be reduced by suitable
design of the tip shape (Yu 1995).
Dynamic Stall
The flow field on the retreating blade side of the rotor is much different from
the flow on the advancing side since the rotor must be trimmed. Because
of the need for higher lift on the retreating blade side due to the lower relative
velocity, the effective blade angle can become relatively large and dynamic stall
can occur. The term dynamic stall refers to the process by which a separation
bubble forms on the body, grows in size and eventually leaves the body as it is
shed downstream. This process induces an additional lift which can be useful
but is eventually lost as the stall vortex leaves the surface. At the same time,
a dramatic increase in drag and moment also occurs. This large increase in
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
546 CONLISK
Figure 14 A comparison of predicted and measured surface pressures for the model OLS rotor
under BVI conditions at two different inboard locations; three BVI interactions occur during this
time range. Here R B A R = Ry and rCv is a dimensionless vortex core radius. The tip Mach number
M = 0.67, advance ratio 0.1632. The distributions are at 3% chord. From Hassan, Tung and
Sankar (1990).
moment induces significant vibratory loads. There have been a large number
of studies of dynamic stall on wings and there have been several reviews on
the subject including Carr (1988) and McCroskey (1981). Consequently, in
deference to existing recent reviews, we limit our discussion of this topic here.
It has long been recognized that the flow past the rotor blade is at least
locally two-dimensional at a given blade section. Two-dimensional dynamic
stall calculations for rotorcraft have been common over the years (Ham and
Garelick 1968). Consequently, much of the current work on dynamic stall is in
the two-dimensional arena (Carr and Chandrasekhara 1995).
From a rotorcraft perspective, the stall vortex is an additional source of
vibration of the rotor blade complicating the calculation of the generation and
evolution of the rotor wake. Substantial additions to the data base on finite
wings oscillating in subsonic flow are the experiments of Lorber et al (1994).
They studied three-dimensional dynamic stall on several rotor blades in order
to determine geometric effects on performance and stall boundaries in both
steady forward flight and descent conditions. They measured blade pressures
and identified a regime where stall occurs followed by the detachment of the
stall vortex identified as an aft propagation of its suction peak. Laminar and
turbulent flow regimes are identified as well. Srinivasan et al (1995) discuss
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
5. EXPERIMENTAL METHODS
So far we have discussed modeling efforts aimed at understanding and com-
puting blade loads and the rotor wake. Experimental measurements are the
foundation of helicopter design and are the main validation tool for modeling
strategies. Many of the modeling results discussed so far have included ex-
perimental comparisons. In this section we examine more closely the various
experimental approaches which have been used in the past, as well as those that
are currently used.
The measurement of the flow field around a rotating blade is a significant
challenge even at relatively low rotational speeds. For chordwise resolution and
to produce the interaction of strong vortices with surfaces, rotors should have
diameters of at least one meter. The facility walls must be at least one meter
away from the hub. Model rotor speeds generally exceed 1000 revolutions
per minute. A piece flying off such a rotor will have momentum in the plane
of the rotor-disk comparable to that of a handgun bullet. Thus stress analysis,
vibration diagnostics, control stations away from the tip-path plane, steel plates,
and bulletproof windows are essential parts of a rotorcraft flow measurement.
Rotor blades flap, teeter, pitch, bend, and twist; these motions pose problems for
mechanical probes as well as in locating a point on a rotor blade in time and space
to the accuracy needed for aerodynamic measurement. As shown in Figure 2,
compressible flow is expected near the tip-vortex core, around the advancing
blade tip, during blade-vortex interaction, and around the leading edge during
dynamic stall. Elsewhere, the flowfield is generally incompressible, limiting the
utility of density-gradient visualization techniques. Light-scattering techniques
encounter the tradeoff between seed particle inertia and the expense of optical
systems. Blade-mounted instrumentation such as sensors of surface pressure,
acceleration, and shear stress must perform under extreme radial acceleration,
and still have sufficient signal-to-noise ratio when transferred from the rotating
to the fixed coordinate system. These considerations limit the applicability
of many diagnostic advances demonstrated in fixed-wing aerodynamics and
bench-scale water tunnels.
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
548 CONLISK
There are roughly four different sets of aerodynamic data which are sought
by the helicopter aerodynamicist. These are (1) flow visualization of the rotor
wake by smoke or some other nonintrusive means such as laser doppler ve-
locimetry (LDV) or wide-field shadowgraphy; (2) direct measurement of the
blade (or airframe) pressure distribution; (3) measurement of rotor inflow ve-
locities (experimentally generated rotor inflow conditions are generally used
as boundary conditions for CFD studies); and (4) direct measurement of the
wake velocity field. Methods for acquisition of this data range from intrusive
methods such as hot wire probes, to nominally nonintrusive means such as
laser doppler velocimetry, wide-field shadowgraphy, and schleiren techniques.
The flow in the rotor wake can be visualized using smoke and other seeding
materials in conjunction with a laser light sheet as well. Lorber (1991) has
described the use of several different visualization and measurement tools for
the comprehensive description of helicopter flow fields.
In the 1950s, flight test experiments were used to validate models of the rotor
wake that were based on momentum theory. Most often these measurements
were done at low disk loadings and in hover or low-speed forward flight under
nominally steady flow conditions. A good discussion of these time-averaged
results and comparison with models based on momentum theory is given in the
review of Gessow (1986).
Much of the “prescribed wake” methodology which powered classical anal-
ysis tools was developed using painstaking human analysis of photographs to
accumulate vortex trajectories. Landgrebe (1972) used smoke emitted from
rakes in a single plane for visualization of the wake of a two-bladed rotor. The
results in Figure 15a show the two-dimensional cross-section of the wake. Both
the tip-vortex and the inboard vortex sheet are visible on this figure. Note the
significant wake contraction which had also been predicted using momentum
theory. Landgrebe also demonstrated the use of conventional schleiren and
laser holography.
The smoke visualizations of the rotor wake were used to develop analytical
expressions for the time-averaged radial and axial positions of the helicopter
wake. The axial coordinate of the tip-vortex could be accurately described by
a formula of the form,
2π
z̄ = k1 ψw for 0 ≤ ψw ≤ ,
b
µ ¶
2π 2π
z̄ = z̄ ψw = 2πb + k 2 ψw − for ψw ≥ , (13)
b b
(a)
(b)
Figure 15 Experimental results from Landgrebe (1972). (a) A sample smoke visualization of a
slice of the rotor wake for a two-bladed rotor. (b) Tip-vortex trajectory for a six-bladed rotor for
several tip Mach numbers.
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
550 CONLISK
and b is the number of blades. k1 and k2 are parameters which are obtained
from experiment; √ k1 is found to vary linearly with CσT while k2 is found to
vary linearly with C T . The radial coordinate of the tip-vortex was found to
correlate well with an equation of the form
r̄ = A + (1 − A)eλw ψw , (14)
where A and λw are constants and A = 0.78 and λw = 0.145 + 27C T ; r̄
is a dimensionless distance also normalized on the rotor radius. Similar wake
coordinates are given for the inboard sheet. The rate of descent of an element of
the tip-vortex is substantially constant prior to the passage beneath the following
blade (ψw ≤ 60◦ , Figure 15b) at which time it jumps to another constant descent
speed. These simple wake coordinate formulas are used in prescribed wake
models to predict the rotor airloads. Similar formulas are discussed by Gray
(1992) and a correction to these formulas for two-bladed rotors is given by
Kocurek and Tangler (1976).
A major data set on which many theoretical and computational schemes are
calibrated is the data set produced by Caradonna and Tung (1981). These tests
involved experiments on a two-bladed untapered and untwisted planform of
aspect ratio six and a NACA-0012 airfoil section. They measured blade pressure
distribution and wake geometry; the pressure distribution was measured using
three pressure transducers along each blade span. The wake properties were
measured with a hot-wire probe oriented with the wire being tangent to the
tip-path plane. Several rotor speeds were tested and the two-bladed rotor was
substantially rigid.
The tip-vortex trajectory resembles the results presented by Langrebe (1972),
and the tip-vortex strength was close to the maximum bound circulation on the
blade sketched in Figure 5. In addition, the vortices were found to closely
resemble the classical Rankine vortex structure. However, the measured blade
loading and that computed using a lifting-surface with a prescribed wake do not
match, leading to the conclusion by the authors that the blade spanwise loading
cannot be predicted without accurate vortex location and strength.
Parthasarathy et al (1985) applied wide-field shadowgraphy to the wake of a
tail rotor. This technique has the major advantages of being nonintrusive and
of not requiring seeding. It uses a point source to project the path-integrated
second spatial derivative of density on a retro-reflective screen. Combined with
pulsed lighting and video, this technique enables quantitative analysis of the
projection of vortex trajectories. Generally, depth perception is not available.
Norman and Light (1987) reported measurements of the tip-vortex geometry in
the wakes of two main-rotor models in hover. They observed some scatter in
the data due to rotor wake instabilities as well as a short-wavelength instability
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
552 CONLISK
Figure 16 Velocity components along a horizontal sweep across the trailing vortex. The time
frame shown corresponds to the motion of the rotor blade from 2.5 to 3.5 chord lengths past the
measurement location as noted in the sketch. The rotor speed is 1100 rpm and the blade tip is at
y = 0. From McAlister (1995).
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
6. INTERACTIONAL AERODYNAMICS
The trend in design of modern helicopters is toward more lightweight and
smaller vehicles while retaining payload capability and performance parame-
ters. This often requires that the rotor be closer to the airframe, and the higher
rotor speeds required often result in higher disk loading. Thus, in recent years,
one of the areas of rotorcraft aerodynamics that has received a large amount of
attention is the interactional aerodynamic phenomena which occur on a heli-
copter. A detailed classification of a variety of aerodynamic interactions among
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
554 CONLISK
major components of the helicopter is given by Sheridan and Smith (1980) and
these interactions are summarized in Figure 17. Two of the most important
interactions are main-rotor tail-rotor interactions and wake-fuselage interac-
tions. The intensity of these interactions depends on the relative proximity,
size, shape, and flight conditions (e.g. low speed, in ground effect, flight with
sideslip) of individual components.
As pointed out by McCroskey (1995), tail rotors are employed to balance the
torque of the main rotor, and for better control in hover and low-speed forward
flight. The performance of the tail rotor may be adversely affected if it interacts
significantly with the main rotor wake. Prouty and Amer (1982) discuss the
difficulties, due to poorly understood interaction effects, in designing the em-
pennage and tail rotor. In addition, this interaction can be a significant source of
noise. Additional references on this interaction are given in McCroskey (1995).
The interaction which has perhaps received the most attention is the vortex
wake–fuselage interaction. This interaction process is especially intense in
hover, climb, and low-speed forward flight since in these flight regimes the
wake is transported almost vertically downward and may directly impinge on
the surface of the fuselage. This can generate additional impulsive and periodic
loading which may adversely influence the handling qualities of the aircraft and
affect the fuselage loading. These loads can also be transmitted to the cabin
in the form of a low-frequency vibration of the fuselage shell and thus can be
annoying to passengers.
Figure 17 Summary of the many interactions among the major components of a helicopter. From
Sheridan and Smith (1980).
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
556 CONLISK
generated. They present results for the pressure on the top of the airframe at
various advance ratios. The accuracy of these measurements was about 10 mm,
or one core diameter; recent experiments by Peterson et al (1995) have achieved
sub-millimeter accuracy using digital imaging. Leishman and Bi (1990) also
report instantaneous surface pressure measurements taken at various advance
ratios and note that these unsteady pressure peaks are very large and can swamp
the steady mean pressure. They focus on the time history of the pressure at
a fixed point on the airframe. Bagai and Leishman (1992b) found that the
presence of the fuselage affects the path of the tip-vortex significantly.
Based on the work of Liou et al (1990), Brand et al (1989, 1990), and
Leishman and Bi (1990), a reasonable picture of the dominant features of the
unsteady pressure on the airframe under direct wake impingement conditions
has emerged. For the given rotor-body separation distance, on the top of the
airframe the blade passage effect generates a large positive pressure rise at
ψ = 0 (and ψ = 2π N
for an N -bladed rotor), followed by a very sharp suction
peak of the same order if the advance ratio is low enough so that the tip-vortex
“collides” with the airframe. For a vortex whose age is 180◦ + ψ (i.e. for a two-
bladed rotor) at the given rotor-body separation distance of Brand et al (1990),
the blade passage effect is dominant up to about ψ = 30◦ with the influence of
the tip-vortex being dominant for the next ∼ 30◦ . Depending on the advance
ratio, the maximum amplitude of the suction peak occurs near ψ = 50◦ ; at this
time, the vortex age is 230◦ .
Despite the fact that the modeling efforts described so far are relatively
sophisticated, all of them demonstrate limited utility in reproducing unsteady
fuselage airloads. In an attempt to better understand the important features of
the pressure distribution on the top of the airframe under vortex-wake collision
conditions, Affes et al (1993a) use a Biot-Savart representation of the tip-vortex
to calculate the tip-vortex path in the time frame from ψ = 0◦ to ψ = 60◦ for two
different advance ratios. At this point in time for the two-bladed rotor studied,
the vortex age is 180◦ + ψ. A typical result for the vortex path for ψ > 0
and the associated pressure, calculated from the unsteady Bernoulli equation
on the top of the airframe, is depicted on Figure 18 (Affes et al 1993a). Note
the good agreement with experiment; no adjustable constants are used in the
analysis and only a crude model of the inboard vortex sheet is used. However,
the computational results begin to deviate from the experimental results for the
pressure distribution at about ψ = 48◦ and the experimental suction peak is
gone by ψ = 60◦ , a time scale of about 0.5 msec. The reason for this is believed
to be the complete transformation of the flow in the vortex core from a high
swirl velocity region characteristic of a “vortex” to a low swirl velocity region
under the action of the axial velocity in the core of the vortex.
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
Additional experimental results (Kim and Komerath 1995) suggest that the
suction peak in the pressure is convected around the retreating side of the
airframe for a long period of time after it has disappeared from the top of the
airframe. This effect is seen in plots of the pressure contours around the airframe
from Kim and Komerath (1995) as shown in Figure 19a at a rotor phase angle
of ψ = 120◦ . Note that on the retreating side (φ < 0; φ refers to the angle
measured from the top of the airframe and not the zero-lift angle of attack as
in Figure 1) the suction peak is still strong while on the advancing side the
pressure on the airframe at a point coinciding with the vortex core is actually
positive. Lee et al (1995) has demonstrated these effects for a model problem;
a summary is sketched in the accompanying Figure 19b. These large-scale
suction peak effects on the retreating side could be important in maneuvering
(a)
(b)
Figure 18 Comparison between computation and experiment of the vortex trajectory and pressure
distribution on the top of a model airframe for advance ratio µ = 0.1. The z-axis is in the direction
of forward flight. (a) Tip-vortex trajectory ψ = 0, 30, 60◦ ; A denotes advancing side and R denotes
retreating side. (b) Pressure at ψ = 42◦ on the left and ψ = 48◦ on the right. From Affes et al
(1993a).
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
558 CONLISK
(a)
Stagnation Suction
Axial Flow
Axial Flow
Γ Γ
Pressure
Figure 19 (a) Pressure distribution around a model airframe as measured by Kim and Komerath
(1994), showing the differences exhibited on the advancing and retreating sides of the airframe.
Arrows denote region below the vortex core. The angle φ measures distance from the top of the
cylinder and φ > 0 denotes the advancing blade side. XRb measures distance along the airframe in
the forward flight direction. (b) A sketch of the tip-vortex structure along the sides of the airframe
as described by Kim and Komerath (1995) and Lee et al (1995). The vortex core radius is greatly
enlarged for clarity.
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
flight and could be a factor in the fatigue life of the airframe. A similar effect
is seen in the computations of Marshall (1994).
Steady flow Navier-Stokes simulations of the full wake-airframe interaction
in forward flight have begun to be analyzed by Zori and Rajagopalan (1995).
The rotor is treated as a source term in the governing equations; the magni-
tudes of the sources are determined from the blade geometry and aerodynamic
characteristics. The detailed boundary layer behavior under the tip-vortex as it
approaches the airframe is discussed by Affes et al (1993b). It is shown that
flow reversal occurs in the form of an eddy which grows in time. Eventually
the boundary layer fluid will wrap around the vortex fluid and this may alter its
effective strength.
The sophistication of large-scale computational techniques has improved to
the point that entire helicopter configurations may be incorporated. Duque and
Dimanlig (1994) and Duque et al (1995) have produced results for the RAH-66
Commanche helicopter. Figure 20 depicts the streakline patterns for the V-22
Tiltrotor operating in the helicopter mode (from Meakin 1993). The V-22 is
a dual purpose machine that can operate in both fixed-wing and rotary-wing
modes. The computation of Meakin (1993) uses an overset grid approach in
which a family of grid systems, each designed for accuracy in a particular
region of the flow, are employed. Data is communicated getween grids by
interpolation. Overset methods are generally viewed as being more accurate
Figure 20 Streaklines for the V-22 Tiltrotor operating in the helicopter mode. From Meakin
(1993).
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
560 CONLISK
than conventional methods at the cost of slightly higher CPU time (Bangalore
and Sankar 1996).
7. SUMMARY
Modern helicopter aerodynamics is challenging because the flow field gener-
ated by a helicopter is extremely complicated and difficult to measure, model,
and predict; moreover, experiments are expensive and also difficult to conduct.
Helicopter aerodynamics encompasses a variety of complicated and inherently
nonlinear flow phenomena including dynamic stall, blade-vortex interaction,
and shock interactions. Fluid-structure interaction problems occur throughout
the rotor cycle because the rotor blades undergo pitching, flapping, and aeroe-
lastic motions dynamically in all regimes of operation. The nature of the rotor
wake depends on a number of factors including the dimensions and shape of
the blades, the number of blades, tip-speed, and regime of operation. Because
of the complexity of the problem, current models and experiments involve a
number of simplifications in order to focus on a specific aspect of the problem.
In spite of this, the past ten years have seen great advances in the under-
standing of isolated rotor wakes. With the advance of computational resources
in recent years and the development of sophisticated experimental techniques
designed to measure the fully three-dimensional character of rotor wakes, rea-
sonable comparisons with experiment may now be obtained through computa-
tional studies at low to moderate Mach number for an isolated rotor in hover at
model scale, provided aeroelastic deformations are negligible. The agreement
between computations and experiment for an isolated rotor in forward flight is
not as good, but advances in this regime are being made at a rapid pace. In addi-
tion, drag and moment are difficult to predict especially under stall conditions,
and there is little experimental data to guide computational approaches.
At the present time, rotor design is accomplished by standard momentum the-
ory for initial estimates of thrust and moment. These estimates are often refined
using vortex methods in conjunction with lifting-line, lifting-surface, or panel
methods to calculate the strength of the wake. Vortex methods are generally
the least computationally demanding of the modeling approaches. CFD tech-
niques, in the form of Navier-Stokes computations, while being substantially
predictive in nature and based on first principles (apart from the specification
of a turbulence model), are at the present time too computationally demanding
and not accurate enough to be used extensively in design. To remedy this,
methods must be sought to improve the speed and computational affordabil-
ity of these programs which can only be run on a supercomputer. The use
of higher-order discretization methods along with adaptive grid schemes can
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
reduce the computational burden and increase accuracy. It is possible that part
of this improvement can be accomplished by the natural evolution of the com-
putational hardware which is continually coming down in price and increasing
in computational speed; however, it is probable that this natural evolution will
not occur fast enough, and future computations will likely benefit from the use
of parallel processing.
Advances in experimental techniques in the past ten years have been sub-
stantial and will continue. These advances will likely help modeling efforts.
Detailed three-dimensional and unsteady measurements of the rotor wake flow
are now beginning to be taken on a routine basis. However, the amount of data
to be acquired and process in this activity is staggering, and efforts must be
taken to process the data in a more timely and efficient manner. This process
is aided by the availablility of post-processing tools, normally designed for
computational studies, which are now being used to reduce experimental data.
The long-range goal of any modeling effort is to be able to compute the
entire unsteady flow field around a helicopter operating in any flight regime.
The objective is to reduce the amount of expensive full-scale testing which
is now required. However, significant barriers to this goal will exist for the
forseeable future. Rotor loads even at model scale in the forward flight and
descent regimes cannot now be calculated. An additional limiting technological
barrier is the issue of noise which is very acute at the higher flight speeds of
modern helicopters. Noise is identified as the major impediment to increased
commercial use of rotorcraft. In order to design quieter helicopters, accurate
computation of BVI and high speed impulsive noise must be made; inevitably
this requires a more accurate computation of the flow field itself.
The extrapolation of model-scale results to the full scale is only partially
successful: Considerable differences between full-scale flight tests and full-
scale wind tunnel tests are common (Tung et al 1995). At the current time,
aeroelastic effects which are significant at full scale, are most often neglected
in the model-scale computations.
Despite these challenges, great advances have been made in modeling the
helicopter rotor wake and the flow around a helicopter during the last ten to
fifteen years. Moreover, with the advances in both experimental and computa-
tional technology, progress is likely to continue at an even greater rate during
the next ten years.
ACKNOWLEDGMENTS
The author is grateful to those who provided references and other information
for this article. Several colleagues have made substantial contributions to the
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
562 CONLISK
Literature Cited
Affes H, Conlisk AT, Kim JM, Komerath NM. ence of a fully interactive wake. Presented at
1993a. Model for rotor tip-vortex-airframe Annu. Forum Am. Hel. Soc., 46th, Washing-
interaction, part 2: comparison with experi- ton D.C.
ment. AIAA J. 31(1)2):2274–82 Betz A. 1915. Die wichtigsten Grundlagen für
Affes H, Xiao Z, Conlisk AT, Kim JM, den Entwurf von Luftschrauben. Z. Mathe-
Komerath NM. 1993b. The three-dimensi- matik 6:97–103.
onal boundary layer flow due to a rotor tip Bliss DB, Teske ME, Quackenbush TR. 1987. A
vortex. Presented at AIAA Fluid Dyn. Conf., new methodology for free wake analysis us-
24th, Orlando. ing curved vortex elements. NASA CR 3958.
Ahmad J, Duque EPN. 1996. Helicopter rotor Bliss DB, Quackenbush TR, Bilanin AJ. 1983.
blade computation in unsteady flows using A new methodology for helicopter free wake
moving embedded grids. J. Aircr. 33:54–60. analyses. Proc. Forum Am. Hel. Soc., 39th.
Bagai A, Leishman JG. 1992a. Experimen- Bliss DB, Miller WO. 1990. Efficient free
tal study of rotor wake/body interactions in wake calculations using analytical numerical
hover. J. Am. Hel. Soc. 37(4):48–57. matching. J. Am. Hel. Soc. 38(2):43–52.
Bagai A, Leishman JG. 1992b. Improved wide- Bliss DB, Dadone LU, Wachspress DA. 1987.
field shadowgraph set-up for rotor wake vi- Computations of rotor wake modeling for
sualization. J. Am. Hel. Soc. 37(3):86–92. high speed applications. Proc. Annu. Forum
Bangalore A, Sankar LN. 1996. Forward flight Am. Hel. Soc., 43rd
analysis of slatted rotors using Navier-Stokes Bramwell ARS. 1976. Helicopter Dynamics.
methods. AIAA paper 96-0675. London: Edward Arnold. 408 pp.
Batchelor GK. 1967. Introduction to Fluid Dy- Brand AG, McMahon HM, Komerath NM.
namics. Cambridge: Cambridge Univ. Press. 1989. Surface pressure measurements on a
615 pp. body subject to vortex wake interaction. AIAA
Beaumier P. 1994. A coupling procedure be- J. 27(5):569–74.
tween a rotor dynamic code and a 3D un- Brand AG, McMahon HM, Komerath NM.
steady full potential code. Presented at Am. 1990. Correlations of rotor wake airframe in-
Hel. Soc. Aeromech. Spec. Conf., San Fran- teraction measurements and flow visualiza-
cisco. tion data. J. Am. Hel. Soc. 35(4):4–15.
Beam R, Warming RF. 1978. An im- Brand AG. 1989. An experimental investigation
plicit factored scheme for the compressible of the interaction between a model rotor and
Navier-Stokes equations. AIAA J. 16(4):393– airframe in forward flight. PhD thesis. Geor-
402. gia Inst. Technol.
Berry JD, Althoff SL. 1990. Inflow velocity per- Bridgeman JO, Steger JL, Caradonna FX. 1982.
turbations due to fuselage effects in the pres- A conservative finite-difference algorithm for
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
the unsteady transonic potential equation in USAAMRDL Tech. Rep. 77-4, Eustis Direc-
generalized coordinates. AIAA paper 82- torate, U.S Army Air Mobility Research and
1388. Development Laboratory, Fort Eustis, VA.
Caradonna FX. 1992. The application of CFD Egolf TA. 1988. Helicopter free-wake pre-
to rotary wing aircraft. NASA TM 102803. diction of complex wake structures under
Caradonna FX, Tung C. 1981. Experimental and blade-vortex interaction operating condi-
analytical studies of a model helicopter rotor tions. Proc. Annu. Forum Am. Hel. Soc.,
in hover. NASA TM 81232. 44th, Washington DC, pp. 819–32.
Caradonna FX, Strawn RC, Bridgeman JO. Egolf TA, Landgrebe AJ. 1983. Helicopter rotor
1988. An experimental and computational wake geometry and its influence in forward
study of rotor-vortex interactions. Proc. Eur. flight. Vol. 1: Generalized wake geometry
Rotorcraft Forum, 14th, Milano. and wake effect on rotor airloads and perfor-
Carr LW. 1988. Progress in analysis and predic- mance. NASA CR 3726.
tion of dynamic stall. J. Aircr. 25(1):6–17. Egolf TA, Sparks SP. 1986. A full potential rotor
Carr LW, Chandrasekhara MS. 1995. An assess- analysis with wake influence using an inner-
ment of the impact of compressibility on dy- outer domain technique. Proc. Annu. Forum
namic stall. AIAA paper 95-0779. Presented Am. Hel. Soc., 42nd, pp. 997–1011.
at Aerosp. Sci., 33rd, Reno, Nevada. Elliott JW, Althoff SL. 1988. Inflow measure-
Castles W, Gray RB. 1951. Empirical relation ment made with a laser velocimeter on a heli-
between induced velocity, thrust, and rate of copter model in forward flight. Volumes. I-V,
descent of a helicopter rotor as determined NASA TM 100541–100545
by wind tunnel tests on four model rotors. Felker FF, Light JS. 1988. Aerodynamic inter-
NACA TN-2474. actions between a rotor and wing in hover. J.
Chang I-C, Tung C. 1985. Numerical solu- Am. Hel. Soc. 33(2):53–61.
tion of the full-potential equation for rotors Funk R, Fawcett PA, Komerath NM. 1993. In-
and oblique wings using a new wake model. stantaneous velocity fields in a rotor wake by
AIAA paper 85-0268. spatial correlation velocimetry. AIAA Paper
Chen RTN. 1990. A survey of nonuniform 93-3081, Fl. Dyn. Conf., 24th, Orlando, Fla.
inflow models for rotorcraft flight dynam- Gallman JM et al. 1995. Effect of wake struc-
ics and control applications. Vertica 14:147– ture on blade-vortex interaction phenomena:
90. acoustic prediction and validation. Presented
Crouse GL, Leishman JG. 1992. Interactional at CEASAIAA Aeroacoustic Conf., 1st, Mu-
aerodynamic effects on rotor performance in nich.
hover and forward flight. Presented at Annu. Gennaretti M, Morino L. 1992. A boundary ele-
Forum Am. Hel. Soc., 48th, Washington D.C. ment method for the potential, compressible
Crouse GL, Leishman JG. 1993. A new method aerodynamics of bodies in arbitrary. Aero. J.
for improved free-wake convergence in hover Roy. Aero. Soc.
and low-speed forward flight. Presented at Gessow A. 1986. Understanding and predicting
Aerosp. Sci. Meet., 21st. helicopter behavior—then and now. J. Am.
Desopper A, Lafon P, Ceroni P, Phillippe JJ. Hel. Soc. 31(1):3–28.
1986. Ten years of rotor flow studies at ON- Gessow A, Myers G. 1952. Aerodynamics of
ERA: state of the art and future studies. Proc. Helicopter. New York: Frederick Ungar
Annu. Forum Am. Hel. Soc., 42nd, Washing- Ghee T, Berry JD, Zori LAJ, Elliott JW. 1995.
ton, D.C. Wake geometry measurements and analyti-
Duque EPN, Dimanlig ACB. 1994. Navier- cal calculations on a small-scale rotor model.
Stokes simulation of the AH-66 (Com- NASA TP 3584.
manche) helicopter. Presented at Am. Hel. Glauert H. 1922. An aerodynamic theory of the
Soc. Aeromech. Spec. Conf., San Francisco. airscrew. British ARC R & M 786.
Duque EPN, Berry JD, Dimanlig ACB, Budge Glauert H. 1928. On horizontal flight of a heli-
AM. 1995. A comparison of computed and copter. R & M 1157.
experimental flowfields of the RAH-66 he- Gorton SA, Poling DR, Dadone L. 1995a.
licopter. Presented at Am. Hel. Soc. Spec. Laser velocimetry and blade pressure mea-
Meet. on Aeromech. Tech. Prod. Des., surements of a blade-vortex interaction. J.
Bridgeport, Conn. Am. Hel. Soc. 40(2):15–23.
Duque EPN, Biswas R, Strawn RC. 1995. A so- Gorton S, Meyers J, Berry JD. 1995b. Veloc-
lution adaptive structured/unstructured over- ity measurements near the empennage of a
set grid solver with applications to helicopter small-scale helicopter model. Proc. Annu.
rotor flows. AIAA paper 95-1766. Presented Forum Am. Hel. Soc., 51st.
at Appl. Aerodyn. Conf., 13th, San Diego, Gray RB. 1955. On the motion of the helical
Calif. vortex shed from a single-bladed hovering
Dvorak FA, Maskew B, Woodward FA. 1977. helicopter rotor and its application to the cal-
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
564 CONLISK
culation of the spanwise aerodynamic load- Kim JM, Komerath NM, Liou SG. 1994. Vor-
ing, Princeton Univ. Aero. Engr. Dept., Re- ticity concentration at the edge of the inboard
port 313. vortex sheet. J. Am Hel. Soc. 39(2):30–4.
Gray RB. 1956. An aerodynamic analysis of Kittleson JK, Yu Y. 1985. Transonic rotor flow
a single-bladed rotor in hovering and low- measurement technique using holographic
speed forward flight as determined from interferometry. J. Am. Hel. Soc. 30(4):3–10.
smoke studies of the vorticity distribution in Klemin A. 1945. Principles of rotary-wing air-
the wake. Princeton Univ. Aero. Engr. Dept., craft. AERO Digest, May 1 and June 1, 1945.
Report 356. Kocurek JD, Tangler JL. 1976. A prescribed
Gray RB. 1992. Vortex modelling for rotor aero- wake lifting surface hover performance anal-
dynamics—the 1991 Alexander A. Nikolsky ysis. Proceedings of the 32nd Forum of the
Lecture. J. Am. Hel. Soc. 37(1):3–14. Am. Hel. Soc., May. Preprint 1001.
Ham ND, Garelick MS. 1968. Dynamic stall Landgrebe AJ. 1972. The wake geometry of a
considerations in helicopter rotors. J. Am. hovering rotor and its influence on rotor per-
Hel. Soc. 13(2):49–55. formance. J. Am. Hel. Soc. 17(4):3–15.
Hardin JC, Lamkin SL 1986. An Euler code Landgrebe AJ, Johnson BV. 1974. Measure-
calculation of blade-vortex interaction noise. ments of model helicopter rotor flow veloci-
Winter Annual Meeting of the American So- ties with a laser Doppler velocimeter. J. Am.
ciety of Mechanical Engineers. Paper 86- Hel. Soc. 19(3):39–43.
WANCA-3. Landgrebe AJ. 1994. New directions in rotor-
Hassan AA, Tung C, Sankar, LN. 1992. Eu- craft computational aerodynamics research
ler solutions for self-generated rotor blade- in the U. S. 75th AGARD Fluid Dynamics
vortex interactions. Int. J. Num. Meth. Fluids. Panel Meeting on Aerodynamics and Aeroa-
15:427–51. coustics of Rotorcraft, Berlin.
Heller H, Schultz K-J , Ahmed SR. 1994. Un- Landgrebe AJ, Moffitt RC, Clark DR. 1977.
steady surface pressure characteristics on he- Aerodynamic technology for advanced rotor-
licopter blades: a key to the physics of rotor craft. J. Am. Hel. Soc. 22 (2, 3)
noise. 19th ICAS Congress, paper 94-2.7.2, Lee J, Xiao Z, Burggraf OR, Conlisk AT,
Anaheim, CA. Komerath NM. 1995. An inviscid analysis
Hess JL, Smith AMO. 1972. Calculation of of vortex-surface collisions. AIAA paper 95-
potential flow about arbitrary bodies. In 2237.
Progress in the Aeronautical Sciences, Vol. Leishman JG, Bi N. 1990. Aerodynamic inter-
8. New York:Pergamon. actions between a rotor and a fuselage in for-
Hess JL. 1990. Panel methods in computational ward flight. J. Am. Hel. Soc. 35(3):22–31.
fluid dynamics. Annu. Rev. Fluid Dynamics. Leishman JG, Baker A, Coyne A. 1995. Mea-
Hoad DR. 1990. Rotor induced-inflow-ratio surements of rotor tip vortices using three-
measurements and CAMRAD calculations. component laser Doppler velocimetry. Am.
NASA TP 2946, January 1990. Hel. Soc. Specialists’ Meeting on Aerome-
Hoad DR, Althoff SL, Elliott JW. 1988. Ro- chanics Technology and Product Design,
tor inflow variability with advance ratio. Pro- Bridgeport, Conn.
ceedings of the 44th Annual Forum of the Leishman JG, Beddoes TS. 1989. A semi-
Am. Hel. Soc., Washington DC, June 16-18, empirical model for dynamic stall. J. Am. Hel.
1988. pp. 57–72. Soc. 34(3):3–17.
Jameson A, Caughey DA. 1977. Numerical cal- Liou SG, Komerath NM, McMahon HM. 1989.
culation of the transonic flow past a swept Velocity measurements of airframe effects on
wing. ERDA Research and Development, a rotor in low-speed forward flight. J. Aircraft
Mathematics and Computing, June 1977. 26(4):340–8.
Johnson W. 1971. A lifting-surface solution for Liou SG, Komerath NM, Lal MK. 1994. Mea-
vortex-induced airloads. AIAA J. 9(4):689– surement around a rotor blade excited in
695. pitch. Part 1: Dynamic inflow. J. Am. Hel.
Johnson W. 1980. Helicopter Theory. Prince- Soc. 39(2):3–12.
ton, NJ: Princeton Univ. Press. Liou SG, Komerath NM, McMahon HM. 1990.
Johnson W. 1986. Recent developments in Measurement of the interaction between a
rotary-wing aerodynamic theory. AIAA J. rotor tip-vortex and a cylinder, AIAA J.
24(8):1219–44. 28(6):975–981.
Katz J, Plotkin A. 1991. Low-Speed Aerody- Lock 1924. Experiments to verify the indepen-
namics from Wing Theory to Panel Methods. dence of the elements of an airscrew blade.
New York: McGraw Hill. British R & M 853.
Kim JM, Komerath NM. 1995. Summary of the Loewy, RG. 1957. A two-dimensional approx-
interaction of a rotor wake with a circular imation to the unsteady aerodynamics of ro-
cylinder. AIAA J. 33(3):470–8. tary wings. J. Aero. Sci. 24(2):81–92.
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
Lorber PF, Egolf TA. 1990. An unsteady rotor- Morino L, Gennaretti M. 1991. Boundary in-
fuselage interaction analysis. J. Am. Hel. Soc. tegral equation methods for aerodynamics.
35(7):32–42. In Computational Nonlinear Mechanics in
Lorber PF. 1991. Aerodynamic results of a Aerospace Engineering, ed. S. N. Atluri.
pressure-instrumented model rotor test at the Prog. Astronaut. Aeronaut. 147:
DNW. J. Am. Hel. Soc., 36(4):66–76. Norman TR, Light JS. 1987. Rotor tip-vortex
Lorber PF. 1993. Tip vortex, stall vortex, and geometry measurements using the wide-field
separation observations on pitching three- shadowgraph technique. J. Am. Hel. Soc.
dimensional wings. AIAA 24th Fluid Dy- 32(2):40–50.
namics Conference, Orlando, Fla. Norman TR, Yamauchi GK. 1991. Full-scale
Lorber PF, Carta FO. 1988. Airfoil dynamic stall investigation of aerodynamic interactions be-
at constant pitch rate and high Reynolds num- tween a rotor and a fuselage. Proceedings of
ber. J. Aircraft 25(6):548–556. the 47th Annual Forum of the Am. Hel. Soc.,
Lorber PF, Stauter RC, Haas RJ, Anderson TJ, pp. 461-86.
Torok MS, Kohlhepp FW. 1994. Techniques Parthasarathy SP, Cho YI, Black LH. 1985.
for comprehensive measurement of model Wide-field shadowgraph flow visualization
helicopter rotor aerodynamics. Proceedings of tip vortices generated by a helicopter rotor.
of the 50th Forum of the Am. Hel. Soc., May AIAA 85-1557, July 1985.
11–14, 1994, pp. 797–814. Peterson K, Frank RB, Mahalingam R,
Mangler KW, Squire HB. 1950. The induced Komerath NK, Conlisk AT. 1995. Recent ex-
velocity field of a rotor. R & M 2642. periments on vortex collision with a cylinder.
Maskew B. 1986. VSAERO theory document. AIAA 95-2236.
NASA CR-4023, November, 1986. Phillipe JJ, Roesch P, Dequin AM, Cler A. 1985.
Marshall JS. 1994. Vortex cutting by a blade. A survey of recent developments in helicopter
Part II: Computations of vortex response. aerodynamics. AGARD-LS-139. Helicopter
AIAA J. 32(7):1428–36. Aeromechanics, April 1985, pp. 2-1–2-40.
Mba MN, Meylan C, Mansca C, Favier D. 1984. Prouty RW, Amer KB. 1982. The YAH-64 em-
Radial distribution of circulation of a rotor in pennage and tail rotor—a technical history.
hover measured by a laser velocimeter. Pro- Proceedings of the 38th Annual Forum of the
ceedings of the Tenth European Rotorcraft Am. Hel. Soc., May 4–7, 1982, pp. 247–61.
Forum, The Hague, Netherlands. Quackenbush TR, Bliss DB, Wachpress, DA.
McAlister KW, Schuler CA, Branum L, Wu JC. 1989. Free wake analysis of rotor configura-
1995. 3-D measurements near a hovering ro- tions for reduced vibratory airloads. Proceed-
tor for determining profile and induced drag. ings of the AHS National Specialists’ Meet-
NASA TP 3577, August 1995. ing on Rotorcraft Dynamics, November.
McCormick BW. 1967. Aerodynamics of Quackenbush TR, Lam C-MG, Bliss DB. 1994.
V/STOL Flight. New York:Academic Press. Vortex methods for the computational analy-
McCroskey WJ. 1995. Vortex wakes of rotor- sis of rotor/body interaction. J. Am. Hel. Soc.
craft. AIAA paper 95-0530. 39(4):14–24.
McCroskey WJ. 1981. The phenomenon of dy- Ramachandran K, Schlechtriem S, Caradonna
namic stall. NASA TM 81264, March. FX, Steinhoff, J. 1993. The application of
Meakin RL. 1993. Moving body overset grid vorticity embedding to the computation of ad-
methods for complete tiltrotor simulations, vancing rotor flows. Presented at Annu. Fo-
AIAA paper 93-3350-CP. Presented at the rum Am. Hel. Soc., 49th, May 19–21, 1993,
11th AIAA Computational Fluid Dynamics pp. 571–84.
Conference, Orlando, FL, July 1993. Ramachandran K, Tung C, Caradonna FX.
Miller WO. 1993. A fast adaptive resolution 1989. Rotor hover performance prediction
method for efficient free wake calculations. using a free-wake, computational fluid dy-
Proceedings of the 49th Forum of the Am. namics method. J. Aircraft 26:1105-10.
Hel. Soc., May 19-21, 1993, pp. 617–30. Rankine, WJ. 1865. On the mechanical princi-
Miller WO, Bliss DB. 1993. Direct periodic so- ples of the action of ship propellers. Trans.
lutions of rotor free wake calculations. J. Am. Inst. Naval Arch. 6:13–39.
Hel. Soc. 38(2):53–60. First presented at the Reichert G. 1985. Helicopter aeromechanics—
Proceedings of the 46th Forum of the Am. introduction and historical review. AGARD-
Hel. Soc., 1990. LS-139. Helicopter Aeromechanics, April
Moore DW. 1972. Finite amplitude waves on 1985, pp. 1-1–1-23.
aircraft trailing vortices. Aero. Quart. 23:307- Sarpkaya, T. 1989. Computational methods
14. with vortices. The 1988 Freeman Scholar lec-
Morino L, Kaprielian Z, Sipcic SR. 1985. Free ture. J. Fluids Engineering 111:5–52.
wake analysis of helicopter rotors. J. Aircraft Scully MP. 1967. A method of computing heli-
9(2):127–40. copter rotor wake distortion. Massachusetts
December 3, 1996 17:28 Annual Reviews AR023-15n AR23-15
566 CONLISK
sively parallel computer. AIAA J. 29(1):58– uct Design, Bridgeport CT, October 11–13,
67. 1995.
Wilson JC, Mineck RE. 1975. Wind-tunnel in- Yu YH, Lee S, McAlister KW, Tung C,
vestigation of helicopter rotor wake effects Wang CM. 1995. Dynamic stall control
on three helicopter fuselage models. NASA for advanced rotorcraft application. AIAA J.
TM X-3185, March 1975. 33(2):289–95.
Yu YH. 1995. Rotor blade vortex interaction Zori LAJ, Rajagopalan RG. 1995. Navier-
noise: generating mechanisms and its control Stokes calculations of rotor-airframe inter-
concepts. Am. Hel. Soc. Specialists’ Meet- action in forward flight. J. Am. Hel. Soc.
ing on Aeromechanics Technology and Prod- 40(2):57–67.