Теплонспользующие установки промышленных предприятий
Расчет температурного режима работы выпарной установки
В технических заданиях на разработку проектов
Многокорпусных выпарных установок приводятся параметры греющего пара первого корпуса и вторичного пара последнего корпуса, начальная и конечная концентрации раствора, производительность по выпаренной воде или выделившейся соли.
Все остальные величины — концентрация раствора в корпусах, температурные депрессии, полезные разности температур, давление, температуры кипения — определяются при расчете. Расчет. многокорпусной выпарной установки (рис. 9.5) ведется методом последовательных приближений. Предварительно определяется общее количество выпаренной воды по всей установке. Дальнейший расчет выполняется в такой последовательности.
Первое приближение:
1. Количество выпаренной воды по корпусам распределяется равномерно. В общем виде, т. е. если имеет место отбор экстрапара, количество выпаренной воды в корпусах без отбора экстрапара получим из уравнения
(9.33) |
2 оВВП£ -ЪЕ1 (=і і=і
В остальных корпусах
(9.34) |
Бвып; — СіВЬШі + Е(.
+ 1 • |
2. Температура вторичного пара
*Г. П)
^В. П,-_| -- ^В. П; + '
3. Концентрацию солей по корпусам определяем при упаривании раствора без выделения твердой фазы по зависимости
По корпусам установки
"К“1)
(9.35)
Тадліянж дан *ае физика - цинические сШст(д
Исходные данные
Расчет
0 бып о&и
/ итерация Г “ |
Ьі = |
Расчет І-Іт, & --- І.-------- Расчет 1/ц і, Ъг, Вабср ^ і Сі, Я <•, /-*1 X Расчет <х, і, лсі, К, м, Ун Расчет I7i. it и а Л {-кипі} tr. nl) X |
(9.36)
«о
В случае упаривания многокомпонентных растворов — в следующей последовательности:
А) находим критерий растворимости (отношение количества содержащейся в растворе воды к количеству некристаллизующейся части его), т. е.
І=і |
Рис. 9.5. Блок-схема расчета многокорпусной выпарной установки |
100
(9.37
Б-Ь,.
Юи
Б) по температуре вторичного пара и параметру Й,- определяем концентрацию некристаллизую - щейся части раствора;
В) по полученной концентрации некристаллизующейся части раствора находим содержание остальных компонентов согласно таблицам или кривым совместной растворимости систем.
3. Количество выделившейся соли и количество поступающего раствора определяют по зависимостям (9.12), (9.14), (9.17), (9.18) Так как раствор переходит из корпуса в корпус, конечная концентрация раствора в первом корпусе Ьс и Ь{ равна начальной концентрации раствора во втором корпусе /7(±1 и Ь±] и т. д.
Расчет а4-,л/, А, Ж, А, *1П, У і
Расчет 'Вс, Ві,
І у ^ і.____
Расчет подогрева - ■ телеа
X
Расчет 6[
X
2 "и 5алее итераций
Г
Результаты расчета ^■кипі, ^ й П і і ЕііВгШІ, Ріі
5. Температурную депрессию раствора выбираем по справочным таблицам для конечных концентраций раствора в выпарных аппаратах. Однако в справочной литературе температурная депрессия приведена для атмосферного давления, а так как корпуса установки обычно работают под давлением, отличающимся от атмосферного, необходим перерасчет температурной депрессии по одному из трех методов.
По закону Бабо
—р—? == согЫ, (9.38)
^ Н
Где Рн — давление пара чистого растворителя при температуре кипения раствора; Рр—давление пара раствора. С помощью зависимости (9.38) по известной температуре кипения раствора, заданной концентрации при каком-нибудь одном давлении определим температуру кипения его при любом другом давлении. Эмпирический закон Бабо справедлив только для разбавленных растворов. Для концентрированных растворов вводится поправка, найденная В. Н. Стабниковым (табл. 9.1).
9.1. Поправочный коэффициент
|
По правилу Дюринга, которое выражает линейность химикотехнических функций,
1 ~ h
; = const, (9.39)
^2 h
Где t—1[, h —12 — разность температур кипения какой-либо
Жидкости при двух произвольно взятых давлениях и какой-либо
Другой жидкости при тех же давлениях. В качестве второй жидкости, для которой обычно известны температуры кипения при различных давлениях, берется вода.
По формуле И. А. Тищенко, выведенной из эмбулиоскопиче- ского закона Рауля и закона Бабо [81!,
/273 + ( n. f
О,. = 0,0162&; ^. (9.40)
/в. п£-
Она справедлива только для водных растворов, имеющих малое значение температурной депрессии (до 15—17 °С) .
6. Температура кипения по корпусам выпарной установки
^КИП; = ^В. П( (9*41)
7. Температура греющего пара по корпусам
^Г. П I — ^В. П[-_________ 1
8. Температуру пленки конденсата греющего пара и перепад температуры в ней представим в виде
(пл[ — 0,25 ^3/г. п,- ^кип(
------ 0,5 |/г П£ /кип; Я )•
( — и. и п< <КИП( Ч хс
9. Коэффициенты теплоотдачи со стороны пара, раствора при принудительной и естественной циркуляции, коэффициент теплопередачи и скорость циркуляции, гидростатическая депрессия, тепловая нагрузка по корпусам определяются по зависимостям.
10. Суммарная полезная разность температур многокорпусной установки
TOC o "1-5" h z "к 1
Е д/п< = /г. п, - ?в. пп - Е (Ь + М - Е Ч (9.4?)
»'=1 к 1=1 С=1
Обычно принимают Е К/ — «к — I )•
4=1 1
11. Распределение общей полезной разности температур между корпусами находим с учетом равенства поверхностей нагрева корпусов. При этом условии обеспечивается взаимозаменяемость выпарных аппаратов, упрощается и удешевляется их изготовление и эксплуатация!
О "к
V —
К к.
Далее составляется таблица, в которую вносятся полученные значения для всех корпусов и по которой определяется температурный режим работы выпарной установки (табл. 9.2).
Показатели
Температура греющего пара * °С Теплосодержание греющего пара Jmi, кДж/кг Температура кипения раствора °С Температурная депрессия vii, °С Гидростатическая депрессия vгi, СС Потери температуры в паропроводе между корпусами ипі, °С
Температура вторичного пара 1Ъ пі, °С Теплосодержание вторичного пара У1 п(, кДж/кг Полезная разность температур Д<п£, °С Скрытая теплота парообразования гъ ш - кДж/кг Энтальпия конденсата греющего пара У2(-, кДж/кг
Так как в техническом задании на проектирование давления вторичного пара, а следовательно, и температура в последнем корпусе известны, то расчет темпера-турного режима работы установки выполняется, начиная с последнего корпуса:
Температура кипения /КИпЛк =
Температура греющего пара ^.п„к = 4ип„к + Ып + 9^;
Температура вторичного пара предпоследнего корпуса
После установления предварительного режима работы выпарной установки определяют коэффициенты самоиспарения и испарения хс и I/,-, а затем количество выпаренной воды и расход греющего пара с учетом этих коэффициентов по зависимостям (9.31), (9.32).
Далее уточняется концентрация раствора по корпусам и производится расчет второго приближения. Во втором приближении рассчитываются депрессии, коэффициенты теплопередачи, скорости циркуляции, тепловые нагрузки и полезные разности температур. по к'орпусам для режима работы, полученного в первом приближении.
По окончании расчетов составляется новый режим работы выпарной установки. Точность расчета контролируется отношением полученных производительностей и параметров работы установки первого, второго и т. д. приближений. Расхождение не должно превышать 5%.
После получения требуемой точности определяют поверхности нагрева корпусов:
<3| <?2
= = и т' д'
Полученное значение поверхности теплообмена увеличивается на 15—20% и округляется до ближайшего большего значения по ГОСТ 11987—81 «Аппараты выпарные трубчатые стальные» или каталогу ЦИНТИхимнефтемаша «Выпарные трубчатые стальные аппараты общего назначения».
Пример. Рассчитать 4-корпусную выпарную установку для электролитических щелоков при следующих исходных данных: схема выпарной установки прямоточная, исходные щелока не подогреваются. Количество корпусов п =-4. Концентрация растворов: исходного по некристаллизующейся части Ьа =
= 9,1 %, по кристаллизующейся части Ь'0 = 17,5 %, по воде *в = 73,4 %; упаренного—по некристаллизующейся части Ьк = 50 %, по кристаллизующейся части *4=2,7°/
Теплота кристаллизации лкр = — 87,9 кДж/кг. Температура греющего пара первого корпуса ег п| = 179 °С, вторичного пара четвертого корпуса <в. п4 = 46 °С> исходного раствора ^ = 70°С, конденсата после подогревателя раствора ^ п = 95 °С.
Размеры грекЭщих трубок — длина I = 7,0 м, внутренний диаметр сУ = 0,034 м, толщина 5СТ = 0,002 м. Коэффициент теплопроводности материала стенки трубки Хст =16,2 Вт/ (м - К).
Свойства соли — плотность? т = 2150 кг/м-1, теплоемкость ст = 0,87 кДж/ (кГ'К), молекулярная масса /Иц = 58, число молекул кристаллизационной воды N = 0. Скорость движения раствора по греющим трубкам и>р = 2,0 м/с.
Температурная депрессия определяется для первого — третьего корпусов по формуле Тищенко И. А., для четвертого—по экспериментальным данным 9,
= 44 °С.
Количество отбираемого экстрапара определяется расчетом, первоначально принимаем Оэ = 20 000 кг/ч.
Коэффициент тепловых потерь Дп1 = &п2 = Дга3 = ДпА = 0,97.
Содержание твердой фазы в циркулирующей суспензии: рт) =0; (Зт2 = 0,05; Ртз = 0,1; Эт4 = 0,05.
Результаты материального и теплового расчета выпарной установки приведены в табл. 9.3.
Результаты |
По корпусам |
|||
Расчетная формула |
— |
|||
Или способ определения |
1 | 2 |
3 |
4 |
19,45 |
20.83 |
Общее количество гаыпа-
"ч
Репной воды 2 Овы,, кг/с
Общее количество выдели
Лившейся СОЛИ
. — 1
Кг/с
Количество упаренного щелока 5К, кг/с
Материальный расчет По формуле (9,16) 74,16
По формуле (9.14)
Тепловой расчет (первое приближение)
|
Расчетная формула или способ определения
3,64 1193,6 179- 175,25 |
3,55 1230 145.5 141.5 |
3,43 1299 112 108,7 |
3,32 1491,5 78,5 |
~Р19 |
*ш,( = 0.25 3<Г>П|.+ |
81,9 |
КДж/(кг К) плотность рп1-, кг/.м3 Температура греющего пара /г. п£> °с Температура пленки конденсата гПЛ I, °С |
Теплоемкость |
9 = 17 500 Вт/м2—задаемся для пер вой итерации
-6,32 |
6,63 |
А^ш.( = °-5 К. п£ — |
Перепад температур в пленке конденсата Д/Пл ( °С |
|
|
|||
^кип I Я
198.5 2018,3 5670 |
194,2 2124.8 5547 1128,3 8588 9288 2430 3.4 1274.9 |
170 2310,7 5180 1128,3 5403 6779 2158 3,17 1524 |
182,9 2224 5520 1746 6705 8187 2344 3,19 1382,3 |
II. — Г„ „ |
0,25 |
10 200 10 200 2 513 3.64 1193,6 |
Коэффициент, учитывающий влияние твердой фазы на теплообмен ус Теплоотдача со стороны щелочи: Без учета влияния твердой фазы а21) Вт/м2- К С учетом твердой фазы аа-, Вт/м* • К Коэффициент теплопередачи К;, Вт/м2 . К Теплоемкость суспензии сс, кДж/кг. К Плотность суспензии р,-, кг/м3 |
Коэффициент конденсации В. Теплота парообразования греющего пара гг П1 кДж/кг Теплоотдача со стороны пара аи, Вт/(м2-К) |
Д*пл4 < 0, поэтому Д^ПЛ4 = ^плЗ [70] [70] А - 2,045,- ;д/ Пл I По формуле (2.12) По формуле (2.4) По формуле (2,12) По формуле (2.13) По формуле (2.13) |
Пок |
Расчетная формула или способ определения
2/к Рг. т~ |
Общий Напор |
0=0 Вы По формуле (9.43) |
По формуле (9.21) |
По формуле (9.22) По формуле (9,24) По формуле (9.25) По формуле (9,32) |
Гидростатическая депрессия иг., СС
Температурный
°с
Теплопгя нагрузка (?,> Вт
Температурный и'пор по корпусам I, °С Температурный режим работы установки, °С. температура кипе - ния - (кип ;
Температура греющего пара
Температура вторич - иого пара 1Я пГ
Параметры пара, кДж/кг: теплосодержание: греющего пара ш- вторичного пара
Конденсата греюще - го пара Укон Коэффициенты:
Испарения пара самоиспарения конденсата а(
Самоиспарения щелочи
О
Концентрации
■'к
V Х(
1
<=!
"к
Расход греюшего пара <7,, кг/с
1,82 |
1,85 |
1,52 |
0.96 |
32,26 |
— |
||
48-106 8,75 |
38-10б 7,15 |
39,6 X хю1 7,7 |
41•10® 8.7 |
167,9 |
142,5 |
117,7 |
90 |
179 |
151,5 |
127 |
99.6 |
152,5 |
128 |
100.6 |
46 |
2790 2750 |
2752 2717 |
2718 2678 |
2678.7 2585 |
758,7 |
6.39 |
533,2 |
417.6 |
0,944 0 |
0,955 0 |
0,96 0.465 |
0,989 0,506 |
0.172 |
0,047 |
0,042 |
—0,032 |
—0,26 |
—0,222 |
—0,172 |
-0.136 |
—0,01 |
—0,01 |
—0,01 |
|
0,0185 |
—0,006 |
—0,004 |
—0,0018 |
0.916 |
0,848 |
0,83 |
0.84 |
18,7 |
15,К |
11.9 |
10,4 |
38.1 |
— |
— |
103а, РЛР^ + Ы1 По формуле (9.42) |
*г. пС ^кип I игч + 4/п I ^В. П!—1 = ^Г. П( + ип/ [70] |
Показатель |
Расчетная формула или способ определен, ИЯ |
Результаты |
По корпусам |
||
- |
2 |
3 |
4 |
||
• |
|||||
Количество выпаренной воды по корпусам О0ЫП м кг/с |
<3ПЫП1 = ЛА +У, |
16,2 |
16,55 |
19,8 |
21.6 |
Удельная тепловая нагрузка qi, Вт/м2 |
?£ = М'ш |
21 986 |
17 374 |
18 050 |
18 775 |
Тепловой расчет (второе приближение)
Расчет начинается с определения Я1 по данным, полученным в перво» приближении. Тепловая нагрузка определяется в зависимости от схемы течения пара и раствора по формулам (9.29) — (9.31). Количество приближений определяется условием
(Зпып / + 1 Овып! ^ о,05 а — количество приближений), бнып /
После получения требуемой точности находим поверхность теплообмена Р. = С}./ (ЬЛ1п<), где Т7, .-= 2194,4 м2; = 2193,4 м2; = 2195 м2; Р^=-
= 2192,9 ма. Принимаем = р2 = = Р4 = 2200 м2.