Chapter 5 Mass, Momentum, and Energy Equations

57:020 Mechanics of Fluids and Transport Processes
Professor Fred Stern Fall 2006
Chapter 5
Chapter 5 Mass, Momentum, and Energy Equations

Flow Rate and Conservation of Mass
1. cross-sectional area oriented normal to velocity vector
(simple case where V A)
U = constant: Q = volume flux = UA [m/s m2 = m3/s]

U constant: Q = A UdA
Similarly the mass flux = m A UdA
2. general case

Chapter 5
Q V ndA
CS
V cos dA
CS
V n dA
m
CS
average velocity:
Q
A
Example:
At low velocities the flow through a long circular tube, i.e. pipe,
has a parabolic velocity distribution (actually paraboloid of
revolution).
i.e., centerline velocity
a) find Q and
Q V ndA udA
A
2 R
udA u (r )rddr
0 0
2 u ( r ) rdr
dA = 2rdr
r
u u max 1

u = u(r) and not

Q=
r
2 u max 1
R
0
1
V u max
2
d 2
Chapter 5
rdr
1
u max R 2
2
Continuity Equation
RTT can be used to obtain an integral relationship expressing
conservation of mass by defining the extensive property B = M
such that = 1.
B = M = mass
= dB/dM = 1
General Form of Continuity Equation
dM
d
0
dV V dA
dt
dt CV
CS
or
V dA
CS
net rate of outflow

of mass across CS
d
dV
dt CV
rate of decrease of
mass within CV
Simplifications:
d
1. Steady flow: dt CV dV 0
2. V = constant over discrete dA (flow sections):

Chapter 5
V dA V A
CS
CS
3. Incompressible fluid ( = constant)

d
V
dA
CS
dV
dt CV
conservation of volume
4. Steady One-Dimensional Flow in a Conduit:

V A 0
CS
1V1A1 + 2V2A2 = 0
for = constant
Q1 = Q 2
Some useful definitions:

Mass flux
V dA
m
Volume flux
Q V dA
Average Velocity
V Q/A
Average Density
1
dA
A
& Q unless = constant

Note: m
Example

Chapter 5
*Steady flow
*V1,2,3 = 50 fps
0
V dA 0
CS
*@ V varies linearly
from zero at wall to
Vmax at pipe center
*find m 4 , Q4, Vmax
*water, w = 1.94 slug/ft3
d
dV
dt CV
4
m
i.e., -1V1A1 - 2V2A2 + 3V3A3 + A V4 dA 4 = 0

= const. = 1.94 lb-s2 /ft4 = 1.94 slug/ft3
4
4 V4 dA 4
m
= V(A1 + A2 A3)
V1=V2=V3=V=50f/s
= 144 50 4 12 2 2 1.52
= 1.45 slugs/s
1.94
Q4 =
=
4 .75
m
ft3/s
V4 dA 4
A4
Q4 =
ro 2
0 0
Vmax
velocity profile
r
1 rddr
ro
dA4

Chapter 5
V4 V4()
ro
r
2 Vmax 1
ro
0
rdr
1 2
r V
Q 3 o max
V4
A
ro2
ro
r2
2 Vmax r dr
ro
0
r2
2 Vmax
2
r0
r
3ro
Vmax =
1
ro 2
3
1
1 1
ro 2 Vmax
3
2 3
2 Vmax ro 2
Q4
3 ro
2.86
fps
1
Vmax
3

Chapter 5
Momentum Equation
RTT with B = MV and = V
FS F B
d
VdV VV R dA
dt CV
CS
V = velocity referenced to an inertial frame (non-accelerating)

VR = velocity referenced to control volume
FS = surface forces + reaction forces (due to pressure and
viscous normal and shear stresses)
FB = body force (due to gravity)
Applications of the Momentum Equation
Initial Setup and Signs
1. Jet deflected by a plate or a vane
2. Flow through a nozzle
3. Forces on bends
4. Problems involving non-uniform velocity distribution
5. Motion of a rocket
6. Force on rectangular sluice gate
7. Water hammer
General form for
moving but
dBsys d
dV V R dA
Recall RTT:
non-accelerating
dt
dt CV
CS
reference frame
VR=velocity relative to CS=V VS=absolute velocity CS
Derivation of the Basic Equation
Subscript not shown in text but implied!
i.e., referenced to CV
Let,
B = MV = linear momentum
=V
V must be referenced to
inertial reference frame
d(M V)
d
F
VdV VV R dA
dt
dt CV
CS
Must be relative to a non-accelerating

