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


ОГЛАВЛЕНИЕ 

стр.

ВВЕДЕНИЕ                     5

1. Определение поверхности теплопередачи выпарных аппаратов            6

2. Расчет концентраций упариваемого раствора                        6

3. Определение температур кипения растворов                        7

4. Расчет полезной разности температур                      8

5. Определение тепловых нагрузок                            8

6. Выбор конструкционного материала                         11

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

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

9. Уточненный расчет поверхности теплопередачи                     19

10. Определение толщины тепловой изоляции                           25

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

12. Определение расхода охлаждающей воды                       27

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

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

15. Расчет производительности вакуум-насоса                       28

ЗАКЛЮЧЕНИЕ                              29

БИБЛИОГРАФИЧЕСКИЙ СПИСОК                      30

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 Определение поверхности теплопередачи выпарных аппаратов      

 

 

2. Расчет концентраций упариваемого  раствора

 

Первое  приближение при расчете.       

Производительность  установки по выпариваемой воде определим  по формуле       

                                              (1)

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

W1 : W2 : W3 = 1,2 : 1,1 : 1,0

Тогда                          

                        (2)

                     (3)

               (4)

Далее рассчитаем концентрацию расходов по корпусам         

 

                             (5)

                          (6)

 

                           (7)

Расчет  выполнен верно.  

 

 

3. Определение температур кипения растворов

 

Общий перепад давлений в установке  равен 

 

                                         (8)

где РГ1 - давление греющего пара в первом корпусе, МПа         

  РБК - давление греющего пара в барометрическом конденсаторе, МПа    

Давление греющих паров в  корпусах              

                                                (9)

 

Расчет  выполнен верно.                    

 

  Таблица 1 Температуры и энтальпии греющих паров  

Давление, МПа 

Температура, °С

Энтальпия, кДж/кг

Теплота парообразования, кДж/кг

РГ1 = 6,797

tГ1 = 284

IГ1 = 2776

rГ1 = 1519

РГ2 = 4,558

tГ2 = 257

IГ2 = 2798

rГ2 = 1675

РГ3 = 2,319

tГ3 = 219

IГ3 = 2800

rГ3 = 1858

РБК = 0,0797

tБК = 93

IБК = 2666

rБК = 2274


 

Определим температуру  вторичных паров по формуле

 

                                               (10)

 

                                                (11)

                                  (12)

где ∆''' = 1°C - величина гидродинамической депрессии на корпус.

 

Сумма гидродинамических депрессий  составляет           

                                                      (13)

  Таблица 2 Давление и теплота парообразования вторичных паров     

Температура, °С

Давление, МПа 

Энтальпия, кДж/кг

Теплота парообразования, кДж/кг

tВП1 =258

РВП1 =4,558

IВП1 =2798

rВП1 =1675

tВП 2 =220

РВП2 =2,319

IВП2 =2800

rВП2 =1858

tВП3 =94

РВП3 =0,080

IВП3 =2666

rВП3 =2274


 

Поверхность теплопередачи первого  корпуса ориентировочная     

                               (14)

где q =80000  Вт/м2 - удельная тепловая нагрузка с естественной циркуляцией  

По ГОСТ 11987-81 для выпарного аппарата с естественной циркуляцией и вынесенной греющей камерой ближайшая поверхность -200  м2. /3/  

Выбираем греющие трубы размером Ø38 мм, толщина стенки δ = 2 мм, длина труб (высота) Н = 4 м.

Плотность раствора по корпусам определим  по справочнику /1/

 

 при этом
кг/м3

 при этом 
кг/м3

 при этом 
кг/м3

Давление в среднем слое кипятильных труб корпусов

                                        (15)

 

где Рвп - давление вторичных паров, МПа                

Н - высота кипятильных труб в аппарате, м         

ε - паронаполнение (объемная доля пара в кипящем растворе), м33

Принимаем ε = 0,5.

Тогда                      

 МПа

 МПа

 МПа

   Таблица 3 Температура кипения и теплота парообразования   

Давление, МПа 

Температура, °С

Энтальпия, кДж/кг

Теплота парообразования, кДж/кг

РСР1 =4,573

tСР1 =258

IСР1 =2797

rСР1 =1668

РСР2 =2,328

tСР2 =220

IСР2 =2800

rСР2 =2853

РСР3 =0,089

tСР3 =75

IСР3 =2670

rСР3 =2266


 

Определим гидростатическую депрессию по корпусам

 

              (16)

              (17)

              (18)

Сумма гидростатических дипрессий составляет

 

                       (19)

Температурную депрессию определим по формуле

 

         (20)

где ∆АТМ '- температурная депрессия при атмосферном давлении, град. /5/

 

Х1 = 97,0 % при этом ∆АТМ1' =12,5 °С

 

Х2 = 17,1 % при этом ∆АТМ2' =1,9 °С

Х3 = 9,7 % при этом ∆АТМ3' =0,8 °С

Тогда                          

 °С

 °С

 °С

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

 

∑∆' = ∆1' + ∆2' + ∆3'                                 (21)

∑∆'= 22,7 + 4,0 + 0,7 = 27,4 °С

Тогда температура кипения растворов  в корпусах составит       

 

tК1 = tГ2 + ∆1' + ∆1'' + ∆1'''                             (22)

