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

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

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

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 = N1/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 = N1/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

+

1x

(11.88)

 

 

g

f

 

 

 

 

Cicchitti et al. (1960) suggested the relationship

 

= x g + (1x ) 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