جريان موائع 2

Chapter 6 Dynamic Similarity and Dimensional Analysis 129
CHAPTER 6
Dynamic Similarity and
Dimensional Analysis
6.1 Definition of Physical Similarity.
Two systems described by the same physics operating under different set
of conditions are said to be physically similar in respect of certain specified
physical quantities, when the ratio of corresponding magnitudes of these
quantities between the two systems is the same everywhere.
There are three types of similarities as in following chart which
constitute the complete similarity between problems of same kind.
Geometric Similarity
Complete similarity Dynamic Similarity
Kinematic Similarity
6.2 Geometric Similarity (G.S).

Geometric similarity implies the similarity of shape such that, the ratio
of any length in one system to the corresponding length in other system is the
same everywhere.
 Prototype:- is the full size or actual scale systems.
 Models:- is the laboratory scale systems.
 The model and prototype may be of identical size, although the two
may then differ in regard to other factors such as velocity and
properties of the fluid.
𝑙𝑒𝑛𝑔𝑡ℎ 𝑜𝑓 𝑚𝑜𝑑𝑒𝑙
𝑙𝑒𝑛𝑔𝑡ℎ 𝑟𝑎𝑡𝑖𝑜 = 𝑙𝑒𝑛𝑔𝑡ℎ 𝑜𝑓 𝑝𝑟𝑜𝑡𝑜𝑡𝑦𝑝𝑒
𝐿𝑚
Or 𝐿𝑟 = 𝐿𝑝
𝑚𝑜𝑑𝑒𝑙 𝑎𝑟𝑒𝑎 𝐿2𝑚
𝑎𝑟𝑒𝑎 𝑟𝑎𝑡𝑖𝑜 = 𝑝𝑟𝑜𝑡𝑜𝑡𝑦𝑝𝑒 𝑎𝑟𝑒𝑎 = 𝐿2𝑝
𝐿2 = 𝐿2𝑟
𝐿𝑟 Known as the model ratio or is the scale factor
6.3 Kinematic Similarity (K.S).
Kinematic similarity refers to similarity of motion
Distance Similarity of length (i.e., G.S.)
Motion
Time Similarity of time intervals.
𝐿𝑚
𝑉𝑚 ⁄𝑇 𝐿𝑚 𝑇𝑚 𝐿𝑟
𝑚
𝑣𝑒𝑙𝑜𝑐𝑖𝑡𝑦 𝑟𝑎𝑡𝑖𝑜 = = 𝐿𝑝 = ÷ =
𝑉𝑝 ⁄𝑇 𝐿𝑝 𝑇𝑝 𝑇𝑟
𝑝
𝐿𝑚
⁄ 2 𝐿𝑚 𝑇2 𝐿
𝑇𝑚
𝑎𝑐𝑐𝑒𝑙𝑎𝑟𝑡𝑖𝑜𝑛 𝑟𝑎𝑡𝑖𝑜 = 𝐿𝑝 = ÷ 𝑇𝑚2 = 𝑇𝑟2
⁄ 2 𝐿𝑝 𝑝 𝑟
𝑇𝑝
𝐿𝑚 3⁄
𝑄𝑚 𝑇𝑚 𝐿3𝑚 𝑇𝑚 𝐿3𝑟
𝑓𝑙𝑜𝑤 𝑟𝑎𝑡𝑒 𝑟𝑎𝑡𝑖𝑜 = = 𝐿𝑝 3
= ÷ =
𝑄𝑝 𝐿3𝑃 𝑇𝑝 𝑇𝑟
⁄𝑇
𝑝
Therefore, geometric similarity is a necessary condition for the kinematic

similarity to be achieved, but not the sufficient one.
6.4 Dynamic Similarity (D.S).

Dynamic similarity is the similarity of forces, in dynamically similar
system, the magnitudes of forces at correspondingly similar point in each
system are in a fixed ratio. In a system involving flow of fluid, different
forces due to different causes may act on a fluid element. These forces are as
follows;
 Viscous force (due to viscosity) 𝐹⃗𝑣
 Pressure force (due to different in pressure) 𝐹⃗𝑝
 Gravity force (due to gravitational attraction) 𝐹⃗𝑔
 Capillary force (due to surface tension ) 𝐹⃗𝑐
 Compressibility force (due to elasticity) 𝐹⃗𝑒
According to Newton's law, the resultant 𝐹𝑅 of all these forces, will

cause the accelartion of a fluid element, hence
𝐹⃗𝑅 = 𝐹⃗𝑣 + 𝐹⃗𝑝 + 𝐹⃗𝑔 + 𝐹⃗𝑐 + 𝐹⃗𝑒 (6.1)
The inertia force 𝐹⃗𝑖 is defined as equal and opposite to the resultant
accelerating force 𝐹⃗𝑅
𝐹⃗𝑖 = −𝐹⃗𝑅
∴ 𝐸𝑞. (6.1) 𝑐𝑎𝑛 𝑏𝑒 𝑒𝑥𝑝𝑟𝑒𝑠𝑒𝑑 𝑎𝑠
𝐹⃗𝑣 + 𝐹⃗𝑝 + 𝐹⃗𝑔 + 𝐹⃗𝑐 + 𝐹⃗𝑒 + 𝐹⃗𝑖 = 0
For dynamic similarity, the magnitude ratios of these forces have to be same
for both prototype and the model. The inertia force 𝐹⃗𝑖 is usually taken as the
common one to describe the ratio
|𝐹⃗𝑣 | |𝐹⃗𝑝 | |𝐹⃗𝑔 | |𝐹⃗ | |𝐹⃗ |
|𝐹⃗𝑖 |
, |𝐹⃗ | , |𝐹⃗ | , |𝐹⃗𝑐 | , |𝐹⃗𝑒|
𝑖 𝑖 𝑖 𝑖
a- Inertia Force.
The inertia force is the force acting on a fluid element is equal in
magnitude to the mass of the element multiplied by its acceleration.
Mass of element ∝ 𝜌𝐿3 ,  is the density, L is the characteristic length,
acceleration of a fluid element is the rate change of velocity in that direction
change with time
𝑉 𝐿
𝑎 ∝ ; 𝑡∝𝑉
𝑡
𝑉2
∴ 𝑎 ∝ 𝐿
𝑉2
The magnitude of inertia force is thus proportional to 𝜌 𝐿3 = 𝜌 𝐿2 𝑉 2
𝐿
This can be written as |𝐹⃗𝑖 | ∝ 𝜌𝐿2 𝑉 2 (6.2)
b- Viscous Force.
The viscous force arises from shear stress in a flow of fluid, therefor,
we can write magnitude of viscous force 𝐹⃗𝑣 = shear stress * surface area.
Shear stress = (viscosity) * rate of shear strain
𝑉
Where rate of shear strain ∝ velocity gradient ∝ 𝐿
2
Surface area ∝ 𝐿
⃗⃗v | ∝ μ V L2 ∝ μVL
|F (6.3)
L
c- Pressure Force.
The pressure force arises due to the difference of pressure in a flow field.
Hence it can be written as
⃗⃗P | ∝ ∆PL2
|F (6.4)
Where Δp is some characteristic pressure in the flow.
d- Gravity Force.
The gravity force on a fluid element is its weight, hence,
⃗⃗
|Fg | ∝ ρL3 g (6.5)
Where g is the acceleration due to gravity (or weight per unit mass)
e- Capillary or Surface Tension Force.

The capillary force arises due to the existence of an interface between
two fluids. It is equal to the coefficient of surface tension 𝜎 multiplied by the
length of a linear element on the surface perpendicular to which the force
acts, therefore,
|𝐹⃗𝑐 | ∝ 𝜎𝑙 (6.6)
f- Compressibility or Elastic Force.

For a given compression (a decrease in volume), the increase in
pressure is proportional to the bulk modulus of elasticity E (∆𝑝 ∝ 𝐸), this
gives rise to a force known as the elastic force .
⃗⃗c | ∝ EL2
|F (6.7)
Note, the flow of fluid in practice does not involve all the forces
simultaneously.
6.4.1 D. S. of Flow Governed by Viscous, Pressure and Inertia Forces.

The ratios of the representative magnitudes of these forces with the help
of Eq's (4.2) to (4.5) as follows:
𝑉𝑖𝑠𝑐𝑜𝑢𝑠𝑓𝑜𝑟𝑐𝑒 |𝐹⃗ | 𝜇𝑉𝐿 𝜇

= |𝐹⃗𝑣| ∝ 𝜌𝑉 2 𝐿2 = 𝜌𝑉𝐿 (6.8)
𝐼𝑛𝑒𝑟𝑡𝑖𝑎𝑓𝑜𝑟𝑐𝑒 𝑖
𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒𝑓𝑜𝑟𝑐𝑒 |𝐹⃗𝑝 | ∆𝑝𝐿2 ∆𝑝

= |𝐹⃗𝑖 |
∝ 𝜌𝑉 2 𝐿2 = 𝜌𝑉 2 (6.9)
𝐼𝑛𝑒𝑟𝑡𝑖𝑎𝑓𝑜𝑟𝑐𝑒
The term LV/ is known as Reynolds number, Re.

𝐼𝑛𝑒𝑟𝑡𝑖𝑎 𝑓𝑜𝑟𝑐𝑒
𝑅𝑒 ∝ 𝑉𝑖𝑠𝑐𝑜𝑢𝑠 𝑓𝑜𝑟𝑐𝑒 is thus proportion to the magnitude ratio.
∆𝑝
The term 𝜌𝑉 2 is known as Euler number, Eu.
∴ 𝑹𝒆 & 𝑬𝒖 Represent the criteria of D.S. for the flows which are affected
only by viscous, pressure and iertia forces. For example are
1- The full flow of fluid in a completely closed conduit

2- Flow of air past a low – speed aircraft
3- The flow of water past a submarine deeply submerged to produce no
waves on the surface.
Hence, Re & Eu for a complete dynamic similarity between prototype and

model must be the same for two. Thus
𝜌𝑝 𝐿𝑝 𝑉𝑝 𝜌𝑚 𝐿𝑚 𝑉𝑚
= (6.10)
𝜇𝑝 𝜇𝑚
∆𝑝𝑝 ∆𝑝𝑚
=𝜌 (6.11)
𝜌𝑝 𝑉𝑝2 2
𝑚 𝑉𝑚
6.4.2 D.S. of Flow Governed by Gravity and Inertia Forces.

A flow of the type in which significant force are gravity and inertia
forces, is found when a free surface is present. For example are
1- The flow of a liquid in an open channel.
2- The wave motion caused by the passage of a ship through water.
3- The flow over weirs and spillways.
The condition for D.S. of such flows requires
 The equality of Eu.
 The equality of the magnitude ratio of gravity to inertia force at
corresponding points in the system.
𝐺𝑟𝑎𝑣𝑖𝑡𝑦𝑓𝑜𝑟𝑐𝑒 |𝐹⃗𝑔 | 𝜌𝐿3 𝑔 𝐿𝑔
= |𝐹⃗𝑖 |
∝ 𝜌𝑉 2 𝐿2 = 𝑉 2 (6.12)
1
(𝐿𝑔)2
The reciprocal the term is known as Froude number, Fr
𝑉
𝑉
∴ 𝐹𝑟 = 1
(𝐿𝑔)2
∴ Dynamic similarity between prototype & model is the equality of Froude
number
√𝐿𝑝 𝑔𝑝 √𝐿𝑚 𝑔𝑚
= (6.13)
𝑉𝑝 𝑉𝑚
6.4.3 D.S. of Flows with Surface Tension as the Dominant Force.

Surface tension forces are important in certain classes of practical
problems such as:
1- Flows in which capillary waves appear.
2- Flows of small jets and thin sheets of liquid injected by nozzle in air.
3- Flow of a thin sheet of liquid over a solid surface.
Dynamic similarity is the magnitude ratio
|𝐹⃗𝑐 | 𝜎𝐿 𝜎
|𝐹⃗𝑖 |
∝ 𝜌𝑉 2 𝐿2 = 𝜌𝑉 2 𝐿 (6.14)
𝜎
The term 𝜌𝑉 2 𝐿 is usually knows as Weber number, 𝑾𝒃 .
For dynamic similarity (𝑊𝑏 )𝑚 = (𝑊𝑏 )𝑝
𝜎𝑚 𝜎𝑝
𝑖. 𝑒. , 𝜌 2 =𝜌 2 (6.15)
𝑚 𝑉𝑚 𝐿𝑚 𝑝 𝑉𝑝 𝐿𝑝
6.4.4 D.S. of Flows with Elastic Force.

The magnitude ratio of inertia to elastic force becomes
Inertia force |𝐹⃗ | 𝜌𝑉 2 𝐿2 𝜌𝑉 2
= |𝐹⃗𝑖 | ∝ = (6.16)
Elastic force 𝑒 𝐸𝐿2 𝐸
𝜌𝑉 2
The parameter is known as Cauchy number.
𝐸
For dynamic similarity flow (Cauchy)m = (Cauchy)p
𝜌 𝑉2 𝜌𝑝 𝑉𝑝2
𝑖. 𝑒. , (𝐸𝑚)𝑚
𝑚
= (𝐸 ) (6.17)
𝑠 𝑠 𝑃
If the flow is isentropic 𝐸 = 𝐸𝑠 is isentropic bulk modulus of elasticity.
𝐸
i = sound wave propagates through a fluid medium =√ 𝜌𝑠
∴ 𝑡ℎ𝑒 𝑡𝑒𝑟𝑚 𝜌𝑉 2 /𝐸𝑠 can be written as 𝑉 2 /𝑖 2
𝑉
The ratio ( 𝑖 ) = 𝑴𝒂 is known as Mach number, in the flow of air past
high-speed aircraft, missiles, propellers and rotary compressors. In these
cases equality of Mach number is a condition of dynamic similarity.
Therefore
(𝑀𝑎 )𝑝 = (𝑀𝑎 )𝑚
𝑉𝑝 𝑉
𝑖. 𝑒 ( 𝑖 ) = ( 𝑖 𝑚 ) (6.18)
𝑝 𝑚
Dimensionless Representation Name Recommended

terms magnitude ratio symbol
of the force
LV/ Inertiaforce Reynolds Re

viscous force number
∆p/V 2 pressureforce Euler number Eu

Inertiaforce
V/(Lg)1/2 Inertia force Froud number Fr

Gravity force
σ surface tension Weber number Wb
V 2 L Inertia force
V/√Es / Inertiaforce Mach number M𝑎

Elastic force
Ex.1
When tested in water at 20C flowing at 2 m/s, an 8-cm diameter sphere
has a measured drag of 5 N. What will be the velocity and drag force on a
1.5m diameter weather balloon moored in sea-level standard air under
dynamically similar condition?
Sol.
𝑘𝑔 𝑘𝑔
For water at 20C 𝜌 ≈ 998 𝑚3 , & 𝜇 = 0.001 𝑚.𝑠
𝑘𝑔 𝑘𝑔
For air at sea level 𝜌 ≈ 1.2255 𝑚3 , 𝜇 = 1.78 ∗ 10−5 𝑚.𝑠
The balloon velocity follows from dynamic similarity, which requires
identical (Reynolds number).
ρVD
𝑹𝒆𝒎 = 𝑹𝒆𝒑 = μ
998∗(2.0)(0.08)
Rem = ≈ 1.6 ∗ 105 = Rep
0.001
𝑉𝑏𝑎𝑙𝑙𝑜𝑛 (1.5) 𝑚
1.6 ∗ 105 = 1.2255 − − − −→ 𝑉𝑏𝑎𝑙𝑙𝑜𝑛 ≈ 1.55
1.78∗10−5 𝑠
Then the two spheres will be identical drag coefficients:
∆𝑝 𝐹 𝐹 F
CD,m = 𝜌𝑉 2 = [ 𝜋𝑑2
] =[ 𝜋𝑑2
] = ρV2 d2
𝜌𝑉 2 ( ) 𝜌𝑉 2 ( )
4 4
𝑚 𝑝
5 𝑏𝑎𝑙𝑜𝑛 𝐹
𝐶𝐷,𝑚 = 𝜋 = 0.4986 = 𝐶𝐷,𝑝 = 1.2255(1.55)2 (1.5)2 𝜋/4
998(2)2 (0.08)2
4
Solve for 𝑓𝑏𝑎𝑙𝑙𝑜𝑛 ≈ 1.296 𝑁
Ex.2 A model of a reservoir having a free water surface within it is drained
in 3 minutes by opening a sluice gate. The geometrical scale of the model is
(1/100). How long would it take to empty the prototype?
Sol.
𝐿3
𝑚
𝑄𝑚 𝑇𝑚 𝐿3𝑟 𝑇𝑚
= 𝐿3
= 𝑇𝑟 =
𝑄𝑝 𝑝 𝑇𝑟 𝑇𝑝
𝑇𝑝
The forces control the flow is

(𝐹𝑔 ) (𝜌𝑔𝐿3 ) 𝑊𝑚 𝐿3𝑚
1- Gravity force = 𝑚
= (𝜌𝑔𝐿3 )𝑚 = = 𝑊𝑟 𝐿3𝑟
(𝐹𝑔 ) 𝑊𝑝 𝐿3𝑃
𝑝 𝑝
(𝐹𝑖 )𝑚 (𝑚 𝑎)𝑚 𝜌 𝐿3𝑚 𝐿 𝐿
2- Inertia force = (𝐹𝑖 )𝑝
= (𝑚 𝑎)𝑝
= 𝜌𝑚 3 (𝑇𝑟2 ) = 𝜌𝑟 𝐿3𝑟 𝑇𝑟2
𝑝 𝐿𝑝 𝑟 𝑟
By equating the two ratio

𝐿𝑟
𝑊𝑟 𝐿3𝑟 = 𝜌𝑟 𝐿3𝑟
𝑇𝑟2
𝐿
𝜌𝑟 𝑔𝑟 𝐿3𝑟 = 𝜌𝑟 𝐿3𝑟 𝑇𝑟2
𝑟
1
𝐿𝑟 2
𝑇𝑚
𝑇𝑟2 = 𝑔 − − − −→ = 100
; 𝑠𝑖𝑛𝑐𝑒 𝑔𝑟 = 1 𝑇𝑚 = 3 𝑚𝑖𝑛
𝑟 𝑇𝑝2 1
∴ 𝑇𝑝 = √900 = 30 min
6.5 The Application of D.S and the Dimensional Analysis.

6.5.1 The concept.
A physical problem may be characterized by a group of dimensionless
similarity parameters or variables rather than by the original dimensional
variables. This gives a clue to reduction in the number of parameters
requiring separate consideration in an experimental investigation.
𝜌𝑉𝐷ℎ
Ex:- 𝑅𝑒 = 𝑅𝑒 2000  4000 by varying V without change in any
𝜇
other independent dimensional variable .
In fact, the variation in the 𝑅𝑒 physically implies the variation in any of the
dimensional parameters defining it.
6.5.2 Dimensional Analysis.

The dimensional analysis is a mathematical technique by which can be
determining many dimensionless parameters and solving several engineering
problems.
There are two existing approaches:-
1- Indicial method.
2- Buckingham's pi theorem.
The dimensional analysis can be explain by the following,
 The Various physical qumtities used in fluid phenomenon can be
expressed in terms of fundamental quantities or primary quantities.
 Fundamental quantities are Mass (M), Length (L), Time (T),
Temperature (𝜃) is used for compressible flow.
 The quantities which are expressed in terms of the fundamental or
primary quantities are called derived or secondary quantities as
(velocity, area, acceleration)
 The expression for a derived quantities in terms of the primary
quantities is called the dimension of the physically quantities.
 A quantity may either be expressed dimensionally in M-L-T or F-L-T
system.
Ex.3
Determine the dimensions of the following quantities.
(i) Discharge.
(ii) Kinematic viscosity.
(iii) Force.
(iv) Specific weight.
Sol.
(i) Discharge = area * velocity
𝐿 𝐿3
= 𝐿2 ∗ = = 𝐿3 𝑇 −1
𝑇 𝑇
(ii) Kinematic Viscosity()=/
𝑑𝑢
Where () given by (τ) =  𝑑𝑦
𝑓𝑜𝑟𝑐𝑒
𝜏 𝑠ℎ𝑒𝑎𝑟𝑠𝑡𝑟𝑒𝑠𝑠
 = 𝑑𝑢 = 𝐿 1 = 𝐴𝑟𝑒𝑎
1
⁄𝑑𝑦 ×
𝑇 𝐿 𝑇
𝐿
𝑚𝑎𝑠𝑠×𝑎𝑐𝑐𝑒𝑙𝑒𝑟𝑎𝑡𝑖𝑜𝑛 𝑀× 2 𝑀×𝐿
𝑇
𝜇= 1 = 1 = 1
𝐴𝑟𝑒𝑎× 𝐿2 × 𝐿2 𝑇 2 ×
𝑇 𝑇 𝑇
𝑀
𝜇 = 𝐿𝑇 = 𝑀𝐿−1 𝑇 −1
𝑚𝑎𝑠𝑠 𝑀
𝑎𝑛𝑑 𝜌 = 𝑣𝑜𝑙𝑢𝑚𝑒 = 𝐿3 = 𝑀𝐿−3
𝜇 𝑀𝐿−1 𝑇 −1
∴=𝜌= = 𝐿2 𝑇 −1
𝑀𝐿 −3
(iii) Force = mass * acceleration
𝑙𝑒𝑛𝑔𝑡ℎ 𝑀𝐿
=𝑀∗ = = 𝑀𝐿𝑇 −2
𝑇𝑖𝑚𝑒 2 𝑇2
𝑀𝐿𝑇 −2
(iv) Specific weight= Weight/volume= force/volume= =
𝐿3
𝑀𝐿−2 𝑇 −2
6.5.3 Dimensions of Physical Quantities.
All physical quantities are expressed by magnitude and units as an
example;
Velocity = 8 m/s; Acceleration = 10 m/s2
( 8 & 10 are the magnitudes but m/s & m/s2 are the dimensions)
SI(system international ) units, in fluid mechanics the primary physical
quantities or (base dimensions) are the [ Mass, Length, Time,
Temperature] are symbolized as [ 𝑴, 𝑳, 𝑻, 𝜽]. Any physical quantity can be
expressed in terms of these primary quantities by using the basic
mathematical definition of the quantity, resulting is known as the dimension
of the quantity, by substitute the mass by force (F)
𝐹 = 𝑀𝐿𝑇 −2 − − − − − −→ 𝑀 = 𝐹𝑇 2 𝐿−1
Let us take some examples.
1- Dimension of stress.
𝑓𝑜𝑟𝑐𝑒
𝑠ℎ𝑒𝑎𝑟 𝑠𝑡𝑟𝑒𝑠𝑠 𝜏 = , 𝑓𝑜𝑟𝑐𝑒 = 𝑚𝑎𝑠𝑠 ∗ 𝑎𝑐𝑐𝑒𝑙𝑎𝑟𝑖𝑜𝑛
𝑎𝑟𝑒𝑎
Dimension of acceleration = dimension of velocity/ dimension of time
𝐿
= 𝐷𝑖𝑚. 𝑜𝑓 𝑑𝑖𝑠𝑡𝑎𝑛𝑐𝑒/(𝐷𝑖𝑚. 𝑜𝑓 𝑡𝑖𝑚𝑒)2 = 𝑇2
𝐷𝑖𝑚. 𝑜𝑓 𝑎𝑟𝑒𝑎 (𝑙𝑒𝑛𝑔𝑡ℎ)2 2
=𝐿
𝑀𝐿
𝑇2
𝜏= = 𝑀𝐿−1 𝑇 −2
𝐿2
2- Dimension of viscosity.
𝑑𝑢
𝜏 = 𝜇 𝑑𝑦
𝜏
𝑂𝑟 , 𝜇= 𝑑𝑢
𝑑𝑦
The dimension of velocity gradient du/dy can be written as

𝐿
du/dy = dimension of u/ dimension of y= 𝑇

= 𝑇 −1
𝐿
Dimension of 
Dim.of τ 𝑀𝐿−1 𝑇 −2
∴  = Dim.of du/dy = = 𝑀𝐿−1 𝑇 −1
𝑇 −1
H.W. Drive the dimensions of the following physical Quantities in [M,L,T],
momentum, work , weight, flow rate , Power.
6.6 Rayleigh's Indicial Method (Method-1).

