Fluid Dynamic
Fluid Dynamic
Fluid Dynamic
FLUID DYNAMICS
Characteristics of Fluid Flow
A B A B
A B A B
A B A B
A B A B
� Ideal Fluids
– the effect of shear stress is ignored
because there is no viscosity.
� Real Fluids
– shear stress occurred between layers of
moving fluid as a result of viscosity.
Compressible and Incompressible Flow
� Laminar Flow
- Re < 2100
- The particle of fluid moving in an orderly
manner and retaining the same relative
position in successive cross sections
or
- Flow is said to be laminar when adjacent
fluid layers move at same velocity and
paths of individual particles of fluid do
not cross each other.
- Occur at low velocities and high viscosity
� Transition Flow
- 2100 < Re < 4000
- The characteristic of this flow is in
between laminar and turbulent flow.
� Turbulent Flow
- Re > 4000
- The particle of fluid no longer moved in
an orderly manner but occupied different
relative position in successive cross
sections
or
- Flow is turbulent when streamlines cross
each other and mixing of fluid flow occur.
- Occur at high velocities and low viscosity
Continuity Equation
Bernoulli’s Equation
Thus
-(p/s)() - gsinq = aS
-(p/s) - gsinq = aS
Divide with
-(1/)(p/s) - gsinq - V/t = 0
Knowing that
V/t = (V/s)(s/t) + v/t
Therefore
-(1/)(p/s) - gsinq - (V/s)(s/t) = 0
-(1/)(p/s) – g(z/s) - (V/s)V = 0
(1/)p + gz + (V)V = 0 Euler
equation
Rearrange
(1/g)P + V2/2g + z = constant
(1/g)P + V2/2g + z = H
where
(1/g)P = pressure head
V2/2g = velocity head
z = elevation head
H = total head
Free Jets
Thus
H = V22/2g Torricelli
Equation
V2 = (2gH)1/2
From experiment
VACT = CVVT
where
VACT = Actual velocity
VT = Theoretical velocity
CV = Velocity coefficient
V2ACT = CV(2gh)1/2
Flowrate
Q = V2A2
= (2gh)1/2A2
Actual Flowrate
QACT = V2ACTA2ACT
where
CD = CCCV = flowrate coefficient
� Vena contracta coefficients
Orifice in Pipe
Assume
D1 = Flow diameter near point A
D2 = Flow diameter near point B
D = Orifice diameter
Bernoulli equation
(1/g)P1 + V12/2g + z1 = (1/g)P2 + V22/2g +
z2
Actual Velocity
V2ACT = CVV2
= CV{(2g[P1 – P2]/g)/(1 – CC2[(D)4/
(D1)4])}1/2
Actual flowrate
QACT = V2ACTA2
= CV{(2g[P1 – P2]/g)/(1 – CC2[(D)4/(D1)4])}1/2
{CCA}
= CD{(2g[P1 – P2]/g)/(1 – CC2[(D)4/
(D1)4])}1/2 A}
Pressure different, P1 – P2
P1 - gh1- gr = P2 - g(h1+h2) - gr +
1gh2
P1 – P2 = gh1 + gr - g(h1+h2) - gr +
1gh2
Divide with g
(P1 – P2)/ g = h1 + r - (h1+h2) - r + s1h2
(P1 – P2)/ g = s1h2 – h2
Therefore
QACT = CD{2g(s1h2 – h2)/(1 – CC2[(D)4/
(D1)4])}1/2 {A}
Venturi Meter
Thus
V12(1 – d14/d24)/2g = (P1 – P2)/g + (z1 – z2)
From manometer
PX = P1 + Ogz1
and
PX = P2 + Og(z2 - h) + mgh
Thus
P2 – P1 = Og(z1+h – z2) - mgh
Therefore
V12(1 – d14/d24)/2g = (Og(z1+h – z2) - mgh)/g
+ (z1 – z2)
Flowrate
Q= VA
= {[(Og(z1+h–z2) - mgh)/g + (z1–z2)]/
[(1 – d14/d24)/2g]}1/2 A1
At point 2
V2 = 0 (stagnation point)
z1 – z2 = 0
Thus
V12/2g = (P2 – P1)/g
V1 = [2(P2 – P1)/]1/2
From manometer
P2 - P1 = 1gR - gR
Therefore
V1 = [2(1gR - gR)/]1/2
Flowrate
Q = VA
= [2(1gR - gR)/]1/2A1
Weir
� Bernoulli equation
(1/g)P1 + z1 + V12/2g = (1/g)P2 + z2 +
V22/2g
Assume point 1 is at the water surface while
point 2 is a small fluid element with thickness
y at depth y from the water surface at the
exit/outlet
At point 1
P1 = 0
V1 = 0
At point 2
P2 = 0
Flowrate
Q = V2A2
= ò V2A
= ò[2g(y)]1/2 (by)
= [2g]1/2 b(y)3/2(2/3)
= (2/3)[2g]1/2 b(h)3/2
Solution :
Bernoulli equation
(1/g)P1 + z1 + V12/2g = (1/g)P2 + z2 +
V22/2g
(V12 –V22)/2g = (P2 – P1)/g + (z2 – z1)
Fluid density
= specific weight/gravity
= 9100/9.81
= 928 kg/m3
Thus
[V12 – (V12d14/ d24)]/2g = (P2 – P1)/g
V12[1 – (d14/ d24)]/2g = (P2 – P1)/g
V12 = 2(P2 – P1)/ [1 – (d14/ d24)]
= 2(550k – 735k)/928[1 –
(0.031 /0.0194)]
4
= (-370k)/(-5648.2)
V1 = 8.09 m/s
Flowrate
Q =V1A1
= 8.09(0.0312)/4
= 6.103 x 10-3 m3/s
Example 2
Solution :
Bernoulli equation
(1/g)P1 + z1 + V12/2g = (1/g)P2 + z2 +
V22/2g
(V12 –V22)/2g = (P2 – P1)/g + (z2 – z1)
At point 2
V2 = 0
And
z2 – z1 = 0
Thus
V12/2g = (P2 – P1)/g
From manometer
P2 – P1 = g(2.5) - mg(2.5)
= 1000(9.81)(2.5) - 900(9.81)(2.5)
= 2452.5 Pa
Therefore
V12/2g = (P2 – P1)/g
V12 = 2(P2 – P1)/
2
V1 = 2(2452.5)/1000
= 4.905
V1 = 2.21 m/s
Flowrate
Q = V1A1
= 2.21(0.082)/4
= 0.011 m3/s
The Momentum Equation
� Examples :
- The lift force on an aircraft is exerted by
the air moving over the wing.
- A jet of water from a hose exerts a force
on whatever it hits.
so at the inlet
volume = area x distance
= A1u1t
where
FR = Force exerted on the fluid by any
solid body touching the control
volume
FB = Force exerted on the fluid body (e.g.
gravity)
FP = Force exerted on the fluid by fluid
pressure outside the control volume
� Step in Analysis
- Draw a control volume
- Decide on co-ordinate axis system
- Calculate the total force
- Calculate the pressure force
- Calculate the body force
- Calculate the resultant force
- In the y-direction :
.
FTY = m (u2Y - u1Y)
= Qu2sinq
- In the y-direction :
FPY = –p2A2sinq
- In the y-direction :
FTY = FRY + FBY + FPY
FRY = FTY - FBY - FPY
= Qu2sinq - (–p2A2sinq)
Resultant force on the fluid is given by
FR = (FRX2 + FRY2)1/2
= Q(u2 - u1)
= - Qu1
- In the y-direction :
.
FTY = m (u2Y - u1Y)
= (Q2/A)sinq
From above u1 = u2 = u3
0 = A2 – A3 – A1sinq
as u2 = u3 = u
Q2 = Q3[(1 + sinq)/(1 - sinq)]
Q1 = Q3[(1 + sinq)/(1 - sinq)] + Q3
= Q3[1 + {(1 + sinq)/(1 - sinq)}]
Example 1
The figure below shows a smooth curved vane
attached to a rigid foundation. The jet of water,
rectangular in section, 75mm wide and 25mm
thick, strike the vane with a velocity of 25m/s.
Calculate the vertical and horizontal
components of the force exerted on the vane
and indicate in which direction these
components act.
Solution :
u1 = 25 m/s
Q = 1.875 x 10-3x 25
= 0.0469 m3/s
So force on vane
RX = - FTX = - 233.44 N
RY = - FTY = -1324.6 N
EXAMPLE 2
Solution :
u1 = 50 – 25
= 25 m/s
Q = 7.85 x 10-3x 25
= 0.1963 m3/s
Calculate the total force using the momentum
equation:
FTX = Q(u2cos20 – u1)
= 1000(0.1963)([-25cos 20] – 25)
= -9.519 kN
FTY = Q(u2sin20)
= 1000(0.1963)(25sin 20)
= 1.678 kN
So force on vane
RX = - FTX = 9.519 kN
RY = - FTY = -1.678 kN