tК1= 257 + 22,7 + 0 + 1= 280,73 °С

tК2 = tГ3 + ∆2' + ∆2'' + ∆2'''                  (23)

tК2= 219 + 4,0 + 0 + 1= 224,03 °С

tК3 = tБК + ∆3' + ∆3'' + ∆3'''                                                 (24)

tК3 = 93 + 0,7 + 2 + 1 = 96,69 °С                      

 

4. Расчет полезной разности температур

 

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

 

∆tП1 = tГ1 - tК1                                      (25)

∆tП1 = 284 - 280,73 = 3,27 °С

∆tП2 = tГ2 - tК2                                 (26)

∆tП2 = 257 - 224,03 = 32,97 °С

∆tП3 = tГ3 - tК3                                                 (27)

∆tП3 = 219 - 96,69 = 122,31 °С

Общая полезная разность температур                

 

∑∆tП = ∆tП1 + ∆tП2 + ∆tП3                                              (28)

∑∆tП = 3,27+32,97+122,31= 158,56 °С

Проверим общую полезную разность температур            

∑∆tП = tГ1 - tБК - (∑∆' + ∑∆'' + ∑∆''')                                        (29)

 

∑∆tП = 284 - 93 - (27,4 + 2,0 + 3) = 158,56 °С

Расчет  выполнен верно.        

 

 

5. Определение тепловых нагрузок

 

Уравнение тепловых балансов по корпусам выглядит так

Q1= D ∙ (IГ1 – i1) = 1,03 ∙ [(GH – W2 – W3) ∙ C2 ∙ (tK1 – tK2) + W1 ∙ (IВП1 – СВ ∙ tK1) + Q1КОНЦ]   (30)

Q2 = W1 ∙ (IГ2 - i2) = 1,03 ∙ [(GH - W3) ∙ C3 ∙ (tK2 - tК3) + W2 ∙ (IВП2 - СВ ∙ tK2) + Q2КОНЦ]        (31)

Q3 = W2 ∙ (IГ3 – i3) = 1,03 ∙ [GH ∙ CH ∙ (tK3 – tH) + W3 ∙ (IВП3 – СВ ∙ tK3) + Q3КОНЦ]      (32)

W = W1 + W2 +W3       (33)

W=25,0

 

где 1,03 - коэффициент, учитывающий потери в окружающую среду       

  С - теплоемкости растворов начальной концентрации и в корпусах     

  D - расход греющего пара в первом корпусе, кг/с            

  QКОНЦ - теплота концентрирования по корпусам (принимаем равным 0).  

  tН - температура кипения исходного раствора в первом корпусе      

 

Принимаем:

(IГ - i) = rГ

(IВП - СВ ∙ tK) = rВП

Окончательно уравнения выглядат так            

D ∙  rГ1 = 1,03 ∙ [(GH -W2 - W3) ∙ C2 ∙ (tK1 - tK2) + W1 ∙ rВП1]                       (34)

