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

Николаев Г.И. Тепловые процессы Учебное пособие

.pdf
Скачиваний:
111
Добавлен:
02.05.2014
Размер:
1.56 Mб
Скачать

3.5. Расчет температурных потерь по корпусам

Общие потери температур по корпусам состоят из трех составляющих:

∆ = ∆1 + ∆11 + ∆111 .

а) 1 - от температурной депрессии,

1 = tn tкип. раств ,

В справочных таблицах (табл XXXVI)[2] находим температуру кипения раствора NaNO3 при атмосферном давлении:

 

Концентрация Температура Температура Депрессия,

 

NaNO3, %

кипения

кипения

град

 

 

 

раствора,

H2O, град

 

 

 

 

 

град

 

 

 

x1

= 15,2

102

100

1

= 2

IIк

x2

= 21,6

103

100

1

= 3

IIIк

x3

= 40,0

107

100

1

= 7

Поправочный коэффициент на искомую температуру

 

определяется по табл. 2.2.

t1

= 130

11

= f1 * 1атм =1,18* 2 = 2,36град;

IIк

t2

= 110

12

= f2 * 1атм =1,07 *3 = 3,21град;

IIIк

t3 = 60

13

= f3 * 1атм = 0,76 * 7 = 5,32град;

_____________________

Итого: ∑∆1 =10,83град.

б) от гидростатического эффекта:

11 = tP tP

Hср = 0,4 м – принимаем расстояние от верхнего уровня жидкости до середины омываемой жидкостью поверхности нагрева.

 

 

 

 

 

 

 

 

 

 

 

P = P + ρ

1

gH

ср

 

= 267,8*103 +1056 *9,81*0,4 = 272 *103 Па = 2,78атм; IIк

1

 

 

 

 

 

 

 

 

 

P = P + ρ

2

gH

ср

=143,7 *103

+1102 *9,81*0,4 =148*103 Па =1,51атм; III

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к

P = P = ρ

3

gH

ср

=19,6 *103

+1255 *9,81* 0,4 = 24,52 *103 Па = 0,25атм.

3

3

 

 

 

 

 

 

 

 

Температура кипения воды при давлении по корпусам и потери температур от гидравлического эффекта:

 

 

t

P1

= f (∆P = 2,78атм)=131,52οС;

 

 

 

 

1

 

t

P1

= f (P = 2,74атм)=130,93οС;

 

 

1

 

 

 

 

 

 

 

 

111

=131,52 130,93 = 0,59град

 

 

 

IIк

t

P2

= f (∆P =1,51атм)=110,9ο

С;

 

 

= f (P

 

 

2

 

t

P2

=1,46атм)=109,9οС;

 

 

2

 

 

 

 

 

 

 

 

112

=110,9 109,9 =1град

 

 

 

IIIк

t

P3

= f (∆P = 0,25атм)= 64,7οС;

 

 

 

 

 

3

 

t

P3

= f (P

= 0,2атм)= 59,7οС;

 

 

3

 

 

 

 

 

 

 

 

113

= 64,7 59,7 = 5град

 

 

 

Итого: ∑∆11 = 6,59град

 

 

 

в) от гидродинамического эффекта (∆111 ).

Потери разности температур на каждом интервале между корпусами принимаем в 1, тогда 111 =1*3 = 3град.

Сумма всех потерь разности температур в целом составит: ∆ =10,83 + 6,59 +3 = 20,42град

1 = 2,36 + 0,59 +1 = 3,59 4град.

3.6.Полезная разность температур

tпол = tп t3 ∆ =143 60 20,42 = 62,58град; tпол1 t1 − ∆1 =143 130 4 = 9град.

82

81

3.7. Температура кипения растворов по корпусам

 

 

 

 

 

tкип = tвп

+ ∆ град

 

 

 

тогда,

 

 

 

 

 

 

 

 

 

 

 

IIIк tкип3

= 60 +5,32 +5 +1 = 71,32град;

 

 

 

 

IIк

tкип2

=110 +3,21+1+1 =115,21град;

 

 

 

tкип1

=130 + 2,36 + 0,59 +1 =133,95град.

 

 

 

Расчетные параметры растворов и паров по корпусам

 

