0% found this document useful (0 votes)
117 views12 pages

Solution: Chapter 2: 2.1 Cavitation

This document discusses several fluid mechanics problems involving pipe flow and pumps. It examines cavitation in a hose, bottling honey with laminar flow, and filling a boiler with turbulent flow. It also analyzes the best syringe size for delivering a given volume over time. Key equations for pressure drop, flow rate, Reynolds number, and filling time are applied to calculate values for each problem.

Uploaded by

서원혁
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
117 views12 pages

Solution: Chapter 2: 2.1 Cavitation

This document discusses several fluid mechanics problems involving pipe flow and pumps. It examines cavitation in a hose, bottling honey with laminar flow, and filling a boiler with turbulent flow. It also analyzes the best syringe size for delivering a given volume over time. Key equations for pressure drop, flow rate, Reynolds number, and filling time are applied to calculate values for each problem.

Uploaded by

서원혁
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 12

Solution: Chapter 2

2.1 Cavitation
(a) Given:
Q = 3.93 × 10−3 m3 /s h2 − h1 = 3 m
D = 5 cm ρ = 1.0 × 103 kg/m3
L = 50 m ν = 1.0 × 10−6 m2 /s

Find: |ᏼ| from hose inlet to pump intake.


Calculate Re and f:
4Q 4(3.93 × 10−3 )
U = = = 2.0 m/s
πD 2 π (5 × 10−2 )2
UD (2)(0.05)
Re = = = 1.0 × 105
ν 1 × 10−6
Flow is turbulent and Blasius eq. [Eq. (2.2-9)] is applicable.

f = 0.0791 Re−1/4 = 0.0791(1 × 105 )−1/4 = 4.45 × 10−3

Find |ᏼ| from Eq. (2.3–6):


2ρU 2 L 2(1000)(2)2 50
|ᏼ| = f = (4.45 × 10−3 )
D 0.05
= 3.56 × 104 Pa

(b) Find P2 (pump intake) for P1 = 1.01 × 105 Pa (nearly empty pool).

ᏼ = −|ᏼ| = P2 − P1 + ρg(h2 − h1 )
P2 = P1 − |ᏼ| − ρg(h2 − h1 )
= 1.01 × 105 − 3.56 × 104 − (1000)(9.81)(3)
= 3.60 × 104 Pa

Although below atmospheric, P2 > PV = 2.34 × 103 Pa, so cavitation will be absent.
(c) Maximum Q without cavitation.
P1 is lowest when the pool is nearly empty, so P2 is also lowest then. Q = Qmax when
P2 = PV and P1 = 1.01 × 105 Pa.

1
Solution: Chapter 2

Find |ᏼ| at Q = Qmax :


|ᏼ| = P1 − P2 − ρg(h2 − h1 )
= 1.01 × 105 − 2.34 × 103 − (1000)(9.81)(3)
= 6.92 × 104 Pa

Re is now unknown, but Re f can be found and the Prandtl–Kármán eq. used
[Eq. (2.2-7)].
D |ᏼ| UD
f = , Re =
2ρU 2 L ν
 1/2
 D D|ᏼ|
Re f =
ν 2ρL
 1/2
0.05 0.05(6.92 × 104 )
= = 9.30 × 103
1 × 10−6 2(1000)(50)

f = [4.0 log(Re f ) − 0.4]−2 = 4.18 × 10−3
 
π D2 π D2 D |ᏼ| 1/2
Qmax = Umax =
4 4 2ρ f L
 1/2
π (0.05)2 0.05(6.92 × 104 )
=
4 2(1000)(4.18 × 10−3 )(50)
= 5.65 × 10−3 m3 /s

2.2 Bottling Honey


Given: Gravity feed with |h| = 2 m.
Tube length = L = 4m.
Honey properties as in Table 1.2.
Mass flow rate ≥ 17 oz/min.
Find: Minimum tube diameter D.
Volume flow rate:
  
