Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Майнор / Книги / Chapter-11---Boiling-and-condensation_2021_Heat-Transfer-Engineering

.pdf
Скачиваний:
0
Добавлен:
08.05.2026
Размер:
2.17 Mб
Скачать

11.3 Pool boiling 361

The boiling Reynolds number Reb is given by

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

q

 

 

 

σ

2

 

ρ f uf Lb

 

 

ρ f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hfg ρ f g

(ρ f ρg )

Reb =

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

µ f

 

 

 

 

 

 

 

 

µ f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

The boiling Nusselt number Nub is given by

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

σ

 

2

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hL

b

 

 

g(ρ f ρg )

 

 

Nub =

 

 

=

 

 

 

 

 

 

 

 

 

 

 

kf

 

 

 

Tsup kf

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(11.18)

(11.19)

The Rohsenow correlation contains an arbitrary constant Csf to account for the influence of liquid-surface combination on the nucleation properties.

The Rohsenow (1952) correlation is

cpf Tsup

 

 

σ

1

n

µ f cpf

 

m+1

q

2

 

 

 

= Csf

 

 

 

 

 

 

 

 

 

hfg

 

 

kf

 

hfg µ f g(ρ f ρg )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Eq. (11.20) can be rearranged as

g(ρ f ρg )

1

 

 

1

2 c

T

n

q = µ f hfg

 

 

 

 

pf sup

 

 

 

m+1

 

σ

 

hfg Csf Prf

 

 

 

(11.20)

(11.21)

where n = 0.33 and m = 0.7. Later, Rohsenow recommended that the value of m be changed to zero for water alone.All fluid properties correspond to the saturated state. The values of the constant Csf for some liquid-surface combinations are given in Table 11.1.

Table 11.1  Values of Csf for the Rohsenow (1952) correlation for various liquid-surface combinations.

Liquid-surface combination

Csf

Water - polished copper

0.0128

n-Pentane - polished copper

0.0154

Carbon tetrachloride - polished copper

0.007

Water - lapped copper

0.0147

Water - scored copper

0.0068

Water - ground and polished stainless steel

0.008

Water - Teflon pitted stainless steel

0.0058

Water - chemically etched stainless steel

0.0133

Water - mechanically polished stainless steel

0.0132

 

 

362CHAPTER 11  Boiling and condensation

11.3.4  Critical heat flux

Zuber (1958) derived an expression for the critical heat flux in pool boiling considering the Taylor and Helmholtz instabilities as

qcrit = qmax = K hfg ρ1/2g σ g(ρ f ρg ) 1/4

(11.22)

where K is a constant that lies between 0.13 and 0.16 based on the experimental data. The recommended value of K is 0.149. Typical qcrit values are in the order of 106 W/m2 for water at near atmospheric pressure.

Example 11.2: Estimate the wall superheat at the onset of departure from nucleate boiling (DNB) or critical heat flux for saturated pool boiling of water on a horizontal surface at atmospheric pressure. Assume the constant Csf to be 0.013 in the Rohsenow correlation.

Solution:

The properties of saturated water at atmospheric pressure are

Tsat = 100 C, ρf = 958 kg/m3 , ρg = 0.598 kg/m3 , hfg = 2257 kJ/kg, σ = 0.05878 N/m, µf = 283.1 × 106 Ns/m2 , cpf = 4.218 kJ/kg K,

Prf = 1.75, Csf = 0.013

Using Eq. (11.22),

qcrit = 0.149 hfg ρ1/2g σ g(ρ f ρg ) 1/4

=0.149 × 2257 × 103 × (0.598)1/2 [0.05878 × 9.81 × (958 0.598)]1/4

=1.26 × 106 W/m2

According to the Rohsenow correlation (this is very widely used and is the workhorse in pool boiling calculations for obtaining the first cut estimates in an engineering problem), Eq. (11.20),

 

q

n

 

σ

 

n/2

 

hfgCsf Prfm+1

 

Tsup =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cpf

 

 

 

 

µ f hfg

g(ρ f

ρg )

 

 

 

 

 

 

 

 

 

1.26

106

 

 

 

 

n

 

0.05878

n/2

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

283.1 10

6

2257 10

3

 

 

 

 

 

 

 

 

 

9.81

