Расчет многокорпусной выпарной установки
Изм.
Лист
№ докум.
Подпись
Дата
Лист
ЭА52.0546.0000 ПЗ
ФГОУ ВПО Башкирский государственный аграрный университет
Факультет: Энергетический
Кафедра: Теплотехники и
энергообеспечения п/п
Специальность: Энергообеспечение п\п
Форма обучения: заочная
Курс, группа: 3,1
Никулина Ксения Владимировна
Курсовая работа
по дисциплине «Теоретические основы теплотехники»
РАСЧЕТ МНОГОКОРПУСНОЙ ВЫПАРНОЙ УСТАНОВКИ
(название работы)
Руководитель:к.т.н. доцент Юхин Д.П
Оценка при защите:
____________________
____________________
(подпись)
«____» __________ 2013 г
Уфа 2013
ОГЛАВЛЕНИЕ
ВВЕДЕНИЕ |
1 ОПИСАНИЕ ТЕХНОЛОГИЧЕСКОЙ СХЕМЫ УСТАНОВКИ 5 |
1.1 Принцип работы двухкорпусной выпарной
установки |
2 ОПРЕДЕЛЕНИЕ ПОВЕРХНОСТИ
ТЕПЛОПЕРЕДАЧИ 2.1 Определение концентрации
упариваемого раствора 2.2 Определение температуры
кипения растворов 2.3 Определение полезной
разности температур 2.4 Определение тепловых нагрузок 2.5 Выбор конструкционного материала 14 2.6 Расчет коэффициентов
теплопередачи 2.7 Распределение полезной разности температур 17 |
3 ОПРЕДЕЛЕНИЕ ТОЛЩИНЫ ТЕПЛОВОЙ ИЗОЛЯЦИИ 18 |
4 РАСЧЕТ БАРОМЕТРИЧЕСКОГО КОНДЕНСАТОРА 19 4.1 Определение расхода 4.2 Определение диаметра 4.3 Определение высоты 4.4 Расчет производительности |
ЗАКЛЮЧЕНИЕ |
БИБЛИОГРАФИЧЕСКИЙ СПИСОК |
ВВЕДЕНИЕ
Выпарные аппараты предназначены для концентрирования жидких растворов практически нелетучих веществ путем частичного удаления растворителя испарением при кипении жидкости. В процессе выпаривания растворитель удаляется из всего объема раствора, в то время как при температурах ниже температур кипения испарение происходит только с поверхности жидкости.
Процесс выпаривания - энергоемкий процесс, особенно если теплота испарения. Как например у воды. Поэтому составляющая на энергозатраты при выпаривании может быть весьма существенной составляющей в себестоимости производства того или иного продукта. Одним из наиболее эффективных способов снижения энергопотребления является применение выпарных батарей- многокорпусных выпарных установок.
Веществом, подлежащим концентрированию в водном растворе, является NaNO3.
Применяется как удобрение; в стекольной, металлообрабатывающей промышленности; для получения взрывчатых веществ, ракетного топлива и пиротехнических смесей для придания огню жёлтого цвета. Получается из природных залежей выщелачиванием горячей водой и кристаллизацией; абсорбцией раствором соды окислов азота; обменным разложением кальциевой или аммиачной селитры с сульфатом, хлоридом или карбонатом натрия.
1 ОПИСАНИЕ ТЕХНОЛОГИЧЕСКОЙ СХЕМЫ УСТАНОВКИ
В химической и смежной с ней отраслях промышленности жидкие смеси, концентрирование которых осуществляется выпариванием, отличаются большим разнообразием как физических параметров (вязкость, плотность, температуря кипения, величина критического теплового потока и др.), так и других характеристик (кристаллизующиеся, пенящиеся, нетермостойкие растворы и др.). Свойства смесей определяют основные требования к условиям проведения процесса (вакуум-выпаривание, прямо- и противоточные, одно- и многокорпусные выпарные установки), а также к конструкциям выпарных аппаратов.
Такое разнообразие требований вызывает определенные сложности при правильном выборе схемы выпарной установки, типа аппарата, числа ступеней к многокорпусной выпарной установке. В общем случае такой выбор является задачей оптимального поиска и выполняется технико-экономическим сравнением различных вариантов с использованием ЭВМ.
Ниже приведен типовой расчет двухкорпусной установки, состоящей из выпарных аппаратов с естественной циркуляцией (с соосной камерой).
- Принцип работы двухкорпусной выпарной установки
Исходный разбавленный раствор из промежуточной емкости центробежным насосом подается в теплообменник, где прогревается до температуры, близкой к температуре кипения, а затем – в первый корпус выпарной установки. Предварительный подогрев раствора повышает интенсивность кипения в выпарном аппарате.
Первый корпус обогревается свежим водяным паром. Вторичный пар, образующийся при концентрировании раствора в первом корпусе, направляется в качестве греющего во второй корпус. Сюда же поступает частично сконцентрированный раствор из 1-го корпуса.
Самопроизвольный
переток раствора и вторичного пара
в следующие корпуса возможен
благодаря общему перепаду давлений,
возникающему в результате создания
вакуума конденсацией вторичного пара
последнего корпуса в барометрическом
конденсаторе смешения, где заданное
давление поддерживается подачей охлаждающей
воды и отсосом неконденсирующихся
газов вакуум-насосом. Смесь охлаждающейся
воды и конденсата выводится из конденсатора
при помощи барометрической трубы
с гидрозатвором. Образующийся во втором
корпусе концентрированный
Конденсат
греющих паров из выпарных аппаратов
выводится с помощью
2 ОПРЕДЕЛЕНИЕ ПОВЕРХНОСТИ ТЕПЛОПЕРЕДАЧИ ВЫПАРНЫХ АППАРАТОВ
Поверхность
теплопередачи каждого корпуса
выпарной установки определяется по
основному уравнению
F = Q/(K·Δtп)
где Q – тепловая нагрузка, кВт;
K – коэффициент теплопередачи, Вт/(м2·K);
Δtп – полезная разность температур, °С.
Для определения тепловых нагрузок Q, коэффициентов теплопередачи К и полезных разностей температур Δtп необходимо знать распределение упариваемой воды, концентраций растворов и их температур кипения по корпусам. Эти величины находятся методом последовательных приближений.
Производительность
установки по выпариваемой воде определяется
из уравнения материального
W = Gн(1
– xн/xк)
где Gн – производительность установки по исходному раствору, кг/с;
xн, xк – массовые концентрации вещества в исходном и упаренном растворе соответственно, %.
W = 12,77·(1 – 6/46) = 11,11 кг/с.
2.1
Определение концентрации
Распределение концентраций раствора по корпусам установки зависит от соотношения нагрузок по выпариваемой воде в каждом аппарате. В первом приближении на основании практических данных принимают, что производительность по выпариваемой воде распределяется между корпусами в соотношении:
w1: w2= 1,0: 1,1
Тогда,
w1 = 1,0W/(1,0 + 1,1) = 11,11/2,1 = 4,37кг/с,
w2 = 1,1W/(1,0 + 1,1) = 1,1∙11,11/2,1 = 5,29 кг/с,
Далее рассчитываются концентрации растворов в корпусах:
x1 = Gнxн/(Gн - w1) = 11,11∙6/(11,11-4,37) =9,9%,
x2 = Gнxн/(Gн - w1 - w2) =11,11∙6/(11,11-4,37-5,92) = 45,9%,
Концентрация раствора в последнем корпусе x2 соответствует заданной концентрации упаренного раствора xк.
2.2
Определение температуры
Общий перепад давлений в установке равен:
ΔPоб
= Pг1 – Pбк
где Pг1 – давление греющего пара, МПа;
Pбк – абсолютное давление в барометрическом конденсаторе, МПа.
ΔPоб =0,546-0,0546=0,4914 МПа.
В первом
приближении общий перепад
Pг1 = 0,546 МПа,
Pг2 = Pг1 - ΔPоб/2 = 0,546 – 0,4914/2 = 0,3003 МПа,
Давление
пара в барометрическом
Pбк = Pг3 - ΔPоб/2 = 0,3003 – 0,4914 = 0,0546 МПа
что соответствует заданному значению Pбк.
По давлениям паров находим их температуры и энтальпии [1]:
Таблица 2.1 Давление паров по температуре и энтальпии
P, МПа |
t, 0C |
I, кДж/кг |
i, кДж/кг |
Pг1 =0,546 |
tг1 = 155,18 |
I1 = 2752 |
i1=654,66 |
Pг2 = 0,3003 |
tг2 = 133,56 |
I2 = 2724,9 |
i2=561,6 |
Pбк = 0,0546 |
tбк = 83,53 |
Iбк = 2648,9 |
iбк=349,75 |
При определении
температуры кипения растворов
в аппаратах исходят из следующих
допущений. Распределение концентраций
раствора в выпарном аппарате с интенсивной
циркуляцией практически
Изменение температуры кипения по высоте кипятильных труб происходит вследствие изменения гидростатического давления столба жидкости.
Температуру
кипения раствора в корпусе принимают
соответствующей температуре
Гидродинамическая
депрессия обусловлена потерей
пара на преодоление гидравлических
сопротивлений трубопроводов
tвп1 = tг2 + Δ1/// = 133,56+1=134,56,
tвп2 = tбк + Δ2/// =83,53+1=84,53.
Сумма гидродинамических депрессий
ΣΔ/// = Δ1/// + Δ2/// = 1+1 = 2 0С.
По температурам вторичных паров определим их давление.
Таблица 2.2 Давление вторичных паров
Температура, °С |
Давление, МПа |
tвп1 =134,56 |
Pвп1 =0,309 |
tвп2 =84,53 |
Pвп2 = 0,057 |
Гидростатическая
депрессия обусловлена
Рср =
Pвп + ρ·g·H (1- ε)/2,
где Н- высота кипятильных труб в аппарате, м;
ρ – плотность кипящего раствора, кг/м3;
ε – паронаполнение (объемная доля пара в кипящем растворе), м3/м3.
Для выбора значения Н необходимо ориентировочно оценить поверхность теплопередачи выпарного аппарата FОР. При кипении водных растворов можно принять удельную тепловую нагрузку аппарата с естественной циркуляцией q = 20000 – 50000 Вт/м2. Примем q = 40000 Вт/м2. Тогда поверхность теплопередачи первого корпуса ориентировочно равна:
Fор
= Q/q = ω1∙r1/q,
где r1 – теплота парообразования вторичного пара, Дж/кг.
Fор = Q/q = ω1∙r1/q = 4,37∙2065∙103 / 40000 = 225,6 м2.
По ГОСТ 11987 – 81 трубчатые аппараты с естественной циркуляцией и соосной греющей камерой состоят из кипятильных труб, высотой 4 и 5 м при диаметре dн = 38 мм и толщине стенки δст = 2 мм. Примем высоту кипятильных труб H = 4 м.
При пузырьковом (ядерном) режиме кипения паронаполнение ε = 0,4 – 0,6.Примем ε = 0,5.
Плотность водных растворов, в том числе NaNO3, при температуре 20 0С и соответствующих концентрациях в корпусах равна:
ρ1 = 1075,7кг/м3, ρ2 = 1393,2 кг/м3.
При определении плотности растворов в корпусах пренебрегаем изменением ее с повышением температуры от 20 0С до температуры кипения ввиду малого значения коэффициента объемного расширения и ориентировочно принятого значения ε.
Давления в среднем слое кипятильных труб корпусов (в Па) равны:
Р1 ср = 30,9∙104 + 1075,7∙9,8∙4∙(1 – 0,5)/2 = 31,95∙104,
Р2 ср =5,7∙104 + 1393,2·9,8∙4∙(1 – 0,5)/2 = 7,07∙104,
Этим
давлениям соответствуют
Таблица
2.3 Зависимость давления от температуры
P, МПа |
t, 0C |
r, кДж/кг |
P1ср = 0,3195 |
t1ср =135,7 |
rвп1 = 2157 |
P2ср = 0,0707 |
t2ср =90,2 |
rвп2 = 2282 |
Определим гидростатическую депрессию по корпусам (в 0C):
Δ1// = t1ср - tвп1 =135,7– 134,56= 1,14 0С,
Δ2// = t2ср - tвп2 = 90,2– 84.53 = 5,67 0С,
Сумма гидростатических депрессий
ΣΔ// = Δ1// + Δ2// = 1,14 + 5,67=6,81 0С.
Температурную депрессию Δ/ определим по уравнению
Δ/
= 1,62∙10-2∙ Δатм/ ∙Т2/
r вп
где Т – температура паров в среднем слое кипятильных труб, К;
Δатм/ - температурная депрессия при атмосферном давлении.
Находим значение Δ/ по корпусам (в 0C):
Δ/1=
1,62∙10-2 ∙
(135,7+ 273)2∙ 1,14 / 2157=1,43,
Δ/2= 1,62∙10-2 ∙ (90,2+ 273)2∙ 5,67 / 2282= 5,31,
Сумма температурных депрессий
ΣΔ/ = Δ1/ + Δ2/ =1,43+5,31=6,74 0С.
Температуры кипения растворов в корпусах равны (в 0C)
tк
= tг + Δ/ +
Δ//
tк1 = tг2 + Δ/1 + Δ//1+Δ///1 = 133,56+1,43+1,14+1=137,13,
tк2 = tбк + Δ/2 + Δ//2+Δ///2 =83,53+5,31+5,67+1=95,51.
В аппаратах с соосной греющей камерой и естественной циркуляцией обычно достигаются скорости раствора v = 0,6 – 0,8 м/с [5,7,8]. Для этих аппаратов масса циркулирующего раствора равна:
M = u∙S∙ρ, (2.8)
М=0,7∙1044∙0,54=394,6
где S- сечение потока в аппарате (м2), рассчитываемая по формуле:
S = Fор ∙dвн/4∙H, (2.9)
S= (225,6 ∙0,038)/(4·4)=0,54 м2
где dвн – внутренний диаметр труб, м;
Н – принятая высота труб, м.
Таким образом, перегрев раствора в j-м аппарате Dtперj равен:
Dtперj = [ωо∙ (Iвп j – св·tкj) - Gнj ∙снj ∙ (tкj-1 - tкj)] / M ∙снj (2.10)
2.3 Определение полезной разности температур
Общая полезная разность температур:
ΔΣtп
= Dtп1 + Dtп2
Полезные разности температур по корпусам (в 0C ) равны:
Dtп1 = tг1 – tк1 = 155,18+137,13=18,05,
Dtп2 = tг2 – tк2 = 133,56-95,51=38,05,
Тогда общая полезная разность температур:
ΣΔtп = 18,05+38,05=56,10С.
Проверим общую полезную разность температур:
ΣΔtп = tг1 - tбк – (ΣΔ/ + ΣΔ/// + ΣΔ// ) (2.12)
ΣΔtп
= 155,18-83,53-(6,74+2+6,81)=56,
2.4 Определение тепловых нагрузок
Расход
греющего пара в 1-й корпус, производительность
каждого корпуса по выпаренной воде
и тепловые нагрузки по корпусам определим
путем совместного решения
Q1 = D∙(Ir1 – i1) = 1,03∙ [Gн ∙cн ∙ (tk1 – tH) + w1∙(Iвп1 – cв∙tk1) + Q1конц] (2.13)
Q2 = w1∙(Ir2 – i2) = 1,03∙[(Gн - w1) ∙c1∙ (tк2 – tк1) + w2∙ (Iвп2 – cв∙tк2) + Q2конц] (2.14)
W = w1
+ w2
где 1,03 – коэффициент, учитывающий 3% потерь тепла в окружающую среду;
сн,с1 – теплоемкости растворов соответственно исходного, в первом, кДж/ (кг·К);
сн= 3,76
с1= 2,41
Q1конц, Q2конц– теплоты концентрирования по корпусам, кВт;
tн – температура кипения исходного раствора при давлении в 1–м корпусе;
tн = tвп1 + Δ/н, (2.16)
где Δ/н – температурная депрессия для исходного раствора.
tн = 134,56+ 1 = 135,560С.
При решении уравнений (2.10) – (2.13) можно принять:
Iвп1 » Iг2; Iвп2 » Iбк
Получим систему уравнений:
Q1 = D∙(2724,9-561,6)
= 1,03∙[12,77∙3,76∙(137,13-135,
Q2
= ω1∙(2648,9-349,75) = 1,03∙[(12,77-ω1)∙2,41∙(95,51-
W = w1 + w2 = 11,11.
Решение этой системы уравнений дает следующие результаты:
D = 5,3 кг/с; Q1 =11465,5 кВт; Q2 =11725,6 кВт;
ω1 = 5,1 кг/с; ω2 = 5,4 кг/с;
Таблица 2.4 Результаты расчета
Параметр |
Корпус | |
1 |
2 | |
Производительность по испаряемой воде w, кг/с |
5,1 |
5,4 |
Концентрация растворов х, % |
9,9 |
46 |
Давление греющих паров Рг, Мпа |
0,546 |
0,3003 |
Температура греющих паров tг, °С |
155,18 |
133,56 |
Температурные потери ΣΔ, град |
3,57 |
11,98 |
Температура кипения раствора tк, °С |
137,13 |
95,51 |
Полезная разность температур, Δtп, град |
18,05 |
38,05 |
Наибольшее отклонение вычисленных нагрузок по испаряемой воде в каждом корпусе от предварительно принятых (ω1=4,37кг/с, ω2=5,29кг/с) не превышает 5%, поэтому нет необходимости заново пересчитывть концентрации, температурные депрессии и температуры кипения растворов.
2.5
Выбор конструкционного
Выбираем конструкционный материал, стойкий в среде кипящего раствора NaNO3 в интервале изменения концентраций от 6 до 46 %. В этих условиях химически стойкой является сталь легированного типа марки Х17. Скорость коррозии её не менее 0,1 мм/год, коэффициент теплопроводности λст = 25,1 Вт/(мК).
2.6
Расчет коэффициентов
Коэффициент
теплопередачи для первого
К1 = 1 / (1/α1 + Σδ/λ + 1/α2) (2.17)
Примем,
что суммарное термическое
Σδ/λ = 0,002/25,1 + 0,0005/2 = 2,87∙10-4 м2∙К/Вт.
Коэффициент теплоотдачи от конденсирующегося пара к стенке α1 равен:
α1 = 2,04·4√(r1·ρ2ж 1·λ3ж 1) / (μж 1·Н·Dt1) (2.18)
где r1 – теплота конденсации греющего пара, Дж/кг;
ρж 1,λж 1,μж 1 – соответственно плотность (кг/м3), теплопроводность (Вт/м*К), вязкость (Па*с) конденсата при средней температуре пленки
tпл = tг 1 - Dt1/2.
где Dt1 – разность температур конденсации пара и стенки, °С.
Расчет α1 ведут методом последовательных приближений. В первом приближении примем
Dt1 =4,25 °С.
Тогда
α1 = 2,04∙4√(2100,5∙103∙9132∙0,613)
Для установившегося процесса передачи тепла справедливо уравнение
q = α1·Δt1
= Δtст / (Σδ/λ) = α2·Δt2
где q – удельная тепловая нагрузка, Вт/м2;
Δtст – перепад температур на стенке, °С;
Δt2 – разность между температурой стенки со стороны раствора и температурой кипения раствора, °С.
Распределение температур в процессе теплопередачи от пара через стенку к кипящему раствору показано на рисунке.
Рисунок 1.1 Распределение температур в процессе теплопередачи от
пара к кипящему раствору через многослойную стенку: 1 – пар; 2 – конденсат; 3 – стенка; 4 – накипь; 5 – кипящий раствор.
Δtст = α1∙Δt1∙Σδ/λ = 7976,85∙4,25∙2,87∙10-4 = 9,7 °С.
Тогда
Δt2 = Δtп 1 - Δtст - Δt = 18-9,7-4,25=4 °С.
Коэффициент теплоотдачи от стенки к кипящему раствору для пузырькового кипения в вертикальных трубок при условии естественной циркуляции раствора равен:
α2
=Аq0.6 =780 q0.6
(λ11.3·ρ10.5·ρп10.06/σ10,5·rв1
Подставив численные значения, получим:
α2
=780q0.6 (0,61 1.3∙10620.5∙3,770.06)/0,0720.5
Проверим правильность первого приближения по равенству удельных тепловых нагрузок:
q/ = α1∙Δt1 =7976,85∙4,25=33901,61 Вт/м2;
q// = α2∙Δt2 =8519,13∙4=34248,93 Вт/м2.
Как видим, q/ ≈ q//
Расхождение между тепловыми нагрузками не превышает 3%, расчет коэффициентов α1 и α2 на этом заканчивается.
Находим К1:
К1 = 1/(1/7965,85+ 2,87·10-4 + 1/8519,13) = 1887,7 Вт/ м2∙К.
Далее рассчитаем коэффициент теплопередачи для второго корпуса К2.
К2=1 / (1/α1 + Σδ/λ + 1/α2)
Σδ/λ = 0,002/25,1 + 0,0005/2 = 2,87∙10-4 м2∙К/Вт.
Расчет α1 ведут методом последовательных приближений. В первом приближении примем Δt1=17,4 °С.
α1 = 2,04∙4√(2166∙103∙933,62∙0,623)
Δtст = 2514,65 ∙17,4∙2,87∙10-4 = 12,56 °С;
Δt2 =38-12,56-17,4= 8,04 °С;
α2 =780q0.6 (0,62 1.3∙11040.5∙3,560.06 )/0,0740.5∙ (2148∙103)0.6 ∙0,5790.66∙24100.3∙
·(0,28∙10-3) 0.3 =5424,19
q/ = α1∙Δt1 = 2514,65∙17,4 = 43754,87 Вт/м2;
q// = α2∙Δt2 = 5424,19∙8,04=43623,23 Вт/м2.
Как видим, q/ ≈ q//
К2 = 1/(1/2514,65 + 2,87∙10-4 + 1/5424,19) = 1150,71 Вт/ м2∙К.
2.7 Распределение полезной
Полезные разности температур в корпусах установки находим из условия равенства их поверхностей теплопередачи: