16.

121 ANALYTICAL SUBSONIC AERODYNAMICS, MASSACHUSETTS INSTITUTE OF TECHNOLOGY

Fundamentals of Fluid Mechanics

1 F UNDAMENTALS OF F LUID M ECHANICS

1.1 A SSUMPTIONS
1. Fluid is a continuum

2. Fluid is inviscid

3. Fluid is adiabatic

4. Fluid is a perfect gas

5. Fluid is a constant-density fluid

6. Discontinuities (shocks, waves, vortex sheets) are treated as separate and serve as boundaries
for continuous portions of the flow

1.2 N OTATION
0
p = pressure (static) V = control volume
0 0
ρ = density S = surface surrounding V
T = temperature (absolute) σ = impermeable body
Q = velocity vector of fluid particles n = normal directed into the fluid
Q = U i + V j + Wk R = gas constant
F = body force per unit mass c p = specific heat at constant pressure
F = ∇Ω c v = specific heat at constant volume
Ω = potential of the force field γ = c p /c v
Gravity field: F = −g k; Ω = −g z e = internal energy per unit mass
p
h = enthalpy per unit mass; h = e + ρ s = entropy per unit mass

1.3 C ONTINUITY E QUATION

∂ρ
+ ∇(ρQ) = 0
∂t
Dρ
+ ρ∇Q = 0
Dt
∂ρ
Ñ Ó
dV 0 + ρ(Qn)d s 0 = 0
V 0 ∂t S 0 +V
∂ρ
Ñ h i
+ ∇(ρQ) dV 0 = 0
V 0 ∂t

1
 1.4 C ONSERVATION OF M OMENTUM
DQ ∇p
=F −
Dt ρ
∂
Ñ Ó
(ρQ)dV 0 + ρQ(Qn)d s 0
X
Fi =
i V0 ∂t S 0 +V

1.5 C ONSERVATION OF T HERMODYNAMIC E NERGY

D h Q2 i ∇ · (pQ)
e+ =− + F ·Q
Dt 2 ρ
D h Q 2 i ∂p
ρ h+ = + ρF · Q
Dt 2 ∂t

1.6 E QUATION OF S TATE

p = RρT (thermally perfect gas)
c p , c v = constants (calorically perfect gas)

2 P RESSURE DISTRIBUTION AND COMPRESSIBILITY

2.1 A SSUMPTIONS
1. Steady flow

2. Inviscid fluid

3. No discontinuities (shocks)

4. Perfect gas

5. One-dimensional motion

6. Adiabatic flow

7. F ≡ 0

8. Isentropic

2.2 N OTATION
( )0 = stagnation conditions, Q = 0
( )∞ = free stream conditions, Q = u c = u ∞ c
( ) = conditions on body surface (airfoil)

Q = u 0 i + u 0 j + ωk

u 0 = u ∞ + γu

2
 2.3 E NERGY E QUATIONS
p
h=e+
ρ
h 1 i
d h + Q2 = 0
2
(Heat content plus kinetic energy is constant)

2.4 P ERFECT G AS R ELATIONS

p = ρRT
pV = RT
1
V≡
ρ

Can show, without effort:

ρV γ = constant
³ 1 ´γ
p = constant
ρ
p
a 2 = γ , a = speed of sound
ρ
r ³ ´
Q= 2c p T0 − T

hT i h ³ p ´ γ−1 i
γ
T0 − T = T0 1 − = T0 1 −
T0 p0
n h ³ p ´ γ−1 io 1
γ 2
Q = 2c p T0 1 −
p0

2.5 M ACH N UMBER

2 2c p (T0 − T ) 2c p (T0 − T )
Q
M2 = = p =
a 2
γρ γRT

2c p ³T
0
´ 2 ³ T0 ´
M2 = −1 = −1
γ(c p − c v ) T (γ − 1) T
T0 h γ − 1 2i
= 1+ M = β(γ, M )
T 2
γ
p 0 ³ T0 ´ γ−1 γ
= = β γ−1
p T
ρ0 ³ T0 ´ γ−1 1 1
= = β γ−1
ρ T

3
 2.6 OTHER USEFUL FORMS , EXPRESSIONS
Q 2 = 2c p (T0 − T )
p0
a 02 = γ = γRT0
ρ0
Q 2 2c p ³ T ´ 2 ³ T ´
= 1 − = 1 −
a 02 γR T0 γ−1 T0
T γ−1 Q 2³ ´
= 1−
T0 2 a0
γ
p h γ − 1 ³ Q ´2 i γ−1
= 1−
p0 2 a0
ρ h γ − 1 ³ Q ´2 i γ−1
1