inertial reference frame

Chapter 5
Newtons 2nd law
where
F = vector sum of all forces acting on the

control volume including both surface and
body forces
= FS + FB
FS = sum of all external surface forces acting at
the CS, i.e., pressure forces, forces
transmitted through solids, shear forces, etc.
FB = sum of all external
body forces, i.e.,
gravity force
Fx = p1A1 p2A2 + Rx
Fy = -W + Ry
free body diagram
R = resultant force on fluid

in CV due to pw and w
i.e., reaction force on fluid
Important Features (to be remembered)

1) Vector equation to get component in any direction must use
dot product
x equation
Fx
d
udV u V R dA
dt CV
CS
y equation
Fy
d
vdV v V R dA
dt CV
CS
Carefully define coordinate

system with forces positive in
positive direction of coordinate
axes

Chapter 5
z equation
Fz
d
wdV w V R dA
dt CV
CS
2) Carefully define control volume and be sure to include all

external body and surface faces acting on it.
For example,
(Rx,Ry) = reaction
force on fluid
(Rx,Ry) = reaction
force on nozzle
3) Velocity V must be referenced to a non-accelerating inertial

reference frame. Sometimes it is advantageous to use a
moving (at constant velocity) reference frame. Note VR = V
Vs is always relative to CS.
4) Steady vs. Unsteady Flow
Steady flow
d
VdV 0
dt CV
i.e., in these cases V used for

B also referenced to CV (i.e.,
V = VR)

Chapter 5
10
5) Uniform vs. Nonuniform Flow

VV R dA
CS
= change in flow of momentum across CS

= VVRA
uniform flow across A
6) Fpres = p ndA
fdV f nds
f = constant, f = 0
= 0 for p = constant and for a closed surface

i.e., always use gage pressure
7) Pressure condition at a jet exit
at an exit into the atmosphere jet
pressure must be pa
Application of the Momentum Equation

1. Jet deflected by a plate or vane
Consider a jet of water turned through a horizontal angle
CV and CS are for

jet so that Fx and Fy
are vane reactions
forces on fluid

x-equation:
Fx Fx
Chapter 5
11
d
udV u V dA
dt
CS
Fx u V A
CS
steady flow
= V1x (V1A1 ) V2 x (V2 A 2 )

continuity equation:
A1V1 = A2V2 = Q
for A1 = A2
V1 = V2
Fx = Q(V2x V1x)
y-equation:
Fy Fy v V A
CS
Fy = V1y( A1V1) + V2y( A2V2)

= Q(V2y V1y)
for above geometry only
where:
note:
V1x = V1 V2x = -V2cos V2y = -V2sin V1y = 0

Fx and Fy are force on fluid
- Fx and -Fy are force on vane due to fluid
If the vane is moving with velocity Vv, then it is convenient to

choose CV moving with the vane
i.e., VR = V - Vv and V used for B also moving with vane
x-equation:
Fx u V R dA
CS
Fx = V1x[-(V Vv)1A1] + V2x[-(V Vv)2A2]

Chapter 5
12
Continuity: 0 = V R dA
i.e., (V-Vv)1A1 = (V-Vv)2A2 = (V-Vv)A
Qrel
Fx = (V-Vv)A[V2x V1x]
Qrel
on fluid
V2x = (V Vv)2x
V1x = (V Vv)1x
Power = -FxVv
For coordinate system

moving with vane
i.e., = 0 for Vv = 0
Fy = Qrel(V2y V1y)
2.
Flow through a nozzle
Consider a nozzle at the end of a pipe (or hose). What force is

required to hold the nozzle in place?
CV = nozzle
and fluid
(Rx, Ry) =
force required to
hold nozzle in
known. place
Then, the
Assume either the pipe velocity or pressure is