представлены в табл. 2.3.

 

 

Таблица 2.3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Параметры

 

Един.изм

 

 

 

Корпуса

 

 

 

 

 

 

.

 

I

 

II

III

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

W-произво-

 

 

кг/с

 

0,295

 

0,324

0,351

0,97

 

 

дительность

 

 

 

 

 

 

 

 

 

 

 

по удаляемой

 

 

 

 

 

 

 

 

 

 

влаге

 

 

 

 

 

 

 

 

 

 

2

X-

 

 

%

 

15,2

 

21,6

40

=

 

концентрация

 

 

 

 

 

 

 

40

 

 

растворов

 

 

 

 

 

 

 

 

 

 

3

P-давление

 

 

МПа/

 

267,8/2,74

 

143,7/1,47

19,6/0,2

Рп =

4

 

паров

 

 

атм

 

 

 

 

 

атм

 

4

tr- температу-

 

град

 

130

 

110

60

tп = 143

 

ра паров

 

 

 

 

 

 

 

 

 

 

5

∑∆ -

 

 

град

 

3,95

 

5,21

11,32

20,53

 

 

температур-

 

 

 

 

 

 

 

 

 

 

 

ные потери

 

 

 

 

 

 

 

 

 

 

6

tкип

 

град

 

133,95

 

115,21

71,32

 

 

 

температура

 

 

 

 

 

 

 

 

 

 

 

кипения

 

 

 

 

 

 

 

 

 

 

 

раствора

 

 

 

 

 

 

 

 

 

 

7

tпол

 

град

 

8,95

 

18,74

43,83

62,42

 

 

полезная

 

 

 

 

 

 

 

 

 

 

 

разность

 

 

 

 

 

 

 

 

 

 

 

температур

 

 

 

 

 

 

 

 

 

 

83

3.8.Расчет коэффициентов теплопередачи по корпусам

Коэффициент теплопередачи определяют по уравнению

аддитивности термических сопротивлений:

K =

1

 

 

 

.

 

1

+δλ

+

1

 

 

α

α

 

1

 

2

 

 

3.8.1.Примем, что суммарное термическое

сопротивление равно термическому сопротивлению стенки δ/ λст и накипи δн / λн . Выбираем конструкционный

материал, стойкий в среде кипящих водных растворов. В этих условиях подходит сталь марки X17 с коэффициентом

теплопроводности

 

λст

= 25,1 Вт/м˚К. Термическое

сопротивление загрязнений со стороны пара не учитываем.

Тогда δ

=

0,002

+

0,0005

= 0,287 *103 , [м2К/Вт].

25,1

 

λ

 

2,42

 

3.8.2. Коэффициент теплоотдачи от конденсирующего пара к стенке первого корпуса.

Определяем по уравнению (2.6) при Рп = 4 атм и Н = 4 м.

α= λ3 ρ2 r= ∆ 0,25 2

12,04 * 4 µ1 Н1t A* t1 , Вт/м К.

11

А= 10650 при Рп = 4 атм из табл. 2.1.

Принимаем t1 = 2 , тогда

α= A * t 0,25 = 1065020,25 = 8955,5Вт/ м2 К .

Для установившегося процесса передачи тепла справедливо уравнение:

q =α1 *t1 =

1

 

*tст =α2 *t2 .

(2.36)

δ

/ λ

 

 

 

 

 

 

84

 

Распределение температур в процессе теплопередачи от пара через стенку к кипящему раствору показано на рис. 2.6.

Рис. 2.6. Распределениетемпературвпроцессетеплопередачи отпарак кипящемурастворучерезстенку:

1 – пар; 2 – конденсат; 3 – стенка; 4 – ржавчина; 5 – кипящий раствор.

3.8.3.Коэффициент теплоотдачи от стенки к

кипящему раствору для

пузырькового кипения

в

вертикальных

трубках

при

условии

естественной

циркуляции раствора определяется по уравнению (2.8).

 

 

 

85

 

 

 

 

 

α

 

= B * q0,6

= 780

 

λ

ρ0,5

ρ0,06

 

 

* q0,6.

 

 

 

 

2

 

2

2

вп

 

 

 

 

 

 

δ 0,5r0,6 ρ0,66с0,3