Based on the fundamental principle of dimensional homogeneity of
physical variables.
Procedure.
1- The dependent variable is identified and expressed as a product of all
the independent variables raised to an unknown integer exponent.
2- Equating the indices of (n) fundamental dimensions of the variables
involved, (n) independent equations are obtained.
3- These (n) equations are solved to obtain the dimensionless groups.
Ex.4
Let us illustrate this method by solving the pipe flow problem with
∆𝑝/𝑙 along the pipe.
Step_1. Here, the dependent variable ∆𝑝/𝑙 can be written as
∆p
= K(V a Dbh ρc μd ) Where, K is constant.
l
Step_2. Inserting the dimension of each variable in the above equation, we
obtain.
ML−2 T −2 = K(LT −1 )a (L)b (ML−3 )c (ML−1 T −1 )d
Equating the indices of M,L and T on both sides , we get
M] c + d = 1
𝐿] 𝑎 + 𝑏 − 3𝑐 − 𝑑 = −2
𝑇] − 𝑎 − 𝑑 = −2
Step-3-:- there are three equations and four unknowns. Solving these
equations in terms of the unknown d, we have
𝑎 =2−𝑑
𝑏 = −𝑑 − 1
𝑐 =1−𝑑
Hence, we can be written
∆𝑝
= 𝐾(𝑉 2−𝑑 𝐷ℎ −𝑑−1 𝜌1−𝑑 𝜇 𝑑 )
𝑙
∆p V2 ρ μ d
= K( D ) (VD ρ)
𝑙 h h
∆pDh 𝜇 𝑑
= 𝐾 (𝑉𝐷 𝜌)
𝑙ρV2 ℎ
Ex.5
Write the equation of displacement for a free fouling body in time T.
Assuming that the displacement dependent on weight, acceleration gravity
and time.
Sol.
𝑆 = 𝐹(𝑊, 𝑔, 𝑇) 𝑑𝑖𝑠𝑝𝑙𝑎𝑐𝑚𝑒𝑛𝑡
𝑆 = 𝐾 𝑊 𝑎 𝑔𝑏 𝑇 𝑐
The equation must be homogenate in dimension
M 0 L1 T 0 = K((MLT −2 )a (LT −2 )b (T)c )
Equating the indices of same dimension of quantities
𝑎=0
𝑎+𝑏 =1 b=1
−2𝑎 − 2𝑏 + 𝑐 = 0 ∴𝑐=2
0 2 )𝑜𝑟
∴ 𝑆 = 𝐾(𝑊 𝑔𝑇 𝑠 = 𝐾𝑔𝑇 2
Ex.6
Find the relation of Reynolds number by dimensional Analysis if
Re=F(,,V,L)
Sol.
𝑅𝑒 = 𝐹(𝜌, , 𝑉, 𝐿)
𝑅𝑒 = 𝐾 𝑎 𝑏 𝑉 𝑐 𝐿𝑑
𝑀0 𝐿0 𝑇 0 = 𝐾 (𝑀𝐿−3 )𝑎 (𝑀𝐿−1 𝑇 −1 )𝑏 (𝐿𝑇 −1 )𝑐 (𝐿)𝑑
𝑎+𝑏 =0 𝑎 = −𝑏
−3𝑎 − 𝑏 + 𝑐 + 𝑑 = 0 𝑑 = −2𝑏 − 𝑐 = −2𝑏 + 𝑏 = −𝑏
−𝑏 − 𝑐 = 0 𝑏 = −𝑐
𝑐 = −𝑏
∴ 𝑅𝑒 = 𝐾 𝜌−𝑏 𝜇 𝑏 𝑉 −𝑏 𝐿−𝑏
𝜇 𝑏
𝑅𝑒 = 𝐾 (𝜌𝑉𝐿) 𝐾 = 1, 𝑏 = −1
Ex.7
Find the dynamic pressure over a submerged body due to the flow of
uncompressible fluid. Assuming the pressure is function of density and
velocity.
Sol.
𝑝 = 𝐹(𝜌, 𝑉)
𝑝 = 𝐾 𝜌𝑎 𝑉 𝑏
𝐹1 𝐿−2 𝑇 0 = (𝐹 𝑎 𝑇 2𝑎 𝐿−4𝑎 )(𝐿𝑏 𝑇 −𝑏 )
From above
1=a , -2=-4a+b, 0=2a-b
∴ 𝑎 = 1 ,𝑏 = 2
𝑝 = 𝐾𝜌𝑉 2
Ex.8
Find the expression for the input power to a fan. By dimension
analysis, assuming the input power depends on the air density, velocity,
viscosity, fan diameter, rotation speed and sound velocity.
Sol.
power = K(ρa db V c ωd μe if )
By using (mass, length, time) as fundamental units.
ML2 T −3 = (ML−3 )a (L)b (L T −1 )c (T −1 )d (ML−1 T −1 )e (L T −1 )f
1=a+e a=1-e
2 = −3a + b + c − e + f then b = 5-2e-c-f
−3 = −c − d − e − f d= 3-c-e-f
Subsititute in power Eqn.
power = K ρ1−e d5−2e−c−f V c ω3−c−e−f μe if
−𝑒
𝜌𝑑2 𝜔 𝑑𝜔 −𝑐 𝑑𝜔 𝑓
𝑝𝑜𝑤𝑒𝑟 = 𝐾 [( ) (𝑉) ( ) ] 𝜔3 𝑑5 𝜌
𝜇 𝑖
The terms between brackets dimension less
1𝑠𝑡 𝑡𝑒𝑟𝑚 = 𝑅𝑒
V= R
2𝑛𝑑 𝑡𝑒𝑟𝑚 = 𝑓𝑎𝑛 𝑟𝑎𝑡𝑖𝑜.
3𝑛𝑑 𝑡𝑒𝑟𝑚 = 𝑀𝑎𝑐ℎ 𝑛𝑢𝑚𝑏𝑒𝑟.
6.7 Buckingham's Pi Theorem. (Method-2)

Assume, a physical phenomenon is described by
n = number of independent variables like x1, x2, x3....xn the phenomenon may
be expressed as
𝐹(𝑥1 , 𝑥2 , 𝑥3 , … … 𝑥𝑚 ) = 0 (6.19)
m = number of fundamental dimensions like mass, time, length and
temperature or force, length, time and temperature.
Buckingham's theorem defining as the phenomenon can be described in terms
of (n-m) independent dimensionless group like 𝜋1 , 𝜋2 , … … … 𝜋𝑚−𝑛 where 𝝅
terms, represent the dimensionless parameters and consist of different
combinations of a number of dimensional variables out of the n independent
variables.
Therefore Eq.(6.19) can be reduced to
𝐹(𝜋1 , 𝜋2 , … … . , 𝜋𝑚−𝑛 ) = 0 (6.20)
6.7.1 Mathematical Description of (𝝅 ) Pi Theorem.

A physical problem described by n number of variable involving m
number of fundamental dimensions (m < n) leads to a system of m linear
algebraic equations with n variables of the from
a11 x1 + a12 x2 + ⋯ a1n xn = b1
a21 x1 + a22 x2 + ⋯ a2n xn = b2 (6.21)
am1 x1 + am2 x2 + ⋯ amn xn = bm
Therefore all feasible phenomena are define with n > m
No physical phenomena is represent
n< m no solution
m=n one solution
All the parameter have fixed value.
In a matrix from
AX=b (6.22)
𝑎11 𝑎12 𝑎1𝑛
Where 𝐴 = [ 𝑎21 𝑎22 𝑎2𝑛 ]
𝑎𝑚1 𝑎𝑚2 𝑎𝑚𝑛
𝑥1 𝑏1
𝑥2 𝑏2
𝑋= . 𝑎𝑛𝑑 𝑏 = .
. .
[𝑥𝑛 ] [𝑏𝑚 ]
6.7.2 Procedure for Determination 𝝅 Terms.

m = Number of fundamental dimensions like mass, (M), time (T), Length
(L), temperature(𝜃)
n = number of independent variables or quantities included in physical
problem such as ( A1,A2,A3-----An) where A1,A2,A3-----An as pressure ,
viscosity and velocity, can also be expressed as
𝐹1(𝐴1 , 𝐴2 , 𝐴3 , − − − − 𝐴𝑛 ) = 0 (6.23)
(n-m) = number of dimensionless parameter (𝜋 ) like 𝜋1 , 𝜋2 , 𝜋3 , … 𝜋𝑛−𝑚
𝜋 , is represent the dimensionless parameters and consist of different
combinitions of a number of dimensional variable. Mathematically, if any
variable A1, depends on independent variable A2.A3…An the function A1 =
𝐹(𝐴2 , 𝐴3 , … 𝐴𝑛 )
According to 𝜋-theorem, Eq. (6.23) can be written in terms of 𝜋- terms
(dimensionless groups). Therefore the above equation can reduced to
𝐹1 (𝜋1 , 𝜋2 , … 𝜋𝑛−𝑚 ) = 0 (6.24)
The method of determining 𝜋 parameters is
- Select (m) of the (A) quantities with different dimensions
- The above selection which contains among them (m) dimensions
- Using the (m) selection as repeating variables together with one of the
other A quantities for each (𝜋). Each 𝜋-term contains (m+1) variables.
Note-1, It is essential that no one of the m selected quantities used as
repeating variable be derived from the other repeating variables.
Note-2, Let A1,A2,A3 contain M,L and T, not necessarily in each one , but
collectively.
Then the first 𝜋 parameter is made up as
𝑦1
𝜋1 = 𝐴2𝑥1 𝐴3 𝐴4𝑧1 𝐴1
𝑦2
The second 𝜋2 = 𝐴2𝑥2 𝐴3 𝐴4𝑧2 𝐴5 (6.25)
𝑥𝑛−𝑚 𝑦𝑛−𝑚 𝑧𝑛−𝑚
And so on until 𝜋𝑛−𝑚 = 𝐴2 𝐴3 𝐴4 𝐴𝑛
In a bove eqn's the exponents are to be determined
- The dimensions of (A) quantities are substituted
- The exponents of M,L and T are set equal to zero in 𝜋 parameters
There produce three equations in three unknowns for each 𝜋 parameter, so
that the x,y,z exponents can be determined , and Hence the 𝜋 parameter.
6.7.3 Selection of Repeating Variables (R.V).

1- (m) repeating variable must contain jointly all the fundamental
dimension
2- The repeating variable must not form the non-dimensional parameter
among them.
3- As far as possible, the dependent variable should not be selected as
repeating variable
4- No two repeating variables should have the same dimensions
5- The repeating variables should be chosen in such a way that one
variable contains geometric property ( e.g , length , L, diameter, d,
height, h) , other variable contains flow property (e.g velocity V,
acceleration a ) and the third variable contains fluid property ( e.g
mass density , weight density W , dynamic viscosity )
(i) L,V,
(ii) d,V,
(iii) l, V,m
(iv) d,V,
Ex.9
Show that the lift force 𝐹𝑙 on airfoil can be express as 𝐹𝐿 =
𝑉𝑑
𝑉 2 𝑑 2 ∅ ( 
, ∝)
Where = mass density , V = velocity of flow
= dynamic viscosity ∝= Angle of incidence
d = A characteristic depth
Sol.
Left force 𝐹𝐿 is function of;, V, d, m, ∝ mathematically, 𝐹𝐿 =
𝑓( , 𝑉, 𝑑, , ∝) − − − −(𝑖)
Or 𝐹1 ( 𝐹𝐿 , 𝜌 𝑉, 𝑑 𝜇, ∝) = 0 − − − − − − − (𝑖𝑖)
∴ Total number of variable, we have n=6
Writing dimensions of each variable
𝐹𝐿 = 𝑀 𝐿 𝑇 −2 , 𝜌 = 𝑀 𝐿−3 , 𝑉 = 𝐿𝑇 −1 , 𝑑 = 𝐿, 𝜇 = 𝑀𝐿−1 𝑇 −1 , ∝= 𝑀0 𝐿0 𝑇 0
Thus, number of fundamental dimensions, m=3
∴ 𝑁𝑢𝑚𝑏𝑒𝑟 𝑜𝑓 𝜋 − 𝑡𝑒𝑟𝑚𝑠 = 𝑛 − 𝑚 = 6 − 3 = 𝟑
Eq. (ii) can be written as : 𝐹1 (𝜋1 , 𝜋2 , 𝜋3 ) = 0 − − − − − (𝑖𝑖𝑖)
Each 𝜋-term contains (m+1) variables, where m=3 and also equal to
repeating variables (R.V). Choosing (d, V, ) as R.V
𝜋1 = 𝑑𝑎1 . 𝑉 𝑏1 . 𝜌𝑐1 . 𝐹𝐿
𝜋2 = 𝑑𝑎2 . 𝑉 𝑏2 . 𝜌𝑐2 . 𝜇
𝜋3 = 𝑑𝑎3 . 𝑉 𝑏3 . 𝜌𝑐3 . ∝
𝜋 − 𝑡𝑒𝑟𝑚:
𝜋1 = 𝑑𝑎1 . 𝑉 𝑏1 . 𝜌𝑐1 . 𝐹𝐿
𝑀0 𝐿0 𝑇 0 = 𝐿𝑎1 (𝐿 𝑇 −1 )𝑏1 (𝑀𝐿−3 )𝑐1 (𝑀𝐿𝑇 −2 )
Equating the exponents of M,L,T respectively, we get
𝑓𝑜𝑟 𝑀: 0 = 𝑐1 + 1 − − − −→ 𝑐1 = −1
𝑓𝑜𝑟 𝐿: 0 = 𝑎1 + 𝑏1 − 3𝑐1 + 1
𝑓𝑜𝑟 𝑇: 0 = −𝑏1 − 2 − −−→ 𝑏1 = −2
∴ 𝑎1 = −𝑏1 + 3𝑐1 − 1 = 2 − 3 − 1 = −2
Substituting the values of 𝑎1 , 𝑏1 𝑎𝑛𝑑 𝑐1 𝑖𝑛 𝜋1 , we get
𝐹
𝜋1 = 𝑑 −2 . 𝑉 −2 . 𝜌−1 . 𝐹𝐿 = 𝜌𝑉 2𝐿𝑑2
𝜋2 − 𝑡𝑒𝑟𝑚:
𝜋2 = 𝑑𝑎2 . 𝑉𝑏2 . 𝜌𝑐2 . 𝜇
𝑀0 𝐿0 𝑇 0 = 𝐿𝑎2 . (𝐿𝑇 −1 )𝑏2 . (𝑀𝐿−3 )𝑐2 . (𝑀𝐿−1 𝑇 −1 )
Equating the exponents of M,L and T respectively , we get
𝑓𝑜𝑟 𝑀: 0 = 𝑐2 + 1 − −−→ 𝑐2 = −1
𝐹𝑂𝑅 𝐿: 0 = 𝑎2 + 𝑏2 − 3𝑐2 − 1
For 𝑇: 0 = −𝑏2 − 1 − −−→ 𝑏2 = −1
∴ 𝑎2 = −𝑏2 + 3𝑐2 + 1 = 1 − 3 + 1 = −1
Substituting the values of 𝑎2 , 𝑏2 , 𝑎𝑛𝑑 𝑐2 𝑖𝑛 𝜋2 , 𝑤𝑒 𝑔𝑒𝑡
𝜇
𝜋2 = 𝑑−1 . 𝑉 −1 . 𝜌−1 . 𝜇 = 𝜌𝑉𝑑
𝜌𝑉𝑑
𝑜𝑟 𝜋2 = 𝜇
𝜋3 − 𝑡𝑒𝑟𝑚
𝜋3 = 𝑑𝑎3 . 𝑉𝑏3 . 𝜌𝑐3 . ∝
𝑀0 𝐿0 𝑇 0 = 𝐿𝑎3 . (𝐿𝑇 −1 )𝑏3 . (𝑀𝐿−3 )𝑐3 . (𝑀0 𝐿0 𝑇 0 )
Equating the exponents of M, L and T respectively , we get
𝑓𝑜𝑟 𝑀: 0 = 𝑐3 + 0 − −−→ 𝑐3 = 0
𝑓𝑜𝑟 𝐿: 0 = 𝑎3 + 𝑏3 − 3𝑐3 + 0
𝑓𝑜𝑟 𝑇: 0 = −𝑏3 + 0 − −−→ 𝑏3 = 0
∴ 𝑎3 = 0
∴ 𝜋3 = 𝑑0 . 𝑣 0 𝜌0 . ∝=∝
Substituting the values of 𝜋1 , 𝜋2 𝑎𝑛𝑑 𝜋3 𝑖𝑛 𝐸𝑞. (𝑖𝑖𝑖), 𝑤𝑒 𝑔𝑒𝑡
𝐹 𝜌𝑉𝑑
𝑓1 (𝜌𝑉 2𝐿𝑑2 , , ∝) = 0
𝜇
𝐹𝐿 𝜌𝑉𝑑
= ∅( ,∝)
𝜌𝑉 2 𝑑2 𝜇
𝜌𝑉𝑑
𝑜𝑟 𝐹𝐿 = 𝜌𝑉 2 𝑑 2 ∅ ( , ∝)
𝜇
Ex.10
The discharge through a horizontal capillary tube is thought to depend
upon the pressure drop per unit length, the diameter and the viscosity. Find
the form of the discharge equation.
Sol.
Quantity Dimensions
Discharge Q L3 T −1
Pressure drop/length M L−2 T −2

Δp/l
Diameter D L
Viscosity  𝑀 𝐿−1 𝑇 −1
∆𝑝
Then 𝐹 (𝑄, , 𝐷 , 𝜇) = 0
𝐿
Three dimension are used, and with four quantities these will be one 𝜋
parameter
𝒏 = 4 , 𝒎 = 3 − − − −→ 𝜋 = 𝑛 − 𝑚 = 4 − 3 = 1
∆𝑝
𝜋 = 𝑄 𝑋1 ( 𝑙 )𝑌1 𝐷 𝑍1 𝜇
Substituting in the dimension gives
𝜋 = (𝐿3 𝑇 −1 ) 𝑋1 (𝑀𝐿−2 𝑇 −2 )𝑌1 (𝐿𝑍1 )(𝑀 𝐿−1 𝑇 −1 ) = 𝑀0 𝐿0 𝑇 0
The exponents of each dimension must be the same on both sides of the
equation.
With L; 3𝑋1 − 2𝑌1 + 𝑍1 − 1 = 0 − − − − − −(𝑖)
M; 𝑌1 + 1 = 0 − − − −→ 𝑌1 = −1
T; −𝑋1 − 2𝑌1 − 1 = 0 − − − − − −→ 𝑋1 = 1
From Eq. (i) 𝑍1 = −4
𝜇
𝜋=𝑄 ∆𝑝
( )𝐷4
𝑙
After solving for Q
∆𝑝 𝐷 4
𝑄=𝐶 𝐿 𝜇
Ex.11
Consider pressure drop in a tube of length l, hydraulic diameter d,
surface roughness ∈, with fluid of density  and viscosity  moving with
average velocity U. Using Buckingham's 𝜋 - theorem obtain an expression
for Δp.
Sol.
This can be expressed as
𝑓(∆𝑝, 𝑈, 𝑑, 𝑙, ∈, 𝜌, 𝜇) = 0
Now n=7 since the phenomenon involves 7 independent parameters.
We select 𝜌, 𝑈, 𝑑 as repeating variables (so that all 3 dimensions are
represented)
Now, 4 𝜋 ----(7-3) parameters are determined as
𝜋1 = 𝜌𝑎1 𝑈𝑏1 𝑑 𝑐1 ∆𝑝
𝜋2 = 𝜌𝑎2 𝑈𝑏2 𝑑𝑐2 𝜇
𝜋3 = 𝜌𝑎3 𝑈𝑏3 𝑑 𝑐3 𝑙
𝜋4 = 𝜌𝑎4 𝑈𝑏4 𝑑 𝑐4 ∈
Now basic units
𝜌 − − − −→ 𝑀𝐿−3
𝑈 − −→ 𝐿𝑇 −1
𝑑 − −−→ 𝐿
∆𝑝 − −−→ 𝑀𝐿−1 𝑇 −2
𝜇 − − − −→ 𝑀𝐿−1 𝑇 −1
∈ − − −→ 𝐿
𝑙 − −−→ . 𝐿
All 𝜋 parameters − − −−→ 𝑀0 𝐿0 𝑇 0
𝑎1 = −1; 𝑏1 = −2 ; 𝐶1 = 0
𝑎2 = −1; 𝑏2 = −1 ; 𝑐2 = −1
𝑎3 = 0; 𝑏3 = 0; 𝑐3 = −1
𝑎4 = 0; 𝑏4 = 0; 𝑐4 = −1
Thus writing 𝜋1 = 𝑓( 𝜋2 , 𝜋3 , 𝜋4 )
∴ The 𝜋 group can be written as follows,
∆𝑝
𝜋1 = 𝑑0 𝑉 −2 𝜌−1 ∆𝑝 = 𝜌𝑉 2 𝑬𝒖 . 𝑵𝒐.
𝜇 𝑑𝑉𝜌
𝜋2 = 𝑑𝑉𝜌 𝑜𝑟 𝑹𝒆 . 𝑵𝒐 .
𝜇
𝑙
𝜋3 = 𝑑−1 𝑉 0 𝜌0 𝑙 = 𝑑
∈
𝜋4 = 𝑑
∴ The new relation can be writing
∆𝑝 𝑑𝑉𝜌 𝑙 ∈
𝑓1 (𝜌𝑉 2 , , 𝑑 , 𝑑) = 0
𝜇
When conclude Δp
∆𝑝 𝑙 ∈
= 𝑓2 (𝑅𝑒 , 𝑑 , 𝑑)
𝜌𝑉 2
𝑊
𝜌= 𝑔
∆𝑝 𝑉2 𝑙 ∈
= 2𝑔 . 𝑓2 (𝑅𝑒 , 𝑑 , 𝑑)
𝑊
The pressure drop is function of (L/d) exponent to(1) in darcy equation
∆𝑝 𝑉2 𝐿 ∈
= 2𝑔 . 𝑑 . 𝑓3 (𝑅𝑒 , 𝑑)
𝑊
Therefore
∆𝑝 𝐿 𝑉2
= (𝐹𝑎𝑐𝑡𝑜𝑟 𝒇) ( ( ))
𝑊 𝑑 2𝑔
Ex.12
Assume the input power to a pump is depend on the fluid weight per
unit volume, flow rate and head produced by the pump. Create a relation by
dimensional analysis between the power and other variables by two methods.
Sol.
Method-1
𝑃 = 𝑓( 𝑊, 𝑄, 𝐻)
𝑃 = 𝐾 𝑊 𝑎 𝑄𝑏 𝐻𝑐
In Dimension analysis
𝐹1 𝐿1 𝑇 −1 = (𝐹𝐿−3 )𝑎 (𝐿3 𝑇 −1 )𝑏 (𝐿)𝑐
Hence
𝑎 = 1 , 1 = −3𝑎 + 3𝑏 + 𝑐
∴ 𝑎 = 1, 𝑏 = 1, 𝑐 = 1
𝑃 = 𝐾𝑊𝑄𝐻
Method-2
𝐹(𝑃, 𝑊, 𝑄, 𝐻) = 0
The variables in dimensions are
P − −−→ FLT −1
Q − −→ L3 T −1
W − −→ FL−3
H − −−→ L
The four variables in 3 fundamental dimensional ∴ 𝜋 𝑔𝑟𝑜𝑢𝑝 𝑖𝑠(4 − 3) = 1
Choice Q, W, H as variable with unknown exponent
∴ 𝜋1 = (𝑄)𝑎1 (𝑊)𝑏1 (𝐻)𝑐1 𝑃
𝑜𝑟 𝜋1 = (𝐿3𝑎1 𝑇 −𝑎1 )(𝐹𝑏1 𝐿−3𝑏1 )(𝐿𝑐1 )(𝐹 𝐿 𝑇 −1 ) = 𝐹 0 𝐿0 𝑇 0
Exponent equality foe F,L,T producing
a1 = −1, b1 = −1, c1 = −1
P
∴ π1 = Q−1 W −1 H −1 P = QWH
𝑃
𝜋1 = 𝐹 (𝑄𝑊𝐻)
𝑃 = 𝐾(𝑄𝑊𝐻)
Ex.13
Assume the input power to a pump is depend on the fluid weight
(W), flow rate (Q) and head produced by pump (H), create a relation by
dimension analysis between the power input and other variables by using
FLT system.
Sol.
Step-1
Quantities Dimensions
Power (P) F L T −1
Flow rate (Q) L3 T −1
Weight(W) per unit volume 𝐹 𝐿−3

Head (H) L
Step-2:- there are four variables & (3) fundamental dimensions

∴ 𝜋 𝑔𝑟𝑜𝑢𝑝 𝑖𝑠 (4 − 3) = 1
Step-3:
Choice Q, W, H as variable with unknown exponent.
𝐹1 (𝑃, 𝑊, 𝑄, 𝐻) = 0−→ 𝐹2 (𝜋1 ) = 0
∴ 𝜋1 = (𝑄)𝑎1 (𝑊)𝑏1 (𝐻)𝑐1 𝑃
Or 𝜋1 = (𝐿3 𝑇 −1 )𝑎1 (𝐹 𝐿−3 )𝑏1 (𝐿)𝐶1 (𝐹𝐿𝑇 −1 )1 = 𝐹 0 𝐿0 𝑇 0
F] 𝑏1 + 1 = 0 − −→ 𝑏1 = −1
L] 3𝑎1 − 3𝑏1 + 𝑐1 + 1 = 0 → 3𝑎1 + 𝑐1 = −4
T] −𝑎1 − 1 = 0 − −→ 𝑎1 = −1
∴ 𝑐1 = −1
𝑃
∴ 𝜋1 = 𝑄 −1 𝑊 −1 𝐻−1 𝑃 = 𝑄𝑊𝐻
𝑃
𝜋1 = 𝑓 (𝑄𝑊𝐻)
𝑃 = 𝐾(𝑄, 𝑊, 𝐻)
Ex.14
Assuming the resistant force for a body submerged in a fluid is function
of (density  , velocity V, viscosity  and characteristic length L). Conclude
a general equation of resistant force by using FLT system.
Sol.
Step.1
Quantities Dimension
Force(F) F
Density() 𝐹 𝐿−4 𝑇 2
Velocity(V) 𝐿 𝑇 −1
Viscosity() 𝐹 𝐿−2 𝑇
Length(L) L
𝐹1 (𝐹, 𝜌, 𝑉, 𝐿, 𝜇) = 0 𝑛 = 5 , 𝑚 = 3
Step-2
We have 5 variables with 3 fundamental dimensions
∴ 𝜋 𝑔𝑟𝑜𝑢𝑝𝑠 = 5 − 3 = 2
We choice 3 repeated variables of unknown exponents
𝜋1 = (𝐿)𝑎1 (𝑉)𝑏1 (𝜌)𝑐1 𝐹 = (𝐿)𝑎1 (𝐿𝑇 −1 )𝑏1 (𝐹 𝑇 2 𝐿−4 )𝑐1 (𝐹) = 𝐹 0 𝐿0 𝑇 0
F] 𝐶1 + 1 = 0 − −→ 𝐶1 = −1
L] 𝑎1 + 𝑏1 − 4𝑐1 = 0 − −→ 𝑎1 + 𝑏1 = −4
T] −𝑏1 + 2𝐶1 = 0−→ 𝑏1 = −2; ∴ 𝑎1 = −2
∴ π1 = L−2 V −2 ρ−1 F = F/L2 V 2 ρ
π2 = (L)a2 (V)b2 (ρ)C2 μ = (L)a2 (𝐿 𝑇 −1 )𝑏2 (𝐹 𝐿−4 𝑇 2 )𝐶2 (𝐹 𝐿−2 𝑇) =
𝐹 0 𝐿0 𝑇 0
F] 𝐶2 + 1 = 0 − −→ 𝐶2 = −1
L] 𝑎2 + 𝑏2 − 4𝐶2 − 2 = 0
a2 + b2 + 2 = −−→ a2 + b2 = −2
T] −𝑏2 + 2𝐶2 + 1 = 0 − −→ 𝑏2 = −1
∴ a2 = −1
𝜇
∴ 𝜋2 = 𝐿−1 𝑉 −1 𝜌−1 𝜇 = 𝐿𝑉𝜌 − −→ 𝜋2−1 = 𝑅𝑒
∴ 𝑓(𝜋1 , 𝜋2 ) = 0
𝐹
𝐹1 (𝜋1 , 𝜋2−1 ) = 0 − −→ 𝐹1 (𝐿2 𝑉 2 𝜌 , 𝑅𝑒 ) = 0
∴ 𝐹 = 𝐿2 𝑉 2 𝜌 𝐹2 (𝑅𝑒 )
Problems.
P6.1 A stationary sphere in water moving at velocity of 1.6 m/s experiences a
drag of 4 N. Another sphere of twice the diameter is placed in a wind
tunnel. Find the velocity of the air and the drag which will give
dynamically similar conditions. The ratio of kinematic viscosities of air
and water is 13, and the density of air is 1.28𝑘𝑔/𝑚3 . [Vair=10.4m/s,
Fd=0.865 N]
P6.2 A 1:80 scale model of an aircraft was tested in air at 20C moves with
speed 40m/s
a- What is the speed model if its submerged in water at 26C,
𝑚2 𝑚2
𝜈𝑎𝑖𝑟 = 14.86 ∗ 10−6 𝑠 ; 𝜈𝑤𝑎𝑡𝑒𝑟 = 0.864 ∗ 10−6 𝑠 [V=2.325 m/s]
b- Determine the air resistance to the prototype if model water resistance
is 7.43N. [Fp=2.643 N]
P6.3 Predicting the general form of input power to a fan , which is depends
on the density, velocity, viscosity of air, diameter, angular velocity of
𝝎𝒅
fan and speed of sound i. [𝑷 = 𝒌 (( 𝑽 ) (𝑹𝒆)(𝑴𝒂)) 𝝆𝑽𝟑 𝒅𝟐 ]
P6.4 Write the equation of displacement (S) for a free failing body with time
(T). Assuming that the displacement depends on weight (W), gravity
(g) and time(T). [S=kgT2]
P6.5 A model spillway has a flow of 100 l/s per meter of width. What is the
actual flow for the prototype spillway if the model scale is 1:20
[qp=8.94 m3/s per m]
P6.6 Researchers plan to test a 1:13 model of a ballistic missile in a high –
speed wind tunnel. The prototype missile will travel at 380m/s through
air at 23C and 95.0 kPa absolute
(a) If the air in the tunnel test section has a temperature -20C at a
pressure of 89 kPa absolute, what its velocity? [Vm=351 m/s]
b) Estimate the drag force on the prototype if the drag force on the
model is 400N. [Fp=72310 N]
P6.7 Derive an expression for the shear stress at the pipe wall when an
incompressible fluid flows through a pipe under pressure. Use diameter
D, flow velocity V, viscosity  and density  of the fluid by using 𝝅
theorem. [𝝉 = 𝝆𝑽𝟐 𝝋(𝑹𝒆)]
P6.8 Derive an expression for the drag on aircraft flying at supersonic speed,
in the form of a function including dimensionless quantities by using
Buckingham's (𝜋) theorem. [𝑭𝒅 = 𝝆𝑳𝟐 𝑽𝟐 𝒇𝟏 (𝑹𝒆, 𝑴𝒂)]
P6.9 Derive an expression for small flow rates over a spillway in the form of
a function including dimensionless quantities. Use dimensional
Analysis with the following parameters height of spillway y, head on
the spillway H, viscosity of liquid , density of liquid , surface tension
𝒉
 and acceleration due to gravity g. [𝒒 = 𝒈𝟏/𝟐 𝑯𝟑/𝟐 𝒇𝟐 (𝒚 , 𝑹𝒆, 𝑾𝒆)]
P6.10 The resisting force F at a plane during flight can be considered as

dependent up on the length of aircraft L velocity V air viscosity, air
density  and black modules of air k. express the functional
relationship between these variables and resisting force using
𝝁 𝒌
dimensional analysis. 𝑭 = 𝑳𝟐 𝑽𝟐 𝝆 ∗ 𝝋 (𝑳𝑽𝝆 , 𝑽𝟐 𝝆)
P6.11 The pressure difference Δp in a pipe of diameter D and length L due to

turbulent flow depends on the velocity V, viscosity , density  and
roughness  using Buckingham's 𝝅 theorem to obtain an expression
𝝁𝑽 𝑳 𝝆𝑽𝑫
for Δp. [∆𝒑 = 𝑫 ∗ 𝑫 ∗ ( 𝝁 )]
P6.12 Prove that the shear stress  in a fluid flowing through a pipe can be
𝜇
expressed by the equation 𝜏 = 𝜌𝑉 2 ∅(𝜌𝐷𝑉 )
Where; D = diameter,  = mass density, V= velocity  = viscosity.
P6.13 A model of a submarine of scale 1/40 is tested in wind tunnel. Find the
speed of air in wind tunnel if the speed of submarine in sea water is
15m/s. Also find the ratio of the resistance between the model and its
prototype. Take the value of kinematic viscosities for sea water and
air as 0.012 stokes and 0.016 stokes respectively. The weight density
of sea water and air are given as 10.1kN/m3 and 0.0122 kN/m3
𝒎 𝑹
respectively. 𝑽𝒎 = 𝟖𝟎𝟎 𝒔 , 𝑹𝒎 = 𝟎. 𝟎𝟎𝟐𝟏𝟒
𝒑
P6.14 A spillway model is to be built to a geometrically similar scale of 1/50

across a flow of 600m width. The prototype is 15m high and
maximum head on it is expected to be 1.5m.
i) What height of model and what head on model should be used.
[Hm=0.3 m]
ii) If flow over the model for a particular head is 12 L/s what flow per
𝑸𝒑 𝒍𝒊𝒕
meter length at prototype is expected. [ 𝑳 = 𝟕𝟎𝟕𝟏 𝒔 ]
𝒑
P6.15 In an airplane model of size 1/50 of its prototype the pressure drop is
4 bar. The model is tested in water. Find the corresponding pressure
drop in the prototype. Take specific weight of air = 0.00124kN/m3. The
viscosity of water as 0.01 poise while the viscosity of air is 0.00018
poise. [p=0.0004114 bar]
2
P6.16 An oil of S.G. 0.9 and viscosity 0.003 N.s/m is to be transported at the
rate of 3000 L/s through a 1.5m diameter pipe. Test was conducted on a
15cm diameter pipe. Using water at 20𝐶 ° if the viscosity at 20𝐶 ° is
0.001 N.s/m2 fined.
i) Velocity of flow in the model. [Vm=5.1 m/s]
ii) Rate at flow in the model. [Q=80.9 lit/s]
P6.17 A geometrically similar model of an air duct is built to 1/25 scale and
tested with water which is 50 times more viscous and 800 times than
air. When tested under dynamically similarity conditions the pressure
drop is 2 bar in the model. Find the corresponding pressure drop in
full scale prototype. [1.024*10-3bar]
P6.18 In a geometrically similar model of spillway the discharge per length is

