MyBook 11
MyBook 11
MyBook 11
Dear Students
I am Dr. Fuad Musleh, I will be teaching you Hydraulics,
CENG321. It is one of the most interesting subjects in Engineering.
You will notice that the very important theories and concepts,
governing fluid motion in pipes, and open channels will be taught
in a practical manner. Which gives you at the end of the course
tools that you need in order to solve the water flow problems.
These scientific notes were collected from different sources:
Several Hydraulics references; notes by Professor Shamsul
Haque Alvi, one of the founders of the civil Engineering
department of University of Bahrain; and finally my own notes
which is supported by continuous additions and updated work. I
believe it will be enough for you to understand the course.
The office hours will be announced soon. However, you are
free to contact me through Microsoft teams and Blackboard at
any time. I will do my best to give clear answers to all your
questions.
Do your best to solve the home works on time, your efforts
will be always appreciated. Using other students work is
completely prohibited, you both might lose the whole mark.
Dr. Fuad Musleh
18/9/2020
1
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
2
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
(2): using the Moody diagram, given a horizontal pipe of 0.1 m diameter, a pressure drop of 0.2 m/m,
a kinematic viscosity of 1×10-6 m2/s, find the flow discharge in the pipe.
……………..………………………….………………………………………………..…………………………………………………………………
……………………………………………………………………………………………..…………………………..……………….
(3): Find the location of the surface in
the manometer.
(5): For the system of pipes in series shown, given that: The flow of water in the system is 1500 L/s,
D1=1.0m, D2 =0.4m, and D3=0.7m.
a. The velocities (m/s) in pipes 1, 2 and 3 are: V1=..............m/s, V2=............. m/s, and V3=............. m/s
b. Doubling the head into the tank will (increase, decrease, not change) the discharge within the pipes.
c. Doubling the head into the tank will (increase, decrease, not change) the velocities within the pipes.
(6): For the system of pipes in series shown, given that: The flow of water in the system is 1000 L/s,
kexit =1, Kentrance =0.5, the Losses due expansion is (V1-V2)2/2g, the kinematic viscosity is =3×10-6 m2/s.
3
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
a. The velocities (m/s) in pipes 1 and 2 are: V1=..............m/s, and V2=............. m/s
b. The Reynolds no. values in pipes 1 and 2 are.............and.............
c. The friction factors in pipes 1 and 2 are.............and.............
d. The difference in head between the two tanks is H =.............m
4
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Definitions:
What is a Water Engineer?
It is also known as: Water Resource Engineer, Wastewater Engineer, Water Resources
Engineer, Water and Sewer Engineer, Water and Waste Engineer.
It is a specialized type of Engineer. A water engineer is someone who deals with the
provision of clean water, disposal of waste water and sewage, and the prevention of
flood damage. He demonstrates a genuine knowledge and interest of the water
industry and environmental issues.
Generally, water resource engineer works on designing pump systems, pipelines and
treatment plants for handling water, wastewater, or black water as it is known in the
industry. Water resource engineers also work on flood control projects and help map
out drainage systems for dealing with water flow. They may also work on fresh water
supply systems, such as those used for drinking water.
In more details, his job involves:
Design water management systems, sewer improvement plans or flood
defence programs, and associated structures, such as pumping stations,
pipework and earthworks and oversee the construction and maintenance of
these systems later.
Repairing, maintaining and building structures that control water resources like
sea defence walls, pumping stations and reservoirs.
He ensure that citizens are provided with a continuous supply of clean,
uncontaminated water for drinking, living, and recreational purposes.
preparing tender documents for construction
reviewing technical submissions
keeping up to date with environmental matters; being aware of policy and
developments
Supervising the operation and maintenance of water and sewage
infrastructure.
Using computer simulations to analyze, for example, potential dam failure.
devising flood defense strategies
monitoring flood levels at times of high risk
managing staff, including other engineers, technicians and site workers.
Fluids (Liquids and gases): a substance which deforms continuously, or flows, when
subjected to shear stress [The fluid particles easily move and change their relative position
and continuously deform (flow) under the action of shear stress even with small forces].
5
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Water resources engineering: is the quantitative study of the hydrologic cycle (the
distribution and circulation of water linking the earth's atmosphere, land and oceans).
Surface runoff is an important factor of water resources, it is measured as the
difference between precipitation and abstractions, such as infiltration (which
replenishes groundwater flow), surface storage and evaporation.
Water resources applications: include the management of the urban water supply,
the design of urban storm-sewer systems, and flood forecasting.
Hydraulic engineering: It is as a sub-discipline of civil engineering concerned with the flow
and conveyance of fluids, principally water and sewage. It consists of the application of fluid
mechanics to water flowing in an isolated environment (pipe, pump) or in an open channel
(river, stream).
Civil engineers are concerned with open channel flow, which is governed by the
interdependent interaction between the water and the channel.
One feature of hydraulic systems in civil engineering is the extensive use of gravity
as the motive force to cause the movement of the fluids. This area of civil
engineering is intimately related to the design of bridges, dams, channels, canals,
and levees, and to both sanitary and environmental engineering.
Hydraulic engineering includes the application of the principles of fluid mechanics
to problems dealing with the collection, storage, control, transport, regulation,
measurement, and use of water.
The hydraulic engineer is concerned with the transport of sediment by the river, the
interaction of the water with its alluvial boundary, and the occurrence of scour and
deposition.
The hydraulic engineer develops conceptual designs for the various features which
interact with water such as spillways and outlet works for dams, culverts for
highways, canals and related structures for irrigation projects, and cooling-water
facilities for thermal power plants.
Shear force: is the force component tangent to the surface. Average shear stress is the shear
force per unit area.
Fluid mechanics: the study of liquids and gases at rest (static) and in motion (dynamics) [the
science that deals with the action of forces on fluids].
Engineering applications: dams and reservoirs; water supply pipelines; groundwater
movement; runoff in parking lots; and pumps, filters, rivers, etc.
Fluids can be classified as ideal fluids and real fluids.
Ideal fluids: are those which are incompressible with zero viscosity and, shear stress is always
zero. Ideal fluid is hypothetical. The main difference between an ideal fluid and a real fluid is
that for ideal flow pressure is constant and for real flow pressure is decreasing along the
flowing fluid.
6
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Real fluids: Fluids with viscosity. Example: Water, Milk, and Honey etc., real fluids are
classified as:
Newtonian and non-Newtonian: Examples of non-Newtonian fluids are Thixotrophic
substance (thixotrophic jelly paints), ideal plastic, Bingham plastic (sewage sludge),
pseudo plastic (clay, milk, cement), dilatant substance (quick sand) etc.
Hydraulic structure: is a structure submerged or partially submerged in any body of water,
which disrupts the natural flow of water.
They can be used to divert, disrupt or completely stop the flow. An example of a
hydraulic structure would be a dam, which slows the normal flow rate of the river in
order to power turbines. A hydraulic structure can be built in rivers, a sea, or any
body of water where there is a need for a change in the natural flow of water.
Hydraulic structures may also be used to measure the flow of water. When used to
measure the flow of water, hydraulic structures are defined as a class of specially
shaped, static devices over or through which water is directed in such a way that
under free-flow conditions at a specified location (point of measurement) a known
level to flow relationship exists. Hydraulic structures of this type can generally be
divided into two categories: flumes and weirs.
7
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Chapter 1: Revision
1.1 Fundamental properties of water:
Mass density: mass of the substance per unit volume (kg/m3).
The density of some fluids is easily changed with a pressure change. Air can easily be
compressed with consequent density change, whereas a very large pressure is required to
effect relatively small density change in water.
Water can be considered incompressible and can be assumed to have constant
density with constant pressure.
Air is a compressible fluid with a variable density as a function of pressure.
In general, densities of liquids are almost constant under constant temperature even
with huge changes in pressure.
Sea water contains dissolved salts (more specific weight) sea water is 4% freshwater.
Water expands when freezes its density decreases and rises up in water.
more volume is needed for water in pipes frozen.
Maximum density of water occurs at 4 oC less dense when heated.
Temperature (oC) Density (water , kg/m3) Specific Weight (, N/m3)
0o (ice) 917 8,996
0o (WATER) 999 9,800
4o 1000 maximum 9,810
20o 998 9,790
30o 996 9,771
60o 983 9,643
Specific weight : the force exerted by the earth’s gravity up on a unit volume of
substance (the weight of the fluid per unit volume). It has the units N/m3, and can be given
as:
weight / volume (mass × g)/Volume ×g
water (@ 4 oC, 1 atm) =9810 N/m3.
water (@ 20 oC, 1 atm) =9790 N/m3.
8
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Specific gravity: The ratio of the specific weight of any liquid at a given temperature to that
of water at 4 oC.
Example1.1: An aquarium holds 0.5 m3 of water. The weight of the aquarium is 5,090 N
when full and 200 N when empty. Determine the temperature of the water.
Solution:
The weight of the water in the aquarium is: Wwater=5,090 N-200 N = 4,890 N
The specific weight of the water is: = 4,890 N/(0.5 m3) = 9,780 N/m3
The temperature of the water can be obtained from the tables to be 25 oC.
Vapour pressure can be defined as the pressure exerted by the water vapour in the
atmosphere independently of other gases.
In still water (static), the pressure will increase as we go deep in the water.
∆𝑃 = × 𝑦 (1.2)
It can be easily concluded that If one travels upward in the fluid (positive z direction),
the pressure decreases; and if one travels downward (negative z direction), the
pressure increases.
A very important equation for the variation of pressure for the case of static fluid with
uniform density can represented as:
(𝑝/ + 𝑧 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡) (1.3)
The left term of the previous equation is called piezometric
head. Z is shown on the following figure.
This equation implies that the piezometric head is constant
throughout an incompressible static fluid.
∆𝑃 = 𝛾𝑦 → 𝑝2 − 𝑝1 = 𝛾(𝑍1 − 𝑍2 ) →
𝒑𝟐 𝒑 𝒑
⁄𝜸 + 𝒁𝟐 = 𝟏⁄𝜸 + 𝒁𝟏 = 𝒄𝒐𝒏𝒔𝒕. = ⁄𝜸 + 𝒁 (1.4)
10
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Pascal's Law:
Pressure added to a confined fluid acts equally in all directions for the fluid itself.
Quiz:
Derive a formula for the variation of hydrostatic pressure between any two points in
water with a specific weight
Solution:
- Consider two arbitrary points A and B along an arbitrary
X-axis with an angle.
- Consider a control volume as shown having a cross
sectional area dA and a length L where the points A and
B lie at the ends of the control volume.
- Let the pressures at A and B be PA and PB.
- Because the prism is at rest, all forces acting on it must
be in equilibrium.
{ ∑𝐹𝑋 = 𝑃𝐴 𝑑𝐴 – 𝑃𝐵 𝑑𝐴 + 𝐿 𝑑𝐴 𝑆𝑖𝑛( ) = 0 }
Note that L sin () = h (the vertical elevation difference
between the two points).
The above equation reduces to:
𝑃𝐵 – 𝑃𝐴 = ℎ Hydrostatic pressure on a prism
11
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
“The difference in pressure between any two points in still
water is always the product of the specific weight of water
and the normal difference in elevation between the two
points.”
If A & B are on the same elevation (h=0) PA = PB.
(for static fluid, the pressure is constant at the same level
12
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
13
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example1.4:
Points 1,2,3, and 4 have equal pressures and their horizontal surface is a surface
of equal pressure.
Points 5 and 6 are on the same horizontal plane but the pressures at them is
different (not connected liquid, different distance from surface).
The horizontal surface through points 7 and 8 that passes through the interface of
the two immiscible liquids is an equal pressure surface (connected liquid between
7 and 8 from downward at the same level).
Points 9 and 10 are not on an equal pressure surface (different liquids, the distance
from them to points7 and 8 are the same with different values of , given that
pressures at 7 and 8 are the same).
Note that: immiscible liquids do not mix under normal conditions).
Quiz:
Two meters of water is equivalent to how many millimetres of mercury?
1.2.3 Manometers
A manometer is a pressure–measurement device consists usually from a
bent tube in the form of U that contains a fluid of known specific gravity.
The difference in pressure at the two ends of the manometer is found by
measuring the difference in elevations of the liquid surfaces.
There are mainly two types of manometers:
1. An open manometer: has one end open to the atmospheric pressure and is
capable of measuring the gauge pressure in a vessel (Figure 1.1).
2. A differential manometer: has each end connected to a different pressure
tap. It is used to measure the pressure difference between two taps
(Figure1.1).
14
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
The liquid used in a manometer
is usually heavier than the fluids
to be measured. It must not mix
with the liquids.
The most frequently used
manometer liquids are mercury
(sp.gr. =13.6), water, alcohol
(sp. gr. =0.9).
The higher the pressure in Figure 1.1 Types of manometers: (a) open manometer
vessel A (Fig. 1.1(b)), the larger (b) differential manometer
the difference, h.
The equations of the two manometers above can be found using the following:
1. Open manometer:
𝑃𝐴 + 𝑦 𝑤𝑎𝑡𝑒𝑟 − ℎ 𝑚𝑒𝑟𝑐𝑢𝑟𝑦 = 0
𝑃𝐴 = ℎ 𝑚𝑒𝑟𝑐𝑢𝑟𝑦 − 𝑦 𝑤𝑎𝑡𝑒𝑟 (1.5)
2. Differential manometer:
𝑃𝐴 + 𝑦 𝑤𝑎𝑡𝑒𝑟 − ℎ 𝑚𝑒𝑟𝑐𝑢𝑟𝑦 – (𝑦 − ℎ) 𝑤𝑎𝑡𝑒𝑟 = 𝑃𝐵
∆𝑃 = 𝑃𝐴 − 𝑃𝐵 = ℎ (𝑚𝑒𝑟𝑐𝑢𝑟𝑦 − 𝑤𝑎𝑡𝑒𝑟 ) (1.6)
15
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example1.6:
A mercury manometer (sp.gr. =13.6
is used to measure the pressure
difference in vessels A and B. As
shown in figure 1.2. Determine the
pressure difference in Pascal’s
(N/m2).
Solution:
Points 3 and 4 (P3 and P4) are on a
surface of equal pressure, the
same for points 1 and 2
P3 = P4 and PA = P1 = P2 .
Figure 1.2 Differential manometer
PA+w (27cm)-Hg. (15cm)-water
(135cm)=PB
∆P=PA-PB=w (135cm-27cm)+
Hg(15cm)
∆P=w [108+(13.6)(15)]cm=9790
N/m3 3.12m =30,500 pascal.
Example1.7:
Find the pressure in pipe A
Solution:
Example 1.8:
If the atmospheric pressure is 0.9 bar abs and a gage
attached to a tank reads 390 mmHg vacuum. What is
the absolute pressure within the tank?
Solution:
Patm = 0.9 bar 100 kPa/bar = 90 kPa, Pgage = [(13.6
m water/m Hg)][390 mm Hg (1m/1000mm)] =
51.9 kPa vacuum=-51.0 kPa
Pabs = 90 +(-51.9) = 38.1 kPa
16
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example1.9:
17
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
1.3 Viscosity:
“It is simply a measure of a fluid resistance to flow.”
A more specific definition, it is a measure of the fluid resistance to gradual deformation
by shear stress or tensile stress. Viscosity is a property of the fluid which opposes the relative
motion between the two surfaces of the fluid in a fluid that are moving at different velocities.
When the fluid is forced through a tube, the particles which compose the fluid generally
move more quickly near the tube's axis and more slowly near its walls (no slip condition) ;
therefore some stress (such as a pressure difference between the two ends of the tube) is
needed to overcome the friction between particle layers to keep the fluid moving. This stress
required to keep the fluid motion is proportional to its viscosity.
A fluid with a low viscosity → flows easily.
A fluid with a high viscosity → flows with difficulty.
→ requires larger force to move it.
→ more resistance to flow.
→ more power consumption due to frictional losses.
The ease with which a fluid [honey, water, oil,...] pours is an Indication of its
Viscosity.
Viscosity is important when we calculate the energy lost from a fluid as it flows in
a pipe.
A fluid that has no resistance to shear stress is known as an ideal or inviscid fluid.
Zero viscosity is observed only at very low temperatures in super fluids. Otherwise,
all fluids have positive viscosity, and are technically said to be viscous or viscid.
One of the most common instruments for measuring kinematic viscosity is the glass
capillary viscometer.
The viscosity of water and liquids increases as the temperature decreases:
→ cold oil pours more slowly than hot oil.
→ a greater force is required to stir cold oil.
In general, the effect of shear stress on a liquid is a continuous yield in angular
deformation in the direction of the shear. Fluid particles easily move and change their relative
position under the action of shear stress even with small values.
→ small values of shear stress fluid motion.
→ the magnitude of shear stress required to cause a specific flow depends
on the viscosity of the fluid.
→ shear stress, τ, is the force per unit area (N/m2).
→ for fluids (water, oil, ..): τ α change of velocity along the depth.
18
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
1.3.1 Viscosity concept:
Assume a thin layer of fluid between two surfaces, one is stationary and the other
is moving.
Note that the fluid has the same velocity as the boundary.
Figure 1.3 Velocity distribution of the flow between a stationary and a moving plate
The fluid in contact with the lower surface has zero velocity (as is the velocity of
the boundary) while that in contact with the upper surface has the velocity V.
As the distance between the two surfaces is small, then the rate of change of
velocity with Y is linear.
𝛥𝑉
Velocity gradient = (1.7)
𝛥𝑌
𝛥𝑉
= 𝜇 𝛥𝑌 (1.8)
19
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
1.3.2 Variation of fluid viscosity with temperature:
As the temperature of liquids is increased, its viscosity will decrease noticeably.
Gases behave differently from fluids in that the viscosity increases as the
temperature increases. Also the amount of change is generally smaller than that
for liquids.
An example of the variation of fluid viscosity with temperature: Engine oil is
generally difficult to pour when it is cold because its viscosity has increased.
A measure of how significant the viscosity of fluid changes with the temperature is
given by the viscosity index (VI). This is more important to lubricating oils and
hydraulic fluid in equipment that operate at wide extremes of temperatures.
{A fluid with a high viscosity index exhibits a small change in viscosity
with temperature}.
Adhesion and cohesion are water properties that affect every water molecule on
earth and also the interaction of water molecules with molecules of other
substances. Essentially, cohesion and adhesion are the "stickiness" that water
molecules have for each other and for other substances. The water drop is
composed of water molecules that like to stick together, an example of the
property of cohesion: The water drop is stuck to the end of the pine needles, which
is an example of the property of adhesion.
{Since cohesion is the predominant cause of viscosity in a liquid, this is
why viscosity decreases with temperature as cohesion does }.
For ordinary pressures viscosity is independent of pressure and depends upon
temperature only. For very great pressures, gases and most liquids shown erratic
variations of viscosity with pressure.
20
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Tutorial. 1 (Revision)
Q2. A liquid when poured into a graduated cylinder is found to weigh 6.5 N when
occupying a volume of 310 ml (milliliters), determine its specific gravity, density, and
specific weight.
Q3. Mercury has a specific weight of 133 kN/m3 when the temperature is 20°C. Determine
its density and specific gravity at this temperature.
21
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q4. If air within the tank is at an absolute pressure of 680 kPa and a temperature of 70°C,
determine the weight of the air inside the tank. The tank has an interior volume of
1.35 m3.
Q5. The tanker carries 1.5(106) barrels of crude oil in its hold. Determine the weight of
the oil if its specific gravity is 0.940. Each barrel contains 42 gallons, and there are
7.48 gal/ft3.
22
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q6. Water in the Swimming pool has a measured depth of 3.03 m, when the Temperature
is 5°C. Determine its approximate depth when the temperature becomes 35°C.
Neglect losses due to evaporation. (Hint: the mass of water does not change)
23
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q8. Oil with a specific gravity of 0.85 forms a layer 0.9 m deep in an open tank that is
filled with water. The total depth of water and oil is 3 m. what is the gage pressure at the
bottom of the tank?
24
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q10.Water, Oil, and Salt water fill a tube as shown, determine the pressure at point 1
(Inside the closed tube).
Q11. Pressure gage B is to measure the pressure at point A in a water flow. If the pressure
at B is 87 kPa, estimate the pressure at A, in kPa. Assume all fluids are at 20oC.
25
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q12. Draw the pressure distribution over the inclined and normal plates
shown.
Q13. The tank contains water and immiscible oil at 20oC. What is h is cm if the density of
the oil is 898 kg/m3?
26
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q14. Determine the elevation difference, Δh, between the water levels of the two open
tanks shown.
Q15. Find the difference in pressure between pipes A and B in the figure shown.
27
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
∑ 𝑸𝒊𝒏 = ∑ 𝑸𝒐𝒖𝒕
𝟓𝒎 𝟐
𝟏𝟐𝒎 𝟐
𝒎𝟑
× 𝟎. 𝟐𝒎 − × 𝟎. 𝟏𝟓𝒎 − 𝟎. 𝟏 + 𝑽𝟐 × 𝟎. 𝟐𝒎𝟐 = 𝟎
𝒔 𝒔 𝒔
28
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
29
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q19.
The oil derrick has drilled 5 km into
the ground before it strikes a crude oil
reservoir. When this happens, the
pressure at the well head A becomes
25 MPa. Drilling “mud” is to be placed
into the entire length of pipe to
displace the oil and balance this
pressure. What should be its density
so that the pressure at A becomes
zero?
30
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q20. The underground storage tank used in a service station contain gasoline filled to
the level A. Determine the pressure at each of the five identified points. Note that
point B is located in the stem, and point C is just below it in the tank. Take
ρgas=730 kg/m3
.•
31
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q21. The pipes connected to the closed tank are completely filled with water. If the
absolute pressure at A is 300 kPa. Determine the force acting on the inside of the
end caps at B and C if the pipe has an inner diameter of 60 mm.
32
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Figure 2.1 Energies within a moving fluid (a) Pipe flow (b) Channel flow
A fluid in motion experiences, in addition to gravity force, pressure forces, viscous and
turbulent shear resistances, boundary resistance and forces due to surface tension
and compressibility effects of the fluid. The presence of all these forces makes the
analysis very complicated.
A simplifying approach to the problem for some cases is made by assuming the fluid to
be ideal, i.e. non-viscous or frictionless and incompressible. Water has a relatively low
viscosity and is incompressible and is found to behave like an ideal fluid.
The governing equation is: The newton second law of motion.
Where ax is the total acceleration in the direction of motion, V is the volume of fluid,
Fx is the sum of all the forces in the direction of motion, and is the mass density of
the fluid.
𝐹𝑥 𝑑𝑣 𝑑𝑣
∑ = 𝜌𝑎𝑥 = 𝜌 ( + 𝑣 ), (2.2)
∀ 𝑑𝑡 𝑑𝑥
𝑑𝑣
{Algebraic sum of all forces per unit volume}, where, , is the local acceleration and
𝑑𝑡
𝑑𝑣
𝑣 , is the convective acceleration.
𝑑𝑥
33
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Where, 𝑭𝑷𝒙 is the pressure force, 𝑭𝒈𝒙 is the gravity force, 𝑭𝝂𝒙 is the viscosity forces,
𝑭𝒕𝒙 is the turbulence force, 𝑭𝑪𝒙 is the compressibility force, and 𝑭𝝈𝒙 is the surface
tension force.
If pressure force and gravity force are present (other forces are insignificant):
If pressure force, gravity and viscous forces are there (other forces are
insignificant):
If pressure force, gravity, viscous, and turbulent forces are there (other forces
are insignificant):
34
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Assume that we have:
1. A moving fluid (pressurized, open channel,..).
2. Equate the sum of all the forces acting on a fluid element to its (mass
acceleration).
3. Consider a cylindrical element in an arbitrary direction S, with a cross sectional
area dA in a flowing fluid (isolated).
4. The element is oriented at an angle 𝜶 with the horizontal plane (x, y).
5. The surrounding fluid is replaced by pressure forces.
6. Viscous forces are assumed to be zero.
7. The element is accelerated in the S direction.
Applying Newton’s second law in the S direction:
∑ 𝐹 = 𝐹𝑃𝑟𝑒𝑠𝑠𝑢𝑟𝑒 + 𝐹𝑔𝑟𝑎𝑣𝑖𝑡𝑦 = 𝑚𝑎𝐿 (2.7)
Let,
dA = Cross-sectional area of the fluid
element.
at B
𝜕𝑃
𝐹𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒 = 𝑃𝑑𝐴 − (𝑃 + ∆𝑆) 𝑑𝐴 = −∆𝑃𝑑𝐴 (2.7)
𝜕𝑆
- Any pressure forces acting on the side of the cylindrical element will not contribute
to a force in S direction.
- The force due to gravity is the component of weight in the S direction:
35
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
- From the previous
∆𝑍
−∆𝑃 × 𝑑𝐴 − 𝜌𝑔∆𝑆 × 𝑑𝐴 × = 𝜌∆𝑆 × 𝑑𝐴 × 𝑎𝑆 (2.9)
∆𝑆
Divide by: 𝑑𝐴 × ∆𝑆
∆𝑃 ∆𝑍
− − 𝜌𝑔 = 𝜌𝑎𝑆 , (2.10)
∆𝑆 ∆𝑆
−𝝏
(𝑷 + 𝝆𝒈𝒁) = 𝝆𝒂𝑺 (2.11)
𝝏𝑺
Euler equation for motion, for incompressible flow with a constant specific gravity
−𝜕
(𝑃 + 𝜌𝑔𝑍) = 𝜌𝑎𝑡 (2.12)
𝜕𝑆
𝑑𝑉
For steady flow, local acceleration is zero ( = 0, pathline streamline)
𝑑𝑡
𝑑 𝑑𝑉 𝑑 𝑉2
− (𝑃 + 𝜌𝑔𝑍) = 𝜌𝑉 =𝜌 ( ) (2.13)
𝑑𝑆 𝑑𝑆 𝑑𝑆 2
𝑽𝟐
𝑷 + 𝜸𝒁 + 𝝆 = 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕 (2.14)
𝟐
This proves the Bernoulli’s equation.
36
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Assumptions for Bernoulli equation:
1) Incompressible fluid (constant density).
2) Ideal fluid (inviscid), viscosity and the energy losses are neglected.
3) No pumps or turbines within the system.
4) Steady flow.
5) No heat transfer.
Bernoulli Equation: “The total energy per unit of weight of fluid at station1 equals
the total energy per unit weight at station 2”
𝑷𝟏 𝑽𝟐𝟏 𝑷𝟐 𝑽𝟐𝟐
+ 𝒁𝟏 + = + 𝒁𝟐 + (2.15)
𝜸 𝟐𝒈 𝜸 𝟐𝒈
𝑃 𝑉2
Energy Grade Line, EGL= ( + 𝑍 + ): the line connecting the total energy at each
𝛾 2𝑔
point within the hydraulic system is horizontal (constant total energy).
𝑃 𝑉2
The values of the terms of energy ( + 𝑍 + ) can change such that the total
𝛾 2𝑔
value will be constant.
𝑃
Hydraulic Grade Line, HGL= ( + 𝑍): sum of elevation and pressure energies,
𝛾
water level or the elevation which water can rise to.
𝑉2
Velocity head ( ) increases as the diameter decreases to keep continuity
2𝑔
equation valid (v1 A1=v2A2).
For constant diameter, as the elevation Z increases (flow up), the pressure
decreases.
37
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Applying Bernoulli equation, it can be assumed that a unit of weight of fluid is flowing
in a hydraulic system carrying three forms of energy:
𝒎𝒈𝒁
i. Potential energy(
𝒎𝒈
): due to the elevation Z of the unit weight relative to some
reference line. (Work is needed to lift the fluid unit weight a distance Z, it can go down
this distance using this energy).
𝒎𝑽𝟐 𝟏
ii. Kinetic energy(
𝟐 𝒎𝒈
): due to the velocity of the fluid unit V.
𝑷 𝒎𝒈 𝟏 𝑷
iii. Flow energy( 𝜸 𝒎𝒈
): due to its pressure, (𝜸 ) is the elevation that the fluid can raise to.
In real flow systems there is friction in the pipe and elsewhere. This produces heat that
is absorbed by the liquid causing a rise in the internal energy and hence the temperature. In
fact the temperature rise will be very small except in extreme cases because it takes a lot of
energy to raise the temperature. If the pipe is long, the energy might be lost as heat transfer
to the surrounding. Since the equation did not include internal energy, the balance is lost and
we need to add an extra term to the right side of the equation to maintain the balance. This
term is the head lost to friction hL.
38
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Quiz
Find QA, QB, and QC assuming inviscid flow (neglect the resistance and apply the energy
equation between points 1, 2, 3 and the outlet). (Pipes are 4 cm-diameter, and containers
are 80 cm). (Hint: Use V=Q/A).
What is the effect of doubling the area of the exit pipe on the velocity and discharge for the
three cases?
39
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example2.1:
Water flows through a pipe reducer as shown in
the figure. The static pressure at (1) and (2) are
measured by the inverted U-tube manometer
containing oil of specific gravity, SG, less than
one. Determine the manometer reading, h.
Solution:
Assuming steady, inviscid, incompressible
flow Bernoulli equation can be used:
1 1
𝑃1 + 𝜌𝑉12 + 𝛾𝑍1 = 𝑃2 + 𝜌𝑉22 + 𝛾𝑍2
2 2
𝟏 𝑨 𝟐
𝑷𝟏 − 𝑷𝟐 = 𝜸(𝒁𝟐 − 𝒁𝟏 ) + 𝝆𝑽𝟐𝟐 [𝟏 − ( 𝟐) ] (2.16)
𝟐 𝑨 𝟏
The value of (𝑷𝟏 − 𝑷𝟐 ) can be measured by the manometer and can be determined by
using the hydrostatic pressure difference equation:
Or
40
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 2.2:
Water at 60 oF is siphoned from a large tank
through a constant hose as shown.
Determine the maximum height of the hill, H,
over which the water can be siphoned
without the occurring of cavitation. The end
of the siphon is 5 ft. below the bottom of the
tank. Atmospheric pressure is 14.7 Pisa.
Solution:
Assuming steady, inviscid, incompressible
flow Bernoulli equation can be used: Figure 2.5 Figure for example 2.2
𝟏 𝟏 𝟏 𝟏
[𝑷𝟏 + 𝝆𝑽𝟐𝟏 + 𝜸𝒁𝟏 = 𝑷𝟐 + 𝝆𝑽𝟐𝟐 + 𝜸𝒁𝟐 ] Or [𝑷𝟏 /𝜸 + 𝑽𝟐𝟏 + 𝒁𝟏 = 𝑷𝟐 /𝜸 + 𝑽𝟐𝟐 + 𝒁𝟐 ] (2.21)
𝟐 𝟐 𝟐𝒈 𝟐𝒈
Let the tank bottom to be the reference, then Z1 = 15ft, Z2 = H, and Z3 = -5ft.
Also, V1=0 (Large tank), P1=0 (open tank), P3=0 (free jet), and from the continuity
equation A2V2 = A1V1. Also it can be said that V2 = V3 (the hose is constant diameter).
Applying Bernoulli between (1) and (3):
1 2
0 + 0 + 15 = 0 + 𝑉 −5
2𝑔 3
𝑓𝑡 𝑓𝑡
𝑉3 = √2𝑔(𝑍1 − 𝑍2 ) = √(2 × 32 2 [15 − (−5)]) = 35.9 ⁄𝑠 (2.22)
𝑠
Applying Bernoulli between (1) and (2) will give the pressure P2 at the top of the hill:
1 1 1 1
𝑃2 = 𝑃1 + 𝜌𝑉12 + 𝛾𝑍1 − 𝜌𝑉22 − 𝛾𝑍2 = 0 + 0 + 𝛾 (𝑍1 − 𝑍2 ) − 𝜌𝑉22 = 𝛾 (𝑍1 − 𝑍2 ) − 𝜌𝑉22 (2.23)
2 2 2 2
(which is negative because 𝒁𝟏 < 𝒁𝟐 )
The vapor pressure (pressure due to vapor only) of water at 60 oF is 0.256 Pisa. Hence,
for incipient cavitation the lowest pressure in the system will be P = 0.256 Pisa. This lowest
pressure will occur at the top of the hill. Thus: 𝑃2 = 0.256 − 14.7 = −14.4 𝑝𝑠𝑖 (gage).
So H can be found from the last equation:
1
𝑃2 = 𝛾(𝑍1 − 𝑍2 ) − 𝜌𝑉22
2
2 1
(−14.4 𝐼𝑏/𝑖𝑛2 ) (144 𝑖𝑛 ⁄ 2 ) = (62.4 𝐼𝑏⁄ 3 ) (15 − 𝐻)𝑓𝑡 − (1.94 𝑠𝑙𝑢𝑔/𝑓𝑡 3 )(35.9𝑓𝑡/𝑠)2
𝑓𝑡 𝑓𝑡 2
𝐻 = 28.2 𝑓𝑡.
The total energy carried by a unit weight of fluid at station (1) plus the energy added
to the fluid by a pump minus the energy removed due to friction minus the energy removed
by a turbine equal the total energy carried by this unit weight at station (2).
𝑃1 𝑉12 𝑃2 𝑉22
+ 𝑍1 + + 𝐻𝑃 − 𝐻𝑇 − 𝐻𝐿 = + 𝑍2 + (2.24)
𝛾 2𝑔 𝛾 2𝑔
Solution of pipe flow problems requires the application of two principles, the law of
conservation of mass (continuity equation) and the law of conservation of energy
(Bernoulli’s or energy equation).
2.2.1.1 The head loss due to viscous effects in the straight pipes
termed the major loss and denoted HL-major.
Friction loss in a pipe is a complex function of: (i) The system geometry, (ii)
The fluid properties and (iii) The flow rate in the system.
42
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Practical observations showed that the head loss is roughly proportional to
the square of the flow rate in most engineering flows (fully developed,
turbulent pipe flow). This observation leads to the Darcy-Weisbach equation
for head loss due to friction (the derivation is shown next chapter):
𝐿 𝑉2
𝐻𝑓 = 𝑓 (2.25)
𝐷 2𝑔
𝜌𝑉𝐷 𝜀
𝑓 = 𝑓𝑢𝑛𝑐𝑡𝑖𝑜𝑛(𝑅𝑒 = , , 𝑃𝑖𝑝𝑒: 𝑐𝑟𝑜𝑠𝑠 − 𝑠𝑒𝑐𝑡𝑖𝑜𝑛, 𝐺𝑒𝑜𝑚𝑒𝑡𝑟𝑦) (2.26)
𝜇 𝐷
The friction factor is insensitive to moderate changes in the flow and might
be considered constant for fully turbulent flow.
For laminar flow, the head loss is proportional to velocity rather than
velocity squared, thus the friction factor is inversely proportional to velocity.
For turbulent flow, Colebrook (1939) found an implicit correlation for the
friction factor in round pipes (Chapter 3).
The familiar Moody Diagram is a log-log plot of the Colebrook correlation on
axes of friction factor and Reynolds number, combined with the f=64/Re
result from laminar flow.
43
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Calculating Head Loss for a Known Flow:
Obtain the allowable head loss from the energy equation, then start by
guessing a friction factor. (0.02 is a good guess if you have nothing better.)
Calculate the velocity from the Darcy-Weisbach equation. From this velocity
and the piping characteristics, calculate Reynolds Number, relative
roughness and thus friction factor. Repeat the calculation with the new
friction factor until sufficient convergence is obtained. 𝑄 = 𝑉𝐴. (more
details will be introduced later in chapters 3 and 5)
2.2.1.2 The head loss in various pipe components termed the minor
loss and denoted HL-minor.
The head loss designations of “major” and “minor” do not necessarily reflect
the relative importance of each type of loss. For a pipe system that contains many
components and a relatively short length of pipe, the minor loss may actually be
larger than the major loss.
For any type of system there are additional so-called minor losses due to (Fig. 2.7):
1. Exit from a pipe into a tank: 𝐻𝑒𝑥𝑖𝑡 = 𝑉 2 /2𝑔
2. Entry to a pipe from a tank: 𝐻𝑒𝑛𝑡𝑟𝑎𝑛𝑐𝑒 = 0.5 𝑉 2 /2𝑔
3. Sudden expansion: 𝐻𝑒𝑥𝑝𝑎𝑛𝑠𝑖𝑜𝑛 = (𝑉1 − 𝑉2 )2 /2𝑔
2
4. Sudden contraction: 𝐻𝑐𝑜𝑛𝑡𝑟𝑎𝑐𝑡𝑖𝑜𝑛 = (0.4~0.5) 𝑉𝑑𝑜𝑤𝑛𝑠𝑡𝑟𝑒𝑎𝑚 /2𝑔
44
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
These components interrupt the smooth flow of the fluid and cause
additional losses because of the flow separation and mixing they induce. In a
typical system with long pipes, these losses are minor compared to the total head
loss in the pipes (the major losses) and are called minor losses. The losses may not
be so minor; for example, a partially closed valve can cause a greater pressure drop
than a very long pipe. The losses are commonly measured experimentally and
correlated with the pipe flow parameter. The data especially for valve are
somewhat dependent on the particular manufacturer’s design.
Figure 2.8 Fittings that causes minor losses within a pipe networks.
45
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
46
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Practice
Study the following figures for hydraulic piping systems and relate each loss in energy to
the its expected reason as minor or major:
47
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Tutorial 2
Q1. Assuming inviscid flow,
Find the velocity of water at points (2) and (3).
If there is no flow at (3), explain why? And find the elevation where the
water will stop.
48
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
49
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q4. Water flows steadily through the pipe shown such that the pressures at
sections (1) and (2) are 300 kPa and 100 k Pa, respectively. Determine the
diameter of the pipe at section (2), D2, if the velocity at section (1) is 20 m/s
and viscous effects are negligible.
Q5. Water flows from a large tank through a large pipe that splits into two smaller pipes
as shown. If viscous effects are not considered, determine:
a) The total flow rate (Q) from the tank, and
b) The pressure at point (1).
Q6. Water is pumped between the two tanks shown. The energy line is as indicated.
b) Which pipe has larger diameter given that both are made of the same material: A
to pump or B to Pump? Explain!
50
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q7. Water is discharged through the drainpipe at B from the large basin at 0.03 m3/s. lf
the diameter of the drainpipe is d = 60 mm, determine the pressure at B just inside
the drain when the depth of the water is h = 2 m.
51
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q8. Water flows from a tank and through the pipe shown. The flow rate is 1.5 m3/s when
the valve is adjusted to produce the water level in the piezometer tube as indicated.
Find the Darcy-Weisbach friction factor value f in the pipe without using the Moody
diagram (consider the minor losses). What do you expect the pipe material to be?
Q9. What is the pressure in tank A, given that the flow rate between the two tanks
through the two pipes shown is 0.2 m3/s, (the diameters and lengths of the pipes
are as shown). Neglect the minor losses and assume the major losses with the pipe
to be:
HL=0.05 (L/D) × V2/(2g)
52
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
A turbulent flow is one in which the fluid particles rapidly mix as they move a long due to
random three dimensional fluctuations.
53
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Figure 3.2 Flow profiles for real and ideal fluid flows
Figure 3.3 Flow profiles for uniform and varied fluid flows
55
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
- Reynolds found that the transition from laminar to turbulent flow in a pipe depends on
the velocity (V), diameter (D), and the viscosity (μ).
- Reynolds suggested an index number ,(Reynolds number), dimensionless:
𝝆𝑽𝑫 𝑽𝑫
𝑹𝒆 = = (3.1)
𝝁 𝝂
- By experiments:
For a flow in a circular pipe, the critical number of Re is 2000, where the flow
starts to change into turbulent flow.
The transition from laminar to turbulent does not happen at Re =2000 exactly but
varies from 2000 to 4000.
Different values were suggested for Reynolds index where the flow can be
considered to be turbulent. For example, Clayton T. Crow guidelines for Re :
56
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
- Almost, a uniform velocity distribution occurs a cross the pipe for turbulent flow
because of the high-velocity fluid at the pipe centre is transported by turbulent
eddies a cross the pipe to the low velocity region near the wall.
In General:
- Laminar pipe flows are associated with low velocities and turbulent flow with high
velocities.
- Laminar flows can occur in small tubes.
- Turbulent flows are the most common.
- Laminar flow occur in circular pipes when fluid flows in orderly laminae (large
number of thin walled concentric tubes, the outer tube adheres to the pipe wall
while the tube next to it moves with a very slight velocity, where the velocity of
each successive tube increases gradually to reach the maximum velocity near the
centre of the pipe.
- The velocity distribution of laminar flow takes the form of a paraboloid of
revolution with mean velocity equal to one-half of the maximum centre velocity.
- In turbulent flow: the turbulent motion causes the slower water particles adjacent
to the pipe wall to mix continuously with high speed particles in the mid-stream.
The low speed particles near the pipe wall will be accelerated because of
momentum transfer.
The velocity distribution in turbulent flow takes the general form of
logarithmic curve in revolution.
Turbulent mixing activities increase with Reynolds number, {velocity
become flatter as Re increases}.
57
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Figure 3.5 Velocity profiles of laminar and turbulent flows in circular pipes
Example 3.1:
Rainwater runoff from a parking lot flows through a 90-cm-diameter pipe, completely filling
it. Would you expect the flow to be laminar or turbulent? Support your answer with
appropriate calculations.
Solution:
𝑉𝐷
𝑅𝑒 = , If Re 4000 the flow is turbulent.
𝜈
𝑅𝑒 𝜈 4000(1×10−6 ) 𝑚 𝑚𝑚
The corresponding velocity is: 𝑉 = = = 0.0044 = 4.44
𝐷 0.9 𝑠 𝑠
Most likely the velocity will be more than that and the flow will be turbulent flow
Quiz:
58
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
3.2 Developing flow and fully developed flow for pipes:
When a fluid is entering a pipe at a uniform velocity, or after passing through
components as bends, valves, pumps, turbines or similar, the fluid particles in the layer in
contact with the surface of the pipe come to a complete stop due to the no-slip condition.
Due to viscosity of the fluid, this layer in contact with the pipe surface, resists the motion
of adjacent layers and slows them down gradually. For the conservation of mass to hold
true the velocity of middle layers of the fluid in the pipe increases (since the layers of fluid
near the pipe surface have reduced velocities). This develops a velocity gradient across the
cross section of the pipe. A fluid needs some length to develop the constant velocity profile.
Boundary layers keep on growing till they meet after some distance downstream from the
entrance region. After this distance velocity profile doesn't change, flow is said to be fully
developed.
𝐿𝑒
= 0.05 𝑅𝑒 𝑓𝑜𝑟 𝑅𝑒 < 2000 (𝐿𝑎𝑚𝑖𝑛𝑎𝑟 𝑓𝑙𝑜𝑤), 𝑤ℎ𝑒𝑟𝑒 𝐷 𝑖𝑠 𝑡ℎ𝑒 𝑝𝑖𝑝𝑒 𝑑𝑖𝑎𝑚𝑒𝑡𝑒𝑟 (3.2)
𝐷
𝐿𝑒
= 50 𝑓𝑜𝑟 𝑅𝑒 > 3000 (𝑇𝑢𝑟𝑏𝑢𝑙𝑒𝑛𝑡 𝑓𝑙𝑜𝑤) (3.3)
𝐷
59
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 3.2:
Consider fluid flowing in a round tube of length 1 m and diameter of 5mm. Classify the flow
as laminar or turbulent and calculate the entrance length for:
(a) air (50 oC) with a speed of 12 m/s, and
(b) water (15 oC) with a mass flow rate of 𝑚̇ = 8 g/s.
The pipe is connected to a reservoir with a smooth and tapered entrance.
Solution:
Air @ (50 oC), will have ν =1.7910-5 m2/s.
Water @(15 oC), will have μ = 1.1410-3 N.s/m2.
a) Air
𝑚
𝑉𝐷 (12 𝑠 ) (0.005𝑚)
𝑅𝑒 = = = 3350 (𝑡𝑢𝑟𝑏𝑢𝑙𝑒𝑛𝑡)
𝜈 −5 𝑚2
1.79 × 10
𝑠
b) Water
4 𝑚̇ 4(0.008 𝑘𝑔/𝑠)
𝑅𝑒 = = = 1787 (𝑙𝑎𝑚𝑖𝑛𝑎𝑟)
𝜋𝐷𝜇 𝜋(0.005𝑚)(1.14 × 1.14 × 10−3 𝑁. 𝑠/𝑚2 )
𝑉𝐷 (𝜌𝑄/𝐴)𝐷 𝑚̇𝐷 4 𝑚̇
[𝑁𝑜𝑡𝑒: 𝑅𝑒 = = =𝜋 = ]
𝜈 𝜇 𝐷2 𝜇 𝜋𝐷𝜇
4
60
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
3.3 Stress distribution in pipe flow: (laminar and turbulent flow)
Required: to derive the equation for the shear-stress distribution on a plane oriented normal
to streamlines for laminar and turbulent flows in circular pipes.
- Consider the flow of a Newtonian fluid in a round tube that is inclined at an angle
with horizontal.
- Assume that the flow is fully developed, steady and laminar.
- A control volume is assumed (cylindrical) of length ∆s and radius r.
Figure 3.6: The forces in a cylindrical control volume within the flow of a round tube.
𝒓 𝒅(𝑷+ 𝜸 𝒁)
𝝉 = [− ] (3.4)
𝟐 𝒅𝒔
Which is the relation of the shear stress distribution on a plane inserted normal to the
direction of flow of pipe.
Shear stress is varying linearly with the radius of the pipe r.
Shear stress is zero at the centre line and maximum at the wall of the pipe.
61
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Figure 3.7 Shear stress distributions over normal area to the flow in a round pipe
“For a fully developed flow in a pipe (laminar or turbulent), the shear stress distribution
on an area that is normal to the stream line of the flow is linear starting from zero at the
centre and increase in the direction of the pipe wall”.
∆(𝑃+𝛾𝑍) ∆(ℎ) 𝑝
𝜏 = 2𝑟 [− ] = 2𝑟 [−𝛾 ] , (𝑤ℎ𝑒𝑟𝑒 ℎ = + 𝑧) (3.5)
∆𝐿 ∆𝐿 𝛾
Quiz:
Hint:
62
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
3.4 Velocity profile for laminar flow in a round tube:
Recall:
Laminar flow is a flow regime in which fluid motion is smooth and occurs in layers,
and the mixing occurs in layers by molecular diffusion (much slower than turbulent
mixing), Re 2000.
Laminar flow is important for flow in small conduits called micro channels, for
lubrication flow, and other flows in which viscous forces are dominant.
Laminar flow in round tubes is called Poiseuille flow or Hagen-Poiseulle flow.
𝑑𝑉 𝑑𝑉 𝑑𝑟 𝑑𝑉
𝜏=𝜇 =𝜇 = −𝜇 (3.7)
𝑑𝑌 𝑑𝑟 𝑑𝑌 𝑑𝑟
𝑟 𝑑(𝑃+ 𝛾 𝑍)
o Using 𝜏 = [− ]
2 𝑑𝑠
2𝜇 𝑑𝑉 𝑑(𝑃+ 𝛾 𝑍)
− ( ) ( ) = [− ] (3.9)
𝑟 𝑑𝑟 𝑑𝑠
To have the left side a function of r and the right side function of s, then both sides
are equal to a constant.
𝑑(𝑃+ 𝛾 𝑍) ∆(𝑃+ 𝛾 𝑍) 𝛾 ∆(𝑃+ 𝛾 𝑍) ∆ℎ
[ = = = 𝛾 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡] , where ∆h is the change
𝑑𝑠 ∆𝐿 𝛾 ∆𝐿 ∆𝐿
𝑃 𝑉2
in fluid energy head ( + 𝑍 + ) over a length ∆L.
𝛾 2𝑔
𝑑𝑉 𝑟 ∆ℎ
= ( ) [ 𝛾 ] , integrate:
𝑑𝑟 2𝜇 ∆𝐿
𝑟2 ∆ℎ
𝑉 = ( )[ 𝛾 ] + 𝐶 (3.10)
4𝜇 ∆𝐿
Apply the no slip condition {V=0 at the wall of the pipe, or V(ro)=0}.
𝑟2 ∆ℎ
0 = ( 𝑜)[ 𝛾 ] + 𝐶
4𝜇 ∆𝐿
63
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
The equation for velocity profile in a pipe for laminar flow:
𝟐 𝟐 𝟐 𝟐
𝑽(𝒓) = − (𝒓𝒐𝟒𝝁
−𝒓
) (𝜸∆𝒉
∆𝑳
) = − (𝒓𝒐𝟒𝝁
−𝒓
) (∆(𝑷+𝜸𝒁)
∆𝑳
) (3.11)
𝒓 𝟐
𝑽(𝒓) = 𝑽𝒎𝒂𝒙 (𝟏 − ( ) ) (3.13)
𝒓𝒐
This indicates that the velocity profile for laminar flow in a pipe is parabolic.
3.4.2 Finding the discharge and the mean velocity for laminar flow in a
round tube:
𝑟 2 2
𝑄 = ∫ 𝑉𝑑𝐴 = ∫0 𝑜 − (𝑟𝑜4𝜇
−𝑟
) (𝛾∆ℎ
∆𝐿
)(2𝜋𝑟𝑑𝑟) (3.14)
2 𝑟𝑜
𝜋 𝛾∆ℎ 𝑟 2 −𝑟𝑜2 𝜋𝑟𝑜4 𝛾∆ℎ
𝑄 = (− ) ( ∆𝐿 ) ( ) | = −( ) ( ∆𝐿 ) (3.15)
4𝜇 2 8𝜇
0
𝑸 𝒓𝟐 𝜸∆𝒉 𝑽𝒎𝒂𝒙
𝑽= = − ( 𝒐) ( )= (3.16)
𝑨 𝟖𝝁 ∆𝑳 𝟐
64
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
3.5 Head loss and friction factor (f) for laminar flow:
o Assume fully developed flow in the pipe.
o Apply energy equation between section (1) and (2).
65
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
̅
32 𝜇 𝐿 𝑉
ℎ𝑓 = (3.19)
𝛾𝐷2
{Laminar, fully developed, steady, Newtonian flow}.
The head loss for laminar flow varies linearly with velocity, it is influenced by
viscosity, pipe length, specific weight of the fluid, and pipe diameter squared.
Extra Reading
The velocity distribution for laminar flow in a pipe is given by the equation:
V = Vmax[1 – (r / r0)2]
Here r0 is the radius of the pipe and r is the radial distance from the center. Determine V in
terms of Vmax and evaluate the kinetic energy correction factor α.
Solution:
A sketch for the velocity distribution is shown in the figure. The discharge is given by
Q = ∫V dA, or
Q = VA = (r)0∫0 Vmax {1 – (r2 / r02)} 2πr dr
or
The mean velocity is one-half the maximum velocity. To evaluate α we apply the equation
α = (1/A) ∫A (V / V)3 dA
α = (1 / r02) (r)0∫0 [{Vmax[1 – (r2 / r02)]3} / {(1/2)3 Vmax3}] 2πr dr
= 2.
66
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝐷2 𝐷2
(𝑃1 − 𝑃2 ) ( ) − 𝑜 ( 𝐷 ∆𝐿) − [( ) ∆𝐿] 𝑠𝑖𝑛() = 0 (3.20)
4 4
67
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝜌𝑉2
𝐿 4𝜏𝑜 4𝜏 𝐿 𝑉2
ℎ𝑓 = ℎ𝑒𝑎𝑑𝑙𝑜𝑠𝑠 𝑖𝑛 𝑝𝑖𝑝𝑒 = ( 𝜌𝑉2 ) (
𝐷 𝛾
2
) = ( 𝜌𝑉𝑜2 ) (𝐷) 2𝑔 (3.25)
2 2
𝜌𝑉 2
Given that o is related to
2
Let
4𝜏𝑜 𝑠ℎ𝑒𝑎𝑟 𝑠𝑡𝑟𝑒𝑠𝑠 𝑎𝑐𝑡𝑖𝑛𝑔 𝑎𝑡 𝑡ℎ𝑒 𝑤𝑎𝑙𝑙
𝑓= = (3.26)
𝜌𝑉 2 /2 𝑘𝑖𝑛𝑒𝑡𝑖𝑐 𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒
𝑳 𝑽𝟐
𝒉𝒇 = 𝒇 (3.27)
𝑫 𝟐𝒈
For fully developed incompressible, and steady flow, laminar or turbulent, the
head loss, or pressure loss due to friction along a given length of pipe depends on
the friction factor, pipe length, diameter, mean velocity squared.
32 𝜇 𝐿 𝑉̅ 𝐿 𝑉2
ℎ𝑓 = =𝑓
𝛾𝐷2 𝐷 2𝑔
32𝜇𝐿𝑉 𝐷 2𝑔 64𝜇
𝑓= = = 64/𝑅𝑒 (3.28)
𝛾𝐷2 𝐿 𝑉2 𝜌𝑉𝐷
68
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
3.7 Power required to maintain a flow:
The power Ρ is the rate of doing work.
For horizontal pipe (Z is constant):
𝑤𝑜𝑟𝑘 𝑓𝑜𝑟𝑐𝑒.𝑑𝑖𝑠𝑡𝑎𝑛𝑐𝑒
Ρ=
𝑢𝑛𝑖𝑡 𝑡𝑖𝑚𝑒
=
𝑢𝑛𝑖𝑡 𝑡𝑖𝑚𝑒
= 𝐹𝑜𝑟𝑐𝑒 × 𝑉𝑒𝑙𝑜𝑐𝑖𝑡𝑦 (3.29)
−𝑑𝑝
P=( ) (𝐴 × 𝐿 × 𝑉 ) = (𝑃 1 − 𝑃 2 )𝑄 (3.30)
𝑑𝑥
For inclined pipe (general equation) , peizometric pressures are used instead of pressure:
𝛲 = 𝛾 (ℎ1 − ℎ2 )𝑄 (3.31)
1 𝑒 2.51
= −2𝑙𝑜𝑔( + ) (3.32)
√𝑓 3.7𝐷 𝑅𝑒√𝑓
69
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
70
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
71
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
72
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
73
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
74
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
75
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
76
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
77
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
The following characteristics are observed from analysing the Moody diagram:
1. For laminar flow, Re < 2000, f = 64/Re, which is independent of relative roughness (e/D
or /D).
2. Turbulent flow is strongly affected by
roughness. For very large Reynolds
numbers (to the right of the dashed line), f
is a function of relative roughness (e/D or
/D) only, which is independent of the
Reynolds number. For such flows,
commonly termed completely turbulent
(fully rough flow), along the wall pipe, the
laminar sub layer is so thin that the surface
roughness completely dominates the
character of the flow near the wall (sub
layer totally broken up and friction
independent of Reynolds number). The
roughness curves are nearly horizontal in
this fully rough regime. Figure 3.11 Viscous sublayer compared to roughness
3. The gap in the figure for which no values of f are given, 2000 < Re < 4000, is a result of
the fact that the flow in this transition range may be laminar or turbulent (or an
unsteady mix of both) depending on the specific circumstances involved. There are no
reliable friction factors in this range.
4. Even for smooth pipes the friction factor is not zero. That is, there is a head loss in any
pipe, no matter how smooth the surface is made. This is a result of the no-slip boundary
condition that requires any fluid to stick to any solid surface it flows over. There is
always some microscopic surface roughness that produces the no-slip behaviour (and
thus f ≠ 0) on the molecular level, even when the roughness is considerably less than
the viscous sub layer thickness. Such pipes are called hydraulically smooth (f = function
of (Re,)).
5. Colebrook (in 1939) combined all data for transition and turbulent flow in smooth as
well as rough pipes into the following relation known as Colebrook equation:
1 𝑒 2.51
= −2𝑙𝑜𝑔( + )
√𝑓 3.7𝐷 𝑅𝑒√𝑓
In 1944, Lewis Ferry Moody plotted the Darcy–Weisbach friction factor into what is now
known as the Moody chart.
6. For smooth pipes (turbulent flow: 4000Re105 ),Blassius experimentally developed the
0.316
following formula: 𝑓 = 0.25
𝑅𝑒
78
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝑃1 −𝑃2 𝑉12 𝑉22
Considering the energy equation: + 𝑍1 − 𝑍2 + − = ℎ𝐿
𝛾 2𝑔 2𝑔
For horizontal, constant diameter pipe:
𝑃1 −𝑃2
∆Z=0, and V1=V2 = ℎ𝐿 → ∆𝑃 = 𝑃2 − 𝑃1 = −𝛾ℎ𝐿
𝛾
Also the Darcy-Weisbach equation for laminar flow:
̅
32 𝜇 𝐿 𝑉
ℎ𝑓 = 𝑄 = 𝑉̅ 𝐴 = 𝜋𝛾𝐷4 ∆𝑃/(8𝜇𝐿)
𝛾𝐷2
These results confirm that for a laminar horizontal pipe flow rate, the flow rate is:
i. Directly proportional to the pressure drop.
ii. Inversely proportional to the viscosity.
iii. Inversely proportional to the pipe length.
iv. Proportional to the pipe diameter to the fourth power, (strong relation).
For inclined constant diameter pipe, the value of ∆P will be:
∆𝑃 = 𝑃2 − 𝑃1 = −𝛾ℎ𝐿 + 𝛾(𝑍1 − 𝑍2 )
More concern is needed to the change in elevation. Cases can be seen were the pressure
could be zero, positive or negative for the the same flow losses depending on the change in
elevation.
Example 3.3:
If the diameter of a pipe is doubled, what effect does this have on the flow rate for a given
head loss? Consider (a) Laminar flow; (b) turbulent flow.
Solution:
̅
32 𝜇 𝐿 𝑉 𝑉
(a) Laminar flow: ℎ𝑓 = = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡 × ,
𝛾𝐷2 𝐷2
V=K1D2 ; Q=AV=D2/4 (KD2)=K2D4
Thus doubling the diameter will increase the flow by a factor of 24=16.
𝐿 𝑉2
(b) Turbulent flow, ℎ𝑓 = 𝑓 ,
𝐷 2𝑔
𝑉2
i) assume fully turbulent (constant f), ℎ𝑓 = 𝐾2 ; 𝑉 = 𝐾3 𝐷1/2 and Q=AV=K3D5/2
𝐷
Thus doubling the diameter will increase the flow rate by a factor of 22/5 = 5.66
for turbulent flow with complete turbulence.
0.316
ii) For smooth pipe with Re100,000, Blasius equation can be used (𝑓 = 0.25
𝑅𝑒
𝐾5
=(𝐷𝑉)0.25);
𝐿 𝑉2 𝑉 7/4
ℎ𝑓 = 𝑓 = 𝐾6 ; thus Q=AV=K8D19/7
𝐷 2𝑔 𝐷5/4
Which means that doubling the diameter will increase the flow rate by a factor
of 219/7=6.56 for turbulent flow in smooth pipes.
79
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 3.4:
Oil (Sp. Gr. =0.94) with a viscosity of 0.0096 N.S/m2 flows in a 90-mm-diameter welded-
steel pipe (e=0.046mm) at a discharge of 7.2 L/s. What is the friction head loss per meter
of pipe?
Solution:
V=Q/A=Q/(D2/4)= 4(0.0072)/( 0.0902)=1.132 m/s
Re=0.09(1.132) (10000.94)/0.0096=9974; e/D=0.046/90=0.000511.
ℎ𝑓 0.032 1.1322
From Moody Diagram, f = 0.032 = = 0.0232 𝑚/𝑚
𝐿 0.09 2(9.81)
Example 3.5:
Crude oil (Sp. Gr. =0.855) at 50oC (ν= 3.810-6 m2/s), flows at 300 L/s through a 450-mm-
diameter pipe (e=0.054 mm) 1500 m long. Find the power loss within the flow.
Solution:
V=0.3/( 0.2252)= 1.886 m/s ; Re=VD/ν =0.45 (1.886)/(3.810-6)= 2.23105 ; e/D=0.054/450
𝐿 𝑉2 1500 1.8862
= 0.00012; from Moody Diagram , f=0.016 ; ℎ𝑓 = 𝑓 = 0.016 = 9.832 𝑚
𝐷 2𝑔 0.45 2(9.81)
Power loss = Q hf = (0.8859807)0.30(9.832) =25.6 kW.
Example 3.6
A fluid flows through a smooth horizontal 2–m long tube of 2 mm diameter with an
average velocity of 2.1 m/s. Determine the head loss and the pressure drop if the fluid is
(a) air, (b) water, (c) mercury. Use the information in the following table.
Solution:
𝑃1 𝑉12 𝑃2 𝑉22
+ 𝑍1 + = + 𝑍2 + + 𝐻𝐿
𝛾 2𝑔 𝛾 2𝑔
𝐿 𝑉2
where, 𝐻𝐿 = 𝑓 , Z1=Z2, and V1=V2
𝐷 2𝑔
𝑚 2
2𝑚 (2.1 )
𝑠
Thus, ℎ𝐿 = 𝑓 ( ) = 225 𝑓
0.002𝑚 2𝑔
And ∆P=P1 – P2 = hL
80
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝑚
𝑉𝐷 (2.1 )(0.002𝑚) 4.2×10−3
𝑠
Also, 𝑅𝑒 = = =
𝜈 𝜈 𝜈
Fluid 𝜈 (m2/s) Re Flow type f hL (m) (N/m3) ∆P (N/m2)
(a)Air 1.4610-5 287 Laminar 64/Re=0.223 50.2 12.0 602
(b)Water 1.1210-6 3750 turbulent 0.0404 9.09 9800 8.91104
(c) Hg 1.1510-7 36500 turbulent 0.0220 4.95 133000 6.58105
Example 3.7:
A viscous fluid flows in a 0.10-m diameter pipe such that its velocity measured 0.012 m away
from the pipe wall is 0.8 m/s. If the flow is laminar, determine the centreline velocity and the
flow rate.
Solution:
𝑟 2
𝑉 (𝑟) = 𝑉𝑚𝑎𝑥 (1 − ( ) ), where D=0.1m.
𝑟𝑜
(0.038𝑚 2
Thus, 0.8 𝑚/𝑠 = 𝑉𝑐𝑒𝑛𝑡𝑒𝑟 (1 − ( ) ) Vc=1.89 m/s.
0.1𝑚/2
𝜋 2 𝜋 2
1 −3
𝑚3
𝑄 = 𝐷 𝑉 = 0.1 ( × 1.89) = 7.42 × 10
4 4 2 𝑠
Example 3.8:
Oil (Sp. Gr.=0.9 , μ=510-1 N.s/m2) flows steadily in a 3-cm pipe. The pipe is vertical, and the
pressure at an elevation of 100 m is 200 kPa. If the pressure at an elevation of 85 m is 250 kPa,
is the flow direction up or down? What is the velocity at the centre of the pipe and at 6 mm
from the centre, assuming that the flow is laminar?
Hints:
The flow will be always from the highest
energy (P/+Z+V2/2g) to the lowest energy.
In case of horizontal constant-diameter pipe,
the flow will be from the highest pressure to
the lowest pressure.
The direction of Z will be always upward Z.
The direction of the flow (s or X), where we
measure L, is at the same line of the pipe and
at the direction from the highest energy to the
lowest energy (flow direction)
For normal pipe, the equation for velocity and
pressure profile for laminar flow can be given Figure
as the following:
81
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 3.9:
Oil (Sp. Gr. =0.85) with kinematic viscosity of 610-4 m2/s flows in a 15 cm pipe at a rate of
0.02 m3/s.
1. What is the head loss per 100 m length of pipe?
2. Find the centre and the average velocity of the flow.
3. Find the wall shear stress.
4. Plot the velocity and the shear profiles.
Solution:
(1).V=Q/A=(0.02 m3/s)/( (0.15m)2/4)= 1.13 m/s
Re=VD/ν =(1.13 m/s)(0.15m)/(610-4 m2/s)=283 2000 flow is laminar
32𝜇𝐿𝑉 32𝜌𝜈𝐿𝑉 32(6×10−4 )×100×1.13
ℎ𝑓 = = = = 9.83 𝑚
𝛾𝐷2 𝜌𝑔𝐷2 9.81×0.152
𝐿 𝑉2 64 𝐿 𝑉 2 64 100 1.132
Or, ℎ𝑓 = 𝑓 = = = 9.83
𝐷 2𝑔 𝑅𝑒 𝐷 2𝑔 283 0.15 2𝑔
82
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 3.10:
A fluid flows through a horizontal 0.25-cm diameter pipe. when the Reynolds number is
1500, the headloss over a 6-m length of the pipe is 192 cm. Determine the fluid velocity.
Solution:
𝐿 𝑉2 64 𝐿 𝑉 2 64 6 𝑉2
ℎ𝐿 = 𝑓 = = = 1.92 𝑚 , (Laminar flow) V=0.606 m/s
𝐷 2𝑔 𝑅𝑒 𝐷 2𝑔 1500 0.0025 2𝑔
Example 3.11:
For laminar flow in a circular pipe, find the velocities at 0.1r, 0.3r, 0.5r, 0.7r, and 0.9r relative
to the center velocity. Plot the velocity profile.
Solution:
Given that the velocity at any diameter of the pipe for laminar flow can be given as:
𝑟 2 𝑉(𝑟) 𝑟 2
𝑉(𝑟) = 𝑉𝑚𝑎𝑥(1 − ( ) ) , or = (1 − ( ) )
𝑟𝑜 𝑉𝑚𝑎𝑥 𝑟𝑜
which can be used to find the relative velocity at any radius of the pipe section:
83
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 3.12:
An oil with a viscosity of μ=0.4 N.s/m2 and density =900 kg/m3 flows in a pipe of diameter
D=0.020 m.
(a) What pressure drop, P1 – P2, is needed to produce a flowrate of Q = 2.010-5 m3/s
for horizontal pipe with x1=0 and x2=10m.
(b) How steep a hill, , must the pipe be on if the oil is to flow through the pipe at the
same rate in part (a), but with P1=P2
(c) If P1=200kPa, what is the pressure at section x3=5m, where x is measured along the
pipe.
Solution:
(a) If the Reynolds number is less than 2000 the flow is laminar.
V=Q/A=(210-5 m3/s)/((0.022 m2)/4)=0.0637 m/s
Re=VD/μ=2.87 2000 the flow is laminar. And the following equation for
̅
32 𝜇 𝐿 𝑉 128 𝜇 𝐿𝑄
laminar flow can be used: ℎ𝐿 = =
𝛾𝐷2 𝛾𝜋𝐷4
The pipe is horizontal, then Z1=Z2. Also, V1=V2 (constant diameter),then the energy
equation can be given as:
𝑃1 𝑃2 𝑉22 𝑉12 𝑃 𝑃2
+ 𝑍1 − 𝑍2 = + ℎ𝐿 + − , or ( 1 − = ℎ𝐿 )
𝛾 𝛾 2𝑔 2𝑔 𝛾 𝛾
128 𝜇 𝐿𝑄
(𝑃1 − 𝑃2 ) = 𝛾ℎ𝐿 =
𝜋𝐷4
128×0.4×10×2×10−5 𝑁
So, ∆𝑃 = 𝑃1 − 𝑃2 = = 20,400 = 20.4 𝑘𝑃𝑎
𝜋(0.02)4 𝑚2
84
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
(b) We have the pipe on a hill of angle such that the pressure drop, ∆P=0, given that
he flow does not change. The approach requires that the value of hL does not
change. the only way to achieve that is to have the change in elevation replace the
pressure drop:
85
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Extra Problems
86
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
87
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
88
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
89
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
90
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Tutorial 3
Q1. Oil having specific gravity 0.89 and µ=0.085 poise flows through a horizontal pipe
of 1.5 cm diameter and 230 cm length. Compute the pressure drop, pressure and
energy gradients in order to maintain a flow of 150 cm3/s.
Q2. It is required to transport 55 L/s of an oil, sp. gr.0.82 and µ=2.1 N.s/m2, to a
distance of 1 km. Permissible energy gradient= 0.04. Determine the size of the
pipe and power required to maintain the flow.
Q3. In a flow through an 8 cm-diameter horizontal pipe with glycerine, ρ=1400 kg/m3
and µ=7.9 Poise, the shear stress at the wall was found to be 300 N/m2.
Determine pressure and energy gradients, mean velocity of flow and Reynolds
number. Also plot the shear stress and velocity distribution diagrams.
Q4. Oil (µ=0.191 N. S / m2, ρ=927.7 kg/m3) flows through an inclined 7.5 cm-diameter
pipe. Find the slope of the pipe that will cause the pressure to be constant along
the pipe's entire length.
Q5. An Oil (µ=8×10-3N.s/m2, sp. gr. =0.88) is flowing through a round tube as shown.
Assuming laminar flow, determine:
a) Direction of flow.
b) Loss of head between 1 and 2.
c) Average velocity of flow
d) Reynolds number.
e) Maximum velocity.
f) Maximum shear stress.
g) Power required to maintain the flow
h) Check that the flow is laminar flow
91
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Solution of Tutorial 3:
(1) Fluid { Oil, S.Gr.=0.89 , µ=0.085 poise= (0.085/10) N.S/m2 }
𝑐𝑚3
150
𝑉=𝜋 𝑠 = 84.9 𝑐𝑚 = 0.85 𝑚/𝑠
2 𝑠
4 (1.5𝑐𝑚)
𝜌𝑉𝐷 (0.89 × 1000) × 0.85 × 0.015
𝑅𝑒 = = = 1335 (𝑙𝑎𝑚𝑖𝑛𝑎𝑟 𝑓𝑙𝑜𝑤)
𝜇 0.0085
𝐿 𝑉 2 64 𝐿 𝑉 2 64 2.3 (0.85)2
ℎ𝑓 = 𝑓 = = = 0.27 𝑚
𝐷 2𝑔 𝑅𝑒 𝐷 2𝑔 1335 0.015 2𝑔
𝑷𝟏 𝑽𝟐𝟏 𝑷𝟐 𝑽𝟐𝟐
+ + 𝒁𝟏 = + + 𝒁𝟐 + 𝒉𝒇
𝜸 𝟐𝒈 𝜸 𝟐𝒈
𝑷𝟏 − 𝑷𝟐 (𝑽𝟐𝟏 − 𝑽𝟐𝟐 )
+ (𝒁𝟏 − 𝒁𝟐 ) + = 𝒉𝒇
𝜸 𝟐𝒈
𝑃1 − 𝑃2
= ℎ𝑓 , (𝑉1 = 𝑉2 𝑎𝑛𝑑 𝑍1 = 𝑍2 )
𝛾
𝑁 𝑁
𝑃1 − 𝑃2 = 𝛾 × ℎ𝑓 = 0.27𝑚 (0.89 × 9810 3 ) = 2357.343 2
𝑚 𝑚
𝟐
∆𝑷 = −𝟐. 𝟑𝟔 𝒌𝑵/𝒎
∆𝑃 −2.36 𝑘𝑁/𝑚2
𝑃𝑟𝑒𝑠𝑠𝑢𝑟𝑒 𝑔𝑟𝑎𝑑𝑖𝑒𝑛𝑡 = =
∆𝐿 2.3 𝑚
ℎ2 − ℎ1 −0.27 𝑚
𝐸𝑛𝑒𝑟𝑔𝑦 𝑔𝑟𝑎𝑑𝑖𝑒𝑛𝑡 = =
∆𝐿 2.3 𝑚
92
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Required: D , Power
𝑄 55 × 10−3 0.07
𝑉= = 𝜋 = 2
𝐴 𝐷2 𝐷
4
ℎ𝐿 32 𝜇𝑉
𝒂𝒔𝒔𝒖𝒎𝒆 𝑳𝒂𝒎𝒊𝒏𝒂𝒓 𝒇𝒍𝒐𝒘 → =
𝐿 𝛾𝐷2
Use the maximum allowable energy gradient (0.04) to find the minimum
allowable diameter of the pipe.
0.07 1 1
0.04 = 32 × 2.1 × 2 2
𝐷 𝐷 0.82 × 9810
0.07
𝜌𝑉𝐷 0.82×1000×[ ]×0.34
(0.34)2
Check , 𝑅𝑒 = = = 9.5 (𝑙𝑎𝑚𝑖𝑛𝑎𝑟)
𝜇 2.1
𝑷𝒐𝒘𝒆𝒓 = 𝜸 (𝒉𝟏 − 𝒉𝟐 )𝑸 = (𝟎. 𝟖𝟐 × 𝟗𝟖𝟏𝟎)(𝟎. 𝟎𝟒 × 𝟏𝟎𝟎𝟎)(𝟓𝟓 × 𝟏𝟎−𝟑 )
= 𝟏𝟕. 𝟔𝟗 𝒌𝒘𝒂𝒕𝒕
93
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
(3) Fluid { Glycerine , ρ =1400 kg/m3, µ=7.9 poise } , D = 8 cm, 𝜏𝑤 =300 N/m2
Required: pressure and energy gradient, average velocity, Re, 𝜏 and velocity
diagrams.
−𝑟𝑜 ∆(𝑃 + 𝛾𝑍)
𝜏𝑤 =
2 ∆𝐿
300𝑁 0.04 ∆(𝑃 + 𝛾𝑍)
=−
𝑚2 2 ∆𝐿
Given that the pipe is horizontal , then ΔZ=0.0 →
∆𝑃 𝑁
= −15000 2
∆𝐿 𝑚 .𝑚
∆ℎ ∆𝑃 1 −15000 𝑁/𝑚3
= = = −1.092 𝑚/𝑚
∆𝐿 ∆𝐿 𝛾 9.81 × 1400 𝑁/𝑚3
𝐷2 ℎ ℎ ∆ℎ
𝑉̅ = 𝛾 , 𝑛𝑜𝑡𝑒 𝑡ℎ𝑎𝑡 = −
32𝜇 𝐿 𝐿 ∆𝐿
(0.08) 2
𝑉̅ = (1400 × 9.81)(1.092) = 3.796 𝑚/𝑠
7.9
32
10
𝜌𝑉𝐷 1400 × 3.796 × 0.08
𝑅𝑒 = = = 538.3 (𝑙𝑎𝑚𝑖𝑛𝑎𝑟)
𝜇 0.79
−𝑟 ∆(𝑃 + 𝛾𝑍) −𝑟
𝜏= = [−15000] = 𝑟[7500]
2 ∆𝐿 2
𝑟 2 𝑟 2
𝑉(𝑟) = 𝑉𝑚𝑎𝑥 (1 − [ ] ) = (2 × 3.796) (1 − [ ] )
𝑟𝑜 0.04
94
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
95
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
(4)
𝑄 0.0085 𝑚3 /𝑠
𝑉= = 𝜋 = 1.92 𝑚/𝑠
𝐴 (0.075) 2
4
𝜌𝑉𝐷 927.7 × 1.92 × 0.075
𝑅𝑒 = = = 699.4 (𝑙𝑎𝑚𝑖𝑛𝑎𝑟)
𝜇 0.191
Applying energy equation between (1) and (2):
𝑃1 𝑉12 𝑃2 𝑉22 32 𝜇 𝑉 𝐿
+ + 𝑍1 = + + 𝑍2 +
𝛾 2𝑔 𝛾 2𝑔 𝛾 𝐷2
𝑃1 𝑃2
= , V1=V2
𝛾 𝛾
𝑍1 − 𝑍2 32 𝜇 𝑉
𝑆𝑜,
=
𝐿 𝛾 𝐷2
𝑄 0.0085 𝑚3 /𝑠
𝑉= = 𝜋 = 1.92 𝑚/𝑠
𝐴 (0.075)2
4
𝑍1 − 𝑍2 32 𝜇 𝑉 32 × 0.191 × 1.92
sin(𝜃) = = = = 0.229
𝐿 𝜌𝑔 𝐷2 927.7 × 9.81 × (0.075)2
→ 𝜃 = 13.2𝑜
96
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
32 𝜇 𝑉 𝐿 𝛾 ℎ𝑓 𝐷 2 (0.88×9810)(19.33)(0.0016)2
c) ℎ𝑓 = → 𝑉= = = 0.236 𝑚/𝑠
𝛾 𝐷2 32 𝜇𝐿 32 (8×10−3 )(5/ sin(45𝑜 ))
𝜌𝑉𝐷 (0.88×1000)(0.236)(0.0016)
d) 𝑅𝑒 = = = 41.536 (𝑙𝑎𝑚𝑖𝑛𝑎𝑟 𝑓𝑙𝑜𝑤)
𝜇 8×10−3
0.0016 2
−𝑟𝑜2 ∆ℎ ( ) 42.06−61.3
2
e) 𝑉𝑚𝑎𝑥 = 𝑉𝑐𝑒𝑛𝑡𝑟𝑒 = 𝛾 =− −3
(0.88 × 9810) =0.47 m/s
4𝜇 ∆𝐿 4(8×10 ) 5/ sin(45𝑜 )
or 𝑉𝑚𝑎𝑥 = 𝑉𝑐𝑒𝑛𝑡𝑟𝑒 = 2𝑉̅ = 2 × 0.236 = 0.47 𝑚/𝑠
0.0016
−𝑟 ∆(𝑃+𝛾𝑍) ( )
2
f) 𝜏𝑚𝑎𝑥 = 𝜏𝑏𝑜𝑢𝑛𝑑𝑎𝑟𝑦 = 𝜏𝑤𝑎𝑙𝑙 = =− (0.88 ×
2 ∆𝐿 2
42.06−61.3
9810) = 11.3 N/m2
5/ sin(45𝑜 )
𝑁 𝜋
g) 𝑃𝑜𝑤𝑒𝑟 = (0.88 × 9810 𝑚3) (61.39 − 42.06)𝑚 [0.236 × (0.0016)2 ] = 0.316 𝑊𝑎𝑡𝑡.
4
97
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Quiz
Answer Sheet
The flow in pipes will always take place in the direction of decreasing
𝑷 𝑽𝟐
energy ( + 𝒁 + ) (decreasing pressure for the case of horizontal pipe).
𝜸 𝟐𝒈
To convert the dynamic viscosity from poise into N.s/m2, divide by 10.
{ µ= 1 Pascal. Second = kg/(m.s)= N.s/m2 = 10 poise }.
Power required to maintain a flow:
𝑃𝑜𝑤𝑒𝑟 = (𝑃1 − 𝑃2 )𝑄 (Horizontal pipe)
𝑝𝑜𝑤𝑒𝑟 = 𝛾 (ℎ1 − ℎ2 )𝑄 (Inclined pipe)
𝑃1 𝑉12 𝑃2 𝑉22
Energy equation: + 𝑍1 + = + 𝑍2 + + ℎ𝐿
𝛾 2𝑔 𝛾 2𝑔
𝑷𝟏 −𝑷𝟐
+ (𝑍1 − 𝑍2 ) = ℎ𝐿 → ∆ℎ = −ℎ𝐿 (constant diameter)
𝜸
𝑃1 − 𝑃2 = 𝛾[(𝑍2 − 𝑍1 ) + 𝛾ℎ𝐿 ] {To find the pressure drop}.
The following equations were derived for the case of circular pipes for the case
of laminar flow.
𝟐 𝟐 𝟐 𝟐 𝟐 𝟐 ∆𝐏 ∆𝐙
𝑽(𝒓) = − (𝒓𝒐𝟒𝝁
−𝒓
) (𝜸∆𝒉
∆𝑳
) = − (𝒓𝒐𝟒𝝁
−𝒓
) (∆(𝑷+𝜸𝒁)
∆𝑳
) = − (𝒓𝒐𝟒𝝁
−𝒓
) (∆𝐋 + 𝜸 ∆𝐋)
𝑟 2
𝑉(𝑟) = 𝑉𝑚𝑎𝑥 (1 − ( ) )
𝑟 𝑜
𝒓𝟐 ∆𝒉 𝟐 𝟐 ∆𝐏 ∆𝐙
𝑽𝒎𝒂𝒙,𝒄𝒆𝒏𝒕𝒆𝒓 = 𝐕(𝐫 = 𝟎) = − ( 𝒐 ) [ 𝜸 ] = − (𝟒𝝁
𝒓𝒐
) (∆(𝑷+𝜸𝒁)
∆𝑳
𝒓𝒐
) = − (𝟒𝝁 ) (∆𝐋 + 𝜸 ∆𝐋)
𝟒𝝁 ∆𝑳
𝜋𝑟𝑜4 𝛾∆ℎ
𝑄 = −( ) ( ∆𝐿 )
8𝜇
𝑄 𝑟2 𝛾∆ℎ 𝑉𝑚𝑎𝑥
The mean velocity: 𝑉 = = −( 𝑜)( ) =
𝐴 8𝜇 ∆𝐿 2
̅
𝟑𝟐 𝝁 𝑳 𝑽 𝟔𝟒 𝑳 𝑽𝟐
𝒉𝒇 = = ; {Laminar, fully developed, steady, Newtonian flow}.
𝜸𝑫𝟐 𝑹𝒆 𝑫 𝟐𝒈
∆(𝑃+𝛾𝑍) ∆(ℎ) ∆P ∆Z
𝜏 = 2𝑟 [− ] = 2𝑟 [−𝛾 ] = −2𝑟 ( + 𝛾 ) , {laminar and turbulent}.
∆𝐿 ∆𝐿 ∆L ∆L
98
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Quiz
Answer Sheet
99
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Quiz
Q. The 6-cm diameter pipe shown contains glycerin, ρ = 1420 kg/m3 and µ=7.8 Poise, flowing at
a flow rate of 6 m3/h.
r (cm) 0 1 2 3
V (m/s)
100
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
3.9 Types of Fluid Flow Problems:
In the design and analysis of piping systems that involve the use of the Moody chart,
we usually encounter three types of problems (figure 3.12):
1. Determining the energy head loss (pressure drop) when the the pipe length and diameter
are given for a specified flow rate ( or velocity).
Use the energy equation.
Find the velocity, Reynolds number, friction factor and solve directly the energy
equation to find the energy loss or pressure difference.
2. Determining the flow rate when the pipe length and diameter are given for a specified
head loss or pressure drop.
The energy equation is used.
Substituting the known values in the energy equation will result in two unknowns (f,
and Q).
Solving the equation requires assuming value for f, substituting this value in the energy
equation to find Q or V.
Reynolds number is found using the last value of V, then using Moody diagram a new
value of f is found.
The new value of f is compared with the previous assumed value. If the difference
between both values is insignificant then we stop and consider the value of Q to be
correct, otherwise we use the new value of f in the energy equation to find Q or V.
the same process is repeated until successive values of f are not significantly different.
3. Determining the pipe diameter when the pipe length and flow rate are given for a specified
pressure drop.
The previous procedure is applied but the unknowns will be the friction factor f and
the diameter.
101
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Case A: Horizontal pipe, constant elevation (Z1=Z2), the pressure must decrease in the
downstream direction.
Case B: Inclined upward, elevation is increasing, the pressure is decreasing.
Case C: Inclined downward, eleveation is decreasing (giving extra value to the pressure).
102
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 3.13:
Oil with =900 kg/m3 and ν =0.00001 m2/s flows at Q=0.2 m3/s through 500 m of 20-cm
diameter cast iron pipe. Determine the head and the pressure drop if the pipe slopes down
at 10o in the flow direction.
Solution: (Type 1 problem to find hf or ∆P)
The head drop or loss can be given as:
𝑷𝟏 𝑽𝟐𝟏 𝑷𝟐 𝑽𝟐𝟐
( + + 𝒁𝟏 ) − ( + + 𝒁𝟐 ) = 𝒉𝒇
𝜸 𝟐𝒈 𝜸 𝟐𝒈
𝑷𝟏 − 𝑷𝟐 (𝑽𝟐𝟏 − 𝑽𝟐𝟐 )
+ (𝒁𝟏 − 𝒁𝟐 ) + = 𝒉𝒇
𝜸 𝟐𝒈
Example 3.14:
When gasoline with a kinematic viscosity of 510-7 m2/s flows in a 200-mm-diameter
smooth pipe the friction head loss is 0.43 m per 100 m. Find the flow rate.
Solution: (Type 2 problem to find Q)
For smooth pipe ( 𝑒 = 0)
The total energy loss can be found as:
ℎ𝑓 = 0.43 = 𝑓 (100/0.2) 𝑉 2 /[2(9.81)]
𝑉 = 0.1299 / 𝑓 1/2
𝑅𝑒 = 𝑉𝐷/𝜈 = 0.2 𝑉 / (510−7 ) = 400000 𝑉
103
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
let the initial assumption: 𝑓 = 0.02 ;𝑉 = 0.918𝑚/𝑠 ; 𝑅𝑒 = 3.67105 𝑓 = 0.014
(Moody diagram)
The following table represents subsequent trials:
Try f V m/s Re New f
0.02 0.918 3.67105 0.014 Try again
0.014 1.101 4.4105 0.0134 Try again
0.0134 1.12 4.48105 0.0134 accepted
𝑄 = (/4)(0.2)2 (1.12) = 0.0352 𝑚3 /𝑠 = 35.2 𝐿/𝑠
Example 3.15:
A steel pipe (e=0.065 mm) 4200 m long is to
convey oil (ν=5.210-5 m2/s) at 300 L/s from
a reservoir with surface elevation 247 m to
one with surface elevation 156 m. What pipe
size is required?
Solution: (Type 3 problem to find D)
2
0.3
( 2)
𝜋𝐷
4200
ℎ𝑓 = 𝑍1 − 𝑍2 = 247 𝑚 − 156 𝑚 = 91 𝑚 = 𝑓 4
𝐷 5 = 0.343 𝑓
𝐷 2𝑔
Tutorial 4.
Q1. Water at 20o C, having =1×10-6 m2/s, flows through a smooth plastic pipe of
16 cm diameter maintaining an average velocity of 2.8 m/s. What is the head
loss per meter of length of pipe?
Q2. Water flows through a 17 cm diameter steel riveted pipe, k = 0.3 mm.
Determine the discharge if there is a loss of 2.8 m in a km length. What is the
power required to maintain the flow? ( = 1×10-6 m2/s).
Q3. What size of steel pipe, k = 0.3 mm, should be used to carry 550 L/s of kerosene,
= 2×10-6 m2/s, with an energy loss of 2.1 m per km.
Solution of tutorial 4:
Q1.
Given:
=1×10-6 m2/s (water) Note: The slope of the pipe is not
important since the requirement is the
Pipe: {smooth plastic pipe ,
total head loss
D = 16 cm}
Flow:{V = 2.8 m/s}
Required:
ℎ𝐿
= ⋯?
𝐿
Solution:
𝑉𝐷 2.8 ×0.16
𝑅𝑒 = = = 448000
𝜈 10−6
→ Turbulent flow
Moody chart {smooth plastic pipe} →
f = 0.0133}
𝒉𝑳 𝟏 𝑽𝟐
=𝒇
𝑳 𝑫 𝟐𝒈
𝟏 𝟐.𝟖𝟐
= 𝟎. 𝟎𝟏𝟑𝟑
𝟎.𝟏𝟔 𝟐𝒈
= 𝟎. 𝟎𝟑𝟑𝟐 𝒎⁄𝒎
105
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q2.
Given:
=1×10-6 m2/s (water)
Pipe: {k = 0.3 mm , D = 17 cm}
ℎ 2.8𝑚
Flow: {Q =…? m/s, 𝐿 = }, Power =… ?
𝐿 1000𝑚
Required:
Q, Power
Solution:
𝑘 0.3 𝑚𝑚
= = 0.00176
𝐷 170𝑚𝑚
ℎ𝐿 1 𝑉2 2.8 1 𝑉2
=𝑓 = =𝑓 → 9.34 × 10−3 = 𝑓𝑉 2
𝐿 𝐷 2𝑔 1000 0.17 2𝑔
1st iteration:
Let f = 0.026
𝑉𝐷 0.599 ×0.17
→ V = 0.599 m/s → 𝑅𝑒 = = = 101890.96 → f = 0.0244
𝜈 10−6
2nd iteration:
Let f = 0.0244
→ V = 0.62 m/s → 𝑅𝑒 = 105394 → f = 0.0243 → V = 0.62 m/s , then
106
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝝅
𝑸 = 𝑽𝑨 = 𝟎. 𝟔𝟐 × (𝟎. 𝟏𝟕)𝟐 = 𝟎. 𝟎𝟏𝟒 𝒎𝟑 /𝒔
𝟒
𝟐. 𝟖𝒎 𝒘𝒂𝒕𝒕
𝑷𝒐𝒘𝒆𝒓 = 𝜸𝑸𝒉 = 𝟗𝟖𝟏𝟎 × 𝟎. 𝟎𝟏𝟒 × = 𝟑𝟖𝟔. 𝟓𝟓
𝟏𝒌𝒎 𝒌𝒎
Q3. Given:
=2×10-6 m2/s (water)
Pipe: {k = 0.3mm , D = …? cm}
ℎ 2.1𝑚
Flow:{Q = 550 l/s , 𝐿 = }
𝐿 1𝑘𝑚
Required:
pipe diameter
Solution:
𝑄 0.55 𝑚3 /𝑠 0.7
𝑉= = 𝜋⁄ 𝐷2 = (1)
𝐴 4 𝐷2
𝑘 0.3
= (2)
𝐷 𝐷
𝑉𝐷 (0.7⁄ 2 )𝐷 3.5×105
𝐷
𝑅𝑒 = = = (3)
𝜈 2×10−6 𝐷
ℎ𝑓 𝑓 𝑉2
=
𝐿 𝐷 2𝑔
2.1 𝑓 0.7 2 1
= ( 2)
1000 𝐷 𝐷 2𝑔
𝑓
= 0.084 (4)
𝐷5
1st iteration:
Let f = 0.020
→ D = 0.75 m (from equation (4))
3.5×105 𝑘 0.03𝑐𝑚
→ 𝑅𝑒 = = 466666.7 , with = = 0.0004
𝐷0.75 𝐷 0.75𝑐𝑚
107
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Class work
The tanks in the figure contains water at 20OC. If the pipe is smooth with L =4500
m and d = 4 cm:
(1) What will the flow rate in m3/h be for ∆𝒛=100 m?
(2) Assuming 5o slope of the pipe down ward, and the pipe entrance of tank (1)
is 85 m below the water surface in the tank. Find the pressure at the middle
distance of pipe (2250 m from tank 1).
𝒌𝒈
(Consider the minor losses, 𝝆 = 𝟗𝟗𝟖 𝒌𝒈/𝒎𝟑 , 𝝁 = 𝟎. 𝟎𝟎𝟏 )
𝒎.𝒔
108
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝐹 𝑢
=𝜇
𝐴 𝑦
𝑑𝑢
𝜏=𝜇
𝑑𝑦
109
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
In case of real fluids, the fluid particles will adhere to the boundary (uboundary =0) even
for small viscosity fluids “No-Slip Condition”.
𝑑𝑢
Farther away from the boundary / wall there is a velocity gradient ( ) which causes
𝑑𝑦
the fluid to exert tangential stress
As we go farther away from the
boundary, the retardation due to the
presence of viscosity is negligible.
Consider a thin flat plate held horizontal (zero slope), that is parallel to a flow in a uniform
stream of velocity U known as a “free stream velocity” or potential velocity”, Figure 4.4.
Boundary Layer can be defined as:
The region adjacent to a surface over which the velocity changes from the free stream
velocity to zero at the surface.
The development and growth of the boundary layer occurs because of the no-slip
condition at the surface where the fluid velocity must be zero (no-slip condition: the fluid
velocity at all fluid–solid boundaries is equal to that of the solid boundary. Particles close to a surface
do not move along with a flow when adhesion is stronger than cohesion. At the fluid-solid interface,
the force of attraction between the fluid particles and solid particles (Adhesive forces) is greater
than that between the fluid particles (Cohesive forces). This force imbalance brings down the fluid
velocity to zero).
The most powerful consequence of the no-slip condition is the appearance of strong
velocity gradients near the wall. This region is often very thin, and it is then called a
boundary layer
Within the boundary layer, adjacent layers of fluid are in relative motion (velocity
gradient), and because all fluids have viscosity, there will be friction between the
layers as they slide over each other. In other words, viscous shear stresses are
produced, with a magnitude given by the viscosity times the velocity gradient.
110
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
The thickness of the boundary layer was defined by Prandtl as δ: The distance from
the plate (boundary) to where 99% of the free stream velocity, U, is achieved again.
{at y = (thickness of B.L.) u = 0.99 U
The fluid particles next to the plate are affected by shear as they pass close to the
leading edge of the plate which slows them. These particles continue to be subject to shear
stress from the plate as they progress downstream continue to decelerate.
In addition, these particles because of their lower velocity retard other particles
adjacent to them farther from the plate. It is noticed downstream that:
The boundary layer becomes thicker.
The velocity gradient at the wall becomes smaller and the local shear stress is reduced
as a result.
Figure 4.4 Development of boundary layer and shear stress along thin flat plate
111
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Transition region:
The region after the critical point where the flow becomes unstable showing
characteristics between laminar and turbulent and the B.L. is unstable is known
as transition region. It is usually small region.
Figure 4.6 Laminar, transition, and turbulent boundary layers over thin flat plate
Downstream. If the flow does not separate from the boundary, the boundary
layer continues to increase to form what is known as:
112
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Turbulent Boundary Layer:
If the plate is smooth, then even though the boundary layer is turbulent, there
is a very thin layer near the boundary in which the flow is still laminar.
This layer is known as “Laminar Sub-Layer”.
Laminar Sub-Layer:
The flow is smooth and steady
The presence of the wall dampens the cross-stream mixing and turbulent
fluctuations.
𝑌
o: wall shear stress is constant. 𝑢 = o
𝑈
Experiments by Prandtl – Karman have provided that the thickness of laminar–sub
layer 𝛿̅ is :
11.6 11.6
𝛿̅ = = (4.1)
√o ⁄𝜌 𝑈∗
Where,
𝑈∗ = √o ⁄𝜌 (4.2)
is the shear velocity, and o is the average shear stress at the boundary.
Figure 4.7 Laminar-Sublayer region for the flow over a thin flat plate
114
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
4.4 Resistance of Smooth and Artificially Roughened Pipes
Recall:
𝐷 𝜕ℎ (𝑃1+𝛾𝑍1)−(𝑃2+𝛾𝑍2)
𝑜 = − 𝛾 =𝐷 (4.3)
4 𝜕𝑠 4𝐿
or
𝜕ℎ (ℎ1−ℎ2)
4𝑜 = −𝐷𝛾 = 𝐷𝛾 ,
𝜕𝑠 𝐿
𝑉2 𝐷 𝜕ℎ 1 (ℎ1−ℎ2) 1
) 4𝑜
1
Divide by ( 𝜌 𝑉2
=𝑓=− 𝛾 𝑉2
= 𝐷𝛾 𝑉2
2 𝜌 4 𝜕𝑠 𝜌 𝐿 𝜌
2 2 2
0.316
𝑓= (4.6)
𝑅𝑒0.25
or
1 𝜀 2.51
= −2𝑙𝑜𝑔 ( + ) (4.8)
√𝑓 14.8𝑅ℎ 𝑅𝑒 √𝑓
where:
- 𝐷ℎ ,is the Hydraulic diameter, (m, ft) – For fluid-filled, circular conduits, Dh = D =
inside diameter.
- 𝑅ℎ = 𝐷/4 is the Hydraulic radius, (m, ft) – For fluid-filled, circular conduits.
115
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Additional, mathematically equivalent forms of the Colebrook equation are:
1 𝑅 𝑅⁄
𝑘
− 2𝑙𝑜𝑔 = 1.74 − 2𝑙𝑜𝑔(1 + 18.57 ) (4.9)
√𝑓 𝑘 𝑅𝑒√𝑓
Where, R: radius of pipe, k: average size of pipe roughness, Re; Reynolds number, f:
friction factor.
This equation was later plotted by Moody.
Moody diagram is frequently used for finding the value of f for commercial pipes.
1 𝜀/𝐷 5.74
= −2𝑙𝑜𝑔 ( + ) (4.10)
√𝑓 3.7 𝑅𝑒0.9
or
0.25
𝑓= 2 (4.11)
𝜀/𝐷 5.74
[𝑙𝑜𝑔( + 0.9)]
3.7 𝑅𝑒
116
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Table 4.1 Roughness coefficient values for some pipes from different sources
ABSOLUTE ROUGHNESS
TYPE OF PIPE
€ (in) € (mm)
GLASS OR PLASTIC SMOOTH SMOOTH
DRAWN TUBING 0.00006 0.0015
COMMERCIAL STEEL OR WROUGHT IRON 0.0018 0.046
ASPHALTED CAST IRON 0.0048 0.12
GALVANIZED IRON 0.006 0.15
CAST IRON 0.0102 0.26
RIVETED STEEL 0.072 1.8
117
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Tutorial 5
Q1. Below are tabulated the actual velocity measurements at various points in the
vertical at a particular station in the wind tunnel:
Y cm 0.20 0.50 0.70 0.90 1.35 1.58 1.75 2.25 2.50 3.00
u m/s 4.93 6.30 6.90 8.01 9.36 10.05 10.35 10.48 10.48 10.48
Plot the velocity profile with exaggerated Y scale and determine the thickness of
(Y/δ) in the relation:
u/U=(Y/δ)m
Q2.In the following cases determine whether the pipe surface will act as smooth,
rough or in transition.
Q3.Find the value of the friction factor f for the following cases:
118
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Solution of tutorial 5:
Q1. Below are tabulated the actual velocity measurements at various points in the
vertical at a particular station in the wind tunnel:
Y cm 0.20 0.50 0.70 0.90 1.35 1.58 1.75 2.25 2.50 3.00
u m/s 4.93 6.30 6.90 8.01 9.36 10.05 10.35 10.48 10.48 10.48
Plot the velocity profile with exaggerated Y scale and determine the thickness of
(Y/δ) in the relation:
u/U=(Y/δ)m
Modified Table:
Y cm 0.20 0.50 0.70 0.90 1.35 1.58 1.75 2.25 2.50 3.00
u m/s
𝒖 𝒎
= [𝒀⁄𝜹]
𝑼
𝒖
𝒍𝒐𝒈 ( ) = 𝒎 𝒍𝒐𝒈[𝒀⁄𝜹]
𝑼
𝒚=𝒎𝒙
𝒖
𝒍𝒐𝒈(𝑼)
𝒎=
𝒍𝒐𝒈(𝒀⁄𝜹)
119
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q2. In the following cases determine whether the pipe surface will act as smooth, rough
or in transition.
Boundary is Transition
120
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q3. Find the value of the friction factor f for the following cases:
(I) By calculation, using laminar flow formula/ Blassius formula/ Colebrook-
White formula
K D
Type of k/D Re f
(mm) (cm)
pipe
Case (1): 1700
Smooth - - -
Case (2): 4800
Rough - - 0.002 Case (3): 1200
Case (4): 106
- 0.3 50 -
Using Colebrook formula
𝟏 𝑹 𝑹⁄
− 𝟐𝒍𝒐𝒈 = 𝟏. 𝟕𝟒 − 𝟐𝒍𝒐𝒈(𝟏 + 𝟏𝟖. 𝟕 𝒌 )
√𝒇 𝒌 𝑹𝒆√𝒇
Complete by assuming values for f until equalizing both sides of the equation
121
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
122
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
The energy loss within the system will be due to: the friction loss within the pipes
𝐿 𝑉2
which can be summed cumulatively (∑ 𝑓 ); the minor loss due to the pipe exit,
𝐷 2𝑔
entrance, valves, contraction and expansion. The energy equation can be applied between
A and B at the up and down streams tanks as:
Given:
𝑉1 𝐴1 = 𝑉2 𝐴2 (continuity equation) (5.3)
𝐷1 2
𝑉2 = ( ) 𝑉1 (5.4)
𝐷 2
used to eliminate 𝑉2 .
Given the values of pipe lengths and diameters known, the previous equation can be
reduced to:
𝑽𝟐𝟏
𝑯= (𝑪𝟏 + 𝑪𝟐 𝒇𝟏 + 𝑪𝟑 𝒇𝟐 ) (5.5)
𝟐𝒈
123
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Where C1, C2, and C3 are known constants.
2
𝑄2 𝐿1 𝑄2 𝐿2 𝑄2 𝑄2 1 1 𝑄2
𝐻 = 0.5 + 𝑓1 + 𝑓2 + ( − ) +
2𝑔𝐴12 𝐷1 2𝑔𝐴12 𝐷2 2𝑔𝐴22 2𝑔 𝐴12 𝐴22 2𝑔𝐴22
2
𝑄2 1 1 1 1 𝐿1 1 𝐿2 1
𝐻= {[0.5 2 + ( 2 − 2 ) + 2 ] + [ 2 ] 𝑓1 + [ ] 𝑓}
2𝑔 𝐴1 𝐴1 𝐴2 𝐴2 𝐷1 𝐴1 𝐷2 𝐴22 2
𝑸𝟐
𝑯= {𝑲 + 𝑲𝟐 𝒇𝟏 + 𝑲𝟑 𝒇𝟐 }
𝟐𝒈 𝟏
124
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 5.1: For the system of pipes in series shown, given that: The flow of water in
the system is 1000 L/s, kexit =1, Kentrance =0.5, the Losses due expansion is (V1-V2)2/2g, the
kinematic viscosity is =3×10-6 m2/s.
𝐻 = 86.97 𝑚
125
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 5.2: Three pipes, as shown, are fitted between two reservoirs. Determine
the discharge if the difference in the water surface elevations is 8 m. Plot the energy
and hydraulic grade lines.
Solution:
Energy equation can be taken between the surfaces of two tanks:
𝑉12 𝐿1 𝑉12 𝐿2 𝑉22 𝐿3 𝑉32 (𝑉1 − 𝑉2 )2 𝑉32 𝑉32
8 + 0 + 0 = 0 + 0 + 0 + 0.5 + 𝑓1 + 𝑓2 + 𝑓3 + + 0.5 +
2𝑔 𝐷1 2𝑔 𝐷2 2𝑔 𝐷3 2𝑔 2𝑔 2𝑔 2𝑔
To simplify the energy equation, The velocities can be introduced related to the discharge:
𝑄 𝑄
𝑉1 = 𝑉3 = 𝜋(0.5)2 = 5.093𝑄 , 𝑉2 = 𝜋(0.6)2 = 3.537𝑄
4 4
8 3
𝑄=√ = 0.233 𝑚 ⁄𝑠
147.206
Notes:
It can be seen how easy it is to solve such problem when the friction factor values
are considered to be fully turbulent.
The final numbers between brackets for the particular losses shows how significant
is the major losses with respect to the minor losses.
The plot of the EGL and HGL is shown in the following figure. It is important to
remember that both lines are parallel since no gradual change in diameter is there.
126
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Also , remember that the difference between these two lines is given by the velocity
head: V2/2g
127
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 5.3: Given the following data for the pipes in series (see the table), find the
discharge through the system for H= 6m.
Given:
ν = 310-6 m2/s
1 = 2 mm, 2 = 0.3 mm
D L V
Pipe /D fultimate Re f V f
(m) (m) (m/s)
1 0.6 300 0.00333 0.026 2.848 5.70105 0.0265 2.819 0.0265
2 1 240 0.0003 0.015 1.025 3.42105 0.0168 1.01484 0.0168
Solution: (The previous table will be used to show the results)
V2 can be given relative to V1 using the continuity equation: V2 = (0.6/1)2 V1 =0.36 V1 .
The first assumption for the friction factor values for both pipes are the ultimate
values.
The energy equation between A and B:
𝝅
𝑸 = 𝑨𝑽 = (𝟎. 𝟔)𝟐 × 𝟐. 𝟖𝟏𝟗 = 𝟎. 𝟕𝟗𝟕 𝒎𝟑 /𝒔
𝟒
128
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
129
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Tutorial 6.
Calculate the value of H for 1250 L/s water flow. Include minor losses in your
calculation.
Solution:
𝑉12 𝐿1 𝑉12 (𝑉1 −𝑉2 )2 𝐿2 𝑉22 𝑉22
𝑍1 + 0 + 0 = 𝑍2 + 0 + 0 + 0.5 + 𝑓1 + +𝑓2 +
2𝑔 𝐷1 2𝑔 2𝑔 𝐷2 2𝑔 2𝑔
1.25
𝑉1 = 𝜋 = 9.947𝑚/𝑠
0.42
4
1.25
𝑉2 = 𝜋 = 3.248 𝑚/𝑠
0.72
4
130
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
a) Compute the flow discharge through the system given that the difference
between the levels of the two reservoirs is H= 6 m.
b) Draw the energy grade line between the two reservoirs.
Solution:
Solving this question requires applying the energy equation EE between the surfaces
of the two tanks (points A and B) as shown in the following figure:
EE between A and B:
𝑉12 𝐿1 𝑉12 (𝑉1 −𝑉2 )2 𝐿2 𝑉22 𝑉22
𝑍1 + 0 + 0 = 𝑍2 + 0 + 0 + 0.5 + 𝑓1 + +𝑓2 +
2𝑔 𝐷1 2𝑔 2𝑔 𝐷2 2𝑔 2𝑔
𝑉12 300 𝑉12 (𝑉1 −𝑉2 )2 240 𝑉22 𝑉22
𝑍1 − 𝑍2 = 6𝑚 = 0.5 + 𝑓1 + +𝑓2 +
2𝑔 0.6 2𝑔 2𝑔 1 2𝑔 2𝑔
(𝑉1 −𝑉2 )2
6𝑚 = 25.484𝑓1 𝑉12 + 12.232𝑓2 𝑉22 + 25.484 × 10−3 𝑉12 + 0.051𝑉22 +
2𝑔
131
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
132
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q3. Determine the flow rate in the horizontal pipeline shown which discharges water
into the air. Draw the EGL and HGL, given that: The kinematic viscosity is = 1×10-6
m2/s.
Solution:
The flow system includes two horizontal pipes connected in series with the same
diameter and different types.
Pipes are conveying water into open area.
Tank (1) is feeding the system from an elevation of 25 m above the existing point
It is required to find the discharge and to draw the EGL and HGL.
The coefficients of pipe entrance and exit are Ke=0.5 and Kexit=1.
Kinematic viscosity is =1×10-6 m2/s.
133
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
134
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
It is required to plot the relation between the system discharge and the level
difference between the two tanks (H-Z).
The data of the pipe is given in the figure.
The coefficients of pipe entrance and exit are Ke=0.5 and Kexit=1.
Kinematic viscosity is =1×10-6 m2/s).
𝑉2 𝐿1 𝐿2
𝐻 − 𝑧 = 10𝑚 = (0.5 + 𝑓1 + 𝑓2 + 1.0 + 1.0)
2𝑔 𝐷 𝐷
2
𝑉 200 180
10𝑚 = (2.5 + 𝑓1 + 𝑓2 )
2𝑔 0.2 0.2
𝑉2
10𝑚 = (2.5 + 1000 𝑓1 + 900 𝑓2 )
2𝑔
135
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q = 0.000 𝑚3 /𝑠
𝑉2 𝐿1 𝐿2
𝐻+0+0=𝑍+0+0+ (0.5 + 𝑓1 + 𝑓2 + 1.0)
2𝑔 𝐷 𝐷
2
𝑉 𝐿1 𝐿2
𝐻 − 𝑍 = 10𝑚 − 𝑍 = (1.5 + 𝑓1 + 𝑓2 )
2𝑔 𝐷 𝐷
2
𝑉
10𝑚 − 𝑍 = (1.5 + 1000 𝑓1 + 900 𝑓2 )
2𝑔
𝑉2 𝐿1 𝐿2
𝐻+0+0=𝑍+0+0+ (0.5 + 𝑓1 + 𝑓2 + 1.0)
2𝑔 𝐷 𝐷
136
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝑉2 𝐿1 𝐿2
𝐻 − 𝑍 = 10𝑚 − 𝑍 = (1.5 + 𝑓1 + 𝑓2 )
2𝑔 𝐷 𝐷
2
𝑉
10𝑚 − 𝑍 = (1.5 + 1000 𝑓1 + 900 𝑓2 )
2𝑔
Q m3/s
0.07
0.06
0.05
0.04
0.03
0.02
0.01
0
0 2 4 6 8 10 12
(H-Z) (m)
Q5.
Three pipes, as shown,
are fitted between two
reservoirs. Determine
the discharge in the
system if the difference
in the water surface
elevations is 9 m. Plot
the energy grade lines
137
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Solution:
𝑉12 𝐿1 𝑉12 𝐿2 𝑉22 𝐿3 𝑉32 (𝑉1 − 𝑉2 )2 𝑉32 𝑉32
9 + 0 + 0 = 0 + 0 + 0 + 0.5 + 𝑓1 + 𝑓2 + 𝑓3 + + 0.5 × +
2𝑔 𝐷1 2𝑔 𝐷2 2𝑔 𝐷3 2𝑔 2𝑔 2𝑔 2𝑔
To simplify the energy equation, the velocities can be introduced related to the discharge:
𝑄 𝑄
𝑉1 = 𝑉3 = 𝜋(0.45)2 = 6.287𝑄 , 𝑉2 = 𝜋(0.55)2 = 4.209𝑄
4 4
9 = Q2 (91.104+35.871+111.922+2.767) = 241.15 Q2
𝟗 𝟑
𝑸=√ = 𝟎. 𝟏𝟗𝟑 𝒎 ⁄𝒔
𝟐𝟒𝟏. 𝟏𝟓
138
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Figure 5.2 Energy losses between branching pipes connecting three tanks through one junction
𝑷𝒋
A very important issue here is the value of the energy at the junction, “𝑯𝒋 = + 𝒁𝒋 ”.
𝜸
Which is the sum of the pressure head and elevation head at the junction. This junction and
its energy elevation head is common between the three tanks!! This means that the flow in
the three pipes is taking place between two different values of energy; one of these values is
common for the three pipes, which is the “head at the junction, “𝑯𝒋 ". The value of “𝑯𝒋 " in
the problem will be the same for the three pipes. So, the Darcy-Weisbach-equation is applied
for three pipes, such that the head loss in the equation is the difference between each tank
level ( 𝑍1 or 𝑍2 or 𝑍3 ) and the junction head “𝐻𝑗 ", with a direction of flow in the pipe “in”
for the highest energy tank, “out” for the lowest tank and “in” or “out” for the intermediate
tank as shown in figure 5.3.
139
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Figure 5.3 Branching pipes with two different possibilities for system feeding scenario.
The conservation of mass principle will have an important role in this problem. Since
we have flow toward and out of the junction, then it should be easily concluded that the total
flow into the junction should equal to the total flow out of the junction.
Problem definition:
The flow through each pipe is required when the reservoir (tank) elevations are given.
The sizes and types of the pipes; and the fluid properties are known.
- The flow into the junction J must be equal to the flow out of the junction.
- Flow will be out of the reservoir (reservoirs) with peizometric head (𝑍) more
than “𝐻𝑗 " toward the junction “in”; on the other hand the flow will be out of
the junction ,”out”, toward the reservoir (reservoirs) with (𝑍) lower than “𝐻𝑗 "
- A value of the head at the junction “𝐻𝑗 " is assumed. It should be more than the
minimum tank elevation and less than the maximum tank elevation.
140
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
- The values of {𝑄1 , 𝑄2 , and 𝑄3 } are calculated by applying the Darcy-Weisbach
equation between the tank levels {𝑍1 , 𝑍2 , and 𝑍3 } and the junction energy
head “𝐻𝑗 ":
𝐿1 𝑉12 𝐿1 𝑄12
𝐻1 − 𝐻𝐽 = 𝑓1 = 𝑓1 (5.7)
𝐷1 2𝑔 𝐷1 2𝑔𝐴21
𝐿2 𝑉22 𝐿2 𝑄22
𝐻𝐽 − 𝐻2 𝑜𝑟 𝐻2 − 𝐻𝐽 = 𝑓2 = 𝑓2 (5.8)
𝐷2 2𝑔 𝐷2 2𝑔𝐴22
𝐿3 𝑉32 𝐿3 𝑄32
𝐻𝑗 − 𝐻3 = 𝑓3 = 𝑓3 (5.9)
𝐷3 2𝑔 𝐷3 2𝑔𝐴23
If the flow into the junction is greater than the flow out of the junction(∆𝑄 > 0),
a higher water elevation head HJ is assumed which will reduce the inflow and
increase the outflow.
Quiz
Can we have flow from the junction toward Tank 3 upward?
How can we express the energies at Tank1, 2, and 3?
What type of energy do we have at the junction J?
141
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 5.4: Find the discharges for water at 60 oF given the following pipe data and
reservoir elevations:
L1=10,000 ft, D1=3ft, 1/D1=0.0002 , Z1=100ft
L2=2000 ft, D2=1.5 ft, 2/D2=0.002 , Z2=60ft
L3=4000 ft, D2=2ft, 3/D3= 0.001 , Z3=30ft
10,000 𝑉12
100 − 65 = 35 = 𝑓1 , 𝑏𝑦 𝑡𝑟𝑖𝑎𝑙𝑠 𝑓1 = 0.014 → 𝑉1 = 6.95 𝑓𝑡/ sec → 𝑄1 = 49.1 𝑐𝑓𝑠 “in”
3 2𝑔
2,000 𝑉22
65 − 60 = 5 = 𝑓2 , 𝑏𝑦 𝑡𝑟𝑖𝑎𝑙𝑠 𝑓2 = 0.024 → 𝑉2 = 3.17 𝑓𝑡/ sec → 𝑄2 = 5.6 𝑐𝑓𝑠 “out”
1.5 2𝑔
4,000 𝑉32
65 − 30 = 35 = 𝑓3 , 𝑏𝑦 𝑡𝑟𝑖𝑎𝑙𝑠 𝑓3 = 0.020 → 𝑉3 = 7.51 𝑓𝑡/ sec → 𝑄3 = 23.6 𝑐𝑓𝑠 “out”
2 2𝑔
So the inflow is greater than the outflow by: ∆𝑄 = 49.1 − 5.60 − 23.6 = 19.9 𝑐𝑓𝑠
It is required to increase the out flow and decrease the inflow, this can be done by increasing
the value of Hj which will decrease the headloss between the feeding reservoirs and the
junction (Hj). This will decrease the inflow, and will increase the energy difference between HJ
and the the reservoir being fed which will increase the outflow.
So, let’s assume [𝐻𝑗 = 𝑍𝑗 + 𝑃𝑗 / = 80 𝑓𝑡], Then,
10,000 𝑉12
100 − 80 = 20 = 𝑓1 , 𝑏𝑦 𝑡𝑟𝑖𝑎𝑙𝑠 𝑓1 = 0.015 → 𝑉1 = 5.07 𝑓𝑡/ sec → 𝑄1 = 35.8 𝑐𝑓𝑠 “in”
3 2𝑔
2,000 𝑉22
80 − 60 = 20 = 𝑓2 , 𝑏𝑦 𝑡𝑟𝑖𝑎𝑙𝑠 𝑓2 = 0.024 → 𝑉2 = 6.35 𝑓𝑡/ sec → 𝑄2 = 11.33 𝑐𝑓𝑠 “out”
1.5 2𝑔
4,000 𝑉32
80 − 30 = 50 = 𝑓3 , 𝑏𝑦 𝑡𝑟𝑖𝑎𝑙𝑠 𝑓3 = 0.020 → 𝑉3 = 8.98 𝑓𝑡/ sec → 𝑄3 = 28.2 𝑐𝑓𝑠 “out”
2 2𝑔
⇒So the inflow is less than the outflow by: 3.73 cfs. (∆𝑄 = 35.8 − 11.33 − 28.2 = −3.73𝑐𝑓𝑠)
⇒Then HJ= 77.6 ft. This value can be used to find the flows in the pipes:
142
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
10,000 𝑉12
100 − 77.6 = 22.4 = 𝑓1 , 𝑏𝑦 𝑡𝑟𝑖𝑎𝑙𝑠 𝑓1 = 0.015 → 𝑉1 = 5.36 𝑓𝑡/ sec → 𝑄1 = 37.83 𝑐𝑓𝑠 “in”
3 2𝑔
2,000 𝑉22
77.6 − 60 = 17.6 = 𝑓2 , 𝑏𝑦 𝑡𝑟𝑖𝑎𝑙𝑠 𝑓2 = 0.024 → 𝑉2 = 5.96 𝑓𝑡/ sec → 𝑄2 = 10.65 𝑐𝑓𝑠 “out”
1.5 2𝑔
4,000 𝑉32
80 − 77.6 = 47.6 = 𝑓3 , 𝑏𝑦 𝑡𝑟𝑖𝑎𝑙𝑠 𝑓3 = 0.020 → 𝑉3 = 8.75 𝑓𝑡/ sec → 𝑄3 = 27.50 𝑐𝑓𝑠 “out”
2 2𝑔
The out flow is 0.32 cfs greater. By extrapolating from the last two HJ values, HJ=77.4ft and
Q1=38.05 cfs, Q2=10.60 cfs, and Q3=27.45 cfs.
143
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 5.5: the three reservoirs shown are connected by pipe to a common junction,
Lengths and diameters of the pipes are as shown on the figure. All pipes are made of
concrete with e=0.6 mm. Determine the discharge in each pipe if the water temperature is
20 oC, (assume minor losses are negligible).
Solution:
Assuming HJ =110 m, and applying Darcy-Weisbach equation between each reservoir and
the junction:
Pipe 1:
1000 𝑉12 𝑒 0.6
120 − 110 = 𝑓1 , 𝑤𝑖𝑡ℎ 𝐷 = 300 = 0.002, 𝑙𝑒𝑡𝑓1 = 0.024(𝑢𝑙𝑡) → 𝑉1 = 1.57 𝑚/ s → 𝑅𝑒 =
0.3 2𝑔
4.71 × 105 → 𝑅𝑒𝑐ℎ𝑒𝑐𝑘𝑖𝑛𝑔 𝑡ℎ𝑒 𝑀𝑜𝑜𝑑𝑦 𝑑𝑖𝑎𝑔𝑟𝑎𝑚 → 𝑓1 = 0.024 𝑖𝑠 𝑐𝑜𝑟𝑟𝑒𝑐𝑡 → 𝑄1 = 0.111𝑚3 /𝑠
“IN”
Pipe 2 :
4000 𝑉22 𝑒 0.6
110 − 100 = 𝑓2 , 𝑤𝑖𝑡ℎ 𝐷 = 500 = 0.0012, 𝑙𝑒𝑡𝑓2 = 0.0205(𝑢𝑙𝑡) → 𝑉1 = 1.09 𝑚/ s → 𝑅𝑒 =
0.5 2𝑔
5.45 × 105 → 𝑅𝑒𝑐ℎ𝑒𝑐𝑘𝑖𝑛𝑔 𝑡ℎ𝑒 𝑀𝑜𝑜𝑑𝑦 𝑑𝑖𝑎𝑔𝑟𝑎𝑚 → 𝑓2 = 0.021 𝑖𝑠 𝑐𝑜𝑟𝑟𝑒𝑐𝑡 → 𝑄2 = 0.212𝑚3 /𝑠
“OUT”
Pipe 3:
2000 𝑉12 𝑒 0.6
110 − 80 = 𝑓3 , 𝑤𝑖𝑡ℎ 𝐷 = 400 = 0.0015, 𝑙𝑒𝑡𝑓3 = 0.022(𝑢𝑙𝑡) → 𝑉1 = 2.31 𝑚/ s → 𝑅𝑒 =
0.4 2𝑔
9.24 × 105 → 𝑅𝑒𝑐ℎ𝑒𝑐𝑘𝑖𝑛𝑔 𝑡ℎ𝑒 𝑀𝑜𝑜𝑑𝑦 𝑑𝑖𝑎𝑔𝑟𝑎𝑚 → 𝑓3 = 0.022 𝑖𝑠 𝑐𝑜𝑟𝑟𝑒𝑐𝑡 → 𝑄3 = 0.29 𝑚3 /𝑠
“OUT”
∑ 𝑄𝑖𝑛 − ∑𝑄𝑜𝑢𝑡 = 0.111 – (0.212 + 0.29 ) = −0.391
144
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
A graph may be constructed of HJ plotted against the corresponding values of ∆Q. The curve
intersects the vertical axes HJ (∆Q=0) at HJ = 97.3 m, which is used to compute the
discharges in pipes as: Q1 = 0.161 m3/s, Q2 = 0.065 m3/s, Q3 = 0.231 m3/s.
145
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Given ultimate friction factor values (constant f for each pipe), a relation can be simply
found between the discharge and the head at the junction HJ for each pipe using Darcy-
Weisbach equation. These relations can be used easily to find the discharge in each pipe
without the need for Moody diagram.
130 𝑄12
Pipe (1): 26 − 𝐻𝐽 = 0.02 = 2829.034 𝑄12
0.15 2𝑔(𝜋×0.152 )2
4
65 𝑄22
Pipe (2): (20 − 𝐻𝐽 ) 𝑜𝑟 (𝐻𝐽 − 20) = 0.02 𝜋 2 = 10741.491 𝑄22
0.1 2𝑔( ×0.12 )
4
200 𝑄32
Pipe (3): 𝐻𝐽 − 0 = 0.02 = 1032.836 𝑄32
0.2 2𝑔(𝜋×0.22 )2
4
Let HJ = 20 m
Q1 = 0.046 m3/s “INFLOW”, Q2 = 0.00 , Q3 = 0.139 m3/s “OUTFLOW”
∆Q = 0.046 – 0.139= - 0.093 m3/s. “Need to decrease HJ”.
Hj (m) Q1 Q2(m3/s) Q3(m3/s) Qin - Qout %ERROR
20 (m3/s)
0.0460529 0 0.1391552 -0.093102307 Decrease HJ
10 0.075204 0.0305118 0.0983976 0.007318229 Decrease HJ 23.98492
11 0.072816 0.028946 0.1032002 -0.001438247 Decrease HJ -4.96872
10.5 0.0740196 0.0297392 0.1008275 0.002931329 Decrease HJ 9.856781
10.6 0.0737804 0.0295823 0.1013065 0.002056241 Decrease HJ 6.95092
10.7 0.0735405 0.0294245 0.1017832 0.001181791 Accepted 4.01635
10.8 0.0732998 0.0292659 0.1022577 0.000307928 Accepted 1.052176
10.9 0.0730583 0.0291064 0.1027301 -0.0005654 Accepted -1.94253
10.85 0.0731791 0.0291862 0.1024942 -0.000128799 Accepted -0.4413
10.84 0.0732033 0.0292022 0.1024469 -4.14641E-05 Accepted -0.14199
10.83 0.0732274 0.0292181 0.1023997 4.58762E-05 Accepted 0.157013
10.835 0.0732153 0.0292102 0.1024233 2.20542E-06 Accepted 0.00755
10.8354 0.0732144 0.0292095 0.1024252 -1.28819E-06 Accepted -0.00441
10.8353 0.0732146 0.0292097 0.1024247 -4.14784E-07 Accepted -0.00142
10.8353 0.0732147 0.0292098 0.1024245 2.19169E-08 Accepted 7.5E-05
From the following table and figure, it can be seen that the values of discharges and Hj are :
The values of P1 and P2 are the values that will raise the total energy in the tanks (2 and 3) into
the same energy as tank (1) such that there will be no difference in energy between the two sides
of each pipe:
𝑃2 9.81𝑘𝑁
Tank (2) : 420 𝑚 + = 775 → 𝑃2 = (775 − 420)𝑚 × = 3482.55 𝑘𝑃𝑎𝑠𝑐𝑎𝑙
𝛾 𝑚3
𝑃3 9.81𝑘𝑁
Tank (3) : 130 𝑚 + = 775 → 𝑃3 = (775 − 130)𝑚 × = 6327.45 𝑘𝑃𝑎𝑠𝑐𝑎𝑙
𝛾 𝑚3
147
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
148
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Tutorial 7
Q1. Determine the flow rates in the three RIVETED STEEL pipes in the figure shown.
(Assume complete turbulence, rough pipes).
149
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q2. Find the flow rate distribution in the CAST IRON pipes in the figure shown. Given
the following data:
Pipe AJ:
1850 m long, 1.25 m diameter.
Pipe BJ:
650 m long, 0.95 m diameter.
Pipe CJ:
2470 m long, 1.6 m diameter.
Q3. Water is pumped from a river to two reservoirs A and B. The water surface in
reservoir A is at the same height as the river whilst that in reservoir B is 20 ft higher.
Pumping from the river takes place by means of a centrifugal pump, the equation
relating flow Q (in cubic ft/sec) and H ft at a constant speed being given by:
{𝐻 = 75 − 10 𝑄2 } . From the river to a junction J a common pipe is used of 8 in
diameter and 500 ft long. The branch J to the reservoir A is 5 in in diameter and 200 ft
long. The branch from J to reservoir B is 6 in in diameter and 200 ft long.
Neglecting all losses other than pipe friction, calculate the discharge to A and B.
Take f as 0.028 throughout.
150
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Assumed Hj Q-
(ft) Q QA QB (QA+QB)
21 2.054794521 1.382235 0.57006 0.1024994
21.1 2.052891051 1.3854228 0.589002 0.0784667
21.2 2.050985816 1.3886033 0.607352 0.05503
21.3 2.049078808 1.3917766 0.625165 0.0321372
21.4 2.047170025 1.3949426 0.642484 0.0097435
21.41 2.046979049 1.3952588 0.64419 0.00753
21.415 2.056979049 1.4052588 0.65419 0.01753
21.42 2.066979049 1.4152588 0.66419 0.02753
21.425 2.076979049 1.4252588 0.67419 0.03753
151
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
5.3 Siphons
Siphon: A fluid flow device commonly used to cause the fluid to flow from one
container in an upward direction over on obstacle (without the aid of a pump) and then flow
into a second lower container.
A siphon consists of a U-tube with one end submerged below the free surface of the
liquid in the container.
The center portion of the U-tube rises above the level of the liquid surface.
The free end of the U-tube lies outside the container at an elevation below the liquid surface
in the container.
Three important points within the
siphon system to perform the analysis
(figure 5.4):
Point 1:
Free liquid surface in the container.
Point 2:
lies in the U-tube at its highest
elevation
Point 3:
In the U-tube at its lowest elevation.
152
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
5.3.2 Analysis to determine the siphon flow rate output and the maximum
allowable head above the liquid:
To determine the volume flow rate, apply Bernoulli’s equation between points (1)
and (3):
(V1 and V3 are velocities at points (1) and (3); V=V3 =velocity in the pipe).
𝑉2 𝐿1−3 𝑉 2 𝑉2 𝑉2
𝑍1 + 0 + 0 = 𝑍3 + 0 + +𝑓 + 0.5 +1× (5.12)
2𝑔 𝐷 2𝑔 2𝑔 2𝑔
𝑉2 𝐿1−3
ℎ= (2.5 + 𝑓 ) (5.13)
2𝑔 𝐷
𝐿1−3
Let 𝐶 = 2.5 + 𝑓
𝐷
1
𝑉 = √2𝑔ℎ × (5.14)
𝐶
𝑄 = 𝐴𝑝𝑖𝑝𝑒 𝑉 (5.15)
To find the maximum allowable pipe elevation above the liquid free surface:
Apply Bernoulli’s equation between points (1) and (2).
𝑃2 𝑉2 𝐿1−2 𝑉 2 𝑉22
0+0+0=𝐻+ + +𝑓 + 0.5 (5.17)
𝛾 2𝑔 𝐷 2𝑔 2𝑔
𝑃2 𝑉2 𝐿1−2
0=𝐻+ + (1.5 + 𝑓 ) (5.18)
𝛾 2𝑔 𝐷
𝑃2
Let ( = −8)
𝛾
𝑉2 𝐿1−2
𝐻 =8− (1.5 + 𝑓 ) , considering the losses. (5.19)
2𝑔 𝐷
𝑉2
𝐻 =8− , neglecting the losses. (5.20)
2𝑔
153
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 5.8: a pipe 30 cm in diameter and
1200 m long connects two reservoirs, one
being 80 m below the other. The pipe has to
cross a ridge, either fully or partially, whose
summit is 10 m above the upper reservoir, as
shown. How high can the pipe go above the
top of the upper reservoir in order that the
pressure at the apex of the siphon does not
fall below 8 m below the atmosphere? The
length of the pipe from the upper reservoir
to the summit of the pipe is 40 m. Take
f=0.018 and consider the minor losses. Find
the flow rate to the lower reservoir. Figure of example 5.8
𝑉2 𝐿1−3 𝑉 2 𝑉2 𝑉2
𝑍1 + 0 + 0 = 𝑍3 + 0 + +𝑓 + 0.5 +1×
2𝑔 𝐷 2𝑔 2𝑔 2𝑔
𝑉2 1200𝑚 1200
𝑍1 − 𝑍3 = 80𝑚 = (2.5 + 0.018 ) , (2.5 + 0.018 = 74.5)
2𝑔 0.3 0.3
1 1
𝑉 = √2𝑔ℎ × = √2 × 9.81 × 80 × = 4.59 𝑚/𝑠
𝐶 74.5
𝑃2
Let ( = −8)
𝛾
154
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 5.8: A siphon consisting of a 3 cm diameter tube is used to drain water from a
tank. The outlet end of the tube is 2 m below the water surface in the tank. Neglecting
friction, calculate the discharge. If the peak point of the siphon is 1.4 m above the water
surface in the tank, estimate the pressure of fluid at the point of siphon.
Data:
Diameter of tube = 3 cm = 0.03 m
Formula needed to solve the question:
Bernoulli's equation for frictionless flow:
𝑃1 𝑉12 𝑃2 𝑉22
+ + 𝑍1 = + + 𝑍2
𝛾 2𝑔 𝛾 2𝑔
𝐻 = 1.4 𝑚
ℎ =2𝑚
Solution:
Applying Bernoulli's equation for the points 1 and 3:
P1 = 0 N/m2(g) , P3 = 0 N/m2(g) , Z1 = 0 m , Z3 = -2 m
Since the rate of fall of liquid level in the tank is almost negligible, V1 = 0 m/sec.
Therefore,
0 + 0 + 0 = 0 + (𝑉32 / 2𝑔) – 2
𝑃2 / 𝛾 = −3.4 m
Let the energies of the water at points A and B in figure 5.5 be hA and hB .
Water will be divided after A into the three differenet pipes in the parallel system.
All the water will meet again at B.
For any part of the water flowing through any of the pipes in the parallel system, the
energy equation can be given between A and B in its path (pipe 1, 2, or 3) as :
𝑃𝐴 𝑉𝐴2 𝑃𝐵 𝑉𝐵2
+ + 𝑍𝐴 = + + 𝑍𝐵 + ℎ𝐿𝐴−𝐵 (5.21)
𝛾 2𝑔 𝛾 2𝑔
156
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Note that the value of hL is the loss of energy in one the parallel pipes between A and
B .This means that water regardless which pipe it will flow in will have the same
amount of headloss: (ℎ𝐿𝑃𝑖𝑝𝑒 1 = ℎ𝐿𝑃𝑖𝑝𝑒 2 = ℎ𝐿𝑃𝑖𝑝𝑒 3 = ⋯ 𝑒𝑡𝑐)
𝑄 = 𝑄1 + 𝑄2 + 𝑄3
For pipes in series, the same fluid flows through all the pipes, and the headlosses are
comulative, while in parallel pipes water will have different routes with equal
headlosses , the discharges will be comulative.
ℎ 𝑇𝑜𝑡𝑎𝑙 = ℎ𝐿𝑃𝑖𝑝𝑒 1 + ℎ𝐿𝑃𝑖𝑝𝑒 2 + ℎ𝐿𝑃𝑖𝑝𝑒 3 , 𝑄 = 𝑄1 + 𝑄2 + 𝑄3 ,
From the previous discussion, the following conditions can be concluded for the flow
in parallel pipes:
Energy losses in all the pipes will be equal. Since the flow through all the pipes starts
from the same point or tank and ends at the same point or tank which means that they
start from the same equal value of energy and end at another equal value of energy,
then the headlosses in all pipes are equal.
For the case where the pipes at the beginning or at the end are connected into a tank,
note that the energy at all the points and elevations within the tank is the same. which
means that all these pipes connected to the tank have the same head at the points of
contact.
Energy @ 1 = Z1
Energy @ 2 = Z2 + (P2/)=Z2+ (Z1-Z2)/ = Z1
Energy @ 3 = Z3 +(P3/)=Z3+(Z1-Z3)/ = Z1
157
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Assuming the pipes
are connected to the same
upstream and downstream
points, the principle of
continuity applies:
Qin =Qout = Q1 + Q2 + Q3
Figure 5.7 Parallel pipes system
Also, the energy equation can be applied between points A and B. This equation can
take three different scenarios:
𝑃𝐴 𝑉𝐴2 𝑃𝐵 𝑉𝐵2
+ 𝑍𝐴 + = + 𝑍𝐵 + + ℎ𝐿−𝑝𝑖𝑝𝑒1 (5.22)
𝛾 2𝑔 𝛾 𝛾
𝑃𝐴 𝑉𝐴2 𝑃𝐵 𝑉𝐵2
+ 𝑍𝐴 + = + 𝑍𝐵 + + ℎ𝐿−𝑝𝑖𝑝𝑒2 (5.23)
𝛾 2𝑔 𝛾 𝛾
𝑃𝐴 𝑉𝐴2 𝑃𝐵 𝑉𝐵2
+ 𝑍𝐴 + = + 𝑍𝐵 + + ℎ𝐿−𝑝𝑖𝑝𝑒3 (5.24)
𝛾 2𝑔 𝛾 𝛾
Recall that the pipes might have different diameters, lengths, and types (roughness).
These parameters are forming the losses terms in equations 5.22, 5,23 and 5.24
𝐿 𝑉2
(ℎ𝐿−𝑝𝑖𝑝𝑒1 = 𝑓1 𝐷1 2𝑔
1
) . But for equations 5.22, 5,23 and 5.24 to be valid, the losses within
1
the pipes must be equal. So, the flow in the three pipes will adjust it self ( distribute it self in
ratios) so that the values of Q1, Q2, and Q3 will achieve:
All the fluid elements will pass at A and B, and will have the same heads at these points (A
and B). This can be real only if hL is the same in all flow routes.
So, 𝑄1 = 0.3667𝑄2
0.3667𝑄2 + 𝑄2 = 4.2 𝑚3 /𝑠 ,
which leads to: 𝑄1 = 1.125 𝑚3 /𝑠, and 𝑄2 = 3.07 𝑚3 /𝑠
Solution:
The relative roughness for all the pipes given that their type is commercial steel and their
diameter is 0.3 m will be:
160
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 5.11:
For a total flow of 12 cfs, determine
the flow through each pipe and the
pressure head at B. Given: =2.0
slug/ft3, ν = 0.00003 ft2/s, PA=80 psi,
ZA = 100 ft, ZB = 80 ft.
Solution:
Pipe Length(ft) Diameter (ft) Q' (m3/s) h'f (ft) Q (ft3/s) hf (ft)
Pipe1 3000 1 ft 0.001
Pipe2 2000 8 in 0.0001
Pipe3 4000 16 in 0.0008
The pressure head at B =
Pipe 1: Assume: 𝑄1′ = 3 𝑐𝑓𝑠;
3.82×1 𝜖
then 𝑉1′ = 3.82 𝑐𝑓𝑠, 𝑅𝑒1′ = = 127,000 , 1 = 0.001, 𝑓1′ = 0.022 (Moody), and
0.00003 𝐷1
′
3000 3.822
ℎ𝑓1 = 0.022 × = 14.97 𝑓𝑡
1 2 × 32.2
′ ′ 2000 𝑉2′2
Pipe 2: ℎ𝑓2 = ℎ𝑓1 = 14.97 = 𝑓2′ ×
8/12 2×32.2
𝜖2 4.11×8/12
= 0.00015, assume 𝑓2′ = 0.020; then 𝑉2′ = 4.11 𝑓𝑡/𝑠 , 𝑅𝑒2′ = = 89,000
𝐷2 0.00003
, 𝑓2′ = 0.019 (Moody), 𝑉2′ = 4.11 𝑓𝑡/𝑠, 𝑄2′ = 1.44 𝑐𝑓𝑠.
′ ′ 4000 𝑉3′2
Pipe 3: ℎ𝑓3 = ℎ𝑓1 = 14.97 = 𝑓3′ ×
16/12 2×32.2
𝜖3 4.01×16/12
= 0.0006, assume 𝑓3′ = 0.020; then 𝑉3′ = 4.01 𝑓𝑡/𝑠 , 𝑅𝑒3′ = = 178,000
𝐷3 0.00003
, 𝑓3′ = 0.020 (Moody) , 𝑄3′ = 5.60 𝑐𝑓𝑠.
→ The total discharge for the assumed conditions is: 𝑄𝑇′ = 3.00 + 1.44 + 5.60 = 10.04 𝑐𝑓𝑠
Correction can be made now to find the discharge in each pipe:
𝑄′ 3.00
𝑄1 = 𝑄1′ × 𝑄 = 10.04 × 12 = 3.58 𝑐𝑓𝑠 𝑉1 = 4.56𝑓𝑝𝑠 ,𝑅𝑒1 = 152,000, 𝑓1 = 0.021, ℎ𝑓1 = 20.4𝑓𝑡
𝑇
𝑄2′ 1.44
𝑄2 = 𝑄′ × 𝑄 = 10.04 × 12 = 1.721 𝑐𝑓𝑠𝑉2 = 4.93𝑓𝑝𝑠, 𝑅𝑒2 = 109,200, 𝑓2 = 0.019, ℎ𝑓2 = 21.6 𝑓𝑡
𝑇
𝑄3′ 5.60
𝑄3 = 𝑄′ × 𝑄 = 10.04 × 12 = 6.70 𝑐𝑓𝑠𝑉3 = 4.80𝑓𝑝𝑠,𝑅𝑒3 = 213,000, 𝑓3 = 0.019, ℎ𝑓2 = 20.4 𝑓𝑡
𝑇
𝑃𝐴 𝑃𝐵 𝑃𝐵 80×144
To find 𝑷𝑩 : + 𝑍𝐴 = + 𝑍𝐵 + ℎ𝑓 = + 100 − 80 − 20.8 = 183.5 𝑓𝑡.
𝛾 𝛾 𝛾 62.4
161
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Solution:
𝐿1 𝑉12 𝐿2 𝑉22 𝐿3 𝑉32
1st trial: 𝑓1 = 𝑓2 = 𝑓3
𝐷1 2𝑔 𝐷2 2𝑔 𝐷3 2𝑔
assume the friction factor f’s values are the same (cancel each other).
Substituting for vales of L and D: 4500𝑉12 = 3750𝑉22 = 1333𝑉32
So, 𝑉2 = 1.095𝑉1 , 𝑉3 = 1.837𝑉1
Using the continuity equatiuon: 𝑄 = 𝑄1 + 𝑄2 + 𝑄3
𝜋
𝑄 = 0.007𝑚3 /𝑠 = [0.022 𝑉1 + 0.042 (1.095𝑉1 ) + 0.062 (1.837𝑉1 )] = 0.00689𝑉1
4
𝑉1 𝐷1
𝑉1 = 1.017 𝑚/𝑠 ; 𝑅𝑒1 = = 5351 , 𝑅𝑒2 = 11723 , 𝑅𝑒3 = 29491 (Turbulent flow).
𝜈
2nd trial: using Moody diagram 𝑓1 = 0.044, 𝑓2 = 0.035, 𝑓3 = 0.029
𝐿1 𝑉12 𝐿2 𝑉22 𝐿3 𝑉32
Using the equation: 𝑓1 = 𝑓2 = 𝑓3
𝐷1 2𝑔 𝐷2 2𝑔 𝐷3 2𝑔
162
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Tutorial 8
Q1. For the given flow piping system, where the total flow is 1.4 m3/s, find:
I- The Gauge Pressure of the water in tank (1), and
II- The flow discharge distribution in each pipe.
III-What should be the diameter of a single pipe from B to C that can replace the 3-
parallel pipes keeping the values of the total system capacity and the
pressure at tank (1) the same?
All the pipes are made of GALVANIZED IRON.
Assume complete turbulence, rough pipes and neglect all the minor losses.
Q2. The total discharge of water in the following parallel pipe system is 0.8 m3/s.
Determine the flow rate in each pipe given that the diameter of each pipe is 0.5 m
and the length of pipes 1 and 3 is 15 m while the length of pipe 2 is 10 m. The pipes
are made of commercial steel and the kinematic viscosity is = 3×10-6 m2/s.
163
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q3. Three pipes A, B, and C are interconnected as shown in the figure. The pipes
characteristics are as follows:
Find:
(1) The rate at which water will flow in each pipe.
(2) The pressure at point P.
Solution of Q1.
Length Diameter Area ε/D 𝒇 𝒇𝑫𝑳 𝟐𝒈𝟏𝑨𝟐
Pipe (m) (m)
(m2) Q' (m3/s) h'f (m) Q (m3/s) hf (m)
AB 280 0.4 0.126 0.00038 0.0157 35.507444 XXXXX XXXXXX 1.4 69.5945902
CD 170 0.5 0.196 0.0003 0.0149 6.7041983 XXXXX XXXXXX 1.4 13.1402287
Pipe1 145 0.35 0.096 0.00043 0.0161 36.7633721 0.25 2.29771076 0.840594328 25.9769555
Pipe2 135 0.25 0.049 0.0006 0.0174 198.950327 0.10746706 2.29771076 0.36134479 25.9769555
Pipe3 140 0.2 0.031 0.00075 0.0183 662.202522 0.05890501 2.29771076 0.198060883 25.9769555
The pressure of water in Tank (1) = ................... Pascal
Diameter of the replacement pipe = ....................m
Pipe (1):
Assume Q'1 =0.25 m3/s
So, ℎ́𝑓1 = 𝑓 𝐷𝐿 2𝑔1𝐴2 × 𝑄́12 = 2.297710𝑚 (table)
Pipe (2):
ℎ́𝑓2 = ℎ́𝑓1 = 2.297710𝑚 = 198.95032𝑄́22 → 𝑄́2 = 0.1074670 𝑚3 /𝑠
Pipe (3):
ℎ́𝑓3 = ℎ́𝑓1 = 2.297710𝑚 = 662.2𝑄́32 → 𝑄́3 = 0.0589 𝑚3 /𝑠
164
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Correction:
𝑄1′ 0.25 𝑚3
𝑄1 = ×𝑄 = × 1.4 = 0.84059 → ℎ𝑓1 =
𝑄𝑇′ (0.25+0.10746+0.0589) 𝑠
25.976 𝑚
Q2 and Q3 are found the same way (See the table)
Find:
(1) The rate at which water will flow in each pipe.
(2) The pressure at point P.
165
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Solution:
Applying energy equation between (1) and (2):
𝑃1 𝑉12 𝑃2 𝑉22
+ + 𝑍1 = + + 𝑍2 + ℎ𝐿 (1−2) (1)
𝛾 2𝑔 𝛾 2𝑔
The losses from the tank to the exit at (2) can be given as the losses through pipe
(A) or (B) plus the losses in pipe (C)
ℎ𝐿 (1−2) = ℎ𝐿 (𝐴) + ℎ𝐿 (𝐶)
𝐿𝐴 𝑉𝐴2 𝐿𝐶 𝑉𝐶2
ℎ𝐿 (1−2) = 𝑓 +𝑓
𝐷𝐴 2𝑔 𝐷𝐶 2𝑔
2000 𝑉𝐴 2
4000 𝑉𝐶2
ℎ𝐿 (1−2) = 0.02 + 0.024
6/12 2(32.2) 8/12 2(32.2)
ℎ𝐿 (1−2) = 1.242𝑉𝐴 + 2.236𝑉𝐶2
2
→ the energy equation between (1) and (2) can be given as:
𝑉𝐶2
0 + 0 + 200 = 0 + + 50 + ℎ𝐿 (1−2)
2(32.2)
𝑉𝐶2
150 = + 1.242𝑉𝐴2 + 2.236𝑉𝐶2
2(32.2)
150 = 1.242𝑉𝐴2 + 2.252𝑉𝐶2 (2)
The continuity equation in the parallel pipes system can be given as:
𝑄𝐴 + 𝑄𝐵 = 𝑄𝐶
The losses in the parallel pipes (A) and (B) are equal:
ℎ𝑓(𝐴) = ℎ𝑓(𝐵)
𝑉𝐵 = 0.7217 𝑉𝐴 (4)
166
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝑉𝐴 = 1.346𝑉𝐶 (5)
𝑉𝐶 = 5.772 𝑓𝑡/𝑠
8 2
𝜋( )
𝑄𝐶 = 𝐴𝐶 𝑉𝐶 = (5.772) = 2.01 𝑓𝑡 3 /𝑠
12
4
𝑉𝐴 = 1.346𝑉𝐶 = 𝑉𝐴 = 1.346(5.772) = 7.769 𝑓𝑡/𝑠
6 2
𝜋( )
𝑄𝐴 = 𝐴𝐴 𝑉𝐴 = 12
(7.769) = 1.53 𝑓𝑡 3 /𝑠
4
To find the pressure at point P, Apply the energy equation between P and (2):
𝑃𝑃 = 0.0135 𝐼𝑏/𝑓𝑡 2
167
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
QUIZ
Q1. Three reservoirs (1), (2), and (3) are connected respectively with pipes 1, 2, and 3 joining
at a common junction J. The elevations of the reservoirs, the lengths of pipes, and the
diameters are all shown on the figure. The friction factor for all the pipes is f = 0.026.
a) The flow from reservoir (1) to junction J is ……………….. m3/s.
b) The flow from reservoir (2) to junction J is ……………….. m3/s.
c) The flow in pipe 3 is ……………………………….. m3/s.
d) The Elevation of water at Tank 3 is…………………………m
𝟓𝟓𝟎 𝑸𝟐𝟏
𝑷𝒊𝒑𝒆 𝟏: 𝟒𝟎 − 𝟐𝟕 = 𝟎. 𝟎𝟐𝟔 𝟐
(𝟏)
𝟎. 𝟐 𝝅 𝟐
𝟐𝒈 ( 𝟒 𝟎. 𝟐 )
𝟒𝟓𝟎 𝑸𝟐𝟐
𝑷𝒊𝒑𝒆 𝟐: 𝟑𝟎 − 𝟐𝟕 = 𝟎. 𝟎𝟐𝟔 𝟐
(𝟐)
𝟎. 𝟏𝟓 𝝅 𝟐
𝟐𝒈 ( 𝟒 𝟎. 𝟏𝟓 )
Q1 = 0.0594 m3/s
Q2 = 0.0154 m3/s
Q3 = 0.0747 m3/s
Z3 17.797 m3/s
168
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q2. given three pipes connected in parallel as shown. The total flow in the system is 7 L/s. All
the pipes are made of galvanized iron. Find the flow rates and the haedlosses in the pipes. Given
that the flow is fully turbulent.
Given:
L1=90m, D1=20mm, (A1=0.000314m2)
L2=150m, D2=40mm,(A2=0.001256m2)
L3=80m, D3=60mm, (A3=0.00282m2)
Crude oil (Sp.Gr.= 0.855)
ν= 3.810-6 m2/s
e = 0.00015 m
The previous tables shows the solution for four different values of
assumed 𝑸 ̃𝟏
It can be seen that all the solutions gave the same results.
169
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Pipe network: Interconnected pipes, could be of large numbers, through which the
water to a given outlets may come from several circuits.
Examples:
Network analysis: it is the analysis of the fluid flow through a hydraulics network,
containing several or many interconnected branches. The aim is to determine the flow rates
and pressure drops in the individual sections of the network.
170
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Considering the laws of conservation of mass and energy, the following conditions can
be made for any network of pipes:
1. The algebric sum of the energy drops around each circuit must be zero. This means
that between any two junctions in the network, the total headloss is independent of
the path taken. For example if we find the losses between the junctions a and g ,in
figure 6.2, as:
This can be given in another approach, if we consider each loop and find the losses in
each pipe , then the sum of the losses in the clockwise direction is the same as that
in the counter-clockwise direction (conservation of energy).
2. The flow into each junction must equal that out of the junction (conservation of mass
or continuity equation).
3. The Darcy-Weisbach equation must be satisfied for each pipe (proper relation between
headloss and discharge for each pipe).
171
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
The previuos coinditions can be summerized as:
The energy drop between any two points is the same and does not depend on the
path taken.
𝒉𝒇 + 𝒉𝒇 = 𝒉𝒇 + 𝒉𝒇
𝒃𝒆 𝒆𝒅 𝒃𝒄 𝒄𝒅
For each junction a,c,..,f; the total inflow must equal the total outflow.
For example:
For junction a,
𝑸𝟏 = 𝑸𝒂𝒃 + 𝑸𝒂𝒊 ,
The previous conditions should be satisfied for each loop and junction.
Based on the previous two basic principles, the pipe networks are generally solved by
the methods of successive approximation. The widely used method of pipe network analysis
is the Hardy-Cross method.
172
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
6.1.1 Finding analysis solution method:
Start by assuming a distribution of flows in each pipe such that the total inflow
equals the total outflow at each junction throughout the whole pipe network
(condition 1 is achieved).
The pipe diameters, lengths, friction cofficients, junction elevations are given.
Based on the assumed values of discharges it is expected that the headlosses
resulting in all the pipes will not achieve the second condition (conservation of
energy).
Values of headlosses in each pipe in a starting loop are computed.
The assumed pipe flows need to be adjusted until headlosses in the clock wise
direction equal headlosses in the counter clock wise direction in the loop.
The headloss in each pipe is:
𝐿 𝑄2 𝐿 1
ℎ𝑓 = 𝑓 = [𝑓 ] 𝑄2 = 𝐾𝑄2 ; (6.4)
𝐷 2𝑔𝐴2 𝐷 2𝑔𝐴2
So, we can use a constant for the losses in each pipe given as:
𝐿 1
𝐾 =𝑓 (6.5)
𝐷 2𝑔𝐴2
In any loop: the total headloss in the clockwise direction (C) is the sum of the
headlosses in all the pipes that carry flow in the clockwise direction, on the other
hand, total headloss in the counter clockwise direction (CC) is the sum of the
headlosses in all the pipes that carry flow in the counter clockwise direction (CC):
2
∑ ℎ𝑓𝐶𝐶 = ∑ 𝐾𝐶𝐶 𝑄𝐶𝐶 (6.7)
Using the assumed flowrates, Q’s , it is not expected that these two values of
equations 6.6 and 6.7 will be equal during the first trial.
The difference {∑ 𝐾𝐶 𝑄𝐶2 − ∑ 𝐾𝐶𝐶 𝑄𝐶𝐶2 }
is the closure error in the first trial.
We are looking for a flow correction ∆𝑸 that when subtracted from 𝑸𝑪 and
added to 𝑸𝑪𝑪 the values {∑ 𝐾𝐶 𝑄𝐶2 = ∑ 𝐾𝐶𝐶 𝑄𝐶𝐶 2 }
will be equalized (second
condition for network of pipes).
Thus ∆Q (the correction) must satisfy:
173
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
∑ 𝐾𝐶 (𝑄𝐶 − ∆𝑄)2 = ∑ 𝐾𝐶𝐶 (𝑄𝐶𝐶 + ∆𝑄)2 (6.8)
ℎ𝑓
Given that: ℎ𝑓 = 𝐾𝑄2 𝐾𝑄 =
𝑄
∑ ℎ𝑓𝐶 −∑ ℎ𝑓𝐶𝐶
∆𝑄 = ℎ𝑓 ℎ𝑓 (6.11)
𝐶
2(∑ +∑ 𝐶𝐶 )
𝑄𝐶 𝑄𝐶𝐶
Example 6.1:
Find the missing flow values in pipes 2, 7, 8; and the input and output values to the
network shown.
174
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 6.2:
A water-supply distribution system for a park is schematically shown. The demands on
different joints of the system are shown on figure 6.4 with flow rates given in litres per
second. Water enters the system at junction A from a water storage tank. All pipes are aged
ductile iron (e=0.26mm) .The water surface elevation in the tank is 50 m above the elevation
of point A. The elevations of the junctions and the lengths of the pipes are provided in the
tables below. Calculate the flow rate in each pipe, the pressure at each junction and the
power required in each pipe. Check if the pressure values will be high enough to satisfy the
customer requirement at each junction (185 kPa). If any of the junctions did not satisfy the
required pressure, suggest modifications to the system to accommodate the customer.
𝐿 1
The values of K were found from the relation: 𝐾 = 𝑓 .
𝐷 2𝑔𝐴2
175
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
A suggested distribution of flow was made as the following:
The direction of flow in the pipes is indicated in the table as the sequence of the letters. For
example the direction of flow in pipe AB is from junction A to B. This direction may change
within the solution ( be reversed), the letters will be reversed then.
The clockwise specifier (C) means only the flows that are moving clockwise in our loop, while
the counter-clockwise specifier (CC) is only the flows that are moving counter-clockwise. It
can be indicated in the table by putting the values for the counterclockwise flow between
brackets. [...].
176
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Using Hardy-Cross method to solve the question requires:
Taking the loops one at a time.
𝑳 𝟏
The losses withen each pipe is found using the Darcy-Weisbach equation 𝑲 = 𝒇 .
𝑫 𝟐𝒈𝑨𝟐
The correction to the flow ∆Q is calculated then.
This value (∆Q) is subtracted from the flow discharges in the clockwise direction and
added to the flows in the counter clockwise direction. This will not change the first
condition (continuity equation) and will make the distribution closer to the correct
distribution that achieve the second condition (energy equation balance) 𝑸 ± ∆𝑸.
Before starting the 2nd loop, the flow in the common pipes should be updated from the
previous corrections.
1st cycle 2nd cycle
Loop
Q hf/Q hf Q QNew
(1) K QNew (m3/s) hf/Q hf
(m3/s) (KQ) (KQ2) (m3/s) (m3/s)
Pipe
AB 0.200 194 38.8 7.76 0.2-(-0.005)=0.205
BG 0.120 678 81.3 9.76 0.12-(-0.005)=0.125
GE 0.000 2,990 0.0 0.0 0-(-0.005)=0.005
AD [0.100] 423 [42.3] [4.23] 0.1+(-0.005)=0.0095
DE [0.100] 1,630 [163.0] [16.3] 0.1+(-0.005)=0.0095
∑ ℎ𝑓𝐶 −∑ ℎ𝑓𝐶𝐶 (7.76+9.76)−(4.23+16.3)
First cycle: ∆𝑄 = ℎ𝑓 ℎ𝑓 = = −0.005 𝑚3 /𝑠
2[(38.8+81.3)+(42.3+163.0)]
2(∑ 𝐶 +∑ 𝐶𝐶 )
𝑄𝐶 𝑄𝐶𝐶
177
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Note that the sign of New Q for pipe GE was negative after the correction in the first cycle,
which requires reversing the direction of flow in the pipe.
1st cycle 2nd cycle
Loop
Q hf/Q hf Q QNew
(1) K QNew (m3/s) hf/Q hf
(m3/s) (KQ) (KQ2) (m3/s) (m3/s)
Pipe
AB 0.200 194 38.8 7.76 0.2-(-0.005)=0.205 0.205 39.8 8.15 0.205
BG 0.120 678 81.3 9.76 0.12-(-0.005)=0.125 0.127 85.8 10.9 0.127
GE 0.000 2,990 0.0 0.0 0-(-0.005)=0.005 [0.008] [23.8] [0.19] [0.008]
AD [0.100] 423 [42.3] [4.23] 0.1+(-0.005)=0.0095 [0.095] [40.2] [3.82] [0.095]
DE [0.100] 1,630 [163.0] [16.3] 0.1+(-0.005)=0.0095 [0.095] [154.7] [14.7] [0.095]
∑ ℎ𝑓𝐶 −∑ ℎ𝑓𝐶𝐶 (8.15+10.9)−(3.82+14.7+0.19)
Second cycle: ∆𝑄 = ℎ𝑓 ℎ𝑓 = = +0.000 𝑚3 /𝑠
2[(39.8+85.8)+(40.2+154.7+23.8)]
2(∑ 𝐶 +∑ 𝐶𝐶 )
𝑄𝐶 𝑄𝐶𝐶
Because the corrections are very small on all three loops, the flow rates can be accepted and
the process can be ended.
178
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Finding the final pressure at each junction:
Q (L/s) hf (m) Power(Kw)
Pipe hf (kPa)
(final) 2
(kQ ) ( Q hf)
AB 205 8.2 80.3 16.4615
AD 95 3.8 37.2 3.534
BC 80 12.2 119.4 9.552
BG 125 10.6 103.8 12.975
GH 33 2.1 20.6 0.6798
CH 30 0.6 5.9 0.177
DE 95 14.7 143.9 13.6705
GE 8 0.2 2.0 0.016
EF 87 14.4 141.0 12.267
HF 63 11.9 116.5 7.3395
Energy equation can be used to find the pressure at each junction be considering two
successive junction as shown:
179
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Joint A B C D E F G H
Elevation(m) 0 0.8 1 0.9 0.8 1 1.1 1.6
Pressure (kPa) 489.5 401.368 280.01 443.489 300.568 157.61 294.631 268.236
The pressure at F is less than 157.61 kPa. Modifications should be made to the pipe
diameters to solve this issue.
180
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 6.3:
Determine the distribution of flow in the pipe network shown using Hardy cross
method.
AB 30 2 60 1800 81.96 81.96 163.92 13434.88 79.786 79.786 159.572 12731.6 77.577
BE 10 3 30 300 61.96 53.074 159.222 8450.548 50.9 56.4676 169.403 9565.77 53.797
DE ]150[ 2 300 45000 98.04 ]98.04[ 196.08 19223.68 100.214 100.214 200.428 20085.7 102.42
AD ]70[ 1 70 4900 18.04 ]18.04[ 18.04 325.44 20.214 20.214 20.214 408.606 22.42
BC 50 1 50 2500 57.88 57.88 57.88 3350.09 57.746 57.746 57.746 3334.6 58.078
CF 60 2 120 7200 67.88 67.88 135.76 9215.39 67.7462 67.7462 135.492 9179.1 68.0782
EF ]70[ 1 70 4900 62.11 62.11 62.11 3857.65 62.2538 62.2538 62.2538 3875.54 61.9318
BE ]61.96[ 3 185.88 11517.125 53.074 50.9 152.7 7772.43 56.4676 53.797 161.391 8682.35 53.475
Tutorial 9
Q1. Perform pipe network analysis and calculate water flow in all branches.
Q2. Perform pipe network analysis and calculate water flow in all branches. Hazen Williams
coefficient for each pipe is provided in table below.
182
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Extra Reading:
An initial flow estimate is made across all pipe branches that satisfies continuity at all
nodes.
Head loss across each pipe is determined using Darcy – Weisbach equation:
hf = K Q²
𝐿 1 8𝑓𝐿
K =𝑓 𝜋 2 = 2 5
𝐷 2𝑔( 𝐷2 ) 𝜋 𝑔𝐷
4
The equation for head loss across each loop is made by summing the loss in each pipe
of the loop assuming the sum of the losses to be zero:
F(I) : K1.Q1² + K3.Q3² - K8.Q8² - K4.Q4² - K2.Q2² = 0
F(II) : K5.Q5² + K7.Q7² - K6.Q6² - K3.Q3² = 0
F(III): K6.Q6² + K10.Q10² - K9.Q9² + K8.Q8² = 0
Newton Raphson method is used to solve above equations for flow correction factor
ΔQ (ΔQI, ΔQII, ΔQIII) across each loop.
In vector form all loops can be written as :
183
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
JL.ΔQ = -F(Qm-1)
Where,
- ΔQ is a vector of loop flow corrections
- JL is the first derivatives of the loop equations evaluated at Qm-1
Coefficient Matrix, JL
𝜕(𝐾1𝑄12 )
= 2𝐾1. 𝑄1
𝜕(△ 𝑄1 )
Diagonal terms are obtained by adding the derivatives of all the pipes in the loop and
always have a positive sign:
𝜕(𝐹𝐼 )
= 2(|𝐾1𝑄1 | + |𝐾3𝑄3 | + |𝐾8𝑄8 | + ⌈𝐾4⌉𝑄4 + |𝐾2𝑄2 |) = 979.03
𝜕(△ 𝑄𝐼 )
𝜕(𝐹2 )
= 2(|𝐾5𝑄5 | + |𝐾7𝑄7 | + |𝐾6𝑄6 | + ⌈𝐾3⌉𝑄3 ) = 794.35
𝜕(△ 𝑄2 )
𝜕(𝐹3 )
= 2(|𝐾6𝑄6 | + |𝐾10𝑄10 | + |𝐾9𝑄9 | + ⌈𝐾8⌉𝑄8 ) = 1408.23
𝜕(△ 𝑄3 )
Off-diagonal terms are gradients for pipes that appear in different loops (common
pipes between two loops) and always have a negative sign.
𝜕(𝐹𝐼 ) 𝜕(𝐹 )
= −𝑛(|𝐾3𝑄3 |) = 𝜕(△𝑄𝐼𝐼 ) = −165.6 (common pipe between loop 1 and 2 (pipe 3))
𝜕(△𝑄𝐼𝐼 ) 𝐼
𝜕(𝐹𝐼 ) 𝜕(𝐹𝐼𝐼𝐼 )
= −𝑛(|𝐾8𝑄8 |) = = −314.81 (common pipe between loop 1 and 3 (pipe 8))
𝜕(△𝑄𝐼𝐼𝐼 ) 𝜕(△𝑄𝐼 )
𝜕(𝐹𝐼𝐼 ) 𝜕(𝐹 )
= −𝑛(|𝐾6𝑄6 |) = 𝜕(△𝑄𝐼𝐼𝐼 ) = −273.72 (common pipe between loop 2 and 3 (pipe 6))
𝜕(△𝑄𝐼𝐼𝐼 ) 𝐼𝐼
In the above derivative change in Loop I due to flow change in Loop II will be due to
common pipe 3. Above gradient is also for δFII/δ(ΔQI).
184
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
The JL matrix is obtained as the following:
- F(Qm-1) is the vector of the sum of energy losses with negative sign given for the flow
counter clockwise
[𝑲𝟏. 𝑸𝟐𝟏 + 𝑲𝟑. 𝑸𝟐𝟑 − 𝑲𝟖. 𝑸𝟐𝟖 − 𝑲𝟒. 𝑸𝟐𝟒 − 𝑲𝟐. 𝑸𝟐𝟐 𝑲𝟓. 𝑸𝟐𝟓 + 𝑲𝟕. 𝑸𝟐𝟕 − 𝑲𝟔. 𝑸𝟐𝟔 − 𝑲𝟑. 𝑸𝟐𝟑 𝑲𝟔. 𝑸𝟐𝟔 + 𝑲𝟏𝟎. 𝑸𝟐𝟏𝟎 − 𝑲𝟗. 𝑸𝟐𝟗 + 𝑲𝟖. 𝑸𝟐𝟖 ]
𝐾1. 𝑄21 + 𝐾3. 𝑄23 − 𝐾8. 𝑄28 − 𝐾4. 𝑄24 − 𝐾2. 𝑄22 = 10.94
𝐾5. 𝑄25 + 𝐾7. 𝑄27 − 𝐾6. 𝑄26 − 𝐾3. 𝑄23 = 6.4
𝐾6. 𝑄26 + 𝐾10. 𝑄210 − 𝐾9. 𝑄29 + 𝐾8. 𝑄28 = −43.62
F matrix is obtained as:
F(Qm-1)T = [-10.94 -6.40 43.62 ]
Below, the matrix is solved in excel to obtain ΔQ.
JL.ΔQ = - F
JL-1 is determined using Minverse function. Select 3×3 cells in excel and type
MINVERSE(Input Array) and press CTRL+SHIFT+ENTER to evaluate inverse of matrix.
Do multiplication of JL-1 and -F vector using MMULT function. Select 3×1 cells in excel
and type MMULT(Inverse Array, F array) and press CTRL+SHIFT+ENTER to evaluate
multiplication. ΔQ array is obtained as following: ΔQ = [ 0.00066 -0.00261 -0.03133]
As flow changes are larg, another iteration is done, with flows adjusted based on ΔQ
values. For example flow rate in Pipe 3 will become.
Q3 = Q30 + ΔQI - ΔQII
= 0.12 + (0.00066) - (-0.00261)
= 0.123
New flowrates are calculated and above steps are repeated.
185
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
To do more iterations copy entire rows and paste them below the above cells to carry
out further calculation till change in flowrates become negligible. For this example ΔQ
becomes negligible in 5 iterations. Final flows in m³/s are as following.
If a pipe’s flow direction changes from the assumed value, the signs for that pipe head loss
terms are switched for all loops containing the pipe during the next iteration in loop
equations. JL matrix signs will remain same as above.
JL.ΔQ = - F
JL-1 is determined using Minverse function. Select 3×3 cells in excel and type
MINVERSE(Input Array) and press CTRL+SHIFT+ENTER to evaluate inverse of matrix.
Do multiplication of JL-1 and -F vector using MMULT function. Select 3×1 cells in excel
and type MMULT(Inverse Array, F array) and press CTRL+SHIFT+ENTER to evaluate
multiplication. ΔQ array is obtained as following.
186
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
As flow changes are larger another iteration is done, with flows adjusted based on ΔQ
values. For example flow rate in Pipe 3 will become.
= 0.123
New flowrates are calculated and above steps are repeated. To do more iterations copy
entire rows and paste them below the above cells to carry out further calculation till
change in flowrates become negligible. For this example ΔQ becomes negligible in 5
iterations. Final flows in m³/s are as following.
If a pipe’s flow direction changes from the assumed value, the signs for that pipe head loss
terms are switched for all loops containing the pipe during the next iteration in loop
equations. JL matrix signs will remain same as above.
187
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example
Perform pipe network analysis and calculate water flow in all branches.
Iteration 1
Head
Pipe Flow Length Diameter e/D Velocity Reynold's Friction Loss
Factor,
m3/s m m m/s Number f K hL nHL/Q
1 0.100 1000 0.40 0.00087 0.80 318310 0.020 162 1.62 32.49
2 0.200 2000 0.45 0.00104 1.26 565884 0.020 183 7.33 73.28
3 0.030 1000 0.30 0.00130 0.42 127324 0.023 779 0.70 46.74
4 0.030 2000 0.30 0.00130 0.42 127324 0.023 1558 1.40 93.48
5 0.070 1000 0.40 0.00130 0.56 222817 0.022 179 0.88 25.06
6 0.080 750 0.30 0.00130 1.13 339531 0.022 556 3.56 89.01
7 0.020 2200 0.30 0.00130 0.28 84883 0.024 1770 0.71 70.80
8 0.170 2000 0.30 0.00173 2.41 721502 0.023 1559 45.05 529.98
9 0.020 500 0.25 0.00130 0.41 101859 0.023 986 0.39 39.42
10 0.020 2200 0.25 0.00173 0.41 101859 0.025 4574 1.83 182.95
Coefficient
Matrix Inverse F DQ
245.99 -46.74 -93.48 4.72E-03 2.76E-04 1.59E-03 -5.00 -0.04182
-46.74 799.10 0.00 2.76E-04 1.27E-03 9.27E-05 -46.57 -0.06073
-93.48 0.00 278.36 1.59E-03 9.27E-05 4.12E-03 -3.38 -0.02618
188
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Iteration2
Head
Pipe Flow Length Diameter e/D Velocity Reynold's Friction Loss
Factor,
m3/s m m m/s Number f K hL nHL/Q
1 0.142 1000 0.40 0.00087 1.13 451438 0.020 160 3.22 45.40
2 0.158 2000 0.45 0.00104 0.99 447548 0.021 184 4.62 58.35
3 0.049 1000 0.30 0.00130 0.69 207547 0.022 757 1.81 74.03
4 0.046 2000 0.30 0.00130 0.65 193718 0.022 1519 3.16 138.68
5 0.096 1000 0.40 0.00130 0.77 306150 0.022 177 1.63 33.98
6 0.054 750 0.30 0.00130 0.76 228420 0.022 565 1.64 60.83
7 0.046 2200 0.30 0.00130 0.65 195993 0.022 1670 3.56 154.24
8 0.109 2000 0.30 0.00173 1.55 463775 0.023 1571 18.76 343.33
9 -0.041 500 0.25 0.00130 0.83 207414 0.022 942 -1.56 76.71
10 -0.041 2200 0.25 0.00173 0.83 207414 0.024 4404 -7.30 358.74
Coefficient
Matrix Inverse F DQ
316.47 -74.03 -138.68 3.84E-03 3.33E-04 1.37E-03 -0.04 -0.00231
-74.03 852.81 0.00 3.33E-04 1.20E-03 1.19E-04 -8.08 -0.00968
-138.68 0.00 387.73 1.37E-03 1.19E-04 3.07E-03 0.39 0.00019
Iteration3
Head
Pipe Flow Length Diameter e/D Velocity Reynold's Friction Loss
Factor,
m3/s m m m/s Number f K hL nHL/Q
1 0.144 1000 0.40 0.00087 1.15 458788 0.020 160 3.32 46.12
2 0.156 2000 0.45 0.00104 0.98 441015 0.021 185 4.48 57.53
3 0.056 1000 0.30 0.00130 0.80 238817 0.022 752 2.38 84.63
4 0.048 2000 0.30 0.00130 0.68 204325 0.022 1515 3.51 145.87
5 0.096 1000 0.40 0.00130 0.76 305545 0.022 177 1.63 33.92
6 0.054 750 0.30 0.00130 0.76 229227 0.022 565 1.65 61.03
7 0.046 2200 0.30 0.00130 0.65 195187 0.022 1670 3.53 153.64
8 0.100 2000 0.30 0.00173 1.41 422705 0.023 1574 15.62 313.57
9 -0.050 500 0.25 0.00130 1.03 256698 0.022 933 -2.37 94.02
10 -0.050 2200 0.25 0.00173 1.03 256698 0.023 4370 -11.10 440.48
Coefficient
Matrix Inverse F DQ
334.15 -84.63 -145.87 3.67E-03 3.33E-04 1.36E-03 -0.03 -0.00003
-84.63 932.69 0.00 3.33E-04 1.10E-03 1.23E-04 0.24 0.00025
-145.87 0.00 394.47 1.36E-03 1.23E-04 3.04E-03 0.00 -0.00001
189
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Iteration4
Head
Pipe Flow Length Diameter e/D Velocity Reynold's Friction Loss
Factor,
m3/s m m m/s Number f K hL nHL/Q
1 0.144 1000 0.40 0.00087 1.15 458879 0.020 160 3.32 46.13
2 0.156 2000 0.45 0.00104 0.98 440934 0.021 185 4.48 57.52
3 0.056 1000 0.30 0.00130 0.79 237634 0.022 752 2.36 84.23
4 0.048 2000 0.30 0.00130 0.68 204414 0.022 1515 3.51 145.93
5 0.096 1000 0.40 0.00130 0.76 305569 0.022 177 1.63 33.92
6 0.054 750 0.30 0.00130 0.76 229195 0.022 565 1.65 61.03
7 0.046 2200 0.30 0.00130 0.65 195218 0.022 1670 3.53 153.66
8 0.100 2000 0.30 0.00173 1.41 423766 0.023 1574 15.69 314.34
9 -0.050 500 0.25 0.00130 1.02 255424 0.022 933 -2.35 93.57
10 -0.050 2200 0.25 0.00173 1.02 255424 0.023 4370 -10.99 438.37
Coefficient
Matrix Inverse F DQ
333.81 -84.23 -145.93 3.67E-03 3.33E-04 1.36E-03 0.00 0.00000
-84.23 930.50 0.00 3.33E-04 1.10E-03 1.23E-04 0.00 0.00000
-145.93 0.00 394.55 1.36E-03 1.23E-04 3.04E-03 0.00 0.00000
190
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example
Perform pipe network analysis and calculate water flow in all branches. Hazen Williams
coefficient for each pipe is provided in table below.
191
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Iteration 1
Head
Pipe Flow Diameter Length C Loss
m3/s mm m HL, m nHL/Q
1 1.000 1000 1000 100 2.13 3.94
2 6.000 750 925 100 219.96 67.82
3 3.000 750 1000 100 65.96 40.68
4 1.000 750 925 100 7.99 14.79
5 3.000 500 350 100 166.31 102.56
6 2.000 500 671 100 150.59 139.30
7 1.000 500 400 100 24.90 46.07
8 1.000 500 650 100 40.47 74.86
9 1.000 750 1000 120 6.17 11.41
10 6.000 1000 1000 120 41.81 12.89
11 6.000 1000 925 120 38.67 11.92
12 6.000 1000 800 120 33.45 10.31
13 1.000 750 763 120 4.71 8.71
14 4.000 1000 650 120 12.84 5.94
15 0.500 750 400 120 0.68 2.53
16 0.500 750 125 120 0.21 0.79
17 0.500 750 400 120 0.68 2.53
18 0.500 500 125 120 1.54 5.70
19 1.000 500 800 120 35.55 65.76
20 2.000 500 125 120 20.02 18.52
21 2.000 500 800 120 128.14 118.53
22 1.000 500 1000 120 44.43 82.20
Coefficient Matrix F
127.22 -67.8 -40.7 0 0 0 0 -280.05
-67.8 104.05 0 -11.41 0 0 0 229.26
-40.7 0 403.47 -74.9 0 0 0 367.30
0 -11.41 -74.86 303.25 -5.94 -118.5 0 -95.69
0 0 0 -5.94 17.17 0 -2.53 8.81
0 0 0 -118.5 0 208.51 -5.70 182.17
0 0 0 0 -2.53 -5.70 11.56 1.33
Inverted Matrix
1.3E-02 8.5E-03 1.5E-03 8.9E-04 3.6E-04 5.1E-04 3.3E-04 ΔQ1 -1.10677
8.5E-03 1.5E-02 1.1E-03 1.1E-03 4.4E-04 6.3E-04 4.1E-04 ΔQ2 1.52773
1.5E-03 1.1E-03 2.8E-03 9.6E-04 3.8E-04 5.5E-04 3.6E-04 ΔQ3 0.87608
8.9E-04 1.1E-03 9.6E-04 4.7E-03 1.9E-03 2.7E-03 1.7E-03 ΔQ4 0.41659
3.6E-04 4.4E-04 3.8E-04 1.9E-03 6.1E-02 1.4E-03 1.4E-02 ΔQ5 0.78179
5.1E-04 6.3E-04 5.5E-04 2.7E-03 1.4E-03 6.4E-03 3.5E-03 ΔQ6 1.13359
3.3E-04 4.1E-04 3.6E-04 1.7E-03 1.4E-02 3.5E-03 9.1E-02 ΔQ7 0.84530
192
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Iteration2
Head
Pipe Flow Diameter Length C Loss
m3/s mm m HL, m nHL/Q
1 -0.107 1000 1000 100 -0.03 0.59
2 3.366 750 925 100 75.47 41.49
3 1.017 750 1000 100 8.92 16.22
4 2.107 750 925 100 31.73 27.86
5 2.124 500 350 100 87.79 76.47
6 1.124 500 671 100 51.85 85.35
7 0.124 500 400 100 0.52 7.81
8 1.459 500 650 100 81.45 103.24
9 -0.111 750 1000 120 -0.11 1.76
10 4.472 1000 1000 120 24.27 10.04
11 7.528 1000 925 120 58.84 14.46
12 6.417 1000 800 120 37.87 10.92
13 1.782 750 763 120 13.70 14.23
14 3.635 1000 650 120 10.75 5.47
15 0.564 750 400 120 0.85 2.80
16 1.345 750 125 120 1.33 1.84
17 -0.345 750 400 120 -0.35 1.85
18 0.788 500 125 120 3.58 8.39
19 -0.134 500 800 120 -0.86 11.88
20 0.866 500 125 120 4.26 9.10
21 1.283 500 800 120 56.36 81.27
22 0.583 500 1000 120 16.40 51.99
Coefficient Matrix F
86.16 -41.5 -16.2 0 0 0 0 -52.63
-41.5 67.75 0 -1.76 0 0 0 40.81
-16.2 0 289.08 -103.2 0 0 0 67.64
0 -1.76 -103.24 254.66 -5.47 -81.3 0 -7.04
0 0 0 -5.47 22.50 0 -2.80 -2.10
0 0 0 -81.3 0 110.64 -8.39 56.19
0 0 0 0 -2.80 -8.39 14.88 1.04
Inverted Matrix
1.7E-02 1.0E-02 1.2E-03 7.5E-04 2.3E-04 5.8E-04 3.7E-04 ΔQ1 -0.35431
1.0E-02 2.1E-02 8.0E-04 6.3E-04 1.9E-04 4.9E-04 3.1E-04 ΔQ2 0.39534
1.2E-03 8.0E-04 4.4E-03 2.4E-03 7.3E-04 1.8E-03 1.2E-03 ΔQ3 0.35149
7.5E-04 6.3E-04 2.4E-03 6.6E-03 2.0E-03 5.1E-03 3.2E-03 ΔQ4 0.38474
2.3E-04 1.9E-04 7.3E-04 2.0E-03 4.6E-02 2.2E-03 9.9E-03 ΔQ5 0.06931
5.8E-04 4.9E-04 1.8E-03 5.1E-03 2.2E-03 1.3E-02 8.0E-03 ΔQ6 0.83237
3.7E-04 3.1E-04 1.2E-03 3.2E-03 9.9E-03 8.0E-03 7.4E-02 ΔQ7 0.55268
193
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Iteration3
Head
Pipe Flow Diameter Length C Loss
m3/s mm m HL, m nHL/Q
1 -0.461 1000 1000 100 -0.51 2.04
2 2.616 750 925 100 47.35 33.49
3 0.311 750 1000 100 1.00 5.93
4 2.461 750 925 100 42.30 31.80
5 1.772 500 350 100 62.82 65.57
6 0.772 500 671 100 25.91 62.05
7 -0.228 500 400 100 -1.61 13.09
8 1.426 500 650 100 78.05 101.24
9 -0.122 750 1000 120 -0.13 1.91
10 4.077 1000 1000 120 20.45 9.28
11 7.923 1000 925 120 64.68 15.10
12 6.801 1000 800 120 42.18 11.47
13 1.851 750 763 120 14.70 14.69
14 3.950 1000 650 120 12.54 5.87
15 1.047 750 400 120 2.69 4.75
16 1.898 750 125 120 2.52 2.46
17 -0.898 750 400 120 -2.02 4.17
18 1.068 500 125 120 6.27 10.87
19 -0.966 500 800 120 -33.34 63.85
20 0.034 500 125 120 0.01 0.58
21 0.835 500 800 120 25.48 56.44
22 0.199 500 1000 120 2.23 20.81
Coefficient Matrix F
73.26 -33.5 -5.9 0 0 0 0 -5.54
-33.5 59.78 0 -1.91 0 0 0 3.00
-5.9 0 247.88 -101.2 0 0 0 10.07
0 -1.91 -101.24 197.73 -5.87 -56.4 0 0.21
0 0 0 -5.87 25.31 0 -4.75 0.52
0 0 0 -56.4 0 131.73 -10.87 -14.12
0 0 0 0 -4.75 -10.87 22.24 -0.96
Inverted Matrix
1.8E-02 1.0E-02 6.4E-04 5.0E-04 1.4E-04 2.2E-04 1.4E-04 ΔQ1 -0.06775
1.0E-02 2.3E-02 4.6E-04 5.3E-04 1.5E-04 2.4E-04 1.5E-04 ΔQ2 0.01168
6.4E-04 4.6E-04 5.3E-03 3.2E-03 9.1E-04 1.4E-03 8.9E-04 ΔQ3 0.03167
5.0E-04 5.3E-04 3.2E-03 7.8E-03 2.2E-03 3.5E-03 2.2E-03 ΔQ4 -0.01795
1.4E-04 1.5E-04 9.1E-04 2.2E-03 4.2E-02 1.8E-03 9.8E-03 ΔQ5 -0.00297
2.2E-04 2.4E-04 1.4E-03 3.5E-03 1.8E-03 9.5E-03 5.0E-03 ΔQ6 -0.12344
1.4E-04 1.5E-04 8.9E-04 2.2E-03 9.8E-03 5.0E-03 5.0E-02 ΔQ7 -0.10401
194
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Iteration4
Head
Pipe Flow Diameter Length C Loss
m3/s mm m HL, m nHL/Q
1 -0.529 1000 1000 100 -0.66 2.29
2 2.536 750 925 100 44.73 32.62
3 0.212 750 1000 100 0.49 4.28
4 2.529 750 925 100 44.48 32.54
5 1.741 500 350 100 60.76 64.57
6 0.741 500 671 100 23.98 59.88
7 -0.259 500 400 100 -2.05 14.62
8 1.476 500 650 100 83.14 104.22
9 -0.151 750 1000 120 -0.19 2.29
10 4.065 1000 1000 120 20.35 9.26
11 7.935 1000 925 120 64.86 15.12
12 6.783 1000 800 120 41.97 11.45
13 1.848 750 763 120 14.66 14.67
14 3.935 1000 650 120 12.45 5.85
15 0.946 750 400 120 2.23 4.35
16 1.794 750 125 120 2.27 2.34
17 -0.794 750 400 120 -1.61 3.75
18 1.049 500 125 120 6.06 10.70
19 -0.843 500 800 120 -25.89 56.85
20 0.157 500 125 120 0.18 2.14
21 0.941 500 800 120 31.75 62.44
22 0.217 500 1000 120 2.62 22.40
Coefficient Matrix F
71.73 -32.6 -4.3 0 0 0 0 -0.08
-32.6 59.30 0 -2.29 0 0 0 0.03
-4.3 0 247.58 -104.2 0 0 0 0.03
0 -2.29 -104.22 208.65 -5.85 -62.4 0 -0.22
0 0 0 -5.85 24.88 0 -4.35 0.02
0 0 0 -62.4 0 132.12 -10.70 -0.02
0 0 0 0 -4.35 -10.70 21.15 -0.05
Inverted Matrix
1.9E-02 1.0E-02 5.0E-04 4.3E-04 1.3E-04 2.2E-04 1.3E-04 ΔQ1 -0.00132
1.0E-02 2.3E-02 4.0E-04 5.3E-04 1.5E-04 2.6E-04 1.7E-04 ΔQ2 -0.00031
5.0E-04 4.0E-04 5.4E-03 3.2E-03 9.3E-04 1.6E-03 1.0E-03 ΔQ3 -0.00062
4.3E-04 5.3E-04 3.2E-03 7.6E-03 2.2E-03 3.8E-03 2.4E-03 ΔQ4 -0.00174
1.3E-04 1.5E-04 9.3E-04 2.2E-03 4.2E-02 1.8E-03 9.6E-03 ΔQ5 -0.00009
2.2E-04 2.6E-04 1.6E-03 3.8E-03 1.8E-03 9.8E-03 5.3E-03 ΔQ6 -0.00116
1.3E-04 1.7E-04 1.0E-03 2.4E-03 9.6E-03 5.3E-03 5.2E-02 ΔQ7 -0.00276
195
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Iteration5
Head
Pipe Flow Diameter Length C Loss
m3/s mm m HL, m nHL/Q
1 -0.530 1000 1000 100 -0.66 2.30
2 2.535 750 925 100 44.69 32.61
3 0.211 750 1000 100 0.49 4.26
4 2.530 750 925 100 44.52 32.55
5 1.741 500 350 100 60.80 64.59
6 0.741 500 671 100 24.01 59.92
7 -0.259 500 400 100 -2.04 14.59
8 1.477 500 650 100 83.26 104.29
9 -0.153 750 1000 120 -0.19 2.31
10 4.066 1000 1000 120 20.35 9.26
11 7.934 1000 925 120 64.85 15.12
12 6.782 1000 800 120 41.95 11.44
13 1.848 750 763 120 14.66 14.67
14 3.934 1000 650 120 12.44 5.85
15 0.943 750 400 120 2.21 4.34
16 1.791 750 125 120 2.27 2.34
17 -0.791 750 400 120 -1.60 3.74
18 1.050 500 125 120 6.08 10.71
19 -0.841 500 800 120 -25.82 56.78
20 0.159 500 125 120 0.18 2.15
21 0.940 500 800 120 31.72 62.41
22 0.218 500 1000 120 2.66 22.55
Coefficient Matrix F
71.73 -32.6 -4.3 0 0 0 0 0.00
-32.6 59.30 0 -2.31 0 0 0 0.00
-4.3 0 247.66 -104.3 0 0 0 0.00
0 -2.31 -104.29 208.85 -5.85 -62.4 0 0.00
0 0 0 -5.85 24.87 0 -4.34 0.00
0 0 0 -62.4 0 132.05 -10.71 0.00
0 0 0 0 -4.34 -10.71 21.14 0.00
Inverted Matrix
1.9E-02 1.0E-02 5.0E-04 4.3E-04 1.3E-04 2.2E-04 1.3E-04 ΔQ1 0.00000
1.0E-02 2.3E-02 4.0E-04 5.3E-04 1.5E-04 2.7E-04 1.7E-04 ΔQ2 0.00000
5.0E-04 4.0E-04 5.4E-03 3.2E-03 9.3E-04 1.6E-03 1.0E-03 ΔQ3 0.00000
4.3E-04 5.3E-04 3.2E-03 7.6E-03 2.2E-03 3.8E-03 2.4E-03 ΔQ4 0.00000
1.3E-04 1.5E-04 9.3E-04 2.2E-03 4.2E-02 1.8E-03 9.6E-03 ΔQ5 0.00000
2.2E-04 2.7E-04 1.6E-03 3.8E-03 1.8E-03 9.8E-03 5.3E-03 ΔQ6 0.00000
1.3E-04 1.7E-04 1.0E-03 2.4E-03 9.6E-03 5.3E-03 5.2E-02 ΔQ7 0.00000
196
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
197
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
198
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Show the suggested flow distribution within the network:
199
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Show that the” inflow=outflow” for all the nodes within the network:
200
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Do at least three cycles of corrections:
201
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
202
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Do at least three cycles of corrections:
203
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Open channel flow assumes that the pressure at the surface is constant (atmospheric)
and the hydraulic grade line matches the surface of the fluid.
The other type of flow within a conduit is a pipe pressurized flow. These two types of
flow are similar in many ways, but differ in one important aspect: open-channel flow has a
free surface, whereas pipe pressurized flow does not.
An open channel is a physical system in which water flows with a free surface. A
channel can be classified as either natural or artificial channel according to its origin:
204
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Artificial channels are those constructed or developed by human effort such as:
Gutters
1. A channel at the edge of a street or road for carrying off surface water.
2. A trough fixed under or along the eaves for draining rainwater from a roof. Also
called regionally eaves trough, rainspout, spouting.
3. A furrow or groove formed by running water.
4. A trough or channel for carrying something off, such as that on either side of a
bowlin galley or that almost level with the water in some swimming pools.
Drainage channel: A channel for conveying storm-water runoff; usually lined with
concrete, grass, riprap, or the like, to reduce erosion of the channel.
Ditches: a long, narrow excavation made in the ground by digging, as for draining or
irrigating land.
Trench.
Floodways.
Tunnels.
log chutes:(chute: an inclined channel, as a trough, tube, or shaft,for conveying water
to a lowerlevel, or a waterfall or steep descent, as in a river)
Navigation channels: A deeper channel cut into the sea or river bed, to enable
larger ships to pass through to a port.
Power canals: a canal used for hydraulic power generation, rather than for transport
of watercraft.
Trough:a channel or conduit for conveying water, as agutter under the eaves of a buil
ding for carryingaway rain water.
spillways
Model channels that are built in the laboratory for experimental investigation studies.
Long distance canals: constructed to achieve the interbasin transfer of water at
National and International levels.
The artificial channel is known by different names, such as " canal "," chute", "culvert",
"drop", "flumes" and "open - flow tunnel", Aqueduct.
The canal is usually a long and mild-sloped channel built in the ground, which may be
lined or unlined with stone masonry, concrete, cement, wood or bituminous materials etc.
205
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
The main differences bewteen open channel flow and pipe pressurized flow are:
- Pipe pressurized flow fills the entire conduit completely, so that the flow boundaries
are fixed by the conduit geometry: Aflow = Aboundary . While the area of the flow in non-
pressurized open channel flow is related to the flow parameters.
- Pressurized pipe flow posseses a hydraulic pressure that varies from one section to
another along the pipe line. While, open channel flow has a free surface under
atmospheric pressure where the depth adjust itself depending on the flow conditions
(Q, bed slope So, boundary friction, geometry). We need to emphasize that the free
surface is subject to atmospheric pressure which remains constant throughout the
entire length of the channel.
206
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
207
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
- The free water surface in figure (b) is Figure 7.4 Comparison between the energy for the
subjected to atmospheric pressure only (zero pressurized pipe flow and the non-pressurized
pressure reference). flow in open channel.
208
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
The following table summerizes the differences between pressurized pipe flow and the free
surface (open channel) flow:
The upper surface of the fluid is exposed to A pipe is a closed conduit and the fluid is
atmosphere. under pressure.
The driving force of flow is the gravity The flow in a pipe is termed as pipe flow only
component at the direction of flow. when the fluid completely fills the cross
Flow conditions are greatly influenced by section with no free surface of fluid.
the channel slope. Hydraulic grade line does not coincide with
Hydraulic grade line coincides with the the water surface.
water surface. The maximum velocity occurs at the pipe
The maximum velocity occurs at a little centre.
distance below the water surface. Velocity Distribution is symmetrical about
The shape of the velocity profile is the pipe axis
dependent on the channel roughness.
Open channel flow: the flow of water in any conduit with the fluid surface exposed to
the atmosphere (free surface).
Open channel hydraulics: the study of the mechanics of water flowing in open
conduits, which generally include channels, streams, and pipes (not flowing under
pressure).
Uniform presmatic channel: the channel maintains a constant shape and slope.
Discharge: volume flow of the fluid passing through a flow section per unit time.
The slope of the channel bottom (So): the vertical fall of the channel bottom divided
by the horizontal run of the bottom.
209
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 7.1
a channel bottom that drops 1 m in a
length of 20 m has a slope of:
𝑑𝑟𝑜𝑝 1𝑚 𝑚
𝑆𝑜 = = = 0.05 = 5%.
𝑙𝑒𝑛𝑔𝑡ℎ 20𝑚 𝑚
Water surface slope (sw): drop of water surface divided by the horizontal run of the
bottom.
Slope of energy grade line (s): drop of energy grade line divided by the horizontal run
of the bottom.
The depth of flow (D)or (y): the vertical distance from the bottom of the channel to
the water surface.
Lining of the channel: the material used to form a channel. Channels are constructed
of many different materials, ranging from grass, lined soil to stone or concrete.
- The type of material forming the channel affects the flow because it increases or
decreases velocity at the interface between the water and the channel lining.
- The effect of lining on the velocity of the channel will affect the cross-section area
and the depth of flow.
- The maximum velocity within the section of flow occurs just below the free surface
in the center of the channel and can reach a magnitude of 2 to 2.5 times the
average velocity.
- The lowest velocity is located along the channel lining.
Figure 7.6: Velocity Distribution - open channel, the number indicates percentage of maximum velocity.
210
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝑻+𝑩
The cross sectional area: the area of a cross section of the flowing fluid: 𝑨 = ( )𝑫
𝟐
Top width (T): width of the channel section at the free surface.
Wetted perimeter (P): the distance along the channel cross section that is in contact
with flowing water.
P= sum of the distance ab, bc, and cd.
Hydraulic radius (Rh): the cross sectional area divided by the wetted perimeter
𝑨
[𝑹𝒉 = ] .
𝑷
𝑸
Average velocity (V): the discharge divided by the flow rate: [𝑽 = ].
𝑨
𝑨
Hydraulic depth (Dh): flow area divided by the top width [𝑫𝒉 = ].
𝑻
Side slope (m, or z): slope of the channel sides defined as vertical over horizontal.
Bottom width (B): width of the channel section at the bottom.
211
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
7.3 Analysis of open channel flow problems:
Open channel hydraulics is the study of the mechanics of water flowing in open
conduits, which generally include channels, streams, and pipes not flowing under pressure.
It should be indicated that the area of the intersection in open channel flow is a
function of depth y which is a function of flow discharge Q, unlike the pressurized flow where
the area of the flow is the same as the area of the pipe. This means that it is required to relate
the area to the depth of flow.
Many ideas can be made about the previous factors to give better understanding of
the relation required to relate them together.
In 1769, the French engineer Antoine Chezy developed the first uniform-flow formula.
This formula was derived based on two assumptions. The first assumption, Chezy assumed
that the force resisting the flow per unit area of the stream bed is proportional to the square
of the velocity (KV2), with K being a proportionality constant. The second assumption was
that the channel was undergoing uniform flow which requires that the depth of flow is
constant.
The Manning formula was introduced later. It is an empirical formula estimating the
average velocity of a liquid flowing in a conduit that does not completely enclose the liquid,
It was first presented by the French engineer Philippe Gauckler in 1867, and later re-
developed by the Irish engineer Robert Manning in 1890. more details will be shown in
section 7.
212
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
213
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Table 7.1 Characteristics of different cross sections Area Wetted Hydraulic Radius Top Width Hydraulic depth
Section (A) Perimeter (P) (Rh) (T) (Dh)
𝑏𝑦
𝑏𝑦 𝑏 + 2𝑦 𝑏 𝑦
𝑏 + 2𝑦
(𝑏 + 𝑚𝑦)𝑦
1 (𝑏 + 𝑚𝑦)𝑦 (𝑏 + 𝑚𝑦)𝑦
[𝑏𝑦 + (𝑚𝑦)𝑦 × 𝑏 + 2𝑦√1 + 𝑚2 𝑏 + 2𝑚𝑦
2 𝑏 + 2𝑦√1 + 𝑚2 𝑏 + 2𝑚𝑦
× 2]
𝑚𝑦
𝑚𝑦 2 2𝑦√1 + 𝑚2 2𝑚𝑦 𝑦/2
2√1 + 𝑚2
1
(2𝜃 − 𝑠𝑖𝑛2𝜃)𝑑𝑜2
8
If the level of water y is
less than the center of
the pipe:
(𝑑𝑜 /2) − 𝑦 1 sin 2𝜃 (sin 𝜃)𝑑0 or 1 2𝜃 − sin 2𝜃
𝜃 = cos−1 𝜃𝑑𝑜 (1 − ) 𝑑𝑜 ( ) 𝑑𝑜
𝑑𝑜 /2 4 2𝜃 2√𝑦(𝑑𝑜 − 𝑦) 8 sin 𝜃
If the leve yl is higher
than the center of the
pipe:
𝑦 − 𝑑𝑜 /2
𝜃 = 90𝑜 + sin−1
𝑑𝑜 /2
214
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Varied flow: if the water depth or discharge change along the length of the channel
can be: gradually varied flow or rapidly varied flow.
Gradually varied flow GVF: occurs when the change of flow depth is small;
and rapidly varied flow: is the case when the change of flow depth is large in
small distance, for example a wave, waterfall, or the rapid transition from a
stream channel into the inlet of a culvert.
215
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Based on the time criterion, open channel flow may be classified into two
categories:
Steady flow: the velocity and water depth at any section do not change with time.
They might be different from one point to another but the whole picture is fixed
with time.
Unsteady flow: the velocity and water depth at any section change with time.
(unsteady uniform flows are very rare in nature).
Summing all the components of these forces in the direction of flow, given
that the case is uniform flow
(∑𝐹 = 0) no acceleration:
∑𝐹 = 0;
𝐹1 + 𝑊𝑠𝑖𝑛() − 𝐹2 − 𝐹𝑓 = 0
𝐹1 = 𝐹2 (Uniform Flow), no change in water depth.
The total weight of the water body is:
𝑊 = 𝐴 𝐿 , (A: cross sectional area, L: length of reach),
Small slope 𝑠𝑖𝑛 = 𝑡𝑎𝑛 = 𝑆𝑜 𝑊𝑠𝑖𝑛() = 𝐴 𝐿 𝑆𝑜 .
The resistance force (Ff ) exerted by the channel boundaries may be
expressed in terms of resistance force per unit area (shear stress)
multiplied by the total channel bed area that is in contact with the flowing
water.
Channel contact area is the product of the wetted perimeter (P) and the
length of the channel reach (L).
In 1769, a French engineer , Antonie Chezy assumed that the shear stress
is proportional to the square of mean velocity
0 V2 the total resistance force Ff = 0 PL = k v2 PL, (P:wetted perimeter,
0 resistance force of channel bed per unit area)
217
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝐹1 + 𝑤 𝑠𝑖𝑛 𝜃 − 𝐹2 − 𝐹𝑓 = 0 ,
𝛾 𝐴
𝑉 = √(( ) 𝑆0 ( )) = 𝑉 = 𝐶 √𝑅ℎ 𝑆𝑒 (7.1)
𝐾 𝑃
Which is known as Chezy’s formula for open channel flow with C =Chezy’s resistance
factor.
2 1
1 3 2
𝑉 = 𝑅ℎ 𝑆𝑒 (S.I units) (7.2)
𝑛
2 1
1 3 2
𝑄 = 𝑅ℎ 𝐴 𝑆𝑒 (7.3)
𝑛
𝑚
Where,[𝑉 ( 𝑠 ) , 𝑅ℎ (𝑚), 𝑄 (𝑚3 /𝑠), 𝐴𝑟𝑒𝑎 (𝐴)(𝑚2 )]
2 1
1.49 3 2
𝑉= 𝑅ℎ 𝑆𝑒 ,(English units) (7.4)
𝑛
2 1
1.49 3 2
𝑄= 𝑅ℎ 𝐴 𝑆𝑒 , (English units) (7.5)
𝑛
𝑓𝑡
Where, [𝑉 ( 𝑠 ) , 𝑅ℎ (𝑓𝑡), 𝑄 (𝑓𝑡 3 /𝑠), 𝐴𝑟𝑒𝑎 (𝐴)(𝑓𝑡 2 )].
218
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Table 7.2 Typical values of Manning n Compiled from tables presented by Chow 1959].
Material n
Metal
Steel 0.012
Non-metal
Lucite 0.009
Glass 0.01
Cement 0.011
Concrete 0.013
Wood 0.012
Clay 0.013
Brickwork 0.013
Gunite 0.019
Masonary 0.025
Natural streams
219
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Rough 0.016
220
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Channels not maintained, weeds and brush Dense weeds, high as flow depth 0.08 – 0.12
uncut
Clean bottom, brush on sides 0.05 – 0.08
Clean bottom, brush on sides, highest stage of flow 0.07 – 0.11
Gutters Manning's n
221
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Rough 0.015
For gutters with small slope, where sediment may accumulate, increase all above values of 0.002
n by
Natural streams - minor streams (top width at floodstage < 100 ft)
1. Main Channels
a. clean, straight, full stage, no rifts or deep pools 0.025 0.030 0.033
b. same as above, but more stones and weeds 0.030 0.035 0.040
c. clean, winding, some pools and shoals 0.033 0.040 0.045
d. same as above, but some weeds and stones 0.035 0.045 0.050
e. same as above, lower stages, more ineffective
0.040 0.048 0.055
slopes and sections
f. same as "d" with more stones 0.045 0.050 0.060
g. sluggish reaches, weedy, deep pools 0.050 0.070 0.080
h. very weedy reaches, deep pools, or floodways
0.075 0.100 0.150
with heavy stand of timber and underbrush
2. Mountain streams, no vegetation in channel, banks usually steep, trees and brush along banks
submerged at high stages
a. bottom: gravels, cobbles, and few boulders 0.030 0.040 0.050
b. bottom: cobbles with large boulders 0.040 0.050 0.070
3. Floodplains
a. Pasture, no brush
1.short grass 0.025 0.030 0.035
2. high grass 0.030 0.035 0.050
b. Cultivated areas
1. no crop 0.020 0.030 0.040
2. mature row crops 0.025 0.035 0.045
3. mature field crops 0.030 0.040 0.050
c. Brush
1. scattered brush, heavy weeds 0.035 0.050 0.070
2. light brush and trees, in winter 0.035 0.050 0.060
3. light brush and trees, in summer 0.040 0.060 0.080
4. medium to dense brush, in winter 0.045 0.070 0.110
222
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
224
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
225
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
226
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 7.2:
a 3-m-wide rectangular channel carries a discharge of 25.3 m3/s at a uniform depth of
1.2 m. Determine the slope of the channel if Manning’s n cofficient is n=0.022.
Solution:
For a rectangular channel, the wetted perimeter and the hydraulic radius are:
𝐴 = 𝑏𝑦 = (3)(1.2) = 3.6 𝑚2
𝑃 = 𝑏 + 2𝑦 = 5.4 𝑚
𝐴 3.6 2
𝑅ℎ = = = = 0.667 𝑚
𝑃 5.4 3
2 1 1
1 1 2/3
Given: 𝑄 = 𝑅ℎ 𝐴 𝑆𝑒 25.3 =
3 2
0.022
0.667 (3.6) 𝑆𝑒 𝑆𝑜 = 𝑆𝑒 = 0.041
2
𝑛
Example 7.3:
If the discharge of the channel in the previous example is increased to 40 m3/s, what is
the normal depth of flow?
Solution:
The slope, the channel width b, and Manning’s n value are the same. Changing the flow
will change only the depth of flow.
The geometric parameters are:
𝐴 = 𝑏𝑦 = (3)𝑦
𝑃 = 𝑏 + 2𝑦 = 3 + 2𝑦
𝐴 3𝑦
𝑅ℎ = =
𝑃 3 + 2𝑦
Substituting these values into the Manning’s equation with n = 0.022 and So = 0.041:
2 1
1 1 3𝑦
𝑄 = 𝑅ℎ 𝐴 𝑆𝑒 40 =
3 2
(3𝑦)( )2/3 (0.041)1/2
𝑛 0.022 3+2𝑦
2
3𝑦 3 (0.022)(40)
(3𝑦) ( ) = = 4.346
3+2𝑦 (0.041)1/2
To solve this problem, it is needed to assume successive values for Y and substitute into
the last equation until the left side is significantly close to the value (4.346) :
2
Y (3𝑦)(3𝑦/(3 + 2𝑦))3 %error
227
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
The solution can be made by Excel to get more accurate value. It can be shown that
2
Y=1.68889 m will give (3𝑦)(3𝑦/(3 + 2𝑦)) = 4.346 with a percentage error
3
0.00022%.
Another approach can be used to solve this question by ploting Q versus Y and use the
plot to find the values of Y that corresponds to Q=40 m3/s.
228
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 7.4:
Find the flow in the following figure if the slope of the channel is 0.0004 and n=0.024.
Side slopes are all 1:1.
A=(3)(3)/2+(3)(20)+(5)(3+3+5)/2+(8)(3+5)+(5)(3+3+5)/2+(3)(25)+(3)(3)/2=263.0 m2
𝑃 = √32 + 32 + 20 + √52 + 52 + 8 + √52 + 52 + 25 + √32 + 32 = 75.63𝑚
R=A/P = 263 / 75.63 = 3.477 m
2
1 1/2 1
𝑄= 𝐴𝑅2/3 𝑆𝑜 = 263.0 (3.477) (0.0004)0.5 = 503 𝑚3 /𝑠
3
𝑛 0.024
229
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Figure 7.14: Normal depth solution procedure; (a) Trapezoidal channels (m=side slops) and (b)
circular channels (do=diameter)
230
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Tutorial 9
Q1. Calculate P, A, and R for the following channel section.
Sol:
231
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Figure 7.15: velocity distribution in a circular pipe flowing part flow. Values are multiple of the
average velocity
Example7.5:
a 6-ft-diameter, concrete pipe is
flowing with a free surface (not under
pressure). If the pipe is laid on a slope
of 0.001 and carries a uniform flow at a
depth of 4 feet, what is the discharge?
Solution:
𝑦−𝑑𝑜 /2
𝜃 = 90𝑜 + sin−1
𝑑𝑜 /2
𝑜 −1
4−3 𝑜
109.5𝑜
𝜃 = 90 + sin = 109.5 = (2𝜋) = 0.608𝜋 𝑟𝑎𝑑𝑖𝑎𝑛𝑠
3 360𝑜
The area of the circular section is:
1 1
𝐴 = (2𝜃 − sin 2𝜃)𝑑𝑜2 = [2(0.608𝜋) − sin 2(0.608𝜋)](6𝑓𝑡)2 = 20.0 𝑓𝑡 2
8 8
𝑃 = 𝜃𝑑𝑜 = (0.608𝜋)(6𝑓𝑡) = 11.5 𝑓𝑡
Then the hydraulic radius is
𝐴 20 𝑓𝑡 2
𝑅ℎ = = = 1.74𝑓𝑡
𝑃 11.5𝑓𝑡
2 1 2
1.49 1.49
using n=0.013 𝑄 = 𝑅ℎ3 𝐴 𝑆𝑒2 = (20)(1.74)3 (0.001)0.5 = 105 𝑓𝑡 3 /𝑠
𝑛 0.013
232
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 7.6:
Find the quantity and velocity of flow in an 18-inch-diameter concrete pipe flowing one
quarter full. The slope of the pipe is 1.00 percent and the Manning’s n is 0.015.
Solution:
(𝑑𝑜 /2) − 𝑦
𝜃 = cos −1
𝑑𝑜 /2
cos −1 (18/2) − 4.5
𝜃= = 60𝑜 = 𝜋/3
18/2
2
1 1 18
𝐴 = 2𝜃 − sin 2𝜃 𝑑𝑜 = 2 𝜋/3 − sin 2 𝜋/3 ( 𝑓𝑡) = 0.345 𝑓𝑡 2
( ) 2 [ ( ) ( )]
8 8 12
18
𝑃 = 𝜃𝑑𝑜 = (𝜋/3) ( 𝑓𝑡) = 1.57 𝑓𝑡
12
𝐴 0.345 𝑓𝑡 2
𝑅ℎ = = = 0.220𝑓𝑡
𝑃 1.57𝑓𝑡
1.49 23 12 1.49 2
𝑄= 𝑅ℎ 𝐴 𝑆𝑒 = (0.345)(0.220)3 (0.01)0.5 = 1.25 𝑓𝑡 3 /𝑠
𝑛 0.015
𝑄 1.25
𝑉= = = 3.62 𝑓𝑡/𝑠
𝐴 0.345
𝑦 − 𝑑𝑜 /2
(𝑑𝑜 /2) − 𝑦 𝜃 = 90𝑜 + sin−1
𝜃 = cos −1 𝑑𝑜 /2
𝑑𝑜 /2
𝜋 𝑦 − 𝑑𝑜 /2
𝜃= + sin−1
2 𝑑𝑜 /2
233
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Given A is constant , what is the sction that will give the maximum Q?
𝑨
It can be easily concluded that ( ) is maximum, or P is minimum (since A is constant).
𝑷
𝒅𝑷
For P to be minimum, [ = 𝟎].
𝒅𝒚
𝑦
𝑅= (7.7)
2
𝑏 = 2𝑦 (7.8)
234
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝐴 = (𝑏 + 𝑚𝑦)𝑦
𝐴
= 𝑏 + 𝑚𝑦
𝑦
𝐴
𝑏 = − 𝑚𝑦
𝑦
Figure 7.16 Cross-section in a trapezoidal channel
𝑃 = 𝑏 + 2𝑦√1 + 𝑚2
𝐴
𝑃 = ( − 𝑚𝑦) + 2𝑦√1 + 𝑚2
𝑦
For minimum P:
𝑑𝑃 𝐴
= 0 = (− 2 − 𝑚) + 2√1 + 𝑚2
𝑑𝑦 𝑦
𝐴
( + 𝑚) = 2√1 + 𝑚2 ,
𝑦2
(𝑏+𝑚𝑦)𝑦
substituting {𝐴 = (𝑏 + 𝑚𝑦)𝑦} in the last equation: ( + 𝑚) = 2√1 + 𝑚2
𝑦2
𝑏 + 2𝑚𝑦 = 2𝑦√1 + 𝑚2
𝑏 = 2𝑦[√1 + 𝑚2 − 𝑚] (7.9)
235
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝑦
𝑅=
2
𝑏 = 2𝑦 [√1 + 𝑚2 − 𝑚]
Among all open-channel shapes, the semicircle has the least perimeter for a given area, so
it is the most hydraulically efficient of all sections. A channel with a semi-circular cross
section, however, has sides that are curved and almost vertical at the water surface
level, which makes the channel expensive to construct (excavation and forming) and
difficult to maintain (bank stability). In practice, semi-circular sections are only used
when pipes are appropriate or in smaller flumes of prefabricated materials.
236
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝑽𝟐
𝑬=𝒚+ (7.11)
𝟐𝒈
Figure 7.18 Specific energy curves for a given channel section and specific discharge
237
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
It can be noted from the plot of the flow depth versus the specific energy (Fig. 7.18)
that:
For a very low velocity with very high depth y, the specific energy E approaches the
value of the depth y (most of the energy is composed of flow depth).
For high velocities , the values of the depths y approach zero (most of the specific
energy is composed of velocity head).
Each value of the specific energy E is corresponding into two values of y,except the
minimum value of specific energy which correponds into one value of depth.
The value of depth y at the minimum specific energy is called critical depth yc. It is
a theoretical concept that depends on the channel shape and the flow quantity
only. Yc does not depend on the channel slope So or roughness.
When the flow depth in the channel is greater than the critical depth yc then the
case represents more tranquil flow called subcritical flow. While when the depth
in the channel is below the critical depth yc then the case is supercritical flow,
which is representing a more rapid case of flow. The case when the depth of flow
concide with the critical depth is critical flow. The velocity of the water at critical
depth is called critical velocity and the slope is critical slope.
The critical depth can be computed using the following relationship which occurs
𝑑𝐸
only at minimum specific energy( = 0 ), (valid for gradually varied flow):
𝑑𝑦
𝑨𝟑 𝑸𝟐
= (7.12)
𝑻 𝒈
Where, A is the cross-sectional arae of the chnnel (m2), T is the top width
of the channel (m), and Q is the flow rate (m3/s).
238
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝑇= 𝑏
𝐴 = 𝑏𝑦
𝐴3 𝑏3 𝑦 3 𝑄2
= =
𝑇 𝑏 𝑔
3 𝑄2
𝑦𝑐 = √ (7.13)
𝑏2 𝑔
Only one particular slope result in a normal depth equal to critical depth.
S= Sc yn = yc
S Sc supercritical flow.
S Sc subcritical flow.
The most common occurrence of critical depth is gradually varied flow, in which
the flow depth varies over a certain distance due to a change in one or more of the
channel attributes. The change could be the channel slope or cross sectional shape.
As the depth varies, yc can be attained within the transition region.
239
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Figure 7.19: Critical depth solution procedure: (a) trapezoidal channels and (b) circular channels
(Fundamentals of Hydraulic Engineering Systems Ned Hwang)
240
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 7.7:
A trapezoidal channel has a bottom width of 5 m and side slopes m=2. If the flow rate is
20 m3/s, what is the critical depth?
Solution:
𝑄2 𝐴3 [(𝑏+𝑚𝑦)𝑦]3
By using the equation = =
𝑔 𝑇 𝑏+2𝑚𝑦
Alternative solution using the figure for critical depth of trapezoidal channels:
𝑄𝑚3/2 20(2)3/2
= = 0.323
𝑔1/2 𝑏 5/2 (9.81)1/2 (5)5/2
𝑚𝑦𝑐 (0.41)(5)
From Fig. 7.19, we obtain = 0.41. Therefore, 𝑦𝑐 = = 1.03𝑚.
𝑏 2
The previous figure shows a profile of channel transition from subcritical flow to
supercritical flow. Depth of flow is at a critical depth near the transition point.
241
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Theoretical critical depth does not change when the channel slope change.
𝑉
𝐹𝑟 = (7.14)
√𝑔𝑦
Where, F is the Froude number, V is the average velocity (m/s), and y is the flow
depth.
When Fr = 1, flow is critical
Fr 1, flow is subcritical
Fr 1, flow is supercritical.
For non-rectangular channels:
The Froude number is defined as:
𝑉
𝐹= (7.15)
√𝑔𝑦ℎ
𝑨
Where yh is the hydraulic depth [𝒚𝒉 = ], A is the cross sectional area and T is the top
𝑻
width.
Example 7.8:
Determine the critical depth and flow velocity in a 10 m wide rectangular channel
carrying 50 m3/s.
Solution:
𝑄2 33 502
𝑦𝑐 = √ 2 = √ = 1.37𝑚
𝑔𝑏 9.81(10)2
242
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Given the discharge Q in a particular channel, the amount of energy head loss
through a jump (∆E) can be determined by measuring the initial and sequent depths and
using the specific energy curves.
Predicting the sequent depth by estimating the energy loss is impractical because
it is difficult to determine the energy loss through the jump.
The relationship between the initial depth and the sequent depth in a hydraulic
jump may be determined by considering the balance of forces and momentum before
and after the jump.
Consider a control volume in the vicinity of a jump as shown:
The balance between the hydro static
forces and the momentum flux through
sections 1 and section 2 per unit width of
the channel is:
𝐹1 − 𝐹2 = 𝜌𝑞(𝑉1 − 𝑉2 )
𝛾 𝛾 𝑞 𝑞 𝑞2 (𝑦1 +𝑦2 )
𝐹1 = 𝑦12 , 𝐹2 = 𝑦22 , 𝑉1 = , 𝑉2 = = 𝑦1 𝑦2
2 2 𝑦1 𝑦2 𝑔 2
243
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝑦2 1 𝑦1 1
= (√1 + 8𝐹𝑟21 − 1) Or = (√1 + 8𝐹𝑟22 − 1) (7.16)
𝑦1 2 𝑦2 2
𝑉1
where, 𝐹𝑟1 = , is the Froude number of the flow.
√𝑔 𝑦1
Example 7.9:
A 10 ft-wide rectangular channel carries 500 cfs of water at a 2-ft depth before
entering a jump. Compute the downstream water depth and the critical depth.
Solution:
𝑄2 3 3 5002
𝑦𝑐 = √ 2 = √ = 4.27𝑓𝑡
𝑔𝑏 9.81(10)2
𝑄 500 𝑓𝑡
𝑉1 =
= = 25
𝑏𝑦 10 × 2 𝑠
𝑉1 25
𝐹𝑟1 = = = 3.12
√𝑔 𝑦1 √32.2 × 2
𝑦2 1
= (√1 + 8𝐹𝑟21 − 1)
𝑦1 2
𝑦2 1
= (√1 + 8(3.12)2 − 1)
2.0 2
𝑦2 = 7.88 𝑓𝑡
{𝐹𝑠 = 𝐹 + 𝜌𝑞} is known as specific force per unit width of the channel.
244
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
For a given discharge, the specific force is a function of the water depth at a given
section.
Plotting Fs vesus water depth will give us a curve similar to the specific energy curve
with a vertex that appears at the critical depth.
A hydraulic jump takes place in a short reach therefore it is reasonable to assume
that through a hydraulic jump, the specific forces immediately before and after the jump
are the same.
The energy head loss through the hydraulic jump ∆E may be estimated by:
V12 V22 1 2 2
q2 1 1
∆E = ( + y1 ) − ( + y2 ) = (V1 − V2 ) + (y1 + y2 ) = ( 2 − 2 ) + (y1 − y2 )
2g 2g 2 2g y1 y2
𝑞2 𝑦1 +𝑦2
Given that
𝑔
= 𝑦1 𝑦2 (
2
)
(y2 −y1 )3
∆E = (7.20)
4y1 y2
245
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
246
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 7.10:
A rectangular channel has a width of b = 3 m, a discharge of 4 m3/s, a bed slope
So=0.001. The channel is lined with Concrete; bottom float-finished, sides are made of
Random stone in mortar. Find the following:
a) The Normal depth yn = ……1.068 m……....
b) The critical depth yc =……..0.566 m…….…
c) The state of flow is: (a) sub-critical, (b) critical, (c) super-critical flow.
d) The minimum specific energy is……0.849 m…….
e) Given that the initial depth in a hydraulic jump is y = 0.1 m, then the sequent
depth is ………1.854 m………….
Solution:
Manning n = 0.018 (Table).
5
1 (3𝑦) ⁄3
(a) 𝑄=4= 0.0010.5 y = 1.068 m
0.0185 (3+2𝑦)2/3
3 42
(b) 𝑦𝑐 = √ = 0.566 𝑚
9.81 (3)2
Example 7.11:
Design a trapezoidal channel having sides slopes with z =1.5, a discharge of 7 m3/s, with a
velocity of 1.8 m/s, and a smooth asphalt lining:
y = .............1.36m......, b =..............0.822m............, S o.......9.1610-4 (0.000916).
Solution:
Unknowns: y, b, S0
Equations:
V=1.8 m/s = (1/0.013) R2/3 S01/2 ........ (1)
A=7/1.8 = ( b + 1.5 y) y ........(2)
Three Unknowns & 2 equations Use most efficient section
𝑏 = 2𝑦 [√1 + 1.52 − 1.5 ] ⇒ 𝑏 = 0.606𝑦
R = y /2
Substitute 𝑏 = 0.606𝑦 in (2) y =1.359 m
R = 1.359 / 2 = 0.6795 substitute in (1) S0 = 0.000916
247
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 7.12:
In a rectangular channel 4 m wide with a flow of 7.53 m3/s the depth is 0.3 m. If a
hydraulic jump occurs, find the depth downstream of it. Calculate the loss of specific
energy through the jump.
Solution:
7.53 2
𝑦1 8𝑞 2 0.3 8( )
𝑦2 = ( ) (√1 + 3 − 1) = ( ) (√1 + 4
− 1) = 1.410 𝑚
2 𝑔𝑦 2 1 9.81(0.3)3
Example 7.13:
A hydraulic jump occurs in a 4-m-wide rrectangular channel carrying 5 m3/s on a slope
of 0.004. The depth after the jump is 1.2 m. Find (a) the depth before the jump, (b) the
losses of energy and power in the jump.
5 2
𝑦2 8𝑞 2 1.2 8( )
𝑦1 = ( ) (√1 + 3 − 1) = ( ) (√1 + − 1) = 0.191 𝑚 4
2 𝑔𝑦 2 9.81(1.2)3
2
248
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example 7.14:
Solution:
The specific energy curve can be constructed using the following relationship:
𝑄2 (15)2 1.27
𝐸 =𝑦+ = 𝑦 + = 𝑦 +
2𝑔(6𝑦)2 2(9.81)(3𝑦)2 𝑦2
3 52
𝑦𝑐𝑟 = √ = 1.37 𝑚 , 𝐸𝑐𝑟 = 2.05𝑚
𝑔
249
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
The total energy head at the inlet as measured with respect to the datum line is:
𝑉𝑖2
𝐻𝑖 = + 𝑦𝑖 + 𝑧𝑖 = 0.1 + 3.6 + 0.4 = 4.1 𝑚
2𝑔
The energy level is shown on the figure by the top horizontal line.
At the exit of the transition, the total energy available is reduced by 0.1 m as given
and as indicated by the EGL in the figure.
𝑉𝑒2
𝐻𝑒 = + 𝑦𝑒 + 𝑧𝑒 = 𝐻𝑖 − 0.1 = 4.0 𝑚
2𝑔
Ee is the specific energy measured with respect to the channel bottom:
Ee=He=4.00m
The value is applied to the specific energy curve to obtain the water depth at the
exit section.
Water surface elevations for four other sections are computed by using the same
method (4, 8, 12, and 16 m from the entrance section). The results for the six
sections are shown:
Section Inlet 4.00m 8.00m 12.00m 16.00m Exit
Specific Energy, E(m) 3.70 3.76 3.82 3.88 3.94 4.00
Water Depth, y(m) 3.60 3.67 3.73 3.79 3.86 3.92
Example 7.15:
A long, rectangular open channel 3m wide carries a discharge of 15m3/sec. The channel
slope is 0.004, and the Manning's coefficient is 0.01. At a certain point in the channel,
flow reaches normal depth.
(b) If a hydraulic jump takes place at this depth, what is the sequent depth?
Solution
3 𝑄2 3 𝑞2 3 (15)2
(a) The critical depth is calculated using the equation 𝑦𝑐 = √ =√ =√ 2=
𝑔𝑏2 𝑔 𝑔(3)
1.37𝑚. The normal depth of this channel can be determined by the Manning equation:
1 23 12
𝑄 = 𝑅ℎ 𝐴 𝑆𝑒
𝑛
250
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝐴1 𝑦1 𝑏
Where 𝐴 = 𝑦1 𝑏, 𝑅ℎ = = , 𝑏 = 3𝑚
𝑃1 2𝑦1 +𝑏
1 3𝑦1 2/3
15 = ( ) , 𝐴(3𝑦1 ) (0.004)1/2
0.01 2𝑦1 + 3
15 4.63𝑚
Solving the equation for y1, we obtain y1 = 1.08 m, 𝑉1 = = , and
3𝑦1 𝑠
𝑉1
𝐹1 = = 1.42
√𝑔𝑦1
𝑦 8𝑞 2
(b) Applying the equation of hydraulic jump: 𝑦2 = ( 1) (√1 + 3 − 1) =
2 𝑔𝑦 1
𝑦
( 21) (√1 + 8𝐹12 − 1) = 1.70𝑚.
(c) The head loss can be estimated by using the following equation:
(𝑦2 − 𝑦1 )3 (0.62)3
△𝐸 = = = 0.032𝑚
4𝑦1 𝑦2 4(1.70)(1.08)
251
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Tutorial. 10
Q1.A sewer line is circular in section has a diameter of 50 cm and is running 0.75 full.
It is laid at a slope of 1:900. If n = 0.02, find the discharge and the mean velocity of
flow.
Q2. Calculate the size of a circular sewer, which is to run 0.8 full, at a gradient of 1
in 250 maintaining a velocity of 1.5 m/s, has n = 0.018. Find also the discharge.
Q3. Prepare a monograph showing Geometric Elements of a circular section giving
relationship between Y/do and {A/Ao, P/Po and AR2/3/AoR02/3}.
Q4.An earthen channel having base width of 5 m and side slopes 1 vertical and 1.5
horizontal carries water with a depth of 1.6 m. The bed slope is 1 in 500. Calculate
the discharge if n = 0.018, what would be the Chezy's C to get the same discharge?
Q5. Calculate the normal depths for the following channels if Q = 6 m3/s, s = 0.0015,
and n = 0.021:
(i) Rectangular channel, b = 8 m
(ii) Trapezoidal channel, b = 7 m, and z = 1 m.
Q6. Design a trapezoidal channel having side slopes with z=1, Q =12 m3/s, s = 1 : 800
and n = 0.02.
Q7. Design a trapezoidal channel carrying water having a velocity of 1.2 m/s and a
discharge of 10 m3/s. n = 0.019 and z = 1.5.
Q8. Determine the most efficient section of a trapezoidal canal to carry discharge of
10 m3/s at a velocity of 1.5 m/s. The canal side slopes are 1V : 0.5H. find the
necessary bed slope also n = 0.02.
252
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Solution of tutorial 10
Q1. A sewer line is circular in section has a diameter of 50 cm and is running 0.75 full. It is
laid at a slope of 1:900. If n = 0.02, find the discharge and the mean velocity of flow.
So= 1/900
n = 0.02
Q=?
V=?
Y=0.7550=37.5 cm
37.5−50/2 2 37.5−50/2
= /2 + sin−1 = 𝜋 , = 90o + sin−1 = 120𝑜
50/2 3 50/2
1
𝐴𝑟𝑒𝑎 (𝐴) = 8(2×23𝜋−sin(2×120𝑜 )(0.5)2 = 0.157 𝑚2
2
𝑃 = 𝜃𝑑𝑜 = 3𝜋 ×0.5=1.0472 𝑚
⁄
1 2⁄ 0.5 1 0.157 2 3 1 0.5
𝑄= 𝐴𝑅 3 𝑆𝑜 = (0.157) ( ) ( )
𝑛 0.02 1.0472 900
Q2. Calculate the size of a circular sewer, which is to run 0.8 full, at a gradient of 1 in 250
maintaining a velocity of 1.5 m/s, has n = 0.018. Find also the discharge.
0.3 𝑑𝑜 0.3 𝑑𝑜
= /2 + sin−1 = 0.704𝜋 , = 90o + sin−1 = 126.87𝑜
0.5 𝑑𝑜 0.5 𝑑𝑜
1
𝐴𝑟𝑒𝑎 (𝐴) = 8(2×0.704𝜋−sin(2×126.87𝑜 )(𝑑𝑜)2 = 0.673(𝑑𝑜 )2
𝑃 = 𝜃𝑑𝑜 = 0.704𝜋𝑑𝑜
2⁄3
1 0.673𝑑02 1
𝑉 = 1.5 𝑚/𝑠 = 𝑛1 𝑅2⁄3𝑆𝑜0.5 = ( ) ( )0.5
0.018 0.704𝜋𝑑0 250
253
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Q4.
𝑄 = 𝐶𝐴√𝑅ℎ 𝑆𝑓
Q5. By trial
1 𝐴5/3 12
3
1 𝐴5/3
𝑄 = 6 𝑚 /𝑠 = 𝑆 = (0.0015)0.5
𝑛 𝑃2/3 𝑒 0.021 𝑃2/3
(i) Rectangular channel, b = 8 m
A =8 × Y ,P=8+2Y
(ii)
1 23 12
𝑄 = 𝑅ℎ 𝐴 𝑆𝑒
𝑛
2
𝑚3 𝐴 𝐴 3 1 0.5
12 = ( ) ( ) (1)
𝑠 0.02 𝑃 800
Two unknowns Y & b
One equation (Manning’s eq.)
→ Need another equation
𝒀
→ Use the most efficient section. {𝑹 = } , {𝒃 = 𝟐𝒀[√𝟏 + 𝒛𝟐 − 𝒛]}
𝟐
254
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝑚3 𝐴 2 1
12 = (𝑅ℎ )3 ( )0.5
𝑠 0.02 800
2
𝑚3 [(𝑏 + 1 × 𝑌) × 𝑌] 𝑌 3 1 0.5
12 = ( ) ( )
𝑠 0.02 2 800
2
𝑚3 [(0.828𝑌 + 𝑌) × 𝑌] 𝑌 3 1 0.5
12 = ( ) ( )
𝑠 0.02 2 800
Q7. Z=1.5 , {V=1.2 m/s , Q=10 m3/s }
𝒀
→ use the most economical section {𝑹 = } , {𝒃 = 𝟐𝒀[√𝟏 + 𝒛𝟐 − 𝒛]}
𝟐
255
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝒀
→ use the most economical section {𝑹 = } , {𝒃 = 𝟐𝒀[√𝟏 + 𝒛𝟐 − 𝒛]}
𝟐
Two equations
𝑄 6 6 1
𝐴= = = (𝑏 + 𝑌 ) × 𝑌 →𝑏= × − 𝑌 (1)
𝑉 1.3 1.3 𝑌
𝑚 1
𝑉 = 1.3 = 𝑅 2/3 (0.002)1/2 (2)
𝑠 0.019
Q10.
5
1 (6𝑦) ⁄3 (6𝑦)5/3
I. 𝑄=8= 0.0020.5 = 3.578 → By trial {y = 0.806 m}
0.02 (6+2𝑦) (6+2𝑦)2/3
3 82
II. 𝑦𝑐 = √ = 0.566 𝑚
9.81 (6)2
Figure-
In uniform flow, the water depth remains a constant value known as the normal. The
energy grade line is parallel to the water surface and the channel bottom. The velocity
distribution also remains unchanged throughout the reach. The computation of one
water depth is sufficient for the entire channel reach. In rapidly varied flow such as a
hydraulic jump, rapid changes in water take place in a short distance. A significant change
in water velocities is associated with the rapid variation of water cross-sectional area. At
this high rate of flow deceleration, the energy loss is very high. The computation of water
depths using the energy principles is not reliable. In this case, computations can only be
carried out applying the momentum principles (𝐹1 − 𝐹2 = 𝜌𝑞(𝑉2 − 𝑉1 ). In gradually
varied flow, velocity changes take place very gradually with distance so that the effects of
acceleration on the flow between two adjacent sections are negligible. Thus, computation
of the water surface profile defined as depth changes along the channel length can be
carried out on energy considerations.
The total energy head at any section in an open channel can be given as:
𝑉2 𝑄2
𝐻= +𝑦+𝑧 = +𝑦+𝑧
2𝑔 2𝑔𝐴2
To compute the water surface profile, we must first obtain the variation of the total
energy head along the channel. Differentiating H with respect to the channel distance x,
we obtain the energy gradient in the direction of the flow:
257
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝑑𝐻 𝑄2 𝑑𝐴 𝑑𝑦 𝑑𝑧 𝑄2 𝑇 𝑑𝑦 𝑑𝑦 𝑑𝑧
=− 3 + + =− 3 + +
𝑑𝑥 𝑔𝐴 𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑔𝐴 𝑑𝑥 𝑑𝑥 𝑑𝑥
Where 𝑑𝐴 = 𝑇(𝑑𝑦). Rearranging the equation gives:
𝑑𝐻 𝑑𝑧
𝑑𝑦 𝑑𝑥 − 𝑑𝑥
=
𝑑𝑥 𝑄2 𝑇
1−
𝑔𝐴3
The term dH/dx is the slope of the energy grade line. It is always a negative quantity
because the total energy head reduces in the direction of the flow, or Se = -dH/dx.
Similarly, the term dz/dx is the slope of the channel bed. It is negative when the elevation
of the channel bed reduces in the direction of the flow; it is positive when the elevation
of the channel bed increases in the direction of the flow. In general, we may write S0 = -
dz/dx. The energy slope in gradually varied flow between two adjacent sections may also
be approximated by using a uniform flow formula. For simplicity, the derivation will be
demonstrated with a wide rectangular channel section where A = by, Q = bq, and Rh = A/P
= by/(b + 2y) == y (for wide rectangular channels because b >> y). Using the Manning
formula we get:
𝑑𝐻 𝑛2 𝑄2 𝑛2 𝑄 2
𝑆𝑒 = − = 4 = 2 10/3
𝑑𝑥 𝑏 𝑦
𝑅ℎ3 𝐴2
The slope of the channel bed may also be expressed in similar terms if uniform flow were
assumed to take place in the channel. Because the slope of the channel bed is equal to
the energy slope in uniform flow, the hypothetical uniform flow conditions are designated
with the subscript n. We have
𝑑𝑧 𝑛2 𝑄 2
𝑆𝑜 = − = ( 2 10/3 )𝑛
𝑑𝑥 𝑏 𝑦
3 𝑄2 3 𝑞2
For rectangular channels, we have 𝑦𝑐 = √ =√
𝑔𝑏2 𝑔
Or
2
𝑔𝐴3𝑐
𝑄 = 𝑔𝑦𝑐3 𝑏 2 =
𝑏
Rearranging the equations, we get:
𝑦𝑛 10/3
𝑑𝑦 𝑆𝑜 [1 − ( 𝑦 )
=
𝑑𝑥 𝑦 3
[1 − ( 𝑐 ) ]
𝑦
258
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝑦𝑛 𝑁
𝑆 [1 − ( ]
𝑑𝑦 𝑜 𝑦)
=
𝑑𝑥 𝑦 𝑀
[1 − ( 𝑐 ) ]
𝑦
where the exponents M and N depend on the cross-sectional shape and the flow
conditions as given by V. T. Chow.
This form of the gradually varied flow equation is very useful for a qualitative analysis,
which helps to understand the gradually varied flow classifications covered in the next
sections. Other forms are often used to compute water surface profiles. Physically, the
term dy/dx represents the slope of the water surface with respect to the bottom of the
channel. For dy/dx = 0, the water depth remains constant throughout the reach or the
special case of uniform flow. For dy/dx < 0, the water depth decreases in the direction of
the flow. For dy/dx > 0, the water depth increases in the direction of the flow. Solutions
of this equation under different conditions will yield the various water surface profiles
that occur in open channels.
Another Approach:
The dynamic equation of GVF can be
obtained by differentiating the energy
equation with respect to the channel
bed (x-direction):H= z + y cos + V2/2g
2
𝑑𝐻 𝑑𝑧 𝑑𝑦 𝑑(𝑉 ⁄2𝑔)
= 𝑑𝑥 + cos 𝜃 +
𝑑𝑥 𝑑𝑥 𝑑𝑥
𝑑𝐻 𝑑𝑧
𝑆𝑓 = − 𝑑𝑥 ; 𝑆𝑜 = sin 𝜃 = − 𝑑𝑥
𝑑𝑦 𝑆𝑜 −𝑆𝑓
= 2
𝑑𝑥 𝑑(𝑉 ⁄2𝑔)
cos 𝜃+
𝑑𝑦
𝑉2 𝑄2
𝑑( ⁄2𝑔) 𝑑( ⁄ ) 2𝑄 2 𝑑𝐴
2𝑔𝐴2
= = − 2𝑔𝐴3 𝑑𝑦
𝑑𝑦 𝑑𝑦
𝑑𝐴
= 𝑇, width of the liquid surface.
𝑑𝑦
2
𝑑(𝑉 ⁄2𝑔) 𝑄2 𝑇
= − 𝑔𝐴3
𝑑𝑦
𝑑𝑦 𝑆𝑜 −𝑆𝑓
For small channel slope, cos = 1.0 = 1−𝑄2 𝑇/(𝐴3 𝑔) , which gives the slope of the water
𝑑𝑥
surface relative to the bed.
Various types of gradually varied flow surface profile can occur. A very important factor is the type
of bed slope, mild, steep, critical, horizontal, or adverse. And the depth of flow relative to the critical
flow and depending on the type of bed slope (mild, steep, ..).
259
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
260
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
3. If the water depth is less than both Yc and Yn' then the water surface profile curves
are type 3. There are S-3, C-3, M-3, H-3, and A-3 curves. Each of these water
surface profile curves is listed and shown schematically in the following figure.
Examples of physical occurrences in open channels are also given. Certain important
characteristics of water surface profile curves can be demonstrated from direct analysis
261
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
𝑦
𝑑𝑦 𝑆𝑜 [1−( 𝑛 )10/3
𝑦
of the gradually varied flow equation ( = 𝑦 3
). By making substitutions into
𝑑𝑥 [1−( 𝑐 ) ]
𝑦
this equation, we note the following:
1. For type 1 curves, y/yc > 1 and y/yn > l. Thus, the value of dy/dx is positive,
indicating that water depth increases in the direction of the flow.
2. For type 2 curves, the value of dy/dx is negative. The water depth decreases in the
direction of the flow.
3. For type 3 curves, the value of dy/dx is again positive. The water depth increases
in the direction of the flow.
Characteristics of water surface profile curves
4. When the actual water depth approaches critical depth, y = yc , the gradually varied
𝑦
𝑑𝑦 𝑆𝑜 [1−( 𝑛 )10/3 𝑑𝑦
𝑦
flow equation ( = 𝑦 3
) yields = ∞, indicating that the slope of the
𝑑𝑥 [1−( 𝑐 ) ] 𝑑𝑥
𝑦
water surface profile curve is theoretically vertical. Likewise, as y approaches yn ,
dy/dx approaches zero, indicating that the water surface profile approaches the
normal depth line asymptotically.
5. A few types of water surface profile curves never approach a horizontal line (S-2,
S-3, M-2, M-3, C-3, H-3, and A-3). Others approach a horizontal line asymptotically,
except for the C-1 curve, which is horizontal throughout the channel reach.
Because Yn = Yc in a critical channel, the gradually varied flow equation yields dy/dx
= SO, indicating the water depth increases at the same rate as the channel bed
elevation decreases, which theoretically, results in a horizontal water surface
profile.
In channels where y < yc, the velocity of water flow is greater than that of the disturbance
wave. For this reason, flow conditions in the downstream channel will not affect those
upstream. The change of water depth resulting from any channel disturbance propagates
in the downstream direction. Thus, computation of the water surface profile should be
carried out in the downstream direction S-2, S-3, C-3, H-3, and A-3).
262
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
In channels where y > yn the speed of wave propagation is greater than the velocity of
water flow. Any disturbance in the downstream channel can travel upstream and affect
the flow conditions upstream as well as downstream. Any change of water depth in the
downstream channel propagates upstream and may also change the water depth in the
upstream channel. Thus, computation of the water surface profile should be carried out
in the upstream direction (M-1, M-2, S-1, C-1, H-2, and A-2). At the break of a channel
from mild to steep slope or a significant drop of the channel bottom, critical depth is
assumed to take place in the immediate vicinity of the brink. At this point, a definite
depth-discharge relationship can be obtained (i.e., control section) and is frequently used
as a starting point for water surface profile computations.
263
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
the next section, either upstream or downstream from the control section. The distance
between sections is critical because the water surface will be represented by a straight
line. Thus, if the depth of flow is changing quickly over short distances, adjacent sections
should be closely spaced to represent accurately the water surface profile. The step-by
step procedure is carried out in the downstream direction for rapid (supercritical) flows
and in the upstream direction for tranquil (subcritical) flows.
Standard Step Method
The standard step method is presented to calculate gradually varied flow water surface
profiles. The method employs a finite difference solution scheme to solve the differential,
𝑦
𝑑𝑦 𝑆𝑜 [1−( 𝑛 )10/3
𝑦
gradually varied flow equation ( = 𝑦 3
). It is the most common algorithm used
𝑑𝑥 [1−( 𝑐 ) ]
𝑦
in computer software packages that solve gradually varied flow profiles. For example, it
is the primary algorithm in the widely used HEC-RAS program developed by the U.S. Army
Corps of Engineers.
The standard step method is derived directly from an energy balance between two
adjacent cross sections (Figure) that are separated by a sufficiently short distance so that
the water surface can be approximated a straight line. The energy relation between the
two sections may be written as:
𝑣22 𝑣12
+ 𝑦2 +△ 𝑧 = + 𝑦1 + ℎ𝐿
2𝑔 2𝑔
where △ 𝑧 is the elevation difference in the channel bottom and is the energy head loss
between the two sections, as shown in the figure. or the equation may be rewritten as:
𝑣22 𝑣12
+ 𝑦2 + 𝑧2 = + 𝑦1 + 𝑧1 + 𝑆𝑒 △ 𝐿
2𝑔 2𝑔
𝐸2′ = 𝐸1′ + 𝑙𝑜𝑠𝑠𝑒𝑠
Where z is the position head (channel bottom elevation with respect to some datum) and
E' is the total energy head (position + depth + velocity). It is important to note that, in
264
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
previous energy equations, the sections 1 and 2 represent downstream and upstream
sections, respectively. If the sections are numbered differently, the losses should always
be added to the downstream side. The computation procedure yields the correct depth
at a cross section that is a distance △ 𝐿 away from a section with a known depth.
Computations begin at a control section and progress upstream (subcritical flow) or
downstream (supercritical flow). For subcritical flow, the water surface profile is
occasionally called a backwater curve because the process moves from downstream to
upstream. Likewise, the profile for supercritical flow is occasionally called a front-water
𝑣22 𝑣12
curve. Equation ( + 𝑦2 + 𝑧2 = + 𝑦1 + 𝑧1 + 𝑆𝑒 △ 𝐿 ) cannot be solved directly for the
2𝑔 2𝑔
unknown depth (e.g., y2) because v2 and Se depend on y2. Therefore, an iterative
procedure is required using successive approximations of y2 until the downstream and
upstream energies balance (or come within an acceptable range). The energy slope (Se)
can be computed by applying the Manning equation, in either SI units:
𝑛2 𝑣 2
𝑆𝑒 = 4/3
𝑅ℎ
Or BG:
𝑛2 𝑣 2
𝑆𝑒 = 4/3
2.22𝑅ℎ
Where, Se is the average of the energy (EGL) slopes at the upstream and downstream
sections. A tabulated computation procedure is recommended as illustrated in the
example problems to follow.
The astute reader may ask why △ 𝐿 is not solved for in the energy equation. By assigning
a depth to the next section instead of the depth, the equation could be used to determine
the distance between sections and avoid the iterative process altogether. This is a
legitimate solution procedure called the direct step but it only works for prismatic
channels (channels of uniform slope and cross section). When water surface profiles are
sought in natural stream channels that are non prismatic, cross sections on these streams
are field surveyed or obtained from geographic information system maps at
predetermined locations, which establishes the distance between sections. Then the
standard step method is used to assess depth of flow at these sections. Fortunately, water
surface profiles are generally solved with the aid of computer software, which takes the
drudgery out of the iterative process
Direct Step Method
In the direct step method, the gradually varied flow equations are rearranged to
determine the distance (△ 𝐿) explicitly between two selected flow depths. This method is
applicable to prismatic channels only because the same cross-sectional geometric
relationships are used for an the sections along the channel.
265
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Replacing sections 1 and 2 with U and D, respectively, and noting that 𝑆𝑜 = (𝑍𝑈 − 𝑍𝐷 )/△
△𝑍 𝑣22 𝑣12
𝐿= ), energy equation ( + 𝑦2 + 𝑧2 = + 𝑦1 + 𝑧1 + 𝑆𝑒 △ 𝐿 )is rearranged as
△𝐿 2𝑔 2𝑔
𝑣𝑑2 𝑣𝑈2
(𝑦𝐷 + 2𝑔) − (𝑦𝑈 + 2𝑔 𝐸𝐷 − 𝐸𝑈
△𝐿 = =
𝑆𝑜 − 𝑆𝑒 𝑆𝑜 − 𝑆𝑒
𝑣2
Where 𝐸 = 𝑦 + is specific energy. In this equation, U and D represent upstream and
2𝑔
downstream sections, respectively. For subcritical flow, the computations begin at the
downstream end and progress upstream. In this case yD and ED would be known. An
appropriate value for yu is selected and the associated Eu is calculated. Then △ 𝐿 is
determined by using Equation last equation.
For supercritical flow, the computations begin at the upstream end and progress
downstream. In this case, yu and Eu would be known. An appropriate value for yD is
selected, and the associated ED is calculated. Then △ 𝐿 is determined by using the last
equation.
Example
A grouted-riprap, trapezoidal channel (n = 0.025) with a bottom width of 4 meters and
side slopes of m = 1 carries a discharge 12.5 m3/sec on a 0.001 slope. Compute the
backwater curve (upstream water surface profile) created a low dam that backs water up
to a depth of 2m immediately behind the dam. Specifically, water depths are required at
critical diversion points that are located at distances of 18 8 m, 423 m, 748 m, and 1,675
m upstream of the dam.
Solution:
Normal depth for this channel can be calculated by using Manning equation (iterative
solution), or the figure for normal depth solution.
𝑛𝑄 (0.025)(12.5)
1/2
= 1/2 8/3
= 0.245
𝑘𝑀 𝑆 𝑏8/3 (1.00)(0.001) (4)
𝑜
From the following figure with m=1,
we obtain yn/b=0.415 therefore, yn=(4m)(0.415)=1.66m.
𝑄2 𝐴3
Critical depth for this channel can be calculated by using the equation = (iterative
𝑔 𝑇
3/2
𝑄𝑚 (12.5)(1)3/2
solution) or by the figure. Using the figure: = 1 = 0.125, from the figure,
𝑔1/2 𝑏5/2 (9.81)2 (4)5/2
𝑚𝑦𝑐 (4)(0.23)
we obtain = 0.23; therefore, 𝑦𝑐 = = 0.92𝑚.
𝑏 1.0
266
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Figure:Critical depth solution procedure: (a) trapezoidal channels, and (b) circular channels.
We will first use the standard step method. Water surface profile computations require
the use of the Manning equation, which contains the variables Rh, and v. Recall from
earlier discussions that Rh = A/P, where A is the flow area and P is the wetted perimeter,
and v = Q/A. The computation procedure displayed in the following table is used to
determine the water surface profile. The depth just upstream from the dam is the control
section, designated as section 1. Energy balance computations begin here and progress
upstream (backwater) because the flow is subcritical (Yc < Yn)· The finite difference
process is iterative; the depth of flow is assumed at section 2 until the energy at the first
𝑣22 𝑣12
two sections match using the equation ( + 𝑦2 + 𝑧2 = + 𝑦1 + 𝑧1 + 𝑆𝑒 △ 𝐿 ). Once
2𝑔 2𝑔
the water depth at section 2 is determined, the depth of flow at section 3 is assumed until
the energies at sections 2 and 3 balance. This stepwise procedure continues upstream
until the entire water surface profile is developed. Because the starting depth of 2.00 m
is greater than the normal depth and normal depth exceeds critical depth, the profile has
an M-1 classification (Figure). The flow depth will approach normal depth asymptotically
as the computations progress upstream, as depicted in the figure.
267
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Once the depth becomes normal, or relatively close, the computation procedure is ended.
The first few standard step computations are displayed in the following table (a);
completion of the problem is left to the student.
Because the channel considered in this example is prismatic, we can also use the direct
step method to calculate the water surface profile. Table (b) is used to determine the
𝑣2 2
𝑣
(𝑦𝐷 + 𝑑 )−(𝑦𝑈 + 𝑈
2𝑔 2𝑔 𝐸𝐷 −𝐸𝑈
profile by setting up and solving Equation △ 𝐿 = = . The
𝑆𝑜 −𝑆𝑒 𝑆𝑜 −𝑆𝑒
calculations in the table are self-explanatory. Like the standard step method, the
computations begin at the downstream end and progress upstream. For the first channel
reach considered, y0 = 2. 00 m is known, and Yu = 1. 91 m is a depth we select based on
the water surface profile (M-1; depths go down) and to compare with the standard step
method solution. Then we calculate the distance between the sections with these two
depths. For the next reach, 1. 91 m becomes the downstream depth, and we select Yu =
1. 82 m. The results are slightly different from those of the standard step method. The
discrepancies result from the iterative nature of the standard step method in which the
results depend on the tolerance limit selected.
268
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
(a) Water surface profile (Backwater) computations using the standard step method
269
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
(b) Water surface profile (backwater) computations using the direct step method.
270
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Example
A rough-concrete trapezoidal channel (n = 0.022) with a 3.5-ft bottom width, side slope
m = 2, and bed slope of 0.012 discharges 185 cfs of fresh water from a reservoir.
Determine the water surface profile in the discharge channel to within 2 percent of
normal depth.
Solution:
Normal depth and critical depth are calculated before solving water surface profiles in
order to determine the gradually varied flow classification.
Normal depth may be determined using the Manning equation in conjunction with the
figure:
𝑛𝑄 (0.022)(18.5)
1/2
= 1/2 8/3
= 0.883
𝑘𝑀 𝑆 𝑏 8/3 (1.49)(0.012) (3.5)
𝑜
From the figure , we obtain yn/b=0.685 therefore, yn=(3.5ft)(0.685)=2.40 ft.
Alternatively, normal depth can be obtained with the Manning equation (using successive
substitution) or by using appropriate computer software.
Critical depth can be computed using the equation and the table:
𝑄2 𝑇 𝑄2 (𝑏 + 2𝑚𝑦𝑐 ) (185)2 [3.5 + 2(2)𝑦𝑐 ]
= = =1
𝑔𝐴3 𝑔[(𝑏 + 𝑚𝑦𝑐 )𝑦𝑐 ]3 32.2[(3.5 + 2𝑦𝑐 )𝑦𝑐 ]3
By successive substitution, we obtain yc=2.76 ft.
Because critical depth exceeds normal depth, the channel is steep with an S-2
classification. Water from the reservoir will enter the channel and pass through critical
depth as depicted in the figure.
Because the control section is at the entrance to the channel and flow is supercritical (Yn
< Yc), computations will proceed in the downstream direction (front-water) starting from
critical depth at the entrance section and approaching normal depth asymptotically. In an
S-2 profile, the water surface elevation changes quickly at first and then approaches
normal depth more gradually. Therefore, use five cross sections, including the control
section, with separation distances (△ 𝐿) of 2, 5, 10, and 40 feet, respectively, moving
downstream. The first few standard step computations are displayed in the following
table (a); completion of the problem is left to the student in Problem 6.8.7.
271
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
We can also use the direct step method to calculate the water surface profile in this
problem because the channel is prismatic. The calculations are summarized in table (b).
Again, the calculations begin from the upstream end and progress downstream. For the
first channel reach considered, Yu = 2.76 ft is known, and YD = 2.66 ft is a depth we select
based on the type of the water surface profile (S-2; depths go down) and to compare with
the standard step method solution. We then calculate the distance between the sections
with these two depths. For the next reach, 2.66 ft is the upstream depth, and we select
YD = 2.58 ft The discrepancies between the standard step and the direct step method
result from the iterative nature of the standard step method in which the results depend
on the tolerance limit selected.
272
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
(a) Water surface profile (front-water) computations using the standard step method
273
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
(b) Water surface profile (front-water) computations using direct step method.
274
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
Another Approach: The gradually varied flow profile can be calculated using numerical
solution of its equation. Computation of the gradually varied flow surface profiles should
proceed upstream from the control section in subcritical flow and downstream from the
control section in supercritical flow.
1. The direct step method:
2. This method is applicable to prismatic channels. Depths of flow are specified and the
distances between successive depths are calculated.
Consider the elements of flow in the figure.
275
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
The depth at the control section is 4.0 m.
A=by+my2 = 3.0 (4.0)+ 1 (4.0)2= 28.0 m2.
𝑅 = 𝐴/𝑃 = 28.0 /(3.0 + 2√4.02 + (1)2 4.02 = 1.956 𝑚
𝑄2 192
Specific energy (𝑦 + ) = 4.0 + 19.62×282 = 4.024 𝑚
2𝑔𝐴2
𝑄 2 𝑛2 𝐸2 −𝐸1
𝑆𝑓 = = 5.44 × 10−5 , ∆𝑥 =
𝐴2 𝑅 4/3 𝑆𝑜 −𝑆𝑓
276
CENG 321 Hydraulics – Course Notes by Dr. Fuad Musleh
4.5
3.5
2.5
Y (m)
1.5
0.5
0
0 200 400 600 800 1000 1200 1400 1600 1800 2000
X (m)
277