µ0,3

 

 

 

 

 

 

 

 

 

2

вп

0

2

2

 

 

 

 

 

 

 

Физические свойства кипящего раствора представлены в

табл. 2.4.

 

 

 

 

 

 

 

 

 

 

Таблица 2.4

 

 

 

 

 

 

 

 

 

 

Корпуса

Параметры

Ед.изм

 

 

 

 

Лит-ра

 

 

 

 

 

 

 

 

 

I

 

 

II

 

III

 

 

 

 

 

 

 

 

 

 

 

 

1

λ – теплопро-

Вт/моК

 

0,61

 

0,62

0,69

[2, 3]

 

 

водность

 

 

 

 

 

 

 

 

 

 

 

 

 

раствора

 

 

 

 

 

 

 

 

 

 

 

 

2

ρ - плотность

кг/м3

1062

 

1104

1399

[3]

 

 

раствора

 

 

 

 

 

 

 

 

 

 

 

 

3

с -

 

 

Дж/

3771

 

3561

2765

[2, 3]

 

 

теплоемкость

кг К

 

 

 

 

 

 

 

 

 

 

 

раствора

 

 

 

 

 

 

 

 

 

 

 

 

4

µ

- вязкость

Па*с

 

0,1*

 

0,29*

0,7*

[2, 3]

 

 

раствора

 

 

 

10-3

 

10-3

10-3

 

 

5

δ – поверхно-

Н/м

 

 

 

 

 

 

 

 

[2, 3]

 

 

стное

 

 

 

 

 

 

 

 

 

 

 

 

 

 

натяжение

 

 

 

 

 

 

 

 

 

 

 

 

 

раствора

 

 

 

 

 

 

 

 

 

 

 

 

6

r

 

- теплота

Дж/кг

2068

 

2148*

2372

[2]

 

 

парообразова-

 

 

 

*103

 

103

*103

 

 

 

ния

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

ρп - плотность

кг/м3

 

3,75

 

2,0

0,098

[2]

 

 

пара

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

ρо

- плотность

кг/м3

0,579

 

0,579

0,579

[2]

 

 

пара

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при Р = 1 атм

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

86

 

 

 

 

 

 

 

α2 = 780*q0,6

0,611,3 *10620,5 *3,750,06

 

=

0,0580,5 *(2068*103 )0,6 *0,5790,66 *37710,3 *(0,1*103 )0,3

 

 

=18,76*q0,6 =18,76*(α1 *t1 )0,6 =18,76*(8955,5*2)0,6

= 6685Вт/ м2 К

В = 18,76.

Рис. 2.7. Зависимость удельной тепловой нагрузки q от разности температур ∆t1

Тепловую нагрузку со стороны кипящего раствора определим по следующей формуле:

q12*∆t2= 6685*2 = 13370 Вт/м2, где ∆t2 = ∆tпол1- ∆tст – ∆t1 = 9-5-2 = 2 0С;

∆tст – разность температур стенки со стороны пара и кипящего раствора, определяется по формуле (2.39):

∆tст= α1* ∆t1* ∑δ/λ = 8955,5*2*0,287*10-3 = 5 град.

Проверим правильность первого приближения по равенству удельных тепловых нагрузок:

q11*∆t1 = 8955,5*2 = 17911 Вт/м2; q22*∆t2 = 6685*2 = 13370 Вт/м2.

Отсюда

∆ =

q1 q2

*100

=

1791113370

*100 = 33,96% .

q2

 

13370

 

 

 

 

 

 

 

 

87

 

Как видим разница между q1 и q2 очень большая, больше 10%, принимаем второе приближение ∆t1=1 град.

Пренебрегая изменением физических свойств конденсата при изменении температуры на 1 град, рассчитываем α1 по

соотношению:

α1 = А* ∆t1-0,25 = 10650*1 = 10650 Вт/м2К; q1 = α1*∆t1 = 10650*1 = 10650 Вт/м2.

Тогда получим:

∆tст = 10650*1*2,87*10-4 = 3 град;

∆t2 = 9-1-3 = 5 град;

α2 = В*q10,6 = 18,76(10650)0,6 = 4893,76 Вт/м2К; q2 = α2*∆t2 = 4893,76*5 = 24469 Вт/м2.

