Fluid Mechanics Presentation
Fluid Mechanics Presentation
Fluid Mechanics Presentation
FLUID MECHANICS
DURATIO
N
Week 1-2
THEME
Fluid properties
(Learning Outcome 1)
Definition of fluid
Introductory
dimensional analysis
Measures of fluid
mass and weight
Ideal gas law
Viscosity
Compressibility of
fluids
Vapour pressure
Surface tension
Solving problems
ASSIGNME
NT/
TEST/PROJ
ECT
Assignment
I
Test I
COMPLETI
ON DATE*
06/08/2010
10/09/2010
Week 3-4
Fluid statics
(Learning
Outcomes 2 )
Pressure at a point
Basic equation for
pressure field
Pressure variation
in a fluid at rest
Measurement of
pressure
Manometry
Hydrostatic force
on a plane surface
Solving problems
Assignment
I
Test I
06/08/2010
10/09/2010
Week
6
5- Flow of
fluids(Learning
Outcomes 3 )
Conservation of
mass
Conservation of
energy
Conservation of
momentum
Assignme
nt II
Test I
03/09/201
0
10/09/201
0
Week 9-10
Flow measurement
(Learning
Outcomes 5 )
Introduction
Rotameter
Venturi meter
Pitot tube
Orifice meter
Vortex flow meter
Assignment
III
Test II
01/10/2010
22/10/2010
Week
15
15/10/2010
22/10/2010
Course objectives
By the end of this module students should be able to:
1. Have a clear understanding of the basic properties of fluid
2. Develop a clear understanding of pressure and its
variation throughout a fluid
3. Perform mathematical calculations for fluid at rest
4. Apply the principles of conservation of mass, momentum and
energy to fluid flow
5. Apply the basic principles of fluid flow to flow of viscous,
incompressible fluids in pipes
6. Perform flow measurements calculations
7. Design centrifugal pumps
Note: If time allows, we will do flow of compressible fluids
6
Area
Deformed
position
Dimensional homogeneity
Now consider an object initially moving at velocity Vo. A
constant acceleration force is applied. Determine the
variation of velocity and distance with time.
Vo
dV Vt Vo
a
dt
t to
Distance,
S
Expression for
acceleration
Vt
Vt Vo
t
Since to=0
Vt = Vo + at
dS
V
dt
ds
Vo at
dt
at 2
S Vo t
2
dS
(Vo at )dt
Example 1.1
A commonly used equation for determining the volume rate flow,
Q, of a liquid through an orifice located inside a tank is
Q 0.61A 2 gh
11
12
Density, kg/m
1000
990
980
970
960
950
0
20
40
60
80
100
Temperature, C
16
17
1.5 Viscosity
The properties of density and specific weight discussed in
section 1.3 are measures of heaviness of a fluid. Another
very important property of a fluid is its viscosity. The viscosity
of a fluid is a measure of resistance to flow under applied
shear stress ,, as show in Fig. 1.3.
Consider
Moving plate, u=V
Moving
plate
Force, F
Velocity
profile
Fluid
Fixed plate
Force
Fig. 1.3 (a) Fluid in shear between
parallel plates
y
V
h
y
u V
h
Eq. 1.7
du V
dy h
19
Newtonian fluids
du
dy
Experimentally determined
du
V
dy
h
Introducing proportionality
constant
Kinematic viscosity
20
oil
H2O, 15 oC
Shearing stress
H2O, 37 oC
1
Air, 15 oC
du
dy
stic
a
l
p
m
a
gh
n
i
B
Shearing stress
ar
e
Sh
g
n
i
nn
i
h
t
n
a
i
on
t
w
Ne
r
ea
h
S
ing
n
e
k
thic
ln a b ln T
D exp( B T )
or
exp(a b ln T )
Andrade Equation
or
CT 2
T S
Sutherland
Equation
23
Dynamic viscosity
Glycerin
Oil
Water
Air
Hydrogen
Temperature
Example 1.3
A dimensional combination of variables that is important in the
study of viscous flow through pipes is called Reynolds number,
NRe, defined as VD/, where is the fluid density, V the mean
fluid velocity, D the pipe diameter and the fluid viscosity. A
Newtonian fluid having viscosity of 0.38 N.s/m 2 and specific
gravity of 0.91 flow through a 25-mm-diametre pipe with a
velocity of 2.6 m/s. Determine the value of NRe. Show how NRe ,
varies with temperature.
25
Solution
N Re
0.38 N .s / m 2
= 156 (kg.m.s-2)/N
But 1 N= 1 kg.m.s-2
Therefore
N Re
NRe= 156
VD
or
D exp( B / T )
N Re
VD
exp(a b ln T )
Example 1.4
The velocity distribution for flow of a Newtonian fluid (see figure)
between two wide, parallel plates is given by the equation:
3
y
u V 1
2
h
Where V is the mean velocity. The fluid has a viscosity of 0.04
lb.s/ft2. When V=2 ft/s and h=0.2 in. Determine: (a) the shearing
stress acting on the bottom wall, and (b) the shearing stress
acting on a plane parallel to the walls and passing through the
centreline (mid plane)
2
27
Solution
du
dy
Newtonian fluid
du 3
2y
V 0 2
dy 2
h
3 Vy
2 h2
du
3V
3V
2 h
dy
h
h
bottom
wall
3
y
u V 1
2
h
3V 0.04 lb.s / ft 2 3 2 ft / s
du
3V
2 0 0
dy
h
mid plane
Given
expression
3V
0
h
Calculation of shearing
stress at the bottom
Calculation of shearing
stress mid plane
28
dp
Ev
d (m / ) /( m / )
E v d
d /
dp
dp
d
d
d d
m( 2 )
2
m/
-Liquids
P
cons tan t
dP
cons tan t
d
Taking derivative
But, we know
dP E v
Ev
P
cons tan t
P
cons tan t
k
Cp
k
Cv
Ev= P
Isothermal
process
Isentropic process:
Frictionless and no heat
exchange with surrounding
31
dp (const.)k k 1 d
Taking derivative
dp
(const.)k k 1
d
P
kP
k 1
k
dp E v kP
E v kP
(isentropic process)
32
Example:
A cubic foot of He at an absolute pressure of 14.7 psi is
compressed isentropically to 0.5 ft3. What is the final pressure?
Solution:
For isentropic compression, we have;
Pf
Pf Pi
i
Pi
k
k
i
f
2
Pf 14.7
1
1.66
Pf = 46.5 psi.
k=1.66
33
Question:
How will the shape of droplets of liquid placed gently on a
surface vary with the size of the droplet?
Answer
For the small droplets such as the
first two on the left hand side, the
droplets adapt a shape that is almost
perfectly spherical, because in this
configuration, there is the least
surface area for a given volume
Capillary action
r 2 h 2r ( cos )
2 cos
h
r
Example
Pressures are sometimes determined by measuring the height of a
column of liquid in a vertical tube. What diameter of clean glass tubing
is required so that the rise of water at 200o C in a tube due to capillary
action is less than 1.00 mm?
Solution
h
2 cos
r
2 cos
r
h
so that
Property of water at 20oC is; = 0.0728 N/m,
= 9.789 kN/m3
2 * 0.0728 * cos 0
r
0.0149
3
3
9.789 * 10 * 1.0 * 10
Ps xs
s
Pz xz
y
X
xyz
Pz xy
40
Assumptions
Fluid at rest
Zero shearing
Fluid moves as rigid body: no relative motion between
adjacent elements
Applying Newtons second Law: F=ma
Y-direction
Sum of forces in the y direction= (pressure)*(small area normal
to pressure)- (force opposing the pressure in the y direction)=
(volume of the element)*(density of the fluid)*acceleration)
xyz
Fy Py xz Ps xs sin 2 a y
41
Z-direction
Sum of forces in the z direction= (pressure)*(small area
normal to pressure)- (force opposing the pressure in the z
direction) (weight of the element) = (volume of the
element)*(density of the fluid)*acceleration)
xyz
xyz
Fz Pz xy Ps xs cos 2 2 a z
It can easily be shown that
z s sin
y s cos
Therefore
Py Ps a y
y
2
Z-direction
Y-direction
z
Pz Ps ( a z )
2
42
Taking limits
y 0, z 0, x 0
Py = Ps
and Pz = Ps
Example
A hydraulic jack consists of an assembly of given dimensions. A
load of 20 kN (P2) must be supported by the jack. What should
be the force applied?
Hydraulic jacks use the
principle of transmission of
pressure.
43
Area A2 is given by
(0.05) 2
A2
0.001963 m 2
5
(0.015) 2
Force applied 10187838 * *
1800 N
4
44
z
P
P
z
z
P P a P a P a
Force acting downwards
due to weight
weight g (z a)
45
P a g (z a )
P g z
P
g
z
Incompressible fluids
P
g
z
P g z
P
g
z
With incompressible fluids
(liquids), density is constant
P Po gz
46
P gz Po
P gz
If we measure pressure
relative to atmospheric
pressure (gauge
pressure, then Po~0
Example
Calculate the pressure (gauge) for a column of
1. 0.2 m water
2. 0.2 gasoline (specific gravity=0.75)
3. 0.2m Mercury (specific gravity=13.6)
47
Solution
gz = 1000*9.81*0.2 = 1962 N/m 2
gz = 0.75*1000*9.81*0.2 = 1472 N/m2
gz = 13.6*1000*9.81*0.2 = 26683 N/m 2
Compressible Fluids
Recall, P
P
RT
P
P
*g
z RT
RT
P
g 2
P
RT z1
Integrating
P2
g
z 2 z1
ln
P1 RT
Isothermal,
T=const.
P1
48
g z 2 z1
P2 P1 exp
RT
g z 2 z1
P2 P1 exp
RT
49
1. Barometer
-Used for measuring atmospheric
pressure
h
-Pressure balance
3. U- Tube Manometer
Using a "U"-Tube enables the pressure of both liquids and gases to be measured with
the same instrument. The "U" is connected as in the figure below and filled with a fluid
called the manometric fluid. The fluid whose pressure is being measured should have a
mass density less than that of the manometric fluid and the two fluids should not be able
to mix readily - that is, they must be immiscible.
-Pressure in a continuous static
fluid is the same at any horizontal
level so
52
If the fluid being measured is a gas, the density will probably be very low in
comparison to the density of the manometric fluid. In this case the term gh1
can be neglected, and the gauge pressure give by
53
6. Choice Of Manometer
Care must be taken when attaching the manometer to vessel, no
burrs must be present around this joint. Burrs would alter the flow
causing local pressure variations to affect the measurement.
Some disadvantages of manometers:
Slow response - only really useful for very slowly varying
pressures - no use at all for fluctuating pressures;
For the "U" tube manometer two measurements must be taken
simultaneously to get the h value. This may be avoided by using
a tube with a much larger cross-sectional area on one side of the
manometer than the other;
It is often difficult to measure small variations in pressure - a
different manometric fluid may be required - alternatively a
sloping manometer may be employed; It cannot be used for very
large pressures unless several manometers are connected in
series;
56
For very accurate work the temperature and relationship
Example of a Piezometer
57
PB = PC
PB = PA + gh1
PB = PAtm + man gh2
We are measuring gauge pressure so
Patm = 0
PA = man gh2 - gh1
58
PC = P D
PC = PA + g ha
PD = PB + g (hb - h) + man g h
PA - PB = g (hb - ha) + hg( man - )
= 990 x9.81x(0.75-1.5) + 0.5x9.81 x(13.60.99) x 103
= -7284 + 61852
= 54 568 N/m2 (or Pa or590.55 bar)
60
61
If the movement in the large cylinder is ignored the term (d/D)2 will disappear
62
Pressure = h
= spec gravity of water
h = depth of water
C = Center of Mass of Gate
CP = Center of Pressure on Gate
Fr = Resultant Force acts at CP
63
Specific weight=
h
FR
P
For a horizontal surface such as shown above, resultant force FR= PA,
where A is the x-sectional area. In terms of fluid depth we have
FR= h A
64
Fr h A y sin A
A
Fr sin y A
A
y A y c A
FR = A YC sin
or FR = A Hc
O
Fr Yr y F sin y A
2
But we know
Fr AYc sin
Yr
2
y
A
A
Yc A
Ix
Yr
Yc A
67
I x I xc AY
2
c
I xc
Yr
Yc
Yc A
Fr X r x F sin xy A
A
Xr
xyA
A
Yc A
I xy
Yc A
Xr
I xyc
Yc A
Xc
68
Parallel axis
theorem
Example
What is the Magnitude and Location of the
Resultant force of water on the door?
W = 62.4 lbs/ft3
Water Depth = 6 feet
Door Height = 4 feet
Door Width = 3 feet
69
Solution
Important variables:
Xc = 1.5
Magnitude of Resultant
Force:
A = 4 x 3 = 12
F R = W A HC
Ixc = (1/12)bh3
HC and Yc = 4
= (1/12)x3x43 = 16 ft4
Location of Force:
YR = (Ixc / YcA) + Yc
YR = (16 / 4x12) + 4 = 4.333 ft
down
XR = Xc (symmetry) = 1.5 ft from
the corner of the door
70
Example
The 4-m-diameter circular gate shown below is located in the inclined
wall of a large reservoir containing water ( =9.80 kN/m3). The gate is
mounted on a shaft along its horizontal diameter. For a water depth of
10 m above the shaft determine the magnitude and location of the
resultant force exerted on the gate by the water.
10 m
c
yr
FR
shaft
71
Solution
To determine the magnitude of the force on
water we apply;
FR = A Hc
= 9.80*10*(3.13*22)
= 1230x103 N= 1.23 MN
Centre of
pressure
Xr
I xyc
Yc A
Xc
I xc
Yr
Yc
Yc A
72
I xc
Yr
Yc
Yc A
I xc
R 4
4
( / 4)(2) 2
10
Yr
73
If the surface is curved the forces on each element of the surface will not be
parallel and must be combined using some vectorial method.
It is most straightforward to calculate the horizontal and vertical components and
combine these to obtain the resultant force and its direction. (This can also be
done for all three dimensions, but here we will only look at one vertical plane).
In the diagram above the liquid is resting on top of a curved base AB.
74
Horizontal forces
Considering the horizontal forces, none can act on CB as there
are no shear forces in a static fluid so the forces would act on
the faces AC and AB as shown below.
76
77
Resultant force
The overall resultant force is found by combining the vertical
and horizontal components vectorialy,
RH
RV
78
Example
A 4 m long curved gate is located in the side of a reservior
containing water as shown below. Determine the magnitude of
the horizontal and vertical components of the force of the
water on the gate. Determine also the resultant force and the
angle it makes with the horizontal.
Solution
Water
9m
A
x
CG
3m
F1
FH
W
C
FV
Free-body diagram
79
Curved surface
Area above
curved surface
=9.81*[(3.14*9/4)+ 3*6)]*4
= 983.6 kN
Length of curved
gate
FV = 983.6 kN
80
FR = ( FH + FV ) 0.5
= (882.9 + 983.6)0.5
=1321.7 kN
= tan-1(FH/FV)
=tan-1(882.9/983.6)
=
81
a x
x
a y
y
For z-axis
p
a z
z
z
y
x
P
P2 1
P3
Constant
pressure lines
ay
p
p
dp
dy dz
y
z
Recall
p
a y
y
p
( g az )
z
dp a y dy ( g a z )dz
Along a line of constant pressure dp= o and therefore
ay
dz
dy
g az
84
ir
pressures=is
from
p
z
The
a function of 2
- z.
variables, r and
o
fr
p
p
w
dp
dr r
dz
2
r
z
o
t
^
-P - k=
a
Axis of
rotation
ar=-rw2
85
dp r 2 dr dz
(i)
pressure is
dr
g
Equation reveals that
w r
z
const. surfaces of constant
2g
pressure are parabolic
2
Pressure distribution
in a rotating liquid
dp 2 rdr dz
w 2 r 2
p
z const.
2
p1
P2
P3
P4
P1
w2 r 2
.
2g
P2
P3
P4
y
x
86
Example
The cross section for the fuel tank of an experimental vehicle is shown
below. The rectangular tank is vented to the atmosphere, and a pressure
transducer is located in its side as illustrated. During testing of the
vehicle, the tank is subjected to a constant linear acceleration, ay (a)
Determine an expression that relates ay and the pressure (in lb/ft2) at the
transducer for a fuel with a SG=0.65. (b) What is the maximum
acceleration that can occur before the fuel level drops below the
transducer?
ay
vent
z
Air
z1
Fuel
(2)
(1)
0.75 ft
0.5 ft
Transducer
0.75 ft
87
Solution
For a constant horizontal acceleration the fuel will move as a rigid body.
The slope of the fuel surface can be expressed as
ay
dz
dy
g
Since az=0
For an arbitrary value of ay, the change in depth z1, of liquid on the right side
of the tank can be found from the Eq.
ay
z1
0.75 ft
g
Since there is no pressure in the vertical z direction, the pressure along the
wall varies hydrostatically with depth. Thus, the pressure at the transducer
is given by
h is the depth of fuel above the transducer
P = h
P= (0.65)(62.4)[0.5-(0.750(ay/g )
= 20.3 30.4 (ay/g)
For z1 0.588 ft
(b) The limiting value for ay ( when the fuel level reaches
the transducer) can be found from the Eq.
z1= (0.75) (ay/g)
For z1 = 0.5, acceleration becomes (ay)max
0.5= (0.75) (ay)max /g)
(ay)max = 2g/3
89
90
Fig. 3.1 (a) System and its surrounding (b) transfers to and from the
system
91
92
93
95
96
Example
Let us now examine the situation illustrated in the figure below.
A tank initially contains 1000 kg of brine containing 10% salt by
mass. An inlet stream of brine containing 20% salt by mass,
flows into the tank at a rate of 20 kg/min. The mixture in the
tank is kept uniform by stirring. Brine is removed from the tank
via an outlet pipe at a rate of 10 kg/min. Find the amount of salt
in the tank at any time t, and the elapsed time 97
when the amount
of salt in the tank is 200 kg.
98
99
100
Recall
Xin Xout = Xsystem
In the present situation, X, refers to energy. Application of this equation over small time
period gives
dMin (e+P/+gz+u2/2)in dMout (e+P/+gz+u2/2)out + dQ-dW= d[M(e+gz+u2/2)]system
(3.3)
Now, let us consider a general system given below
Steady-state
energy balance
101
102
103
104
105
106
107
108
109
m1
m2
m3
m4
m5
mn
Momentum= mu
u1
u2
u3
u4
u5
u dm
un
110
111
112
D 2
d
( p1 p 2 ) wDL ()
4
dt
Momentum balance in
the right direction
113
Two scenarios
(a)Control surface moving with trolley. When the control surface
moves at the same velocity as the trolley, an observation
perceives water entering the system across BC not with velocity
u but with a relative velocity (u-v) so that the convection of mass
into the control volume is A(u-v).
Stationery
Observermoving system
115
A
2
a
(u v)
m
116
dm d
dm
Au
( AL)
Av
dt dt
dt
since dL/dt= v
Momentum balance gives
Rate of addition = Rate of increase
Au 2
d
dv
dm
dL
(mv ALu ) m v
Au
dt
dt
dt
dt
dm
Auv =
= ma v
since dv/dt=a and dL/dt=v
dt
From mass balance equation dm/dt = A(u-v)
118
Force, F
Solution
Taking the right direction as
positive and assuming fluid
mass flow as m the rate of
transfer of momentum into
the system by convection is
mu
120
121
124
125
u12 p1 u 22 p 2
5.1
126
127
Pitot tube
The Pitot tube (named after Henri Pitot in 1732) measures
a fluid velocity by converting the kinetic energy of the
flow into potential energy. The conversion takes place at
the stagnation point, located at the Pitot tube entrance
(see the schematic below). The device is based on
Bernoulli principle and is used in finding the velocity of a
moving craft such as a boat and airplane.
128
129
Pitot-static tube
For pipe flow, we make use of pitot-static tube as shown
below for measuring velocity at
130
3. Venturi meter
A venturi meter is a tube with a constricted throat that increases
velocity and decreases pressure. They are used for measuring the
flowrate of compressible and incompressible fluids in pipeline
Contracting
section
Expanding
section
Pipe
Pipe
Throat
V1
D1
High-pressure tap
p1
D2
V2
D1
Low-pressure tap
p2
131
132
133
134
Problem:
Water is flowing in a schedule 40 pipe having nominal D1 = 4 inch with a flowrate of
325 gal/min at 60F. The pressure differential in the manometer connected to the two
taps of the venturi meter is 50 inch Hg. Calculate (a) throat diameter D2 of the venturi
meter to the nearest 1/8 inch and (b) power consumed by the venturi meter (assume
10% pressure loss).
135
136
137
4. Rota Meter
Annular
area a
p2
Scale
z2
Tapered
tube
Float
Flow
Mg
z1
p1
139
140
141
2 Mg 1
Q A u1 A
2
A
A 1
a
142
Buoyancy Force:
Vdisp-W = 4/3 x x R3
FB = W x Vdisp-w
FB = 62.4 x 33.51
FB = 2091.024 lbs up
Sum of the Forces:
Fy = 0 = 500 - 2091.024 + T
144
FR = A YC sin
or FR = A Hc
YR = (Ixc / YcA) + Yc
XR = (Ixyc / YcA) + Xc
but for a rectangle or
circle: XR = Xc
For 90 degree walls:
F145
R = A Hc
Surrounding S
System V
Loss through
outlet(s)
Accumulation
(or depletion)
Boundary B
Addition through
inlet(s)
V=1m
Po= 1 bar
To vacuum pump
V=0.001 m3/s
(independent of
pressure)
146
dMin
Work done,
dW
system
dMout
Heat added, dQ
147
W (Work)
system
Outlet
q (heat)
148
Vent
Storage
tank
40 ft
flow
Vent
Supply
tank
1
3
4.5 ft
4 ft
pump
149
P1
Direction
of flow
P2
150
System is moving
with velocity v
B
Water
supply
jet
v
Stationery
observer
C
151
System is
stationery
B
Water
supply
jet
L
m
Stationery
observer
152
shield
m
u
Force, F
153
A2
A1
1
154
water
1
Stagnation point
155
movable
Pipe wall
u1
Velocity
profile
liquid
Stagnation point
156
157
20 kg/min
Salt content
20% by mass
Tank, initial
content 100kg
10 kg/min
158