JOURNAL OF PROPULSION AND POWER

Vol. 10, No. 5, Sept.-Oct. 1994

Design of Optimum Propellers

Charles N. Adkins*
Falls Church, Virginia 22042
and
Robert H. Liebeckt
Douglas Aircraft Company, Long Beach, California 90846

Improvements have been made in the equations and computational procedures for design of propellers and
wind turbines of maximum efficiency. These eliminate the small angle approximation and some of the light
loading approximations prevalent in the classical design theory. An iterative scheme is introduced for accurate
calculation of the vortex displacement velocity and the flow angle distribution. Momentum losses due to radial
Downloaded by ORTA DOGU TEKNIK UNIVERSITESI on March 10, 2019 | http://arc.aiaa.org | DOI: 10.2514/3.23779

flow can be estimated by either the Prandtl or Goldstein momentum loss function. The methods presented here
bring into exact agreement the procedure for design and analysis. Furthermore, the exactness of this agreement
makes possible an empirical verification of the Betz condition that a constant-displacement velocity across the
wake provides a design of maximum propeller efficiency. A comparison with experimental results is also
presented.

Nomenclature £ = displacement velocity ratio, v'lV

a = axial interference factor 17 = propeller efficiency
a' = rotational interference factor A = speed ratio, V/flR
B = number of blades g = nondimensional radius, r/R = AJC
b = axial slipstream factor \e = nondimensional Prandtl radius
Cd = blade section drag coefficient £0 = nondimensional hub radius
C, = blade section lift coefficient p = fluid density
Cp - power coefficient, P/pn3D5 a = local solidity, Bc/2/nr
CT = thrust coefficient, T/pn2D4 <f> — flow angle
Cx = torque force coefficient <f)t = flow angle at the tip
Cy = thrust force coefficient fl = propeller angular velocity
c = blade section chord
D = propeller diameter, 2R Superscript
D' = drag force per unit radius ' = derivative with respect to r or £, unless otherwise
F = Prandtl momentum loss factor noted
G = circulation function
/ = advance ratio, VlnD Introduction
K = Goldstein momentum loss factor
L'
n
=
=
lift force per unit radius
propeller rps
I N 1936, a classic treatise on propeller theory was authored
by H. Glauert.1 In this work, a combination of momentum
theory and blade element theory, when corrected for mo-
P = power into propeller mentum loss due to radial flow, provides a good method for
Pc = power coefficient, 2P/pV37rR2 analysis of arbitrary designs even though contraction of the
Q = torque propeller wake is neglected. Although the theory is developed
R = propeller tip radius for low disc loading (small thrust or power per unit disc area),
r = radial coordinate it works quite well for moderate loading, and in light of its
T = thrust simplicity, is adequate for estimating performance even for
Tc = thrust coefficient, 2T/pV27rR2 high disc loadings. The conditions under which a design would
V = freestream velocity have minimum energy loss were stated by A. Betz2 as early
v' = vortex displacement velocity as 1919; however, no organized procedure for producing such
W = local total velocity a design is evident in Glauert's work. Those equations which
wn = velocity normal to the vortex sheet are given by Betz make extensive use of small-angle approx-
wt = tangential (swirl) velocity imations and relations applicable only to light loading con-
x = nondimensional distance, £lrlV ditions. Theodorsen3 showed that the Betz condition for min-
a = angle of attack imum energy loss can be applied to heavy loading as well.
j8 = blade twist angle In 1979, E. Larrabee4 resurrected the design equations and
F = circulation presented a straightforward procedure for optimum design.
s = drag-to-lift ratio However, there are still some problems: first, small angle
approximations are used; second, the solution for the dis-
placement velocity is accurate only for vanishingly small val-
Presented as Paper 83-0190 at the AIAA 21st Aerospace Sciences ues (light loading), although an approximate correction is
Meeting, Reno, NV, Jan. 10-13, 1983; received July 23, 1992; re-
vision received July 6, 1993; accepted for publication Dec. 15, 1993. suggested for moderate loading; and third, there are viscous
Copyright © 1994 by the American Institute of Aeronautics and terms missing in the expressions for the induced velocities.
Astronautics, Inc. All rights reserved. These viscous terms must be included in the design equations
*Consulting Engineer. Member AIAA. if they are to be consistent with the classical propeller analysis.
tMcDonnell Douglas Fellow. Fellow AIAA. This approach is given later.
676
 ADKINS AND LIEBECK: DESIGN OF OPTIMUM PROPELLERS 677