Очевидно, что q1≠q2.

Для расчета в третьем приближении строим графическую зависимость удельной тепловой нагрузки q от разности температур между паром и стеной в первом корпусе (рис. 2.7.) и определяем ∆tопт = 1,7 град.

Отсюда получим:

α1 = А* ∆tопт-0,25 = 10650/1,70,25 = 9326 Вт/м2К; q1 = α*∆tопт = 9326*1,7 = 15856 Вт/м2;

∆tст= α1* ∆tопт* ∑δ/λ = 9326*1,7*0,287*10-3 = 4,55 град;

∆t2 = 9-4,55-1,7 = 2,75 град;

α2 = В(q)0,6 = 18,76(15856)0,6 = 6214 Вт/м2К; q2 = α2*∆t2 = 6214*1,7 = 17087 Вт/м2.

Определим процент отклонения между q1 и q2

∆ =

q2 q1

*100

=

17087 15856

*100 = 7,76% .

q

 

15865

 

 

 

 

1

 

 

 

 

 

 

 

 

 

88

 

Процент отклонения допустимого < 10% для

технических расчетов, тогда принимаем

α1 = 9326 Вт/м2К и α2 = 6214 Вт/м3К.

Следовательно,

К1 = 93261 + 0,2871*103 + 62141 =1828Вт/ м2 К .

Ориентировочное соотношение коэффициентов теплопередачи по корпусам при выпаривании водных растворов солей К1 : К2 : К3 = 1 : 0,58 : 0,34.

На основании таких ориентировочных соотношений коэффициенты теплопередачи составят:

IК

К1

= 1828 Вт/м2К;

IIК

К2 = 1060 Вт/м2К;

IIIК

К3

= 621 Вт/м2К.

3.9. Составление тепловых нагрузок по корпусам

По условию, раствор подается на выпарку подогретым до температуры кипения.

Тогда раствор тепла составит по корпусам: IК

Q1 =W1 * r1 = 0,295* 2179 *103 = 643000Вт

 

 

 

 

IIК

 

G *C (t

 

 

 

 

)= 0,324*2233*103

1,09*3,69*103

 

 

Q

2

=W

2

*r

1

t

2

*

г

 

 

 

 

2

 

 

1

 

 

1

 

 

 

 

 

 

*(134 115)

= 647000Вт,

 

 

 

 

 

 

де

G1 = GH W1

=1,39 0,295 =1,09кг/ с

 

 

 

 

 

 

 

IIIК

 

 

 

 

 

(t

 

 

 

 

)= 0,351* 2359*103 0,771*3,43*103 *

 

 

Q

3

=W

3

* r G

2

*C

2

2

t

3

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

*(115 73)

= 722000Вт,

 

 

 

 

 

 

 

 

где G2

 

= GH W1 W2 =1,39 0,295 0,324 = 0,77кг/ с.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

89

 

 

 

3.10. Расход греющего пара в I корпусе:

D =

Q1

=

643000

= 0,3кг/ с.

 

 

 

 

 

 

r

2140

 

 

 

 

1

 

 

 

 

 

Удельный расход пара:

d =

D

=

 

0,3

= 0,31 кг пара D / кг удаляемой влаги W.

W

0,97

 

 

 

 

3.11. Распределение полезной разности температур по корпусам

Распределение полезной разности температур по корпусам сделаем в двух вариантах из условий минимальной общей поверхности и из условий равной поверхности корпусов,

т. е. пропорционально Q

и

Q

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

K

 

 

 

 

 

 

 

 

 

Найдем факторы пропорциональности:

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

Q

*10

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

K

 

IК

 

Q1

K1

=

643

 

= 0,593

 

 

Q1

 

=

643000

= 351 или 0,351

 

 

 

1828

 

 

 

K1

 

1828

 

 

IIК

Q2

K2

=

647

0,781

 

 

Q2

K 2

=

647000

= 610 или 0,61

 

 

 

 

1060

 

 

 

 

 

 

 

1080

 

 

III

К

Q3

K

 

=

722

=1,078

 

 

Q3

 

=

722000

= 1162 или 1,16

 

 

 