= 1−
ρ0 2 a0
γ−1 2
a 2 = a 02 − Q
2

2.7 P RESSURE , VELOCITY RELATIONS IN ISENTROPIC FLOW

With some effort, one may show:
γ
p h γ−1 2 ³ Q 2 ´i γ−1
= 1+ M∞ 1 − 2
p∞ 2 u∞

Expanding the right-hand side:

p γ³ Q2 ´ 2 γ ³ Q 2 ´2 4 γ(2 − γ) ³ Q 2 ´3 6 γ(2 − γ)(3 − 2γ) ³ Q 2 ´4 8

= 1+ 1− 2 M ∞ + 1− 2 M ∞ + 1− 2 M ∞ + 1− 2 M ∞ +...
p∞ 2 u∞ 8 u∞ 48 u∞ 384 u∞

Obtain an expression for

p − p∞
cp = 1
2 ρ ∞ u∞
2

Let
γV
Q = u ∞ + γV, ¿1
U∞
Find c p and discuss its limitations.

3 S IMILARITY OF FLOWS
3.1 R EQUIREMENTS FOR SIMILARITY OF FLOWS
1. Similarity in boundary geometry
Boundary of one flow can be made to coincide with that of another if its linear dimensions are
multiplied by a constant

2. Dynamic constraint
Dependent variables of one flow are proportional to those of another at the corresponding
points.

4
Example Problem - Illustration

Consider the dynamics of an incompressible fluid flow with constant.

Equation of incompressibility:
D p ∂ρ ∂ρ
= + ui =0
Dt ∂t ∂x i
Equation of continuity:
∂u i
=0
∂x i
Introduce dimensionless variables:

0 ui 0 ρ 0 p 0 xi 0 tU
ui = , ρ = , p = , xi = , t =
U ρ0 ρ 0U 2 L L

U , ρ 0 , L − reference quantities

3.2 L INEAR MOMENTUM

0 0
0
³∂ 0 ∂ ´ 0 ∂p ρL γ ∂2 0
ρ + uα 0 ui = − 0 + 2 Fi + 0 ui
∂t
0
∂x α ∂x U U L ∂x α ∂x a
i
0 0 0
∂ρ 0 ∂ρ ∂u α
+ uα =0 =0
∂t
0 0 0
∂x α ∂x a
Froude no: F = pU −→ inertia forces
gravity force
gL
Reynolds no: Re = UγL −→ inertia force
viscous force

F and Re must be the same for both flows. This is sufficient for dynamic similarity along with similar
boundary geometry.

U , ρ 0 , L may be different for both flows.

4 E QUATIONS GOVERNING IRROTATIONAL FLOWS OF A HOMENTROPIC GAS

For this class of flows the simplification is through the introduction of the velocity potential, φ, where

Q = ∇φ

or
∂φ
ui =
∂x i
and the vorticity is zero: ω = ∇ ×Q = ∇ × ∇φ = 0 where ω is the vorticity vector.

The unsteady Bernoulli equation may be written, for this class of flows:

∂Q ³1 ´ 1
+ ∇ Q 2 − Qxω = − ∇p
∂t 2 ρ

5
since, p = p(ρ), ω=0
∂Q ³1 ´ 1
+ ∇ Q 2 + ∇p = 0
∂t 2 ρ
or ³ ∂φ 1 ∂p ´
Z
∇ + Q2 + =0
∂t 2 ρ
therefore
∂φ 1 2 dp
Z
+ Q + = f (t )
∂t 2 ρ
Absorb f (t ) into φ and obtain
∂φ 1 2 dp
Z
+ Q + = constant
∂t 2 ρ
Differentiate above equation with respect to time, t :

∂2 φ ∂Q 1 ∂ρ
+Q · + a2 =0
∂t 2 ∂t ρ ∂t

Expressing the continuity equation in terms of φ:

1 ∂p 1
+ ∇2 φ + Q · ∇ρ = 0
ρ ∂t ρ

Linear momentum equation rewritten yields

1 1 1 1 n ∂Q o
Q· ∇ρ = 2 Q · ∇p = 2 Q − − (Q · ∇)Q
ρ a ρ a ∂t

Combining the above three equations yields:

1 ∂2 Φ 2 ∂Q 2 1 h i
+ Q · = ∇ φ − Q · (Q · ∇)Q
a 2 ∂t 2 a 2 ∂t a2
∂Φ
since u i = ∂x i , the above equation may be written:

