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

01 POWER ISLAND / 01 CCPP / V. Ganapathy-Industrial Boilers and HRSG-Design (2003)

.pdf
Скачиваний:
84
Добавлен:
23.06.2022
Размер:
4 Мб
Скачать
☆

8.48

Q:

Discuss the simplified approach to designing fire tube boilers.

A:

Engineers often must estimate the size of heat transfer equipment such as heat exchangers, gas coolers, boilers, and economizers for preliminary costing and to check space requirements. With the approach presented here, one can quickly determine one or more configurations to accomplish a certain amount of heat transfer. One can also size equipment so as to limit the pressure drop without performing lengthy calculations. Life-cycle costing can then be applied to select the optimum design.

Two situations will be discussed [8].

1.The tube-side heat transfer coefficient governs the overall heat transfer. Examples: Fire tube boilers; gas coolers; heat exchangers in which a medium such as air or flue gas flows on the tube side and a fluid with a high heat transfer coefficient flows on the outside. Phase changes can also occur on the outside of the tubes.

2.The shell side governs. Examples: Water tube boilers, steam–air exchangers, and gas–liquid heat transfer equipment. See Q8.49.

Tube-Side Transfer Governs

In a fire tube boiler, gas flows inside the tubes and a steam–water mixture flows on the outside. The gas heat transfer coefficient is small, about 10–20 Btu=ft2 h F, compared to the outside coefficient of 2000–3000 Btu=ft2 h F. The metal resistance is also small; hence the gas-side coefficient governs the overall coefficient and the size of the equipment.

The energy transferred is given by

Q ¼ UA DT ¼ WiCp ðT1 T2Þ

ð89Þ

The overall heat transfer coefficient is obtained from Eq. (4),

1

¼

do

1

 

d

 

do

 

do

 

 

 

 

þ

 

þ

 

ln

 

þ ff i

 

þ ffo

U

hidi

ho

24Km

di

di

Because the inside coefficient governs U, we can rewrite Eq. (4) as follows (neglecting lower order resistances, such as ho, metal resistance, and fouling factors, which contribute to about 5% of U):

di

ð90Þ

U ¼ 0:95hi do

Copyright © 2003 Marcel Dekker, Inc.

The value of the tube-side coefficient is obtained from the familiar Dittus– Boelter equation, Eq. (8),

Nu ¼ 0:023 Re0:8 Pr0:4

where

 

 

 

 

Nu ¼

hidi

;

Re ¼ 15:2

w

12k

dim

The fluid transport properties are evaluated at the bulk temperature. Substituting Eqs. (8)–(11) into Eq. (90) and simplifying, we

following expression [Eq. (12)]:

hi ¼ 2:44w0:8F1=di1:8

where

F1 ¼ ðCp=mÞ0:4k0:6

have the

ð91Þ

Combining Eqs. (89)–(91) we have, after substituting A ¼ 3:14diLN=12 and for flow per tube w ¼ Wi=N,

 

Q

 

¼ 0:606

LN0:2

ð92Þ

D

T F

W 0:8

d0:8

1

i

i

 

This simple equation relates several important variables. Given Q; DT; Wi, and F1, one can try combinations of L; di, and N to arrive at a suitable configuration. Also, for given thermal data, LN0:2=di0:8 is constant in Eq. (92).

F1 is shown in Table 8.35 for flue gas and air. For other gases, F1 can be computed from Eq. (91).

When a phase change occurs, as in a boiler, DT is written as

D

T

ðT1 tsÞ ðT2 tsÞ

ð

93

Þ

¼ ln½ðT1 tsÞ ðT2 tsÞ&

 

 

Combining Eqs. (92) and (93) and simplifying, we arrive at the expression

ln

T1

ts

 

0:606

 

F1

 

N0:2

 

L

 

94

 

T2 ts

¼

Cp

 

Wi0:2di0:8

ð

Þ

 

 

 

 

Factor F1=Cp is also given in Table 8.35.

Equation (94) relates the major geometric parameters to thermal performance. Using this method, one need not evaluate heat transfer coefficients.

Copyright © 2003 Marcel Dekker, Inc.

TABLE 8.35 Factors F1=Cp; F2=Cp; F2, and F3 for Air and

Flue Gasa

Temp ( F)

F1=Cp

F2

F2=Cp

F3

Air

 

 

 

 

100

0.6660

