Shell Momentum Balances

Outline

1.Flow Through a Vertical Tube

2.Flow Through an Annulus
3.Exercises
Flow Through a Vertical Tube

The tube is oriented

vertically.

What will be the

velocity profile of a
fluid whose direction
of flow is in the +z-
direction
(downwards)?
Flow Through a Vertical Tube

Same system, but

this time gravity will
also cause
momentum flux.
 Flow Through a Vertical Tube

 rate of momentum   rate of momentum 

     force of gravity 
 in by molecular 
  out by molecular    acting on system   0
 transport   transport   
   

pressure : PA z  0  PA z  L
net momentum flux :  rz A1 r   rz A2 r  r
gravity : + gV (why positive?)

Adding all terms together:

P  2 r r  z  0  P  2 r r  z  L   rz  2 rL  r   rz  2 rL  r r
  g (2 r rL)  0
 Flow Through a Vertical Tube
P  2 r r  z 0  P  2 r r  z L   rz  2 rL  r   rz  2 rL  r r   g (2 r rL)  0

Dividing by 2 Lr :
 P z  0  P z  L   rz r r   rz r r  r
 r    gr  0
 L  r

Let r  0 :
 P0  PL  r  d  r   gr  0
   rz 
 L  dr
 Flow Through a Vertical Tube
 P0  PL  d
  r   rz r    gr  0
 L  dr
Rewriting:
d  P0  PL   P0  PL  g (0)   gL 
 rz r     g r    r
dr  L   L L 
d  P0   g (0) PL   gL 
 rz r     r
dr  L L 

d  0 L 
We let: z  Pz   gz  rz r    r
dr  L 
 Flow Through a Vertical Tube

d  P0  PL 
 rz r    r
Flow through a
dr  L  circular tube

d  0 L  Flow through a

 rz r    r
dr  L  vertical tube
 Flow Through a Vertical Tube

0 L  2 2
vz    R  r 
 4 L 

 0 L  2 Hagen-Poiseuille
vave  D
 32L  Equation
 Outline

1.Flow Through a Vertical Tube

2.Flow Through an Annulus
3.Exercises
Flow Through an Annulus

Liquid is flowing upward

through an annulus (space
between two concentric
cylinders)

Important quantities:
R : radius of outer cylinder
κR : radius of inner
cylinder
Flow Through an Annulus

Assumptions:
1. Steady-state flow
2. Incompressible fluid
3. Only Vz component is
significant
4. At the solid-liquid interface,
no-slip condition
5. Significant gravity effects
6. Vmax is attained at a
distance λR from the
center of the inner cylinder
(not necessarily the center)
 Flow Through an Annulus

 rate of momentum   rate of momentum 

     force of gravity 
 in by molecular 
  out by molecular    acting on system   0
 transport   transport   
   

pressure : PA z  0  PA z  L
net momentum flux :  rz A1 r   rz A2 r  r
gravity :   gV (why negative?)

Adding all terms together:

P  2 r r  z  0  P  2 r r  z  L   rz  2 rL  r   rz  2 rL  r r
  g (2 r rL)  0
 Flow Through an Annulus
 P0  PL  d
  r   rz r    gr  0
 L  dr
Rewriting:
d  P0  PL   P0  PL  g (0)   gL 
 rz r     g r    r
dr  L   L L 
d  P0   g (0) PL   gL 
 rz r     r
dr  L L 

d  0 L 
We let: z  Pz   gz  rz r    r
dr  L 
 Flow Through an Annulus
d  0 L 
 rz r    r
dr  L 

Solving: BOUNDARY CONDITION!

d  0 L  At a distance λR from the center of
 rz r    r the inner cylinder, Vmax is attained in
dr  L  the annulus, or zero momentum flux.
  L  2
 rz r   0  r  C1  0 L  C1
 2L  0  R 
 2L  R
 0 L  C1
 rz   r 
 2L  r   L 
 
2
 C1    0  R
 2L 
 Flow Through an Annulus

0 L   R 
2
0 L 
 rz   r   
 2L   2L  r

Rewriting:

 rz 