³ u 2 ´ ∂2 Φ ³ v 2 ´ ∂2 Φ ³ w 2 ´ ∂2 Φ uv ∂2 Φ v w ∂2 Φ uw ∂2 Φ
∗ 1− 2 + 1 − + 1 − − 2 − 2 − 2 =
a ∂x 2 a 2 ∂y 2 a 2 ∂z 2 a 2 ∂x∂y a 2 ∂y∂z a 2 ∂x∂z

1 ∂2 φ u ∂2 φ v ∂2 φ w ∂2 φ
+ 2 + 2 + 2
a 2 ∂t 2 a 2 ∂x∂t a 2 ∂y∂t a 2 ∂z∂t
where
∂Φ ∂Φ ∂Φ
u= v= w=
∂x ∂y ∂z
For steady flow of a calorically perfect gas:

h 0 = constant

Q2
cp T +
= c p T0
2
γ−1 2
a 2 = a 02 − (Φx + Φ2y + Φ2z )
2
Equation * is the potential-flow equation.

6
 5 S MALL P ERTURBATION T HEORY
We will consider a slender body immersed in a uniform flow, viz.,

in the uniform flow:

Q = U∞ i
in the perturbed flow:
Q = ui + v j + w k
0 0 0
Q = (U∞ + u )i + v + w
j k

Q = ∇φ
Now define a perturbation velocity potential, φ(x, y, z), where

0 ∂φ
u =
∂x
0 ∂φ
v =
∂y
0 ∂φ
w =
∂z
∴ Φ(x, y, z) = U∞ x + φ(x, y, z)
Using the notation in eqn(*):
0 ∂Φ ∂φ
u = U∞ + w = = U∞ +
∂x ∂x
0 ∂Φ ∂φ
v =v = =
∂y ∂y
0 ∂Φ ∂φ
w =w = =
∂z ∂z
∂2 φ
Φxx = = φxx
∂x 2
∂2 φ
Φy y = = φy y
∂y 2
∂2 φ
Φzz = = φzz
∂z 2
∂2 φ
Φx y = = φx y
∂x∂y
∂2 φ
Φy z = = φy z
∂y∂z

7
 ∂2 φ
Φxz = = φxz
∂x∂z
Substituting Φ = U∞ x + φ and multiplying eqn(*) by a 2 we obtain the perturbation equation or per-
turbation velocity potential equation, for steady flow:
∗∗[a 2 −(U∞ +φx )2 ]φxx +[a 2 −(φ y )2 ]φ y y +[a 2 −(φz )2 ]φzz −2(U∞ +φx )φ y φx y −2(U∞ +φx )φz φxz −2φ y φz φ y z = 0
Note a 2 may be expressed as:
γ−1 0 0 0 0
a 2 = a∞
2
− (2u U∞ + u 2 + v 2 + w 2 )
2
γ − 1
a 2 = a∞2
− (2φx U∞ + (φx )2 + (φ y )2 + (φz )2 )
2
Also, note that eqn(**) is exact! It is also non-linear.

5.1 P ERTURBATIONS
Assume the perturbations are small, viz.,
0 0 0
u v w
¿ 1; ¿ 1; ¿1
U∞ U∞ U∞
In the limit of small perturbations, we may neglect the terms containing squares of the perturbation
velocities in comparison to those containing first powers. Eqn(**) with a 2 substituted becomes
2 2 φx 2 φx 2
φy 2 φz
(1−M ∞ )φxx +φ y y +φzz = M ∞ (γ+1) φxx +M ∞ (γ−1) (φ y y +φzz )+2M ∞ φx y +2M ∞ φxz
U∞ U∞ U∞ U∞
Note that each term on the right-hand side is non-linear. Each term on the right-hand side contains
a perturbation velocity (φx , φ y , or φz ). Hence, we may neglect the right-hand side in comparison to
the left-hand side. We obtain
2
(1 − M ∞ )φxx + φ y y + φzz = 0

5.2 QUESTIONS
1. What is the equation where M ∞ → 1?

2. What is the equation where M ∞ À 1?

6 B OUNDARY C ONDITIONS
1. The body surface is a stream line. (inviscid, irrotational flows)

2. Flow velocity must be tangent to body surface

3. Velocity vector has to be orthogonal to the unit normal of the body surface