0.2m3/s. If the scale of the model is 1/36, find the discharge per meter
run of the prototype. [qp=43 m3/s]
P6.19 The force required to tow a 1:30 scale model of a modern boat in a
lake at a speed of 2m/s is 0.5N. Assuming that the viscous resistance
due to water and air is negligible in comparison with the wave
resistance calculate the corresponding speed of the prototype for
dynamically similar conditions. What would be the force required to
propel the prototype at that velocity in the same lack? [Fp=13500 N]
P6.20 In an airplane model at size 1/40 of its prototype the pressure drop is
7.5 kN/m2 the model is tested in water. Find the corresponding
pressure drop in the prototype. Take density of air = 1.24 kg/m3,
density of water = 100 kg/m3, Viscosity of air = 0.00018 poise,
Viscosity of water = 0.01 poise. [pp=1.225 N/m2]
Chapter 7 Viscous Incompressible Flows in Pipes 155
CHAPTER 7
Viscous Incompressible Flows
in Pipes
Part-One
(Laminar Flow)
7.1 Introduction.
Real fluids possess viscosity, while ideal fluid is inviscid. The viscosity
of fluid introduce resistance to motion by developing shear or frictional stress
between the fluid layers and between fluid layers and the boundary, which
causes the real fluid to a adhere to the solid boundary and hence no relative
motion between fluid layer and solid boundary.
Viscosity causes the flow to occur in two modes laminar and turbulent flow.
Reynolds number < 2000, the flow is always laminar through a pipe which is
critical value of Re for circular pipe. Flow between parallel plates based on
mean velocity and distance between the plates.
𝑅𝑒 = 𝑣𝑖𝑠𝑐𝑜𝑢𝑠𝑓𝑜𝑟𝑐𝑒
The flow is laminar when one of the conditions occurs
i) Viscosity is very high.
ii) Velocity is very low.
iii) The passage is very narrow.
7.2 Relationship between Shear Stress and Pressure Gradient.

The shear stress is maximum at the boundary and gradually decreases
with increase in distance from the solid boundary where the velocity is zero
at the boundary. A pressure gradient exists which overcome the shear
resistance and causes the fluid to flow. Due to non uniform distribution of
velocity, the fluid at any layer moves at a higher velocity than the layer
below.
The motion of the fluid element will be resisted by shearing or frictional
force which must be overcome by maintaining a pressure gradient in the
direction of flow, from fig 7.1,
Let  = shear stress on the lower face ABCD of the element

