Table 11.3 Properties of saturated water (Meyer et al., 1993; Harvey et al., 2000; IAPWS, 1994).
T |
p |
ρf |
hfg |
cpf |
f × 106 |
kf × 103 |
Prf |
σ × 103 |
βf × 106 |
(K) |
(bar) |
(kg/m3) |
(kJ/kg) |
(kJ/kg K) |
(Ns/m2) |
(W/mK) |
|
(N/m) |
(K–1) |
273.16 |
0.006177 |
999.8 |
2502 |
4.22 |
1790.6 |
561 |
13.47 |
75.65 |
-68 |
275 |
0.006985 |
999.9 |
2497 |
4.214 |
1681.8 |
564.5 |
12.55 |
75.39 |
-35.5 |
280 |
0.009918 |
999.9 |
2485 |
4.201 |
1433.9 |
574 |
10.63 |
74.68 |
43.6 |
285 |
0.01389 |
999.5 |
2473 |
4.193 |
1239.4 |
583.5 |
8.91 |
73.95 |
112 |
290 |
0.0192 |
998.8 |
2461 |
4.187 |
1083.7 |
592.7 |
7.66 |
73.21 |
172 |
295 |
0.02621 |
997.8 |
2449 |
4.183 |
957.9 |
601.7 |
6.66 |
72.46 |
226 |
300 |
0.03537 |
996.5 |
2438 |
4.181 |
853.8 |
610.3 |
5.85 |
71.69 |
275 |
305 |
0.04719 |
995 |
2426 |
4.18 |
766.9 |
618.4 |
5.18 |
70.90 |
319 |
310 |
0.06231 |
993.3 |
2414 |
4.179 |
693.5 |
626 |
4.63 |
70.11 |
361 |
320 |
0.1055 |
989.3 |
2390 |
4.181 |
577.0 |
639.6 |
3.77 |
68.47 |
436 |
340 |
0.2719 |
979.5 |
2342 |
4.189 |
422.0 |
660.5 |
2.68 |
65.04 |
565 |
360 |
0.6219 |
967.4 |
2291 |
4.202 |
326.1 |
673.7 |
2.03 |
61.41 |
679 |
373.15 |
1.014 |
958.3 |
2257 |
4.216 |
281.7 |
679.1 |
1.75 |
58.91 |
751 |
400 |
2.458 |
937.5 |
2183 |
4.256 |
218.6 |
683.6 |
1.36 |
53.58 |
895 |
420 |
4.373 |
919.9 |
2123 |
4.299 |
186.7 |
682.5 |
1.18 |
49.41 |
1008 |
440 |
7.337 |
900.5 |
2059 |
4.357 |
162.8 |
678 |
1.05 |
45.10 |
1132 |
460 |
11.71 |
879.5 |
1989 |
4.433 |
144.3 |
670.2 |
0.955 |
40.66 |
1273 |
480 |
17.9 |
856.5 |
1912 |
4.533 |
129.7 |
659 |
0.892 |
36.11 |
1440 |
500 |
26.39 |
831.3 |
1825 |
4.664 |
117.7 |
643.9 |
0.853 |
31.47 |
1645 |
520 |
37.69 |
803.6 |
1730 |
4.838 |
107.6 |
624.6 |
0.833 |
26.78 |
1909 |
540 |
52.37 |
772.8 |
1622 |
5.077 |
98.8 |
600.1 |
0.835 |
22.08 |
2266 |
560 |
71.06 |
738 |
1499 |
5.423 |
90.8 |
570.1 |
0.864 |
17.40 |
2783 |
580 |
94.48 |
697.6 |
1353 |
5.969 |
83.4 |
534.6 |
0.931 |
12.80 |
3607 |
600 |
123.4 |
649.4 |
1176 |
6.953 |
75.7 |
495.3 |
1.06 |
8.38 |
5141 |
620 |
159 |
586.9 |
941 |
9.354 |
67.3 |
454.1 |
1.39 |
4.27 |
9092 |
640 |
202.7 |
481.5 |
560 |
25.94 |
55.3 |
414.9 |
3.46 |
0.81 |
39710 |
647 |
220.4 |
357.3 |
47 |
3905 |
46.9 |
1323 |
138 |
0.00 |
7735000 |
|
|
|
|
|
|
|
|
|
|
381 drop pressure phase-Two 8.11
382 CHAPTER 11 Boiling and condensation
Table 11.4 Properties of saturated steam (Meyer et al., 1993; Harvey et al., 2000).
T |
p |
ρg |
cpg |
g × 106 |
kg × 103 |
Prg |
βg × 106 |
(K) |
(bar) |
(kg/m3) |
(kJ/kg K) |
(Ns/m2) |
(W/m K) |
|
(K–1) |
273.16 |
0.006177 |
0.004855 |
1.884 |
9.216 |
17.07 |
1.02 |
3681 |
275 |
0.006985 |
0.005507 |
1.886 |
9.26 |
17.17 |
1.02 |
3657 |
280 |
0.009918 |
0.007681 |
1.891 |
9.382 |
17.44 |
1.02 |
3596 |
285 |
0.01389 |
0.01057 |
1.897 |
9.509 |
17.73 |
1.02 |
3538 |
290 |
0.0192 |
0.01436 |
1.902 |
9.641 |
18.03 |
1.02 |
3481 |
295 |
0.02621 |
0.01928 |
1.908 |
9.778 |
18.35 |
1.02 |
3428 |
300 |
0.03537 |
0.02559 |
1.914 |
9.92 |
18.67 |
1.02 |
3376 |
305 |
0.04719 |
0.0336 |
1.92 |
10.06 |
19.01 |
1.02 |
3328 |
310 |
0.06231 |
0.04366 |
1.927 |
10.21 |
19.37 |
1.02 |
3281 |
320 |
0.1055 |
0.07166 |
1.942 |
10.52 |
20.12 |
1.02 |
3195 |
340 |
0.2719 |
0.1744 |
1.979 |
11.16 |
21.78 |
1.01 |
3052 |
360 |
0.6219 |
0.3786 |
2.033 |
11.82 |
23.69 |
1.01 |
2948 |
373.15 |
1.014 |
0.5982 |
2.08 |
12.27 |
25.1 |
1.02 |
2902 |
400 |
2.458 |
1.369 |
2.218 |
13.19 |
28.35 |
1.03 |
2874 |
420 |
4.373 |
2.352 |
2.367 |
13.88 |
31.13 |
1.06 |
2914 |
440 |
7.337 |
3.833 |
2.56 |
14.57 |
34.23 |
1.09 |
3014 |
460 |
11.71 |
5.983 |
2.801 |
15.26 |
37.66 |
1.13 |
3181 |
480 |
17.9 |
9.014 |
3.098 |
15.95 |
41.45 |
1.19 |
3428 |
500 |
26.39 |
13.2 |
3.463 |
16.65 |
45.67 |
1.26 |
3778 |
520 |
37.69 |
18.9 |
3.926 |
17.38 |
50.44 |
1.35 |
4274 |
540 |
52.37 |
26.63 |
4.54 |
18.15 |
56.1 |
1.47 |
4994 |
560 |
71.06 |
37.15 |
5.41 |
19.01 |
63.34 |
1.62 |
6091 |
580 |
94.48 |
51.74 |
6.76 |
20.02 |
73.72 |
1.84 |
7904 |
600 |
123.4 |
72.84 |
9.181 |
21.35 |
91.05 |
2.15 |
1135 |
620 |
159 |
106.3 |
14.94 |
23.37 |
126.7 |
2.76 |
2000 |
640 |
202.7 |
177.1 |
52.59 |
27.94 |
250 |
5.88 |
7995 |
647 |
220.4 |
286.5 |
53340 |
39.72 |
1573 |
135 |
9274000 |
|
|
|
|
|
|
|
|
model overestimating the acceleration pressure drop, as it neglects the slip between the phases and assumes equal acceleration for both liquid and vapor phases. It may be noted that the homogenous model gives accurate results for dispersed flow – bubbly flow with small bubbles dispersed in liquid and spray flow or drop flow with small droplets dispersed in vapor, as the vapor velocity is approximately equal to the liquid velocity. The homogenous model also gives accurate results for higher pressures where vapor properties approach liquid properties. For accurate estimates of pressure drop over a wide range of operating conditions and vapor qualities, the correlations based on the separated flow model are recommended (Collier and Thome, 1994; Stephan, 1992; Jayaramu et al., 2019).
11.8 Two-phase pressure drop 383
In summary, the preceding sections gave a step by step approach to quickly estimate heat transfer and pressure drop in boiling and condensation based on the empirical correlations reported by many researchers in the past. These should be good enough for commonly encountered engineering applications, though the reader is encouraged to take the cognizance of the assumptions under which these are valid. Should situations arise where these assumptions do not hold, more advanced correlations need to be used. If these too do not work for a problem, numerical modeling or in-house experimentation may be required to determine heat transfer coefficient and pressure drop.
Example 11.3: Water flows upward in a vertical uniformly heated tube of length 0.9 m. The mass flux and the inlet subcooling are 500 kg/m2s and 0 K, respectively. The inner diameter of the tube is 25 mm. The system pressure is 10 bar. Determine the critical heat flux (given F1 = 0.478, F2 = 0.662, F3 = 0.4, F4 = 0.0166 for p = 10 bar) and the corresponding local quality at the exit.
Solution:
Using Eq. (11.44),
|
|
D G ( hsub ) |
|
|
|
|
|
|
|
|
|
||||||||||
|
|
A' + |
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
||||||
qcrit = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
C' + z |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
= |
|
A' |
(since ( |
Tsub ) = 0) |
|
|
|
|
|
||||||||||||
C ' + z |
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
A' = |
2.317 |
D G hfg |
|
|
|
|
|
F1 |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
4 |
|
|
|
|
|
|
|
|
1/ 2 |
G |
|
|||||||||
|
|
|
|
|
1 + 0.0143 F2 D |
|
|
|
|||||||||||||
= |
|
|
|
0.025 500 2014.9 103 0.478 |
|
||||||||||||||||
2.317 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
(1 + |
0.0143 0.662 |
|
|
|
|
|||||||||||||||
|
|
|
4 |
0.0251/ 2 500) |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
= 3.988 106 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
n = 2 − 0.00725 p |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
= 2 − 0.00725 × 10 (since p = 10 bar) |
|
||||||||||||||||||||
= 1.9275 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
C' = |
|
0.077 F3 D G |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
G n |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
1 + 0.347 F4 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
1356 |
|
|
|
|
|
|
|
|||||||
= |
|
0.077 0.4 |
0.025 500 |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
500 1.9275 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
1 + 0.347 0.0166 |
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1356 |
|
|
|
|
|
||
= 0.384
384 CHAPTER 11 Boiling and condensation
Therefore,
= 3.988 ×106 qcrit 0.384 + 0.9
= 3.104 ×106 W/m2
From energy balance, the exit quality can be determined as follows:
G |
π D2 |
(hf + xe hfg − hf ) = qcrit π D z |
|||
|
4 |
|
|
4 qcrit z |
|
|
|
xe = |
|
||
|
|
|
G Dh |
||
|
|
|
|
fg |
|
|
|
= |
|
4 × 3.104 × 106 × 0.9 |
|
|
|
500 × 0.025 × 2014.9 × 103 |
|||
|
|
xe = |
0.4436 |
|
|
It is instructive to note the order of magnitude of critical heat flux for flow boiling. As a rule of thumb it goes as MW/m2, for water.
Example 11.4: For the problem in Example 11.3, with heat flux equal to the critical heat flux, determine the local heat transfer coefficients and the corresponding wall temperatures at two locations: z = 0.4 m and 0.8 m from the inlet. Use the Chen correlation.
Solution:
Saturated water properties corresponding to 10 bar are:
k f = 0.677 W/mK, Prf = 0.967, σ = 0.04226 N/m, cpf = 4.403 × 103 J/kg K,
ρg = 5.16 kg/m3 , ρf = 887 kg/m3 , hfg = 2014.9 × 103 J/kg, µf = 148.5 × 10−6 Ns/m2 , µg = 15.02 × 10−6 Ns/m2
Using Eq. (11.34),
1 |
|
|
|
|
|
x 0.9 |
|
|
|
ρ f |
|
0.5 |
µg |
|
0.1 |
|
|
|
|
|
||||||||
|
|
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Xtt |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
1 |
− x |
|
|
|
ρg |
|
|
µ f |
|
|
|
|
|
|
|
||||||||||
|
q = q |
|
= 3.104 106 W/m2 |
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
crit |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
At z = 0.4 m, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
x = |
4 q z |
= |
|
|
|
4 × 3.104 × 106 × 0.4 |
|
= 0.197 |
|||||||||||||||||||
|
|
|
500 × 0.025 × 2014.9 × 103 |
|||||||||||||||||||||||||
|
|
|
|
|
G Dhfg |
|
|
|
|
|||||||||||||||||||
1 |
|
|
|
|
0.197 |
|
|
0.9 |
|
887 0.5 |
|
15.02 × 10−6 0.1 |
||||||||||||||||
|
|
|
|
= |
|
|
|
|
|
|
|
|
× |
|
|
|
× |
× 10 |
−6 |
|
||||||||
|
Xtt |
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
1 |
− 0.197 |
|
|
|
|
5.16 |
|
148.5 |
|
|
|||||||||||||||
|
|
|
= 2.944 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Re |
|
= |
G D(1 − x) |
= |
500 × 0.025 × (1 − 0.197) |
= 67612 |
||||||||||||||||||||||
|
|
|
||||||||||||||||||||||||||
|
|
f |
|
|
|
|
|
|
f |
|
|
|
|
|
|
|
148.5 × 10−6 |
|
|
|
|
|
||||||
11.8 Two-phase pressure drop 385
Using Eq. (11.35), and since |
1 |
> 0.1, |
|||
|
|||||
|
|
|
|
Xtt |
|
|
1 |
0.736 |
|
|
|
F = 2.35 |
|
|
+ 0.213 |
|
|
|
|
|
|||
|
Xtt |
|
|
|
|
=2.35(2.944 + 0.213)0.736
=5.477
Using Eq. (11.32),
0.8 |
0.4 |
kf |
|
||
hc = 0.023Ref |
Prf |
|
|
|
F |
|
|||||
|
|
|
D |
|
|
=0.023 (67612)0.8 (0.967)0.4 0.677 5.477
0.025
=24555W/m2 K
=24.55kW/m2 K
Using Eq. (11.33),
ReTP = Ref × F1.25
= 67612 × 5.4771.25
= 5.665 ×105
Using Eq. (11.41),
S = |
1 |
|
|
1+ 2.53 ×10−6 × ReTP1.17 |
|
= |
0.0683 |
|
Using Eq. (11.40),
hNcB = 0.00122 |
|
|
k0.79f c0.45pf ρ0.49f |
|
(Tw − Tsat ) |
0.24 |
(psat (Tw )− psat (Tsat )) |
0.75 |
(S) |
|
||||||
|
|
|
|
|
|
|
|
|
|
|||||||
σ |
0.5 |
0.29 |
0.24 |
0.24 |
|
|
|
|||||||||
|
|
|
µ f |
hfg |
ρg |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.734 × 43.62 × 27.83 |
|
|
0.24 |
(psat (Tw )− psat (Tsat )) |
|
|||||||
= 0.00122 |
|
|
|
|
|
|
|
|
|
|
(Tw − Tsat ) |
|
0.75 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
0.2055 × 0.0776 × 32.58 × 1.4826 |
|
|
|
|
|
|||||||||
×0.0683
=0.0963 × (Tw − Tsat )0.24 (psat (Tw )− psat (Tsat ))0.75
Using Eq. (11.30),
hTP = hc + hNcB
q = hTP (Tw − Tsat )
3.104 × 106 = 24555 + 0.0963 × (Tw − Tsat )0.24 (psat (Tw )− psat (Tsat ))0.75 (Tw − Tsat ) Since q = qcrit , (Tw − Tsat ) must be high.
We now face a special difficulty here. The heat flux on the left-hand side of the equation is known. The right-hand side contains Tw, which we are seeking. Unfortunately we cannot solve for Tw, as the expression on the right-hand side has the saturation pressure
psat (Tw), which is as yet unknown in view of the fact that Tw is unknown. Hence we have no option but to go for an iterative solution with an assumed value of Tw and iterating it
to convergence.
Let (Tw − Tsat ) = 50 C
386 CHAPTER 11 Boiling and condensation
Tw = Tsat (10 bar) + 50 C
=180 + 50 C
=230 C
psat (Tw ) = 27.98 bar = 27.98 × 105 N/m2
psat (Tw ) − psat (Tsat ) = (27.98 − 10) × 105 = 17.98 × 105 N/m2
q = hTP (Tw − Tsat )
=24555 + 0.0963 × (Tw − Tsat )0.24 (psat (Tw ) − psat (Tsat ))0.75 (Tw − Tsat )
=24555 + 0.0963 × 500.24 × (17.98 × 105 )0.75 × 50 =1.832 × 106 < 3.104 × 106 W/m2
Let (Tw − Tsat ) = 65 C
Tw = Tsat (10 bar) + 65 C =180 + 65 C
= 245 C
psat (Tw ) = psat (245 C) = 36.62 × 105 N/m2
psat (Tw ) − psat (Tsat ) = (36.62 − 10) × 105 = 26.62 × 105 N/m2
q = hTP (Tw − Tsat )
q= 24555 + 0.0963 × 650.24 × (26.62 × 105 )0.75 × 65 = 2.719 × 106 < 3.104 × 106 W/m2
Let (Tw − Tsat ) = 70 C
Tw = 180 + 70 C = 250 C psat (Tw ) = 39.776 × 105 N/m2
psat (Tw ) − psat (Tsat ) = 29.776 × 105 N/m2
q = hTP (Tw − Tsat )
= 3.058 × 106 ≈ 3.104 × 106 W/m2
Therefore, at z = 0.4 m,
|
hTP = |
3.058 × 106 |
= 43685 W/m2K |
|
|
|
|
|
|||||||||||||
|
Tw = |
|
|
|
70 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
250 C |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
At z = |
0.8 m, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
x = |
|
|
4 q z |
= |
|
|
|
4 × 3.104 × 106 × 0.8 |
= 0.394 |
|||||||||||
|
|
|
500 × 0.025 × 2014.9 × 103 |
||||||||||||||||||
|
|
|
|
|
G Dhfg |
|
|
|
|
|
|||||||||||
1 |
|
|
|
0.394 |
|
|
0.9 |
|
887 0.5 |
15.02 × |
10−6 |
0.1 |
|||||||||
|
|
|
= |
|
|
|
|
|
|
|
× |
|
× |
|
|
−6 |
|
= 7.07 |
|||
|
Xtt |
|
− 0.394 |
|
10 |
||||||||||||||||
|
|
1 |
|
|
|
|
|
5.16 |
148.5 × |
|
|
|
|||||||||
Re |
|
= |
|
G D(1 − x) |
= |
500 × 0.025 × (1 − 0.394) |
= 51041 |
||||||||||||||
|
|
|
|
||||||||||||||||||
|
|
f |
|
|
|
|
f |
|
|
|
|
|
|
148.5 × 10−6 |
|
|
|
|
|
||
11.8 Two-phase pressure drop 387
|
|
1 |
0.736 |
F = 2.35 |
|
|
+ 0.213 |
|
|||
|
Xtt |
|
|
=2.35 × (7.07 + 0.213)0.736
=10.13
hc = 0.023 Re0f .8 Prf0.4 kDf F
=36349W/m2K
ReTP = Ref × F1.25
=9.244 × 105
S = |
|
|
|
1 |
|
= 0.0398 |
1 |
+ 2.53 |
×10 |
−6 1.17 |
|||
|
ReTP |
|||||
Using S = 0.0398 in the expression for hNcB,
hNcB = 0.0561× (Tw − Tsat )0.24 (psat (Tw )− psat (Tsat ))0.75
hTP = hc + hNcB
q = hTP (Tw − Tsat ) = 36349 + 0.0561× (Tw − Tsat )0.24 (psat (Tw )− psat (Tsat ))0.75 (Tw − Tsat )
Let(Tw − Tsat ) = 65 C
Tw = Tsat + 65 C = 180 + 65 C = 245 C
psat (Tw )− psat (Tsat ) = 26.62 × 105 N/m2
Substituting the above values in the expression for q,
q = 3.017 ×106 ≈ 3.104 ×106 W/m2
Therefore, at z = 0.8 m,
hTP = 3.017 ×106 = 46415W/m2K 65
Tw = 245 C
Note that this sample problem does not ask for the calculation of the heat transfer coefficient at z = 0.9 m (tube exit). It is intentional, as the heat flux is equal to the critical heat flux, which indicates dryout at the exit. The Bowring CHF correlation presented in the chapter (or for that matter any other CHF correlation) estimates the CHF that corresponds to the onset of dryout (or DNB) at the tube exit. The Chen correlation should not be used to calculate the two-phase heat transfer coefficient for dryout and post-dryout (spray flow or drop flow) heat transfer, though the correlation gives some values of the heat transfer coefficient. The heat transfer coefficient calculated at the tube exit using the Chen correlation can at best be taken as the value prior to the onset of dryout at the tube exit.
388CHAPTER 11 Boiling and condensation
Example 11.5: For the problem in Example 11.3, with heat flux equal to the critical heat flux, evaluate the pressure drop across the tube using the homogeneous model.
Solution:
Fluid properties corresponding to 10 bar are:
f = 148.5 × 10−6 Ns/m2 , g = 15.02 × 10−6 Ns/m2 Specific volumes v f = 1.127 × 10−3 m3 /kg, vg = 0.1937m3 /kg
v fg = vg − v f = 0.1926 m3 /kg
The tube exit quality, xe = 0.4436
Fanning friction factor at the tube inlet:
Re |
|
= |
GD = |
500 × 0.025 |
= 84175 |
|
f |
|
f |
148.5 ×10−6 |
|
Using the Blasius equation,
f |
|
= |
0.079 = 0.00464 |
|
f |
|
Re1/4f |
Fanning friction factor at the tube exit:
Two-phase viscosity can be calculated from the relation
1 |
|
= |
|
|
x |
+ 1− x |
|
||
|
|
|
|
|
|
|
g |
f |
|
|
|
|
|
||||||
|
|
|
|
= |
0.4436 |
+ (1− 0.4436) |
|||
|
|
|
|
|
|
15.02 × 10−6 |
148.5 × 10−6 |
||
|
|
|
= 30.4 × 10−6 Ns/m2 |
||||||
Re |
= |
|
|
GD = 500 × 0.025 = 4.16 × 105 |
|||||
TP |
|
|
|
|
30.4 × 10−6 |
||||
Using the Blasius equation,
= 0.079 =
fTP ReTP1/4 0.0031
The average value of fTP over the tube length is given by
fTP, avg = |
ff + fTP |
= 0.00387 |
|
2 |
|||
|
|
The total pressure drop can be determined using Eq. (11.91). The frictional pressure drop from the homogeneous model is given by
p |
|
= 2 fTP,avg G2 vf |
z |
|
1 + vfg xe |
|
|
|
|
||||||||
|
f |
|
|
D |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
vf 2 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 0.00387 5002 1.127 |
10−3 0.9 |
0.1926 0.4436 |
||||||||||||
|
|
= |
|
|
|
0.025 |
|
|
|
|
1 |
+ |
1.127 10 |
−3 |
|
||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
||||
= 3054 N/m2
11.8 Two-phase pressure drop 389
The acceleration pressure drop from the homogeneous model is given by
pa = G2 v fg xe
=5002 × 0.1926 × 0.4436
=21359 N/m2
The gravitational head pressure drop from the homogeneous model is given by
pg = |
g sin |
(θ ) z |
|
+ |
vfg |
|
|
|
|
ln 1 |
|
xe |
|||
v fg |
xe |
v f |
|||||
|
|
|
|
||||
|
|
|
|
|
|
|
|
9.81 1 0.9 |
|
|
0.1926 |
|
|
|
|
|
= |
|
|
ln 1 |
+ |
|
|
|
0.4436 |
|
0.1926 |
0.4436 |
|
-3 |
||||||
|
|
1.127 10 |
|
|
|
|
|||
= 449 N/m2
The total pressure drop
ptotal = pf + pa + pg
=3.054 + 21.359 + 0.449 kN/m2
=24.862 kN/m2
Example 11.6: Saturated water at 0.123 bar with a mass flux of 400 kg/m2s enters a vertical tube and flows upward. There is saturated steam at 1.013 bar condensing on the tube outer surface. The length and diameter of the tube are 0.6 m and 0.025 m, respectively. Assume laminar film condensation. Determine
(a)The average heat transfer coefficient on the condensation side.
(b)The average heat transfer coefficient on the boiling side, which can be assumed to be approximately equal to the mean of the local heat transfer coefficients at the inlet and outlet of the tube. You may use the Chen correlation.
(c)The average tube wall temperature.
(d)The rate of condensation (total), kg/s.
(e)The exit quality on the boiling side.
Solution:
The average heat transfer coefficient relationship for laminar film condensation on a vertical plate, given by Eq. (11.63), can also be used for condensation on a vertical tube.
|
|
|
|
ρ f (ρ f − ρg )g hfg |
k3f |
1/ 4 |
|
|
|
||||
h f |
= 0.943 |
|
|
|
|
|
µ f z (Tsat − Tw ) |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Though all liquid properties in the above expression correspond to the film tem-
perature (Tw + Tsat)/2, and hfg corresponds to Tsat, let all the properties be evaluated at Tsat since Tw is unknown.
390 CHAPTER 11 Boiling and condensation
|
|
= 0.943 |
|
958 × (958 − 0.598) × 9.81 × 2257 × 103 × 0.6813 |
1/4 |
||
h f |
|
283.1 × 10−6 × 0.6 × (100 − Tw ) |
|
||||
|
|
|
|
|
|||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
||
|
|
= 0.943 × 13939.67 × |
|
|
|
||
|
|
(100 − Tw )1/4 |
|
||||
The two-phase heat transfer coefficient for saturated boiling is obtained using Eq. (11.30),
hTP = hc + hNcB
Using Eq. (11.34) and the saturated fluid properties at 0.123 bar,
1 |
ρ f 0.5 µg |
|
0.1 |
x 0.9 |
|
|
|
|
||||||||||
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Xtt |
|
|
|
|
|
|
|
|
|
|||||||||
|
ρg µ f |
|
|
|
1− x |
|
|
|
|
|
||||||||
|
|
988 |
0.5 |
|
|
10.1 10−6 |
0.1 |
|
x 0.9 |
|||||||||
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
547.8 10 |
−6 |
|
||||||||||
|
|
0.083 |
|
|
|
|
|
|
|
1− x |
||||||||
|
= 73.17 |
|
x |
0.9 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
− x |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||||
From Eq. (11.35),
|
|
|
1 |
|
|
|
|
||
F = 1 for |
|
|
|
≤ 0.1 |
|
|
|
||
|
|
|
|
|
|||||
|
|
Xtt |
|
|
|
|
|||
|
1 |
|
|
0.736 |
|
1 |
|
||
= 2.35 |
|
|
+ 0.213 |
for |
|
> 0.1 |
|||
|
|
||||||||
|
Xtt |
|
|
|
|
Xtt |
|
||
From Eq. (11.32),
h = 0.023 |
G (1 − x) D 0.8 |
cp 0.4 |
kf |
(F) |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
c |
|
f |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
k f |
D |
|
|
|||||||
= 0.023 |
400 |
0.025 |
0.8 |
(1 − x) |
0.8 |
(3.56) |
0.4 |
|||||||||
|
|
−6 |
|
|
|
|
|
|||||||||
|
547.8 10 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
= 2521.53 (1 − x)0.8 F
From Eq. (11.40),
0.643 (F)0.025
hNcB = 0.00122 |
|
|
k0.79f |
c0.45pf ρ0.49f |
|
(Tw |
− Tsat ) |
0.24 |
(psat (Tw ) − psat |
(Tsat )) |
0.75 |
(S) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||
σ |
0.5 |
0.29 |
0.24 |
0.24 |
|
|
|
|
||||||||||||
|
|
|
µ f |
|
hfg |
ρg |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= 0.00122 |
|
|
|
0.7055 42.16 29.34 |
|
|
|
T 0.24 |
|
p0.75 |
S |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
sup |
|
sup |
|
|
|
||
|
|
0.2606 0.1133 33.92 0.55 |
|
|
|
|
|
|
|
|
||||||||||
From Eq. (11.33),
ReTP = Ref × F1.25
=G(1 − x) D × F1.25
f
=18254.8 × (1 − x) × F1.25
From Eq. (11.41),
S = |
1 |
|
|
|
|
1 + 2.53 × 10 |
−6 Re1.17 |
|
|
|
TP |