621

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

K3

 

621

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q K = 2,452

 

 

 

 

 

 

 

 

 

Q K = 2,121

90

Полезные разности температур по корпусам:

Вариант минимальной Вариант равной поверхности поверхности

Itпол1 = tпол * Q1 K1 =

К Q K

 

=

62,58*0,593 =15,13град

 

tпол *Q1

K1

 

 

 

 

62,58* 0,351

 

 

 

2,452

 

 

 

 

 

 

 

 

=

 

 

=

IIК

 

 

 

 

 

 

Q

 

 

 

2,121

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q2 K2

 

 

 

 

K

 

 

 

 

 

 

 

 

tпол2

= tпол *

=

=10,35град

 

 

 

 

 

 

 

 

 

Q K

 

 

62,58* 0,61

=18,0град

 

=

62,58* 0,781

=19,93град

 

 

 

 

2,121

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,452

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

IIIК

 

 

 

 

 

62,58*1,16

 

 

 

 

 

 

 

 

tпол *

Q3

K3

 

 

= 34,23град

 

tпол3

 

 

2,121

 

 

 

=

Q K

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

62,58*1,078

= 27,51град

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,452

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

91

3.12. Определение поверхности нагрева

F =

 

 

Q

=

 

 

 

643000

= 23м2

 

F1 =

 

 

Q1

=

 

 

643000

 

= 33,9м

2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K1tпол1

1828 *10,35

 

K1

tпол1

1828*15,13

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F =

 

 

Q2

 

=

647000

 

= 33,9м2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q2

 

 

 

 

 

647000

 

 

 

2

 

 

K2 tпол2

 

 

 

 

1060 *18,0

 

 

F =

 

 

 

=

 

 

 

 

= 30м

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q3

 

 

 

 

722000

 

 

 

 

K2 tпол2

 

1060 *19,93

 

F =

 

 

=

 

= 33,9м2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

K3tпол3

 

 

 

 

621* 34,23

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F =

 

 

Q3

=

 

 

 

722000

= 30м2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

K3tпол3

 

 

621*38,68

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F = 83м2

 

 

 

 

 

 

 

F =101,7м2

 

 

 

В связи с небольшой

 

разницей поверхностей двух

 

 

 

 

 

вариантов, принимаем вариант равной поверхности,

 

обеспечивающей однотипность аппарата.

 

 

 

 

 

 

 

 

 

 

 

 

Сравнение распределенных из условий равенства

 

поверхностей теплопередачи и предварительно рассчитанных

 

значений полезных разностей температур представлено ниже:

 

Распределение

 

 

 

в

 

первом

 

 

IК

 

 

 

 

 

 

IIК

 

IIIК

 

 

 

 

 

10,35

 

 

 

 

 

18,0

34,23

 

приближении значения tпол , оС

 

 

8,95

 

 

 

18,74

43,89

 

Предварительно

рассчитанные

 

 

 

 

 

 

 

 

 

 

 

 

 

 

значения tпол,

оС

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как видно, значения полезных разностей температур отличаются незначительно, что позволяет не производить уточненный расчет.

Учитывая большую разность температур в третьем корпусе, произведем уточненный расчет (как пример).

92

3.13.Уточненный расчет теплопередачи

Всвязи с тем, что существенное изменение давлении, по сравнению с рассчитанным в первом приближении, происходит только в 1-м и 2-м корпусах, где суммарные потери незначительны, во втором приближении принимаем

такие же значения 1 , 11 и 111 как в первом приближении.

Полученные после перераспределения температур (давлении) параметры растворов и паров по корпусам представим в табл. 2.5.

Таблица 2.5

Параметры

Един.

 

Корпуса

 

 

 

измер.

I

II

III

1

W – количество

кг/с

0,295

0,324

0,351

 

удаляемой влаги

 

 

 

 

2

X – концентрация

%

15,2

21,6

40

 

раствора

 

 

 

 

3

tn – температура

град

130

 

 

 

греющего пара

 

 

 

 

4

tпол - полезная

град

10,35

18,0

34,23

 

разность температур

 

 

 

 

5

tкип – температура

град

119,65

94,1

54,66

 

кипения раствора

 

 

 

 

 