The purpose of this article is to correct these difficulties VORTEX FILAMENT

and bring the design method into exact agreement with the AFTER TIME
INCREMENT, At
analysis. It is then possible to verify empirically the optimality
of the design. This work was initiated at McDonnell Douglas
in 1980 in response to a requirement for simple estimates of
propeller performance. In-house methods, if they existed, had AXIS
been irretrievably archived. An early version was presented
as AIAA Paper 83-0190. Continuous requests for copies of
the paper plus some refinements to the method have moti-
vated its publication in the Journal.
VORTEX FILAMENT (t = 0)
Momentum Equations
Detailed axial and general momentum theory is described Fig. 2 Definition of displacement velocity v' in the propeller wake.
by Glauert,1 and only a brief summary is given here to em-
phasize several important features. Consider a fluid element
of mass dm, far upstream moving toward the propeller disc develops the Betz condition for heavy loading by including
in a thin, annular stream tube with velocity V. It arrives at the contraction of the wake. He shows that sufficiently far
the disc with increased velocity, V(l + 0), where a is the downstream in the contracted wake, the vortex sheet must
Downloaded by ORTA DOGU TEKNIK UNIVERSITESI on March 10, 2019 | http://arc.aiaa.org | DOI: 10.2514/3.23779

axial interference factor. At the disc, dm exists in the annulus be the same regular screw surface for a propeller of minimum
27rr dr, and the mass rate per unit radius passing through the induced energy loss. This optimum vortex sheet acts as an
disc is 2irrpV(l + a), neglecting radial flow. The element dm Archimedean screw, pumping fluid aft between rigid spiral
moves downstream into the far wake, increasing speed to the surfaces.
value V(l 4- &), where b is the axial slipstream factor. Axial At the blade station, r, the total lift per unit radius is given
momentum theory determines b to be exactly 20, whereas the by
general theory (which includes rotation of the flow) deter-
mines b to be approximately 2a. Using the axial approxi-
mation, which is generally accepted, the overall change in L> = f = BpWT (2)
momentum of the element is 2VaF dm where F, the momen-
tum loss factor, accounts for radial flow of the fluid. The and in the wake, the circulation in the corresponding annulus
thrust per unit radius T', acting on the annulus can now be is
expressed as
i rri BY = 27rrFwf (3)
T = — = 27rrpV(l + a)(2VaF) (la)
Setting the circulation F in Eq. (2) equal to that in Eq. (3)
will ultimately determine that circulation distribution F(r) that
By similar arguments, the torque per unit radius Q' is given
minimizes the induced power of the propeller.
by In order to obtain F(r), it is necessary to relate vvr to a more
measurable quantity. Figure 2 shows the wake vortex filament
Q'lr = 2irrpV(l + a)(2Slra'F) (Ib) at station r and the definition of the various velocity com-
ponents there. The motion of the fluid must be normal to the
Flow geometry about a blade element at the disc is shown in local vortex sheet, and this normal velocity is wn. Therefore,
Fig. 1, where W acts on the blade element with a, and acts the tangential velocity is given by
on the disc at <j>. F goes from about 1 at the hub (where the
radial flow is typically negligible), to 0 at the tip, and is not
unlike the spanwise loading of a wing. The functional form Wt = wn sin </>
of this factor was first estimated by Prandtl1-2 and a more
accurate, though more complex, form was determined by However, for a coordinate system fixed to the propeller disc,
Goldstein5 and Lock.6"8 the axial velocity of the vortex filament would be

Circulation Equations v' = H^/COS (/>