1 lb 0.454 kg
17 oz = 0.482 kg
16 oz lb
  
kg 1 min
w = mass flow rate = 0.482 = 8.04 × 10−3 kg/s
min 60 s
w 8.04 × 10−3
Q = vol. flow rate = = = 5.66 × 10−6 m3 /s
ρ 1.42 × 103
The high viscosity suggests laminar flow. Poiseuille’s equation:
π D4 |ᏼ|
Q=
128 μ L
 
128μLQ 128μLQ 1/4
D4 = , D=
π |ᏼ| π |ᏼ|
|ᏼ| = |P| + ρg|h| = ρg|h| (atm P at both ends)

2
Solution: Chapter 2

 1/4  1/4
128νLQ 128(1.34 × 10−2 )(4)(5.66 × 10−6 )
D= =
π g|h| π (9.81)(2)
= 2.82 × 10−2 m = 2.82 cm.
This is a very reasonable diameter, indicating that the bottling process is feasible.
Check Re to confirm that the flow is really laminar:
4Q 4 (5.66 × 10−6 )
U = = = 9.08 × 10−3 m/s
πD 2 π (2.82 × 10−2 )2
UD (9.08 × 10−3 )(2.82 × 10−2 )
Re = = = 1.91 × 10−2
ν 1.34 × 10−2
This is much smaller than 2100, which confirms the applicability of Poiseuille’s equation.

2.3 Filling a Boiler


Given: V = 7.57 m3 (volume to be transferred)
Copper tubing with D = 2.60 cm and L = 15 m.
h = 3 m
|P| = 4.0 × 105 Pa
(a) Filling time at 20◦ C (ρ = 998 kg/m3 , ν = 1.00 × 10−6 m2 /s)
Assume that the√flow is turbulent and might be affected by roughness. Although Re
is unknown, Re f can be found from what is given, making the Colebrook – White
eq. [Eq. (2.5.2)] a good choice.
D |ᏼ| UD
f = , Re =
2ρU 2 L ν
 
 D D|ᏼ| 1/2
Re f =
ν 2ρL
ᏼ = P + ρgh = −|P| + ρgh.
|ᏼ| = −ᏼ = |P| − ρgh
= 4.0 × 105 − (998)(9.81)(3) = 3.71 × 105 Pa
 1/2
 0.026 0.026(3.71 × 105 )
Re f = = 1.48 × 104
1.00 × 10−6 2(998)(15)
k 1.5 × 10−6
= = 1.56 × 10−5 (k for “drawn metals”, Table 2.2)
3.7D 3.7(0.026)
  −2
1.26 k
f = 4 log √ +
Re f 3.7D
  −2
1.26 −5
= 4 log + 1.56 × 10
1.48 × 104
= 3.91 × 10−3
   1/2
D |ᏼ| 1/2 0.026(3.71 × 105 ) m
U = = −3
= 9.07
2ρ f L 2(998)(3.91 × 10 )(15) s

3
Solution: Chapter 2

π D2 π (0.026)2 m3
Q= U = (9.07) = 4.82 × 10−3
4 4  s
V 7.57 1 min
t= = = 1.57 × 103 s = 26.2 min
Q 4.82 × 10−3 60 s
Check Re to confirm turbulent flow:
UD 9.07(0.026)
Re = = = 2.36 × 105 > 2100.
ν 1.00 × 10−6
(b) Filling time at 52◦ C (ρ = 987 kg/m3 , ν = 5.29 × 10−4 /987 = 0.536 × 10−6 m2 /s).

Re f ∝ ν −1 ρ −1/2
  
 1.00 998 1/2
Re f = 1.48 × 10 4
= 2.78 × 104
0.536 987
k
= 1.56 × 10−5 as before
3.7D
  −2
1.26 −5
f = 4 log + 1.56 × 10 = 3.52 × 10−3
2.78 × 104
U ∝ (ρ f )−1/2
 