tкип = t1 − ∆tпол

 

 

 

 

6

tвп - температура

град

116,7

89,89

44,34

 

вторичного пара

 

 

 

 

 

tвп = tкип − (∆1 + ∆11 )

 

 

 

 

7

Рвп – давление

атм

1,814

0,715

0,0977

 

вторичного пара

 

 

 

 

8

tr = tвп − ∆111

град

 

115,7

88,89

 

 

93

 

 

 

а) Расчет тепловых нагрузок.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q = D * r

 

=1,03*W

* r =1,03* 0,295* 2217 *103 = 674кВт,

 

 

1

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

где «1,03» - 3% потерь тепла в окружающую среду.

 

 

r = f (P

 

=1,814)= 2217 *103 Дж/ кг ;

 

 

 

1

 

 

 

 

 

 

вп1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r = f

(Р = 4атм)= 2141*103 Дж / кг;

 

 

 

D

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D =

 

674

 

= 0,31кг / с;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2141

 

 

 

 

 

 

tкип2 )]=

 

 

 

 

 

 

 

 

 

 

 

 

 

G1 * C1 (tкип1

 

 

 

 

Q2

=1,03[W2 * r2

 

 

 

 

=1,03[0,324 * 2285 1,09 * 3,651(119,68 94,1)]658кВт,

 

 

Q3

=1,03[W3 * r3

G2 * C2 (tкип2 tкип3 )]=

 

 

 

=1,03[0,351* 2390 0,771* 2,765(94,1 54,66)]= 777кВт.

 

 

б) Расчет коэффициентов теплопередачи.

 

 

 

К1

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

674 *10

3

 

=1920Вт/ м2 К;

 

 

=

 

 

 

 

 

 

 

 

1

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

33,9 *10,35

 

 

 

 

 

 

F * tпол1

 

 

 

 

 

 

 

 

 

 

 

 

К2

=

 

 

 

 

 

Q

2

 

 

 

 

=

 

658 *103

=1078Вт/ м2 К;

 

 

 

 

F * t

 

 

 

33,9 *18

 

 

 

 

 

 

 

пол2

 

 

 

 

 

 

 

 

 

 

 

 

 

К3

=

 

 

 

 

 

 

Q3

 

=

 

777

 

 

= 669Вт/ м2 К .

 

 

 

 

 

 

 

 

 

 

 

 

33,9*34,23

 

 

 

 

 

 

F * tпол3

 

 

 

 

 

 

 

 

 

 

 

в) Распределение полезной разности температур:

 

 

 

 

 

 

 

 

 

 

Q1

 

K1

 

 

 

 

 

 

 

 

 

 

 

 

674

 

 

 

 

 

 

t

 

= ∑∆t

 

 

 

 

 

 

= 62,58

 

 

 

 

 

1920

 

 

 

 

 

=10,32град;

 

 

 

 

 

 

 

 

 

 

 

 

674

 

 

658

 

 

777

 

 

 

пол1

 

 

 

пол Q K

 

 

 

 

+

 

+

 

 

 

 

 

 

 

 

 

 

 

 

0,61

 

 

 

 

 

 

 

 

 

 

 

1920

 

1078

669

 

 

 

tпол2

= 62,58

 

=18град;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tпол3 = 62,58

1,16

 

= 34,24град.

 

 

 

 

 

 

 

 

 

 

2,12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

94

 

 

 

 

 

 

 

 

Проверка суммарной полезной разности температур:

tпол =10,32 +18 +34,24 = 62,56град.

Сравнение полезных разностей температур, полученных во 2-ом и 1-ом приближениях:

Значение ∆tпол во 2-м

I к

II к

III к

10,32

18,0

34,24

приближении, 0С

Значение ∆tпол в 1-м

10,35

18,0

34,23

приближении, 0С

Различия между полезными разностями температур по корпусам в 1-ом и 2-ом приближениях не превышает 1 %. Если же разница превышает 5 %, необходимо выполнить следующее, взяв за основу расчет ∆tпол из 2-го приближения и т.д., за совпадения полезных разностей температур.

г) Поверхность теплопередачи выпарных аппаратов:

F =

 

 

 

Q

 

 

=

 

674 *10

3

 

