Lab Manual
Lab Manual
Lab Manual
UNIVERSITY OF OTTAWA
CVG2116
INTRODUCTORY FLUID MECHANICS LABORATORY
LABORATORY MANUAL
1.
2.
3.
4.
5.
6.
2010 EDITION
Winter 2010
TABLE OF CONTENTS
A.
Objective ...............................................................................................2
Laboratory Procedure ..........................................................................2
Submission of Reports ........................................................................3
Format ...................................................................................................3
EXPERIMENTAL SETUPS........................................................................8
1.0
2.0
3.0
4.0
5.0
6.0
B.
C.
ii
A.
Winter 2010
1.0
Winter 2010
OBJECTIVE
It should be noted that although results do not always seem to fit the prescribed
theory, they are not necessarily incorrect. It is your task and obligation during
working on this course and certainly in your future engineering career to
determine "what went wrong". Remember that there is always an explanation, be
it your fault or not. Because this is such an important aspect of engineering, any
"cooking" of the experimental results and calculations is unacceptable.
You should also indicate in your repots the practical lessons and applications of
the theory or of the results that you have obtained. This tool is in preparation for
the "real world".
2.0
LABORATORY PROCEDURE
Winter 2010
No horseplay will be tolerated in the laboratory. Students are not to tamper with
experimental apparatus beyond adjustments specified in this handout or by the
Instructor. It is not necessary to wear a lab coat for these exercises.
3.0
SUBMISSION OF REPORTS
All data obtained will be shared among the group members, whereas each group
will submit one laboratory report. However, it is strongly advised that each
student participate in every laboratory report! Duties can be divided based on
data analysis, presentation of results, theoretical and explanations and
experimental procedures, discussions and/or conclusion and recommendations.
It is also suggested that each member of the group read the report before it is
handed in. That way, all group members will endorse the report collectively and
give it their seal of approval.
The text of the report must be produced using a word processor on a computer.
Though recommended, it is not always necessary to use the word processor to
enter equations or graphical material in the report; these may be hand-written.
Unless otherwise stated, the laboratory report should be handed in two weeks
after the experiment has been performed unless otherwise stated. A good
practice for report writing is to always assume that the reader is not well informed
on the report topic! Therefore it is your responsibility be as clear and detailed as
possible, especially with respect to how you layout the wording and phrasing of
your sentences.
4.0
FORMAT
4.1
Winter 2010
Title Page
Abstract
It is advised that this section is written last, since it is a brief summary of your
report. It should include the objective of the experiment, the method by which you
have conducted the experiment, the key findings of your experiment and your
main conclusion.
4.3
Objective
In your own words, write one brief paragraph (3 to 4 lines) outlining the
purpose(s) of the laboratory experiment. This section may be separate or
included towards the end of the introduction after presenting the reader with the
purpose behind carrying out such an experiment.
4.4
Introduction
This section is intended to present the reader the purpose behind your
experiment or investigation. It begins by familiarizing the reader with historical
developments in this particular field and by mentioning the purpose or
importance of conducting such an investigation. The authors must present the
scope of the experiment in the context of the "real world" and give examples of
practical applications.
4.5
Winter 2010
Theory
This section should briefly provide the reader with the background in order to
understand the framework of the experiment and the principles by which the
theory has been derived and simplified.. Derivations should include equations
numbered in sequence, where all new terms must be identified clearly with the
appropriate symbols and units (if applicable).
4.6
List of Equipment/Apparatus/Materials
This section should be in point form. At least one sketch or photograph of the
main features of the laboratory apparatus is required. Do not forget to label and
number your figure captions accordingly!
4.7
Experimental Procedure/Methodology
This section should be in point form. You may also make a reference to the lab
handout or manual (e.g., As outlined in the laboratory manual/handout). Note
and describe any steps taken that were different from the procedure in the
handout.
4.8
In tables (numbered and captioned), show data collected during the experiment.
Columns/row headings should be indicated using the appropriate titles, symbols
and units. This section may be included in an appendix, especially if it contains a
large set of data.
4.9
Calculations
Winter 2010
specified. Make tables for repeated calculations in the "Results" section. Each
experiment in this manual contains a section describing the procedure for
performing the calculations. Do not forget units and always use appropriate and
consistent significant figures!
4.10
Results
Present all of your results in neat and clear tabular form where graphs should be
used for plots. A plot of the results will usually help you identify the pattern/trend
or the mathematical-physical relationship between the different variables. Give
the tables and graphs concise headings with table/figure numbers. Do not write
numerical results with insignificant digits.
Example:
Measured flume width
= 12.7 cm
= 6.5 cm
Area
= 12.7cm x 6.5 cm
= 82.55 cm2 (incorrect)
= 82.6 cm2 (correct)
4.11
Discussion
Refer to the heading "Report" in the individual laboratory handout. You may also
include one paragraph summarizing the theory (no equations are necessary) in
this section. After a complete understanding of the theory behind the experiment,
you should be able to analyze the numerical and graphical results. You should
refer to the results in this section (using the corresponding figure/table number)
and discuss any trend or behavior that you have observed. Remember that the
results may not always fit the theory. You should be able to surmise non-trivial
reasons for deviations. You should indicate the accuracy of the results, as well as
both major and minor causes of error or inaccuracy. You must also state possible
Winter 2010
Discuss whether the objective(s) have been achieved and state what indication
prompted you to conclude that. Summarize your results and percentage errors,
and comment on the main source of error (if any). If you have an original idea or
concept, you have the full right to claim it as such. However, if you use an idea or
results (not raw data!) from someone else, for example after a discussion with a
classmate, or from a published article, you must indicate that this is the case.
You will not lose marks if you do this judiciously. Simply remember that you are
expected to be able to research and think matters out; copying material from
another group is not acceptable and will not be tolerated! This section of the
report is not to be longer than one page. Keep it brief and concise.
4.13
References
Any literature sources that you have consulted in writing this report (e.g.
laboratory manual, course textbooks, scientific papers, reports, etc) must be
included in this section. Please note that Wikipedia and similar web based
encyclopedias are not considered official references! Please refer to the Guide
for writing Laboratory Reports of the Civil Engineering Department with respect to
article or book citations.
4.14
Appendices
This section (if applicable) should comprise of additional information not included
in the report such as: raw data, computer codes, sample and other calculations,
other tables and figures, other hardcopies, etc.
B.
Winter 2010
EXPERIMENTAL SETUPS
Gate
1.0
Winter 2010
Introduction
The flow rate will depend upon the velocity of the flow. If, for example, there is a
flow of water from the end of a pipe with the same velocity, v0, in all parts of the
flow, an increment of volume, V, will flow out of the pipe in a time increment, t,
and is equal to:
V = v 0 At
(1.1)
The velocity, v0, is the velocity normal to the cross-sectional area, A, of the pipe.
The flow rate, Q, can be expressed as:
Q=
dV
V
= limt 0
= v0 A
dt
t
(1.2)
Suppose that the normal velocity varies from point to point throughout the flow.
It would be reasonable to find the flow rate by dividing the area into increments,
multiplying each area by the normal velocity through it, and summing these
partial flow rates. The general expression for the flow rate in this case is the
integral of the normal velocity times the elemental area, dA, over the total crosssectional area, A.
Winter 2010
Q = v n dA
(1.3)
v=
Q
A
(1.4)
The Pitot tube, or stagnation tube, which was invented by Henri de Pitot (16951771), can be used to determine the velocity of flow at different points throughout
the flow by measuring the total or stagnation pressure head, pstagnation/, and the
static pressure head, pstatic/, at a point. The difference between the two pressure
heads is the velocity head, vn2/2g, as shown in the expression below,
v n2 pstagnation pstatic
=
2g
(1.5)
The most common type of Pitot tube was developed by Ludwig Prandtl (18751973). For this type, the tubes for the total (or stagnation) pressure and the static
pressure measurement are combined into one piece, whereas the openings for
the static pressure measurement are located at the proper point along the body
of the probe to give the same static pressure reading, as if no tube was present.
Proper use of the instrument requires that it be correctly aligned to point directly
(or parallel) into the oncoming flow.
10
1.2
Winter 2010
Objective
Experimental Procedure
1.
2.
Measure the depth of the flow using the point gauge provided.
3.
Lower the Pitot tube to the bottom of the channel and record the vertical
position.
4.
Raise the Pitot tube such that it is right below the surface (by about
10mm) and record its vertical position.
5.
Subtract the two positions obtained in step 3 and 4, and divide the
difference by 4 to obtain the vertical increment for each flow depth to be
used.
6.
Place the Pitot tube on the channel bottom, as close to the wall as
possible (left or right), starting at the lowest position. Keep the tube
aligned with the flow.
7.
Record the manometer readings for both the static and total pressure
head.
8.
Slide the Pitot tube towards the opposite wall by adding the horizontal
increment specified on the data sheet to the horizontal position, while
recording both manometer readings.
9.
Once the manometer readings have been recorded four times (i.e., wallto-wall), repeat readings of the Pitot tube at 1/4 depth, 1/2 depth, 3/4
depth, and near the surface by adding the vertical increment to the
vertical location. Please refer to the figure in the data sheet for more
clarification.
10. Measure the flow rate using the tank on the scale at the end of the
channel and a stop watch. Use at least 3 discharge measurements.
11
1.4
Winter 2010
Calculations
1. Calculate the velocity, v, using Eq. (1.5) at each point where the static and
stagnation pressure heads were measured.
2. Plot each point on a cross section of the flow and label it with the velocity
found there.
3. Draw lines of equal velocity (isovelocity or isovels) as you would draw
contour lines on a topographic map. It is likely that interpolation between
the data points will be required to plot the isovelocity lines. Note: There
are several programs that can do this task (e.g. Surfer, Golden Software).
4. Calculate the mean velocity between each adjacent isovelocity lines and
measure the area between them. Determine the discharge through
adjacent isovelocity lines by multiplying the mean velocity by the area
between them. Then determine the total discharge in the cross-sectional
area by summing the discharge through adjacent isovelocity lines as
indicated in Eq. (1.3).
5. Calculate the flow's mean velocity by dividing the discharge by the total
cross-sectional area.
1.5
Report
12
Winter 2010
1/3
2/3
3/4
Depth
Width W : _________
Depth D : _________
Pitot tube diameter: 4mm
N.B. Flow going through
the sheet
1/2
1/4
Left wall
Width
Stagnation (or total) and static pressure head (cm or mm) for the cross-section:
Left wall
1/3 width
2/3 width
Right wall
Position: 247mm
Position: 288mm
Position: 328mm
Position: 370mm
Bottom
Position:
depth
Position:
depth
Position:
depth
Position:
Surface
Position:
____________
13
2.0
Winter 2010
BERNOULLIS EQUATION
2.1
Introduction
The form of the equation most often used in hydraulics, expresses each
parameter in terms of "heads" as shown below,
+z+
v2
= constant
2g
where,
(2.1)
p/ is pressure head.
z is elevation head.
v2/2g is velocity head.
The first two heads are combined to give the piezometric head, h, as shown
below,
h=
+z
(2.2)
In the calculations for this experiment some dimensionless plots are required.
An advantage of dimensionless plots is that the results for more than one set of
conditions may be used for other calculation based on similarity laws. Reduced
to non-dimensional terms, the results from many different situations may be able
to be expressed in a single graphical relation. This principle of similarity is very
important in physical modeling in which results from geometrically small scale
tests are applied to a similar situation on a larger scale (also known as
prototype).
14
2.2
Winter 2010
Objective
The objective of this experiment is to apply and validate the Bernoulli equation to
flow in a closed conduit and to observe the energy loss in a Venturi metre.
2.3
Experimental Procedure
1. Establish steady flow in the apparatus when the flow rate is near full
capacity.
2. Make sure that there is no air in the system.
3. Record the water levels in the piezometer tubes (hA to hL). Note: There is
no piezometer/station at I. Also, the zero mark on the piezometers scale
is 167mm above the Venturi centerline.
4. Measure the flow rate by weighing as directed by the lab demonstrator.
5. Repeat the process at three lower flow rates.
15
2.4
Winter 2010
Calculations
Area
Spacing
(mm)
(mm2)
(mm)
25.99
530.5
23.97
451.3
20
18.21
260.4
12
15.75
194.8
14
16.01
201.3
15
17.65
244.7
15
19.29
292.2
15
20.94
344.4
15
22.50
397.6
15
24.22
460.7
15
26.29
542.8
20
Station
H A = hA +
vA2
2g
(2.3)
4. For each of the four runs, compute h/HA. Plot these values (h/HA) for all
stations against the distance along the conduit, all on one graph. Use a
different plotting symbol (legend and/or colour) for each run.
16
Winter 2010
All dimensions in mm
Direction of Flow
A
B
20
C
12
D
14
E
15
F
15
G
15
H
15
J
15
15 20
Report
17
Winter 2010
BERNOULLIS EQUATION
hE
[mm]
Q1
Time:
Mass:
Q2
Time:
Mass:
Q3
Time:
Mass:
Q4
Time:
Mass:
18
hF
[mm]
hG
[mm]
hH
[mm]
hJ
[mm]
hK
[mm]
hL
[mm]
3.0
Winter 2010
Introduction
Method 1: Integration of the fluid pressure exerted on a plane surface gives the
component of the fluid force normal to the plane.
Fn = pdA
(3.1)
Similarly, the force of fluid pressure on a curved surface can be obtained from a
known pressure distribution. The determination of the force by its individual
components may, in this case, be the best procedure; in some cases, the use of
a coordinate system adapted to the curved surface may be advantageous.
The method of integration is a basic technique for finding pressure forces which
can be used whenever the pressure distribution is known, either by
measurement, or by analysis. Also, the point of application of the force (or centre
of pressure) can be found by integrating the moment of the force about the
subject axis.
19
Winter 2010
For steady flow, the momentum flux equation states that the sum of the external
forces on the fluid within a certain volume, called a control volume, is equal to the
difference between the momentum flux out of the control volume and the
momentum flux into the control volume (i.e., the net rate of change of
momentum). In many cases where the flow rate, Q, is constant, the velocity can
be represented by an average velocity, v, so that this statement can be
expressed as,
JG
JJG JJG
2 v1
F = Q (v
where,
(3.2)
The control volume is usually chosen in a way that simplifies the problem,
keeping in mind what are the known and required quantities.
The momentum equation does not give any information on the point of
application of any of the external forces. Also, when necessary, body forces,
such as the weight of the fluid in the control volume, must be included in the
summation of forces.
3.2
Objective
20
3.3
Winter 2010
Experimental Procedure
Note: Perform the experiment with the conditions set by the lab instructor. Do not
attempt to change any of these settings!
1. Measure the head over the V-notch weir using the manometer on the side.
2. Ask the lab instructor to adjust the height of the gate opening for the first
condition.
3. Measure the gate opening, , using a metallic ruler.
4. Measure the depth of water upstream, Y1, using a metallic ruler. Ensure
that the measurement is recorded as close to the gate as possible.
5. Measure the depth of the water, Y2, downstream of the gate using the
point gauge. These measurements should be made where the depth is
uniform. Use the average of three measured values.
6. Record the manometer readings for the piezometer taps connected to the
gate. Ensure that there is no air in the lines.
7. Ask the lab instructor to set conditions for a second run using a larger gate
opening and then repeat steps 3 through 6.
8. It is convenient to use the channel bottom as a reference datum. To
establish the zero of the manometers, determine the manometer reading
for the tap in the channel bottom and, at the same time, measure the
depth of at one of the tapping points on the channel (i.e., downstream the
gate). Subtract the depth of water from the manometer reading to get the
zero reading of the manometers.
9. Measure the channel width using a metallic ruler.
10. From the demonstrator, obtain the zero reading for the V-notch weir. The
height of the tapping points on the gate with respect to the bottom of the
gate is provided on the data sheet.
3.4
Calculations
1. Find both the stagnation and hydrostatic pressure heads at each tapping
point on the gate.
21
p
=hz
stag
p
= Y1 z
hydr
where,
Winter 2010
(3.3)
(3.4)
2. Plot the pressure head vs. distance above channel bottom for the two
runs. Each run (i.e. gate opening condition) should be on different graphs,
however, a comparison between the stagnation and hydrostatic pressure
head for the same gate opening condition should be carried out on the
same graph. Assume that the stagnation pressure head is zero at the
water surface and at the bottom of the gate (Note: The hydrostatic
pressure head is zero at the water surface only). Draw smooth curves
showing the pressure variation for the stagnation pressure head plot.
3. Find the area under the pressure head curves (for both stagnation and
hydrostatic pressure calculations), apply any scaling factors, and calculate
the normal forces on the gate for the two runs (i.e. for both gate openings).
p
Fn = dA
A
(3.5)
4. Use the standard V-notch weir formula to find the flow rate, Q,
Q = 1.365H 2.5
(SI units)
(3.6)
22
Winter 2010
5. Using a control volume that includes the water between the uniform flow
sections where the depths were measured, find the force on the gate
using the momentum equation. Be sure to include the hydrostatic forces
on the water at the two ends of the control volume.
3.5
Report
1. Comment on the pressure distribution for the two conditions and compare
the plot obtained from the stagnation pressure head to the hydrostatic
pressure head.
2. Discuss why the pressure distribution shown in your plot differs from the
pressure distribution when the gate is closed. Explain how the gate
opening affects the pressure distribution on the gate.
3. Compare the values of the forces calculated by the two methods.
4. Comment upon the advantages and disadvantages of the two methods.
5. Discuss possible practical applications in calculating forces on structures.
6. List possible sources of errors involved with respect to each method when
determining the force on the gate, as well the sources of error generally
noticed while conducting this experiment.
23
Winter 2010
24
4.0
Winter 2010
IMPULSE TURBINE
4.1
Introduction
The impulse turbine operates at atmospheric pressure by the force of a highspeed jet acting on vanes mounted on a wheel. The type of impulse water turbine
used in this experiment is called the Pelton turbine and is essentially a
development of the water wheel. Pelton turbines are used where the head is high
and the flow rate is relatively low.
The jet of an impulse turbine is produced by a needle valve which controls the
flow rate. In prototype turbines, the opening of the valve is adjusted automatically
by speed control devices as the load on the turbine varies. The load on the
turbine is the electrical generator which must rotate at constant speed while
generating electrical current.
25
Winter 2010
3. Head at the needle valve, H. This head represents the energy available to
the turbine. Available power or power input is expressed as:
Pin = HQ
(4.1)
4. Power output or brake power. This is the power transmitted to the shaft by
the turbine. In a test setup it is usually measured by a Prony brake which
is a device used to create and measure the torque on a rotating shaft.
Then the power output can be expressed as,
Pout = T
Where,
(4.2)
T is the torque.
4.2
Pout
Pin
(4.3)
Objective
Experimental Procedure
Winter 2010
Calculations
T = Fd
(4.4)
Report
27
Winter 2010
IMPULSE TURBINE
1st head at needle valve: ________
Break #1
Break #2
Break #3
Break #4
Break #5
Break #3
Break #4
Break #5
Break #3
Break #4
Break #5
Brake Force
[lbf]
Discharge
[cfm]
Shaft speed
[rpm]
2nd head at needle valve: ________
Break #1
Break #2
Brake Force
[lbf]
Discharge
[cfm]
Shaft speed
[rpm]
3rd head at needle valve: ________
Break #1
Break #2
Brake Force
[lbf]
Discharge
[cfm]
Shaft speed
[rpm]
28
Winter 2010
100
90
80
70
60
50
40
Example:
Observed: 43 ft
Corrected: 36.3 ft
30
20
10
0
0
10
15
20
25
30
Pressure, P [PSI]
Observed Reading
True Reading
29
35
40
5.0
Winter 2010
Introduction
L v2
hf = f
D 2g
(5.1)
The Moody diagram (refer to the Fluid Mechanics textbook) relates f to the
Reynolds number and the relative roughness (e/D or ks/D) in the most convenient
way. Two commonly used formulas for certain types of flow are shown below,
Laminar flow:
f =
64
Re
(Hagen-Poiseuille equation)
30
(5.2)
Winter 2010
5.2
0.316
Re1 4
(Blasius equation)
(5.3)
Objective
The objective of this experiment is to measure the head loss in pipes of different
diameters and at different flow rates, as well as to compute experimental values
of the Darcy-Weisbach coefficient, f, and then compare them with theoretical or
accepted values.
5.3
Experimental Procedure
1. After checking with the lab instructor that the apparatus is ready, open the
manometer valves for the first pipe. Be sure the other manometer valves
are closed.
2. Slowly open the discharge valve of the pipe until the maximum flow
possible is established in the pipe. The maximum flow will be determined
by the maximum height of the manometer, or by the fully open valve
position. Allow a few minutes after setting the valve to ensure steady
conditions.
3. Read the water level in the left and right manometers and determine the
flow rate using the level gauge on the volumetric tank at the outlet.
4. Take 4 additional readings at lower flow rates, repeating the manometer
readings and the volumetric rate measurements.
5. Repeat steps 1 through 4 for the other pipes.
6. Record the average temperature of the water during the experiment.
Note: Pipe internal diameters are: 7.6mm, 14.5mm, 20.6mm and 26.8 mm,
volumetric tank area is 0.207 m2 and the distance between the manometer
taps, L, is 2.44 m.
31
5.4
Winter 2010
Calculations
Report
Note: The transition zone is very unstable and may extend beyond the region
indicated on the Moody diagram.
32
Winter 2010
Left manometer
reading (mL) [cm]
Right manometer
reading (mR) [cm]
Left manometer
reading (mL) [cm]
Right manometer
reading (mR) [cm]
Left manometer
reading (mL) [cm]
Right manometer
reading (mR) [cm]
Left manometer
reading (mL) [cm]
Right manometer
reading (mR) [cm]
33
6.0
Winter 2010
FORCED VORTEX
6.1
Introduction
There are two types of vortices distinguished in the dynamics of the motion and
the resulting velocity distributions. These are forced and free vortices. The forced
vortex is caused by external forces on the fluid such as the impeller of a pump,
whereas the free vortex naturally occurs in the flow and can be observed in a
drain or in the atmosphere in the form of a tornado.
An equation for the forced vortex can be derived (refer to the textbook) by
applying Newton's law to a fluid element and assuming there are no shear
stresses acting on the fluid (i.e., no relative motion between adjacent particles).
The resulting equation can be expressed between two points:
h2 h1 =
where,
(r
2g
2
2
r12
(6.1)
In this experiment, a forced vortex is created by a rotating plate fitted with blades;
this is very similar to the action of the impeller of a pump.
34
6.2
Winter 2010
Objective
Experimental Procedure
1. Record the initial level in the piezometer tubes. Check whether the same
reading is obtained at both ends of the tube bank. If not, the leveling
screws should be adjusted so that the readings are the same.
2. Turn on the transformer switch.
3. Set the dial at 20. Allow enough time for the levels in the tubes to become
steady. Record the levels on the piezometers.
4. Determine the rotational speed with the tachometer.
5. Repeat steps 5 and 6 for the remainder dial settings. For the highest dial
setting, ensure that none of the piezometer levels are off scale.
6. Set the dial to zero and turn the switch off.
7. Record the distances from the axis of rotation to the piezometer tapping
points.
6.4
Calculations
1. For each run, plot the piezometric head reading against the radial
distance. Draw a smooth curve through the data points, extending it to r =
0 (i.e. y-axis).
2. For each run, plot on graph paper the difference in piezometric head
between each measured value and the value at r = 0, against r2. Draw a
straight line through the points for each run.
3. Find the slope of each line and compare it to its theoretical value from the
given equation:
p
2 2
r
+ z + z =
0 2g
p
where,
(6.2)
35
Winter 2010
Report
36
Winter 2010
FORCE VORTEX
Tube
4
Tube
5
Tube
6
Tube
7
Tube
8
Tube
9
Tube
10
Head
pressure
[in]
Dial 30
Head
pressure
[in]
Dial 40
Head
pressure
[in]
Dial 50
Head
pressure
[in]
Dial
MAX:60
37
Tube
7
Tube
8
Tube
9
Tube
10
C.
Winter 2010
REFERENCES
Beatriz, M.P. (2006). Guide for Writing Laboratory Reports, Department of Civil
Engineering, University of Ottawa, Ottawa, Canada, September 2003 (rev.
August 2006).
Daugherty, R.L. and Franzini, J.B. (1979). Fluid Mechanics with Engineering
Applications, 7th Edition, McGraw-Hill Inc, New York, N.Y.
Streeter, V.L. and Wylie, E.B. (1985). Fluid Mechanics, 8th Edition, Mcgraw-Hill
College, New York, N.Y.
Vennard, J.K. and Street, R.L. (1982). Elementary Fluid Mechanics, 6th Edition,
John Wiley & Sons Inc, New York, N.Y.
38