998(3.91) 1/2 m
U = 9.07 = 9.61
987(3.52) s
1
t∝
U  
9.07
t = 26.2 = 24.7 min
9.61
Only 1.5 min is saved.

2.4 Syringe Pump


(a) Best syringe size from among 1 ml, 5 ml, and 50 ml.
Total volume to be delivered is
   
100 µL 0.1 mL
V = (30 min) = (30 min) = 3 mL.
min min
The 1 ml syringe is too small.
Pusher velocity needed for other syringes:
4Q
Us =
π D2s
4(0.1 cm3 /min) cm mm
Us (5 ml) = = 8.70 × 10−2 = 0.87
π (1.21 cm)2 min min
4(0.1 cm3 /min) cm mm
Us (50 ml) = = 1.79 × 10−2 = 0.18
π (2.67 cm)2 min min
Either would work, but 5 ml is preferable because its Us is not as close to the limits
of the pump.

4
Solution: Chapter 2

(b) |P| in horizontal tubing


 
UD 4Q D 4Q
Re = = =
ν πD 2 ν π Dν
   
L m3 1 min m3
Q = 100 × 10−6 10−3 = 1.67 × 10−9
min L 60s s
−9
4(1.67 × 10 )
Re = = 2.47 (laminar flow)
π (8.6 × 10−4 )(1 × 10−6 )
128 µLQ
|P| = |ᏼ| = (Poiseuille)
π D4
128(1 × 10−3 )(0.5)(1.67 × 10−9 )
= = 62.2 Pa
π (8.6 × 10−4 )4
As shown, |P| depends only on Q and the tubing dimensions. The syringe size does
not matter.
(c) |P| for outlet above inlet, h2 − h1 = 10 cm.
|ᏼ| (from Poiseuille eq.) will be the same as before.
ᏼ = P + ρgh
−P = −ᏼ + ρgh = |ᏼ| + ρgh (ᏼ < 0)
ρgh = (1 × 103 )(9.81)(0.1) = 981 Pa
−P = 62 + 981 = 1043 Pa = |P|
Moving the liquid upward requires a larger |P|. Because Poutlet is fixed at atmo-
spheric, Pinlet must be larger.

2.5 Flue Gases


Given: D, L, and |P|
Find: w (mass flow rate)
Assume that the flow is turbulent, in √ which case it might be affected by roughness.
Although Re is not known initially, Re f can be calculated, making the Colebrook –
White equation [Eq. (2.5-2)] a good choice.
 
1 1.26 k
√ = −4 log √ + Eq. (2.5-2)
f Re f 3.7D
   3 1/2
 U D D|P| 1/2 D |P|
Re f = =
ν 2ρU 2 L 2ν 2 ρL
−5
μ 2.6 × 10 m2
ν= = = 4.13 × 10−5
ρ 0.63 s
 
 (6.4 × 10−2 )3 (63)
1/2
Re f = = 1.12 × 103
2(4.13 × 10−5 )2 (0.63)(6.1)
k = 4.6 × 10−5 m for steel (Table 2.2)
k 4.6 × 10−5
= = 7.19 × 10−4
D 6.4 × 10−2

5
Solution: Chapter 2

  −2
1.26 k
f = 4 log √ +
Re f 3.7D
  −2
1.26 7.19 × 10−4
= 4 log +
1.12 × 103 3.7
= 7.53 × 10−3
From Fig. 2.6, Re ∼
= 104 , confirming that the flow is turbulent.
Find U from f:
 1/2  1/2
D |P| (6.4 × 10−2 )(63) m
U = = = 8.35
2ρ f L 2(0.63)(7.53 × 10−3 )(6.1) s
Check Re:
UD (8.35)(6.4 × 10−2 )
Re = = = 1.29 × 104 (turbulent, as assumed)
ν 4.13 × 10−5
Find w:
π D2 π kg
w = ρQ = U ρ = (6.4 × 10−2 )2 (8.35)(0.63) = 1.69 × 10−2
4 4 s