0.0897

0.3730

0.5920

200

0.6870

0.0952

0.3945

0.6146

300

0.7068

0.1006

0.4140

0.6350

400

0.7225

0.1056

0.4308

0.6528

600

0.7446

0.1150

0.4591

0.6810

1000

0.7680

0.1220

0.4890

0.6930

1200

0.7760

0.1318

0.5030

0.7030

 

 

0.1353

 

0.7150

Flue gasa

 

 

 

 

200

0.6590

0.0954

0.3698

0.5851

300

0.6780

0.1015

0.3890

0.6059

400

0.6920

0.1071

0.4041

0.6208

600

0.7140

0.1170

0.4300

0.6457

800

0.7300

0.1264

0.4498

0.6632

1000

0.7390

0.1340

0.4636

0.6735

1200

0.7480

0.1413

0.4773

0.6849

aFlue gas is assumed to have 12% water vapor by volume.

Gas Pressure Drop

Now consider gas pressure drop. The equation that relates the geometry to tubeside pressure drop in in.H2O is

DPi ¼ 9:3 10 5

f

W

 

 

2

 

n

 

 

 

 

i

ðL þ 5diÞ

 

 

 

N

di5

ð95Þ

¼ 9:3 10 5

 

W

 

2

 

 

 

i

 

K2n

 

 

 

N

 

 

 

 

where

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K2 ¼ f ðL þ 5diÞ=di5

 

 

 

 

 

 

ð96Þ

Combining Eqs. (94)–(96) and eliminating N,

 

 

 

ln

T1 ts

0:24

 

 

F1

 

K

 

 

n0:1

 

 

97

 

T2 ts ¼

Cp

1 DPi0:1

ð

Þ

 

 

 

 

where

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K1 ¼ ðL þ 5diÞ0:1Lf 0:1=di1:3

 

 

 

 

 

ð98Þ

Copyright © 2003 Marcel Dekker, Inc.

TABLE 8.36 Values of K1 as a Function of Tube Diameter and Length

 

 

 

 

di

(in.)

 

 

 

 

L (ft)

1.00

1.25

1.50

1.75

2.00

2.25

2.50

2.75

3.00

 

 

 

 

 

 

 

 

 

 

8

7.09

5.33

4.22

3.46

2.92

2.52

2.20

1.95

1.75

10

8.99

6.75

5.34

4.38

3.70

3.17

2.78

2.46

2.21

12

10.92

8.20

6.48

5.31

4.48

3.85

3.36

2.98

2.67

14

12.89

9.66

7.63

6.25

5.27

4.53

3.95

3.50

3.14

16

14.88

11.14

8.80

7.21

6.07

5.21

4.55

4.02

3.61

18

16.89

12.65

9.98

8.17

6.88

5.91

5.15

4.56

4.10

20

18.92

14.16

11.17

9.14

7.70

6.60

5.76

5.10

4.56

22

20.98

15.70

12.38

10.12

8.52

7.31

6.37

5.64

5.05

24

23.05

17.24

13.59

11.11

9.35

8.02

6.99

6.19

5.54

26

25.13

18.80

14.81

12.11

10.19

8.74

7.61

6.74

6.03

28

27.24

20.37

16.05

13.11

11.00

9.46

8.74

7.30

6.52

 

 

 

 

 

 

 

 

 

 

TABLE 8.37 Values of K2 as a Function of Tube Diameter and Length

 

 

 

 

 

 

di (in.)

 

 

 

 

L (ft)

1.00

1.25

1.50

1.75

2.00

2.25

2.50

2.75

3.00

 

 

 

 

 

 

 

 

 

 

8

0.2990

0.1027

0.0428

0.0424

0.0109

0.0062

0.0037

0.0024

0.0016

10

0.3450

0.1171

0.0484

0.0229

0.0121

0.0069

0.0041

0.0027

0.0018

12

0.3910

0.1315

0.0539

0.0252

0.0134

0.0075

0.0045

0.0029

0.0019

14

0.4370

0.1460

0.0595

0.0277

0.0146

0.0082

0.0049

0.0031

0.0021

16

0.4830

0.1603

0.0650

0.0302

0.0158

0.0088

0.0053

0.0033

0.0022

20

0.5750

0.1892

0.0760

0.0350

0.0183

0.0101

0.0060

0.0038

0.0025

22

0.6210

0.2036