unknown (velocity or pressure) and the exit velocity V 2 can be
obtained from combined use of the continuity and Bernoulli
equations.
Bernoulli:
1
1
p1 z1 V12 p 2 z 2 V22
2
2
z1=z2

Chapter 5
13
1
1
p1 V12 V22
2
2
Continuity:
A1V1 = A2V2 = Q
V2
A1
D
V1
A2
d
V1
1
D
p1 V12 1
2
d
Say p1 known:
1/ 2
2p1
V1
4
1 D d

To obtain the reaction force Rx apply momentum equation in xdirection

Fx
d
udV u V dA
dt CV
CS
u V A
= CS
steady flow and uniform

flow over CS
Rx + p1A1 p2A2 = V1(-V1A1) + V2(V2A2)
= Q(V2 - V1)
Rx = Q(V2 - V1) - p1A1
To obtain the reaction force Ry apply momentum equation in ydirection
Fy v V A 0
since no flow in y-direction

Ry Wf WN = 0 i.e., Ry = Wf + WN
CS
S under
Numerical Example: Oil with S = .85 flows in pipe
1.65
g
pressure of 100 psi. Pipe diameter is 3 and nozzle tip diameter
D/d = 3
is 1
1
Q = 4 12 V
= .716 ft3/s
2

Chapter 5
14
V1 = 14.59 ft/s
V2 = 131.3 ft/s
Rx = 141.48 706.86 = 569 lbf
Rz = 10 lbf
This is force on nozzle
3. Forces on Bends
Consider the flow through a bend in a pipe. The flow is
considered steady and uniform across the inlet and outlet
sections. Of primary concern is the force required to hold the
bend in place, i.e., the reaction forces Rx and Ry which can be
determined by application of the momentum equation.

Chapter 5
15
Rx, Ry = reaction force on

bend i.e., force
required to hold
bend in place
0 V A V1A1 V2 A 2
Continuity:
i.e., Q = constant =
V1A1 V2 A 2
Fx u V A
p1A1 p 2 A 2 cos R x V1x V1A1 V2 x V2 A 2
x-momentum:
= Q V2 x V1x
y-momentum:
Fy v V A
p 2 A 2 sin R y w f w b V1y V1A1 V2 y V2 A 2
= Q V2 y V1y
4. Problems involving Nonuniform Velocity Distribution

See text pp. 215 216

Chapter 5
16
5. Force on a rectangular sluice gate

The force on the fluid due to the gate is calculated from the xmomentum equation:
Fx u V A
F1 FGW Fvisc F2 V1 V1A1 V2 V2 A 2

FGW F2 F1 Q V2 V1 Fvisc
usually can be neglected
y2
y
y 2 b 1 y1b Q V2 V1
2
2
1
FGW b y 22 y12 Q V2 V1
2
Q 2 1
1

b y 2 y1
Moment of Momentum Equation

See text pp. 221 229
Energy Equations
V1
Q
y1b
V2
Q
y2b

Chapter 5
17
Derivation of the Energy Equation

The First Law of Thermodynamics
The difference between the heat added to a system and the work
done by a system depends only on the initial and final states of
the system; that is, depends only on the change in energy E:
principle of conservation of energy
E = Q W
E = change in energy
Q = heat added to the system
W = work done by the system
E = Eu + Ek + Ep = total energy of the system
potential energy
kinetic energy
Internal energy due to molecular motion
The differential form of the first law of thermodynamics

expresses the rate of change of E with respect to time
dE
QW
dt
rate of work being done by system

rate of heat transfer to system

Chapter 5
18
Energy Equation for Fluid Flow

The energy equation for fluid flow is derived from Reynolds
transport theorem with
Bsystem = E = total energy of the system (extensive property)
= E/mass = e = energy per unit mass (intensive property)
= u + ek + ep
dE
d
edV CS e V dA
dt
dt CV
d
Q& W& (u ek e p )dV (u ek e p )V dA
CS
dt CV
This can be put in a more useable form by noting the following:

Total KE of mass with velocity V MV 2 / 2 V 2
mass
M
2
Ep
Vz
ep
gz
(for Ep due to
M V
ek
V2 V
gravity only)
V2
V2
d
&
&
Q W
gz u dV
gz u V dA
Cs
dt CV 2
2
rate of work
done by system
rate of change
of energy in CV
rate of heat
transfer to sysem
Rate of Work Components:
W
s W
f
W
flux of energy
out of CV
(ie, across CS)