At each radial position along the blade, infinitesimal vor-
tices are shed and move aft as a helicoidal vortex sheet. Since where the increase in magnitude of v' over wn is due to ro-
these vortices follow the direction of local flow, the helix angle tation of the filament. This is analogous to a barber pole where
of the spiral surface is ((>, shown in Fig. 1. The Betz condition it appears that the stripes are translating in spite of the fact
for minimum energy loss, neglecting contraction of the wake, that only a rotational velocity exists. It will become clear that
requires the vortex sheet to be a regular screw surface; i.e., it is convenient to use v', and the corresponding displacement
rtan <t> must be a constant independent of radius. Theodorsen3 velocity ratio, £ = v'lV. The tangential velocity is then

Wt = V£ sin cf> cos 4>

and the circulation of Eq. (3) can be expressed as

F - 27rV2£G/(Bty (4)

G = F x cos (/> sin </> (5)

and x is the speed ratio given by

DISC PLANE
Fig. 1 Flow geometry for blade element at radial station r.
 678 ADKINS AND LIEBECK: DESIGN OF OPTIMUM PROPELLERS

Constraint Equations
For design, it is necessary to specify either 7, delivered by
the propeller or the power P, delivered to the propeller. The
nondimensional thrust and power coefficients used for design
are

Tc = 2T/(pV2irR2) (9a)
3 2 3 2
Pc = 2P/(pV 7rR ) = 2Q(l/(pV >7TR ) (9b)

and using these definitions, Eq. (6) can be written as

T'c = I(£ - Itf (lOa)

P'C=J[{ + J& (lOb)
where the primes denote derivatives with respect to f, and
DISC PLANE-
Fig. 3 Force diagram for a blade element. /; = 4£G(1 - e tan 0) (lla)
Downloaded by ORTA DOGU TEKNIK UNIVERSITESI on March 10, 2019 | http://arc.aiaa.org | DOI: 10.2514/3.23779

The circulation equations for thrust T , and torque Q', per 1'2 = A(/;/2£)(l + e/tan </>)sin <j> cos </> (lib)
unit radius can be written by inspection of Fig. 3 as /; - 4fG(l + e/tan 0) (lie)
2
T = L' cos 0 - D' sin 4> = L' cos 0(1 - e tan <£) (6a) J'2 = (/;/2)(l - s tan </>)cos 0 (lid)

Since £ is constant for an optimum design, a specified thrust

Q'lr = L' sin <£ + D' cos 0 = L' sin </>(! + e/tan gives the constraint equations
(6b)
£ = (A/2/2) - [(A/2/,)2- TC/I2]»2 (12)
where e is the drag-to-lift-ratio of the blade element. Next,
using Eq. (2), L' can be replaced by F(r) which, in turn, is PC = J^ + / 2 f 2 (13)
related to conditions in the wake by Eq. (3). Based on the
flow in the wake, F(r) is given by Eqs. (4) and (5), and T Similarly, if power is specified, the constraint relations are
and Q'lr are reduced to being functions of <f> and the dis-
placement velocity, £ = v'lV. The local flow angle <f> will £ = ~(JJ2J2) + [(V2/2)2 + PJJ^'2 (14)
clearly be a function of the radius; however, at this stage of
the analysis, the optimum distribution £(r) is not yet deter- Tc = U - I2£2 (15)
mined. Several diagrams and an excellent photograph of the
vortex sheet can be found in a 1980 work by Larrabee.9 where the integration has been carried out over the region
f = & t o f = 1.
Condition for Minimum Energy Loss Blade Geometry
At this point, a departure from Larrabee's4 design proce- For the element dr of a single blade at radial station r, let
dure is made, and the momentum equations, Eqs. (1), and c be the chord and Cl the local lift coefficient. Then, the lift
the circulating equations, Eqs. (6), are required to be equiv- per unit radius of one blade is
alent. This condition results in the interference factors being
related to £ by the equations = Pwr
2
a = (£/2)cos (£(l - £ tan </>) (7a) where F is given by Eq. (4). It follows directly that

a' '= (£/2*)cos (/> sin <£(! + e/tan </>) (7b) We = 47rXGVR£/(ClB) (16)

Assume for the moment that £ is known; then the local value
where Eqs. (4) and (5) have been used to express L' in terms of cj) is known from Eq. (8), and the above relation is a
of £, and the terms in epsilon correctly describe the viscous function only of the local lift coefficient. Since the local Rey-
contribution. Equations (7), together with the geometry of nolds number is We divided by the kinematic viscosity, Eq.
Fig. 1, lead to the important simple relation (16) plus a choice for C/ will determine the Reynolds number
and £, from the airfoil section data. The total velocity is then
tan <£ = £/2)/x - (1 + £/2)X/( (8) determined by Fig. 1 as

Here, A is a constant, and £ varies from £0 at the hub to unity W = V(l + fl)/sin (17)
at the edge of the disc. The relation between the two non-
dimensional distances and the constant speed ratio is where a is given by Eq. (7), and the chord is then known from
Eq. (16). If the choice for C/ causes £ to be a minimum, then
=
viscous as well as momentum losses will in most cases be
x = (r/R)/\ = f/A minimized, and overall propeller efficiency will be the highest
possible value. For preliminary considerations, it is usually
Recalling the Betz2 condition, r tan </> = const, Eq. (8) sufficient to choose one Ch the design Ch for determining
proves that for the vortex sheet to be a regular screw surface, blade geometry. (Any Cl specification is permissible as long
£ must be a constant independent of radius. This is the con- as the optimum blade loading distribution, cC/(r), is main-
dition for minimum energy loss. It should be noted that Eq. tained.) Since a is known from C, and Reynolds number, the
(8) results from Eq. (7) whether viscosity is included or not. blade twist with respect to the disc is j8 = a. + <f). G is zero
 ADKINS AND LIEBECK: DESIGN OF OPTIMUM PROPELLERS 679

at the edge of the disc, and the tip chord is therefore always
zero for a finite lift coefficient.

Design Procedure
Either F or K, relation for the momentum loss function can
be selected. For the sake of simplicity, only the Prandtl re-
lation is described as

F = (2/77)arc (18)

where

f= (19)

and </>, is the flow angle at the tip. From Eq. (8)
DISC PLANE-
tan <t>t = A(l + £/2) (20) Fig. 4 Force coefficients for propeller blade element analysis.
Downloaded by ORTA DOGU TEKNIK UNIVERSITESI on March 10, 2019 | http://arc.aiaa.org | DOI: 10.2514/3.23779

so that a choice for £ determines the function F as well as <f> and the relations for the thrust 7" and torque Q' per unit
by radius are then
tan $ = (tan <£, (21) T = ($)pW2BcCy (22a)
2
which is simply the condition that the vortex sheet in the wake Q'lr = ($)PW BcCx (22b)
is a rigid screw surface (r tan <£ = const). For an initial value,
£ = 0 will suffice. Again, it is required that the loading Eqs. (22) be exactly
The design is initiated with the specified conditions of power equal to the momentum result Eqs. (1). With the use of the
(or thrust), hub and tip radius, rotational rate, freestream flow geometry in Fig. 1, this requires the interference factors
velocity, number of blades, and a finite number of stations to be
at which blade geometry is to be determined. Also, the design
lift coefficient—one for each station if it is not constant— a = o-KI(F - <rK) (23a)
must be specified. The design then proceeds in the following a' = <rK'/(F + aK') (23b)
steps:
1) Select an initial estimate for £ (£ = 0 will work). where
2) Determine the values for F and <j) at each blade station
by Eqs. (18-21).
3) Determine the product We, and Reynolds number from K = Cy/(4 sin2 (/>) (24a)
Eq. (16). K' = Cx/(4 cos 0 sin </>) (24b)
4) Determine e and a from airfoil section data.
5) If e is to be minimized, change C, and repeat Steps 3 and cr is given by
and 4 until this is accomplished at each station.
6) Determine a and a' from Eq. (7), and Wfrom Eq. (17). er = Bc/(2m)
7) Compute the chord from step 3, and the blade twist
)3 = a + (f>. Equations (23) correct the placement of the factor F used by
8) Determine the four derivatives in / and / from Eq. (11) Glauert in his equations (5.5) of Chapter VII as identified by
and numerically integrate these from £ = £0 to £ = 1. Larrabee.4
9) Determine £ and Pc from Eqs. (12) and (13), or £ and The relation for the flow angle is obtained from Fig. 1 and
Tc from Eqs. (14) and (15). Eqs. (23) as
10) If this new value for £ is not sufficiently close to the
old one (e.g., within 0.1%) start over at step 2 using the
new £. tan <t> = [V(l + a)]/[nr(l - a')] (25)
11) Determine propeller efficiency as TJPC, and other fea-
tures such as solidity. For determining the function, F, in Eq. (18), Glauert suggests
The above steps converge rapidly, seldom taking more than the relation sin $, = f sin 4> be used in Eq. (19). It is rec-
three or four cycles. An accurate description of viscous losses ommended that Eq. (21) be used instead, i.e.,
can be obtained by creating another design with e equal to
zero and noting the difference in propeller efficiency. tan </>, = £ tan </>

which is exact for the analysis of an optimally designed pro-

Analysis of Arbitrary Designs peller at the design point.
The analysis method is outlined here in order to discuss The analysis procedure requires an iterative solution for
problems of convergence for off design and for square-tipped the flow angle </> at each radial position, £. An initial estimate
propellers in general, and to point out two minor errors in for </> can be obtained from Eq. (8) by setting £ equal to zero.
Glauert's work. Figure 4, which is simply an alternate version Since j3 is known, the value for a in Fig. 3 is /3 — <£, and the
of Fig. 3, shows the relation between the propeller force coef- airfoil coefficients are known from the section data. The
ficients, Cy and Cx, and the airfoil coefficients, C, and Cd. The Reynolds number is determined from the known chord and
equations are W, which is obtained from Fig. 1 and Eq. (23a), and the new
estimate for <£ is then found from Eq. (25). A direct substi-
Cy = C, cos <f> - Cd sin <£ = C,(cos </> - e sin </>) tution of the new <f> for the old value will cause adequate
convergence for an optimum design which is being analyzed
Cx = C, sin </> + Cd cos </> = C7(sin </> + e cos </>) at the design point. However, for analysis off-design and for
 680 ADKINS AND LIEBECK: DESIGN OF OPTIMUM PROPELLERS

nonoptimum designs, some recursive combination of the old mation of propeller performance, and the extent of airfoil
and new values for <f> is required to cause adequate conver- data required for propeller analysis exceeds that which is typ-
gence. Under some conditions (usually near the tip), con- ically used in analysis of wings. Between takeoff, climb, and
vergence may not be possible at all due to large values for cruise, propellers typically operate over a relatively wide range
the interference factors, a and a', in Eq. (23). Since Fis zero of blade section Reynolds numbers—a variation by factor of
at the tip and a is not for a square tip propeller, the value 5 is not uncommon. Also, during takeoff, climb, and wind-
for a is - 1 and a' is +1. Such values are physically impossible milling, some portion of a propeller blade is likely to be stalled,
since the slipstream factors are approximately twice the values either positively or negatively. An example set of character-
at the rotor plane. Wilson and Lissaman10 suggest empirical istics for the NACA 4415 airfoil is given in Fig. 5. These
relations for resolving this problem, whereas Viterna and would be supplemented by additional drag data in the un-
Janetzke11 give empirical arguments for clipping the magni- stalled region for a range of Reynolds numbers.
tude of a and a' at the value 0.7 (a/F at the tip is finite at the
design point for an optimum propeller). Empirical Optimality
For analysis, the conventional thrust and power coefficients In Chapter VII of Glauert's work, his equation (2.20) shows
are that when blade friction is neglected, the most favorable dis-
tribution of circulation is where the displacement velocity is
CT = TI(pn2D4) constant across the wake. Here, the term x2/(l + ;c2) is the
small-angle approximation of G, given by Eq. (5) in this ar-
Downloaded by ORTA DOGU TEKNIK UNIVERSITESI on March 10, 2019 | http://arc.aiaa.org | DOI: 10.2514/3.23779