0.0816

0.0375

0.0195

0.0108

0.0064

0.0040

0.0027

24

0.6670

0.2180

0.0870

0.0400

0.0207

0.0114

0.0067

0.0042

0.0028

26

0.7130

0.2320

0.0926

0.0423

0.0219

0.0121

0.0071

0.0045

0.0030

28

0.7590

0.2469

0.0982

0.0447

0.0231

0.0217

0.0075

0.0047

0.0031

 

 

 

 

 

 

 

 

 

 

 

K1 and K2 appear in Tables 8.36 and 8.37 respectively, as a function of tube ID and length. In the turbulent range, the friction factor for cold-drawn tubes is a function of inner diameter.

Using Eq. (97), one can quickly figure the tube diameter and length that limit tube pressure drop to a desired value. Any two of the three variables N; L, and di determine thermal performance as well as gas pressure drop. Let us discuss the conventional design procedure:

Copyright © 2003 Marcel Dekker, Inc.

1.Assume w, calculate N.

2.Calculate U, using Eqs. (4) and (90).

3.Calculate L after obtaining A from Eq. (89).

4.Calculate DPi from Eq. (95).

If the geometry or pressure drop obtained is unsuitable, repeat steps 1–4. This procedure is lengthy.

Some examples will illustrate the simplified approach. The preceding equations are valid for single-pass design. However, with minor changes one can derive the relationships for multipass units (e.g., use length ¼ L=2 for twopass units).

Example 1

A fire tube waste heat boiler will cool 66,000 lb=h of flue gas from 1160 F to 440 F. Saturation temperature is 350 F. Molecular weight is 28.5, and gas pressure is atmospheric. If L is to be limited to 20 ft due to layout, determine N and DPi for two tube sizes: (1) 2 1.77 in. (2 in. OD, 1.77 in. ID) and (2) 1.75 1.521 in.

Solution. Use Eq. (92) to find N. Use 2 in. tubes. F1=Cp from Table 8.35

is 0.73 for flue gas at the average gas temperature of 0.5

 

(1160

þ

440)

¼

800 F.

 

 

 

 

440 350

 

¼

 

 

 

 

 

 

 

 

 

 

ln

1160

350

 

 

2:197

 

 

 

 

 

 

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2:197 ¼

0:606 0:73 N0:2

 

 

 

 

 

 

ð66;000Þ0:2 ð1:77Þ0:8

 

 

 

 

 

¼

0:6089N0:2;

N ¼ 611

 

 

 

 

 

 

 

Compute DPi using Eq. (95). From Table 8.37, K2 is 0.035. Compute the gas specific volume.

Density ðrÞ ¼ 28:5

 

 

 

492

 

 

 

¼ 0:031 lb=ft3

 

 

 

 

 

 

359

ð

460

þ

800

Þ

 

 

 

 

 

 

 

n ¼ 32:25 ft3=lb

 

 

 

 

 

 

Substituting into Eq. (95), we have

 

0:035 32:25

DPi ¼ 9:3 10 5

611

 

 

66;000

 

2

 

 

 

 

¼ 1:23 in:H2O

Repeat the exercise with 1.75 in. tubes; length remains at 20 ft. From Eq. (92) we note that for the same thermal performance and gas flow, N0:2L=di0:8 ¼

Copyright © 2003 Marcel Dekker, Inc.

0:076 32:25

a constant. The above concept comes in handy when one wants to quickly figure the effect of geometry on performance. Hence,

6110:2

20

¼ N0:2

20

 

 

 

 

 

0:8

ð1:521Þ

0:8

 

ð1:77Þ

 

 

 

 

N ¼ 333

 

 

With smaller tubes, one needs fewer tubes for the same duty. This is due to a higher heat transfer coefficient; however, the gas pressure drop would be higher. From Table 8.37, K2 ¼ 0.076 for 1.521 in. tubes. From Eq. (95),

DPi ¼ 9:3 10 5 66;000 2 333

¼ 8:95 in:H2O

Example 2

Size the heat exchanger for 2.0 in. tubes with a pressure drop of 3.0 in.H2O. For the same thermal performance, determine the geometry.

Solution. The conventional approach would take several trials to arrive at the right combination. However, with Eq. (97), one can determine the geometry

rather easily:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

1160 350

 

2:197

 

 

0:24

 

 

F1

 

 

K1n0:1

 