(958 − 0.598)

2257 103 0.013 (1.75)m+1

4.218 103

=(1971.96)n (6.258 106 )n/2 6.955 (1.75)m+1

=4.933n 6.955 1.75m+1

The original Rohsenow correlation had values of n = 0.33 and

m = 0.7.

Tsup = 4.9330.33 × 6.955 ×1.751.7

=30.49 K

11.3 Pool boiling 363

For n = 0.33 and m = 0 (recommended by Rohsenow, later for only water),

Tsup = 4.9330.33 × 6.955 ×1.75

= 20.60 K

If one looks at Examples 11.1 and 11.2, it is evident that Tsup for critical heat flux (or DNB) is higher than that is required for nucleation, as it should be.

11.3.5  Film boiling

Berenson (1960) derived an expression for the minimum heat flux (Leidenfrost point E in Fig. 11.2):

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

g(ρ f ρg ) 2

 

σ

4

q

 

= 0.09ρ

h'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

min

g

 

fg

 

ρ f + ρg

 

 

g(ρ f ρg )

 

 

 

 

 

 

 

 

 

 

 

 

 

where,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h'

= h

fg

+ 0.68 c

pg

T

 

 

 

 

fg

 

 

 

 

sup

 

 

(11.23)

(11.24)

Heat transfer in stable film boiling occurs by convection (conduction) through the vapor film and also by radiation across the film.

Berenson (1960) proposed the following expression for the convective heat transfer coefficient for film boiling on flat surfaces:

 

 

 

 

 

 

 

 

 

kg3hfg' ρg g(ρ f ρg )

 

 

 

 

 

h

= 0.425

 

 

 

 

 

 

1/ 2

c

 

 

σ

 

 

 

µg (Tw Tsat )

 

 

 

 

 

 

 

 

 

 

 

 

 

g(ρ f ρg )

1

4

(11.25)

Bromley (1950) proposed the following expression for film boiling on horizontal cylinders:

 

 

 

 

1

 

3 '

4

h = 0.62

kg hfg ρg g(ρ f ρg )

 

(11.26)

µg (Tw Tsat )D

 

c

 

 

 

 

 

 

 

 

 

 

For film boiling on spheres, Eq. (11.26) is used with the constant 0.62 replaced by 0.67.

Note that in Eq. (11.25) and Eq. (11.26), the vapor properties are evaluated at pressure pf and temperature (Tw + Tsat ) / 2 and liquid properties at Tsat .

364 CHAPTER 11  Boiling and condensation

Radiation heat transfer across the vapor film can be calculated using the following expression:

hr =

 

 

1

 

σ (Tw4 Tsat4 )

 

 

 

 

 

 

 

 

1

+

 

1

1 (Tw Tsat )

(11.27)

 

εw

α f

 

 

 

 

 

 

 

where εw is the emissivity of solid (wall), and αf is the absorptivity of liquid, which is equal to the emissivity of the liquid.

Bromley (1950) proposed the following expression for the combined effects of convection and radiation:

hc

1

 

 

3

+ hr

(11.28)

h = hc

 

 

 

h

 

 

 

For hr < hc , Bromley (1950) suggested the following approximation:

h = hc + 0.75hr

(11.29)

11.4  Flow boiling

Flow boiling (or forced convection boiling) occurs when a liquid is forced through a tube or over a surface that is maintained at a temperature greater than the saturation temperature of the liquid. The applications of flow boiling include steam generators for thermal and nuclear power plants. The fluid dynamics and heat transfer mechanisms of flow boiling are more complex than in pool boiling, as the forced flow influences the formation of vapor bubbles, their separation from the heating surface, and the two-phase flow patterns along the tube.

11.4.1  Flow boiling regimes

Fig. 11.5 shows the variation of the wall and fluid temperatures, flow patterns, and heat transfer regions for the flow boiling in a vertical uniformly heated tube (constant wall heat flux). Heat transfer to the subcooled liquid takes place by single-phase convection alone until the wall temperature is higher than the saturation temperature, at which point the boiling commences while the liquid may still be subcooled. There is a decrease in the wall temperature due to the increased heat transfer coefficient as a result of boiling. The size of the vapor bubbles formed increases as the bulk liquid attains saturation temperature. At lower qualities, there are small bubbles dispersed in the liquid (bubbly flow), and the heat transfer mechanism is mainly nucleate boiling. The local heat transfer coefficient mainly depends on the heat flux and is not affected by the quality in nucleate boiling. As the quality (or the void fraction) increases along the flow, bubbles coalesce and form larger bubbles separated by liquid slugs, and the flow pattern is termed as a slug flow. With further increase in the vapor quality along the flow, a vapor core is formed, which is irregular, termed as a churn-annular flow.