Cp = P/(pn3D5) ticle. The effect of profile drag is shown by Glauert in his

equation (3.5), which states that the optimum distribution for
Using Eqs. (22) and (24), the differential forms with respect the displacement velocity ratio is
to f are given by
f = Io - (26)
C'T = <rK')cos <t>}2
where the effect of profile drag on thrust has been ignored.
c; = In order to study this problem empirically, consider a general
function H(x) and the two first-order terms of its Laurent
When these have been integrated from the hub to the tip, the series, IIx and x, and describe the displacement velocity dis-
propeller efficiency is tribution as
77 = CTJ/CP
£= £0 32/x (27)
where/ = V/(nD) is the advance ratio. Propeller performance
is typically described by plots of Cr, Cp, and 17 vs /. which includes the case of Eq. (26). It is desired to find values
for 8l and 82 which maximize propeller efficiency subject to
Airfoil Section Data the constraints of Eqs. (10). To solve this problem, £ in Eqs.
(10) is replaced by Eq. (27). Then, a choice for 5, and 52 will
The proper accounting for blade section (airfoil) charac- enable a determination of £0 and a calculation of overall pro-
teristics has proven essential for accurate and reliable esti- peller efficiency. A systematic study of various propeller con-
2.0-.
ditions was undertaken using the design and analysis proce-
dures of this article. Nonzero values for 5t and 82 that caused
an increase in propeller efficiency could not be found under
any conditions. Therefore, it was concluded that a constant
displacement velocity is at least locally optimum whether pro-
file drag is considered or not. Momentum and viscous losses
are then uncoupled; former is minimized by constant dis-
placement velocity, the latter by choosing a C, distribution so
that the drag-to-lift ratio is a minimum everywhere.
Some may criticize the authors for including viscous terms
in the development of the optimized circulation distribution
-30 -20 70 80 for the propeller design problem. Certainly, this is contrary
to classical lifting-line wing theory where the analogous elliptic
lift distribution is obtained inviscidly. However, it has been
NACA4415
shown that the combined momentum-circulation equations,
RN = 0.5 x10 6
-1.0 J
Eqs. (7), produce the screw surface equations, Eqs. (8), whether
viscosity is included or not. It is evident in Eq. (11) that the
momentum and viscous terms are directly separable for use
0.2-, in Eq. (10). The viscous terms, when integrated in Eq. (10),
"C/4
account for the difference in power and thrust between an
0.1- inviscid and viscous propeller design. If one accepts the clas-
sical blade element analysis equations as a measure of per-
-30 -20 10 20 30 40 50 70 80 formance, the momentum and viscous losses are indeed un-
a (DEG) coupled, and the viscous design equations produce a propeller
with minimal losses.
Alternatively, if one is still concerned about including vis-
cosity in the design procedure, a propeller could be designed
inviscidly and the blade section drag simply added to the
-.03-
inviscid design. This would require that the analysis procedure
be applied with drag set equal to zero when solving for the
-.04 J
induction velocities. The performance of the propeller would
Fig. 5 Example of airfoil section characteristics including normal and then be obtained by adding drag to this inviscid solution with-
chord force coefficients. out altering the induction velocities. However, this would
 ADKINS AND LIEBECK: DESIGN OF OPTIMUM PROPELLERS 681