𝜕𝜏
𝜏 + 𝜕𝑦 𝛿𝑦 = Shear stress on the upper face 𝐴́𝐵́ 𝐶́ 𝐷́ of the element.
For 2-dimensional steady flow there will be no shear stress on the
vertical faces ABB'A' & CDD'C' as in Fig. 7.1. Thus the only forces acting
on the element in the direction of flow (x-axis) will be the pressure and shear
forces. Let 𝛿𝑥, 𝛿𝑦 𝑎𝑛𝑑 𝛿𝑧 are element thickness in x, y and z directions.
Net shearing force on the element in y-direction is equal to
𝜕𝜏 𝜕𝜏
= (𝜏 + 𝜕𝑦 𝛿𝑦) 𝛿𝑥. 𝛿𝑧 − 𝜏. 𝛿𝑥𝛿𝑧 = 𝜕𝑦 𝛿𝑥. 𝛿𝑦. 𝛿𝑧 (7.1)
Net pressure force on the element in x-direction is equal to
𝜕𝑝 𝜕𝑝
= 𝑝. 𝛿𝑦. 𝛿𝑧 − (𝑝 + 𝜕𝑥 𝛿𝑥) 𝛿𝑦. 𝛿𝑧 = − 𝜕𝑥 𝛿𝑥. 𝛿𝑦. 𝛿𝑧 (7.2)
For the flow to be steady and uniform these begin no acceleration, the sum of
the forces must be zero, from (7.1 & 7.2)
𝜕𝜏 𝜕𝑝
. 𝛿𝑥. 𝛿𝑦. 𝛿𝑧 − 𝜕𝑥 . 𝛿𝑥. 𝛿𝑦. 𝛿𝑧 = 0
𝜕𝑦
The relationship between shear stress and pressure gradient is
𝜕𝑝 𝜕𝜏
= 𝜕𝑦 (7.3)
𝜕𝑥
Eq. (7.3) indicates that the pressure gradient in the direction of flow is equal
to the shear gradient in the direction normal to the direction of flow. This is
applicable for laminar and turbulent flow.
𝜕𝜏
𝜏+ 𝛿𝑦
` 𝜕𝑦
D' C'
y 𝜕𝑝
B' 𝑝+ 𝛿𝑥
A'
p 𝜕𝑥
D C
𝜏
A B
z
Figure 7.1: Pressure and Shear stress Forces on a Fluid Element.
7.3 Laminar Flow between Parallel Plates.

Consider two parallel plates with (h) distance apart. For steady flow
between them a pressure gradient 𝜕𝑝/𝜕𝑥 exist which related shear stress in y-
direction.
𝑑𝑢
Since 𝜏 = 𝜇 𝑑𝑦
𝜕2 𝑢 1 𝜕𝑝
Then 𝜕𝑦 2 = 𝜇 𝜕𝑥 y
Integration gives
1 𝜕𝑝
𝑢 = 2𝜇 (𝜕𝑥 ) 𝑦 2 + 𝐶1 𝑦 + 𝐶2 Umax
𝑢 = 0 𝑎𝑡 𝑦 = 0 𝑎𝑛𝑑 𝑦 = ℎ
h
1 𝜕𝑝 . x
𝐶2 = 0 , 𝐶1 = − 2𝜇 𝜕𝑥 ℎ
And Figure 7.2: Laminar flow between parallel plates.
1 𝜕𝑝 2
𝑢 = 2𝜇 (− 𝜕𝑥 ) (ℎ𝑦 − 𝑦 )
This equation is a parabola shape with vertex at center line (y=h/2) when
the maximum velocity occurs
1 𝜕𝑝 ℎ2 ℎ2
𝑈𝑚𝑎𝑥 = 2𝜇 (− 𝜕𝑥 ) ( 2 − )
4
1 𝜕𝑝
Or𝑢𝑚𝑎𝑥 = 8𝜇 (− 𝜕𝑥 ) ℎ2 (7.4)
The (-ve) of pressure gradient is the pressure drop in the direction of flow.
The discharge dq through a small area of depth dy per unit width is
dq = udy
1 𝜕𝑝
𝑑𝑞 = (− ) (ℎ𝑦 − 𝑦 2 )𝑑𝑦
2𝜇 𝜕𝑥
ℎ 1 𝜕𝑝
𝑄= ∫0 2𝜇 (− 𝜕𝑥 ) (ℎ𝑦 − 𝑦 2 )𝑑𝑦
ℎ
1 𝜕𝑝 𝑦2 𝑦3
𝑄 = 2𝜇 (− 𝜕𝑥 (ℎ − ))
2 3
0
1 𝜕𝑝
𝑄 = 12𝜇 (− 𝜕𝑥 ) ℎ3 (7.5)
𝑄 𝑄
Mean velocity of flow u̅ = 𝑓𝑙𝑜𝑤 𝑎𝑟𝑒𝑎 = ℎ∗1
1 𝜕𝑝
(− )ℎ3 1 𝜕𝑝
12𝜇 𝜕𝑥
u̅ = = 12𝜇 (− 𝜕𝑥 ) ℎ2 (7.6)
ℎ
The above velocity (𝑢̅) may be used to calculate the pressure drop
12𝜇𝑢̅ 𝜕𝑝 ̅
12𝜇𝑢
−𝜕𝑝 = ℎ2 𝜕𝑥 → 𝜕𝑥 = ℎ2 (7.7)
From Eq's (7.4 & 7.6)
3
𝑈𝑚𝑎𝑥 = 2 𝑢̅
The pressure drop between two sections with distance 𝑥1 𝑎𝑛𝑑 𝑥2 from origin
is
𝑝2 𝑥2 12𝜇𝑢
̅
∫𝑝 (−𝑑𝑝) = ∫𝑥 ℎ2 𝑑𝑥
1 1
̅
12𝜇𝑢
𝑝1 − 𝑝2 = (𝑥2 − 𝑥1 )
ℎ2
If L is the length between the sections

̅𝐿
12𝜇𝑢
𝑝1 − 𝑝2 = ℎ2 (7.8)
The variation of shear stress in the y-direction is
𝜕 1 𝜕𝑝
𝜏 = 𝜇 [ (− ) (ℎ𝑦 − 𝑦 2 )]
𝑑𝑦 2𝜇 𝜕𝑥
1 𝜕𝑝
𝜏 = 2 (− 𝜕𝑥 ) (ℎ − 2𝑦 )
𝜕𝑝 ℎ
𝜏 = − 𝜕𝑥 (2 − 𝑦 ) (7.9)
Shearing stress varies linearly with y, it is maximum at the boundary, y=0
and y=h
ℎ 𝜕𝑝 ℎ 12𝜇 ̅
6𝜇𝑢
𝑎𝑡 𝑦 = 0 𝜏0 = 2 [ − 𝜕𝑥 ] = 2 ( ℎ2 𝑢̅) = (7.10)
ℎ
ℎ 𝜕𝑝 ℎ 12𝜇 ̅
6𝜇𝑢
𝐴𝑡 𝑦 = ℎ 𝜏0 = − 2 (− 𝜕𝑥 ) = − 2 ( ℎ2 𝑢̅) = − ℎ (7.11)
Shearing stress at the center y = h/2 is zero
𝜕𝑝 ℎ ℎ
𝜏 = − 𝜕𝑥 (2 − 2 ) = 0
Ex.1
Incompressible fluid flows through a rectangular passage of width (b),
small depth (t) and length L, in the direction of flow. If the pressure drop
between the two ends is p calculate the shear stress at the wall of the passage
in terms of mean velocity and the coefficient of viscosity
Sol.
̅𝐿
12𝜇𝑢
𝑝 = 𝑡2
𝑡 𝜕𝑝 𝑡𝑝
𝜏0 = 2 (− 𝜕𝑥 ) = 2𝐿
̅
6𝜇𝑢
Then 𝜏0 = 𝑡
Ex.2
Water at 20C flows between two large parallel plates separated by a
distance of 16mm. calculate
i) Max. velocity
ii) Shear stress at the wall if the average velocity is (0.4 m/s) (take 
for water = 0.01 Poise)
Sol.
𝑡2 𝜕𝑝 3 3 𝑚
i) 𝑈𝑚𝑎𝑥 = 8𝜇 (− 𝜕𝑥 ) = 2 𝑢̅ = 2 ∗ 0.4 = 0.6 𝑠
ii) Shear stress at the wall
6𝜇𝑢̅
𝜏0 = − 𝑡 = −0.15 𝑛/𝑚2
7.4 Couette Flow.
Couette flow is the flow between two parallel plates as in Fig. 7.3
one plate is at rest and the other is moving with a velocity U, assuming
infinitely large in z-direction
The governing Equation is

𝑑𝑝 𝑑2 𝑢
= 𝜇 𝑑𝑦 2
𝑑𝑥
Flow is independent of any variation in z-direction the boundary condition
are
i) At y = 0 , u = 0 y
ii) At y = h , u = 𝑈 Moving plate
After integration twice we get U
1 𝑑𝑝 2
𝑢 = 2𝜇 𝑑𝑥 𝑦 + 𝑐1 𝑦 + 𝑐2
h
At y = 0, u = 0, then 𝑐2 = 0 x
1 𝑑𝑝
𝑢 = 2𝜇 𝑑𝑥 𝑦 2 + 𝑐1 𝑦
Fixed plate
At y=h , u= 𝑈
𝑈 1 dp Figure 7.3: Couette flow between parallel plates
∴ 𝐶1 = h − 2μ dx h
Then the expression for u becomes
𝑦 1 𝑑𝑝
𝑢 = ℎ 𝑈 − 2𝜇 𝑑𝑥 (ℎ𝑦 − 𝑦 2 ) (7.12)
Multiply and divided by(ℎ2 )
𝑦 ℎ2 𝑑𝑝 𝑦 𝑦
Or 𝑢 = ℎ 𝑈 − 2𝜇 . 𝑑𝑥 . ℎ (1 − ℎ ) (7.13)
Eq. 7.13 can also be expressed in the form
𝑢 𝑦 ℎ2 𝑑𝑝 𝑦 𝑦
= ℎ − 2𝜇𝑈 . 𝑑𝑥 . ℎ (1 − ℎ)
𝑈
𝑢 𝑦 𝑦 𝑦
Or = ℎ + 𝑃 ℎ (1 − ℎ) (7.14)
𝑈
ℎ2 𝑑𝑝
Where 𝑃 = − 2𝜇𝑈 (𝑑𝑥 )
P is known as the non- dimensional pressure gradient. When P=0, the
velocity distribution across the channel reduced to
𝑢 𝑦
𝑈
= ℎ 𝑖𝑠 𝑘𝑛𝑜𝑤𝑛 𝑎𝑠 𝑠𝑖𝑚𝑝𝑙𝑒 𝑐𝑜𝑢𝑒𝑡𝑡𝑒 𝑓𝑙𝑜𝑤 .
 When P>0 , i.e for a negative pressure gradient (-dp/dx ) in the
direction of motion , the velocity is positive over the whole gap.
 When P < 0, these is positive or adverse pressure gradient in the
direction of motion and the velocity over a portion of channel width
can become negative and back flow may occur near the wall, which is
at rest
7.4.1 Maximum and Minimum Velocities.

The variation of maximum and minimum velocity in the channel is
found out by setting du/dy = 0 from Eq. 7.14, we can write
𝑑𝑢 𝑈 𝑃𝑈 𝑦
= ℎ + ℎ (1 − 2 ℎ)
𝑑𝑦
Setting du/dy = 0 gives
𝑦 1 1
= 2 + 2𝑃 (7.15)
ℎ
By studying Eq. 7.15 we conclude that
1- The maximum velocity for P=1 occurs at y/h =1 and equal to U. For P>1
, the maximum velocity occurs at a location y/h <1.
2- i.e that with P>1, the fluid particles attain a velocity higher than that of
the moving plate.
3- for P= -1 , the minimum velocity occurs at y/h =0 for P< -1 , the minimum
velocity occurs at a location y/h>1
4- This means that these occurs a back flow near the fixed plate. The values
of maximum and minimum velocities can be determined by substituting
the value of y from Eq. 7.15 into Eq. 7.14 as
𝑈(1+𝑃)2
𝑈𝑚𝑎𝑥 = for P ≥ 1
4P
(7.16)
𝑈(1+𝑃)2
𝑈𝑚𝑖𝑛 = 4P 𝑓𝑜𝑟 𝑃 ≤ 1
The expression for shear stress can be obtained by substituting the value of u
in Newton's equation of viscosity
𝑑𝑢 𝑑 𝑦𝑈 1 𝑑𝑝
𝜏= 𝜇 = 𝜇 [ − (ℎ𝑦 − 𝑦 2 )]
𝑑𝑦 𝑑𝑦 ℎ 2𝜇 𝑑𝑥
𝑈 𝑑𝑝 ℎ
𝜏 = 𝜇 ℎ + (− 𝑑𝑥 ) (2 − 𝑦) (7.17)
The shear stress at the center is
𝑈
𝜏 = 𝜇ℎ (7.18)
Ex.3
Laminar flow takes place between parallel plates 10 mm apart. The plates
are inclined at 45 with the horizontal. For oil of viscosity 0.9 kg/m.s and
mass density is 1260 kg/m3, the pressure at two points 1.0 m vertically apart
are 80 𝑘𝑁/𝑚2 and 250 𝑘𝑁/𝑚2 when the upper plate moves at 2.00 m/s
velocity relative to the lower plate but in opposite direction to flow determine
i) velocity distribution
ii) max. velocity Flow
iii) shear stress on the top plate p1=250 kPa
10 mm Plate V=2m/s
Sol.
Consider section 1&2 from Bernoulli's Eqn. p2=80 kPa
𝑝1 𝑝2 1.0 m
H1 − H2 = − ( 𝛾 + 𝑍1 ) + ( 𝛾 + 𝑍2 )
250000 80000
= − (9.806∗1260 + 1) + (9.806∗1260 + 0)
𝐻1 − 𝐻2 = −21.234 + 6.475 = −14.759𝑚 𝑖𝑛 1.414𝑚 𝑙𝑒𝑛𝑔𝑡ℎ
Since 𝐻1 is greater than 𝐻2 , flow will be in down word direction .
𝜕𝐻 14.759
= − 1.414 = −10.438
𝜕𝑥
𝑘𝑁⁄
𝜕𝑝 𝜕𝐻 𝑚2
And 𝜕𝑥 = 𝛾 𝜕𝑥 = −10.438 ∗ 1260 ∗ 9.806 = −128.97 𝑚
𝑢 𝑦 1 𝜕𝑝
= ℎ − 2𝜇𝑈 𝜕𝑥 (𝑦ℎ − 𝑦 2 )
𝑈
𝑚 0.9𝑘𝑔
𝑈 = −2 , ℎ = 0.01𝑚, 𝜇 =
𝑠 𝑚𝑠
2 1
∴ 𝑢 = − 0.01 𝑦 − 2∗0.9 (−128967.33)(0.01𝑦 − 𝑦 2 )
i) 𝑢 = 516.4364 𝑦 − 71648.5𝑦 2
To find y at which u is max. set du/dy=0 = 516.486 − 143297.2𝑦 or
𝑦 = 3.604 ∗ 10−3 𝑚
ii) ∴ 𝑢𝑚𝑎𝑥 = (516.486 ∗ 0.003604) − (71648.2 ∗ 0.0036042 ) =
𝑚
0.9308 𝑠
𝜕𝑢
iii) 𝜏0 = 𝜇 (𝜕𝑦)𝑦=0.01 = 0.9(516.486 − 143297.2 ∗ 0.01) =
−824.837 𝑁/𝑚2
7.5 Pipe of Circular Cross-Section.

7.5.1 Hagen-Poiseuille Flow.
Consider fully developed laminar flow through a straight tube of circular
cross – section as in Fig. 7.4 . Rotational symmetry is considered to make the
flow two – dimensional axisymmetry. Let us take x-axis as the axial of the
tube along which all the fluid particles travel, i.e.
𝑉𝑥 ≠ 0, 𝑉𝑟 = 0, 𝑉𝜃 = 0
Now from continuity equation, we obtain
𝜕𝑉𝑟 𝑉 𝜕𝑉 1 𝜕𝑣
+ 𝑟𝑟 + 𝜕𝑥𝑥 = 0 [𝑓𝑜𝑟 𝑟𝑜𝑡𝑎𝑡𝑖𝑜𝑛𝑎𝑙 𝑠𝑦𝑚𝑚𝑒𝑡𝑟𝑦, 𝑟 . 𝜃𝜃 = 0]
𝜕𝑟
This means 𝑉𝑥 = 𝑉𝑥 (𝑟, 𝑡)
𝜕𝑉 𝜕
Invoking [𝑉𝑟 = 0, 𝑉𝜃 = 0 𝜕𝑥𝑥 = 0, 𝑎𝑛𝑑 𝜕𝜃 ( 𝑎𝑛𝑦 𝑞𝑢𝑎𝑛𝑡𝑖𝑡𝑛𝑔) = 0]
With Navier–Stokes equation, we obtain in the x-direction
𝜕𝑉𝑥 1 𝜕𝑝 𝜕2 𝑉 1 𝜕𝑉𝑥
= − 𝜌 . 𝜕𝑥 + 𝜈 ( 𝜕𝑟 2𝑥 + 𝑟 . ) (7.19)
𝜕𝑡 𝜕𝑟
For steady flow, the governing equation becomes
𝜕2 𝑉𝑥 1 𝜕𝑉𝑥 1 𝑑𝑝
+ . = (7.20)
𝜕𝑟 2 𝑟 𝜕𝑟 𝜇 𝑑𝑥
The boundary conditions are
𝜕𝑉
i) At 𝑟 = 0 , 𝑉𝑥 𝑖𝑠 𝑓𝑖𝑛𝑖𝑡 & 𝜕𝑟𝑥 = 0
ii) At r = R, Vx=0 yield Eq. 7.20 can be written after multiplying by r
𝑑2 𝑉𝑥 𝑑𝑉𝑥 1 𝑑𝑝
𝑟 + = 𝜇 . 𝑑𝑥 𝑟
𝑑𝑟 2 𝑑𝑟
𝑑 𝑑𝑉𝑥 1 𝑑𝑝
𝑜𝑟 (𝑟 ) = 𝜇 𝑑𝑥 𝑟 by integration
𝑑𝑟 𝑑𝑟
𝑑𝑉𝑥 1 𝑑𝑝
𝑟 = 2𝜇 . 𝑑𝑥 𝑟 2 + 𝐴
𝑑𝑟
𝑑𝑉𝑥 1 𝑑𝑝 𝐴
= 2𝜇 . 𝑑𝑥 𝑟 + 𝑟 by integration
𝑑𝑟
1 𝑑𝑝
𝑉𝑥 = 4𝜇 . 𝑑𝑥 𝑟 2 + 𝐴 𝑙𝑛 𝑟 + 𝐵
𝑑𝑉𝑥
At 𝑟 = 0 𝑉𝑥 = 𝑓𝑖𝑛𝑖𝑡𝑒 & =0→𝐴=0
𝑑𝑟
𝑎𝑡 𝑟 = 𝑅 , 𝑉𝑥 = 0
1 𝑑𝑝
𝐵 = − 4𝜇 . 𝑑𝑥 . 𝑅 2
𝑅2 𝑑𝑝 𝑟2
∴ 𝑉𝑥 = 4𝜇 (− 𝑑𝑥 ) [1 − 𝑅2 ] (7.21)
This shows that the axial velocity profile in a fully developed laminar
pipe flow is having parabolic variation along r.
At 𝑟 = 0 , 𝑎𝑠 𝑠𝑢𝑐ℎ , 𝑉𝑥 = 𝑉𝑥 𝑚𝑎𝑥
𝑅2 𝑑𝑝
𝑉𝑥 𝑚𝑎𝑥 = 4𝜇 (− 𝑑𝑥 ) (7.22)
r R
x Vmax
Figure 7.4: Flow in circular pipe.

7.5.2 Volumetric Flow Rate.
The average velocity in pipe is
𝑹
𝑄 ∫𝟎 2πr Vx (r)dr
𝑉𝑎𝑣. = 𝜋𝑅2 = substitute Eq. 7.21
πR2
2𝜋𝑅2 𝑑𝑝 𝑅2 𝑅4
(− )[ − 2 ]
4𝜇 𝑑𝑥 2 4𝑅
𝑜𝑟 𝑉𝑎𝑣 . = πR2
𝑅2 𝑑𝑝 1
𝑉𝑎𝑣. = 8𝜇 (− 𝑑𝑥 ) = 2 𝑉𝑥 max → 𝑉𝑥 𝑚𝑎𝑥 = 2𝑉𝑎𝑣 (7.23)
Now, the discharge Q through a pipe is given by
𝑄 = 𝜋𝑅 2 𝑉𝑎𝑣 (7.24)
𝑅2 𝑑𝑝
𝑄= 𝜋𝑅 2 8𝜇 (− 𝑑𝑥 )
𝜋𝑑4 𝑑𝑝
𝑜𝑟 𝑄 = − 128𝜇 (𝑑𝑥 ) (7.25)
From Eq's 7.22 & 7.23
𝑝1 − 𝑝2 𝜇 𝑉𝑎𝑣
= 4 𝑉𝑚𝑎𝑥 2
= 32𝜇 (7.26)
𝐿 𝑅 𝑑2
Eq. 7.26 is known as the Hagen- Poiseuille equation.
Ex.4
Oil mass density is 800 kg/m3 and dynamic viscosity is 0.002 kg/m.s
flow through 50mm diameter, pipe length is 500 m and the discharge flow
rate is 0.19*10-3 m3/s determine
i) Reynolds number of flow.
ii) Center line velocity.
iii) Loss of pressure in 500 m length.
iv) Pressure gradient.
v) Wall shear stress.
Sol.
4𝑄 4∗0.19∗10−3 𝑚
𝑉𝑎𝑣. = 𝜋𝑑2 = = 0.0968
𝜋∗(0.05)2 𝑠
𝑉𝑑𝜌 0.0968∗0.05∗800
i) 𝑅𝑒 = = = 1936.0
𝜇 0.002
𝑚
ii) 𝑉𝑥 𝑚𝑎𝑥 = 2𝑉𝑎𝑣. = 2 ∗ 0.0968 = 0.1936 𝑠
iii) From Eq. 7.26
𝑝1 − 𝑝2 𝜇 𝑉
= 4 𝑉𝑚𝑎𝑥 𝑅2 = 32𝜇 𝑑𝑎𝑣2
𝐿
32𝜇𝑉𝑎𝑣 𝐿 32∗0.002∗0.0968∗ 500 𝑁
∴ 𝑝1 − 𝑝2 = = = 1239.04 𝑚2
𝑑2 (0.05)2
𝑁
𝑑𝑝 𝑝1 −𝑝2 1239.04 2.478 2
𝑚
iv) = = = = 2.478 𝑝𝑎/𝑚
𝑑𝐿 𝐿 500 𝑚
(𝑝1 −𝑝2 )𝑑 0.05 𝑁
v) 𝜏0 = = (1239.04) ∗ = 0.03098 , 𝐸𝑞. 7.28
4𝐿 4∗500 𝑚2
7.5.3 Shear Stress in Horizontal Pipe.

A force balance for steady flow in horizontal pipe as in Fig. 7.5 yields
𝑝1 (𝜋𝑟 2 ) − 𝑝2 (𝜋𝑟 2 ) − 𝜏(2𝜋𝑟𝐿) = 0
(𝑝 −𝑝 )𝑟
𝑜𝑟 𝜏 = 1 2𝐿 2 (7.27)
From Eq. 7.27
𝑎𝑡 𝑟 = 0 𝜏 = 0
𝑟 = 𝑅 𝜏 = 𝜏0
(𝑝 −𝑝 )𝑑
𝜏0 = 1 4𝐿 2 (7.28)
Eq. 7.27 is valid for laminar & turbulent flow.
𝑝 −𝑝
( 1𝜌𝑔 2 ) Represent the energy drop per unit weight (ℎ𝐿 ) multiply Eq. 7.27 by
(g/g) yields
𝜌𝑔𝑟 𝑝 −𝑝 𝜌𝑔ℎ
𝜏 = 2𝐿 ( 1𝜌𝑔 2 ) = 2𝐿 𝐿 𝑟 (7.29)
2𝜏0 𝐿 4𝜏0 𝐿
∴ ℎ𝐿 = = (7.30)
𝜌𝑔𝑅 𝜌𝑔𝑑
𝜏 = 𝜏0 𝑎𝑡 𝑟 = 𝑅
τo
τAs R
r
p1A p2A Vmax
Figure 7.5: Forces on element in horizontal pipe.
7.5.4 Shear Stress in Inclined Pipe.

The energy equation may be written in pipe when related the loss to
available energy reduction
𝑝1 𝑉2 𝑝 𝑉2
1
+ 2𝑔 + 𝑧1 = 𝜌𝑔2 + 2𝑔
2
+ 𝑧2 + ℎ𝑓
𝜌𝑔 1−2
𝑉2
Since the velocity head (2𝑔) is the same
𝑝1 −𝑝2
ℎ𝑓 = + 𝑧1 − 𝑧2 (7.31)
𝜌𝑔
∆𝑝
∴ ℎ𝑓 = 𝜌𝑔 + ∆𝑧 (7.32)
Applying the linear – momentum eqn. in the L-direction
∑ 𝐹𝑙 = 0 = (𝑝1 − 𝑝2 )𝐴 + 𝛾𝐴𝐿 𝑠𝑖𝑛𝜃 − 𝜏0 𝐿𝑃 = 𝑚̇(𝑉2 − 𝑉1 ) = 0
(P) is the wetted perimeter of the conduit ,i.e , the portion of the perimeter
where the wall is in contact with the fluid when the conduit not circular pipe.
𝐿 𝑠𝑖𝑛𝜃 = 𝑧1 − 𝑧2
𝑝1 − 𝑝2 τ0 𝐿𝑃
+ 𝑧1 − 𝑧2 = (7.33)
𝜌𝑔 𝜌𝑔𝐴
From Eq. 7.31& 7.33
𝜏0 𝐿𝑃
ℎ𝑓 = 𝜌𝑔𝐴 (7.34)
From experiment
𝜌
𝜏0 =  2 𝑉 2 (7.35)
𝜌 𝐿𝑃 𝐿 𝑉2
∴ ℎ𝑓 =  2 𝑉 2 𝛾𝐴 =  𝑅 2𝑔 (7.36)
Rh=A/P
Rh= hydraulic Radius of the conduit
For a pipe Rh=d/4 ; =f/4
Where  is the non-dimensional factor, the ℎ𝑓 head loss due to friction can be
written as follows,
𝑓 𝐿 4 𝑉2 𝐿 𝑉2
∴ ℎ𝑓 = 4 = 𝑓 𝑑 2𝑔 (7.37)
𝑑 2𝑔
Eq. 7.37 is the Darcy – Weisbach equation, valid for duct flows of any cross-
section and for laminar and turbulent flow, f is the friction factor f=4 
By equating Eq's 7.30 & 7.37

4𝜏0 𝐿 𝐿 𝑉2
= 𝑓 𝑑 2𝑔
𝜌𝑔𝑑
𝑓𝜌𝑉 2
∴ 𝜏0 = 8 (7.38)
In Hagen-Poiseuille eqn.
∆𝑝𝑑2
𝑉𝑎𝑣 = From Eq. 7.26
32𝜇𝐿
∆𝑝
∆𝑝 = 𝜌𝑔 ℎ𝑓 − −−→ ℎ𝑓 = 𝜌𝑔
𝜌𝑔 ℎ𝑓 𝑑2
∴ 𝑉𝑎𝑣 = 32𝜇𝐿
32𝑉𝑎𝑣 𝜇𝐿 𝐿 𝑉2
ℎ𝑓 = = 𝑓 𝑑 2𝑔
𝜌𝑔𝑑2
64
64𝑉𝑎𝑣 𝜇𝐿 𝜌𝑑𝑉𝑎𝑣 𝐿 𝑉𝑎𝑣 2 64 𝐿 𝑉𝑎𝑣 2
= ( 2𝜌𝑔𝑑2 ) = =𝑅
𝜇 𝑑 2𝑔 𝑒 𝑑 2𝑔
𝐿 𝑉𝑎𝑣 2 64 𝐿 𝑉𝑎𝑣 2
ℎ𝑓 = 𝑓 𝑑 =𝑅 (7.39)
2𝑔 𝑒 𝑑 2𝑔
64
∴ 𝑓 = 𝑅𝑒 (7.40)
It applies to all roughness and may be used for the solution of laminar flow
problems in pipes.
From above equations the laminar head loss as followes
2
64 𝐿 𝑉𝑎𝑣 32𝜇𝐿𝑉𝑎𝑣 128𝜇𝐿𝑄
ℎ𝑓(𝑙𝑎𝑚𝑖𝑛𝑎𝑟) = 𝑅𝑒 = = (7.41)
𝑑 2𝑔 𝜌𝑔𝑑2 𝜋𝜌𝑔𝑑4
From Eq. 7.22
4𝑉 𝜇𝐿 32𝑉 𝜇𝐿
𝑝1 − 𝑝2 = 𝑚𝑎𝑥 = 𝑑𝑎𝑣2
𝑅2
Pressure drop per unit weight
∆𝑝 32𝜇𝐿𝑉
ℎ𝑓 = 𝜌𝑔 = 𝜌𝑔 𝑑2𝑎𝑣 for laminar flow (7.42)
Ex.5
An oil of viscosity 0.9 𝑁𝑠/𝑚2 and S.G. 0.9 is flowing through a
horizontal pipe of 60 mm diameter. If the pressure drop in 100 m length of
the pipe is 1800𝑘𝑁/𝑚2 , determine:
(i) The rate of flow of oil.

(ii) The center-line velocity.
(iii) The total friction drags over 100 m length.
(iv) The power required to maintain the flow.
(v) The velocity gradient at the pipe wall.
(vi) the velocity and shear stress at 8mm from the wall
Sol.
Area of the pipe,
𝜋
𝐴 = 4 ∗ (0,06)2 = 2.827 ∗ 10−3 (𝑚2 ) Pressure drop in (100m) length of the
pipe, ∆𝑝 = 1800 𝑘𝑁/𝑚2
i) the rate of flow,Q
32𝜇𝑉 𝐿
𝑝1 − 𝑝2 = ∆𝑝 = 𝑑2𝑎𝑣
∆𝑝 𝑑2
𝑉𝑎𝑣 = 32𝜇𝐿
1800∗103 ∗(0.06)2 𝑚
∴ 𝑉𝑎𝑣 = = 2.25
32∗0.9∗100 𝑠
𝜌𝑉𝑑 0.9∗1000∗2.25∗ 0.06
Reynolds number, 𝑅𝑒 = = = 135
𝜇 0.9
As Re is less than 2000, the flow is laminar and the rate of flow is,
𝑚3 𝑙𝑖𝑡
𝑄 = 𝐴 ∗ 𝑉𝑎𝑣 = 2.827 ∗ 10−3 ∗ 2.25 = 6.36 ∗ 10−3 = 6.36
𝑆 𝑠
ii) the center-line velocity , 𝑉𝑚𝑎𝑥
𝑚
𝑉𝑚𝑎𝑥 = 2𝑉𝑎𝑣 = 2 ∗ 2.25 = 4.5 𝑠
iii) the total frictional drag over (100m) length
(𝑝 −𝑝 )𝑑
From 𝜏0 = 1 4𝐿 2
0.06
∴ 𝜏0 = 1800 ∗ 103 ∗ 4∗100 = 270 𝑁/𝑚2
∴ Friction drag for (100m) length
𝐹𝑑 = 𝜏0 ∗ 𝐴𝑠 = 𝜏0 ∗ 𝜋𝑑𝐿 = 270 ∗ 𝜋 ∗ 0.06 ∗ 100
𝐹𝑑 = 5089 𝑁
(iv) The power required to maintain the flow, P,
𝑃 = 𝐹𝑑 ∗ 𝑉𝑎𝑣 = 5089 ∗ 2.25 = 11451 𝑊
= 15.35 h.p
Alternatively,
𝑃 = 𝑄. ∆𝑝 = 0.00636 ∗ 1800 ∗ 103 = 11448 𝑊
𝑑𝑢
(v) The velocity gradient at the pipe wall, (𝑑𝑦) ;
𝑦=0
𝜕𝑢
𝜏0 = 𝜇. (𝜕𝑦)
𝑦=0
𝜕𝑢 𝜏0 270
𝑜𝑟 (𝜕𝑦) = = = 300 𝑠 −1
𝑦=0 𝜇 0.9
(vi) the velocity and shear stress at (8mm) from the wall,
𝑅2 𝜕𝑝 𝑟2
𝑉 = 4𝜇 (− 𝜕𝑥 ) (1 − 𝑅2 )
1 𝜕𝑝
Or 𝑉 = − 4𝜇 . 𝜕𝑥 (𝑅 2 − 𝑟 2 )
Here, 𝑦 = 8𝑚𝑚 = 0.008𝑚
But y = R-r
∴ 0.008 = 0.03 − 𝑟 − −−→ 𝑟 = 0.022𝑚
1 1800∗ 103 𝑚
∴ 𝑉(8𝑚𝑚) = + 4∗0.9 ∗ (0.032 − 0.0222 ) = 2.08
100 𝑠
𝜏 𝜏0 𝜏0 270
For linear relation = − −−→ 𝜏(8𝑚𝑚) = 𝑟 ∗ = 0.022 ∗ 0.03 =
𝑟 𝑅 𝑅
2
198 𝑁/𝑚
∆𝑝
Or 𝜏 = 2𝐿 ∗ 𝑟 𝑓𝑟𝑜𝑚 𝐸𝑞. 7.27
0.022 𝑁
𝜏 = 1800 ∗ 103 ∗ 2∗100 = 198 𝑚2
Table 7.1: Summary of used equations in pipe
𝑅2 𝜕𝑝 𝑟2
Velocity in circular pipe. 𝑉𝑥 = (− ) [1 − 2 ]
4𝜇 𝜕𝑥 𝑅
𝑉𝑚𝑎𝑥 (max. velocity) 𝑉𝑚𝑎𝑥 = 2𝑉𝑎𝑣
𝑅2 𝑑𝑝 1
𝑉𝑎𝑣 (Average velocity) 𝑉𝑎𝑣 = (− ) = 𝑉𝑚𝑎𝑥
8𝜇 𝑑𝑥 2
∆𝑝 𝜇 32𝜇𝑉𝑎𝑣
Pressure loss along pipe = 4𝑉𝑚𝑎𝑥 2 =
𝐿 𝑅 𝑑2
(𝑝1 − 𝑝2 )𝑑
Wall shear stress 𝜏0 =
4𝐿
(𝑝1 − 𝑝2 )𝑟
Shear stress at any r 𝜏=
2𝐿
4𝜏0 𝐿
Energy losses ℎ𝑓 =
𝜌𝑔𝑑
Energy loss by friction 𝐿 𝑉2
factor ℎ𝑓 = 𝑓
𝑑 2𝑔
4 𝐴𝑟𝑒𝑎
Hydraulic diameter 𝑑ℎ =
𝑤𝑒𝑡𝑡𝑒𝑑 𝑝𝑟𝑖𝑚𝑒𝑡𝑒𝑟
2
64 𝐿 𝑉𝑎𝑣 32𝜇𝐿𝑉
Energy loss in Laminar ℎ𝑓 𝑙𝑎𝑚𝑖𝑛𝑎𝑟 = =
flow 𝑅𝑒 𝑑 2𝑔 𝛾𝑑 2
= 128𝜇𝐿𝑄/𝜋𝜌𝑔𝑑 4
Part-2
Turbulent Flow
7.6 Friction Factor Calculations.
Experimentation shows the following to be true in turbulent flow.
1- The head loss varies directly as the length of the pipe.
2- The head loss varies almost as the square of the velocity.
3- The head loss varies almost inversely as the diameter.
4- The head loss depends upon the surface roughness of the interior pipe
wall.
5- The head loss depends upon the fluid properties of density and
viscosity.
6- The head loss is independent of the pressure.
𝐿 𝑉2
ℎ𝑓 = 𝑓 𝑑 . 2𝑔
𝑓 = 𝑓(𝑉, 𝑑, 𝜌, 𝜇, 𝜖, 𝜖́ , 𝑚)
∈ is a measure of the size of the roughness projection and has the dimension
of a length.
∈́ is a measure of the arrangement or spacing of the roughness elements.
m is a form factor.
For smooth ∈=∈ ′ = 𝑚 = 0 → 𝑓 = 𝑓(𝑉, 𝐷, 𝜌, 𝜇) averaged into non-
𝜌𝑑𝑉
dimensionless group namely 𝜇 = 𝑅𝑒
For rough pipes the terms ∈,∈' may be made dimensionless by dividing by d
𝜌𝑑𝑉 ∈ ∈′
∴ 𝑓 = 𝐹 ( 𝜇 , 𝑑 , 𝑑 , 𝑚) Proved by experimental plot of friction factor
aganst the 𝑅𝑒 on a log-log chart. Blasius presented his results by an empirical
formula is valid up to about Re =100000
0.316
𝑓= 1
𝑅𝑒 4
In rough pipe ∈/d is called relative roughness.
∈
𝑓 = 𝐹 (𝑅𝑒, 𝑑) is limited and not permit variation of ∈'/d or m.
Moody has constructed one of the most convenient charts for determining
friction factors. In laminar flow, the straight line masked “laminar flow” and
the Hangen-Poiseuille equation is applied and from which 𝑓 = 64/𝑅𝑒
𝐿 𝑉2 ∆𝑝𝑅 2
ℎ𝑓 = 𝑓 ; 𝑉𝑎𝑣 =
𝑑 2𝑔 8𝜇𝐿
The Colebrook formula provides the shape of ∈/d = constant curves in the
transient region
∈
1 𝑑 2.51
= −0.86 ln ( 3.7 + ) (7.43)
√𝑓 𝑅𝑒 √𝑓
7.7 Simple Pipe Problem.

Six variables enter into the problem for incompressible fluid, which are
Q, L, d, ℎ𝑓 , V, ∈. Three of them are given (L, V, ∈) and three will be find.
Now, the problems type can be solved as follows,
Problem Given To find (unknown)

I 𝑄, 𝐿, 𝑑, 𝑉, ∈ ℎ𝐹
II ℎ𝒇 , 𝐿, 𝑑, 𝑉, ∈ 𝑄
III ℎ𝒇 , 𝑄, 𝐿, 𝑉, ∈ 𝑑
In each of the above problem the following are used to find the unknown
quantity
(i) The Darcy – Weisbach Equation.
(ii) The Continuity Equation.
(iii) The Moody diagram.
In place of the Moody diagram Fig. 7.6, the following explicit formula
for f may be utilized with the restrictions placed on it
1
∈ 106 3
𝑓 = 0.0055 [1 + (2000. 𝑑 + 𝑅𝑒 ) ] Moody equation
∈
4 ∗ 103 ≤ 𝑅𝑒 ≤ 107 & 𝐷 ≤ 0.01
1.325 ∈
𝑓= 2 10−6 ≤ 𝐷 ≤ 10−12 , 5000 ≤ 𝑅𝑒 ≤ 108 (7.44)
∈ 5.74
[ln( + 0.9 )]
3.7𝑑 𝑅𝑒
1% yield diff-with Darcy equation

The following formula can be used without Moody chart is
1 6.9 ∈⁄𝑑 1.11
≈ −1.8 log [𝑅𝑒 + ( 3.7 ) ] (7.45)
𝑓 1/2 𝑑
Eq. 7.45 is given by Haaland which varies less than 2% from Moody chart.
7.7.1 Solution Procedures.

I- Solution for 𝒉𝒇 .
With Q,∈, and d are known
𝑉𝑑 4𝑄
𝑅𝑒 =  = 𝜋𝑑
And f may be looked up in Fig. 7.6 or calculated from Eq. 7.44. Substitution
of (f) in Eq. 7.37 yields ℎ𝑓 the enrgy loss due to flow through the pipe per
unit weight of fluid.
Ex.6
Determine the head (energy) loss for flow of 140 l/s of oil =0.00001
2
𝑚 /𝑠 through 400 m pipe length of 200 mm- diameter cost–iron pipe
Sol.
4𝑄 4(0.14)
𝑅𝒆 = 𝜋𝐷 = 𝜋(0.2)(0.00001) = 𝟖𝟗𝟏𝟐𝟕
The relation roughness is ∈/D= 0.25/200= 0.00125 from a given diagram by
interpolation f= 0.023 by solution of Eq. 7.44, f=0.0234; hence
2
𝐿 𝑉2 400 0.14 1
ℎ𝑓 = 𝑓 𝑑 2𝑔 = 0.023 0.2 [𝜋(0.2)2 ] 2(9.81)
4
ℎ𝑓 = 46.58𝑚.
II- Solution for Discharge Q.

𝑉 & 𝑓 Are unknown then Darcy – Weisbach equation and moody diagram
must be used simultaneously to find their values.
1- Givens ∈/d
f value is assumed by inspection of the Moody diagram
2- Substitution of this trail f into the Darcy – equation produce a trial value
of V.
3- From 𝑽 a trial 𝑹𝒆 is computed.
4- An improved value of f is found from moody diagram with help of 𝑹𝒆
5- When f has been found correct the corresponding 𝑽 and Q is determined
by multiplying by the area.
Ex.7
Water at 15 C flow through a 300mm diameter riveted steel pipe,
∈=3mm with a head loss of 6 m in 300 m. Determine the flow rate in pipe.
Sol.
The relative roughness is ∈/d = 0.003/0.3=0.01, and from diagram a trial
f is taken as (0.038). By substituting into Eq. 7.37 Darcy equation
300 𝑉2
6 = 0.038 0.3 2(9.81)
𝑚
∴ 𝑉 = 1.76 𝑠
𝑚2
At T=15C   = 1.13 ∗ 10−6 𝑠
𝑉𝑑 1.715∗0.3
∴ 𝑅𝑒 = = 1.13∗10−6 = 467278

∈
From the Moody diagram f=0.038 at (𝑅𝑒 & 𝐷)
ℎ𝑓 .𝑑.2.𝑔 6∗0.3∗2∗9.81 𝑚
And from Darcy  𝑉𝑎𝑣 = √ = √ = 1.76 𝑆
𝑓.𝐿 0.038∗300
(6∗0.3)(2)(9.81) 𝑚3
∴ 𝑄 = 𝐴𝑉 = 𝜋 (0.15)2 √ (0.038)(300)
= 0.1245 𝑠
III- Solution for Diameter d.

Three unknown in Darcy-equation f, V , d, two in the continuity
equation V, d and three in the Re number equation
To element the velocity in Eq. 7.37 & in the expression for Re, simplifies the
problem as follows.
𝐿 𝑄2
ℎ𝑓 = 𝑓 𝑑 2
𝑑2 𝜋
2𝑔( )
4
8𝐿𝑄 2
𝑂𝑟 𝑑5 = ℎ 2 𝑓 = 𝐶1 𝑓 (7.46)
𝑓 𝑔𝜋
8𝐿𝑄 2
In which C1 is the known quantity ℎ 2
𝑓 𝑔𝜋
4𝑄
From continuity 𝑉𝑑2 = 𝜋
𝑉𝑑 4𝑄 1 𝐶2
𝑅𝑒 = = = (7.47)
 𝜋 𝑑 𝑑
4𝑄
𝐶2 is the known quantity 𝜋 the solution is now effected by the following
procedure
1- Assume the value of f.
2- Solve Eq. 7.46 for d.
3- Solve Eq. 7.47 for Re.
4- Find the relative roughness ∈/d.
5- With 𝑅𝑒 and ∈/d, Look up new f from a diagram.
6- Use the new f, and repeat the procedure.
7- When the value of f does not change in the two significant steps, all
equations are satisfied and the problem is solved.
Ex.8
Determine the size of clean wrought-iron pipe required to convey 4000
𝑓𝑡 2
gpm oil, =0.0001 𝑠 , 10000 ft pipe length with a head loss of 75 ft .lb/lb.
Sol.
4000
The discharge is 𝑄 = 448.4 = 8.93 𝑐𝑓𝑠
8𝐿𝑄 2 8∗10000∗8.932
From Eq. 7.46, 𝑑5 = ℎ 𝑓= 𝑓 = 𝟐𝟔𝟕. 𝟔𝟓 𝒇
𝑓 𝑔𝜋 2 75∗32.2∗𝜋 2
And from Eq. 7.47,
4𝑄 1 4∗8.93 1 113700
𝑅𝑒 = 𝜋 𝑑 = 𝜋∗0.0001 𝑑 = 𝑑
And from Table 7.2 ∈= 0.00015 ft
If f=0.02 (assumed value) ,  d = 1.35 ft
𝑅𝑒 = 81400
∈/d= 0.00011 from Moody chart f = 0.0191
In repeating the procedure, d = 1.37 ft  𝑅𝑒 = 82991 − −→ 𝑓 = 0.019

Therefore
𝑑 = 1.382 ∗ 12 = 16.6 𝑖𝑛
Figure 7.6: The Moody chart for pipe friction with smooth and rough
walls [1].
Table 7.2: Recommended roughness values [1].

7.7.2 General Applications
Application-1
A pump delivers water from a tank (A)( water surface elevation=110m)
to tank B (water surface elevation= 170m). The suction pipe is 45m long and
35cm in diameter the delivered pipe is 950m long 25cm in diameter. Loss
head due to friction 𝒉𝒇𝟏 = 5𝑚 and hf2 = 3m If the piping are from
pipe(1)= steel sheet metal
pipe(2)= stainless – steel
Calculate the following
i) The discharge in the pipeline
ii) The power delivered by the pump.
Sol.
Given
𝑚2
𝜈𝑤 = 1.007 ∗ 10−6 𝑠
𝑑1 = 35 𝑐𝑚 = 0.35𝑚 ; 𝑑2 = 25𝑐𝑚 = 0.25𝑚
𝐿1 = 45𝑚 ; 𝐿2 = 950 𝑚
From table 7.2 ∈1 = 0.05 𝑚𝑚
∈2 = 0.002 𝑚𝑚
∈1 0.05
= 350 = 1.428 ∗ 10−4
𝑑1
∈2 0.002
= = 8 ∗ 10−6
𝑑2 250
Assume 𝑓1 = 0.013; 𝑓2 = 0.008
𝐿 𝑉2
ℎ𝑓1 = 𝑓1 𝑑1 . 2𝑔
1
1
1 45 𝑉2 𝑚 𝑉𝑑 7.66∗0.35
5 = 0.013 0.35 . 2∗9.81 − −−→ 𝑉1 = 7.66 −→ 𝑅𝑒1 = = 1.007∗10−6
𝑠 𝜈
𝑅𝑒1 = 2662363 = 2.66 ∗ 106
𝐿 𝑉2 950 𝑉2
ℎ𝑓2 = 𝑓2 𝑑2 2𝑔
2 2
= 0.008 0.25 . 2∗9.81 = 3.0𝑚
2
𝑚 1.39∗0.25
𝑉2 = 1.39 𝑅𝑒2 = 1.007∗10−6 = 3.45 ∗ 105
𝑠
1st Trail
∈
(𝑅𝑒1 & 𝑑1 ) − −→ 𝑓1 = 0.0138
1
∈2
(𝑅𝑒2 & 𝑑 ) −→ 𝑓2 = 0.014
2
45 𝑉2 𝑚
1
ℎ𝑓1 = 5 = 0.0138 0.35 . 2∗9.81 −→ 𝑉1 = 7.435 − −−→ 𝑅𝑒1 = 2.58 ∗ 106
𝑠
950 𝑉2 𝑚
2
ℎ𝑓2 = 3 = 0.014 0.25 2∗9.81 − −→ 𝑉2 = 1.051 − −→ 𝑅𝑒2 = 2.6 ∗ 105
𝑠
𝑛𝑑
2 𝑡𝑟𝑖𝑎𝑙
∈
(𝑅𝑒1 & 𝑑1 ) 𝑓1 = 0.0165 , 𝑓2 = 0.015
1
From f1 &f2
1 45 𝑉2
ℎ𝑓1 = 5 = 0.0165 0.35 2∗9.81 → 𝑉1 = 6.8 𝑚/𝑠
2950 𝑉2
ℎ𝑓2 = 3 = 0.015 0.25 2∗9.81 → 𝑉2 = 1.01 𝑚/𝑠
Re1=2.36*106
Re2=2.52*105
3𝑟𝑑 𝑡𝑟𝑖𝑎𝑙
∈ ∈
(𝑅𝑒1 & 𝑑1 ) , (𝑅𝑒2 & 𝑑2 ) → 𝑓1 = 0.0169, 𝑓2 = 0.015
1 2
From Darcy-equation gives V1=0.6.72 m/s, V2=1.016 m/s.
Q=A1*V1= 0.6462 m3/s
From energy equation
𝑝1 𝑉12 𝑝2 𝑉12
+ + 𝑧1 + ℎ𝑝 = + + 𝑧2 + ℎ𝑓
𝛾 2𝑔 𝛾 2𝑔
(6.72)2 (1.016)2
+ 110 + ℎ𝑝 = 2∗9.81 + 170 + 8 Since p1=p2
2∗9.81
hp = 65.75 m
𝑃 = 𝛾𝑄ℎ𝑝 = 9810 ∗ 0.6462 ∗ 65.75 = 416.8 𝑘𝑊 The power delivered by
the pump.
Application-2
In a pipeline of diameter 350mm and length 75m, water is flowing at a
velocity of 2.8 m/s. Find the head lost due to friction, using Darcy-Eq.&
Moody chart, pipe material is Steel–Riveted kinematic viscosity = 0.012
stoke
Sol.
𝐿 𝑉2 𝑚
ℎ𝑓 = 𝑓 𝑑 . 2𝑔 ; 𝑑 = 0.35𝑚 , 𝐿 = 75𝑚; 𝑉 = 2.8 𝑠
From table 7.2 for steel riveted ∈=3.0 mm
∈ 0.003
= 0.35 = 8.57 ∗ 10−3
𝑑
𝑚2 𝑚2
1 = 104 𝑠𝑡𝑜𝑘𝑒 ∴  = 0.012 ∗ 10−4
𝑠 𝑠
𝑉𝑑 2.8∗0.35
𝑅𝑒 = = 0.012∗10−4 = 816666 = 8.1 ∗ 105

∈
𝑎𝑡 (𝑅𝑒 & ) −→ 𝑓 = 0.0358
𝑑
75 2.82
∴ ℎ𝑓 = 0.0358 0.35 2∗9.81 = 3.0𝑚
By determine the value of f by Eq. 7.45.
∈ 1.11
1 6.9 𝑑
1 ≈ −1.8 log [𝑅𝑒 + ( 3.7 ) ]
𝑓2 𝑑
1 6.9 8.57∗10−3 1.11
1 = −1.8 log(8.16∗106 + ( ) ) = 5.2646
3.7
𝑓2
f=0.036 Δf = 0.0002
Application-3
Oil having absolute viscosity 0.1 Pa.s and relative density 0.85 flow
through an iron pipe with diameter 305mm and length 3048 m with flow rate
𝑚3
44.4 ∗ 10−3 . Determine the head loss per unit weight in pipe.
𝑠
Sol.
𝑄 44.4∗10−3 𝑚
𝑉= =1 = 0.61
𝐴 𝜋(0.305)2 𝑆
4
𝑉𝑑𝜌 0.61∗0.305∗850
𝑅𝑒 = = = 1580
𝜇 0.1
i.e the flow is laminar .
64 64
𝑓 = 𝑅 = 1580 = 0.0407
𝑒
𝐿 𝑉 3048 (0.61)2
∴ ℎ𝑓 = 𝑓 𝑑 2𝑔 = 0.0407 ∗ 0.305 ∗ = 7.71 𝑚
2𝑔
7.8 Minor Losses.

The losses which occur in pipelines because of bend, elbows, joints,
valves, etc, are called minor losses hm. the total losses in pipeline are
ℎ𝐿 = ℎ𝑓 + ℎ𝑚 (7.48)
Also, others minor losses can be explain as follows,
A. Losses due to sudden expansion in pipe.

From momentum equation
∑ 𝐹𝑥 = 𝜌𝑄 (𝑉2 − 𝑉1 )
𝑄 = 𝐴2 𝑉2
𝑝1 𝐴1 − 𝑝2 𝐴2 = 𝜌𝐴2 (𝑉22 − 𝑉1 𝑉2 )
Divided by 𝛾𝐴2 since 𝐴1 = 𝐴2
𝑝1 −𝑝2 𝑉22 −𝑉1 𝑉2
= − − − (𝑎)
𝛾 𝑔
From Bermonlli's equation between section 1&2 as in Fig.7.7
𝑝1 1𝑉2 𝑝2 2 𝑉2
+ 2𝑔 = + 2𝑔 + ℎ𝑚
𝛾 𝛾
𝑝1 −𝑝2 𝑉22 −𝑉12
= + ℎ𝑚 − − − (𝑏)
𝛾 2𝑔
Equating Eq. (a & b)
𝑉22 −𝑉1 𝑉2 𝑉22 −𝑉12
= + ℎ𝑚
𝑔 2𝑔
𝑉22 −𝑉1 𝑉2 𝑉22 −𝑉12 2𝑉22 −2𝑉1 𝑉2 −𝑉22 +𝑉12
ℎ𝑚 = − =
𝑔 2𝑔 2𝑔
𝑉12 −2𝑉1 𝑉2 +𝑉22 (𝑉1 −𝑉2 )2

ℎ𝑚 = = − − − (𝑐)
2𝑔 2𝑔
𝑄 = 𝐴1 𝑉1 = 𝐴2 𝑉2
𝐴
𝑉2 = 𝐴1 𝑉1 𝑠𝑢𝑏𝑠𝑡𝑖𝑡𝑢𝑡𝑒 𝑖𝑛 𝐸𝑞. (𝑐)
2
1 2
V1 V2
V1 V2
p1 A1 p2 A2
Figure 7.7: Sudden expansion.
𝐴 𝐴 2
𝑉12 −2𝑉1 1 𝑉1 +( 1 𝑉1 )
𝐴2 𝐴2
∴ ℎ𝑚 = 2𝑔
𝐴 𝐴 2
𝑉12 (1−2 1 + ( 1 ) )
𝐴2 𝐴2
ℎ𝑚 = 2𝑔
𝑉12 1 𝐴1 2 𝑉2
ℎ𝑚 = 2𝑔 (1 − 𝐴 ) = 𝐾 2𝑔 (7.49)
2
2 2
𝐴 𝑑2
𝐾 = (1 − 𝐴1 ) = (1 − 𝑑12 )
2 2
𝑑
If sudden expansion from pipe to a reservoir 𝑑1 = 0
2
𝑉12
∴ ℎ𝑚 = 2𝑔
B. Head loss due to a sudden contraction in the pipe cross section.
𝑄 = 𝐴2 𝑉2 = 𝐴𝑐 𝑉𝑐
c 2
𝐶𝑐 = 𝐴𝑐 /𝐴2
𝑝𝑐 𝐴2 − 𝑝2 𝐴2 = 𝜌𝑄(𝑉2 − 𝑉𝑐 )
𝑝𝑐 𝐴2 − 𝑝2 𝐴2 = 𝜌𝐴2 (𝑉22 − 𝑉2 𝑉𝑐 )
Divided by 𝛾 𝐴2
𝑝𝑐 −𝑝2 𝑉22 −𝑉𝑐 𝑉2
= − − − (𝑎)
𝛾 𝑔
Applying B.E between sections c & 2 c 2
𝑝𝒄 𝑉𝒄𝟐 p2 V22
+ 𝟐𝒈 = + 2g + hm Figure 7.8: Sudden contraction.
𝜸 γ
pc −p2 V22 −V2c

= + hm − − − −(b)
γ 2g
A
Vc = A 2 𝑉2 , Equating Eq (a& b)
C
𝑉22 −𝑉𝑐 𝑉2 𝑉22 −𝑉𝐶2
= + ℎ𝑚
𝑔 2𝑔
𝑉22 −𝑉𝐶 𝑉2 𝑉22 −𝑉𝑐2 2𝑉22 −2𝑉𝑐 𝑉2 −𝑉22 +𝑉𝑐2
ℎ𝑚 = −( )=
𝑔 2𝑔 2𝑔
𝑉𝑐2 −2𝑉𝑐 𝑉2 +𝑉22 𝑉22 𝐴2 2 2𝐴2
ℎ𝑚 = = 2𝑔 ((𝐴 ) − + 1)
2𝑔 𝑐 𝐴𝑐
2
2 𝑉2 1 2
ℎ𝑚 = 2𝑔 [ (𝐶 ) − 1 ]
𝑐
𝑉22
ℎ𝑚 = 𝐾 2𝑔 (7.50)
1 2 𝑑2
𝐾 = (𝐶 − 1) , From experimental, 𝐾 ≈ 0.42 (1 − 𝑑22 )
𝑐 1
𝐴𝑐 is the cross – sectional area of the vena-contracts
𝐴
𝐶𝑐 is the coefficient of contraction is defined by 𝐶𝑐 = 𝐴𝑐
2
C. The head loss at the entrance to a pipe line.
From a reservoir the head losses is usually taken as
0.5𝑉 2
ℎ𝑚 = is the opening is square-edged
2𝑔
0.01𝑉 2 0.05𝑉 2
ℎ𝑚 = ~ if the rounded entrance
2𝑔 2𝑔
Table 7.3 lists the loss coefficient K for four types of valve, three angles of
elbow fitting and two tee connections. Fitting may be connected by either
internal screws or flanges, hence the two are listings.
Table 7.3: Head loss coefficients K for typical fittings.
Entrance losses are highly dependent upon entrance geometry, but exit
losses are not. Sharp edges or protrusions in the entrance cause large zones of
flow separation and large losses as shown in Fig. 7.9. As in Fig. 7.10, a bend
or curve in a pipe, always induces a loss larger than the simple Moody
friction loss, due to flow separation at the walls and a swirling secondary
flow arising from centripetal acceleration.
Figure 7.9: Entrance and exit loss coefficients,(a) reentrant inlets,

(b) rounded and beveled inlets. Exit losses are K=1.
Figure 7.10: Resistance coefficients for 90 bends.

Table 7.3 gives the losses coefficients for the fully open condition. In
case of partially open valve the losses can be much higher. Fig. 7.11 gives
average losses for three valves as a function of percentage open. The opining
distance ratio h/D as the x-axis in Fig. 7.11 is shown by Fig.7.12 of valve
geometry.
Figure 7.11: Average-loss coefficients for partially open valves.
Figure 7.12: Typical commercial valve; (a) gate valve, (b) globe valve,
(c) angle valve, (d) swing-check valve, (e) disk-type gate valve.
Application-4
𝑓𝑡 2
Water,  = 1.94 slugs/ft3, 𝑎𝑛𝑑 𝜈 = 1.1 ∗ 10−5 , is pumped between
𝑠
𝑓𝑡 3
two reservoir at 0.2 through 400 ft of 2 in diameter pipe and several minor
𝑠
∈
losses, as shown in figure. The roughness ratio is 𝑑 = 0.001. Compute the
horse power required.
Sol.
Write the steady- flow energy equation between section 1 &2 the two
reservoir surface:
𝑝1 𝑉12 𝑝2 𝑉22
+ + 𝑧1 = ( + + 𝑧2 ) + ℎ𝑓 + ∑ ℎ𝑚 − ℎ𝑝
𝜌𝑔 2𝑔 2𝑔 2𝑔
Where ℎ𝑝 is the head increase across the pump 𝑝1 = 𝑝2 , 𝑉1 = 𝑉2 ≈ 0,
solve for the pump head
𝑉2 𝐿
ℎ𝑝 = 𝑧2 − 𝑧1 + ℎ𝑓 + ∑ ℎ𝑚 = 120 − 20 + 2𝑔 ( 𝑓 𝑑 + ∑ 𝐾)
𝑄 0.2 𝑓𝑡
𝑉=𝐴=1 2 2
= 9.17
𝜋( ) 𝑠
4 12
Now list and sum the minor loss coefficients
Loss K
Sharp entrance (fig, 7.9) 0.5
Open globe value (2 in Table 7.3) 6.9
𝑅 ∈ 0.15
12-in bend Fig. 7.10 𝑑 = 6, 𝑑 = 0.001
Regular 90 elbow ( 𝑇𝑎𝑏𝑙𝑒 7.3) 0.95
Half – closed gate value ( Fig. 7.11) 3.8
Sharp exit (Fig. 7.9) 1.0
∑𝐾 13.3
Calculate the 𝑅𝑒 and pipe friction factor

2
𝑉𝑑 9.17( )
𝑅𝑒 = = 1.1∗1012−5 = 139000
𝜈
∈
For 𝑑 = 0.001 , from the Moody chart read f = 0.0216
9.172 0.0216(400)
∴ ℎ𝑝 = 100 + 2(32.2) [ 2 + 13.3]
( )
12
ℎ𝑝 = 100 + 84 = 184 𝑓𝑡 𝑝𝑢𝑚𝑝 ℎ𝑒𝑎𝑑
𝑓𝑡 𝑙𝑏𝑓
𝑃 = 𝜌𝑔𝑄ℎ𝑝 = 1.94(32.3)(0.2)(184) = 2300 𝑠
𝑓𝑡 𝑙𝑏𝑓 2300
1 ℎ. 𝑝. = 550 𝑠 ∴ 𝑃 = 550 = 4.2 ℎ. 𝑝.
For an efficiency 70 to 80% , a pump is needed with an input power about 6
h.p.
Application-5
Water is to be supplied to the inhabitants of a college campus through a
supply main. The following data is given,
Distance of the reservoir from the campus = 3000 m
Number of inhabitance = 4000
Consumption of water per day of each inhabitant= 180 liters
Loss of head due to friction = 18 m,
Co-efficient of friction for the pipe, f = 0.007. If the half of the daily supply is
pumped in 8 hr, determine the size of the supply main (d) .
Sol.
180 𝑚3
Total supply per day = 4000 ∗ 1000 = 720 𝑑𝑎𝑦
720 𝑚3
𝑄 = 2∗8∗3600 = 0.0125 𝑠
ℎ𝑓 = 18𝑚
8∗𝐿 𝑄 2 8∗3000∗(0.0125)2
Assume f= 0.03 ∴ 𝑑5 = ℎ ∗𝑓 = ∗ 0.03 = 6.455 ∗
𝑓 𝑔 𝜋2 18∗9.81∗𝜋 2
10−5
∴ 𝑑 = 0.1452𝑚
4𝑄 4∗0.0125
𝑅𝑒 = 𝜋𝜈𝑑 = 𝜋∗1.007∗10−6 ∗0.1452 = 108904 = 1.08 ∗ 105
∈ 0.046∗10−3 ∈
= 0.1452 = 3.168 ∗ 10−4 From ( 𝑅𝑒 & ) − −→ 𝑓 = 0.0195
𝑑 𝑑
5 −5
𝑑 = 4.2 ∗ 10 − −→ 𝑑 = 0.13352 𝑚 .
Minor losses may be expressed in terms of the equivalent Le of pipe that

has the same head loss in m.N /N
𝐿𝑒 𝑉 2 𝑉2
𝑓 = 𝐾 2𝑔
𝑑 2𝑔
𝐾𝑑
K is the sum of several losses, solving for 𝐿𝑒 gives 𝐿𝑒 = (7.51)
𝑓
Application-6
I) Find the discharge through the pipeline as in below figure for H=10m,
II) determine the head loss hL for Q=60 l/s. III) compare the result of
discharge with equivalent length .
Sol.
The energy equation applied between points 1 & 2, including all the
losses, may be written
2 𝑉2 2 1 𝑉2
2 102 𝑉 2
2 𝑉2 5.3 𝑉22
I) H1 + 0 + 0 = 2𝑔 + 0 + 0 + 2 2𝑔 + 𝑓 0.2032 2𝑔 + 2 ∗ 0.26 2𝑔 + 2𝑔
Loss coefficients (K):-
Entrance = 0.5
Each elbow = 0.26
Globe valve (partially open h/d=0.6) = 5.3
2 𝑉2
∴ 𝐻1 = 2𝑔 (7.32 + 502 𝑓) − − − − − (𝐴)
When the head is given, this problem is solved as the second type of simple
∈ 0.26
pipe problem. If 𝑑 = 203.2 = 1.28 ∗ 10−3 , f = 0.0205
2 𝑉2 𝑚
10 = 2𝑔 ( 7.32 + 502 ∗ 0.0205) − −→ 𝑉2 = 3.337
𝑠
𝑚2
𝜈 = 1.01 ∗ 10−6 𝑠
∈ (3.337∗0.2032)
= 0.00128 ; 𝑅𝑒 = = 6.7 ∗ 105
𝑑 (1.01∗10−6 )
∈
From Moody chart at { 𝑅𝑒 & 𝑑} − −→ 𝑓 = 0.0208
𝑚 ∈
Repeating the procedure gives 𝑉2 = 3.32 𝑠 , 𝑅𝑒 = 6.6 ∗ 105 , 𝑎𝑛𝑑 𝑑 , 𝑓 =
0.0209. from eq. A gives V2=3.31 m/s. The discharge is
𝜋
𝑄 = 𝑉2 𝐴2 = (3.31) ( 4 ) (0.2032)2 = 107.34 𝑙/𝑠
II) For the second part , with Q is known, the solution is straight forward;
𝑄 0.06 𝑚 ∈
𝑉2 = 𝐴 = 𝜋 2
= 1.85 𝑠 𝑅𝑒 = 3.7 ∗ 105 𝑎𝑛𝑑 𝑑 , 𝑓 = 0.0212
( )(0.2032)
4
From Eq. A
(1.85)2
∴ ℎ𝐿 = 2(9.806) (6.32 + 502 ∗ 0.0212) = 2.959 𝑚
III) With equivalent lengths Eq. 7.51 the value of f is an approximated, say
f = 0.0205. The sum of minor losses is K= 6.32
𝐾𝑑 6.32∗0.2032
𝐿𝑒 = 𝑓 = 0.0205 = 62.64 𝑚
The total length of pipe is
62.64 + 102 = 164.64 𝑚
𝐿+𝐿𝑒 𝑉22 164.64 𝑉 2
2
By Darcy equation. 10 = 𝑓 𝑑 2𝑔
= 𝑓 0.2032 2𝑔
𝑚 ∈
If 𝑓 = 0.0205, 𝑉2 = 3.43 , 𝑅𝑒 = 6.9 ∗ 105 𝑎𝑛𝑑 𝑑 , 𝑓 = 0.0203
𝑠
𝑚
From eq.A , 𝑉2 = 3.347 𝑎𝑛𝑑 𝑄 = 108.5 𝑙/𝑠
𝑠
H=10
203.2 mm –diam. Clean cast-iron
2
m
Water at
Globe valve
20 C (Partially open)
12 m
Standard elbows flanged type

90 regular
30 60 m
Square- m
edge
entrance
7.9 Pipe in Series.

In this typical series–pipe system as in Fig. 7.13, the H (head) is
required for a given Q or the Q wanted for a given H. Applying the energy
equation from A to B including all losses gives
𝑉2 𝐿 𝑉2 (𝑉1 −𝑉2 )2 𝐿 𝑉2 𝑉2
1
𝐻 + 0 + 0 = 0 + 0 + 0 + 𝐾𝑒 2𝑔 + 𝑓1 𝑑1 2𝑔
1
+ + 𝑓2 𝑑2 2𝑔
2 2
+ 2𝑔 (7.52)
1 2𝑔 2
From continuity eqn.
𝑉1 𝑑12 = 𝑉2 𝑑22
𝑉2 is eliminated from eqn , so that
𝑉2 𝑓1 𝐿1 𝑑 2 2 𝑓2 𝐿2 𝑑 4 𝑑 4
1
𝐻 = 2𝑔 { 𝐾𝑒 + + [ 1 − (𝑑1 ) ] + (𝑑1 ) + (𝑑1 ) } (7.53)
𝑑1 2 𝑑2 2 2
For known lengths and sizes of pipes this reduces to

𝑉2
1
𝐻 = 2𝑔 ( 𝐶1 + 𝐶2 𝑓1 + 𝐶3 𝑓2 ) (7.54)
𝐶1 , 𝐶2 & 𝐶3 𝑎𝑟𝑒 𝐾𝑛𝑜𝑤𝑛
 Q is given then 𝑅𝑒 is computed and f's may be looked up in Moody
diagram then H is found by direct substitution.
 For a given H, the values of V1,f1,f2 are unknown in Eq. 7.54.
By Assuming values of 𝑓1 & 𝑓2 (may be equaled) then 𝑉1 𝑖𝑠 𝑓𝑜𝑢𝑛𝑑 1𝑠𝑡 𝑡𝑟𝑖𝑎𝑙
and from 𝑉1 − −→ Re's are determined and values of 𝑓1 & 𝑓2 look up from
Moody diagram.
And at these value , a better 𝑉1 is computed from Eq. 7.53 since f varies so
slightly with 𝑅𝑒 the trial solution converges very rapidly. The same
procedures apply for more than two pieces in series.
A
H
B
Water
2 Water
K
e
Figure 7.13: Series pipe.
system
7.10 Equivalent pipes.
Two pipe system (in series) are said to be equivalent when the same
head loss produces the same discharge in both system. From Darcy equation.
𝐿 𝑄12 𝐿 8𝑄12
ℎ𝑓1 = 𝑓1 𝑑1 2 = 𝑓1 𝑑15
1 (𝑑2 𝜋) 2𝑔 1 𝜋2 𝑔
14
And for a second pipe
𝐿 8𝑄22
ℎ𝑓2 = 𝑓2 𝑑25
2 𝜋2𝑔
For the two pipes to be equivalent
ℎ𝑓1 = ℎ𝑓2 𝑄1 = 𝑄2
After equating ℎ𝑓1 = ℎ𝑓2 and simplifying
𝑓1 𝐿1 𝑓2 𝐿2
=
𝑑15 𝑑25
Solving for 𝐿2 gives
𝑓 𝑑 5
𝐿2 = 𝐿1 𝑓1 (𝑑2 ) (7.55)
2 1
Ex.9
From Fig. 7.13, 𝐾𝑒 = 0.5, 𝐿1 = 300𝑚 , 𝑑1 = 600𝑚𝑚, ∈1 = 2𝑚𝑚, 𝐿2 =
𝑚2
240𝑚, 𝑑2 = 1𝑚, ∈2 = 0.3𝑚𝑚 , 𝜈 = 3 ∗ 10−6 𝑠
𝑎𝑛𝑑 𝐻 = 6𝑚. Determine
the discharge through the system.
Sol.
From the energy Eq. 7.53
𝑉2 300 240
1
6 = 2𝑔 [ 0.5 + 𝑓1 + (1 − 0.62 )2 + 𝑓2 0.64 + 0.64 ]
0.6 1.0
After simplifying
𝑉2
1
6 = 2𝑔 (1.0392 + 500𝑓1 + 31.104𝑓2 )
∈1 ∈
From = 0.0033, 𝑑2 = 0.0003, and Moody diagram values of f 's are
𝑑1 2
assumed for the fully turbulent range.
𝑓1 = 0.026 𝑓2 = 0.015
𝑚 𝑚
By solving for 𝑉1 with these value , 𝑉1 = 2.848 𝑠 , 𝑉2 = 1.025 𝑆
2.848∗0.6
𝑅𝑒1 = = 569600
3∗10−6
1.025∗1.0
𝑅𝑒2 = 3∗10−6 = 341667
At these Re's and from Moody diagram, 𝑓1 = 0.0265, 𝑓2 = 0.0168, by
𝑚 𝑚3
solving again for 𝑉1 , 𝑉1 = 2.819 and 𝑄 = 0.797
𝑠 𝑆
Ex.10
Solve Ex.9 by means of equivalent pipes.
Sol.
First by expressing the minor losses in terms of equivalent length, for
pipe 1, since K=Ke+(1-(d1/d2)2)2
𝐾 𝑑 0.91∗0.6
𝐾1 = 0.5 + (1 − 0.62 )2 = 0.91 − −−→ 𝐿𝑒1 = 1𝑓 1 = 0.026 = 21𝑚
1
𝐾2 𝑑2 1∗1
For pipe-2 𝐾2 = 1 − −−→ 𝐿𝑒2 = = 0.015 = 66.7 𝑚
𝑓2
After selecting𝑓1 & 𝑓2 . The problem is reduced to 321m of 600-mm diam. &
306.7 m of 1-m pipe diam.
By expressing the 1-m pipe terms of an equivalent length of 600-mm pipe, by
Eq. 7.55
𝑓 𝑑 5 0.015 0.6 5
𝐿𝑒 = 𝑓2 𝐿2 (𝑑1 ) = 306.7 0.026 ( 1 ) = 13.76 𝑚
1 2
Now, by adding to the 600-mm pipe, the problem is reduced to 334.76m of
600-mm for finding the discharge through it, ∈1 = 2𝑚𝑚, 𝐻 = 6 𝑚.
334.76 𝑉 2
6=𝑓 0.6 2𝑔
0.6
𝑓 = 0.026, −−→ 𝑉 = 2.848 −→ 𝑅𝑒 = .848 ∗ = 569600
(2 ∗ 10−6 )
∈
From Re and = 0.0033 from Moody diagram f=0.0265
𝑑
from above equation
𝑚 𝑚3
𝑉 = 2.821 → 𝑄 = 𝜋 (0.3)2 (2.821) = 0.798
𝑠 𝑠
7.11 Pipes in Parallel.

The second type of pipe- system is the parallel flow type, in this case as
shown in Fig. 7.14 the head losses are same in any of the lines and the total
flow is the sum of flow rate in each pipe.
The minor losses are added into the lengths of each pipe as equivalent
lengths. From Fig. 7.14 the conditions to be satisfied are
𝑝 𝑝
ℎ𝑓1 = ℎ𝑓2 = ℎ𝑓3 = 𝛾𝐴 + 𝑧𝐴 − ( 𝛾𝐵 + 𝑧𝐵 ) (7.56)
𝑄 = 𝑄1 + 𝑄2 + 𝑄3
Figure 7.14: Parallel pipes system.
𝑧𝐴 & 𝑧𝐵 are the elevations of point A&B.

Q is the discharge through the approach pipes.
Two types of problems occur
1) The elevations of HGL(hydraulic grid line) at A&B are known, to
find the discharge.
2) Q is known, to find the distribution of flow and the head loss, size of
pipe, fluid properties and roughness's are assumed to be known.
Case-1. as the simple pipe problem. Since, head loss is the drop in HGL.
These discharges are added to determine the total discharge.
Case-2. The recommended procedure is as follows.
1- Assume a discharge Q', through pipe 1.
2- Solve forℎ′𝑓1 , using assumed discharge.
3- Using ℎ′𝑓1 , find 𝑄′2 , 𝑄′3 .
4- With the three discharges for a common head loss, now assume that
the given Q is split up among the pipes in the same proportion as
𝑄′1 , 𝑄′2 , 𝑄′3 ; thus
Q′ Q′ Q′
Q1 = ∑ Q1′ 𝑄; 𝑄2 = ∑ Q′2 𝑄; 𝑄3 = ∑ Q′3 𝑄 (7.57)
5- Check the correctness of these discharges by computing ℎ𝑓1 , ℎ𝑓2 , ℎ𝑓3

for the computing 𝑄1 , 𝑄2 , 𝑄3 .
Ex.11
In Fig. 7.14, 𝐿1 = 3000 𝑓𝑡 , 𝑑1 = 1 𝑓𝑡 , ∈1 = 0.001 𝑓𝑡
𝐿2 = 2000 𝑓𝑡 , 𝑑2 = 8𝑖𝑛 , ∈2 = 0.0001 𝑓𝑡
𝐿3 = 4000 𝑓𝑡 , 𝑑3 = 16 𝑖𝑛 , ∈3 = 0.0008 𝑓𝑡
𝑆𝑙𝑢𝑔𝑠 𝒇𝒕𝟐
𝜌 = 2.00 , 𝝂 = 0.00003
𝑓𝑡 3 𝒔
𝑝𝐴 = 80 𝑝𝑠𝑖, 𝑧𝐴 = 100 𝑓𝑡 , 𝑧𝐵 = 80 𝑓𝑡 .
For a total flow of 12 cfs, determine the flow through each pipe and the
pressure at B.
Sol.
For pipe-1-
𝑓𝑡
Assume 𝑄1′ = 3 𝑐𝑓𝑠 ; 𝑡ℎ𝑒𝑛 𝑉1′ = 3.82 𝑠
3.82∗1
∴ 𝑅𝑒1′ = 0.00003 = 127000 .
∈1
= 0.001 From Moody chart 𝑓1′ = 0.022
𝑑1
′ 3000 3.822
∴ ℎ𝑓1 = 0.022 = 14.97 𝑓𝑡
1.0 64.4
For pipe-2-
2000 𝑉 ′2
14.97 = 𝑓2′ 0.667 2𝑔
2
− − − − − (𝑎)
Assume 𝑓2′ = 0.020 (Recomended fully turbulent flow)
2
𝑓𝑡 4.01∗ ∗1
𝑡ℎ𝑒𝑛 𝑉2′ = 4.01 − −→ 𝑅𝑒2′ = 3
= 89000
𝑠 0.00003
∈2
= 0.00015
𝑑2
From Moody chart  𝑓2′ = 0.019
𝑓𝑡
Then from Eq. (a) 𝑉2′ = 4.11 𝑠 − −→ 𝑄2′ = 1.44 𝑐𝑓𝑠
For pipe -3-
4000 𝑉3′2
14.97 = 𝑓3′ 1.333 − − − − − (𝑏)
2𝑔
𝑓𝑡 4.01∗1.333
Assume 𝑓3′ = 0.019 𝑡ℎ𝑒𝑛 𝑉3′ = 4.01 − −→ 𝑅𝑒3′ = = 178000
𝑠 0.00003
∈3
= 0.0006
𝑑3
∈
From Moody chart. At (𝑅𝑒3′ & 𝑑3 ) − −→ 𝑓3′ = 0.02 from Eq.(b).
3
𝑓𝑡
𝑉3′ = 4.0 𝑠 & 𝑄3′ = 5.6 𝑐𝑓𝑠
The total discharge for the assumed conditions is
∑ 𝑄 ′ = 3.0 + 1.44 + 5.6 = 10.04 𝑐𝑓𝑠
From Eq. 7.57
𝑄′ 3.00
𝑄1 = ∑ 𝑄1 ′ . 𝑄 = 10.04 12 = 3.58 𝑐𝑓𝑠
1.44
𝑄2 = 10.04 12 = 1.72 𝑐𝑓𝑠
5.6
𝑄3 = 10.04 12 = 6.7 𝑐𝑓𝑠
Check the values of ℎ𝑓1 , ℎ𝑓2 , ℎ𝑓3
3.58 𝑓𝑡 ∈1
𝑉1 = = 4.56 − −→ 𝑅𝑒1 , = 152000}−→ 𝑓1 = 0.021
1 𝑠 𝑑1
𝜋∗4
hf1=20.4ft
1.72 𝑓𝑡 ∈2
𝑉2 = 𝜋 = 4.93 − −→ 𝑅𝑒2 = 109200}−→ 𝑓2 = 0.019
𝑠 𝑑2
9
hf2=21.6 ft
6.7 𝑓𝑡 ∈
𝑉3 = 4𝜋 = 4.8 −→ 𝑑 3 𝑅𝑒3 = 213000}−→ 𝑓3 = 0.019𝑓2
𝑠 3
9
hf3=20.4ft
is about midway between 0.018 & 0.019. to satisfy the condition ℎ𝑓1 =
ℎ𝑓2 = ℎ𝑓3
∴ 𝑖𝑓 𝑓2 = 0.018 𝑡ℎ𝑒𝑛 ℎ𝑓2 = 20.4 𝑓𝑡 𝑠𝑎𝑡𝑖𝑠𝑓𝑦𝑖𝑛𝑔 .
To find pB
𝑝𝐴 𝑝
+ 𝑧𝐴 = 𝛾𝐵 + 𝑧𝐵 + ℎ𝑓
𝛾
𝑝𝐵 144
𝑜𝑟 = 80 ∗ 64.4 + 100 − 80 − 20.4 = 178.1𝑓𝑡
𝛾
In which the average head loss was taken. Then
178.1∗64.4
𝑝𝐵 = 144 = 79.6 𝑝𝑠𝑖
Problems.
P7.1 Two parallel plates kept 100mm a part have laminar flow of oil between
them with a maximum velocity of 1.5 m/s Calculate
i) The discharge per meter width.
ii) The shear stress at the plates.
iii) The difference in pressure between two points 20m a part.
iv) The velocity gradient of the plates.
v) The velocity at 20mm from the plate.
Assume viscosity of oil 𝜇 = 24.5 𝑝𝑜𝑖𝑠𝑒.
P7.2 A liquid of viscosity oil 0.1N.s/m2 is filled between two parallel plates
10mm a part. If the upper plate is moving at 2m/s and the pressure
difference between two sections 10m apart is 9.81kN/m2. Find the
relation of velocity and also determine the shear stress on the moving
plate.
P7.3 What is the pipe diameter which must be used for oil flow rate
0.0222m3/s at 15.6 C the energy loss per unit weight ℎ𝑓 = 22.0𝑚 for a
horizontal pipe with length 1000m. Using =0.00021m2/s & relative
density= 0.912kg/m3.
P7.4
a) Calculate the shear stress at the wall of pipe if pipe diameter is
305mm and water head loss is 15m per weight of 300m length.
b) Calculate the shear stress at a point 51mm from centerline axis of
pipe.
c) Determine the average velocity when f=0.05.
P7.5 The dynamic viscosity of oil is 0.1Pa.s and relative density is 0.85. The
oil flow rate is (44.4 ∗ 10−3 )m3/s through a pipe with 305mm diameter
and 3048m length. Determine the energy loss per unit weight.
P7.6 A heavy oil flow from A to B through a horizontal pipe of cast–iron

with 153mm diameter and 104.4m length. The pressure at section A is
1.069MPa and at B is 34.48kPa. The kinematic viscosity of oil
412.5*10-6m2/s and relative density 0.918. Determine the flow rate.
P7.7 The distance between two point A&B along pipe line is 1224m, pipe
diameter is 153mm. Point B level is above point A about 15.39m . The
pressures at A&B are 848kPa and 335kPa respectively. If the pipe
manufacturing from wrought iron determines the discharge of oil flow
between A&B, take =3.83*10-6m2/s, =854kg/m3.
P7.8 Calculate the energy loss per unit weight in cast-iron pipe with 278mm
diameter and 450m length when the flow are
a) Water at 15.6C and velocity is 1750mm/s.
b) Oil at 15.6C and same velocity, take, water=1.13*10-6m2/s,
oil=4.41*10-6m2/s.
P7.9 Determine the diameter of wrought iron pipe when the water flow rate is
1.25m3/s and pipe length is 3225m. The drop in hydraulic grad line is
68.8m.
P7.10 The parallel galvanized–iron pipe system as in figure delivers water at

20C with total flow rate of 0.036m3/s. If the pump is wide open and
not running, with a loss coefficient k= 1.5, determine
a) The flow rate in each pipe.
b) The overall pressure drop.
P7.11 Show that the discharge per unit width between two parallel plate
distance (t) a part. When one plate is moving at velocity U while the
other one is held stationary for the condition of zero shear stress at the
fixed plate is q=Ut/3.
P7.12 An oil of viscosity 9poise and S.G.= 0.9 it is flowing through a

horizontal pipe of 60mm diameter. If the pressure drop in 100m
length of the pipe is 1800kN/m2 determine.
i. The rate of flow oil.
ii. The centerline velocity.
iii. Total frictional drag over 100m length.
iv. The power required to maintain the flow.
v. The velocity gradient at the pipe wall.
vi. The velocity and shear stress at 8mm from the wall.
P7.13 Oil flow in the pipe as shown in the figure has S.G.= 1.26 and
viscosity 1.5N.s/m2 and consider the head of pump is 40m. Find the
flow rate in the pipe, Hint (the head loss after the pump is enter in
calculation).
Pump
L=120m 8m
d=0.15
m
P7.14 Find the head loss of flow through a pipe the flow rate in the pipe is
200L/s, and the pipe made from stainless steel. The fluid has S.G. = 0.9
and dynamic viscosity 0.07 N.s/m2. The diameter of pipe is 300mm and
length of pipe is 200m.
P7.15 Oil flow through a pipe has 300mm diameter commercial steel pipe. If
the head loss for 400m pipe length is about 10m. Determine the flow
rate in the pipe. Consider S.G.=0.9, =0.09 N.s/m2.
P7.16 The system in the following figure consists a tank volume 20m3 and is
filled after 1000s. If the friction in pipes is 0.03 and ∈=
0.000075 𝑚 find,
i) The pressure developed by the pump on its delivery side.
ii) The power delivered to water by the pump.
d=150mm
R
R=0.45m 20m
Pump R
D=300mm
50m 50m
Flow
P7.17 Three pipes of diameters 300mm, 200mm and 400mm and lengths
450m, 255m and 315m respectively are connected in series between
two tanks. The difference in water surface levels in two tanks is 18m.
Determine the rate of flow of water if the coefficients of friction for
three pipes are 0.03, 0.0312 and 0.0288 respectively considering
i) Minor losses.
ii) Neglecting minor losses.
P7.18 Solve the above problem by the equivalent pipe technique consider the
diameter of equivalent pipe is 0.2m and f=0.0312. Assume with
minor losses.
P7.19 Two pipes have diameters 200mm, 250mm and lengths 200m and
150m respectively connected in parallel as shown in figure. If the
total flow rate is 0.1m3/s and the points (1) and (2) have same
elevation and same diameters calculate
i) the flow rate in each pipe.
ii) pressure drop between (1) of (2) assume the fluid is oil has
S.G.=0.9 and dynamic viscosity = 0.002N.s/m2 and ∈ for two
pipe is 0.25mm
Pipe-1 , d1, L1
1 2
Pipe-2 , d2, L2
P7.20 the following figure consists of 1250m of 5.8cm cast-iron pipe, two
45 regular screwed type and four 90screwed long-radius elbows, a
fully open screwed globe valve and a sharp exit into reservoir. If the
elevation at point 1 is 425m, what gage pressure is required at point 1
to deliver 0.0045m3/s of water at 25C into reservoir?
Chapter 8 Introduction to Boundary Layer 193
CHAPTER 8
Introduction to Boundary
Layer
8.1 Boundary Layer Definitions and Characteristics.
Boundary Layer is a small region developing around a boundary surface
of body, in which the velocity of the flowing fluid increase rapidly from zero
at the boundary surface and approaches the velocity of main stream. The
layer adjacent to the boundary surface is known as boundary Layer (B.L.).
Firstly Introduced by L.Prandtl in 1904. Fluid medium around bodies moving
in fluids, can be divided into following two regions
(i) A thin layer adjoining the boundary called B.L where the viscous
shear takes place.
(ii) A region outside the boundary layer where the flow behavior is
quite like that of an ideal fluid and the potential flow theory is
applicable.
𝑈∞ : is the velocity at the outer edge of the B.L.
𝛿: is called the dynamic B.L thickness where u=0.99U as shown in Fig. 8.1.
𝑇𝑤 : is the wall temperature where the fluid immediately at the surface is
equal to the temperature of the surface.
𝛿𝑇 : is the thermal B.L thickness, where the temperatures are changing as
𝑇 = 𝑇𝑤 𝑎𝑡 𝑦 = 0 𝑇 = 𝑇∞ 𝑎𝑡 𝑦 = 𝛿𝑇
In dynamic B.L. u=u(y) u=0 at y=0; U=𝑈∞ 𝑎𝑡 𝑦 = 𝛿,
𝛿 = 𝛿(𝑥) 𝑎𝑛𝑑 𝛿𝑇 = 𝛿𝑇 (𝑥)
𝜕𝑢
𝜏0 = 𝜇( )𝑤 is the shear stress at the wall. (8.1)
𝜕𝑦
The displacement of the streamlines ( 𝛿𝑑 ) in the free stream as a result
of velocity deficits in the B.L is known the displacement thickness. The
momentum layer thickness ( 𝜃 ) is the equivalent thickness of a fluid layer
with velocity U with momentum equal to the momentum lost due to friction
, and is defind as the momentum thickness 𝜃.
𝜕𝑇
𝑞𝑤 = −𝑘 ( ) is the heat flux at the wall. (8.2)
𝜕𝑦 𝑤
Both ( 𝜏𝑤 , 𝑞𝑤 ) are function of distance from the leading edge

𝜏𝑤 = 𝜏𝑤 (𝑥), 𝑞𝑤 = 𝑞𝑤 (𝑥)
Figure 8.1: Growth of a boundary layer on a flat plate.
8.2 Boundary Layer Theory (Flow over Flat Plate).

For flow over a flat plate with zero pressure gradient, the transition
process occurs when
𝑈∞ 𝑥𝐿 𝑚
𝑅𝑒 = = 3 ∗ 105 . We assume that 𝑈∞ = 1.0 over a thin plate 10
𝜈 𝑠
long 𝑅𝑒𝑐𝑟𝑖𝑡𝑖𝑐𝑎𝑙 = 3 ∗ 105
𝑈∞ 𝑥𝐿 1.0∗𝑥𝐿
𝑅𝑒 = = = 3 ∗ 105
𝜈 1.6∗10−5
𝑥𝐿 = 4.8𝑚 . From Fig. 8.2, the following cases can be show.
(i) The thickness of B.L. (𝛿) increases with distance from leading edge
𝑥, 𝛿 − −→ 𝛿(𝑥)
(ii) 𝛿 Decreases as U increases.
(iii) 𝛿 Increases as kinematic viscosity  increases.
𝑈
(iv) 𝜏0 ≈ 𝜇 ( ) ; ℎ𝑒𝑛𝑐𝑒 𝜏0 decreases as x-increases. However, when B.L
𝛿
becomes turbulent.
𝑈∞ 𝑥
v) If 𝑅𝑒 = < 5 ∗ 105 B.L is laminar (Velocity distribution is parabolic)
𝜈
𝑈∞ 𝑥
If 𝑅𝑒 = > 5 ∗ 105 B.L is turbulent on that portion (Velocity
𝑉
distribution follows log law as a power law).
𝑈𝑥
vi) The critical value of at which B.L change from laminar to turbulent
𝑣
depends on:
 Turbulence in ambient flow.
 Surface roughness.
 Pressure gradient.
 Plate curvature.
 Temperature difference between fluid and boundary.
Figure 8.2: Comparison of flow past a sharp flat plate at low and high
Reynolds numbers
(a) Laminar, low –Re flow; (b) high-Re flow.
8.3 Displacement Thickness (𝜹𝒅 ) of B.L.

Consider the mass flow per unit depth across the vertical line y=0 and
𝑦 = 𝑦1 as shown in Fig. 8.3.
𝑦
A= actual mass flow between 0 and 𝑦1 = ∫0 1 𝜌𝑢 𝑑𝑦
B= theoretical mass flow between 0 and 𝑦1 if B.L were not present
𝑦
B = ∫0 1 𝜌∞ 𝑈∞ 𝑑𝑦
B-A= decrement in mass flow due to presence of B.L, i.e (missing mass
𝑦
flow)= ∫0 1( 𝜌∞ 𝑈∞ − 𝜌𝑢)𝑑𝑦 − − − (𝑎)
Express this missing mass flow as the product of 𝜌∞ 𝑈∞ and a heigh (𝛿𝑑 )
that is
Missing mass flow = 𝜌∞ 𝑈∞ 𝛿𝑑 − − − (𝑏)
Equating Eq's ( a & b)
𝑦
𝜌∞ 𝑈∞ 𝛿𝑑 = ∫0 1(𝜌∞ 𝑈∞ − 𝜌𝑢)𝑑𝑦
𝑦 𝜌𝑢
𝛿𝑑 = ∫0 1 (1 − ) 𝑑𝑦 𝑖𝑓 𝜌 = 𝜌∞ & 𝑦1 = 𝛿
𝜌 ∞ 𝑈∞
𝛿 𝑢
∴ 𝛿𝑑 = ∫0 ( 1 − ) 𝑑𝑦 (8.3)
𝑈 ∞
Physically;
1) Missing mass flow
2) Deflected the streamline upward throught a distance 𝛿𝑑
∴ 𝛿𝑑 is the distance through which the external inviscid flow is displaced by
the presence of the B.L.
External Streamline
y1 Streamline 𝛿𝑑
U U
Edge of
B.L.
u
𝛿

(a) (b)
Figure 8.3: (a) Hypothetical flow with no B.L.(inviscid case).
(b)Displacement thickness in actual flow with B.L.
8.4 Momentum Thickness ().

To understand the physical interpretation of  as the momentum
thickness. Consider the mass flow across a segment dy as in Fig. 8.4, given
by 𝑑𝑚 = 𝜌𝑢 𝑑𝑦 then
A= momentum flow across 𝑑𝑦 = 𝑑𝑚 𝑢 = 𝜌𝑢2 𝑑𝑦 − − − (𝑎)
If this same elemental mass flow were associated with the free- stream
velocity, where the velocity is 𝑈∞ , then
B= momentum flow at free stream velocity associated with dm
= 𝑑𝑚𝑈∞ = (𝜌𝑢 𝑑𝑦)𝑈∞ − − − − − (𝑏)
Hence, B-A= decrement in momentum flow (missing momentum flow)
associated with mass dm
= (𝜌𝑢 𝑈∞ − 𝜌𝑢2 )𝑑𝑦 = 𝜌𝑢 (𝑈∞ − 𝑢)𝑑𝑦 − − − (𝑐)
The total decrement in momentum flow across the vertical line from y=0
to 𝑦 = 𝑦1 is the integral of Eq. (c)
𝑦
∴ 𝑚𝑖𝑠𝑠𝑖𝑛𝑔 𝑚𝑜𝑚𝑒𝑛𝑡𝑢𝑚 𝑓𝑙𝑜𝑤 = ∫0 1 𝜌𝑢( 𝑈∞ − 𝑢)𝑑𝑦 − − − −(𝑑)
Assume that the missing momentum flow is the product of 𝜌∞ 𝑈∞2 and a
height . Then
2
𝑚𝑖𝑠𝑠𝑖𝑛𝑔 𝑚𝑜𝑚𝑒𝑛𝑡𝑢𝑚 𝑓𝑙𝑜𝑤 = 𝜌∞ 𝑈∞ 𝜃 − − − −(𝑒)
Equating Eq's( d& e)
𝑦
𝜌∞ U∞2 𝜃 = ∫0 1 𝜌𝑢( 𝑈∞ − 𝑢)𝑑𝑦
𝑦 𝜌𝑢 𝑢
∴ 𝜃 = ∫0 1 (1 − ) 𝑑𝑦 𝑖𝑓 𝜌 = 𝜌∞ 𝑎𝑛𝑑 𝑦1 = 𝛿
𝜌 ∞ 𝑈∞ U∞
𝛿 𝑢 𝑢
∴ 𝜃 = ∫0 (1 − ) 𝑑𝑦 (8.4)
U∞ 𝑈∞
∴ 𝜃 is an index that is proportional to the decrement in momentum flow due
to the presence of the B.L. It is the height of an ideal stream tube which is
carrying the missing momentum flow at free stream conditions.
8.5 Energy Thickness(𝜹𝒆 ).

Energy thickness is defined as the distance, measured perpendicular to
the boundary of the solid body, by which the boundary should be displaced to
compensate for the reduction in kinetic energy K.E. of the flowing fluid.
The mass of flow per second through the elementary strip= 𝜌𝑢 𝑑𝑦
1 1
K.E. of this fluid inside the B.L. = 𝑚̇𝑢2 = ( 𝜌𝑢 𝑑𝑦)𝑢2
2 2
K.E of the same mass of fluid before entering the B.L.
1
= ( 𝜌𝑢 𝑑𝑦)U∞2 ; Loss of K.E. through elementary strip is equal to
2
1 1
= (𝜌𝑢 𝑑𝑦)𝑈∞2 − (𝜌𝑢 𝑑𝑦)𝑢2
2 2
1
= (𝜌𝑢)(U∞2 − 𝑢 2 )𝑑𝑦
2
𝛿1
∴ 𝑇𝑜𝑡𝑎𝑙 𝑙𝑜𝑠𝑠 𝑜𝑓 𝐾. 𝐸. 𝑜𝑓 𝑓𝑙𝑢𝑖𝑑 = ∫0 𝜌𝑢( U∞2 − 𝑢2 )𝑑𝑦 − − − (𝑖)
2
Let 𝛿𝑒 = distance by which the plate is displaced to compensate for the
reduction in K.E. then loss of K.E. through 𝛿𝑒 of fluid flowing with velocity
U as follows
1
= (𝜌 𝑈∞ 𝛿𝑒 )U∞2 − − − − − (𝑖𝑖)
2
Equating Eq's (i) and (ii), we have
1 𝛿1
(𝜌𝑈 𝛿𝑒 )U∞2 = ∫0 𝜌𝑢 (U∞2 − 𝑢)𝑑𝑦
2 2
1 𝛿
𝑜𝑟 , 𝛿𝑒 =
U3∞
∫0 𝑢 (U 2 − u2 )dy
𝛿 𝑢 𝑢2
𝛿𝑒 = ∫0 (1 − ) 𝑑𝑦 (8.5)
U∞ U2∞
Ex.1
𝑢 𝑦
The velocity distribution in the B.L. is given by : = , where u is the
U∞ 𝛿
velocity at a distance y from the plate and 𝑢 = 𝑈 𝑎𝑡 𝑦 = 𝛿, 𝛿 𝑏𝑒𝑖𝑛𝑔 𝐵. 𝐿.
thickness. Find
(i)The displacement thickness (ii) the momentum thickness (iii) the energy
thickness and (iv) the value of 𝛿𝑑 /𝜃
Sol.
𝑢 𝑦
Velocity distribution =
U∞ 𝛿
(i) the displacement thickness 𝛿𝑑
𝛿 𝑢
𝛿𝑑 = ∫0 (1 − ) dy
U ∞
𝛿 𝑦
= ∫0 (1 − ) 𝑑𝑦
𝛿
𝑦2
= [𝑦 − ]
2𝛿
𝛿2 𝛿 𝛿
𝛿𝑑 = (𝛿 − )=𝛿− =
2𝛿 2 2
(ii) The momentum thickness, :
𝛿 𝑢 𝑢
𝜃 = ∫0 (1 − ) 𝑑𝑦
U∞ U∞
𝛿
𝛿 𝑢 𝑢 𝛿 𝑦 𝑦2 𝑦2 𝑦3
= ∫0 (1 − ) 𝑑𝑦 = ∫0 ( − 2 ) 𝑑𝑦 = [ − 2]
U∞ U∞ 𝛿 𝛿 2𝑦 3𝛿 0
𝛿2 𝛿3 𝛿 𝛿 𝛿
𝜃=[ − ]= − =
2𝛿 3𝛿 2 2 3 6
(iii) The energy thickness , 𝛿𝑒 :
𝛿 𝑢 𝑢2
𝛿𝑒 = ∫0 (1 − ) 𝑑𝑦
U∞ U2∞
𝛿𝑦 𝑦2 𝛿 𝑦 𝑦3
= ∫0 (1 − 2 ) 𝑑𝑦 = ∫0 ( − ) 𝑑𝑦
𝛿 𝛿 𝛿 𝛿3
𝛿
𝑦2 𝑦4 𝛿2 𝛿4 𝛿 𝛿 𝛿
𝛿𝑒 = [ − 3] = − = − =
2𝛿 4𝛿 0 2𝛿 4𝛿 3 2 4 4
𝛿
𝛿𝑑 2
(iv) The value of = 𝛿 =3
𝜃
6
Ex.2
The velocity distribution in the boundary layer is given by
𝑢 3𝑦 1 𝑦2
= − Calculate the following,
𝑈 2𝛿 2 𝛿2
𝛿𝑑
(i) The ratio of displacement thickness to B.L thickness ( ),
𝛿
𝜃
(ii) The ratio of momentum thickness to B.L thickness ( ).
𝛿
Sol.
𝛿 𝑢 𝛿 3𝑦 1 𝑦2
i) 𝛿𝑑 = ∫0 (1 − ) 𝑑𝑦 = ∫0 (1 − + 2 ) 𝑑𝑦
𝑈 2𝛿 2𝛿
𝛿
3 𝑦2 1 𝑦3 3 𝛿2 1 𝛿3
𝛿𝑑 = [𝑦 − + 2 ] = [𝛿 − + ∗ ]
2 2𝛿 2 3𝛿 0 4 𝛿 2 3𝛿 2
3 𝛿 5 𝛿𝑑 5
𝛿𝑑 = (𝛿 − 𝛿 + ) = 𝛿 =
4 6 12 𝛿 12
𝛿𝑢 𝑢
ii) 𝜃 = ∫0 (1 − ) 𝑑𝑦
𝑈 𝑈
𝛿 3𝑦 1 𝑦2 3𝑦 1 𝑦2
𝜃 = ∫0 ( − 2 ) (1 − + 2 ) 𝑑𝑦
2𝛿 2𝛿 2𝛿 2𝛿
𝛿 3𝑦 9 𝑦2 3 𝑦3 1 𝑦2 3 𝑦3 1 𝑦4
𝜃 = ∫0 ( − 2 + 3 − 2 + 3 − 4 ) 𝑑𝑦
2𝛿 4𝛿 4𝛿 2𝛿 4𝛿 4𝛿
𝛿 3𝑦 9 𝑦2 1 𝑦2 3 𝑦3 3 𝑦3 1 𝑦4
𝜃 = ∫0 [ − ( 2 + 2 ) + ( 3 + . 3 ) − 4 ] 𝑑𝑦
2𝛿 4𝛿 2𝛿 4𝛿 4 𝛿 4𝛿
𝛿 3𝑦 11 𝑦 2 3 𝑦3 1 𝑦4
𝜃 = ∫0 [ − + − ] 𝑑𝑦
2𝛿 4 𝛿2 2 𝛿3 4 𝛿4
𝛿
3 𝑦2 11 𝑦3 3 𝑦4 1 𝑦5
𝜃=[ ∗ − ∗ 2 + ∗ 3 − ∗ 4]
2 2𝛿 4 3𝛿 2 4𝛿 4 5𝛿 0
3 𝛿2 11 𝛿3 3 𝛿4 1 𝛿5
=[ ∗ − ∗ + ∗ − ∗ ]
2 2𝛿 4 3𝛿 2 2 4𝛿 3 4 5𝛿 4
3 11 3 1 19 𝜃 19
=( 𝛿− 𝛿+ 𝛿− 𝛿) = 𝛿−→ =
2 12 8 20 120 𝛿 120
8.6 Von Karman Integral Equation.

From the 2nd law of Newton's and momentum equation for a C.V as
shown in Fig.8.4 the following equations can be written
𝑑
∑𝐹 = ∫ 𝜌 𝑉 ⃗ 𝑑∀ + ∫ 𝜌𝑉 ⃗ (𝑉
⃗ . 𝑛)𝑑𝐴.
𝑑𝑡 𝑐.𝑣 𝑐.𝑠
𝜕𝑢 2 𝜕𝑣 𝜕𝑢
∑ 𝐹𝑥 = 𝜌 (𝑢 + 𝑑𝑥) 𝑑𝑦 − 𝜌𝑢2 𝑑𝑦 + 𝜌 (𝑣 + 𝑑𝑦) (𝑢 + 𝑑𝑦) 𝑑𝑥 − 𝜌𝑢 𝑣𝑑𝑥
𝜕𝑥 𝜕𝑦 𝜕𝑦
𝜕𝑢 𝜕𝑢 2
∑ 𝐹𝑥 = 𝜌 [𝑢2 + 2𝑢 𝑑𝑥 + ( 𝑑𝑥) ] 𝑑𝑦 − 𝜌𝑢2 𝑑𝑦 + 𝜌𝑣𝑢𝑑𝑥 +
𝜕𝑥 𝜕𝑥
𝜕𝑢 𝜕𝑣 𝜕𝑣 𝜕𝑢
𝜌𝑣 𝑑𝑦𝑑𝑥 + 𝜌 𝑢 𝑑𝑦 𝑑𝑥 + 𝜌 . (𝑑𝑦)2 ∗ 𝑑𝑥 − 𝜌𝑢𝑣 𝑑𝑥
𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦
𝜕𝑢 𝜕𝑣 𝜕𝑢
∑ 𝐹𝑥 = 𝜌 (2 𝑢 + 𝑢 + 𝑣 ) 𝑑𝑦 𝑑𝑥
𝜕𝑥 𝜕𝑦 𝜕𝑦
𝜕𝑣 𝜕𝑢
After disregarding second – order terms & =− from C.E.
𝜕𝑦 𝜕𝑥
𝜕𝑢 𝜕𝑢
∴ ∑ 𝐹𝑥 = 𝜌 ( 𝑢 + 𝑣 ) 𝑑𝑥𝑑𝑦 (8.6)
𝜕𝑥 𝜕𝑦
The summation of forces on C.V may be written as
𝜕𝑝 𝜕𝜏𝑥
∑ 𝐹𝑥 = 𝑝 𝑑𝑦 − (𝑝 + 𝑑𝑥) 𝑑𝑦 + (𝜏𝑥 + 𝑑𝑦) 𝑑𝑥 − 𝜏𝑥 𝑑𝑥
𝜕𝑥 𝜕𝑦
∑ 𝐹𝑥 = 𝑝 𝑑𝑦 − 𝑝𝑑𝑦 − 𝑑𝑥𝑑𝑦 + 𝜏𝑥 𝑑𝑦 + 𝑑𝑦𝑑𝑥 − 𝜏𝑥 𝑑𝑥
𝜕𝑥 𝜕𝑦
∑ 𝐹𝑥 = (− + ) 𝑑𝑥𝑑𝑦 (8.7)
𝜕𝑥 𝜕𝑦
Equating Eq's (8.6 & 8.7)
𝜕𝑢 𝜕𝑢 𝜕𝑝 𝜕𝜏𝑥
𝜌 (𝑢 + 𝑣 ) 𝑑𝑥 𝑑𝑦 = (− + ) 𝑑𝑥𝑑𝑦 (8.8)
𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑦
𝜕𝑣
𝜌 (𝑣 + 𝑑𝑦)
𝜕𝑦
𝜕𝜏𝑥
𝜌 (𝜏𝑥 + 𝑑𝑦)
𝜕𝑦
3 2
𝜕𝑢
𝜌 (𝑢 + 𝑑𝑥)
𝜕𝑥
u
p 𝜕𝑝
(𝑝 + 𝑑𝑥)
v 𝜕𝑥
4 1
τx
Figure 8.4: Distribution of pressure forces on control volume.

The shear stress is very nearly equal to

𝜕𝑢
𝜏𝑥 = 𝜇 ; Substituting in Eq.8.8 gives
𝜕𝑦
𝜕𝑢 𝜕𝑢 1 𝜕𝑝 𝜇 𝜕2𝑢
𝜌(𝑢 + 𝑣 )=− +( ) (8.9)
𝜕𝑥 𝜕𝑦 𝜌 𝜕𝑥 𝜌 𝜕𝑦 2
Table 8.1 Masses and momentum fluxes on control volume faces

Surface Mass flux Flux of x-momentum
(1-2) 𝜕𝑢 2
𝜕𝑢
𝜌 (𝑢 + 𝑑𝑥) 𝑑𝑦 𝜌 (𝑢 + 𝑑𝑥) 𝑑𝑦
𝜕𝑥 𝜕𝑥
(3-4) 𝜌𝑢 𝑑𝑦 𝜌𝑢2 𝑑𝑦
(2-3) 𝜕𝑣 𝜕𝑣 𝜕𝑢
𝜌 (𝑣 + 𝑑𝑦) 𝑑𝑥 𝜌 (𝑣 + 𝑑𝑦) (𝑢 𝑑𝑦) 𝑑𝑥
𝜕𝑦 𝜕𝑦 𝜕𝑦
(4-1) 𝜌𝑣 𝑑𝑥 𝜌𝑣 𝑢𝑑𝑥
From the Integral method of momentum equation for Von Karman

integral as follows
ℎ 𝜕𝑢 𝜕𝑢 ℎ 1 𝑑𝑝 ℎ 𝜕2𝑢
∫0 (𝑢 𝜕𝑥 + 𝑣 𝜕𝑦) 𝑑𝑦 = ∫0 − 𝜌 𝑑𝑥 𝑑𝑦 + ∫0 𝑣 𝜕𝑦 2
𝑑𝑦 (8.10)
Where h is an undefined distance from the wall to outside the boundary layer.
Integrating the second term in the integral by parts, we have
ℎ 𝜕𝑢 ℎ 𝜕𝑣
[∫0 𝑣 𝑑𝑦] = [𝑣𝑢]ℎ0 − ∫0 𝑢 𝑑𝑦 − − − −(𝑎)
𝜕𝑦 𝜕𝑦
From C.E, v at y=h is given by ;
ℎ 𝜕𝑣 ℎ 𝜕𝑢
𝑣 = ∫0 𝑑𝑦 = − ∫0 𝑑𝑦 − − − −(𝑏)
𝜕𝑦 𝜕𝑥
Substitute Eq. (b in a) and 𝑢 = 𝑈∞ 𝑎𝑡 𝑦 = ℎ
ℎ 𝜕𝑢 ℎ 𝜕𝑢 ℎ 𝜕𝑢
[∫0 𝑣 𝑑𝑦] = −U∞ ∫0 𝑑𝑦 + ∫0 𝑢 𝑑𝑦 − − − −(𝑐)
𝜕𝑦 𝜕𝑥 𝜕𝑥
𝜕𝑝
CASE (A) =0
𝜕𝑥
𝜕𝑝
Substitute Eq. (c) in Eq's (8.9& 8.10) and neglecting =0
𝜕𝑥
ℎ 𝜕𝑢 𝜕𝑢 𝜕𝑢
∫0 𝜌 (𝑢 𝜕𝑥 + 𝑢 𝜕𝑥 − 𝑈∞ 𝜕𝑥 ) 𝑑𝑦 = −𝜏0
ℎ 𝜕𝑢 𝜕𝑢
∫0 𝜌 ( 2𝑢 𝜕𝑥 − U∞ 𝜕𝑥 ) 𝑑𝑦 = −𝜏0 − − − − − (𝑒)
ℎ 𝜕𝑢 𝜕𝑢
∫0 𝜌 ( U∞ 𝜕𝑥 − 2𝑢 𝜕𝑥 ) 𝑑𝑦 = 𝜏0 − − − − − (𝑓)
𝜕𝑢2 𝜕𝑢
Since = 2𝑢
𝜕𝑥 𝜕𝑥
ℎ 𝜕𝑢 𝜕
∫0 𝜌 (( U∞ 𝜕𝑥 ) − 𝜕𝑥 (𝑢2 )) 𝑑𝑦 = 𝜏0
𝑑 ℎ
𝜏0 = 𝜌
𝑑𝑥
∫0 𝑢 ( U∞ − 𝑢)𝑑𝑦 𝑠𝑖𝑛𝑐𝑒 ℎ = 𝛿
𝑑 𝛿
𝜏0 = 𝜌
𝑑𝑥
∫0 𝑢 ( U∞ − 𝑢)𝑑𝑦
𝜏0 𝑑 𝛿 𝑢 𝑢
= ∫0 (1 − ) 𝑑𝑦 (8.11)
𝜌U2∞ 𝑑𝑥 U∞ 𝑈∞
𝜏0 𝑑𝜃
= (8.12)
𝜌U2∞ 𝑑𝑥
𝛿 𝑢 𝑢
Since 𝜃 = ∫0 (1 − ) 𝑑𝑦 (8.13)
U∞ U∞
Eq. 8.13 is the Von Karman equation without pressure gradient dp/dx=0
𝜕𝑝
CASE (B) ≠0
𝜕𝑥
When the pressure gradient in Eq's (8.9 or 8.10) is included and from adding
it in Eq. (e) as followes
1 𝑑𝑝 𝑑𝑈∞
− = U∞ 𝑓𝑟𝑜𝑚 𝐵𝑜𝑢𝑟𝑛𝑎𝑙𝑙𝑖′𝑠 𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛
𝜌 𝑑𝑥 𝑑𝑥
ℎ 𝜕𝑢 𝜕𝑢 𝑑𝑈∞
Now∫0 𝜌 (− U∞ + 2𝑢 − U∞ ) 𝑑𝑦 = −𝜏0 (8.14)
𝜕𝑥 𝜕𝑥 𝑑𝑥
𝜕 𝜕𝑢 𝜕𝑈∞
(𝑢 U∞ ) = 𝑈∞ + 𝑢
𝜕𝑥 𝜕𝑥 𝜕𝑥
𝜕𝑢 𝜕 𝜕𝑈∞
𝑈∞ = (𝑢 U∞ ) − 𝑢 Substituted in Eq. 8.14 after multiply by (-1)
𝜕𝑥 𝜕𝑥 𝜕𝑥
ℎ 𝜕𝑢 𝜕U∞ 𝜕U∞ 𝜕 𝜏0
∫0 [ 2𝑢 𝜕𝑥 − U∞ + 𝑢 − ( 𝑢 U∞ )] 𝑑𝑦 = −
𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜌
𝜕𝑢 𝜕𝑢2
In above eqn. 2𝑢 =
𝜕𝑥 𝜕𝑥
ℎ 𝜕 𝑑U∞ ℎ 𝜏
∫0 𝜕𝑥 [ 𝑢 ( 𝑢 − U∞ ] 𝑑𝑦 −
𝑑𝑥
∫0 ( U∞ − 𝑢 )𝑑𝑦 = − 𝜌0
𝛿 𝑢 𝑢
𝜃 = ∫0 (1 − ) 𝑑𝑦
U∞ U∞
𝛿 𝑢
𝛿𝑑 = ∫0 (1 − ) 𝑑𝑦
U∞
𝜏0 𝑑 2 𝑑U
∴ = (U∞ 𝜃) + U∞ 𝛿𝑑 ∞
𝜏0 2 𝑑𝜃 𝑑U∞ 𝑑U∞
= U∞ + 𝜃 ∗ 2U∞ + U∞ 𝛿𝑑
𝜌 𝑑𝑥 𝑑𝑥 𝑑𝑥
𝜏0 2 𝑑𝜃 𝑑U∞
= U∞ + (2𝜃 + 𝛿𝑑 )U∞ (8.15)
We assume
𝜏0 𝛿𝑑
𝐶𝑓 = 1 𝑎𝑛𝑑 𝐻 =
2
𝜌 U2∞ 𝜃
𝐶𝑓 is the friction factor and H is the shape factor.

𝐶𝑓 𝑑𝜃 1 𝑑U∞
∴ = + 𝜃(2 + 𝐻 ) (8.16)
2 𝑑𝑥 U∞ 𝑑𝑥
𝑑𝑝 𝑑U∞
From Eq's (8.15 and 8.16) the case of =0& =0
𝑑𝑥 𝑑𝑥
No gradient of pressure along the x-axis
𝜏0 𝑑𝜃 1
= = 𝑐𝑓 (8.17)
𝜌U2∞ 𝑑𝑥 2
8.7 Approximate Solution to the Laminar B.L.

We have four conditions that proposed velocity profile should satisfy on
flat plat with zero pressure gradients.
𝑢 = 0 𝑎𝑡 𝑦 = 0
𝑢 = U∞ 𝑎𝑡 𝑦 = 𝛿
𝜕𝑢
= 0 𝑎𝑡 𝑦 = 𝛿
𝜕𝑥
𝜕2 𝑢
= 0 𝑎𝑡 𝑦 = 0
𝜕𝑦 2
𝑢
Let = 𝐴 + 𝐵𝑦 + 𝐶𝑦 2 + 𝐷𝑦 3 is a cubic polynomial will satisfy the four
U∞
conditions.
From above conditions
At y=0 , u=0
𝑢
=0=𝐴 − − − − − −(𝑎)
U∞
At 𝑦 = 𝛿 𝑢 = U∞
𝑈
= 1 = 𝐵𝛿 + 𝐶𝛿 2 + 𝐷𝛿 3 − − − − − (𝑏)
U∞
𝜕𝑢
At 𝑦 = 𝛿 =0
𝜕𝑦
𝜕𝑢
= 𝐵 + 2𝑐𝑦 + 3𝐷𝑦 2
𝜕𝑦
𝐵 + 2𝐶𝛿 + 3𝐷𝛿 2 = 0 − − − − − (𝑐)
At y=0
𝜕2𝑦
= 0 = 2𝐶 + 6𝐷𝑦 − −−→ 𝐶 = 0 − − − −(𝑑)
𝜕𝑦 2
From Eq. (c)  𝐵 = −3𝐷𝛿 2 − − − − − (𝑒)
1
From Eq. (b) 𝐵 = − 𝐷𝛿 2 − − − −(𝑓)
𝛿
Equating (e & f )
1 1 3
−3𝐷𝛿 2 = − 𝐷𝛿 2 − − − −→ 𝐷 = − ; 𝐵=
𝛿 2𝛿 3 2𝛿
Hence a good approximation for the velocity profile in a laminar flow is

𝑢 3 𝑦 1 𝑦 3
= ( )− ( )
𝑈∞ 2 𝛿 2 𝛿
Let us now use this velocity profile to find 𝛿(𝑥)𝑎𝑛𝑑 𝜏0 (𝑥). From Von
Karman's integral Eq's (8.11 & 8.12)
= ∫0 (1 − ) 𝑑𝑦
𝜌U2∞ 𝑑𝑥 U∞ U∞
𝑢 3𝑦 𝑦3
Assume the velocity take the following form = − in the
𝑈∞ 2𝛿 2𝛿 3
B.L. Calculate the thickness and the wall shear.
= ∫ (1 − ) 𝑑𝑦
𝜌U2∞ 𝑑𝑥 0 U∞ U∞
𝜏0 𝑑 𝛿 3𝑦 𝑦3 3𝑦 𝑦3
= ∫0 (2 𝛿 − 2𝛿 3) (1 − 2𝛿 + 2𝛿 3) 𝑑𝑦
𝜌U2∞ 𝑑𝑥
𝜏0 𝑑 𝛿 3𝑦 9𝑦 2 3𝑦 4 𝑦3 3𝑦 4 𝑦6
= ∫0 (2 𝛿 − 4𝛿 2 + 4𝛿 4 − 2𝛿 3 + 4𝛿 4 − 4𝛿 6 ) 𝑑𝑦
𝜌U2∞ 𝑑𝑥
𝜏0 𝑑 𝛿 3𝑦 9𝑦 2 6𝑦 4 𝑦3 𝑦6
2 = ∫ ( − 4𝛿2 + 4𝛿 4 − 2𝛿 3 − 4𝛿 6) 𝑑𝑦
𝜌𝑈∞ 𝑑𝑥 0 2 𝛿
𝛿
𝜏0 𝑑 3𝑦 2 9𝑦 3 6𝑦 5 𝑦4 𝑦7
= [ − 2 + 4 − 3 − 6]
𝜌U2∞ 𝑑𝑥 4𝛿 12𝛿 20𝛿 8𝛿 28𝛿 0
𝜏0 𝑑 3 9 6 𝛿 𝛿
= [ 𝛿− 𝛿+ 𝛿− − ]
𝜌U2∞ 𝑑𝑥 4 12 20 8 28
𝜏0 𝑑𝛿 𝑑𝛿
= (0.1392) = 0.1392
𝜌U2∞ 𝑑𝑥 𝑑𝑥
2
∴ 𝜏0 = 0.1392 𝜌U∞ 𝑑𝛿/𝑑𝑥 (8.18)
𝜕𝑢
At the wall we know that 𝜏0 = 𝜇 | or using the cubic profile
𝜕𝑦 𝑦=0
𝜕𝑢 3
| = 𝐵 + 2𝐶𝑦 + 3𝐷𝑦 2 = 𝐵 = ( )U∞
𝜕𝑦 𝑦=0 2𝛿
3
∴ 𝜏0 = 𝜇 (𝑈∞ ) (8.19)
2𝛿
Equating the foregoing expression (8.18 & 8.19) for 𝜏0 (𝑥), we find that
2 𝑑𝛿 3
0.139𝜌U∞ = 𝜇 (U∞ )
𝑑𝑥 2𝛿
3
2
𝜇 𝜈
𝛿 𝑑𝛿 = 𝑑𝑥 = 10.8 𝑑𝑥 (8.20)
0.139𝜌𝑈∞ 𝑈∞
From using at 𝛿 = 0 𝑎𝑡 𝑥 = 0 (the leading edge) Eq. 8.20 is integrated to
give
√𝑥
𝛿 = 4.65√𝜈𝑥/𝑈∞ 𝑚𝑢𝑙𝑡𝑖𝑝𝑙𝑦 𝑏𝑦
√𝑥
𝑥
𝛿 = 4.65 (8.21)
√𝑅𝑒𝑥
Where 𝑅𝑒𝑥 is the local Reynolds number. Substituted Eq. 8.21 in Eq.
2 𝜈
8.19 giving the wall shear as 𝜏0 = 0.323 𝜌U∞ √ 𝑥 U∞
0.323𝜌U2∞
𝜏0 = (8.22)
√𝑅𝑒𝑥
1 2
The shearing stress is made dimensionless by dividing by 𝜌U∞ . The local
2
skin friction coefficient 𝐶𝑓 its
𝜏0 0.646 0.646
𝐶𝑓 = 1 = = (8.23)
2
𝜌U2∞ 𝑥
√𝑈∞ 𝜈 √𝑅𝑒𝑥
If the wall shear is integrated over the length L, the result per unit width
is the drag force.
𝐿
𝐹𝐷 = ∫0 𝜏0 𝑑𝑥 = 0.646 𝜌U∞ √U∞ 𝐿𝜈 Where 𝜏0 from Eq. 8.22
0.646 𝜌U2∞ 𝐿
𝐹𝐷 = (8.24)
√𝑅𝑒𝐿
𝐿 2 1
Or 𝐹𝐷 = ∫0 𝜏0 𝑑𝑥 = 𝜏0 . 𝐿 = 𝐶𝐹 . 𝜌U∞ .𝐿 since 𝜏0 from Eq. 8.23
2
𝐹 2 .𝐿
0.646 𝜌𝑈∞ 1.292
∴ 𝐶𝐹 = 1 𝐷2 =1 = (8.25)
𝜌U∞ .𝐿 𝜌 U2∞ .𝐿√𝑅𝑒𝐿 √𝑅𝑒𝐿
2 2
Where 𝑅𝑒𝐿 is the Reynolds number at the end of flat plate.
Ex.3
Assume a parabolic velocity profile and calculate the B.L thickness
and the wall shear. Compare with those calculate above.
Sol.
The parabolic velocity profile is assumed to be
𝑢
= 𝐴 + 𝐵𝑦 + 𝐶𝑦 2
U∞
With three conditions
𝜕𝑢
𝑢 = 0 𝑎𝑡 𝑦 = 0; 𝑢 = U∞ 𝑎𝑡 𝑦 = 𝛿; = 0 𝑎𝑡 𝑦 = 𝛿 ∴ 𝐴 = 0
𝜕𝑦
1 = 𝐴 + 𝐵𝛿 + 𝐶𝛿 2 = 𝐵𝛿 + 𝐶𝛿 2 ;
0 = 𝐵 + 𝐶 ∗ 2𝛿
2 1
Then A=0; 𝐵 = ; 𝐶=− . The velocity profile is
𝛿 𝛿2
𝑢 𝑦 𝑦2
= 2 − − − − − − −(𝑎)
U∞ 𝛿 𝛿2
This is substituted into Von Karman's integral equation (8.11 ) to obtain
𝜏0 𝑑 𝛿 2𝑦 𝑦2 2𝑦 𝑦2
= ∫ ( 𝛿 − 𝛿 2 ) [1 − + ] 𝑑𝑦
𝜌U2∞ 𝑑𝑥 0 𝛿 𝛿2
𝑑 𝛿 2𝑦 4𝑦 2 2𝑦 3 𝑦2 2𝑦 3 𝑦4
= ∫ [ − + − 2 + − ] 𝑑𝑦
𝑑𝑥 0 𝛿 𝛿2 𝛿3 𝛿 𝛿3 𝛿4
𝛿
𝑑 2𝑦 2 4𝑦 3 2𝑦 4 𝑦3 2𝑦 4 𝑦5
= [ − 2 + 3 − 2 + 3 − 4]
𝑑𝑥 2𝛿 3𝛿 4𝛿 3𝛿 4𝛿 5𝛿 0
𝑑 4 1 1 1 1
= [𝛿 − 𝛿 + 𝛿 − 𝛿 + 𝛿 − 𝛿]
𝑑𝑥 3 2 3 2 5
𝑑 4 1 1 1 1
= [1 − + − + − ] 𝛿
𝑑𝑥 3 2 3 2 5
𝑑 30−25−3 𝑑 2
= ( )𝛿 = 𝛿
𝑑𝑥 15 𝑑𝑥 15
2 2 𝑑𝛿
∴ 𝜏0 = 𝜌U∞ − − − − − (𝑏)
15 𝑑𝑥
𝜕𝑢
We also use 𝜏0 = 𝜇 |
𝜕𝑦 𝑦=0
From Eq. (a)
2
𝜏0 = 𝜇 U ∞ − − − − − (𝑐)
𝛿
Equating Eq's ( b & c) we obtain
2 2 𝑑𝛿 2
𝜌U∞ = 𝜇𝑈∞
15 𝑑𝑥 𝛿
𝜈
𝛿 𝑑𝛿 = 15 𝑑𝑥 Using 𝛿 = 0 𝑎𝑡 𝑥 = 0 𝑎𝑓𝑡𝑒𝑟 𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑡𝑖𝑜𝑛
U∞
𝛿 𝑥 𝜈
∫0 𝛿𝑑𝛿 = ∫0 15 U 𝑑𝑥
∞
𝛿2 𝜈
= 15 𝑥
2 𝑈∞
𝜈𝑥
∴ 𝛿 = 5.48 √ − − − −(𝑑)
𝑈∞
This is 18% higher than the value using the cubic profile, the wall shear is
found to be
2𝜇U∞
𝜏0 = Substitute Eq. (d)
𝛿
2𝜇U∞ U∞
𝜏0 = √
5.43 𝜈𝑥
𝜈 2
0.365𝜌𝑈∞
2
𝜏0 = 0.365 𝜌𝑈∞ √ =
𝑥 𝑈∞ √𝑅𝑒𝑥
This is 13% higher than the value using the cubic velocity profile.
8.8 Solution of Turbulent B.L. Power – Law Form.

The power– law form is
7 𝑅𝑒𝑥 < 107
𝑢 𝑦 1/𝑛
= ( ) n= 8 107 < 𝑅𝑒𝑥 ≤ 108 (8.26)
U∞ 𝛿
9 108 < 𝑅𝑒𝑥 ≤ 109
The Balsius formula is an empirical relation for the local friction coefficient
1
𝜈 4
𝐶𝑓 = 0.046 ( ) (8.27)
U∞ 𝛿
𝜏0 21
We have 𝐶𝑓 = 1 − −→ 𝜏0 = 𝐶𝑓 𝜌U∞
2
𝜌U2∞ 2
𝜈 1/4
2
∴ 𝜏0 = 0.023 𝜌U∞ ( ) (8.28)
U∞ 𝛿
Substitute the velocity profile in Von Karman's integral Eq. 8.11 with 𝑅𝑒𝑥 <
167
1
𝜏0 𝑑 𝛿 𝑦 1/7 𝑦 7
= ∫ [( ) (1 − ( ) )] 𝑑𝑦
𝜌U2∞ 𝑑𝑥 0 𝛿 𝛿
7 2 𝑑𝛿
𝜏0 = 𝜌U∞ (8.29)
72 𝑑𝑥
Equating Eq's (8.28 & 8.29) for 𝜏0 , we find that
1
1
𝜈 4
𝛿 𝑑𝛿 = 0.237 (
4 ) 𝑑𝑥
U∞
Assuming a turbulent flow from the leading edge L>>𝑥𝐿 , from integration
the above eqn.
5
𝛿4 𝜈 1/4
5 = 0.237 ( ) 𝑥
U∞
4
After taking the root (4/5) and multiplied by (x/x)1/5 gives the following
𝛿 = 0.38 𝑥 (𝑅𝑒𝑥 )−1/5 (8.30)
Substituting Eq. 8.30 for 𝛿 in to Eq. 8.27 we find that
𝐶𝑓 = 0.059(𝑅𝑒 )−1/5 𝑓𝑜𝑟 𝑅𝑒 < 107 (8.31)
1 2 (𝐿𝑊)
The drag force = 𝐹𝐷 = 𝜏0 𝐴 = 𝐶𝑓 𝜌U∞ (8.32)
2
Where W is the width of plate and Cf from equation 8.31.
Ex.4
Estimate the boundary layer thickness at the end of a 4-m-long flat
𝑚
surface if the free-stream velocity U∞ = 5 . Use atmospheric air at 30𝐶 ° .
𝑠
And predict the drag force if the surface is 5m wide
a) Neglect the laminar portion of the flow
b) Account for the Laminar portion using 𝑅𝑒𝑐𝑟𝑖𝑡. = 5 ∗ 105 .
Sol.
a) Let us first assume turbulent flow from the leading edge. The B.L
thickness is given by Eq. 8.30. It is
𝛿 = 0.38 𝑥 𝑅𝑒𝑥 −1/5

5∗4
𝛿 = 0.38 ∗ 4 ∗ ( )−1/5 = 0.092𝑚
1.6∗10−5
U∞ .𝐿
The drag force is using Eq's (8.31 & 8.32); 𝑅𝑒 =
𝜈
1 2 (𝐿𝑊) −1/5 1 2
𝐹𝐷 = 𝐶𝑓 ∗ 𝜌U∞ = 0.059(𝑅𝑒 ) ∗ 𝜌U∞ (𝐿. 𝑊)
2 2
5∗4 −1/5 1
𝐹𝐷 = 0.059 ( ) ∗ ∗ 1.16 ∗ 52 ∗ 4 ∗ 5 = 1.032𝑁
1.6∗10−5 2
5∗4 6
𝑅𝑒𝐿 = = 1.25 ∗ 10 . Hence, the calculation is acceptable
1.6∗10−5
The distance is found as follows
U∞ 𝑥𝐿
𝑅𝑒𝑐𝑟𝑖𝑡. = 5 ∗ 105 =
𝜈
5∗105 ∗1.6∗10−5
∴ 𝑥𝐿 = = 1.6 𝑚
5
The B.L thickness at 𝑥𝐿 is found from Eq. 8.21 with
𝑥𝐿 4.65∗1.6
𝛿 = 4.65 = = 0.0105 𝑚
√𝑅𝑒𝑥 √5∗105
To find the origin of turbulent flow, using equation of B.L thickness in
turbulent as
4
−1/5 0.38𝑥′ 0.38𝑥 5
𝛿 = 0.38 𝑥 (𝑅𝑒𝑥 ) = 1 = 1
U 𝑥 5 𝑈
( ∞
𝜈
) ( 𝜈∞ )5
𝛿 U∞ 1/5
∴ 𝑥′4/5 = ( ) where 𝛿 = the same at end of L.B
0.38 𝑉
0.0105 5/4 5 1/4
𝑥′ = ( ) ( −5 ) = 0.2663𝑚
0.38 1.6∗10
The distance 𝑥𝑡 as in figure is then
𝑥𝑡 = 𝐿 − 𝑥𝐿 + 𝑥′ = 4 − 1.6 + 0.266 = 2.666𝑀
To find the B.L thickness at the end of plate using Eq. 8.30
5∗2.666 −1/5
𝛿 = 0.38 𝑥 𝑅𝑒 −1/5 = 0.38 ∗ 2.666 ∗ ( )
1.6∗10−5
𝛿 = 0.0662𝑚
The value of part (a) is 28% to high when compared with this more accurate
value. T.B.L
L.B.L
xL x'
xt
Problems.
P8.1 How far from the leading edge can turbulence be expected on an airfoil
traveling at 100 m/s if the kinematic viscosity at different elevation as
followes,
𝜈1 = 1.169 ∗ 10−5 𝑚/𝑠 2 𝑎𝑡 𝑇 = −20𝐶 °
𝜈2 = 1.087 ∗ 10−5 𝑚/𝑠 2 𝑎𝑡 𝑇 = −30°
𝜈3 = 1.008 ∗ 10−5 𝑚/𝑠 2 𝑎𝑡 𝑇 = −40° use 𝑅𝑒𝑐𝑟𝑖𝑡 = 6 ∗ 105 and
assume flat plate with zero pressure gradient.
[xL1=0.1014 m, xL2=0.065 m, xL3=0.06]
𝜋𝑦
P8.2 Assume that 𝑢 = U∞ sin ( ) in a zero pressure gradient boundary
2𝛿
Layer. Calculate; (𝑎) 𝛿(𝑥), (𝑏) 𝜏0 (𝑥).
P8.3 A boundary Layer profile is approximated with

𝑦
𝑢 = 3𝑈∞ 0 < 𝑦 ≤ 𝛿/6
0
𝑦 1 𝛿 𝛿
𝑢 = 𝑈∞ ( + ) <𝑦≤
𝛿 3 6 2
𝑦 2 𝛿
𝑢 = 𝑈∞ ( + ) <𝑦≤𝛿
3𝛿 3 2
Determine 𝛿(𝑥)𝑎𝑛𝑑 𝜏0 (𝑥).
P8.4 A Laminar flow is maintained in a B.L on a 6m-long 5m-wide flat plate

with 20C atmospheric air flowing at 4 m/s. Assuming a parabolic
profile. Calculate
(a) 𝛿 at x=6m. [5.84 mm]
(b) 𝜏0 𝑎𝑡 𝑥 = 6𝑚. [0.025 N/m2]
(c) The drag force on one side. [0.75 N]
P8.5 The velocity profile at a given x-location in B.L is assumed to be

2𝑦 𝑦2
𝑢(𝑦) = 10 ( − ). A stream is 2 cm from the flat plate at the
𝛿 𝛿2
leading edge. How far is it from the plate when x=3m (i.e what is h)?
Also, calculate the displacement thickness at x =3m. Compare the
displacement thickness to (h=2 cm).
[𝟐. 𝟖𝟗 𝒄𝒎, 𝜹𝒅 = 𝟎. 𝟖𝟗𝟒 𝒎𝒎]
P8.6 Atmospheric air at 20C flows at 10 m/s over a 2m-long 4m-wide flat
plat. Calculate the maximum B.L thickness and drag force on one side
assuming
a) Laminar flow over the entire length. [𝜹 = 𝟐. 𝟏𝟒 𝒎𝒎, 𝑭𝒅 = 𝟏. 𝟑𝟖 𝑵]

b) Turbulent flow over the entire length. [𝜹 = 𝟒𝟓 𝒎𝒎, 𝑭𝒅 = 𝟏. 𝟔𝟖𝟕 𝑵]
P8.7 Fluid flows over a flat plat at 20 m/s. Neglecting the laminar portion
then determine  and 𝜏0 𝑎𝑡 𝑥 = 6 𝑚 𝑖𝑓 𝑡ℎ𝑒 𝑓𝑙𝑢𝑖𝑑 𝑖𝑠:
a- Atmospheric air at 20C [=95 mm, 0=17.64 N/m2]
b- Water at 20C [=59 mm, 0=279.44 N/m2]
P8.8 Find the shear stress and the thickness of the boundary layer
i- at the trailing edge of smooth flat plate 3.0 m wide and 0.6 m long
parallel to flow immersed in 15 C° water flowing at the an undisturbed
velocity of 0.9 m/s. Assume a laminar boundary layer over the whole
plate. [=4 mm, 0=0.38 N/m2]
ii- At the center. [=2.86 mm, 0=0.537N/m2]
iii- Find the total friction drag on one side of the plate. [Fd=0.684 N]
2
P8.9 i- Given the general equation for a parabola u=ay +by+c, drive the
dimensional velocity distribution equation (u/U).
ii- Determine the shear stress at 150 mm and 300 mm back from the
leading edge of the plate placed longitudinally in oil (S.G=0.925,
=0.73*10-4m2/s) flowing with undisturbed U=0.6 m/s.
[0=3.052 N/m2, 0=2.16 N/m2]
P8.10 A 80 m long streamlined train has 2.3 m high sides and 2.4 m wide
top. Determine the power required to overcome the skin – friction drag
when the train is traveling at 15 m/s through standard atmosphere at
sea level, assuming the drag on the sides and top to be equal to that on
one side of a flat plate 5m wide and 80 m long. [P=3.825 kW]
Chapter 9 Turbomachinery 211
CHAPTER 9
Turbomachinery
9.1 Introduction.
The pump family is a machine which is designed to add energy to the
fluid, but the turbines family which extract energy from the fluid. Both types
are usually connected to a rotating shaft, hence this is the turbomachinery.
Machine which deliver liquid are called pump, if machine delivers gases can
be classified into the categories as follows
Fan, pressure rise very

small (few inches of water)
Three types of Blower, up to 1 atm.

machines
Compressor above 1 atm.
If machine delivers liquid can be classified into two type as follows
Positive–displacement pumps
(PDP's) are limited by (low
pressure, small volume).
pumps are used
today
Rotodynamic or momentum
change pumps. (From small
volume to very large volume
and higher pressures).
Reciprocating piston or plunger
External gear pump
Double-screw pump
PDP's Sliding vane pump
Three lobe pump
Double circumferential pump
Flexible – tube squeegee pump
Centrifugal pumps
Rotodynamic Jet pumps

Axial pumps
pumps
Electromagnetic
Special designs
pumps
Fluid actuated
pump
All types of pumps can be shows in Fig. (9.1), but this chapter is concern to
the centrifugal pump type, which are explain in details in the next sections.
Figure (9.1) Schematic design of positive-displacement pumps: (a) reciprocating piston or

plunger, (b) external gear pump, (c) double-screw pump, (d) sliding vane, (e) three lobe
pump, ( f ) double circumferential piston, (g) flexible-tube squeegee.
9.2 Centrifugal Pump.

The main components of centrifugal pumps are
The impeller
The casing
Centrifugal Pump The drive shaft with gland and packing
Suction pipe with one way valve (foot

valve)
Delivery pipe with delivery valve

The liquid enters the eye of the impeller axially due to the suction created by
the impeller motion as shown in Fig (9.2). The impeller blades guide the fluid
and impart momentum to the fluid, which increase the total head (or pressure)
of the fluid. The casing can be simple volute type or diffuser can be used as
desired as in Fig. (9.3).
Figure 9.2 Volute type centrifugal pump.
Figure 9.3 Diffuser pump.
The volute is a spiral casing of gradually increasing cross section. A

part of the kinetic energy in the fluid is converted to pressure in the casing.
9.2.1 Impeller.
The impeller consists of two disc plates with blades mounted
perpendicularly on its surface in between. The blades of the rotating impeller
transfer energy to the fluid there by increasing pressure and velocity. The fluid
enters the impeller eye through fluid sucking then flow through the impeller
channels formed by the blades between hub plate and shroud plate. Fig.(9.4)
shows the impeller components and the flow direction relatively to the
impeller.
Figure (9.4) the impeller components and the flow directions.
The blades may be of three different orientations
Radial type 2=90°
Blade types Backward type 290°
Forward type 290°
Fig. (9.5) shows the blade type orientations can be used.
Figure (9.5) Different blade arrangements.

Also the impeller can be classified as follows,
Open impeller type (when the blades

mounted perpendicularly on disc).
Shrouded impeller type (when another

Impeller types
disc is used to cover the blades).
Blades spreading out from the shaft.
9.2.2 Classification of centrifugal pump.

Centrifugal pump may be classified in several ways as follows,
Low speed
Basis of speed Medium speed
High speed
Radial flow
Direction of flow Mixed flow
Axial flow
Low head (10m and below)
Basis of head Medium head (10m – 50m)
High head (above 50m)

Single entry
Basis of entry
Double entry
Fig. (9.6) shows the single and double entry of centrifugal pump.
Figure (9.6) (a) Single entry, (b) Double entry pumps
9.3 Pump Head.

The main assumption in this analysis the flow is to be steady, the pump
basically increases the Bernoulli head of the flow between the eye, point (1),
and the exit, point (2) as show in Fig. (9.7)
Figure (9.7) schematic of a typical centrifugal pump.

The change in head is denoted by (H).
𝑝 𝑉2 𝑝 𝑉2
𝐻 = (𝜌𝑔 + 2𝑔 + 𝑍) − (𝜌𝑔 + 2𝑔 + 𝑍) = ℎ𝑠 − ℎ𝑓 (9.1)
2 1
Where hs is the pump head supplied and hf is the losses. Now, the primary
output parameter for any turbomachine is the net head (H).
Usually V2 V1 and (Z2 - Z1) less than one meter, then, the net pump head is
𝑝 −𝑝 ∆𝑝
𝐻 = 2𝜌𝑔 1 = 𝜌𝑔 (9.2)
The power delivered to the fluid simply equals the specific weight times the
discharge times the net head change.
𝑃𝑤 = 𝜌𝑔𝑄𝐻 (9.3)
This is called the water horsepower. The power required to drive the pump is
the brake horsepower
Pbhp =  Tsh (9.4)
Where  is the shaft angular velocity.
Tsh is the shaft torque.
9.4 Pump Theory.

To predict actually the head, power, efficiency and flow rate of a pump
two theoretical approaches are possible
Simple one-dimensional formula.
Pump Theory
Complex digital – computer models.
One-dimensional flow.
To construct an elementary
theory of pump performance Idealized fluid – velocity
in this chapter, we assume vector through the impeller.
the following
Angular momentum
theorem for C.V.
To analyze the centrifugal pump, we choose the annular region that

encloses the impeller section as the control volume C.V. as in Fig. (9.8), which
shows the idealized velocity diagram and the impeller geometry.
Figure (9.8) Inlet and exit velocity diagrams for

an idealized pump impeller.
Since, r1 inside radius of impeller.
r2 outside radius of impeller.
 The angular velocity of shaft impeller blades.
At inlet,
u1 = r1 is the circumferential speed of the tip impeller at r1.
w1 velocity component of fluid tangent or parallel to the blade angle 1.
V1 is the absolute velocity of fluid at entrance, is the vector sum of
w1 and u1.
At exit,
u2 = r2 is the circumferential speed of the tip impeller at r2.
w2 velocity component outlet of fluid tangent or parallel to the blade angle
2.
V2 is the absolute velocity of fluid at exit, is the vector sum of w2 and u2.
The angular – momentum theorem to a radial flow devices (sec. 4.7),
we arrived to a result for the applied torque Tsh
𝑇𝑠ℎ = 𝑚̇ (𝑟2 𝑉𝑡2 − 𝑟1 𝑉𝑡1 ) = 𝜌𝑄(𝑟2 𝑉𝑡2 − 𝑟1 𝑉𝑡1 ) (9.5)
Where Vt1 and Vt2 are the absolute circumferential or tangent velocity
components of the flow at inlet and exit. The power delivered to the fluid is
thus,
𝑃𝑤 = 𝜔𝑇𝑠ℎ = 𝜌𝑄(𝑢2 𝑉𝑡2 − 𝑢1 𝑉𝑡1 ) (9.6)
𝑃𝑤 1
Or 𝐻 = 𝜌𝑔𝑄 = 𝑔 (𝑢2 𝑉𝑡2 − 𝑢1 𝑉𝑡1 ) (9.7)
Eq's (9.5, 9.6 and 9.7) are the Euler turbomachine equations, showing that the
torque, power and ideal head are functions only of the rotor-tip velocities u1,
u2 and the absolute fluid tangential velocities Vt1 and Vt2 independent of the
axial velocities.
We can rewriting these relation in other form, from the geometry of
Fig. (9.6)
𝑉 2 = 𝑢2 + 𝑤 2 − 2𝑢𝑤𝑐𝑜𝑠𝛽 ; 𝑤 𝑐𝑜𝑠𝛽 = 𝑢 − 𝑉𝑡
1
Or 𝑢𝑉𝑡 = 2 (𝑉 2 + 𝑢2 − 𝑤 2 ) (9.8)
Substituting Eq. (9.8) into Eq. (9.7) gives
1
𝐻 = 2𝑔 [(𝑉22 − 𝑉12 ) + (𝑢22 − 𝑢12 ) − (𝑤22 − 𝑤12 )] (9.9)
Substituting for (H) from its definition in Eq. (9.1) and rearranging, we obtain
the classic relation
𝑝 𝑤2 𝑟 2 𝜔2
+ 𝑍 + 2𝑔 − = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡 (9.10)
𝜌𝑔 2𝑔
This is the Bernoulli equation in rotating coordinates and applies to either
two or three-dimensional ideal incompressible flow. For a centrifugal pump,
the power can be related to the radial velocity Vn=Vt tan and the continuity
equation
𝑃𝑤 = 𝜌𝑄(𝑢2 𝑉𝑛2 𝑐𝑜𝑡𝛼2 − 𝑢1 𝑉𝑛1 𝑐𝑜𝑡𝛼1 ) (9.11)
Where
𝑄 𝑄
𝑉𝑛2 = and 𝑉𝑛1 = (9.12)
2𝜋𝑟2 𝑏2 2𝜋𝑟1 𝑏1
Where b1 and b2 are the blade widths at inlet and exit, with the pump
parameters (r1, r2, 1, 2 and ) are known. Eq. (9.7) or Eq. (9.11) is used to
compute idealized power and head versus discharge.
The 'design' flow rate Q* is commonly estimated by using that the flow enters
exactly normal to the impeller; 1=90 ; Vn1=V1
Ex.1
Given are the following data for a commercial centrifugal water pump:
r1=4 in, r2= 7 in, 1=30, 2 =20, speed N=1440 r/min. Estimate
a- The design – point discharge.
b-The water horsepower.
c- The head if b1=b2=1.75 in.
Sol.
2𝜋𝑁 2𝜋(1440)
(a) 𝜔 = 60 = 60 = 150.8 𝑟𝑎𝑑/𝑠
u1=r1=150.8 (4/12) = 50.3 ft/s
u2=r2=150.8 (7/12) = 88.0 ft/s
For design point Vn1=V1 and 1=90 then the inlet velocity diagram as
follows
Vn1=V1
u1=50.3 ft/s
𝑉𝑛1
tan 30 = → 𝑉𝑛1 = 𝑢1 ∗ 𝑡𝑎𝑛30 = 29.0 𝑓𝑡/𝑠
𝑢1
Hence the discharge is

4 1.75 𝑓𝑡 3
𝑄 = 2𝜋𝑟1 𝑏1 𝑉𝑛1 = 2𝜋 ( ) ( ) (29.0) = 8.87
12 12 𝑠
(b) The outlet radial velocity follows from Q
𝑄 8.87
𝑉𝑛2 = = = 16.6 𝑓𝑡/𝑠
2𝜋𝑟2 𝑏2 2𝜋( 7 )(1.75)
12 12
This enable us to construct the outlet – velocity diagram as in below
Vn2=16.6 ft/s
u2=88.0 ft/s
The tangential component is

𝑉𝑡2 = 𝑢2 − 𝑉𝑛2 cot 𝛽2 = 88.0 − 16.6 𝑐𝑜𝑡20° = 42.4 𝑓𝑡/𝑠
16.6
∝2 = 𝑡𝑎𝑛−1 = 21.4°
42.4
The power is then computed as follows, with Vt1 = 0 at the design point.
𝑃𝑤 = 𝜌𝑄𝑢2 𝑉𝑡2 = (1.94 ∗ 8.87 ∗ 88.0 ∗ 42.4) = 64100 𝑓𝑡 𝑙𝑏/𝑠
64100
𝑃𝑤 = = 117 ℎ𝑝
550
(c) Finally, the head is estimated by the following
𝑃 64100
𝐻 ≈ 𝑤 = (62.4)(8.87) = 116 𝑓𝑡
𝜌𝑔𝑄
9.5 Pressure Developed by the Impeller.
Fig. (9.9) shows the general system arrangement of a centrifugal
pump.
Zs Suction level above the water level.
Zd Delivery level above the impeller outlet.
h fs, hfd Friction losses head in both suction and delivery sides.
Vs, Vd Fluid velocities in pipes for both sides.
Applying B.E. between the water level and pump suction.
𝑝𝑠 𝑠 𝑉2 𝑝𝑎
+ 𝑍𝑠 + 2𝑔 + ℎ𝑓𝑠 = (9.13)
𝛾 𝛾
𝑝𝑠 𝑝𝑎 𝑠 𝑉2
∴ = − 𝑍𝑠 − 2𝑔 − ℎ𝑓𝑠 (9.14)
𝛾 𝛾
Similarly applying B.E. theorem between the pump delivery and the
delivery at the tank,
𝑝𝑑 𝑑 𝑉2 𝑝𝑎 𝑑 𝑉2
+ 2𝑔 = + 𝑍𝑑 + 2𝑔 + ℎ𝑓𝑑
𝛾 𝛾
𝑝𝑑 𝑝𝑎
𝑜𝑟 = + 𝑍𝑑 + ℎ𝑓𝑑 (9.15)
𝛾 𝛾
Where pd is the pressure at the pump delivery, from Eq.(9.14 &9.15)
𝑝𝑑 𝑝𝑠 𝑝𝑎 𝑝𝑎 𝑉𝑠2
− =( + 𝑍𝑑 + ℎ𝑓𝑑 ) − ( − 𝑍𝑠 − − ℎ𝑓𝑠 )
𝛾 𝛾 𝛾 𝛾 2𝑔
𝑝𝑑 −𝑝𝑠 𝑉𝑠2 𝑉𝑠2
= 𝑍𝑑 +𝑍𝑠 + 2𝑔 + ℎ𝐿 = 𝐻𝑒 + 2𝑔 (9.16)
𝛾
Where He is the effective head and hL is the total friction head.
Zd
Zs
Figure (9.9) Centrifugal pump system
9.6 Manometric Head.

The official code defines the head on the pump as the difference in total
energy head at the suction and delivery flanges. This head is defined as
manometric head. The total energy at suction inlet or suction side expressed as
suction head of fluid and termed as Hs
𝑝𝑠 𝑠𝑉2
𝐻𝑠 = + 2𝑔 + 𝑍𝑠𝑑
𝛾
Where Zsd is the height of suction gauge from datum. The total energy at the
delivery side of the pump expressed as delivery head and termed as Hd
𝑝𝑑 𝑑𝑉2
𝐻𝑑 = 𝛾
+ 2𝑔 + 𝑍𝑑𝑑
Where Zdd is the height of delivery gauge from datum. Know, the difference in
total energy of fluid is defined as Hm
𝑝 𝑝 𝑉𝑑2 −𝑉𝑠2
𝐻𝑚 = 𝐻𝑑 − 𝐻𝑠 = ( 𝛾𝑑 − 𝛾𝑠 ) + ( ) + (𝑍𝑑𝑑 − 𝑍𝑠𝑑 ) (9.17)
2𝑔
Substituting Eq.(9.16) in Eq. (9.17) and rearranging will be as follows
𝑑𝑉2
𝐻𝑚 = 𝐻𝑒 + 2𝑔 + (𝑍𝑑𝑑 − 𝑍𝑠𝑑 ) (9.18)
𝑉𝑑2
As (𝑍𝑑𝑑 − 𝑍𝑠𝑑 ) is small and 2𝑔 is also small as the gauges are fixed as close as
possible.
∴ 𝐻𝑚 = 𝑠𝑡𝑎𝑡𝑖𝑐 ℎ𝑒𝑎𝑑 + 𝑎𝑙𝑙 𝑙𝑜𝑠𝑠𝑒𝑠 (9.19)
9.7 Pump Efficiency.

When the losses are exists in the pump, then (Pw ) is
actually less than (Pbhp). The goal of the pump designer is to make the
efficiency () as high as possible over as broad a range of discharge (Q) as
possible. The efficiency can be divided to three types as follows
Manometric Efficiency
Pump Efficiency Mechanical Efficiency
Volumetric Efficiency
A- Manometric efficiency.
The ideal head (H) can be expressed by Eq. (9.7)
1
𝐻𝑖𝑑𝑒𝑎𝑙 = (𝑢2 𝑉𝑡2 − 𝑢1 𝑉𝑡1 )
𝑔
From Fig. (9.10), the inlet tangent velocity is generally equal to zero due to
no guide vanes at inlet. So the inlet triangle is right angle as shown in
Fig.(9.10.a). Therefor
V1 = Vn1 and are vertical
𝑉 𝑉
𝑡𝑎𝑛𝛽1 = 1 , 𝑜𝑟 , = 𝑛1
𝑢1 𝑢1
Since, Vt1 = 0.0
𝑢 𝑉
𝐻𝑖𝑑𝑒𝑎𝑙 = 2𝑔 𝑡2 (9.20)
From the outlet triangle
2𝜋𝑁 𝜋𝑑 𝑁 Vr1
𝑢2 = 𝜔𝑟2 = 60 𝑟2 = 602 𝛽1
V1=Vn1
𝑉𝑛2
𝑡𝑎𝑛𝛽2 = 𝑢
2 −𝑉𝑡2 u1
𝑉𝑛2
𝑢2 − 𝑉𝑡2 = 𝑡𝑎𝑛𝛽 (a) Inlet Triangle
2
𝑉𝑛2
∴ 𝑉𝑡2 = 𝑢2 −
𝑡𝑎𝑛𝛽2
From Eq. (9.20) Vn2
𝑢2 𝑉
∴ 𝐻𝑖𝑑𝑒𝑎𝑙 = 𝑛2
(𝑢2 − 𝑡𝑎𝑛𝛽 ) 2
𝑔 2
Now, manometric efficiency is defined (b) Outlet Triangle

as the ratio of manometic head and ideal
head
Figure (9.10) Velocity triangles for
backward curved bladed pump.
𝐻𝑚 𝐻 𝑔 𝐻𝑚 𝑔
𝜂𝑚 = 𝐻 = 𝑢 𝑚𝑉 = 𝑉𝑛2 (9.21)
𝑖𝑑𝑒𝑎𝑙 2 𝑡2 𝑢2 (𝑢2 − )
𝑡𝑎𝑛𝛽2
B- Mechanical Efficiency.
The mechanical efficiency is defined as
𝐸𝑛𝑒𝑟𝑔𝑦 𝑡𝑟𝑎𝑛𝑠𝑓𝑒𝑟 𝑡𝑜 𝑡ℎ𝑒 𝑓𝑙𝑢𝑖𝑑
𝜂𝑚𝑒𝑐ℎ = 𝑊𝑜𝑟𝑘 𝑖𝑛𝑝𝑢𝑡
𝜌𝑄(𝑢 𝑉 ) 𝜌𝑄(𝑢2 𝑉𝑡2 )

𝜂𝑚𝑒𝑐ℎ = 𝑃𝑜𝑤𝑒𝑟2𝑖𝑛𝑝𝑢𝑡
𝑡2
= (9.22)
𝑃𝑏ℎ𝑝
C- Volumetric Efficiency.
There are always some leakage after being
imparted energy by the impeller.
Volum deliverd (Q)
Volumetric efficiency = Volume passing throught impeller (Q+Q )
L
Where (QL) is the losses of fluid due to leakage, thus, the total efficiency
is simply the product of its three parts.
𝜂𝑜 = 𝜂𝑚 𝜂𝑚𝑒𝑐ℎ 𝜂𝑣𝑜𝑙. Or from basic definition can be defined as
𝑃 𝜌𝑔𝑄𝐻
𝜂𝑜 = 𝑃 𝑤 = 𝑃 (9.23)
𝑏ℎ𝑝 𝑏ℎ𝑝
Ex.2
The following details refer to a centrifugal pump. Outer diameter 30cm,
eye diameter 15 cm. blade angle at inlet (1=30) blade angle at outlet
(2=25). Impeller speed (N=1450 rpm). The velocity remains constant. The
whirl (tangent) velocity at inlet is zero. Determine (a)- the torque applied if the
(m=0.76%). (b)- The power when the width of blades at outlet is 2 cm and
equal the width at inlet. (c)- The head.
Sol.
2𝜋𝑟 𝑁 𝜋𝑑 𝑁 𝜋∗0.3∗1450
𝑢2 = 602 = 602 = = 22.78 𝑚/𝑠
60
2𝜋𝑟 𝑁 𝜋𝑑 𝑁 𝜋∗0.15∗1450
𝑢1 = 601 = 601 = = 11.39 𝑚/𝑠
60
V1=Vn1 since Vt1=0
From inlet velocity diagram
𝑉
tan 30 = 𝑢𝑛1 → 𝑉𝑛1 = 𝑢1 ∗ 𝑡𝑎𝑛30 = 11.39 ∗ 𝑡𝑎𝑛30 = 6.58 𝑚/𝑠
1
Vt2 = u2-x
𝑉𝑛2 V1=Vn1
𝑡𝑎𝑛𝛽2 = 1=30
𝑋
𝑉𝑛2 u1
∴𝑋=
𝑡𝑎𝑛𝛽2
𝑉𝑛2 6.58
∴ 𝑉𝑡2 = 𝑢2 − = 22.78 − = 8.67 𝑚/𝑠 Vn2=6.58 m/s
𝑡𝑎𝑛𝛽2 tan 25 x
Power delivered to the fluid as
u2=22.78 m/s
𝑃𝑤 = 𝜔𝑇𝑠ℎ = 𝜌𝑄(𝑢2 𝑉𝑡2 − 𝑢1 𝑉𝑡1 )
𝑄 = 2𝜋 ∗ 𝑟2 ∗ 𝑏 ∗ 𝑉𝑛2 = 𝜋 ∗ 𝑑2 ∗ 𝑏 ∗ 𝑉𝑛2 = 𝜋 ∗ 0.3 ∗ 0.02 ∗ 6.58
= 0.124 𝑚3 /𝑠
(a).
𝑇𝑠ℎ = 𝑚̇ (𝑟2 𝑉𝑡2 − 𝑟1 𝑉𝑡1 ) = 𝜌𝑄(𝑟2 𝑉𝑡2 − 𝑟1 𝑉𝑡1 )
𝑇𝑠ℎ = 𝜌 ∗ 𝜋 ∗ 𝑑2 ∗ 𝑏 ∗ 𝑉𝑛2 (𝑟2 𝑉𝑡2 − 𝑟1 𝑉𝑡1 )
= 1000 ∗ 𝜋 ∗ 0.3 ∗ 0.02 ∗ 6.58 ∗ (0.15 ∗ 8.67 − 0.0)
= 161 𝑁. 𝑚
(b)
𝑃𝑤 = 1000 ∗ 0.124 ∗ (22.78 ∗ 8.67 − 11.39 ∗ 0.0) = 24490 𝑊
𝑃𝑤 24490
𝑃𝑤 = = ≈ 33 ℎ𝑝
746 746
𝑃𝑤 𝑃𝑤 33
𝜂𝑚 = → 𝑃𝑏ℎ𝑝 = = ≈ 44 ℎ𝑝
𝑃𝑏ℎ𝑝 𝜂𝑚 0.76
(c)
Head (H) as
𝑃𝑤 1 1
𝐻 = 𝜌𝑔𝑄 = 𝑔 (𝑢2 𝑉𝑡2 − 𝑢1 𝑉𝑡1 ) = 9.81 (22.78 ∗ 8.67 − 0.0) = 20 𝑚
9.8 Pump Performance Curves.

Performance charts are always plotted for constant shaft-rotation speed
N (rpm). The basic independent variable is taken to be discharge Q in (gal /min
or m3/s). The dependent variables or output, are taken to be head H for liquid
(pressure rise p for gases), brake horsepower 𝑃𝑝ℎ𝑝 and efficiency. Fig.(9.11)
shows typical performance curves for centrifugal pump.
Figure (9.11) Centrifugal pump characteristics at constant speed.
The efficiency  is always zero at no flow and Qmax, and it reaches a maximum
perhaps 80 to 90 percent at about 0.6Qmax. This is the design point flow rate
Q* or best efficiency point (BEP) and max. . The head and horsepower at BEP
will be termed H* and P*bhp.
9.9 Net Positive – Suction Head (NPSH).

NPSH, which is the head required at the pump inlet to keep the liquid
from cavitation or boiling. The NPSH is defined as
𝑝 𝑉2 𝑝
𝑁𝑃𝑆𝐻 = 𝜌𝑔𝑠 + 2𝑔
𝑠
− 𝜌𝑔𝑣 (9.24)
Where (ps &Vs) are the pressure and velocity at the pump inlet and pv is the
vapor pressure of the liquid. The right–hand side is equal or greater than
(NPSH) in the actual system to avoid cavitation. Applying the energy equation
between sump surface and the pump suction level.
𝑝𝑠 𝑉2
𝑠 𝑝𝑎
+ 2𝑔 + 𝑍𝑠 + ℎ𝑓𝑠 = (9.25)
𝛾 𝛾
From Eq's (9.24&9.25)
𝑝 𝑝
𝑁𝑃𝑆𝐻 = 𝛾𝑎 − 𝛾𝑣 − 𝑍𝑠 − ℎ𝑓𝑠 (9.26)
Thoma cavitation parameter is defined by
𝑝 𝑝
(𝑁𝑃𝑆𝐻) ( 𝑎⁄𝛾)−( 𝑣⁄𝛾)−𝑍𝑠 −ℎ𝑓𝑠
𝜎= 𝐻 = 𝐻
At cavitation conditions
𝑝 𝑝
𝜎 = 𝜎𝑐 𝑎𝑛𝑑 𝛾𝑠 = 𝛾𝑣
𝑝 𝑝
( 𝑎⁄𝛾)−( 𝑣⁄𝛾 )−𝑍𝑠 −ℎ𝑓𝑠 𝑉2
𝑠
𝜎𝑐 = = 2𝑔𝐻 (9.27)
𝐻
The height of suction, the frictional losses in the suction line play an important
role for avoiding cavitation at a location.
9.10 Outlet Blade Angle and Specific Speed.

Different blade arrangements can be used in design of centrifugal
pumps. Three possible orientations of outlet blade angle. Forward curved blade
(β2 >90º), radial curved blade (β2=90º) and backward curved blade (β2 < 90º).
Fig. (9.12) are shown the velocity triangles for three arrangements.
Vt2 u2=Vt2 u2
u2 Vt2
Vr2
V2 Vn2 V2 Vn2 Vn2
V2 Vr2
Figure (9.12) Blade shape arrangements with outlet velocity triangles
The dimensionless parameter (Ns) is known as the specific speed of pumps. In

practice is used as
𝑁√𝑄
𝑁𝑠 = 3 (9.28)
𝐻 ⁄4
The (Ns) quantity have been derived from dimensional analysis technique, and
defined as the speed at which the pump will operate to deliver unit flow under
unit head.
Ex.3
The outer diameter and width of a centrifugal pump impeller are 55 cm and 3
cm. The pump runs at 1300 rpm. The suction head is 7 m and the delivery head
is 45 m. The frictional drops in suction side is 2.5 m and in the delivery side is
9 m. The blade angle at out let is 33°. The manometric efficiency is 83% and
the overall efficiency is 77%. Determine
a- The power required to drive the pump.
b- Calculate the pressures at the suction and delivery side of the pump.
Sol.
Assume the inlet tangential velocity is zero, the total head against the pump is
H = 45+7+2.5+9= 63.5 m
𝜋𝑑 𝑁 𝜋∗0.55∗1300 𝑚
𝑢2 = 602 = 60
= 37.4 𝑠 ; 𝑠𝑖𝑛𝑐𝑒 𝑢2 = 𝜔𝑟2
𝑔𝐻 9.81∗63.5
𝜂𝑚 = 𝑢 = 0.83 = 𝑠𝑜𝑙𝑣𝑖𝑛𝑔 𝑓𝑜𝑟 𝑉𝑡2 = 20.0 𝑚/𝑠
2 𝑉𝑡2 37.4∗𝑉𝑡2
The outlet velocity triangle is used to calculate Vn2,
𝑉𝑛2
𝑡𝑎𝑛𝛽2 = 𝑢 −𝑉
2 𝑡2
∴ 𝑉𝑛2 = 𝑡𝑎𝑛𝛽2 (𝑢2 − 𝑉𝑡2 ) = 𝑡𝑎𝑛33(37.4 − 20.0) = 11.36 𝑚/𝑠
Flow rate Q is
𝑄 = 𝜋 ∗ 𝑑2 ∗ 𝑏2 ∗ 𝑉𝑛2 = 𝜋 ∗ 0.55 ∗ 0.03 ∗ 11.36 = 0.58886 𝑚3 /𝑠
𝑃𝑤 𝑃𝑤 𝜌𝑔𝑄𝐻
𝜂𝑜 = 𝑆𝑜𝑙𝑣𝑖𝑛𝑔 𝑓𝑜𝑟 𝑃𝑜𝑤𝑒𝑟 (𝑃𝑏ℎ𝑝 ) = =
𝑃𝑏ℎ𝑝 𝜂𝑜 𝜂𝑜
1000 ∗ 9.81 ∗ 0.58886 ∗ 63.5
∴ 𝑃𝑏ℎ𝑝 = = 476391.56 𝑊
0.77
Now, consider the water level and the suction level are (1 & 2)
𝑝𝑎 𝑉1 2 𝑝2 𝑉2 2
+ + 𝑍1 = + + 𝑍2 + 𝑙𝑜𝑠𝑠𝑒𝑠
𝛾 2𝑔 𝛾 2𝑔
101000 𝑝2 11.362 𝑝
+0+0 = + + 7 + 2.5 𝑠𝑜𝑙𝑣𝑖𝑛𝑔 𝑓𝑜𝑟 2⁄𝛾
9810 9810 2 ∗ 9.81
𝑝2
⁄𝛾 = 1.218 𝑚 𝑎𝑏𝑠𝑜𝑙𝑢𝑡𝑒
Consider suction side and delivery side as 2 and 3,
𝑝2 𝑉2 2 𝑝3 𝑉3 2 𝑢2 𝑉𝑡2
+ + 𝐻𝑖𝑑𝑒𝑎𝑙 = + + 𝑍3 ; 𝑠𝑖𝑛𝑐𝑒 𝐻𝑖𝑑𝑒𝑎𝑙 =
𝛾 2𝑔 𝛾 2𝑔 𝑔
2 2 2
𝑉3 = √𝑉𝑡2 + 𝑉𝑛2
= √20 + = 23 𝑚/𝑠 11.362
11.362 37.4 ∗ 20 𝑝3 232
1.218 + + = + ; 𝑠𝑜𝑙𝑣𝑖𝑛𝑔 𝑓𝑜𝑟
2 ∗ 9.81 9.81 𝛾 2 ∗ 9.81
𝑝3
= 57 𝑚 𝑎𝑏𝑠𝑜𝑙𝑢𝑡𝑒.
𝛾
Ex.4
A centrifugal pump was tested for cavitation initiation. Total head was
40 m and flow rate was 0.06 m3/s. Cavitation started when the total head at the
suction side was 3 m. The atmospheric pressure was 760 mm Hg and the vapor
pressure at this temperature was 2 kPa. It was proposed to install the pump
where the atmospheric pressure is 700 mm Hg and the vapor pressure at the
location temperature is 1 kPa. If the pump develops the same total head and
flow, can the pump be fixed as the same height as the lab setup? What should
be the new height?
Sol.
At the suction point
Total head = Vapour pressure + velocity head.
∴ Velocity head = Total head – Vapor pressure in head of water
𝑉𝑠2 2∗103
= 3 − 103 ∗9.81 = 2.796 𝑚
2𝑔
Cavitation parameter (σ ) is defined by
𝑠 𝑉2 2.796
𝜎 = 2∗𝑔∗𝐻 = = 0.0699
40
Applying B.E. between water level and suction point,
𝑝1 𝑉2
1 𝑝𝑎𝑡𝑚
+ 2𝑔 + 𝑍1 = − ℎ𝑓
𝛾 𝛾
𝑝𝑎𝑡𝑚 𝑝𝑣 760∗13.6 2∗103
(𝑍1 + ℎ𝑓 ) = − 𝜎𝐻 − =( ) − 2.796 − 103 ∗9.81 = 7.336 𝑚
𝛾 𝛾 1000
At the new location (head and flow rate being the same, friction loss will be
the same)
𝑝𝑎𝑡𝑚 𝑝𝑣 700∗13.6 1∗103
(𝑍1 ′ + ℎ𝑓 ′ ) = − 𝜎𝐻 − =( ) − 2.796 − 103 ∗9.81 = 6.622 𝑚
𝛾 𝛾 1000
ℎ𝑓 = ℎ𝑓′
∴ (𝑍1 − 𝑍1′ ) = 0.714 𝑚
The pump should be lowered by 0.714 m, since the new height is 6.622 m.
Problems.
P9.1 Determine the discharge flow rate from centrifugal pump at 1000
rev./min, the head being 14.5 m. the vane angle at outlet is 30º to the
periphery. The impeller diameter is 0.3 m and width is 0.05 m. The
manometric efficiency of the pump is 0.85%. [ 1.14 m3/s]
P9.2 A centrifugal pump impeller is 0.5 min diameter and delivers 2 m3/min
of water. The peripheral velocity is 10 m/s and the flow velocity is 2 m/s.
the blade outlet angle is 35º. Whirl at inlet is zero. Determine the power
and torque delivered by the impeller. [2.18 kW, 54.5 N m]
P9.3 A centrifugal pump with 2.3 m diameter impeller running at 327 rev./min.
delivers 7.9 m3/s of water. The head developed is 72.8 m. The width of
the impeller at outlet is 0.22 m. if the overall efficiency is 91.7% determine
the power to drive the pump. Also determine the blade angle at exit.
[61.52 kW, 13º]
P9.4 A centrifugal pump running at 1000 rev./min. works against a head of 80

m delivering 1 m3/s. The impeller diameter and width are 80 cm and 8
cm respectively. Leakage loss 3 percent of discharge. Hydraulic efficiency
is 80%. External mechanical loss is 10kW. Calculate the blade angle at
outlet, the power required and the overall efficiency.
[15.5º, 1020kW, 76.9%]
P9.5 The diameters of a centrifugal pump impeller is 600 mm and that of the
eye is 300 mm. The vane angle at inlet is 30° and that at outlet is 45°. If
the absolute velocity of water at inlet is 2.5 m/s determine the speed and
manometric head. The whirl at inlet is zero. [275.8 rpm, 5.44 m]
P9.6 The speed of a centrifugal pump was 240 rpm and it is required to develop
22.5 m head when discharging 2 m3/s of water. The flow velocity at outlet
is 2.5 m/s. The vanes at outlet are set back at 30° to the tangential
direction. Determine the manometric efficiency and the power required to
drive the pump. Impeller diameter is 1.5 m. [81%, 545 kW]
P9.7 A centrifugal pump delivers 50 l/s when running at 1500 rpm. The inner
and outer diameters are 0.15 m and 0.25 m respectively. The blades are
curved at 30° to the tangent at the outlet. The flow velocity is 2.5 m/s and
is constant. The suction and delivery pipe diameters are 15 cm and 10 cm,
respectively. The pressure head at suction is 4 m below atmosphere. The
pressure at the delivery is 18 m above atmosphere. The power required
was 18 kW. Determine the vane angle at inlet for zero whirl at inlet. Also
find the manometric efficiency and overall efficiency.
[12°, 77.2%, 64.6%]
P9.8 A centrifugal pump running at 1450 rpm delivers 0.11 m3/s of water
against a head of 23 m. The impeller diameter is 250 mm and the width is
50 mm. The manometric efficiency is 75% determine vane angle at outlet.
[42.1°]
P9.9 A six stages centrifugal pump running at 660 rpm is to deliver 0.8 m3/s
of water of water against a manometric head of 88 m. The vanes are curved
back at 38° to the tangent at outer periphery. The velocity of flow is 35% of
the peripheral velocity at outlet. The hydraulic losses are 33% of the velocity
head at the outlet of the impeller. Determine the diameter of the impeller at
outlet and the manometric efficiency. [0.5 m, 87.4%]
P9.10 A multi stage pump is required to deliver of 2.5 l/s water against a
maximum discharge head of 250 m. The diameter of radial bladed impeller
should not be more than 18 cm. Assume a speed of 3000 rpm. Determine the
impeller diameter, number of stages and power. [160 mm, 4, 8.76 kW]

جريان موائع 2

Uploaded by

Copyright:

Available Formats

جريان موائع 2

Uploaded by

Document Information

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

جريان موائع 2

Uploaded by

Copyright:

Available Formats

Chapter 6 Dynamic Similarity and Dimensional Analysis 129

Complete similarity Dynamic Similarity

6.2 Geometric Similarity (G.S).

Therefore, geometric similarity is a necessary condition for the kinematic

6.4 Dynamic Similarity (D.S).

According to Newton's law, the resultant 𝐹𝑅 of all these forces, will

Where Δp is some characteristic pressure in the flow.

e- Capillary or Surface Tension Force.

f- Compressibility or Elastic Force.

6.4.1 D. S. of Flow Governed by Viscous, Pressure and Inertia Forces.

𝑉𝑖𝑠𝑐𝑜𝑢𝑠𝑓𝑜𝑟𝑐𝑒 |𝐹⃗ | 𝜇𝑉𝐿 𝜇

𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒𝑓𝑜𝑟𝑐𝑒 |𝐹⃗𝑝 | ∆𝑝𝐿2 ∆𝑝

The term LV/ is known as Reynolds number, Re.

1- The full flow of fluid in a completely closed conduit

Hence, Re & Eu for a complete dynamic similarity between prototype and

6.4.2 D.S. of Flow Governed by Gravity and Inertia Forces.

6.4.3 D.S. of Flows with Surface Tension as the Dominant Force.

6.4.4 D.S. of Flows with Elastic Force.

Dimensionless Representation Name Recommended

LV/ Inertiaforce Reynolds Re

∆p/V 2 pressureforce Euler number Eu

V/(Lg)1/2 Inertia force Froud number Fr

V/√Es / Inertiaforce Mach number M𝑎

The forces control the flow is

By equating the two ratio

6.5 The Application of D.S and the Dimensional Analysis.

6.5.2 Dimensional Analysis.

The dimension of velocity gradient du/dy can be written as

du/dy = dimension of u/ dimension of y= 𝑇

6.6 Rayleigh's Indicial Method (Method-1).

6.7 Buckingham's Pi Theorem. (Method-2)

6.7.1 Mathematical Description of (𝝅 ) Pi Theorem.

6.7.2 Procedure for Determination 𝝅 Terms.

6.7.3 Selection of Repeating Variables (R.V).

Pressure drop/length M L−2 T −2

Flow rate (Q) L3 T −1

Weight(W) per unit volume 𝐹 𝐿−3

Step-2:- there are four variables & (3) fundamental dimensions

P6.10 The resisting force F at a plane during flight can be considered as

P6.11 The pressure difference Δp in a pipe of diameter D and length L due to

P6.14 A spillway model is to be built to a geometrically similar scale of 1/50

P6.18 In a geometrically similar model of spillway the discharge per length is

7.2 Relationship between Shear Stress and Pressure Gradient.

Let  = shear stress on the lower face ABCD of the element

7.3 Laminar Flow between Parallel Plates.

If L is the length between the sections

The governing Equation is

7.4.1 Maximum and Minimum Velocities.

7.5 Pipe of Circular Cross-Section.

Figure 7.4: Flow in circular pipe.

7.5.3 Shear Stress in Horizontal Pipe.

Figure 7.5: Forces on element in horizontal pipe.

7.5.4 Shear Stress in Inclined Pipe.

By equating Eq's 7.30 & 7.37

(i) The rate of flow of oil.

Table 7.1: Summary of used equations in pipe

7.7 Simple Pipe Problem.

Problem Given To find (unknown)

1% yield diff-with Darcy equation

7.7.1 Solution Procedures.