8
The body surface is described by f (x, y, z)
f (x, y, z) = 0
Boundary condition is expressed as
Q · ∇ f (x, y, z) = 0
or
∂f
ui =0
∂x i
Introducing the perturbation velocities
0
u = U∞ + u
0
v =v
0
w =w
Substituting,
0 ∂f ∂f 0 ∂f
(U∞ + u ) +v +w =0
∂x ∂y ∂z
0
Since u ¿ U∞ , we may write:
∂f 0 ∂f 0 ∂f
U∞ +v +w =0
∂x ∂y ∂z
This equation must be satisfied on the surface of the body. Consider the two-dimensional case:
0
w =0
∂f
=0
∂z
We obtain: 0
v ∂ f /∂x d y
=− =
U∞ ∂ f /∂y d x
u −1
Therefore U∞ is the slope of the body (approximately) the slope of the streamline. Recall that
0 ∂φ dy ´
u = = U∞
∂y d x BODY
0
Now for thin bodies, a small angle of attack, y BODY ≈ 0: this suggests an expansion of v (x, y) in a
powers of y:
0 0
³ ∂v 0 ´
v (x, y) = v (x, 0) + y + ...
∂y y=0
0
³d y ´
∴ v(x, y) = v (x, 0) ∼ = U∞
d x BODY
For three-dimensional planar flows
∂f ∼
=0
∂z
and the boundary condition becomes
0
³ ∂y ´
v (x, 0, z) = U∞
∂x BODY
at infinity:
0
u →0
0
v →0
0
w →0
0 0 0
or w , v , and w are finite.

9
 7 L INEARIZED P RESSURE C OEFFICIENT
Let’s revisit the pressure coefficient, c p :
p − p∞
cp ≡ 1
2 ρ ∞U∞
2

where p is the pressure (local, static) at the location or point of interest in the flow field. Note that c p
is dimensionless.
Since,
1 1 γp ∞ γ U2 γ
ρ ∞U∞ 2
= ρ ∞U∞2
= p ∞ 2∞ = p ∞ M ∞ 2
2 2 γp ∞ 2 a∞ 2
then
2 h p i
cp = − 1
γM ∞
2 p
∞
for an inviscid, adiabatic, isentropic, steady flow and
0 0 0
Q = (U∞ + u )i + v j + w k

We show that
1 1 2
h + Q 2 = h ∞ + U∞
2 2
which for a calorically perfect gas leads directly to
2
T γ − 1 U∞ −Q2
= 2
T∞ 2 a∞
γ−1 0 0 0 0
= 1 − 2 [2u U∞ + u 2 + v 2 + w 2 ]
2a ∞
Isentropic flow conditions lead to:
p h T i γ
γ−1
=
p∞ T∞
γ
p h γ − 1 2 ³ 2u 0 u 02 + v 02 + w 02 ´i γ−1
= 1− M∞ + 2
p∞ 2 U∞ U∞
In the case of small velocity perturbations,
u0
¿1
U∞
³ u 0 ´2
n1
U∞
³ v 0 ´2
n1
U∞
³ w 0 ´2
n1
U∞
Using the binomial expansion, we show that

p γ 2 ³ u0 u 02 + v 02 + w 02 ´
= 1 − M∞ 2 + 2
+ ...
p∞ 2 U∞ U∞
therefore:
2u 0
cp = −
U∞
Discuss the limitations implied in the above expression for c p .

10
 8 C ROCCO ’ S T HEOREM
Consider the motion of a fluid element. The fluid element may both translate and rotate.
Let:
v = translational velocity
w = rotational velocity
ω = angular velocity
where
1
w = ∇×v
2
∇ × v ≡ vorticity
Combine Euler’s equation, first and second laws of thermodynamics:

∂v
ρ + ρ(v · ∇)∇ = −∇p
∂t
∇p
T ∇s = ∇h − v∇p = ∇h −
ρ
v2
h = h0 −
2
We obtain:
∂v n
T ∇s = ∇h 0 − v × (∇ × v) + Crocco’s Theorem
∂t
For steady flow, we obtain
T ∇s = ∇h 0 − v × (∇ × v)
or
v × (∇ × v) = ∇h 0 − T ∇s
For two-dimensional, steady flows:
1 ³ ∂s ∂h 0 ´
2w = T −
v ∂n ∂n
Vorticity =⇒ rates of change of entropy and stagnation enthalpy normal to the streamlines

Flow over a supersonic blunt body:

11
For this flow,
h 0 = constant
∂h 0
=0
h
∂s
6= 0 (why?)
∂n

12
MIT OpenCourseWare
https://ocw.mit.edu/

16.121 Analytical Subsonic Aerodynamics

Fall 2017

For information about citing these materials or our Terms of Use, visit: https://ocw.mit.edu/terms.