440 350

¼

¼

Cp

DPi0:1

 

 

 

 

From Table 8.35, F1=Cp ¼ 0.73; DPi ¼ 3, n ¼ 32.25. Then

ln

1160 350

¼

2:197

¼

0:24K

1 ð

32:25

Þ

0:1

 

 

 

440

 

350

 

 

 

 

 

 

 

 

 

 

 

 

0:73

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¼ 0:222K1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30:1

 

 

 

 

 

 

 

 

 

 

K1 ¼

9:89

 

 

 

 

 

 

 

 

 

 

 

 

From Table 8.36, we can obtain several combinations of tube diameter and length that have the same K1 value and would yield the same thermal performance and pressure drop. For the 1.77 in. ID tube, L is 21.75 ft. Use Eq. (92) to calculate the number of tubes.

N0:2

2:197 ¼ 0:606 0:73 21:75 ð66;000Þ0:2 ð1:77Þ0:8

N ¼ 402

Thus, several alternative tube geometries can be arrived at for the same performance, using the preceding approach. One saves a lot of time by not calculating heat transfer coefficients and gas properties.

Copyright © 2003 Marcel Dekker, Inc.

Life-Cycle Costing

Such techniques determine the optimum design, given several alternatives. Here, the major operating cost is from moving the gas through the system, and the installed cost is that of the equipment and auxiliaries such as the fan. The lifecycle cost is the sum of the capitalized cost of operation and the installed cost:

LCC ¼ Cco þ Ic

The capitalized cost of operation is

C

co ¼

C

Y

1 Y T

 

a

 

1

 

Y

 

 

 

 

 

 

where Y ¼ ð1 þ eÞ=ð1

þ iÞ.

The annual cost

of operating the fan is estimated as

Ca ¼ 0:001 PHCe

where the fan power consumption in kW is

P ¼ 1:9 10 6 Wi DPi rZ

The above procedure is used to evaluate LCC. The alternative with the lowest LCC is usually chosen if the geometry is acceptable. (Ce is cost of electricity.) and H is the number of hours of operation per year.

8.49

Q:

Discuss the simplified approach to designing water tube boilers.

A:

Whenever gas flows outside a tube bundle—as in water tube boilers, economizers, and heat exchangers with high heat transfer coefficients on the tube side— the overall coefficient is governed by the gas-side resistance. Assuming that the other resistances contribute about 5% to the total, and neglecting the effect of nonluminous transfer coefficients, one may write the expression for U as

U ¼ 0:95ho

ð99aÞ

where the outside coefficient, ho, is obtained from

Nu ¼ 0:35 Re0:6 Pr0:3

ð99bÞ

Copyright © 2003 Marcel Dekker, Inc.

where, using Eqs. (16)–(18) and (21),

 

 

hodo

 

Gd

 

 

mCp

Nu ¼

 

 

;

Re ¼

 

 

 

; Pr

¼

 

12k

 

12 m

k

G ¼ 12

 

 

Wo

 

 

 

 

NwLðST doÞ

 

 

 

Equation (99) is valid for both in-line (square or rectangular pitch) and staggered (triangular pitch) arrangements. For bare tubes, the difference in ho between inline and staggered arrangements at Reynolds numbers and pitches found in practice is 3–5%. For finned tubes, the variation is significant.

Substituting Eqs. (17)–(21) into Eq. (99a) and (99b) and simplifying,

ho ¼ 0:945G0:6F2=do0:4

ð100Þ

U ¼ 0:9G0:6F2=do0:4

ð101Þ

where

 

F2 ¼ k0:7ðCp=mÞ0:3

ð102Þ

F2 is given in Table 8.35. Gas transport properties are computed at the film temperature.

A ¼ pdoNwNd L=12

Combining the above with Eq. (89) and simplifying gives

Q=DT ¼ UA ¼ p0:9G0:6F2doNwNd L=12do0:4 ¼ 0:235 F2G0:6NwNd Ldo0:6

Substituting for G from Eq. (21),

Q

0:6

Nw0:4L0:4Nd

 

 

¼ 1:036F2Wo

 

ð103Þ

DT

ðST =do 1Þ0:6

The above equation relates thermal performance to geometry. When there is a phase change, as in a boiler, further simplification leads to

ln

T1

ts

 

2:82

 

F2

 

Nd

 

104

 

T2 ts