11.4 Flow boiling 365

FIGURE 11.5

Wall and fluid temperatures, flow patterns, and heat transfer regions for flow boiling in a vertical heated tube.

A continuous vapor core is then formed, surrounded by the liquid, and is termed as an annular flow. Because of the higher velocities caused by the increased vapor fraction, the nucleation from the heating surface is suppressed. This is due to the thinner thermal boundary layer resulting from the increased velocities. The heat transfer is mainly due to the convective evaporation at the liquid-vapor interface. The vapor core contains entrained liquid droplets formed as a result of the interaction between the liquid and vapor phases. The thickness of the liquid film decreases due to evaporation, and the velocity increases due to the increased quality, both of which result in an increase in the heat transfer coefficient (decrease in the wall temperature) along the flow. At some point, the liquid film completely dries out due to evaporation, and the vapor comes in contact with the wall, which results in a sudden decrease in the

366CHAPTER 11  Boiling and condensation

heat transfer coefficient (increase in the wall temperature). This point is called “dryout”, which is very important from the viewpoint of the design and operating conditions. Further downstream, the evaporation of the entrained liquid droplets increases the velocity, and hence there is a slight increase in the heat transfer coefficient. The vapor flow with entrained droplets is called a spray flow or mist flow, with the wall being liquid deficient. Eventually the tube is completely filled with vapor, at which point a single-phase convection starts. The quality (x) shown in the figure indicates thermodynamic quality obtained from energy balance. The actual quality can be different from the thermodynamic quality, mainly in the subcooled boiling and spray flow conditions.

In actual applications, one may not see all the flow patterns described here, as the flow patterns depend on the inlet condition, mass flux, heat flux, tube diameter, and tube length. For the case of the wall temperature controlled heating, the local heat flux depends on the local two-phase heat transfer coefficient. Fig. 11.6 shows the influence of heat flux on the heat transfer coefficient variation along the flow, for

a constant mass flux. It can be seen that with the increase in heat flux (q2), the heat transfer coefficient increases for the same vapor quality, but the dryout occurs early (at a lower quality). With further increase in the heat flux (q3 and q4), there occurs no

FIGURE 11.6

Heat transfer coefficient variation with quality for constant mass flux and heat flux with q1 < q2 < q3 < q4 < q5.

11.4 Flow boiling 367

transition to annular flow, and there is DNB in the saturated nucleate boiling, which causes a sudden reduction in the heat transfer coefficient. The DNB occurs even in the subcooled boiling at higher heat flux (q5), as shown in the figure. Note that the DNB indicates film boiling similar to that in pool boiling and “dryout” indicates the dryout of liquid film. Normally the wall temperature increase caused by DNB is much higher than that by dryout. Note that q (DNB-subcooled) > q (DNB-saturated) > q (dryout).

It must be noted that boiling (both pool boiling and flow boiling) heat transfer is highly sensitive to the heating surface characteristics. Boiling characteristics change with surface ageing, especially for water (Collier and Thome, 1994; Jayaramu et al., 2019), and hence it is very important that the experimental runs be repeated until the repeatable data is obtained, which can be used for the heat transfer analysis, comparison with the literature data, or development of correlations.

11.4.2  The Chen correlation

By far the most successful and widely used two-phase correlation for flow boiling of water and organic fluids in conventional size tubes is the one proposed by Chen (1963). The correlation predicts the data for both the saturated nucleate boiling region and the convective vaporization region. It assumes that both nucleate boiling and convection occur in some proportion over the entire range of correlation, and their contributions are additive.

hTP = hNcB + hc

(11.30)

where hTP , hNcB, and hc are the local two-phase heat transfer coefficient, the nucleate boiling contribution, and the convection contribution, respectively; hc can be evaluated by a Dittus-Boelter type equation that was presented in Chapter 5.