2.6 Hydraulic Fracturing


Given Q, D, L, h, and P2 (outlet), find P1 (inlet). Start with Re:
 
UD 4Q D 4Q
Re = = =
ν π D2 ν π νD
0.467 × 10−3 −7 m
2
ν= = 4.75 × 10 for 60◦ C water.
0.983 × 103 s
4 0.21
Re = = 3.75 × 106 (turbulent flow)
π (4.75 × 10−7 )(0.15)
See if roughness is important:
k = 4.6 × 10−5 m for steel (Table 2.2)
k 4.6 × 10−5
= = 3.1 × 10−4
D 0.15
Inspection of Fig. 2.6 indicates that roughness significantly increases f at this Re and k/D.
For turbulent flow with roughness, Eq. (2.5–3) is convenient for finding f:
   1.11
−2
6.9 k
f = 3.6 log +
Re 3.7D
   1.11
−2
6.9 3.1 × 10−4
= 3.6 log + = 3.81 × 10−3
3.75 × 106 3.7

Find |ᏼ| from f:


2ρU 2 L f 4Q m
|ᏼ| = U = = 11.9
D π D2 s

6
Solution: Chapter 2

2(0.983 × 103 )(11.9)2 (2.3 × 103 )(3.81 × 10−3 )


=
0.15
= 1.62 × 10 Pa = 16.2 MPa
7

Find P1 (top of pipe):

ᏼ = P + ρgh
−P = P1 − P2 = −ᏼ + ρgh = |ᏼ| − ρg(h) (ᏼ < 0, h < 0)
ρg|h| = (0.983 × 103 )(9.81)(2.3 × 103 ) = 2.22 × 107 Pa = 22.2 M Pa
P1 = P2 + |ᏼ| − ρg|h|
= 44 + 16.2 − 22.2 = 38 M Pa

Note that even though the downward flow encounters a large resistance in the pipe, P1 <
P2 due to the large ρg|h|.

2.7 Drag Reduction


Given:

U = 3.0 m/s ρ = 850 kg/m3


D = 1.0 m ν = 1.06 × 10−5 m2 /s
L = 100 km

Find: |P|(= |ᏼ|) for maximum drag reduction.


In Example 2.3–2 it was found that Re = 2.83 × 105 and roughness was assumed (implic-
itly) to be negligible. Using k for steel in Table 2.2, k/D = 4.6 × 10−5 . Figure 2.6 con-
firms that roughness effects are indeed small at this Re and k/D.
Approach: Find f iteratively from Eq. (P2.7-1)
Find |P| from f and the given data.
For Re = 2.83 × 105 , Fig. P2.7 gives f = ∼ 6 × 10−4 . An improved value from
Eq. (P2.7-1) is

Re f = (2.83 × 105 )(6 × 10−4 )1/2 = 6.93 × 103

f = [19 log(Re f ) − 32.4]−2 = [19 log(6.93 × 103 ) − 32.4]−2
= 6.07 × 10−4

A second refinement is

Re f = (2.83 × 105 )(6.07 × 10−4 )1/2 = 6.97 × 103
f = [19 log(6.97 × 103 ) − 32.4]−2 = 6.06 × 10−4

which is adequate convergence. Then

2ρU 2 L
|P| = f ∝ f
D

7
Solution: Chapter 2

|P| = |P|0 = 5.59 × 106 Pa and f = f0 = 3.66 × 10−3 without additives (Example
2.3-2).
   −4 
fMDR 6 6.06 × 10
|P|MDR = |P|0 = 5.59 × 10
f0 3.66 × 10−3
= 9.26 × 105 Pa
This is an 83% reduction. Polymeric additives have been used, for example, in the Trans-
Alaska Pipeline.

2.8 Economic Pipe Diameter


(a) Pipe diameter D∗ that minimizes total cost
C = K1 Dn L + K2 Q|ᏼ|
Relate Q|ᏼ| to D:
 