W1 ∙ rГ2 = 1,03 ∙ [(GH -W3) ∙ C3 ∙ (tK2 - tК3) + W2 ∙ rВП2]

W2 ∙ rГ3 = 1,03 ∙ [GH ∙ CH ∙ (tK3 - tH) + W3 ∙ rВП3]

W = W1 + W2 + W3 = 25,0

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

С = 4,19 ∙ (1 - Х)                   (33)

С2 = 4,19 ∙ (1 - 0,171) = 3,479 кДж/кг∙К

С3 = 4,19 ∙ (1 - 0,097) = 3,787 кДж/кг∙К

СН = 4,19 ∙ (1 - 0,070) = 3,902 кДж/кг∙К

Температура кипения исходного раствора в  первом корпусе      

tН = tВП1 + ∆H'                  (34)

tН = 258 + 4,63 = 262,63 °С

Температурную депрессию определим по формуле

 

                                        (35)

где ∆АТМ' - температурная дипрессия при атмосферном давлении, град.     

ХН = 7,0 % при этом ∆НАТМ' = 1,4 °С

Тогда  

 °С

Подставим известные величины в  уравнения и решим систему  уравнений относи-

тельно производительности по испаряемой воде

D ∙ 1519 = 1,03 ∙ [ (26,9 - W2 - W3) ∙ 3,479 ∙ (280,73 - 224,03) + W1 ∙ 1675]

W1 ∙ 1675 = 1,03 ∙ [ (26,9 - W3) ∙ 3,787 ∙ (224,03 - 96,69) + W2 ∙ 1858]

W2 ∙ 1858 = 1,03 ∙ [26,9 ∙ 3,902 ∙ (96,69 - 262,63) + W3 ∙ 2274]

                                                               W1 + W2 + W3 =25                                 (36)              

Упростим уравнения              

1519 ∙ D = 1725,3 ∙ W1 - 203,16 ∙ W2 - 203,16 ∙ W3 + 5474,1

1675 ∙ W1 = 1913,7 ∙ W2 - 496,7 ∙ W3 + 13384,4

1858 ∙ W2 = 2342,2 ∙ W3 - 17971,4

W1 + W2 + W3 = 25,0

Далее 

W3 = 0,7933 ∙ W2 + 7,6728

W1 = 0,9073 ∙ W2 + 5,7152

1519 ∙ D - 1725,3 ∙ W1 + 203,16 ∙ W2 + 203,16 ∙ W3 = 5474,1

0,9073 ∙ W2 + 5,7152 + W2 + 0,7933 ∙ W2 + 7,6728 = 25,0

Отсюда

W1 = 9,616 кг/с

W2 = 4,300 кг/с

W3 = 11,084 кг/с

D = 12,468 кг/с

Проверим правильность решения  уравнения

W1 + W2 + W3 = 25,0 кг/с                                                  (37)

9,616 + 4,300 + 11,084 = 25,0 кг/с

Расчет выполнен верно.                     

 

Определим тепловые нагрузки

Q1 = D ∙ rГ1                   (38)

Q1= 12,468 ∙ 1519 = 18939 кВт

Q2 = W1 ∙ rГ2                                          (39)

Q2= 9,616 ∙ 1675 = 16107 кВт

Q3 = 4,300 ∙ 1858 = 7989,1 кВт

Полученные параметры сведем в  таблицу 4.        

           Таблица 4 Параметры растворов и паров по корпусам         

Параметр

Корпус

1

2

3

Производительность по испаряемой воде W, кг/с

9,616

4,300

11,084

Концентрация растворов Х, %

97,0

17,09

9,7

Давление греющих паров РГ, МПа

6,797

4,558

2,319

Температура греющих паров tГ, °С

284

257

219

Температурные потери ∑∆, °С

23,73

5,03

3,69

Температура кипения раствора tК, °С

280,73

224,03

96,69

Полезная разность температур ∆tП, °С

3,27

32,97

122,31


 

 

6. Выбор конструкционного материала

 

Выбираем конструкционный материал, стойкий в среде кипящего раствора в интер-вале изменения концентрации от 5% до 40%. В этих условиях лучше всего подходит 

сталь марки Х17, имеющая скорость корозии менее 0,1 мм в год, коэффициент теплопро-