hc D = 0.023 ReTP0.8 PrTP0.4

(11.31)

kTP

 

Chen (1963) argued that heat is essentially transferred from wall to a liquid film in both dispersed and annular flows, so it is appropriate to use the liquid thermal conductivity in Eq. (11.31). Furthermore, the magnitudes of the liquid and vapor Prandtl numbers are approximately the same, and the two-phase value should not be very different from this value (Prf or Prg ).

Eq. (11.31) can be written as

G (1

x) D

0.8

 

cp 0.4

kf

 

(F)

 

hc = 0.023

 

 

 

 

 

 

 

 

 

 

 

 

 

(11.32)

 

f

 

 

 

 

 

 

 

 

 

 

 

k f

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

where,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.8

 

 

 

ReTP

 

0.8

 

 

 

 

 

 

ReTP

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F =

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

(11.33)

Ref

 

 

 

G (1 − x) D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

368 CHAPTER 11  Boiling and condensation

The parameter F is expressed as a function of the Lockhart-Martinelli parameter, Xtt .

 

 

ρg

0.5 µ f

 

0.1

 

1

x 0.9

Xtt =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ f

µg

 

 

 

 

x

 

F = 1 for

 

1

≤ 0.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xtt

 

 

 

 

 

 

 

 

 

 

1

 

 

0.736

 

 

1

 

F = 2.35

 

 

 

+ 0.213

 

for

 

> 0.1

 

 

 

 

Xtt

 

 

 

 

 

 

Xtt

 

(11.34)

(11.35)

The quantity Xtt2 is the ratio of the frictional pressure gradient with liquid phase alone in the channel to that with vapor (gas) phase alone in the channel, with both phases assumed to be turbulent. Based on the experimental data, the parameter Xtt is related to the parameter F as in Eq. (11.35). With increase in the quality x, the parameter Xtt decreases, which causes an increase in the parameter F. The parameter F can be understood as an “enhancement factor” that enhances the forced convective evaporation with increase in the quality x. As discussed earlier in the Flow boiling regimes section, increasing the quality increases the velocity, which, together with the thinning of liquid film (caused by evaporation), increases the heat transfer coefficient in the annular regime.

Figure 11.4 shows that the liquid superheat is not constant across the thermal

boundary layer. The mean superheat of the liquid ( Te ) in which the vapor bubbles

grow is lower than the superheat at the wall (

Tsup ). Since the thermal boundary

layer is thick in pool boiling, the difference

between Te and Tsup is small.

However, the thinner thermal boundary layer and the steeper temperature gradient in the forced convection result in a significant difference between Te and Tsup.

Chen uses the Forster and Zuber (1955) analysis with Te and the corresponding vapor pressure difference ( pe ) to write the equation for hNcB as

h

= 0.00122

 

 

k0.79f c0.45pf ρ0.49f

 

(

T )0.24

(

 

 

0.5

0.29

0.24

0.24

 

NcB

 

σ

 

e

 

 

 

 

 

µ f

hfg

ρg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A suppression factor, S, is then defined as the ratio of

 

T

0.99

S =

e

 

 

 

Tsup

 

 

 

 

T

0.24

 

T

0.75

S =

e

 

 

e

 

 

 

 

 

Tsup

 

Tsup

 

 

 

 

 

 

 

 

T

0.24

 

p

0.75

S =

e

 

 

e

 

 

 

 

 

 

Tsup

 

psup

 

 

 

 

 

 

 

pe )0.75

Te to Tsup

(11.36)

(11.37)

(11.38)

Eq. (11.38) is based on the Clasius-Clapeyron equation,

11.4 Flow boiling 369

pe

=

psup

=

 

hfg

(11.39)

T

T

v

fg

T

 

 

 

e

 

sup

 

 

sat

 

Substituting Eq. (11.38) in Eq. (11.36),

hNcB = 0.00122

 

 

k0.79f c0.45pf ρ0.49f

 

(

Tsup )

0.24

( psup )

0.75

(S)

(11.40)

 

 

 

 

 

 

 

 

 

σ

0.5

0.29

0.24

0.24

 

 

 

 

 

µ f

hfg

