CE 336 Lab 5 Report Friction in Pipes
CE 336 Lab 5 Report Friction in Pipes
Fluid Mechanics
CE 336
1
Table of Contents
Abstract………………………………………………………………………….5
Group Tasks
2
Purpose
When a fluid flows in a closed conduit pipe, its energy loss can be
demonstrated in this experiment. Frictional head losses of the fluid flow in a
pipe may be calculated using experimental and theoretical coefficients
used in the three empirical formulas listed in subsequent theory section.
3
Abstract
4
The frictional losses resulting from the fluid’s internal shear forces, are
added to the losses due to the surface imperfections of the material that
form the conduit. Real surfaces further restrict the movement of the
dynamic fluid.
Magnified surface of an
aluminum pipe
Real pipe
surfaces are
not perfectly
smooth
Theory
Fundamentally, the system concurs with the first law of thermodynamics for
opened systems. Hence, the conservation of energy principle applies:
Δ E=E¿ −( E¿¿ out + losses)=0 ¿
Consequently,
E¿ =E out +losses
and
P 1 V 12 P 2 V 22
z 1+ + =z 2+ + + headlosses Equation 1
γ 2g γ 2g
5
Z1 CONTROL Z2
P1 VOLUME P2
V1 V2
slow
fast
stream
stream
Figure 2 provides the basis for the energy analysis. The continuity
equation correlates the fluid speed and the pipe geometry to the flow.
Q1 = V1 * Area1
Q2 = V2 * Area2 Equation 2
Area1= Area2
Q1 = Q2
Assuming that Z1 = Z2 and that once a steady flow has been established,
the average speed of the water within the pipe is constant in the
longitudinal direction and so V1 average = V2 average.
6
The calculation for the actual head loss through the pipe can now be
determined from the pressure head differences at the entrance and exit of
the control volume:
length pipe∗frictioncoefficient∗ρ V 2
P1−P2=
2∗Diameter pipe
Knowing the fluid speed and viscosity, the flow’s Reynolds number can be
determined using:
The actual head loss through the pipe can be measured experimentally,
using a manometer. Then the coefficient of friction can be determined using
equation 4. Alternatively, the coefficient can also be determined using the
Moody chart, if the conduit’s roughness ratio is known.
The roughness
can be
estimated
knowing the
pipe material
7
The friction
coefficient will be
estimated with
the knowledge of
the pipe’s
roughness and
the flow’s
Reynolds Reynolds number will be determined
number
with equation 5
Equipment
The main components were the pipes in which the water was flowing for
this experiment. Below, figure shows the flow pipe system.
low pressure
The ¾ The ¾ head The
‘old’ pipe ‘new’ pipe
pipe
system
The ½
High pressure pipe
head
8
The pressure drop through the system was being monitored. Below, figures
1 and 2 show the capillary system which can divert a fraction of the flow to
the manometer.
The high
pressure
side
capillary
tube
valves
The low
pressure
side
capillary
9
tube
valves
The pipe
discharge
valves
The capillary
tubes route
the flow to
the
manometer
The
manometer
determines
the pressure
drop
The pipe
discharge
hose to the
10
weighing
station
The scale
weighs the
water
discharge
and
determines
the actual
flow
The water pressure on the high side of the pipe system was kept constant
during each run. During a run, an observer was posted at the input valve to
continuously control the system pressure according to the pressure gauge.
Below, the pressure control center is shown.
The gauge
shows the
high side
pressure
11
Procedure Part I
1. Open the main water valve. Then open all the capillary tubes valves to
verify an adequate flow to the manometer. Verify the flow. Then close all
the capillary tubes.
2. Place the discharge collection vessel on the scale and reset zero for
scale (tare weight).
3. Open the incoming and outgoing valves of the ¾” ‘old’ pipes were
opened so that the flow could move within its volume.
5. Set the inbound system pressure to 10 PSI and then allow the flow to
discharge into the receptacle for a few minutes, until flow is steady.
6. Simultaneously divert the flow to the weighing vessel and start timing
using a chronometer to measure and record how much time it takes to
12
collect ½ cubic foot of water. Then calculate the run’s actual volumetric flow
rate. Then discard the weighed water.
Procedure Part II
13
Raw Data
Part I:
14
31.2 43.12 0.0116 23,6 - 7.80 1.2 16.0125
10 5.63 81 7.45 7
31.2 48.09 0.0104 21,2 - 7.60 1.2 15.6975
8 5.05 34 7.35 5
31.2 53.86 0.0093 18,9 - 6.15 1.0 12.5475
6 4.51 59 5.80 0
½ 31.2 66.19 0.0076 15,4 - 3.65 0.6 7.5600
new 4 3.67 27 3.55 0
Part II:
Pipe Pressure Act. Hydraulic T.E.L
15
Head Darcy
Moody Radius Hazen Manning
(psi) Loss h2 W. Slope S
f_theo Rh(ft) W. C n
(ft) f_exp
10 15.70 0.1718 0.1189 0.01633 0.805 50.61 0.0193
¾ inch 8 14.81 0.1878 0.1300 0.01633 0.759 48.51 0.0202
old 6 10.82 0.1864 0.1290 0.01633 0.555 49.31 0.0201
4 6.51 0.1895 0.1312 0.01633 0.334 49.91 0.0203
10 19.16 0.1736 0.1202 0.01761 0.983 49.47 0.0197
¾ inch 8 16.85 0.1942 0.1344 0.01761 0.864 47.02 0.0208
new 6 12.08 0.1892 0.1310 0.01761 0.619 48.27 0.0205
4 7.61 0.1866 0.1292 0.01761 0.390 49.51 0.0204
10 16.01 0.0854 0.0591 0.01281 0.821 73.97 0.0131
8 15.70 0.1041 0.0721 0.01281 0.805 67.04 0.0145
½ new
6 12.55 0.1044 0.0723 0.01281 0.643 67.55 0.0145
4 7.56 0.0950 0.0658 0.01281 0.388 72.26 0.0138
16
Table 5: Data for Part II of the experiment
Sample Calculation
For Part I:
Convert the diameter to feet to maintain consistency of the measurement
units:
0.784 inches
D= =0.0653 foot
12
2
' π 2
Pip e s Area= ( 0.0653 ft ) =0.00335 ft
4
Calculate the actual flow rate of the water, knowing the time interval it took
17
0.5 ft 3 ft 3
Q= =0.0229
21.8 sec sec
Calculate the flow’s velocity, knowing the flow rate quantity and the cross
section of the pipe in which this flow is moving. Now, we can determine the
ft 3
0.0149
Q sec ft
Velocity of flow = = 2
=4.44
Area 0.00335 ft sec
ft
'
Reynol d s Number =
Velocity∗Diameter
=
( 4.44
sec )
( 0.0653 ft )
=23,827
kinematic viscosity ft 2
1.217E-5
s
The mercury’s delta height was then converted to feet so that the head loss
Δ h ft .=¿ ¿
18
Calculate the head loss, knowing the density difference between mercury
and water.
p1 p2 w
hL m 12.6 * ∆h
w w
Absolute roughness for a galvanized steel equals 0.005 ft. Knowing the
K 0.005 ft
Relative Roughness= = =7.66E-3 ft
D 0.0653 ft
Using this quantity, the Moody diagram was used to calculate a theoretical
value for the pipe’s frictional coefficient. According to the diagram, this
theoretical value is
f theoretical=0.045
ft
f experimental=
(
( 15.70 ft ) ( 0.0653 ft ) ( 2 ) 32.2
sec2 ) =0.1717
2
ft
(
( 19.5 ft ) 4.44
sec )
The experimental coefficient of friction for the pipe can be calculated and
equation,
19
To determine the dimensionless pipe characteristics, the hydraulic radius of
the pipe and the energy line slope were first calculated.
Area 0.00335 ft 2
Hydraulic Radius ( Rh ) = = =0.0163 ft
Perimeter π ( 0.0653 ft )
H Loss 15.70 ft
Slope(T . E . L .)= = =0.805
Pipe length 19.5 ft
ft
4.44 =1.318C 0 ( .0163 ft ).63 ( 0.805 ).54∧so C=50.61
sec
ft 1.486
4.44 = ( 0.0163 ft )2 /3 ( 0.805 )1/ 2∧so n=0.0193
sec n
Below, table 6 shows the pipe dimensions and their respective estimated
20
0.0705 0.0039 0.00709
4 0.846 0 0.00050
0.0512 0.0020 0.00977
10 0.614 6 0.00050
0.0512 0.0020 0.00977
8 0.614 6 0.00050
½ new
0.0512 0.0020 0.00977
6 0.614 6 0.00050
0.0512 0.0020 0.00977
4 0.614 6 0.00050
Next, table 7 and table 8 show the numerical analysis of the first part of the
experiment. They are the values which are required in the lab manual’s
21
31.2 54.8 0.0091 14,6 - 3.15 0.5 6.5100
4 2 2.72 09 3.05 2
31.2 25.2 0.0198 29,3 - 9.30 1.5 19.1625
10 9 5.07 67 8.95 2
31.2 28.5 0.0175 26,0 - 8.20 1.3 16.8525
8 2 4.50 41 7.85 4
31.2 33.2 0.0150 22,3 - 5.90 0.9 12.0750
¾ 6 6 3.85 30 5.60 6
inch 31.2 41.6 0.0120 17,8 - 3.75 0.6 7.6125
new 4 3.08 53 3.50 0
31.2 43.1 0.0116 23,6 - 7.80 1.2 16.0125
10 2 5.63 81 7.45 7
31.2 48.0 0.0104 21,2 - 7.60 1.2 15.6975
8 9 5.05 34 7.35 5
31.2 53.8 0.0093 18,9 - 6.15 1.0 12.5475
6 6 4.51 59 5.80 0
½ 31.2 66.1 0.0076 15,4 - 3.65 0.6 7.5600
new 4 9 3.67 27 3.55 0
Act.
Hydraulic
Pressure Head Darcy T.E.L
Pipe Moody Radius Hazen Manning
(psi) Loss W. Slope S
f_theo Rh(ft) W. C n
h2 (ft) f_exp
10 15.70 0.1718 0.1189 0.01633 0.805 50.61 0.0193
¾ 14.81 0.1878 0.1300 48.51 0.0202
8 0.01633 0.759
inch
6 10.82 0.1864 0.1290 0.01633 0.555 49.31 0.0201
old
4 6.51 0.1895 0.1312 0.01633 0.334 49.91 0.0203
¾ 10 19.16 0.1736 0.1202 0.01761 0.983 49.47 0.0197
22
8 16.85 0.1942 0.1344 0.01761 0.864 47.02 0.0208
inch 12.08
6 0.1892 0.1310 0.01761 0.619 48.27 0.0205
new
4 7.61 0.1866 0.1292 0.01761 0.390 49.51 0.0204
10 16.01 0.0854 0.0591 0.01281 0.821 73.97 0.0131
½ 8 15.70 0.1041 0.0721 0.01281 0.805 67.04 0.0145
new 6 12.55 0.1044 0.0723 0.01281 0.643 67.55 0.0145
4 7.56 0.0950 0.0658 0.01281 0.388 72.26 0.0138
Part II:
Again, the experimental volumetric flow rate of the water is calculated with
the knowledge of the time interval it took to collect ½ cubic foot of water in
0.5 ft 3 ft 3
Q= =0.0219
22.79 sec sec
23
Next, it was necessary to obtain the speed in the parallel pipes so that a
calculated head loss could be determined and then compared to the actual
head loss which was measured using the manometers. The actual head
The speed was calculated based on the set of two equations with 2
unknowns that were given with the Darcy-Weysbach relation and the fact
that the arithmetic sum of the 2 flow rates would have to equal the total flow
f 1 L1 V 21 f 1 L1 V 21
=
2 g D1 2 g D1
D 1 f 2 V 22
¿ so V 1=
√
D2 f 1
D1 f 2
¿ finally , Qtotal = A1 (√ ) √D2 f 1
V 22 + A2 V 2
Q total A2 V 2
Therefore , =V 2 +
D1 f 2 D1 f 2
A1 (√ )
D2 f 1
A1 (√ )
D2 f 1
ft
Ultimately , V 2=3.02
sec
24
ft
V 1=3.32
sec
ft 2 (
Head losscalculated =
( 3.32
sec ) 15.70 ft ) (0.0708)
=1.69 ft
ft
(
( 0.0653 ft )( 2 ) 32.2
)
sec
2
This number was compared to that of the actual head loss and a relative
error was calculated in percentage.
25
Manometer Measure The
Reading d Head f1 Head
f2 %Error
Combination Pressur Left Right Delta Loss hL Loss
Pipes e (psi) (in) (in) in (ft) (ft) hL
- 1.75 0.30
10 1.65 4.08 0.086 0.076 1.69 58.7
- 1.60 0.24
8 1.70 3.29 0.087 0.075 1.37 58.4
D1 0.75" - 1.35 0.20
and D2: 0.5" 6 1.80 2.72 0.083 0.077 1.19 56.2
Results
Below provides the friction coefficients that were derived from the
experimental data. The theoretical coefficients were obtained on a Moody
diagram with the known theoretical pipe roughness and the experimental
Reynolds numbers. The water’s energy line was calculated from the
experimental data. In an effort to correlate the experimental data to popular
26
methods of quantifying the friction that exists within closed pipes, the
Hanzen-Williams and the Manning coefficients were used for calculations.
Part I:
27
31.2 54.82 0.009 14,6 - 3.15 0.52 6.5100
4 1 2.72 09 3.05
31.2 25.29 0.019 29,3 - 9.30 1.52 19.1625
10 8 5.07 67 8.95
31.2 28.52 0.017 26,0 - 8.20 1.34 16.8525
8 5 4.50 41 7.85
31.2 33.26 0.015 22,3 - 5.90 0.96 12.0750
6 0 3.85 30 5.60
¾
inch 31.2 41.6 0.012 17,8 - 3.75 0.60 7.6125
new 4 0 3.08 53 3.50
31.2 43.12 0.011 23,6 - 7.80 1.27 16.0125
10 6 5.63 81 7.45
31.2 48.09 0.010 21,2 - 7.60 1.25 15.6975
8 4 5.05 34 7.35
31.2 53.86 0.009 18,9 - 6.15 1.00 12.5475
6 3 4.51 59 5.80
½ 31.2 66.19 0.007 15,4 - 3.65 0.60 7.5600
new 4 6 3.67 27 3.55
28
6 0.846 0.0705 0.00390 0.00050 0.00709
4 0.846 0.0705 0.00390 0.00050 0.00709
10 0.614 0.0512 0.00206 0.00050 0.00977
8 0.614 0.0512 0.00206 0.00050 0.00977
½ new
6 0.614 0.0512 0.00206 0.00050 0.00977
4 0.614 0.0512 0.00206 0.00050 0.00977
Act.
Hydraulic
Pressure Head Darcy T.E.L
Pipe Moody Radius Hazen Manning
(psi) Loss h2 W. Slope S
f_theo Rh(ft) W. C n
(ft) f_exp
10 15.70 0.1718 0.1189 0.01633 0.805 50.61 0.0193
¾ inch 8 14.81 0.1878 0.1300 0.01633 0.759 48.51 0.0202
old 6 10.82 0.1864 0.1290 0.01633 0.555 49.31 0.0201
4 6.51 0.1895 0.1312 0.01633 0.334 49.91 0.0203
10 19.16 0.1736 0.1202 0.01761 0.983 49.47 0.0197
¾ inch 8 16.85 0.1942 0.1344 0.01761 0.864 47.02 0.0208
new 6 12.08 0.1892 0.1310 0.01761 0.619 48.27 0.0205
4 7.61 0.1866 0.1292 0.01761 0.390 49.51 0.0204
29
10 16.01 0.0854 0.0591 0.01281 0.821 73.97 0.0131
8 15.70 0.1041 0.0721 0.01281 0.805 67.04 0.0145
½ new
6 12.55 0.1044 0.0723 0.01281 0.643 67.55 0.0145
4 7.56 0.0950 0.0658 0.01281 0.388 72.26 0.0138
30
Table 15: Part II results
Graphs
31
Darcy Weisbach Friction vs. Pressure (psi)
0.2
0.18
0.16
HalfNewPipe
Friction
0.14
QuarterNewPipe
QuarteOldPipe
0.12
0.1
0.08
4 5 6 7 8 9 10
Pressure (psi)
32
3/4 New Friction Coefficient vs. Reynolds Number
0.5
0.4
0.3
0.2
0.1
Friction Coefficient
-0.1
-0.2
-0.3
-0.4
QuarterNewPipe
-0.5
1.6 1.8 2 2.2 2.4 2.6 2.8 3
Reynolds Number 4
x 10
33
3/4 Old Friction Coefficient vs. Reynolds Number
0.5
0.4
0.3
0.2
0.1
Friction Coefficient
-0.1
-0.2
-0.3
-0.4
QuarterOldPipe
-0.5
1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4
Reynolds Number 4
x 10
34
1/2 New Friction Coefficient vs. Reynolds Number
0.5
0.4
0.3
0.2
0.1
Friction Coefficient
-0.1
-0.2
-0.3
-0.4
HalfNewPipe
-0.5
1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4
Reynolds Number 4
x 10
35
Discussion
The study, given its high level of experimental difficulty, gave results that
which were obtained from a Moody diagram. The ¾-inch old pipe had the
best correlation with a divergence of less than 6% while the 2 other pipes’
The parallel flow analysis gave results that were highly conclusive: The
was about 13%. The closest correlation data was obtained from the values
that were obtained while testing with a 10 PSI pipe pressure – the
divergence was less than 7%. The 4 PSI test was nearly as correlative with
36
Conclusion
shown that pipe flow velocity was a function of pipe internal pressure, as
relation was almost linear. Water speed decreased about 10% for every 2
PSI drop in pressure. However, there was a hidden energetic cost to the
elevated pressures. In all cases, the energy line slope was more
pressure because the pipe’s ability to impede flow is greater while the
have come from using a stopwatch to time flow quantity with the
mechanical scale. A digital unit would have given a better accuracy and
repeatability.
37
Bibliography
38