водности λ = 25,1 Вт/м∙К.                       

 

 

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

 

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

                                                       (40)

где   α1, α2 - коэффициенты теплоотдачи от конденсирующего пара к стенке и от      кипящего раствора к стенке, Вт/м2∙К                 

         λ - коэффициент теплопроводности, Вт/м∙К     

 

Суммарное термическое сопротивление  определим по формуле       

 м2∙К/Вт                (41) 

 

где   δСТСТ, δНН - термическое сопротивление стенки и накипи, сопротивление загряз-

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

        δСТ = 0,002 м - толщина стенки трубы                

        δН = 0,0005 м - толщина накипи                   

        λН = 2,42 Вт/м∙К - теплопроводность накипи

 

 

           


 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

раствору  через многослойную стенку: 1 - пар; 2 – конденсат; 3 – стенка; 4 – накипь;   

5 –  кипящий раствор.                         

 

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

Формуле

                     


(42)  

      (42)

 

где   r1 - теплота конденсации греющего пара, Дж/кг                     

        ρж1, λж1, μж1 - соответствено плтность (кг/м3), теплопроводность (Вт/м∙К), вязкость

конденсата  при средней температуре пленки (Па∙с) /1/           

 

Средняя температура пленки

                                                        43)

где   ∆t1 - разность температур конденсата пара и стенки, град.          

 

Расчет  α1 ведут методом последовательных приближений.         

 

1-е  приближение                       

Принимаем ∆t1 = 0,5 °С и проверим правильность первого приближения по равен-ству удельных тепловых нагрузок.      

 

Средняя температура пленки

 °С

 °С

 °С                  

              Таблица 5 Физические свойства конденсата при средней температуре     

Параметр

Корпус

1

2

3

Теплота конденсации греющего пара rГ, кДж/кг

1512

1668

1858

Плотность конденсата ρж, кг/м3

743,2

786,9

840,3

Теплопроводность конденсата λж, Вт/м∙К

0,572

0,608

0,648

Вязкость конденсата μж, Па∙с ∙ 10-3

0,092

0,103

0,122


 

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

        


 

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

                           (44)     

где q - удельная тепловая нагрузка, Вт/м2                 

  ∆tСТ - перепад температур на стенке, град               

  ∆t2 - разность между температурой стенки со стороны раствора и температурой  

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

 

Перепад температур на стенке определим  по формуле

                                                (45)

Тогда разность температур составляет

∆t2 = ∆tП1 - ∆tСТ - ∆t1                       (46)

 

∆t2= 3,27 - 1,58 - 0,5 =1,20 °С

Коэффициент теплопередачи от стенки к кипящему раствору определим по формуле

                         (47)

где ρж, ρп, ρ0 - соответственно плотность жидкости, пара и пара при абсолютном давле-

нии р = 1 атм., кг/м3                       

  ρ0 = 0,579 кг/м3 - плотность пара при абсолютном давлении        

  σ - поверхностное натяжение, Н/м                  

 

 

  μ - вязкость раствора, Па∙с                     

 

Физические свойства раствора в  условиях кипения запишем в таблицу 6

              Таблица 6 Физические свойства кипящих растворов и их паров       

Параметр

Корпус

1

2

3

Плотность раствора ρ, кг/м3

1291

960

903

Вязкость раствора μ, Па∙с ∙ 10-3

0,684

0,256

0,194

Теплоемкость раствора с, Дж/кг∙К

2802

3857

4180

Теплопроводность раствора λ, Вт/м∙К

0,641

0,668

0,677

Поверхностное натяжение σ, Н/м

0,11

0,079

0,077

Теплота парообразования rВ, кДж/кг

1675

1858

2274

Плотность пара ρп, кг/м3

22,7

11,5

0,48


 

Подставим известные величины в  уравнение

Коэффициент теплопередачи от стенки к кипящему раствору составляет

                                         (48)

Проверим правильность первого  приближения

q1 = α1 ∙ ∆t1

q1= 11013 ∙ 0,5 = 5506,6 Вт/м2

q2 = α2 ∙ ∆t2