Chapter 5
19
For convenience of analysis, work is divided into shaft work Ws

and flow work Wf
Wf = net work done on the surroundings as a result of
normal and tangential stresses acting at the control
surfaces
= Wf pressure + Wf shear
Ws = any other work transferred to the surroundings
usually in the form of a shaft which either takes
energy out of the system (turbine) or puts energy into
the system (pump)
Flow work due to pressure forces Wf p (for system)
Note: here V uniform over A
CSCV
Work = force distance

at 2 W2 = p2A2 V2t (on surroundings)
System at time t + t
rate of work W 2 p 2 A 2 V2 p 2 V 2 A 2
neg. sign since pressure

at 1 W1 = p1A1 V1t
System at time t
1 p1 V 1 A 1
W
force on surrounding fluid
actsIn
ingeneral,
a direction opposite
to the motion of the system
boundary
fp p V A
W
for more than one control surface and V not necessarily uniform
over A:

Chapter 5
20
fp p V dA p V dA
W
CS
CS

f W
fp W
fshear
W
Basic form of energy equation
p
Q& W&s W&fshear V dA
CS

V2
V2
d

gz u dV
gz u V dA
CS
dt CV 2
2
Q& W&s W&fshear
V2
d

gz u dV
dt CV 2
Usually this term can be eliminated

by proper choice of CV, i.e.
CS normal to flow lines.
Also, at fixed boundaries the
velocity is zero (no slip
condition) and no shear
stress flow work is done.
Not included or discussed in
text!
V2
p

gz u V dA
CS
2
h=enthalpy
Simplified Forms of the Energy Equation

Energy Equation for Steady One-Dimensional Pipe Flow
Consider flow through the pipe system as shown

Chapter 5
21
Energy Equation (steady flow)

V2
p
&
&
Q Ws
gz u V dA
CS

2
3
V
1
1
1
Q& W&s
gz1 u1 1V1 A1
dA1
A1
A
1
2

p2
2V23

gz2 u2 2V2 A2
dA2
A2
A2
*Although the velocity varies across the flow sections the

streamlines are assumed to be straight and parallel;
consequently, there is no acceleration normal to the streamlines
and the pressure is hydrostatically distributed, i.e., p/ +gz =
constant.
*Furthermore, the internal energy u can be considered as
constant across the flow sections, i.e. T = constant. These
quantities can then be taken outside the integral sign to yield
p1
V13
&
&
Q Ws
gz1 u1 V1dA1
dA1
A
A
1
1
2

p2
V23
gz2 u2 V2 dA2
dA2
A
A
2
2
2


22
Q V dA VA
Recall that
So that
V dA VA m
Define:
mass flow rate
V3
V A
V
dA
m
2
2
2
A
K.E. flux
3
i.e.,
Chapter 5
1 V
dA
A A V
K.E. flux for V= V =constant across pipe
= kinetic energy correction factor
2
p
p
22
V
V
1
Q& W& 1 gz1 u1 1 m& 2 gz2 u2 2 m&

1 & & p1
V 1 p2
V2
Q W gz1 u1 1
gz2 u2 2
m&
= 1 if V is constant across the flow section

> 1 if V is nonuniform
Nnote that:
laminar flow = 2
turbulent flow = 1.05 1 may be used
Shaft Work
Shaft work is usually the result of a turbine or a pump in the
flow system. When a fluid passes through a turbine, the fluid is
doing shaft work on the surroundings; on the other hand, a pump
does work on the fluid
W
W
t and W
W
where W
are
work
magnitudes of power
time
Using this result in the energy equation and deviding by g
results in
s

Chapter 5
23
W&p p1
V12 W&t p2
V22 u2 u1 Q&
z1 1
z2 2
&
&
&
mg
2 mg
2
g
mg
mechanical part
Note: each term has dimensions of length
Define the following:
hp
ht
p
W
g
m
p
W
Qg
p
W
Q
t
W
g
m
u2 u1 Q&
hL
head loss
&
g
mg
thermal part