Table 1 Propeller design solution

r c ft * RN a a'
0.5000 0.3424 58.3125 54.8118 0.4449 0.0348 0.0633
0.8958 0.4605 41.8645 38.3637 0.8104 0.0644 0.0365
1.2917 0.4269 32.2669 28.7661 0.9834 0.0804 0.0219
1.6875 0.3569 22.2978 22.7927 1.0295 0.0890 0.0142
2.0833 0.2796 18.7971 18.7971 0.9740 0.0938 0.0098
2.4792 0.1913 15.9619 15.9619 0.7830 0.0968 0.0072
2.8750 0.0000 13.8552 13.3552 0.0000 0.0000 0.0000
Input: brake horsepower = 70, 2 blades; hub diam = 1 ft, tip diam = 5.75 ft; blade section: NACA
4415, C, = 0.7, velocity = 110 mph, rpm = 24001.
Output: thrust = 207.61 Ib, 77 = 0.86996.
Note: a and a' have been set equal to zero at the tip.

Table 2 Propeller analysis solution

r * c, RN a a'
Downloaded by ORTA DOGU TEKNIK UNIVERSITESI on March 10, 2019 | http://arc.aiaa.org | DOI: 10.2514/3.23779

0.5000 54.8116 0.7000 0.4449 0.0348 0.0633

