11.4 Flow boiling 371
Table 11.2 Functions F1, F2, F3, and F4 for the Bowring critical heat flux correlation, for only water (Bowring, 1972).
Pressure (bar) |
F1 |
F2 |
F3 |
F4 |
1 |
0.478 |
1.782 |
0.4 |
0.0004 |
5 |
0.478 |
1.019 |
0.4 |
0.0053 |
10 |
0.478 |
0.662 |
0.4 |
0.0166 |
15 |
0.478 |
0.514 |
0.4 |
0.0324 |
20 |
0.478 |
0.441 |
0.4 |
0.0521 |
25 |
0.48 |
0.403 |
0.401 |
0.0753 |
30 |
0.488 |
0.39 |
0.405 |
0.1029 |
35 |
0.519 |
0.406 |
0.422 |
0.138 |
40 |
0.59 |
0.462 |
0.462 |
0.1885 |
45 |
0.707 |
0.564 |
0.538 |
0.2663 |
50 |
0.848 |
0.698 |
0.647 |
0.3812 |
60 |
1.403 |
0.934 |
0.89 |
0.7084 |
68.9 |
1.000 |
1.000 |
1.000 |
1.000 |
70 |
0.984 |
0.995 |
1.003 |
1.03 |
80 |
0.853 |
0.948 |
1.033 |
1.322 |
90 |
0.743 |
0.903 |
1.06 |
1.647 |
100 |
0.651 |
0.859 |
1085 |
2.005 |
110 |
0.572 |
0.816 |
1.108 |
2.396 |
120 |
0.504 |
0.775 |
1.129 |
2.819 |
130 |
0.446 |
0.736 |
1.149 |
3.274 |
140 |
0.395 |
0.698 |
1.168 |
3.76 |
150 |
0.35 |
0.662 |
1.186 |
4.227 |
160 |
0.311 |
0.628 |
1.203 |
4.825 |
170 |
0.277 |
0.595 |
1.219 |
5.404 |
180 |
0.247 |
0.564 |
1.234 |
6.013 |
190 |
0.22 |
0.534 |
1.249 |
6.651 |
200 |
0.197 |
0.506 |
1.263 |
7.32 |
|
|
|
|
|
11.4.4 A brief overview of flow boiling in micro-channels
Rapid miniaturization of electronic devices poses challenges to the thermal management of electronics. Data center cooling and cooling of rocket motors and laser diodes involve high heat fluxes. As micro-channels have high heat transfer surface area to fluid flow volume ratio, and flow boiling heat transfer coefficients are high, heat sinks based on the flow boiling in micro-channels is one of the promising techniques to efficiently dissipate high heat fluxes. Evaporating vapor bubbles in microchannels are confined by the channel cross-sectional dimensions (Gedupudi et al., 2011), and hence the flow patterns during flow boiling in micro-channels can be different from those in the conventional size tubes. The confined vapour bubble growth
372 CHAPTER 11 Boiling and condensation
leads to pressure fluctuations, which can result in transient flow reversal if upstream compressibility is present (Gedupudi et al., 2011; Jain et al., 2019). Another characteristic is the possibility of occurrence of intermittent dryout caused by the cyclic passage of liquid slugs and long confined bubbles (Jain et al., 2019). Enhancing the two-phase heat transfer coefficient and the critical heat flux and reducing the flow instabilities and pressure drop are the research challenges in this area (Kandlikar et al., 2006; Karayiannis and Mahmoud, 2017).
11.5 Condensation
The condensation process is the reverse of the boiling process and involves the change of a vapor phase to a liquid phase. Just as liquid superheat is required to induce the nucleation of bubbles in boiling, vapor subcooling is required to induce the nucleation of droplets in condensation. If the condensate forms a continuous film, then it is called film-wise condensation (or simply, film condensation), which occurs on wetted surfaces. In drop-wise condensation (or simply, drop condensation), the vapor condenses into small liquid droplets of different sizes that coalesce and fall down the cooled surface, and it occurs on non-wetted surfaces. Drop condensation, as seen from the above definition, offers less thermal resistance, and this results in heat transfer coefficients as much as 5 to 10 times the values in film condensation. Though drop condensation would be preferred to film condensation, it is difficult to achieve or maintain drop condensation. Much of the recent research focuses on promoting drop condensation using microstructured surfaces, which is beyond the scope of this textbook. Additionally, there is a lack of reliable theories of drop condensation.
11.6 Film condensation on a vertical plate
Nusselt (1916) made the following simplifying assumptions to solve the film condensation problem analytically, which is often considered a classic treatise on the subject of condensation:
1.The condensate flow in the film is laminar.
2.The plate temperature, Tw, is uniform and less than the vapor saturation temperature, Tsat.
3.The fluid properties are constant.
4.Negligible shear stress at the liquid-vapor interface.
5.Momentum changes in the film are negligible.
6.The heat transfer across the film is by conduction only, and the condensate temperature profile is linear.
Fig. 11.7 shows the typical growth of condensate film as well as velocity and temperature profiles. The analysis is made per unit width of the plate in the direction
11.6 Film condensation on a vertical plate 373
FIGURE 11.7
Film condensation on a vertical plate.
perpendicular to the plane of the paper. Consider a force balance on a differential element of the film at any axial location, x, shown in Fig. 11.7.
ρ f gdxdy + τ (y + dy)dx + p(x)dy = τ (y)dx + p(x + dx)dy |
(11.48) |
||||||||
|
∂τ |
|
∂ p |
|
|
|
(11.49) |
||
|
|
= −ρ f g + |
|
|
|
|
|
||
|
∂y |
∂x |
|
|
|
|
|
||
Using the boundary layer type approximation |
∂ p |
|
, |
|
|||||
|
= 0 |
|
|||||||
|
|
|
|
|
|
∂y |
|
|
|
|
|
dp |
|
|
|
|
|
|
(11.50) |
|
|
dx = ρg g |
|
|
|
||||
|
|
|
|
|
|
||||
Therefore, Eq. (11.49) changes to |
|
|
|
|
|
|
|
||
∂τ |
− ρg )g |
|
|
|
|
||||
∂y = − (ρ f |
|
|
|
(11.51) |
|||||
|
|
|
|
|
|
|
|
|
|
According to Newton’s law of viscosity, |
|
|
|
|
|
|
|
||
|
|
τ = µ f |
∂u |
|
|
|
(11.52) |
||
|
|
∂y |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
Substituting Eq. (11.52) in Eq. (11.51),
∂2 u2 |
= |
−g (ρ f − ρg ) |
(11.53) |
∂y |
|
µ f |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
374 |
CHAPTER 11 Boiling and condensation |
|
|
|
|
|
|
|
|
|||||||||||
|
Integrating Eq. (11.53) twice and applying the boundary conditions u = 0 at y = 0, |
|||||||||||||||||||
|
and ∂u = 0 at y = δ, the velocity profile becomes |
|
|
|
|
|
|
|
|
|||||||||||
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u(y) = |
(ρ f |
− ρg )gδ 2 y |
− |
|
y2 |
|
(11.54) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
µ f |
|
|
2δ |
2 |
|||||||||
|
|
|
|
|
|
|
|
|
δ |
|
|
|
|
|
||||||
|
The condensate mass flow rate at any location x is given by |
|
||||||||||||||||||
|
|
|
δ (x) |
ρ f u |
(y)dy = |
ρ f (ρ f |
|
− ρg )gδ 3 |
(11.55) |
|||||||||||
|
m(x) = ∫0 |
|
|
|
|
|
3µ f |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
dm |
= |
ρ f |
(ρ f − ρg )gδ |
2 |
d δ |
|
|
|
(11.56) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
µ f |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
From the heat transferred by conduction, |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
dm = |
kf (Tsat − Tw ) |
dx |
|
|
|
|
|
(11.57) |
|||||||||
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
δ hfg |
|
|
|
|
|
|
|
|
|
|
|
|
From Eq. (11.56) and Eq. (11.57), |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
3 |
|
|
|
|
|
k f (Tsat − Tw )µ f |
|
|
|
|
|
(11.58) |
||||||
|
|
δ dδ = ρ f (ρ f − ρg )ghfg |
|
dx |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
||||||||||||||
On integrating Eq. (11.58) with the boundary condition δ = 0 at x = 0,
|
4k f (Tsat − Tw )µ f |
1/4 |
|
δ (x) = |
|
x |
|
ρ f (ρ f − ρg )ghfg |
|||
|
|
To obtain the local heat transfer coefficient,
hx dx (Tsat − Tw ) = |
k f (Tsat − Tw ) |
|||
|
δ (x) |
|||
|
|
|
||
hx = |
k f |
|
||
δ (x) |
||||
|
||||
Substituting Eq. (11.61) in Eq. (11.59),
|
|
ρ f (ρ f − ρg )ghfg k3f |
1/4 |
|
h |
= |
|
|
|
4µ f (Tsat − Tw )x |
||||
x |
|
|
||
|
||||
|
|
|
|
(11.59)
(11.60)
(11.61)
(11.62)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
11.6 Film condensation on a vertical plate |
375 |
||||
The average heat transfer coefficient is |
|
|
|
|
||||||
|
|
|
1 L |
|
|
ρ f (ρ f − ρg )ghfg k3f |
1/4 |
|
|
|
|
|
|
|
|
|
|||||
hL = |
|
|
h dx = 0.943 |
|
|
|
|
|
||
L ∫0 |
µ f (Tsat − Tw )L |
(11.63) |
|
|||||||
|
|
|
x |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
All liquid properties in Eq. (11.62) and Eq. (11.63) are best evaluated at the mean
film temperature (Tsat + Tw ) / 2, and hfg and ρg at Tsat. Although Eq. (11.63) has been derived for a vertical flat plate, the expression can also be used for film condensation
on the inside or outside surfaces of a vertical tube provided the radius is large compared with the film thickness at the bottom of the tube.
For an inclined flat plate that makes an angle θ with the horizontal, it can be readily seen that Eq. (11.63) modifies to
|
|
|
ρ f (ρ f − ρg )gsin (θ )hfg k3f |
1/4 |
|
|
|
|
|
||||
hL = 0.943 |
|
|
|
(11.64) |
||
µ f (Tsat − Tw )L |
||||||
|
|
|
|
|
||
|
|
|
|
|
|
|
The average heat transfer coefficient can also be expressed as a function of the condensate flow rate at the bottom of the plate.
From energy balance,
|
|
|
. |
|
|
||
|
h |
L L (Tsat − Tw ) = mL hfg |
|||||
|
|
|
. |
|
|
(11.65) |
|
|
L (Tsat − Tw ) = |
mL |
hfg |
||||
|
|
|
|
|
|
||
|
hL |
||||||
|
|
|
|||||
Substituting Eq. (11.65) in Eq. (11.63),
|
|
|
ρ f (ρ f − ρg )gk3f |
1/3 |
|
|
hL = 0.925 |
|
|
(11.66) |
|||
. |
||||||
|
|
|
µ f mL |
|
|
|
|
|
|
|
|
||
The film Reynolds number at the bottom of the plate is
|
|
|
|
|
|
|
|
. |
|
|
ρ f (ρ f |
− ρg )gδ 3 |
||
|
ρ f u(L) Dh |
|
ρ f u(L)4δ |
|
4 mL |
|
4 |
|||||||
Reδ = |
|
|
|
= |
|
|
|
= |
|
= |
|
|
|
(11.67) |
|
|
µ f |
|
|
µ f |
µ f |
|
3µ2f |
|
|||||
Substituting Eq. (11.67) in Eq. (11.66) and assuming (ρ f − ρg ) ≈ ρ f ,
|
L |
µ |
2 |
1/3 |
|
||
h |
|
||||||
|
|
|
f |
|
= 1.47Reδ−1/3 |
(11.68) |
|
|
|
2 |
|||||
k f |
|
ρ f |
g |
|
|
||
|
|
|
|
|
|
|
|
Experimental data suggests that the wave-free laminar regime exists at the bottom of the plate for Reδ ≤ 30.
376 CHAPTER 11 Boiling and condensation
It may be noted that Reδ in Eq. (11.68) is unknown. Reδ |
can be expressed as |
||||||||||
|
|
|
|
|
|
µ2f |
1/3 |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
4hL |
|
|
|
|
|
|
|||
|
2 |
|
|
|
|||||||
Reδ = |
|
|
|
ρ f g |
φ |
(11.69) |
|||||
|
|
|
k f |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
||
where parameter φ is |
|
|
|
|
|
|
|
|
|
|
|
φ = |
k f |
L (Tsat − Tw ) |
|
|
|
||||||
µ f |
|
|
|
µ2f |
|
1/3 |
|
|
(11.70) |
||
|
hfg |
|
|
|
|
|
|
|
|||
|
2 |
|
|
|
|
||||||
|
|
|
|
|
|
ρ f g |
|
|
|
||
Substituting Eq. (11.69) in Eq. (11.68),
|
L |
µ |
2 |
1/3 |
|
|
||
h |
φ−1/4 |
|
||||||
|
|
|
f |
|
= 0.943 |
(11.71) |
||
k f |
2 |
|||||||
|
ρ f |
g |
|
|
|
|||
|
|
|
|
|
|
|
|
|
The following equations are recommended for different flow regimes at the bottom of the plate.
Wave-free laminar:
|
L |
µ |
2f |
1/3 |
|
|
||
h |
|
|
||||||
|
|
|
|
|
|
= 0.943 φ−1/ 4 |
, φ ≤ 15.8 |
(11.72) |
k f |
2 |
|
||||||
|
ρ f |
g |
|
|
|
|||
|
|
|
|
|
|
|
|
|
Laminar-wavy (Kutateladze, 1963):
|
L |
µ |
2f |
1/3 |
|
1 |
|
|
|
|
|
|
h |
= |
(0.68 φ + 0.89) |
0.82 |
, 15.8 |
≤ φ ≤ 2530 |
|
||||||
|
|
|
|
|
|
φ |
|
(11.73) |
||||
k f |
2 |
|
|
|||||||||
|
ρ f |
g |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Turbulent (Labuntsov, 1957):
|
L |
µ2f |
1/3 |
|
1 |
|
|
|
|
|
|
|
|
||
h |
|
|
1/ 2 |
|
|
4 /3 |
|
|
|
||||||
|
|
|
|
|
|
= |
|
|
φ − 53)Prf |
+ 89 |
|
|
, |
φ ≥ 2530, Prf ≥ 1 |
(11.74) |
k |
|
ρ2 g |
|
|
|||||||||||
|
φ (0.024 |
|
|
||||||||||||
|
|
f |
|
f |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
As an improvement to the Nusselt theory, Rohsenow (1956) considered thermal advection effects and non-linear temperature profile across the film and showed that the latent heat of vaporization, hfg, should be changed to
h'fg = hfg + 0.68 cpf (Tsat − Tw ) |
(11.75) |
It may be noted that the presence of non-condensable gases in vapor degrades the heat transfer performance, as they increase the thermal resistance at the liquid-vapor interface. In view of this, the general practice in the design of a condenser is to vent the non-condensable gas.
11.7 Condensation on horizontal tubes 377
11.7 Condensation on horizontal tubes
In surface condensers, condensation occurs on a bank of horizontal tubes, as shown in Fig. 11.8.
Nusselt derived the following mean heat transfer coefficient relationship for laminar film condensation on the outside surface of a single horizontal tube:
|
|
ρ f (ρ f − ρg )g h'fg k3f |
1/ 4 |
|
|
h |
= 0.725 |
|
|
(11.76) |
|
Dµ f (Tsat − Tw ) |
|||||
D1 |
|
|
|||
|
|
||||
|
|
|
|
|
For a vertical column of N horizontal tubes, Jakob (1936) proposed the following mean heat transfer coefficient correlation:
hDN = N−1/4
hD1
where hDN is the mean heat transfer coefficient for N tubes.
The mean heat transfer coefficient for the Nth tube in the column,
(11.77)
hDN , is given by
hDN = N3/4 − (N − 1)3/4 |
|
||
hD1 |
(11.78) |
||
where hD1 is the mean heat transfer coefficient for the first tube. |
|
||
Kern (1958) recommended the following correlations for N ≥ 10 : |
|
||
|
|
DN = N−1/6 |
|
|
h |
(11.79) |
|
|
h |
|
|
|
|
D1 |
|
FIGURE 11.8
Film condensation on a bank of horizontal tubes.
378 CHAPTER 11 Boiling and condensation
hDN |
= N5/6 − (N − 1)5/6 |
(11.80) |
h |
|
|
D1 |
|
|
The above equations clearly show that the Nth tube will have a much lower heat transfer coefficient than the first, which is intuitively apparent.
Note that the Nusselt theory neglects the shear at the liquid-vapor interface. The interfacial shear becomes significant for condensation with flowing vapor and hence is an important factor that influences the heat transfer coefficient for condensation in tubes. Readers can refer to advanced texts (Collier and Thome, 1994; Stephan, 1992) for condensation in vertical and horizontal tubes.
11.8 Two-phase pressure drop
Pressure drop is an important consideration in the design of heat transfer equipment and more so for two-phase flows. Two-phase drop involves the pressure drop due to wall and interfacial shear stresses, the acceleration pressure drop for the case of evaporating flows, and the pressure drop due to gravitational head. The homogenous model and the separated flow model are usually used to estimate the two-phase pressure drop. The separated flow model, which is more complex, makes better prediction than the homogenous model but is beyond the scope of the present study. The interested reader can refer to advanced texts, such as Collier and Thome (1994) and Stephan (1992).
The homogenous model is based on the following assumptions:
(i)the vapor and liquid velocities are equal (ug = uf = u )
(ii)the thermodynamic equilibrium exists between the phases, and
(iii)an appropriately defined single-phase friction factor can be used for two-phase flow.
The total pressure gradient is the sum of the frictional, acceleration, and gravitational head pressure gradients,
dp |
dp |
dp |
dp |
||||
|
|
= |
|
+ |
|
+ |
|
dz |
dz f |
dz a |
dz g |
||||
The frictional pressure gradient is given by
dp |
|
2 fTP G2 |
|
|
2 fTP G2 |
(vf + xvfg ) |
|
|
v |
|
|||||
− |
|
= |
|
|
= |
|
|
D |
|
D |
|||||
dz f |
|
|
|
||||
(11.81)
(11.82)
where G, v, and fTP are the total mass flux (the product of density and velocity), specific volume, and the two-phase Fanning friction factor, respectively. Note that the Fanning friction factor is expressed as τ w / (ρu2 / 2), where τ w is the wall shear stress. The Darcy friction factor is four times the Fanning friction factor.
11.8 Two-phase pressure drop 379
The acceleration pressure gradient is given by
dp |
1 |
|
d |
. |
|
|
|
|
|
1 d |
|
2 d |
v |
|
|||||
− |
= |
|
|
|
|
(m |
u |
)= |
|
|
|
|
|
(AG u ) = G |
|
|
(11.83) |
||
|
|
|
|
|
|
|
|
|
|
|
|||||||||
dz a |
|
A dz |
|
|
|
|
|
|
A dz |
|
dz |
|
|||||||
|
− |
dp |
= G |
2 |
|
d |
(v f |
+ x v fg ) |
|
|
(11.84) |
||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
dz |
|
|
||||||||||||||
|
|
dz a |
|
|
|
|
|
|
|
|
|||||||||
Assuming vfg ≈ vg and neglecting the compressibility of the liquid phase,
dp |
|
2 |
|
dvg |
|
|
dx |
|||
− |
|
= G |
|
x |
|
|
+ v fg |
|
|
|
|
dz |
dz |
||||||||
dz a |
|
|
|
|
|
|
||||
|
|
= G |
2 |
|
dvg dp |
+ v fg |
||||
|
|
|
x |
|
|
|
|
|||
|
|
|
|
|
||||||
|
|
|
|
|
|
dp dz |
|
|||
The gravitational head pressure gradient is given by
dx
(11.85) dz
dp |
|
|
|
g sin (θ ) |
|
g sin (θ ) |
|
||||
= ρ g sin (θ ) = |
|
|
|||||||||
− |
|
|
|
|
|
= |
|
(11.86) |
|||
|
|
|
|
(v f + x v fg ) |
|||||||
dz g |
|
|
|
v |
|
|
|||||
where θ is the angle made by the inclined pipe with the horizontal. Substituting Eqs. (11.82), (11.85) and (11.86) in Eq. (11.81),
|
|
|
|
2 fTP G2 |
v |
|
+ xv |
|
+ G |
2 v |
|
dx |
+ |
g sin (θ ) |
|
|
|
|
|
|
|
f |
fg |
|
(vf + xvfg ) |
|
|||||||
dp |
|
D |
|
|
fg |
|
|
|
dz |
|
|
|||||
− |
|
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
(11.87) |
|
|
|
|
|
|
|
dvg |
|
|
|||||||
dz |
|
|
|
|
1 + G |
2 |
|
|
|
|||||||
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
dp |
|
|
|
||
The Fanning friction factor ( fTP ) can be evaluated using a mean two-phase viscosity, .
McAdams et al. (1942) suggested the following relationship for the mean twophase viscosity.
1 |
= |
x |
+ |
1− x |
(11.88) |
|
|
|
g |
f |
|||
|
|
|
|
|||
Cicchitti et al. (1960) suggested the relationship
|
= x g + (1− x ) f |
(11.89) |
For turbulent flow, the two-phase Fanning friction factor can be evaluated using the Blasius equation
fTP = |
0.079 |
= |
|
0.079 |
(11.90) |
||||
1/4 |
|
GD |
1/4 |
||||||
|
ReTP |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|||
380 CHAPTER 11 Boiling and condensation
Eq. (11.90) can be plugged into Eq. (11.87) to obtain the pressure gradient (dp/dz), and one needs to know the quality x and its axial gradient (dx/dz).
11.8.1 Total pressure drop
To obtain the total pressure drop, Eq. 11.87 needs to be integrated with respect to length. For simplicity, the following assumptions can be made.
1.The properties v f and vg remain constant over the length considered. The friction factor fTP can be an average value, which remains constant over the length considered.
2.The term G2 x (dvg /dp) in the denominator of Eq. (11.87) is much less than 1. The assumption of negligible compressibility of the gaseous phase is valid for many cases, especially at high pressures.
For the case of evaporation with the constant wall heat flux condition, which results in dx/dz = constant, the integration of Eq. (11.87) over a length z with the inlet and outlet qualities being zero and xe, respectively, leads to
p = |
2 f |
TP,avg |
G2 v |
f |
z |
|
+ |
x |
e |
v |
fg |
|
+ G2 |
|
|
|
|
||||
|
|
|
1 |
|
|
|
|
vfg xe |
|
|
|||||||||||
|
D |
|
|
|
|
vf |
|
|
|
||||||||||||
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
g sin (θ )z |
|
|
|
vfg |
|
(11.91) |
||||||||
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
ln 1+ xe |
|
|
|
||||
|
|
|
|
|
|
vfg xe |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
vf |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Note that Eq. (11.91) is based on the inlet condition being the saturated liquid (i.e., inlet quality x = 0). If the inlet is subcooled (i.e., the inlet temperature is less than saturation temperature), then the subcooled length can be determined and single-phase (liquid phase) pressure drop equation can be used over the subcooled length, neglecting the effect of any subcooled boiling on the pressure drop. There are models even for the accurate determination of pressure drop in the subcooled boiling condition (Ramesh and Gedupudi, 2019).
Saturated water and steam properties are given in Tables 11.3 and 11.4, respectively. It is important to note that liquid properties mainly depend on temperature, and the effect of pressure on liquid properties is negligible and can be ignored. There are property tables even for the subcooled (compressed) liquid water, but they are not absolutely necessary. One can determine the liquid properties for the specified (or considered) temperature, making use of the saturated water properties (Table 11.3). For example, for water at 1 bar and 300 K, the saturated water properties corresponding to 300 K can be taken as the water properties.
Gravity head pressure drop is zero for horizontal flow. At higher qualities, the acceleration pressure drop will be high. But then the frictional pressure drop also increases proportionately, some of which comes from the interfacial shear that the homogenous model will not be able to predict (as it assumes equal velocities for both phases). Hence, there is a possibility of the homogenous model underestimating the frictional pressure drop. On the other hand, there is a possibility of the homogenous