ρg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S approaches unity at low flow rates and zero at high flow rates. Based on the experimental data, S can be expressed as

S =

1

(11.41)

 

1+ 2.53 ×106 ReTP1.17

 

The experimental data used for the development of the Chen correlation considered the two-phase flow in vertical tubes, and the fluids were water (0.55 to 34.8 bar, upward flow and downward flow), benzene, cyclohexane, heptane, methanol, and pentane, all at 1.03 bar with upward flow. The average deviation between the predicted (calculated) and measured data from the experimental cases was ±12%. Though the Chen correlation finds extensive applications for different fluids and operating conditions, it cannot be used for liquid metals or for horizontal tubes due to the asymmetrical nature of flow patterns, such as the stratified flows in horizontal tubes. Many variations of the Chen correlation were proposed later to extend its applicability to a wider range of fluids and operating conditions, and even for different channel geometries (e.g., offset strip fins). In fact, the Chen correlation’s basic philosophy of the nucleate boiling suppression factor S and the liquid-phase convection enhancement factor F has been used even in the development of some of the correlations (called Chen-type correlations) for flow boiling in mini/micro-channels (Jayaramu et al., 2019).

11.4.3  Critical heat flux in flow boiling

The occurrence of the critical heat flux (CHF) and the nature of the CHF (dryout and DNB) in flow boiling were discussed earlier. This section presents the estimation of CHF for flow boiling in vertical uniformly heated round tubes. Flow is in the upward direction.

For uniformly heated tubes, the onset of the CHF occurs first at the tube (or channel) exit, and hence overheating first occurs at the tube exit, in the absence of flow and thermal instabilities.

Boiling cannot occur if the wall temperature is less than the saturation temperature. Therefore, the lower limiting value of the critical heat flux corresponds to the condition of the tube exit wall temperature being equal to the saturation temperature. The fluid temperature at the tube exit can be determined using the inlet condition, the mass flux, and the heat flux. The tube exit wall temperature can be determined using the heat flux, the single phase heat transfer coefficient (hf 0 ), and the fluid

370 CHAPTER 11  Boiling and condensation

temperature at the exit. For a tube of diameter D, length z, mass flux G and the inlet subcooling ( Tsub )i , the lower limiting value of the critical heat flux is

qcrit ,min =

( Tsub )

 

4z

 

i

(11.42)

+ 1

 

 

 

 

 

Gcpf D

 

hf 0

 

The upper limiting value of the critical heat flux corresponds to the condition of the exit quality being equal to 1(x (z) = 1),

 

 

=

GDh

fg

 

+

cpf (

Tsub )

 

 

q

crit ,max

 

1

 

i

 

(11.43)

 

 

 

 

 

 

 

4z

 

 

 

 

hfg

 

 

 

 

 

 

 

 

 

 

 

 

 

Based on the extensive experimental data on the CHF for water in vertical uniformly heated round tubes, the following correlation was proposed by Bowring (1972):

 

 

 

 

 

 

(

 

Tsub )

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A' + DGcpf

 

 

 

 

 

 

 

 

 

qcrit =

 

4

 

 

 

 

(11.44)

 

 

 

 

 

 

 

 

 

 

 

where,

 

 

C' + z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A' =

DG hfg

 

 

F1

 

 

(11.45)

 

 

2.317

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 0.0143 F2 D

1/ 2

 

 

 

 

 

4

 

 

1.0

 

G

 

 

 

0.077F3 DG

 

 

 

 

(11.46)

 

 

C' = 1.0 + 0.347F

(G / 1356)n

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

D, G, and h

fg

are in m, kg/m2s, and J/kg, respectively. The exponent n is given by

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n = 2.0 0.00725 p

 

 

 

 

(11.47)

with the operating pressure p in bar. The pressure functions F1, F2, F3, and F4 are given in Table 11.2.

The empirical correlation was obtained based on the experiments with water in the following parameter ranges:

p : 2 190 bar

D: 0.002 0.045 m z : 0.15 3.7 m

G :136 18600 kg/m2s

The RMS error is 7%, and 95% of the test data lie within ±14%. The empirical correlation does not differentiate between the types of CHF - dryout and DNB (film boiling). Experimental data on the variation of CHF with different parameters is available in Collier and Thome (1994) and Stephan (1992).