0 L  R  r  2  R 

2L       
 R   r 

From the definition of flux:

dv z 0 L  R  r  2  R 
dv
 rz   z
      
dr
dr 2 L  R   r 
 Flow Through an Annulus

dv z 0 L  R  r  2  R 
      
dr 2 L  R   r 

Solving:
0 L  R 1  r 
 2

vz        R ln  r   C2 
2

2 L R  2  
 Flow Through an Annulus

0 L  R  1  r 2  
vz       2R ln  r   C2 
2L R  2  
Rewriting:
  L R 2
 r
2
 R  Take out R/2
vz      2 ln  r    C2 
0 2
Multiply r in log term
4 L  R   R  by R/R (or 1)

0 L  R  r 
2
 2
2  r   
vz      2 ln   ln(R)   C2  Expand log term
4 L  R   R  
  L R 2
 r
2
r  Lump all constants
vz      2 ln   C2  into C2
0 2

4 L  R  R 
 Flow Through an Annulus

0 L  R  r 
2
 2
r 
vz      2 ln   C2 
2

4 L  R  R 

We have two unknown constants: C2 and λ

We can use two boundary conditions:

No-slip Conditions
At r = κR, vz = 0
At r = R, vz = 0
 Flow Through an Annulus

  L R 2
 r
2
r 
vz      2 ln   C2 
0 2

4 L  R  R 
Using B.C. #1:

0
0 L  R2 2
  2 2 ln    C2 
4 L
0   2  2 2 ln    C2
 C 2  1
Using B.C. #2:
0 L  R 2  2
1
0 1  C2   2 
2

0  1  C2
4 L
ln  
 Flow Through an Annulus
 C 2  1
  L R 2
 r
2
r 
 2
1
vz      2 ln   C2 
0 2

4 L  2 
2
 R  R  ln  


0 L  R 2
 r    2  1   r  
2

vz      ln   1 
4 L  R   ln   R  
 Shell Balances

1. Identify all the forces that influence the flow

(pressure, gravity, momentum flux) and their
directions. Set the positive directions of your axes.
2. Create a shell with a differential thickness across the
direction of the flux that will represent the flow
system.
3. Identify the areas (cross-sectional and surface areas)
and volumes for which the flow occurs.
4. Formulate the shell balance equation and the
corresponding differential equation for the
momentum flux.
 Shell Balances

5. Identify all boundary conditions (solid-liquid, liquid-

liquid, liquid-free surface, momentum flux values at
boundaries, symmetry for zero flux).
6. Integrate the DE for your momentum flux and
determine the values of the constants using the BCs.
7. Insert Newton’s law (momentum flux definition) to
get the differential equation for velocity.
8. Integrate the DE for velocity and determine values of
constants using the BCs.
9. Characterize the flow using this velocity profile.
 Shell Balances

Important Assumptions*
1. The flow is always assumed to be at steady-
state.
2. Neglect entrance and exit effects. The flow is
always assumed to be fully-developed.
3. The fluid is always assumed to be
incompressible.
4. Consider the flow to be unidirectional.

*unless otherwise stated

Design Equations for Laminar
and Turbulent Flow in Pipes
 Outline

1.Velocity Profiles in Pipes

2.Pressure Drop and Friction Loss (Laminar
Flow)
3.Friction Loss (Turbulent Flow)
4.Frictional Losses in Piping Systems
 Velocity Profiles in Pipes

Recall velocity profile in a circular tube:

 P0  PL  2 2
vz    R  r 
 4 L 
1. What is the shape of this profile?
2. The maximum occurs at which region?
3. What is the average velocity of the fluid
flowing through this pipe?
Velocity Profiles in Pipes
 Velocity Profiles in Pipes
Velocity Profile in a Pipe:

 P0  PL  2 2
vz    R  r 
 4 L 
Average Velocity of a Fluid in a Pipe:

 P0  PL  2
vave  D
 32L 
Maximum vs. Average Velocity
 Outline

1.Velocity Profiles in Pipes

2.Pressure Drop and Friction Loss (Laminar
Flow)
3.Friction Loss (Turbulent Flow)
4.Frictional Losses in Piping Systems
 Recall: Hagen-Poiseuille
