Fundamental of Aerodynamic 5th Ed
Fundamental of Aerodynamic 5th Ed
Fundamental of Aerodynamic 5th Ed
PART 1
CHAPTER 1
Fundamental Principles
dp
= -gpdy
(1.52)
p
or
+ pgh =
Pt + pgh1
constant
(1.54)
= P2 + pgh2
Such equations govern, among other things, the operation of a manometer, and
also lead to Archimedes' principle that the buoyancy force on a body immersed
in a fluid is equal to the weight of the fluid displaced by the body.
edge, and moments about the quarter chord, all per unit span. Also,
calculate the location of the center of pressure.
1.s Consider an airfoil at 12 angle of attack. The normal and axial force
coefficients are 1.2 and 0.03, respectively. Calculate the lift and drag
coefficients.
1.6 Consider an NACA 2412 airfoil (the meaning of the number designations
for standard NACA airfoil shapes is discussed in Chapter 4). The
following is a tabulation of the lift, drag, and moment coefficients about
the quarter chord for this airfoil, as a function of angle of attack.
a (degrees)
- 2.0
0
2.0
4.0
6.0
8.0
10.0
12.0
14.0
1.16 PROBLEMS
1.1 For most gases at standard or near standard conditions, the relationship
among pressure, density, and temperature is given by the perfect gas
equation of state: p = pRT, where R is the specific gas constant. For air
at near standard conditions, R = 287 J/ (kg K). (More details on the
perfect gas equation of state are given in Chapter 7.) Using the above
information, consider the following two cases:
a. At a given point on the wing of a Boeing 727, the pressure and
temperature of the air are 1.9 x 104 N/m 2 and 203 K, respectively.
Calculate the density at this point.
b.
c,
Cd
Cm,c / 4
0.05
0.25
0.44
0.64
0.85
1.08
1.26
1.43
1.56
0.006
0.006
0.006
0.007
0.0075
0.0092
0.0115
0.0150
0.0186
- 0.042
- 0.040
- 0.038
- 0.036
- 0.036
-0.036
- 0.034
- 0.030
- 0.025
From this table, plot on graph paper the variation of Xcp/c as a function
ofa.
1.7 The drag on the hull of a ship depends in part on the height of the water
waves produced by the hull. The potential energy associated with these
waves therefore depends on the acceleration of gravity g. Hence, we can
state that the wave drag on the hull is D = f(p 00 , V00 , c , g) where c is a
length scale associated with the hull, say, the maximum width of the hull.
Define the drag coefficient as C O = D / q00 c2 . Also, define a similarity
parameter called the Froude number, Fr= V / .Jgc. Using Buckingham's
pi theorem, prove that Co = f (Fr).
1.8 The shock waves on a vehicle in supersonic flight cause a component of
drag called supersonic wave drag Dw. Define the wave-drag coefficient as
Co,w = Dw/q00 S, where Sis a suitable reference area for the body. In
supersonic flight, the flow is governed in part by its thermodynamic
properties, given by the specific heats at constant pressure c P and at
constant volume Cv. Define the ratio cp/cv = y. Using Buckingham's
pi theorem, show that Co,w = f(M00 , y). Neglect the influence of friction.
1.9 Consider two different flows over geometrically similar airfoil shapes, one
airfoil being twice the size of the other. The flow over the smaller airfoil
has freestream properties given by T00 = 200 K, p00 = 1.23 kg/m3 , and
V00 = 100 mis. The flow over the larger airfoil is described by
T00 = 800 K, p 00 = 1.739 kg/m3 , and V00 = 200 mis. Assume that both,
and a are proportional to T 112 Are the two flows dynamically similar?
99
100
PART 1
Fundamental Principles
1.10 Consider a Lear jet flying at a velocity of 250 mis at an altitude of 10 km,
where the density and temperature are 0.414 kg/m3 and 223 K,
respectively. Consider also a one-fifth scale model of the Lear jet being
tested in a wind tunnel in the laboratory. The pressure in the test section of
the wind tunnel is 1 atm = 1.01 x 105 N/m2 . Calculate the necessary
velocity, temperature, and density of the airflow in the wind-tunnel test
section such that the lift and drag coefficients are the same for the
wind-tunnel model and the actual airplane in flight. Note: The relation
among pressure, density, and temperature is given by the equation of state
described in Problem 1.1.
1.11 AU-tube mercury manometer is used to measure the pressure at a point
on the wing of a wind-tunnel model. One side of the manometer is
connected to the model, and the other side is open to the atmosphere.
Atmospheric pressure and the density of liquid mercury are
1.01 x 105 N/m2 and 1.36 x 104 kg/m3 , respectively. When the
displacement of the two columns of mercury is 20 cm, with the high
column on the model side, what is the pressure on the wing?
1.12 The German Zeppelins of World War I were dirigibles with the following
typical characteristics: volume= 15,000 m3 and maximum diameter=
14.0 m. Consider a Zeppelin flying at a velocity of 30 mis at a standard
altitude of 1000 m (look up the corresponding density in Appendix D).
The Zeppelin is at a small angle of attack such that its lift coefficient is
0.05 (based on the maximum cross-sectional area). The Zeppelin is flying
in straight-and-level flight with no acceleration. Calculate the total weight
of the Zeppelin.
1.13 Consider a circular cylinder in a hypersonic flow, with its axis
perpendicular to the flow. Let <P be the angle measured between radii
drawn to the leading edge (the stagnation point) and to any arbitrary point
on the cylinder. The pressure coefficient distribution along the cylindrical
surface is given by C P = 2 cos 2 <P for O ::: <P ::: 1r / 2 and 3n / 2 ::: ::: 2n
and CP = 0 for 1r / 2 ::: ::: 3n / 2. Calculate the drag coefficient for the
cylinder, based on projected frontal area of the cylinder.
1.14 Derive Archimedes' principle using a body of general shape.
1.15 Consider a light, single-engine, propeller-driven airplane similar to a
Cessna Skylane. The airplane weight is 13000 N and the wing reference
area is 16.533 m2 . The drag coefficient of the airplane CD is a function of
the lift coefficient CL for reasons that are given in Chapter 5; this function
for the given airplane is Cv = 0.025 + 0.054Cz.
a. For steady, level flight at sea level, where the ambient atmospheric
density is 1.225 kg/m3, plot on a graph the variation of CL, CD, and
the lift-to-drag ratio L / D with flight velocity ranging between 21 mis
and 75 mis.
b. Make some observations about the variation of these quantities with
velocity.
CHAPTER 1
1.16 Consider a flat plate at zero angle of attack in a hypersonic flow at Mach
10 at standard sea level conditions. At a point 0.5 m downstream from the
leading edge, the local shear stress at the wall is 282 N/m2 . The gas
temperature at the wall is equal to standard sea level temperature. At this
point, calculate the velocity gradient at the wall normal to the wall.
1.17 Consider the Space Shuttle during its atmospheric entry at the end of a
mission in space. At the altitude where the Shuttle has slowed to Mach 9,
the local heat transfer at a given point on the lower surface of the wing is
0.03 MW/m2 . Calculate the normal temperature gradient in the air at this
point on the wall, assuming the gas temperature at the wall is equal to the
standard sea-level temperature.
t.18 The purpose of this problem is to give you a feel for the magnitude of
Reynolds number appropriate to real airplanes in actual flight.
a. Consider the DC-3 shown in Figure 1.1. The wing root chord length
(distance from the front to the back of the wing where the wing joins
the fuselage) is 4.3 m. Consider the DC-3 flying at 320 kilometers per
hour at sea level. Calculate the Reynolds number for the flow over the
wing root chord. (This is an important number, because as we will see
later, it governs the skin-friction drag over that portion of the wing.)
b. Consider the F-22 shown in Figure 1.5, and also gracing the cover of
this book. The chord length where the wing joins the center body is
6.5 m. Consider the airplane making a high-speed pass at a velocity of
396 mis at sea level (Mach 1.2). Calculate the Reynolds number at the
wing root.
1.19 For the design of their gliders in 1900 and 1901, the Wright brothers used
the Lilienthal Table given in Figure 1.65 for their aerodynamic data. Based
on these data, they chose a design angle of attack of 3 degrees, and made
all their calculations of size, weight, etc., based on this design angle of
attack. Why do you think they chose three degrees?
Hint: From the table, calculate the ratio of lift to drag, UD, at 3 degrees angle of
attack, and compare this with the lift-to-drag ratio at other angles of attack. You
might want to review the design box at the end of Section 1.8, especially Figure
1.36, for the importance of UD.
101
198
PART 1
CHAPTER 2
Fundamental Principles
The velocity potential <I> is defined for irrotationa l flows only, such that
(2. 154)
2.7 The velocity field given in Problem 2.3 is called source flow, which will be
In cartesian coordinate s
a<1>
u= -
ax
a<1>
w= -
a<1>
v= -
az
ay
(2.156)
2.8
In cylindrical coordinate s,
a<1>
V, = -
ar
1 a<1>
Ve = - -
ae
(2. 157)
In spherical coordinate s,
a<1>
V, = -
ar
1 a<1>
Ve= - -
ae
1 a<1>
V<1> = - - rsin () act>
(2. 158)
A line of constant <I> is an equipotential line. Equipoten tial lines are perpendicular to streamlines (for two-dimen sio nal irrotationa l flows).
2.19 PROBLEMS
2.9
2.10
2.11
2.12
2.1 Consider a body of arbitrary shape. If the pressure distributio n over the
2.2
2.3
2.4
2.5
2.6
surface of the body is constant, prove that the resultant pressure force on
the body is zero. [Recall that this fact was used in Equation (2.77).]
Consider an airfoil in a wind tunnel (i.e., a wing that spans the entire test
section). Prove that the lift per unit span can be o btained from the pressure
distributio ns on the top and bottom walls of the wind tunnel (i.e., from the
pressure distributio ns o n the walls above and below the airfoil).
Consider a velocity field w here the x and y componen ts of velocity are
2
2
2
given by u = ex/ (x + y ) and u = ey / (x + y2) where e is a constant.
Obtain the equations of the streamline s.
Consider a velocity field where the x and y components of velocity are
2
2
2
2
given by u = cy / (x + y ) and v = -ex / (x + y ), where c is a constant.
Obtain the equations of the stream lines.
Consider a velocity field where the radial and tangential componen ts of
velocity are V, = 0 and Ve = er , respectively, where c is a constant.
Obtain the equations of the streamline s.
Consider a velocity field where the x and y componen ts of velocity are
given by u = ex and v = -ey, where c is a constant. Obtain the equations
of the streamlines.
2.13
2.14
199
308
PART 2
CHAPTER 3
Elementary flows
Surface velocity:
(3.53)
1/f
(3.55)
= -lnr
1/f =
V, = -
Ve=O
21tr
<P
K sine
21t - r -
(3 .87)
(3.112)
1/f = - lnr
(3.114)
</J
=-
21t
(3.140)
D=O
Kutta-Joukowski theorem
for a closed two-dimensional body of arbitrary shape, the lift per unit span is
L' = Poo Voof.
(3 .88)
21t _r_
(3.62)
K Cose
(3.125)
21t R
(3.72)
(3.67)
2Jt
A
21t e
2Jt
r
Ve = - -
_.:. + I: - 1
2
1-1
2Jt
an;
00
cos,B; = 0
(i
= 1, 2, ... , n)
(j ,'1 )
V, =0
21tr
(3.152)
(3.105)
3.22 PROBLEMS
Note: All the following problems assume an inviscid, incompressible flow. Also,
standard sea level density and pressure are 1.23 kg/m3 and 1.01 x 105 N/m2 ,
respectively.
3.1 For an irrotational flow, show that Bernoulli's equation holds between any
(3.92)
(3.100)
C P = 1 - 4 sin2 e
(3.101)
(1 - R:)
r
_!::_ In ..c_
21t
(3.118)
(continued)
3.2 Consider a venturi with a throat-to-inlet area ratio of 0.8, mounted on the
side of an airplane fuselage. The airplane is in flight at standard sea level.
If the static pressure at the throat is 1.006 x 105 N/m2 , calculate the
velocity of the airplane.
3.3 Consider a venturi with a small hole drilled in the side of the throat. This
hole is connected via a tube to a closed reservoir. The purpose of the
venturi is to create a vacuum in the reservoir when the venturi is placed in
an airstream. (The vacuum is defined as the pressure difference below the
outside ambient pressure.) The venturi has a throat-to-inlet area ratio of
0.85. Calculate the maximum vacuum obtainable in the reservoir when the
venturi is placed in an airstream of 90 mis at standard sea level conditions.
3.4 Consider a low-speed open-circuit subsonic wind tunnel with an
inlet-to-throat area ratio of 12. The tunnel is turned on, and the pressure
309
CHAPTER 3
310
PART 2
difference between the inlet (the settling chamber) and the test section is
read as a height difference of 10 cm on a U-tube mercury manometer.
(The density of liquid mercury is 1.36 x 104 kg/m 3 .) Calculate the
velocity of the air in the test section.
3.5 Assume that a Pi tot tube is inserted into the test-section flow of the wind
tunnel in Problem 3.4. The tunnel test section is completely sealed from
the outside ambient pressure. Calculate the pressure measured by the Pitot
tube, assuming the static pressure at the tunnel inlet is atmospheric.
3.6 A Pitot tube on an airplane flying at standard sea level reads 1.07 x
105 N/m2 . What is the velocity of the airplane?
3. 7 At a given point on the surface of the wing of the airplane in Problem 3.6,
the flow velocity is 130 mis. Calculate the pressure coefficient at this point.
3.8 Consider a uniform flow with velocity V00 Show that this flow is a
physically possible incompressible flow and that it is irrotational.
3.9 Show that a source flow is a physically possible incompressible flow
everywhere except at the origin. Also show that it is irrotational
everywhere.
3.10 Prove that the velocity potential and the stream function for a uniform
flow, Equations (3.53) and (3.55), respectively, satisfy Laplace's equation.
3.11 Prove that the velocity potential and the stream function for a source flow
Equations (3.67) and (3.72), respectively, satisfy Laplace's equation.
'
3.12 Consider the flow over a semi-infinite body as discussed in Section 3.11. If
V00 is the velocity of the uniform stream, and the stagnation point is 1 m
upstream of the source:
a. Draw the resulting semi-infinite body to scale on graph paper.
b. Plot the pressure coefficient distribution over the body; that is, plot Cp
versus distance along the centerline of the body.
3.13 Derive Equation (3.81). Hint: Make use of the symmetry of the flow field
shown in Figure 3.18; that is, start with the knowledge that the stagnation
points must lie on the axis aligned with the direction of V00
3.14 Derive the velocity potential for a doublet; that is, derive Equation (3.88).
Hint: The easiest method is to start with Equation (3.87) for the stream
function and extract the velocity potential.
3.15 Consider the nonlifting flow over a circular cylinder. Derive an expression
for the pressure coefficient at an arbitrary point (r , 8) in this flow, and
show that it reduces to Equation (3.101) on the surface of the cylinder.
3.16 Consider the nonlifting flow over a circular cylinder of a given radius,
where V00 = 6 mis. If V' is doubled, that is, V00 = 12 mis, does the shape
of the streamlines change? Explain.
3.17 Consider the lifting flow over a circular cylinder of a given radius and with
a given circulation. If V00 is doubled, keeping the circulation the same,
does the shape of the streamjjnes change? Explain.
3.19 A typical World War I biplane fighter (such as the French SPAD shown in
Figure 3.50) has a number of vertical interwing struts and diagonal
bracing wires. Assume for a given airplane that the total length for the
vertical struts (summed together) is 7.5 m, and that the struts are
cylindrical with a diameter of 0.05 m. Assume also that the total length
of the bracing wires is 24 m, with a cylindrical diameter of 0.0023 m.
Calculate the drag contributed by these struts and bracing wires when
the airplane is flying at 192 km/h at standard sea level. Compare this
component of drag with the total zero-lift drag for the airplane, for which
the total wing area is 21.4 m2 and the zero-lift drag coefficient is 0.036.
3.20 The Kutta-Joukowski theorem, Equation (3.140), was derived exactly
for the case of the lifting cylinder. In Section 3.16 it is stated without
proof that Equation (3.140) also applies in general to a two-dimensional
body of arbitrary shape. Although this general result can be proven
mathematically, it also can be accepted by making a physical argument as
well. Make this physical argument by drawing a closed curve around the
body where the closed curve is very far away from the body, so far away
that in perspective the body becomes a very small speck in the middle of
the domain enclosed by the closed curve.
3.21 Consider the streamlines over a circular cylinder as sketched at the right of
Figure 3.26. Single out the first three streamlines flowing over the top of
the cylinder. Designate each streamline by its stream function, o/1, o/2, and
,jl3 . The first streamline wets the surface of the cylinder; designate ,jf1 = 0.
The streamline above that is \11 2, and the next one above that is o/3
Assume the streamlines start out in the free stream equally spaced. Hence,
the volume flow rates between the streamlines are the same. The
streamline ,jJ2 passes through the point (l.2R,7t/2) directly above the top
of the cylinder. Calculate the location of the point directly above the top of
the cylinder through which the streamjjne ,jJ3 flows. Comment on the
spacing between the streamlines directly above the top.
3.22 Consider the flow field over a circular cylinder mounted perpendicular to
the flow in the test section of a low-speed subsonic wind tunnel. At
standard sea-level conditions, if the flow velocity at some region of the
flow field exceeds about 400 km/h, compressibility begins to have an
effect in that region. Calculate the velocity of the flow in the test section of
the wind tunnel above which compressibility effects begin to become
important, i.e., above which we can not accurately assume totally
incompressible flow over the cylinder for the wind tunnel tests.
3.23 Prove that the flow field specified in Example 2.1 is not incompressible;
i.e., it is a compressible flow as stated without proof in Example 2.1.
311
406
PART 2
c HA PT EA
T_he Kutta condition is an observation that for a lifting airfoil of given shape at a
give~ an~Je o~ attack, nature adopts that particular value of circulation around
th~ ~1rfoil which results in the flow leaving smoothly at the trailing edge. If the
tra~l~ng-edge _angle is finite, then the trailing edge is a stagnation point. If the
trailing e??e 1s cusped, the_n the velocities leaving the top and bottom surfaces
at the trrulmg edge are firute and equal in magnitude and direction. In either
case,
y(TE)
=0
(4. 10)
The vortex panel method is an important numerical technique for the solution
of the inviscid, incompressible flow over bodies of arbitrary shape, thickness,
and angle of attack. For panels of constant strength, the governing equations
are
y ;; oO
~ ds
I: -'
on;
2Jr
n
j= I
and
=o
(i
= 1,2, . .. , n)
Y;
= -Yi - I
which is one way of expressing the Kutta condition for the panels immediately
above and below the trailing edge.
Thin airfoil theory is predicated on the replacement of the airfoil by the mean
ca~ber line. A vortex sheet is placed along the chord line, and its strength
~dJusted. such that, in conjunction with the uniform freestream, the camber
hne becom~~ a streamline of the flow while at the same time satisfying the
Kutta cond1t:Ion. The strength of such a vortex sheet is obtained from the
fundamental equation of thin airfoil theory:
d~ = Voo (a - dxdz)
_I ( y(~)
2rr Jo x - ~
(4. 18)
Symmetric airfoil
1.
2.
3.
4.
c1 = 2rra.
4.1 Consider the data for the NACA 2412 airfoil given in Figure 4.10.
Calculate the lift and moment about the quarter chord (per unit span)
for this airfoil when the angle of attack is 4 and the freestream is at
standard sea level conditions with a velocity of 15 mis. The chord of the
airfoil is 0.6 m.
4.2 Consider an NACA 2412 airfoil with a 2-m chord in an airstream with a
velocity of 50 mis at standard sea level conditions. If the lift per unit span
is 1353 N, what is the angle of attack?
4.3 Starting with the definition of circulation, derive Kelvin's circulation
theorem, Equation (4.1 1).
4.4 Starting with Equation (4.35), derive Equation (4.36).
4.5 Consider a thin, symmetric airfoil at 1.5 angle of attack. From the results
of thin airfoil theory, calculate the lift coefficient and the moment
coefficient about the leading edge.
4.6 The NACA 4412 airfoil has a mean camber line given by
Cm ,ac
= 0.
Cambered airfoil
1.
c,
= 2rr [a +
dz
Lift slope
dc1/da = 2rr.
3. The aerodynamic center is at the quarter-chord point.
4. The center of pressure varies with the lift coefficient.
2.
4.17 PROBLEMS
(4.57)
407
408
PART 2
4.9 Starting with Equations (4.35) and (4.43), derive Equation (4.62).
4.10 For the NACA 2412 airfoil, the lift coefficient and moment coefficient
about the quarter-chord at - 6 angle of attack are -0.39 and - 0.045,
respectively. At 4 angle of attack, these coefficients are.0.65 and -0.037
respectively. Calculate the location of the aerodynamic center.
'
4.11 Consider again the NACA 2412 airfoil discussed in Problem 4.10. The
airfoil is flying at a velocity of 60 mis at a standard altitude of 3 km (see
Appendix D). The chord length of the airfoil is 2 m. Calculate the lift per
unit span when the angle of attack is 4.
4.12 For the -airfoil in Problem 4.11, calculate the value of the circulation
around the airfoil.
4.13 In Section 3.15 we studied the case of the lifting fl.ow over a circular
cylinder. In real life, a rotating cylinder in a flow will produce lift; such
real flow fields are shown in the photographs in Figures 3.34(b) and (c).
Here, the viscous shear stress acting between the flow and the surface of
the cylinder drags the flow around in the direction of rotation of the
cylinder. For a cylinder of radius R rotating with an angular velocity w in
an otherwise stationary fluid, the viscous flow solution for the velocity
field obtained from the Navier-Stokes equations (Chapter 15) is
R2w
Vo = -r
where V0 is the tangential velocity along the circular streamlines and r is
the radial distance from the center of the cylinder. (See Schlichting,
Boundary-Layer Theory, 6th ed., McGraw-Hill, 1968, page 81.) Note that
Vo varies inversely with r and is of the same form as the inviscid flow
velocity for a point vortex given by Equation (3. 105). If the rotating
cylinder has a radius of 1 m and is flying at the same velocity and altitude
as the airfoil in Problem 4.11 , what must its angular velocity be to produce
the same lift as the airfoil in Problem 4.11? (Note: You can check your
results with the experimental data for lift on rotating cylinders in Hoerner,
Fluid-Dynamic Lift, published by the author, 1975, p. 21-4, Fig. 5.)
4.14 The question is often asked: Can an airfoil fly upside-down? To answer
this, make the following calculation. Consider a positively cambered
airfoil with a zero-lift angle of -3. The lift slope is 0.1 per degree.
(a) Calculate the lift coefficient at an angle of attack of 5. (b) Now
imagine the same airfoil turned upside-down, but at the same 5 angle of
attack as part (a). Calculate its lift coefficient. (c) At what angle of attack
must the upside-down afrfoil be set to generate the same lift as that when it
is right-side-up at a 5 angle of attack?
4.15 The airfoil section of the wing of the British Spitfire of World War II fame
(see Figure 5.19) is an NACA 2213 at the wing root, tapering to an NACA
2205 at the wing tip. The root chord is 2.5 m. The measured profile drag
coefficient of the NACA 2213 airfoil is 0.006 at a Reynolds number of
c HA PT ER
409
484
PART 2
CHAPTER 5
ro
w= - -
(5.35)
2b
CL
a - -' - nAR
(5.42)
Cz
(5.43)
C0 . ' -
nAR
ao
Q= - --
--
1 + a 0 /nAR
(5.69)
General wing:
c2
c2
n AR
neAR
ao
a = -- -- - - --
1 + (ao/n AR)(I
+ r)
(5.70)
For low-aspect-ratio wings, swept wings, and delta wings, lifting-surface theory must be used. In modern aerodynamics, such lifting-surface theory is
implemented by the vortex panel or the vortex lattice techniques.
5.10 PROBLEMS
5.1 Consider a vortex filament of strength r in the shape of a closed circular
loop of radius R. Obtain an expression for the velocity induced at the
center of the loop in terms of r and R .
5.2 Consider the same vortex filament as in Problem 5.1. Consider also a
straight line through the center of the loop, perpendicular to the plane of
the loop. Let A be the distance along this line, measured from the plane of
the loop. Obtain an expression for the velocity at distance A on the line, as
induced by the vortex filament.
The measured lift slope for the NACA 23012 airfoil is 0.1080 degree- ' ,
and a L = O = -1 .3. Consider a finite wing using this airfoil, with AR = 8
and taper ratio = 0.8. Assume that o = r. Calculate the lift and induced
drag coefficients for this wing at a geometric angle of attack = 7.
CY
Downwash is constant:
485
486
PART 2
C?nsider t~e Spitfir~ in Problem 5.9 on its landing approach at sea level
with a landmg velocity of 112 km/h. Calculate the induced drag
coefficient for this low-speed case. Compare your result with the
high~speed case in Problem 5.9. From this, what can you conclude about
the relative importance of the induced drag coefficient at low speeds
compared to that at high speeds?
Three-Dimensional
Incompressible Flow
Treat nature in terms of the cylinder, the sphere, the cone, all in perspective.
Paul Cezanne, 1890
PREVIEW BOX
I:
6.1 INTRODUCTION
To this point in our aerodynamic discussions, we have been working mainl y in
a two-dimensional world; the flows over the bodies treated in Chapter 3 and the
airfoils in C hapter 4 involved only two dimensions in a single plane-so-called
planar flows. In Chapter 5, the analyses of a finite wing were carried out in the
plane of the wing, in spite of the fact that the detailed flow over a finite wing is truly
three-dimensional. The relative simplicity of dealing with two dimensions (i.e.,
having only two independent variables) is self-evident and is the reason why a
large bulk of aerodynamic theory deals with two-dimensional flows. Fortunately,
487
CHAPTER 7
546
PART. 3
#f
qpdV-
pV dS+
#f
p(f : V)dV
(7 .43)
D(e+V 2 /2)
.
Dt
=pq -V pV+p(f V)
(7.44)
~:the flow is steady and adiabatic, Equations (7.43) and (7.44) can be replaced
ho = h
v2
+ -2
= const
p =pRT
(7.1)
(7 .6a)
Total te~pe~a~re To and total enthalpy ho are defined as the properties that
would exist 1f (m our imagination) we slowed the fluid element at a poi t th
flow to zero velocity adiabatically. Similarly, total pressure Po and totatd~:sit e
~o are define~ as the properties that would exist if (in our imagination) w~
s owed the fluid ele_ment at a point in the flow to zero velocity isentropicall If
~!e~e~l flow_field 1s a?iab~tic, ho is ~onstant throughout the flow; in contra~-. if
o field 1s nonadiabattc, h 0 vanes from one point to another Similar! 1f
general fl~w field is isentropic, Po and Po are constant throughout th / m contrast, 1f the flow field is nonisentropic, Po and Po vary from one p:in~:~
another.
s:ock waves ~e very thin regions in a supersonic flow across which the
P ess~re, density, temperature, and entropy increase; the Mach number, flow
velocity, and total pressure decrease; and the total enthalpy stays the same.
7.8 PROBLEMS
Noted: All problems deal with calorically perfect air as the gas, unless otherwise
note . Also, recall that 1 atm = 1.01 x 10s N/m2.
7.3 Just upstream of a shock wave, the air temperature and pressure are 288 K
and 1 atm, respectively; just downstream of the wave, the air temperature
and pressure are 690 K and 8.656 atm, respectively. Calculate the changes
in enthalpy, internal energy, and entropy across the wave.
7.4 Consider the isentropic flow over an airfoil. The freestream conditions are
T = 245 Kand p 00 = 4.35 x 104 N/m2. At a point on the airfoil, the
00
pressure is 3.6 x 104 N/m2. Calculate the density at this point.
7.5 Consider the isentropic flow through a supersonic wind-tunnel nozzle. The
reservoir properties are To = 500 K and p 0 = 10 atrn. If p = l atm at the
nozzle exit, calculate the exit temperature and density.
7.6 Consider air at a pressure of 0.2 atm. Calculate the values of r r and r 5
Express your answer in SI units.
7.7 Consider a point in a flow where the velocity and temperature are 400 mis
and 267 K, respectively. Calculate the total enthalpy at this point.
7.8 In the reservoir of a supersonic wind tunnel, the velocity is negligible, and
the temperature is 1000 K. The temperature at the nozzle exit is 600 K.
Assuming adiabatic flow through the nozzle, calculate the velocity at the
exit.
= 0.819 kg/m3,
and V = 300 mis. At a point on the airfoil surface, the pressure is
00
0.5 atm. Assuming isentropic flow, calculate the velocity at that point.
7.10 Calculate the percentage error obtained if Problem 7 .9 is solved using
(incorrectly) the incompressible Bernoulli equation.
7.11 Repeat Problem 7.9, considering a point on the airfoil surface where the
7.9 An airfoil is in a freestream where p00 = 0.61 atm, p00
547
598
PART 3
P2
CHAPTER 8
U1
-= -=
P2
P1
[
S2 -
s,
(y
8.9 PROBLEMS
+ l)M f
2 + (y - l)Mf
= 1 + ~(M21 y +1
1)
+1
(y
l)M f
+ ])M?
(8.61)
(8.65)
the test section, the flow velocity is 420 mis. Calculate the test-section
Mach number. Assume the tunnel flow is adiabatic.
8.3 At a given point in a flow, T = 300 K, p = l.2 atm, and V = 250 mis. At
this point, calculate the corresponding values of p0 , To, p*, T *, and M *.
(8.67)
8.4 At a given point in a flow, T = 389 K, p = l.6 atm, and V = 900 mis. At
this point, calculate the corresponding values of p 0 , To, p*, T *, and M *.
8.5 Consider the isentropic flow through a supersonic nozzle. If the
test-section conditions are given by p = l atm, T = 230 K, and M = 2,
8.6
- Rln
[1 + ~(Mt-1)]
y+l
Po,2 =
Po.1
(8.68)
e -(s2-s,) / R
(8.73)
8.7
8.8
For a calorically perfect gas, the total temperature is constant across a normal
shock wave:
8.9
To,2 = To.,
However, there is a loss in total pressure across the wave:
8.10
For subsonic and supersonic compressible flow, the freestream Mach number
is determined by the ratio of Pi tot pressure to freestream static pressure. However, the equations are different:
Subsonic flow:
Mf
= _ 2_
[(Po.,) (y Y- 1
P1
2
(y + l) M t ]
4 yMf- 2(y - 1)
=[
l )/ y _
y/(y- l )
i]
8.12
8.13
(8.74)
l -y +2yMf
y+l
(8.80)
8.14
8.15
599
600
PART 3
8.16
8.17
8.18
8.19
e.
Calcu late the tempe rature at a stagnation point on the vehicl
the flow
in
tube
Pitot
a
,
In the test sectio n of a supersonic wind tunnel
a
(from
t
remen
reads a pressure of 1: 13 atm. A static pressure measu
late
Calcu
atm.
l
pressure tap on the sidewall o~ the test section) yields O.
n.
sectio
the Mach numb er of the flow m the test
the moo
When the Apollo comm and modu le returned to earth from
Using the r:~~;ts
entere d the earth' s atmos phere at a Mach numb er of 36.
the ratio of
with
gas
t
perfec
cally
calori
a
for
fro~ the presen t chapt er
stagnation
the
at
rature
tempe
gas
the
t
predic
sp~c1fic heats equal to 1.4,
ream
freest
the
where
e
altitud
an
at
36
Mach
at
pomt of the Apoll o
r.
answe
your
of
ty
validi
the
on
ent
Comm
temperature is 300 K.
as it entered
36
Mach
at
e
The stagnation tempe rature on the Apollo vehicl
ted in
predic
than
value
the atmosphere was 11 ,000 K, a much different
of
ratio
a
with
gas
t
perfec
Prob_lem 8.17 for the case of a calori cally
ons that
reacti
cal
chemi
to
spec1fi_c h~ats equal t~ 1.4. The differe nce is due
The
tion.
ioniza
and
on
occur m ~If at_these high tempe rature s-diss ociati
nt
consta
with
gas
t
anaJ~ses m this book assuming a caloricaJJy perfec
flows. However
specific h~ats ~e not valid for such chemically reacting
gas results are '
t
perfec
aJJy
as an ~ngm eenng appro ximat ion, the caloric
ic heats a
specif
of
ratio
the
somet imes appli~d with a lower value of
effects 'of hi h
the
te
simula
to
try
so-cal led "effective gamm a," in order to
in this g
stated
ion
condit
the
For
tempe rature chemi cally reacting flows.
yield a
to
ary
necess
a
gamm
proble m, calculate the value of the effective
eam
freestr
the
e
Assum
point.
tempe rature ?f 11,00 0 K at the stagnation
K.
300
1s
tempe rature
isentropic flow.
Prove that the total pressu re is constant throughout an
CH AP TE R
Obliq ue Sh ock an d
Expa nsi on Waves
shock wave is extrem ely
In the case of air (and the same is true for all gases) the
are still applic able for
thin so that calculations based on one-dimensional.flow
g through it, even
determining the chang es in velocity and density on passin
sion, provided that
dimen
one
to
d
limite
not
is
when the rest of the flow system
ered.
only the velocity compo nent norma l to the wave is consid
Pj EVIE W BOX
for hypers onic airplane:,. T he perform ance of such
Tai..e a look at Figure 9. 1. What you ,ee b the
e ngine~ depend ~ in part on the nature of the oblique
co 1puted wave pallem -both shock and expan~ion
wave pattern~ both upstrea m of and ins ide the en~ , es-ge nerated by a generi c super-, onic transp ort
gi ne. The material in thi !, c hapter is vital to you if
config uration flying at Mach 1.7 at an ahillld e of
you are interes ted in designing ,upers onic and hyperin
l.1 km. All the\c waves are obliqu e 10 the flow,
~onic wind tunne ls. where the oblique wave patte rn:,
co tra,1to the norma l shock wave\ discw, \ed in Chapcreated by mode b in the tunne l and in the diffuser
te r 8. The present c hapter is a ll about obliqu e shock
down~ tream of the model!> affect the perfor mance
and ex pansio n waves.
the tunne l. In fact, the materia l in this c hapter b
The material in thi~ chapte r i:. vita l to a funda- of
vital to a w hole ho:,t of applica tion, in ~uper-,onic
me ntal under~ tanding of \ uper~o ni c flow. Moreo ver,
flow.
it 1, vital to you if you a re inte rested in de~ign ing
By now you get the me~~ age-th e material in
accept
lly
nta
me
nviron
e
and
ble
fea,i
ically
an econom
thi s chapte r is simply \i tal to your 1-tudy of superahle supers onic transpo rt. The shod. wave~ in Fig1,onic flows. It is the bread and butler of such flow!..
u c 9.1 create the major ~ource of drag {wave drag) on
So lay out the bread, ~pread it with butler. and consume
the: a irplane . and the waves. when they propag ate to
thi ~ c hapter. And on top of everyth ing ebe. learning
the ground . cause the muc h discu~sed ~onic boom:
oblique :,hock and expans ion waves is excitin g.
1 he materi al in this chapte r i~ abo , ita l to you if about
thi ~.
Mjet engine!> I predict that you are going to e njoy
, l U are intere!>ted in design ing SCRA
601
706
PART 3
CHAPTER 1 O
10.7 SUMMARY
Continuity:
(10.l)
Momentum:
(10.5)
Energy:
(10.9)
dA
- = (M 2 A
du
u
1)-
(10.25)
tells us that
In a supersonic wind tunnel, the ratio of second throat area to first throat area
should be approximately
A1,2
Po.1
(10.39)
A1,1 Po.2
If A,, 2 is reduced much below this value, the diffuser will choke and the tunnel
will unstart.
10.8 PROBLEMS
10.1 The reservoir pressure and temperature for a convergent-divergent nozzle
are 5 atm and 289 K, respectively. The flow is expanded isentropically to
supersonic speed at the nozzle exit. If the exit-to-throat area ratio is 2.193,
calculate the following properties at the exit: Me, Pe, Te, Pe, Ue, Po,e , To,e
10.2 A flow is isentropically expanded to supersonic speeds in a
convergent-divergent nozzle. The reservoir and exit pressures are 1 and
0.3143 atm, respectively. What is the value of Ae I A*?
10.3 A Pitot tube inserted at the exit of a supersonic nozzle reads
8.92 x 104 N/m2 . If the reservoir pressure is 2.02 x 105 N/m2, calculate
the area ratio AeI A* of the nozzle.
10.4 For the nozzle flow given in Problem 10.1, the throat area is 25 cm2
Calculate the mass flow through the nozzle.
10.5 A closed-form expression for the mass flow through a choked nozzle is
.
PoA*
m= - -
./To
(
y-1 z)] Cy+ t)/Cy -t )
A ) 2- 1 [ 2
- - - 1+ - -M
( -A*
2
M
y +1
2
(10.32)
This tells us that the Mach number in a duct is governed by the ratio of local
duct area to the sonic throat area; moreover, for a given area ratio, there are
two values of Mach number that satisfy Equation (10.32)-a subsonic value
and a supersonic value.
!.
(-2-)
y+l
(y+ l)/(y- 1)
707
708
PART 3
CHAPTER 10
709