q2= 2183,5 ∙ 1,20 = 2616,4 Вт/м2

q1 ≠ q2

 

2-е  приближение                      

Принимаем ∆t1 = 0,26 °С и проверим правильность первого приближения по равен-ству удельных тепловых нагрузок.     

 

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


 

 (49)

 

 

 

Перепад температур на стенке определим  по формуле

∆tСТ = 12969 ∙ 0,26 ∙ 0,00029 = 0,97 °C

Тогда разность температур составляет

∆t2 = 3,27 - 0,97 - 0,26 = 2,0 °С

Коэффициент теплопередачи от стенки к кипящему раствору составляет

α2 = 12,434 ∙ (12969,2 ∙ 0,26)0,6 = 1626,8 Вт/м∙К

Проверим правильность второго  приближения

q1 = α1 ∙ ∆t1

q1= 12969 ∙ 0,26 = 3372 Вт/м2

q2 = α2 ∙ ∆t2

q2= 1626,8 ∙ 2,0 = 3334,1 Вт/м2

q1 ≈ q2

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

 Вт/м2∙К

 

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

 

1-е  приближение                      

Принимаем ∆t1 = 2 °С и проверим правильность первого приближения по равенству удельных тепловых нагрузок.

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


                                                                       

                                                                                                                   Вт/м∙К    

 

 

Перепад температур на стенке определим  по формуле

Тогда разность температур составляет

∆t2 = ∆tП2 - ∆tСТ - ∆t1

∆t2 = 32,97 - 4,79 – 2 = 26,19 °С

Подставим известные величины в  уравнение

Коэффициент теплопередачи от стенки к кипящему раствору составляет

α2 = 14,694 ∙ (α1 ∙ ∆t1)0,6

α2= 14,694 ∙ (8357,7 ∙ 2)0,6 = 5023,5 Вт/м∙К

Проверим правильность первого  приближения            

q1 = α1 ∙ ∆t1

q1= 8357,7 ∙ 2 = 16715 Вт/м2

q2 = α2 ∙ ∆t2

q2 = 5023,5 ∙ 26,19 = 131556 Вт/м2

q1 ≠ q2

 

2-е  приближение                      

Принимаем ∆t1 = 10,7 °С и проверим правильность первого приближения по ра-

венству удельных тепловых нагрузок.                  

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

                                                                                           


 

 

Перепад температур на стенке определим  по формуле 

∆tСТ = 5495,4 ∙ 10,7 ∙ 0,00029 = 16,8 °C

Тогда разность температур составляет

∆t2 = 32,97 - 16,8 - 10,7 = 5,44 °С

Коэффициент теплопередачи от стенки к кипящему раствору составляет

 

α2 = 14,694 ∙ (5495,4 ∙ 10,7)0,6 = 10685 Вт/м∙К

Проверим правильность второго  приближения

q1 = α1 ∙ ∆t1

q1= 5495,4 ∙ 10,7 = 58800 Вт/м2

q2 = α2 ∙ ∆t2

q2= 10685 ∙ 5,44 = 58119 Вт/м2

q1 ≈ q2

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

Вт/м2∙К

 

Коэффициент теплопередачи третьего корпуса 

 

1-е  приближение                      

Принимаем ∆t1 = 2 °С и проверим правильность первого приближения по равенству удельных тепловых нагрузок.       

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


 

                                              Вт/м∙К    

 

Перепад температур на стенке определим  по формуле

Тогда разность температур составляет

∆t2 = ∆tП2 - ∆tСТ - ∆t1

∆t2 = 32,97 - 5,11 – 2 = 25,87 °С

Подставим известные величины в  уравнение

Коэффициент теплопередачи от стенки к кипящему раствору составляет

α2 = 11,408 ∙ (α1 ∙ ∆t1)0,6

α2= 11,408 ∙ (8921,3 ∙ 2)0,6 = 4055,9 Вт/м∙К

Проверим правильность первого  приближения 

q1 = α1 ∙ ∆t1

q1= 8921,3 ∙ 2 = 17843 Вт/м2

q2 = α2 ∙ ∆t2

q2= 4055,9 ∙ 25,87 = 104906 Вт/м2

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