= 34м2 ;

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

k1

tпол1

 

 

1920 *10,32

 

 

 

 

 

F =

 

 

 

Q

2

 

 

=

658*103

 

= 33,9м2 ;

 

 

 

 

 

 

 

 

 

 

 

2

 

k2

tпол2

 

 

 

1078*18

 

 

 

 

 

 

 

 

 

 

 

F =

 

 

 

Q

 

 

=

 

777 *103

 

= 33,9м2 .

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

k3

tпол3

 

 

 

669 *34,24

 

 

 

 

 

 

 

95

По ГОСТ 11987-81 (главаV)[4] выбираем аппарат со следующими характеристиками:

Номинальная поверхность

FН =40 м2

теплообмена

d = 38*2 мм

Диаметр труб

h = 4000 мм

Высота труб

D = 600 мм

Диаметр греющей камеры

D1 = 1200 мм

Диаметр сепаратора

D2 = 400 мм

Диаметр циркуляционной

На = 15000 мм

трубы

Ма = 3000 кг

Общая высота аппарата

 

Масса аппарата

 

В дальнейшем этот аппарат представить с чертежом в пояснительной записке – раздел подбор и принцип работы основного оборудования.

3.14. Расчет барометрического конденсатора

Для создания вакуума в выпарных установках обычно применяют конденсаторы смешения с барометрической трубой. В качестве охлаждающего агента используют воду, которая подается в конденсатор с перфорированными полочками, при температуре окружающей среды (около 20 0С). Вода с конденсатом выводится из аппарата через барометрическую трубу высотой до 10 м, т.к. 1 атм соответствует 10 м. вод. ст. для поддержания постоянства вакуума в системе, одновременно откачивают неконденсирующиеся газы вакуум-насосом.

Необходимо рассчитать расход охлаждающей воды, основные размеры барометрического конденсатора и барометрической трубы, производительность вакуум-насоса.

96

3.14.1. Определение расхода охлаждающей воды. Расход охлаждающей воды определяют из теплового

баланса конденсатора:

Gв =

Wв (Iбκ Cв *tk )

,

(2.37)

Cв (tk tн )

 

 

 

где Iδк – энтальпия паров в барометрическом конденсаторе, Дж/кг;

tк – конечная температура смеси воды и конденсата,

0С;

tн – начальная температура охлаждающей воды, 0С.

Разность температур между паром и жидкостью на выходе из конденсатора должна быть 3÷5 град.

Поэтому конечную температуру воды на выходе из конденсатора принимают на 3÷5 град ниже температуры

конденсации паров:

tk = tбk – 3 = 60-3 = 57 0С.

Тогда

= 0,351(2608,3 4,18*57) =

Gв 5,38кг/ с. 4,18(57 20)

3.14.2. Расчет высоты барометрической трубы.

В соответствии с нормами внутренний диаметр барометрической трубы dбт=300 мм.

Скорость воды в барометрической трубе Uв равна:

Uв =

4(Gв +W3 )

=

 

4(5,38 + 0,381)

 

= 0,08м/ с

 

983*3,14 * 0,32

 

 

ρвπd ебт

 

 

 

 

 

 

 

 

Высоту барометрической трубы определяют по

формуле:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В

 

 

 

 

 

 

 

 

2

 

 

Нбт =

 

 

 

+ ζ +

Hбт

 

 

Uв

 

,

 

+ 1

 

 

*

 

 

+ 0,5

ρв * g

dбт

 

2g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

97

 

 

 

 

 

 

 

где В – вакуум в барометрическом конденсаторе, Па; 0,5 м – запас высоты на возможное изменение

барометрического давления.

ζ =ζвх +ζвых = 0,5 +1,0 =1,5; B=Pатм-Pбк=9,8*104-2*104=7,8*104Па;

Re =

4 * dбт * ρв

=

U * dбт

=

0,08* 0,3

= 60000 ;

 

 

 

 

0,4 *106

 

 

 

 

µв

 

 

 

 

ϑв

 

 

 

 

 

 

 

 

 

 

68

 

0,25

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

Re 560

 

λ = 0,11

 

 

 

+

 

 

 

 

при

10

 

 

 

 

; (2.38)

 

 

 

 

 

 