Equation

 P0  PL  2
vave  D
 32L 

Describes the pressure drop and flow of

fluid (in the laminar regime) across a
conduit with length L and diameter D
 Hagen-Poiseuille Equation

32Lv ave
P0  PL  2
D

Pressure drop / Pressure loss (P0 – PL):

Pressure lost due to skin friction
 Friction Loss

32Lv ave
P0  PL  2
D
In terms of energy PO  PL 32Lvave
Ff  
lost per unit mass:  D 2

Mechanical energy lost due to friction in

pipe (because of what?)
 Friction Factor

Definition: Drag force per wetted surface

unit area (or shear stress at the surface)
divided by the product of density times
velocity head

 P0  PL  AC  AS
S
f 

 v 2
2
  v 22
 
 Friction Factor

2
Ff L v
 4 fF
gc D 2gc
Frictional force/loss head is proportional
to the velocity head of the flow and to
the ratio of the length to the diameter of
the flow stream
 Friction Factor for Laminar Flow

Consider the Hagen-Poiseuille equation

(describes laminar flow) and the
definition of the friction factor:
 P0  PL  2 Ff PO  PL L v2
vave  D   4 fF
 32L  gc g D 2gc

16
Prove: fF  Valid only for laminar flow
NRe
 Outline

1.Velocity Profiles in Pipes

2.Pressure Drop and Friction Loss (Laminar
Flow)
3.Friction Loss (Turbulent Flow)
4.Frictional Losses in Piping Systems
 Friction Factor for Turbulent
Flow
2
Ff L v
 4 fF
gc D 2gc

1. Friction factor is dependent on NRe and

the relative roughness of the pipe.
2. The value of fF is determined
empirically.
 Friction Factor for Turbulent
Flow
How to compute/find the value of the friction factor for
turbulent flow:

1. Use Moody diagrams.

- Friction factor vs. Reynolds number with a series of
parametric curves related to the relative roughness

2. Use correlations that involve the friction factor f.

- Blasius equation, Colebrook formula, Churchill
equation (Perry 8th Edition)
 Moody Diagrams
Important notes:
1. Both fF and NRe are plotted in logarithmic scales.
Some Moody diagrams show fD (Darcy friction
factor). Make the necessary conversions.
2. No curves are shown for the transition region.
3. Lowest possible friction factor for a given NRe in
turbulent flow is shown by the smooth pipe line.
 Friction Factor Correlations
1. Blasius equation for turbulent flow in smooth
tubes:
0.079
fF  4000  NRe  10 5
0.25
NRe
2. Colebrook formula

1   2.51 
 2log10   
fD  3.7D NRe fD 
 
 Friction Factor Correlations
3. Churchill equation (Colebrook formula explicit in fD)

1  0.27  7  
0.9

 2log10     
fD  D  NRe  
 
4. Swamee-Jain correlation
0.25
fD 
  5.74 
2log10   0.9 
 3.7D NRe 
 Equivalent Roughness, ε

Materials of Construction Equivalent Roughness (m)

Copper, brass, lead (tubing) 1.5 E-06
Commercial or welded steel 4.6 E-05
Wrought iron 4.6 E-05
Ductile iron – coated 1.2 E-04
Ductile iron – uncoated 2.4 E-04
Concrete 1.2 E-04
Riveted Steel 1.8 E-03
 Frictional Losses for Non-Circular
Conduits
Instead of deriving new correlations for f, an approximation
is developed for an equivalent diameter, Deq, which may be
used to calculate NRe and f.

S
Deq  4RH  4
Pw
where RH = hydraulic radius
S = cross-sectional area
Pw = wetted perimeter: sum of the length
of the boundaries of the cross-section
actually in contact with the fluid
 Equivalent Diameter (Deq)

S
Deq  4RH  4
Pw
Determine the equivalent diameter of the
following conduit types:
1. Annular space with outside diameter Do and
inside diameter Di
2. Rectangular duct with sides a and b
3. Open channels with liquid depth y and liquid
width b