0.8958 38.3638 0.7000 0.8104 0.0644 0.0365
1.2917 28.7661 0.7000 0.9834 0.0804 0.0219
1.6875 22.7927 0.7000 1.0295 0.0890 0.0142
2.0833 18.7971 0.7000 0.9740 0.0938 0.0098
2.4792 15.9619 0.7000 0.7830 0.0968 0.0072
2.8750 12.5862 0.7000 0.0000 0.0000 0.0000
Input: propeller geometry from Table 1; r, C, and /3 at the same radial locations; velocity =110 mph,
rpm = 2400.
Output: brake horsepower = 70, thrust = 207.61 Ib, 17 = 0.86996.

1.6

1.4

1.2
1.0
I
0.8
"7 DESIGN

0.6

0.4

DISC PLANE 0.2

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
rR
/ T.P
Fig. 6 Force coefficients for windmill blade element.
Fig. 8 C, distributions for example propeller.

0.08 0.8

0.06 0.6 0.20 1.0

CP.CT
0.8
0.04 0.15

0.02

0.05
0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
REF: NACA TN 1834, PROP MODEL 5

Fig. 7 Example of propeller performance.

require an additional layer of iteration to achieve a specified Fig. 9 Comparison of theory and experiment.
design thrust or power. In light of the favorable agreement
between the present theory and the experimental results given geometry for a windmill is shown in Fig. 6, where the primary
later in this article, it is argued that such an increase in com- distinction is that the blade section is inverted (as compared
plexity is not justified. with a propeller), and the local angle of attack is measured
from below the local velocity vector. Corresponding relations
for the angles are
Windmills
All of the analyses described in this article are directly windmill a = </> - /3
applicable to the windmill problem after a minor adjustment
in the angle definitions of Fig. 1. The corresponding flow propeller a = ft — $
 682 ADKINS AND LIEBECK: DESIGN OF OPTIMUM PROPELLERS