¼

Cp

G0:4ðST =do 1Þdo0:4

ð

Þ

 

 

 

If the tube diameter and pitch are known, one can estimate Nd or G for a desired thermal performance.

Let us now account for gas pressure drop. The equation that relates the gas pressure drop to G is Eq. (28) of Chapter 7:

DPo ¼ 9:3 10 10 G2f Nd r

Copyright © 2003 Marcel Dekker, Inc.

For in-line arrangements, the friction factor is obtained from Eq. (29) of Chapter 7:

f ¼ Re 0:15 X

where

X ¼ 0:044 þ

0:08SL=do

ðST =do 1Þ0:43þ1:13do=SL

Another form of Eq. (28) of Chapter 7 is

D

P

 

1:34

 

 

10 7

 

 

Wo1:85nNd2:85m0:15X

 

105

 

o ¼

 

Nw1:85L1:85do0:15ðST doÞ1:85

ð

Þ

 

 

 

 

 

 

Substituting for f

 

in Eq. (28) of Chapter 7 and combining with Eq. (104) we can

relate DPo to performance in a single equation:

 

 

 

DPo ¼ 4:78 10 10 G2:25 ðST doÞ

 

 

 

 

 

 

 

ln

 

 

T1 ts

 

 

X

 

 

106

 

 

 

 

 

T2 ts

do0:75 F3r

ð

Þ

 

 

 

 

 

where

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F3 ¼ ðF2=CpÞm 0:15

 

 

 

 

ð107Þ

F3 is given in Table 8.35. With Eq. (107), one can easily calculate the geometry for a given tube bank so as to limit the pressure drop to a desired value. An example will illustrate the versatility of the technique.

Example

In a water tube boiler, 66,000 lb=h of flue gas is cooled from 1160 F to 440 F. Saturation temperature is 350 F. Tube outside diameter is 2 in., and an in-line arrangement is used with ST ¼ SL ¼ 4 in. Determine a suitable configuration to limit the gas pressure to 3 in.H2O.

Let us use Eq. (106). Film temperature is 0.5 (800 þ 350) ¼ 575 F. Interpolating from Table 8.35 at 475 F, F3 ¼ 0:643. Gas density at 800 F is 0.031 lb=ft3 from Example 1.

DPo ¼ 4:78 10 10 G2:25 ð4 2Þ

ln 1160 350 440 350

ð0:044 þ 0:08 2Þ

20:75 0:643 0:031

¼128 10 10 G2:25 ¼ 3

Copyright © 2003 Marcel Dekker, Inc.

Hence, G ¼ 5200 lb=ft2 h. From Eq. (21) one can choose different combinations of Nw and L:

NwL ¼ 66;000 12=ð2 5200Þ ¼ 76

If Nw ¼ 8, then L ¼ 9.5 ft.

Calculate Nd from Eq. (104):

ln

1160 350

¼

2:197

 

 

440 350

 

 

 

 

F2

 

Nd

 

 

¼

2:82

 

 

 

Cp

G0:4ðST do 1Þdo0:4

2:197 ¼ 2:82 0:426Nd =ð52000:4 1 20:4Þ or Nd ¼ 74

Thus, the entire geometry has been arrived at.

8.50

Q:

How is the bundle diameter of heat exchangers or fire tube boilers determined?

A:

Tubes of heat exchangers and fire tube boilers are typically arranged in square or triangular pitch (Fig. 8.20). The ratio of tube pitch to diameter could range from

1.25to 2 depending on the tube size and the manufacturer’s past practice. Looking at the triangular pitch arrangement, we see that half of a tube area

is located within the triangle, whose area is given by

Area of triangle ¼ 0:5 0:866p2 ¼ 0:433p2

If there are N tubes in the bundle, then

Total area occupied ¼ 0:866Np2

If the bundle diameter is D, then 3:14 D2=4 ¼ area of bundle ¼ 0.866Np2, or

D ¼ 1:05pN0:5

 

 

 

ð108Þ

Similarly, for

the square pitch, the area occupied by one tube

¼

p2

. Hence bundle

2

2

 

 

area ¼ 3:14 D

=4 ¼ Np

, or

 

 

 

D ¼ 1:128pN0:5

 

 

 

ð109Þ

In practice, a small clearance is added to the above number for manufacturing purposes.

Copyright © 2003 Marcel Dekker, Inc.