Chapter 5
24
Head Loss
In a general fluid system a certain amount of mechanical energy
is converted to thermal energy due to viscous action. This effect
results in an increase in the fluid internal energy. Also, some
heat will be generated through energy dissipation and be lost
(i.e. - Q
). Therefore the term
from 2nd law
u2 u1 Q&
hL
0
g
gm&
represents a loss in
mechanical energy due
to viscous stresses
Note that adding Q

to system will not make hL = 0 since this
also increases u. It can be shown from 2nd law of
thermodynamics that hL > 0.
Drop over V and understand that V in energy equation refers

to average velocity.
Using the above definitions in the energy equation results in
(steady 1-D incompressible flow)
p1
V2
p
V2
1 1 z1 h p 2 2 2 z 2 h t h L
2g
2g
form of energy equation used for this course!

Chapter 5
25
Comparison of Energy Equation and Bernoulli Equation

Apply energy equation to a stream tube without any shaft work
Infinitesimal stream tube
Energy eq :
1=2=1
p1 V12
p V2
z1 2 2 z 2 h L
2g
2g
If hL = 0 (i.e., = 0) we get Bernoulli equation and

conservation of mechanical energy along a streamline
Therefore, energy equation for steady 1-D pipe flow can be
interpreted as a modified Bernoulli equation to include viscous
effects (hL) and shaft work (hp or ht)
Summary of the Energy Equation
The energy equation is derived from RTT with
B = E = total energy of the system
= e = E/M = energy per unit mass
= u +
1 2
V +gz
2

internal
KE
Chapter 5
26
PE
dE d
W
edV eV dA Q
dt dt CV
CS
from 1st Law of

Thermodynamics
heat add work

p
V
z
ed done
2g
1
2
1
Neglected in text presentation

W
s W
p W
v
W
shaft work
done on or by pressure work Viscous stress
system (pump done on CS work on CS
or turbine)
p pV dA p V dA
W
CV
CS
s W
t W
p
W
W
t W
p d edV e p e V dA
Q
dt CV
CS
1
e u V 2 gz
2
For steady 1-D pipe flow (one inlet and one outlet):
1) Streamlines are straight and parallel
p/ +gz = constant across CS
2) T = constant u = constant across CS
3)
define
1 V
dA
A CS V
= KE correction factor

Chapter 5
27
V
V
3
A
m
V dA
2
2
2
mechanical energy
p1
V2
p
V2
1 1 z1 h p 2 2 2 z 2 h t h L
2g
2g
p m
g
hp W
t m
g
ht W
u2 u1 Q&
hL
head loss
&
g
mg
Thermal
energy
Note: each term has
units of length
V is average velocity
(vector dropped) and
corrected by
> 0 represents loss in mechanical energy due to viscosity

Chapter 5
28
Concept of Hydraulic and Energy Grade Lines

p1
V2
p
V2
1 1 z1 h p 2 2 2 z 2 h t h L
2g
2g
Define
HGL =
p
z
p
V2
z
2g
point-by-point
application is
EGL =
graphically
displayed
HGL corresponds to pressure tap measurement + z
EGL corresponds to stagnation tube measurement + z
EGL = HGL if V = 0
EGL1 = EGL2 + hL
for hp = ht = 0
hL =
L V2
D 2g
i.e., linear variation in L for D,

V, and f constant
f = friction factor
f = f(Re)
p2
h
p2
V22
tube: 2g h
pressure tap:
stagnation
h = height of fluid in
tap/tube
EGL1 + hp = EGL2 + ht + hL
EGL2 = EGL1 + hp ht hL
abrupt change
due to hp or ht
L V2
D 2g
Helpful hints for drawing HGL and EGL

1. EGL = HGL + V2/2g = HGL for V = 0

2.&3.
p1
z1
2
V1
2g
p2
L V2
D 2g
in pipe means EGL and HGL will slope

downward, except for abrupt changes due to ht or hp
hL f
z2
2
V2
2g
hL
HGL2 = EGL1 - hL
hL
2g
Chapter 5
for abrupt expansion
29

Chapter 5
30
4. p = 0 HGL = z
5. for
hL f
L V2
D 2g
= constant L
i.e., linearly increased for increasing
f V2
L with slope D 2g
EGL/HGL slope downward