2ρU 2 f L 2ρ f L 4Q 2 32 ρ f LQ2
|ᏼ| = = =
D D π D2 π 2 D5
 −b π b  νD b
4 Q
f ∼
= aRe −b
= a = a
π νD 4 Q

32a π ρLQ ν D
b 3 b b
Q|ᏼ| = 2
π 4 D5 Qb

32a π b
Let c =
π2 4
Q|ᏼ| = cρLν b Q3−b Db−5
Total cost is now
C(D) = K1 Dn L + K2 cρLν b Q3−b Db−5
To minimize C, set dC/dD = 0:
dC
= 0 = nK1 Dn−1 L + K2 cρLν b Q3−b (b − 5)Db−6
dD
nK1 Dn−1 = (5 − b)(K2 cρν b Q3−b Db−6
(5 − b)K2 c b 3−b
Dn−1−b+6 = D5+n−b = ρν Q
nK1
(5 − b)K2 c 1−b b 3−b
= ρ μQ
nK1
 
∗ (5 − b)K2 c 1/(5+n−b) (1−b)/(5+n−b) b/(5+n−b) (3−b)/(5+n−b)
D = ρ μ Q
nK1
(b) Dependence of D∗ on ρ and μ
Use n = 1.5 and b = 0.16.
1−b 0.84
ρ exponent is = = 0.132
5+n−b 6.34
b 0.16
μ exponent is = = 0.025
5+n−b 6.34

8
Solution: Chapter 2

Both are small, implying insensitivity of D∗ to property values.


(c) Dependence of U ∗ on Q

D∗ ∝ Q(3−b)/(5+n−b) = Q2.84/6.34 = Q0.45


Q Q
U ∗ ∝ ∗ 2 ∝ 2(3−b)/(5+n−b)
(D ) Q
2(3 − b) 5 + n − b − 6 + 2b −1 + n + b
1− = =
(5 + n − b) 5+n−b 5+n−b
−1 + 1.5 + 0.16
= = 0.10
5 + 1.5 − 0.16
U ∗ ∝ Q0.10

Whereas D∗ is somewhat sensitive to Q, U ∗ is relatively insensitive.

2.9 Microfluidic Device


(a) Flow-rate ratios.
From flow in = flow out, Q1 = Q5 .
From L2 = L4 , Q2 = Q4
From L3 = L2 /2, Q3 = 2Q2 .

Q1 = Q2 + Q3 + Q4 = Q2 + 2Q2 + Q2 = 4Q2
Q2 1 Q3 1 Q4 1 Q5
= , = , = , =1
Q1 4 Q1 2 Q1 4 Q1
(b) Calculate Q1 = overall flow rate for PA − PD = 3000 Pa, L1 = 1 mm, and square
channels with side length 100 µm.
Let ri j = resistance between nodes i and j
rk = resistance of channel k
There are 3 parallel channels between B and C. With ri ∝ Li ,
1 1 1 1 1 2 1 4
= + + = + + =
rBC r2 r3 r4 r2 r2 r2 r2
With L2 = 4L1 , r2 = 4r1 and
1 1
= or rBC = r1 .
rBC r1
With L5 = L1 , r5 = r1 . Adding the series resistances,

rAD = r1 + rBC + r5 = 3r1

Friction factor for laminar flow:


DH |P| c
f = 2
= c = const.
2ρU L ReH

9
Solution: Chapter 2

For square channels, c = 14.23, DH = a = side length (Table 2.1 with a = b), and
Re4 = U DH /ν = U a/ν. Thus
a |P| cν a2 |P|
= or U =
2ρU 2 L Ua 2cμ L
Apply to segment 1:
|P|AD
|P|1 = = 1000 Pa because r1 = rAD /3.
3
L1 = 1 × 10−3 m
a = 1 × 10−4 m
(1 × 10−4 )2 1000 m
U1 = −3 −3
= 0.351
2(14.23)(1 × 10 ) 1 × 10 s
m3
Q1 = U1 a2 = (0.351)(1 × 10−4 m)2 = 3.51 × 10−9
s

2.10 Murray’s Law


(a) Radius that minimizes total energy expended (Et ) for one vessel.

Et = Q|ᏼ| + bπ R2 L
8μLQ
|ᏼ| = (Poiseuille)
π R4
8μLQ2
Et = + bπ R2 L ≡ AR−4 + BR2
π R4
Minimize Et with respect to R, for constant A = 8μLQ2 /π and B = 6π L:
dEt
= −4AR−5 + 2BR = 0
dR
 
2A 16μLQ2 1 16μQ2
R6 = = =
B π bπ L bπ 2
  1/6
16μQ2
R= (independent of L)
bπ 2
(b) Obtain Murray’s law by applying result from part (a) to each vessel at a bifurcation.
Q1 = Q2 + Q3 (vessels labeled as in Fig. P2.10)
Ri ∝ Q1/3
i or Qi ∝ R3i (optimum found above)
R31 = R32 + R33
(c) Extend Murray’s law to generations j and k
Let Qi j = flow rate in vessel i at generation j
Same total flow rate at each generation:


Nj

Nk
Qi j = Qik
i=1 i=1

10
Solution: Chapter 2

Energy minimization:
Qi j ∝ R3i j
Combine:

Nj

Nk
R3i j = R3ik
i=1 i=1

2.11 Open-Channel Flow


(a) Dimensional analysis
Governing quantities with dimensions:
[U ] = LT−1 [g] = LT−2
[H ] = L [ρ] = ML−3
[W ] = L [μ] = ML−1 T−1
[k] = L [S] = dimensionless (slope)
7 quantities − 3 dimensions = 4 dimensionless groups
For N1 choose U (the unknown) for the numerator and g and H for the denominator:
U
N1 =
ga H b
[U ] = LT−1 = [ga H b ] = (LT−2 )a Lb = La+b T−2a
L: 1 = a + b
T: −1 = −2a
a = 1/2 = b
U
N1 = (= Fr1/2 )
(gH )1/2
Note: The reason for preferring H over W or k is that U may be expected to be inde-
pendent of W for W/H → ∞ (as in a parallel-plate channel) and independent of k for
k/H → 0 (smooth surfaces), whereas H will always affect U. That g is a key quantity
in a gravity-driven flow is obvious.
A known group with ρ and μ is Re:
UHρ
N2 = = Re
μ
The remaining two groups must include W and k. Let
W
N3 = = φ aspect ratio
H
k
N4 = roughness parameter
H
Conclusion:
U
= F (S, Re, φ, k/H )
(gH )1/2

11
Solution: Chapter 2

U will certainly increase with increasing slope and decrease with increasing rough-
ness, but the form of F is largely unknown.
(b) Force balance for length L of channel
Let τw be the shear stress averaged over the solid surfaces. Then
shear force = τw (side area + bottom area) = τw (2H + w)L
gravitational force = (fluid mass) · (grav, accel. in flow direction)
= (ρLHW )g sin θ = ρLHW gS
ρLHW gS = τw (2H + W )L
 
HW ρgH Sφ
τw = ρgS =
2H + W 2+φ
(c) Force balance using hydraulic diameter
4A
Let DH = with A = cross-sectional area
C
C = wetted perimeter (solid only).
For the rectangular channel,
4HW 4H φ
DH = =
2H + W 2+φ
The general expression for τw is
ρALgS = τwCL
ρgSA ρgSDH
τw = =
C 4
(d) Friction factor independent of flow rate (f independent of Re)
 
2τw 2 ρgSDH gSDH
f = = =
ρU 2 ρU 2 4 2U 2
Arrange as in part (a), but with DH instead of H:
 1/2
U S
1/2
=
(gDH ) 2f
It is seen now that U ∝ S 1/2 . Although independent of ReH and presumably φ, in a
rough channel f will depend on k/DH [Eq. (P2.11–2)].

12

You might also like