0.3 Conclusions and Recommendations

J = 0.914
THEORY WITH PRANDTL F
- — — THEORY WITH GOLDSTEIN K
0 NACATEST
The propeller theory of Glauert has been extended to im-
prove the design of optimal propellers and refine the calcu-
lation of the performance of arbitrary propellers. Extensions
of the theory include 1) elimination of the small angle as-
sumptions in the optimal design theory; 2) accurate calcula-
tion of the vortex displacement velocity which properly ac-
counts for the blade section drag; and 3) elimination of the
small angle assumptions in the Prandtl momentum loss func-
tion for both design and analysis. These extensions bring the
0.1 design and analysis procedures to exact numerical agreement
within the precision of computer analysis.
The primary approximation remaining in both procedures
is the use of the axial momentum equations which require the
Downloaded by ORTA DOGU TEKNIK UNIVERSITESI on March 10, 2019 | http://arc.aiaa.org | DOI: 10.2514/3.23779

REF: NACA TN 1834, PROP MODEL 5

increase in wake velocities to be twice those at the disc. Under
certain conditions this approximation is not good and gives
0 0.2 0.4 0.6 0.8 1.0 rise to the unnatural conditions and convergence problems
r/R-np
described in the analysis section. Improvements might be made
Fig. 10 Comparison of propeller analyses thrust coefficient. by replacing the axial momentum equations with relations
more closely aligned with the general theory, particularly in
as shown in Figs. 6 and 1, respectively. In these figures, C, those differential stream tubes in which "heavy loading" ex-
for the windmill is negative with respect to that for the pro- ists. Such conditions appear to be more prevalent in the anal-
peller, and this sign change together with the angle definition yses at off-design conditions than in the design itself, and,
will convert the propeller methods to the windmill application. when combined with poststall misknowledge, can lead to large
For the design case, the input Pc value should be negative, errors in analysis. However, for design and analysis within
and the resulting values of v' (and the interference factors a the conventional operating regime, both procedures are sim-
and a') and Tc will also be negative. (Thrust is of less interest ple, accurate, and reliable. This method has been extended
for a windmill since it typically represents the tower load and by Page and Liebeck13 to the design and analysis of dual-
is not a main performance parameter.) Similarly, the analysis rotation propellers. A favorable comparison between theory
results for a windmill rotor will yield negative values for both and experiment was also observed.
Pc and Tc.
Examples
References
As a sample calculation, the design of a propeller for a light
airplane is considered. The design conditions and the resulting l
Glauert, H., "Airplane Propellers, "Aerodynamic Theory, edited
design are described in Table 1 which gives for each radial by W. Durand, Div. L, Vol. 5, Peter Smith, Gloucester, MA, 1976,
station: blade chord, blade pitch angle, local flow angle, local pp. 169-269.
2
Reynolds number, and the interference coefficients a and a'. Betz, A., with appendix by Prandtl, L., "Screw Propellers with
This propeller geometry has, in turn, been analyzed at the Minimum Energy Loss," Gottingen Reports, 1919, pp. 193-213.
3
design condition and the result is given in Table 2. Agreement Theodorsen, T., Theory of Propellers, McGraw-Hill, New York,
is virtually exact. Analysis over a range of values of the ad- 1948.
4
vance ratio, J = V/(nD), provides the typical propeller per- Larrabee, E., "Practical Design of Minimum Induced Loss Pro-
pellers," Society of Automotive Engineers, Business Aircraft Meet-
formance plots which are shown in Fig. 7, and Fig. 8 gives ing and Exposition, Wichita, KS, April 1979.
the blade lift coefficient distribution over a range of /s where 5
Goldstein, S., "On the Vortex Theory of Screw Propellers," Pro-
the design condition is the C, — 0.7 and is a constant line at ceedings of the Royal Society of London, Series A, Vol. 123, 1929,
J = 0.7. pp. 440-465.
6
A calibration of the method is given by comparing its the- Lock, C., "The Application of Goldstein's Airscrew Theory to
oretical prediction with experimental results. Reid12 has eval- Design," British Aeronautical Research Committee, RM 1377, Nov.
uated several conventional propellers extensively by experi- 1930.
7
ment, and one of these has been chosen for comparison. Lock, C., "An Application of Prandtl's Theory to an Airscrew,"
Figure 9 gives Cp, CT, and 77 vs / for both Reid's experiments British Aeronautical Research Committee, RM 1521, Aug. 1932.
8
Lock, C., "Tables for Use in an Empirical Method of Airscrew
and the corresponding theoretical prediction. The agreement Strip Theory Calculations," British Aeronautical Research Commit-
here is quite good, with most of the disparity occurring after tee, RM 1674, Oct. 1934.
the blade is stalled. This propeller uses NACA 16-series air- 9
Larrabee, E., "The Screw Propeller," Scientific American, Vol.
foils, and no poststall data were available. 243, No. 1, 1980, pp. 134-148.
10
Figure 10 gives the comparison of the blade thrust coeffi- Wilson, R., and Lissaman, P., "Applied Aerodynamics of Wind
cient distribution as measured by Reid and calculated by the Power Machines," Oregon State Univ., NSF/RA/N-74-113, PB-
method. Two theoretical results are shown: one using F, and 2318595/3, Corvallis, OR, July 1974.
H
the other using the more complex (from a calculation point Viterna, A., and Janetzke, D., "Theoretical and Experimental
of view) K. In principle, the accuracy of the method should Power from Large Horizontal-Axis Wind Turbines," Proceedings from
the Large Horizontal-Axis Wind Turbine Conference, DOE/NASA-
be better with the Goldstein factor for a propeller with few LeRC, July 1981.
blades—this example had three blades—and the two factors 12
Reid, E. G., "The Influence of Blade-Width Distribution on
should give similar results as the number of blades is in- Propeller Characteristics," NACA TN 1834, March 1949.
creased. The results of Fig. 10 confirm this trend, and the 13
Page, G. S., and Liebeck, R. H., "Analysis of Dual-Rotation
overall comparison for both factors is regarded as quite good. Propellers," AIAA Paper 89-2216, Aug. 1989.