6. for change in D change in V
i.e.
V1A1 = V2A2
D 2
D 22
V1 1 V2
4
4
V1D12 V1D 22
change in distance between

HGL & EGL and slope
change due to change in hL

7. If HGL < z then p/ < 0
Chapter 5
31
i.e., cavitation possible
condition for cavitation:

p p va 2000
gage pressure
N
m2
p va ,g p A p atm p atm 100,000

p va ,g
N
m2
10 m
9810 N/m3

Chapter 5
32

Chapter 5
33

Chapter 5
34
Application of the Energy, Momentum, and

Continuity Equations in Combination
In general, when solving fluid mechanics problems, one should
use all available equations in order to derive as much
information as possible about the flow. For example, consistent
with the approximation of the energy equation we can also apply
the momentum and continuity equations
Energy:
p1
V2
p
V2
1 1 z1 h p 2 2 2 z 2 h t h L
2g
2g
Momentum:
2
2
Fs V2 A 2 V1 A1 Q V2 V1
Continuity:
A1V1 = A2V2 = Q = constant
one inlet and

one outlet
= constant

Chapter 5
35
Abrupt Expansion
Consider the flow from a small pipe to a larger pipe. Would like
to know hL = hL(V1,V2). Analytic solution to exact problem is
extremely difficult due
to the occurrence of
flow separations and
turbulence. However, if
the assumption is made
that the pressure in the
separation region
remains approximately
constant and at the
value at the point of
separation, i.e, p1, an approximate solution for hL is possible:
Apply Energy Eq from 1-2 (1 = 2 = 1)
p1
V2 p
V2
z1 1 2 z 2 2 h L
2g
2g
Momentum eq. For CV shown (shear stress neglected)

Fs p1A 2 p 2 A 2 W sin u V A
A 2 L
z
L
= V1 (V1A1 ) V2 (V2 A 2 )
= V22 A 2 V12 A1
W sin
next divide momentum equation by A2
A2
p1 p 2
V2 V2 A V2 A A
z1 z 2 2 1 1 1 1 1 1

g
g A2 g A2 A2
from energy equation

Chapter 5
36
V22
V12
V22 V12 A1
hL
2g 2g
g
g A2
V22 V12
2A
1 1
hL
2g 2g
A2
hL
continutity eq.
V1A1 = V2A2
A
1 2
V2 V12 2V12 1
2g
A2
2V1V2
hL
1
V2 V1 2
2g
If V2 = V1 ,
1 2
hL =
V1
2g
A1 V2
A 2 V1

Chapter 5
37
Forces on Transitions
Example 7-6
Q = .707 m3/s
V22
.
1
head loss = 2g
(empirical equation)
Fluid = water
p1 = 250 kPa
D1 = 30 cm
D2 = 20 cm
Fx = ?
First apply momentum theorem

Fx u V A
Fx + p1A1 p2A2 = V1(V1A1) + V2(V2A2)

Fx = Q(V2 V1) p1A1 + p2A2
force required to hold transition in place

Chapter 5
38
The only unknown in this equation is p2, which can be obtained

from the energy equation.
p1 V12 p 2 V22
hL
2g
2g
V22 V12
p 2 p1
hL
2 g 2g
note: z1 = z2 and = 1
drop in pressure
V22 V12
Q
V
A
p
h
2
1
2 1
L
x
2g 2g
p1A1
p2
In this equation,
(note: if p2 = 0 same as nozzle)
continuity
V1 = Q/A1 = 10 m/s
V2 = Q/A2 = 22.5 m/s
A1V1 = A2V2
V2
A1
V1
A2
i.e. V2 > V1
V22
h L .1
2.58m
2g
Fx = 8.15 kN
is negative x direction to hold

transition in place

Chapter 5 Mass, Momentum, and Energy Equations

Uploaded by

Copyright:

Available Formats

Chapter 5 Mass, Momentum, and Energy Equations

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Chapter 5 Mass, Momentum, and Energy Equations

Uploaded by

Copyright:

Available Formats

57:020 Mechanics of Fluids and Transport Processes

Professor Fred Stern Fall 2006