е

e

 

 

dбт

60000

 

 

 

 

 

 

 

 

 

λ =

64

 

при Re<2320;

 

 

 

 

 

 

 

 

 

 

(2.39)

Re

 

 

0,316

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ =

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Re0,25

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для гладких труб при 232 Re 10

1

 

 

 

 

(2.40)

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ=0,11e0,25;

для турбулентного режима при Re > 560 1e , (2.41)

отсюда при Re = 60000

 

0,2

 

68

 

0,25

λ = 0,11

 

+

 

 

= 0,023 ,

300

60000

 

 

 

 

где ∆ = 0,2мм для стальных труб умеренного загрязнения

[4, глава 1]:

Нбт =

7,8*104

 

 

 

Нбт

0,08

 

= 8,6м.

 

+ 1

+1,5

+0,023

 

*

 

+ 0,5

983*9,8

0,3

2 *9,8

 

 

 

 

 

 

 

 

 

 

 

 

98

 

 

 

 

3.14.3. Диаметр барометрического конденсатора. Диаметр барометрического конденсатора Dбк

определяют из уравнения расхода:

Dбк =

4W3

 

,

ρn *U n

*π

 

 

гдеρn плотность пара, кг/м3;

Un – скорость пара в конденсаторе принимают 15 – 25

м/с при остаточном давлении порядка 104.

Тогда Dбк =

4 * 0,351

= 0,42м .

0,128* 20 *3,14

 

 

Выбираем Dбк = 500 мм [4, глава V, приложение 5]. 3.14.4. Расчет производительности вакуум-насоса. Производительность вакуум-насоса Gвозд определяется

количеством газа (воздуха), который необходимо удалить из барометрического конденсатора:

G

возд

= 2,5*105 (W +G

в

)+ 0,01W ,

(2.42)

 

 

3

 

3

 

где 2,5*10-5 – количество газа, выделяющегося из 1 кг

воды;

0,01 – количество газа, подсасываемого в

 

 

конденсатор через неплотности на 1 кг паров.

 

Тогда

 

 

 

 

 

Gвозд = 2.5*105 (0,351+5,38)+ 0,01* 0,351 = 3,65*103 кг/ с

Объемная производительность вакуум-насоса равна:

Vвозд

=

R(273 +tвозд )*Gвозд

,

 

 

 

где

 

M возд * Рвозд

 

 

 

 

 

tвозд = tn + 4 + 0,1(tк tн )= 60 + 4 + 0,1(57 20)= 67,7οС

тогда

 

 

 

 

 

Vвозд =

8310(273 + 67,7)*3,65 *103

= 0,018м3 / с =1,09м3

/ мин(2.43)

 

 

 

 

29 *19,6 *103

 

 

99

 

 

 

 

 

 

 

 

 

Согласно ГОСТа 1867-57 [4, глава V, приложение 6] принимаем вакуум-насос производительностью 1,5 м3/мин.

4. Пример расчета выпаривания пищевых продуктов

(2-46)

Пример 4.1. Рассчитать двухкорпусную выпарную установку непрерывного действия для сгущения томатной массы в количестве G=1,5кг\с от Хн =5% до Хк,=20%. Масса поступает на выпаривание подогретой до температуры кипения. Давление пара, греющего I корпус p=0,13 Мпа, остаточное давление вторичного пара, поступающего в барометрический конденсатор, pк=12КПа. Схема установки приведена на рис.2.7.

Решение. Количество воды выпаренной в двух корпусах установки, определяем по уравнению (2.2):

 

 

W=1,5(1-5/20)=1,125кг/с.

 

На

основании

практических

данных

примем

следующее

соотношение массовых

количеств

воды,

выпариваемой по корпусам:

 

 

 

 

I : II=1,0 : 1,1.

 

 

Тогда количество воды выпариваемой по корпусам составит в I корпусе: W1=1,125 1,0 /(1+1,1)=0,536кг/с;

во I корпусе: W2=1,125 1,1/(1+1,1)=0,589кг/с,

________________

ИтогоW=1,125кг/с.

Из первого корпуса во второй переходит томатной

массы:

G1=G-W1=1,5-0,536=0,964 кг/с.

100