Расчёт и проектирование трёхкорпусной выпарной установки непрерывного действия

Размещено на http://www.allbest.ru/

Содержание

Введение

1. Расчет трехкорпусной выпарной установки

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

1.1.1 Концентрации упариваемого раствора

1.1.2 Температуры кипения растворов

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

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

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

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

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

1.1.8 Уточненный расчет поверхности теплообмена

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

2. Расчет теплообменных аппаратов

2.1 Расчет кожухотрубчатого теплообменника

2.2 Расчет пластинчатого теплообменника

2.3 Выбор оптимального нормализованного теплообменного аппарата

Заключение

Введение

Тепловыми называются технологические процессы, протекающие при условии подвода или отвода тепла.

К тепловым процессам относятся: нагревание, охлаждение, конденсация, теплообмен и испарение. Частным случаем испарения является процесс выпаривания [1].

Нагревание - процесс повышения температуры перерабатываемых материалов путем подвода к ним тепла. Нагревание применяется в химической технологии для ускорения массообменных и химических процессов.

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

Конденсация - процесс сжижения паров вещества путем отвода от них тепла.

Теплообмен - процесс распространения тепла из одной части пространства в другую.

Испарение - перевод в парообразное состояние какой-либо жидкости путем подвода к ней тепла. Частным случаем испарения является весьма широко распространенный в химической технике процесс выпаривания.

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

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

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

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

Выпаривание проводят при атмосферном давлении, под вакуумом или под давлением, большим атмосферного. Образующийся при выпаривании растворов пар называется вторичным, или соковым [3].

Выпаривание под вакуумом имеет ряд преимуществ по сравнению с атмосферной выпаркой: снижается температура кипения раствора, что дает возможность использовать этот способ для выпаривания растворов термически нестойких веществ; повышается полезная разность температур, что ведет к снижению требуемой поверхности теплопередачи выпарного аппарата; несколько снижаются потери теплоты в окружающую среду (так как снижается температура стенки аппарата); появляется возможность использования теплоносителя низкого потенциала. К недостаткам выпаривания под вакуумом относятся удорожание установки (так как требуется дополнительное оборудование - конденсатор, вакуум-насос и др.), а также несколько больший расход греющего пара на 1 кг выпариваемой жидкости (вследствие снижения давления над раствором происходит увеличение теплоты испарения растворителя).

При выпаривании под повышенным давлением (выше атмосферного) вторичный пар может быть использован в качестве греющего агента для различных технологических нужд.

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

При больших производительностях (от нескольких кубических метров выпариваемого раствора в час и выше), что характерно для промышленности, выпаривание проводят по непрерывному принципу. В аппаратах непрерывного действия обычно создают условия для интенсивной циркуляции раствора, т.е. в таких аппаратах гидродинамическая структура потоков близка к модели идеального смешения. Поэтому концентрация раствора в таких аппаратах ближе к конечной, что приводит к ухудшению условий теплопередачи (так, с повышением концентрации раствора увеличивается его вязкость и, следовательно, снижается коэффициент теплоотдачи от стенки к раствору) [6].

Принципиальная технологическая схема трехкорпусной вакуум-выпарной установки непрерывного действия представлена на рисунке 1.1.

Исходный раствор подается из емкости 1 центробежным насосом 2 через теплообменник 3 в первый корпус выпарной установки 4. В теплообменнике 3 исходный раствор нагревается до температуры близкой к температуре кипения раствора в первом корпусе выпарной установки.

Первый корпус установки обогревается свежим (первичным) паром. Вторичный пар, образующийся при кипении раствора в первом корпусе, направляется в качестве греющего пара во второй корпус 5; сюда же поступает частично сконцентрированный раствор из первого корпуса. Аналогично упаренный раствор из второго корпуса подается в третий корпус 6 , обогреваемый вторичным паром второго корпуса. Упаренный до конечной концентрации в третьем корпусе готовый продукт поступает из него в емкость 10. По мере прохождения из корпуса в корпус давление и температура пара понижаются, и из последнего (третьего) корпуса пар с низким давлением отводится в барометрический конденсатор смешения 7, в котором при конденсации пара создается вакуум. Раствор и вторичный пар перемещаются из корпуса в корпус самотеком благодаря общему перепаду давления, возникающего в результате избыточного давления в первом корпусе и вакуума в последнем. Воздух и неконденсирующиеся газы, поступающие в установку с охлаждающей водой (в конденсаторе) и через не плотности трубопроводов, отсасываются через ловушку 8 вакуум-насосом.

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

Задание на курсовое проектирование

Рассчитать и спроектировать трехкорпусную выпарную установку непрерывного действия для концентрирования G_н=20000 кг/ч водного раствора Na₂CO₃ от начальной концентрации х_н=3% до конечной х_к=35 % при следующих условий:

1. обогрев производится насыщенным водяным паром давлением Р_г1=0,98 МПа;

2. давление в барометрическом конденсаторе Р_бк=0,0036 МПа;

3. выпарной аппарат - тип-2, исполнение-2;

4. взаимное направление пара и раствора - прямоток;

5. отбор экстра-пара не производится;

6. раствор поступает в первый корпус подогретым на 3-5 градусов ниже температуры кипения раствора.

теплопередача выпарной аппарат

1. Расчет трехкорпусной выпарной установки

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

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

(1)

Для определения тепловых нагрузок Q, коэффициентов теплопередачи К и полезных разностей температур Дt_п необходимо знать распределение упариваемой воды, концентраций растворов и их температур кипения по корпусам. Эти величины находят методом последовательных приближений.

Первое приближение

Определяем производительность установки в кг/с

G_н=20000/3600=5,55 кг/с (2)

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

(3)

кг/с

1.1.1 Концентрации упариваемого раствора

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

(4)

кг/с (5)

кг/с (6)

кг/с (7)

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

% (8)

% (9)

% (10)

Концентрация раствора в последнем корпусе соответствует заданной концентрации упаренного раствора .

1.1.2 Температуры кипения растворов

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

^МПа (11)

В первом приближении общий перепад давлений распределяются между корпусами поровну. Тогда давления греющих паров в корпусах равны:

Р_г1=0,98 МПа;

 МПа (12)

МПа (13)

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

МПа (14)

что соответствует заданному значению Р_6к.

По давлениям паров находим их температуры и энтальпии (табл.П57 [9]) полученные данные введем в таблицу 1.

Таблица 1 - Параметры в корпусах и барометрической камере

Давление, МПа	Температура, оС	Энтальпия, кДж/кг
Рг1=0,98	tг1=177,2	Iг1=2782,4
Рг2=0,65453	t г2=160,5	Iг2 =2768,4
Рг3=0,32906	tг3=134,9	Iг3=2732,8
Рбк=0,0036	tбк=26,1	Iбк=2543,5

При определении температуры кипения растворов в аппаратах исходят из следующих допущений. Распределение концентраций раствора в выпарном аппарате с интенсивной циркуляцией практически соответствует модели идеального перемешивания. Поэтому концентрацию кипящего раствора принимают равной конечной в данном корпусе и, следовательно, температуру кипения раствора определяют при конечной концентрации.

Изменение температуры кипения по высоте кипятильных труб происходит вследствие изменения гидростатического давления столба жидкости. Температуру кипения раствора в корпусе принимают соответствующей температуре кипения в среднем слое жидкости. Таким образом, температура кипения раствора в корпусе отличается от температуры греющего пара в последующем корпусе на сумму температурных потерь ?Д от температурной (), гидростатической (Д") и гидродинамической (Д"') депрессий (?Д=+Д"+Д"').

Гидродинамическая депрессия обусловлена потерей давления пара на преодоление гидравлических сопротивлений трубопроводов при переходе из корпуса в корпус. Так в расчетах принимают Д"'=1,0-1,5 ⁰С на корпус. Примем для каждого Д"'=1 °С. Тогда температуры вторичных паров в корпусах будут равны:

^°С (15)

^°С (16)

^°С (17)

Сумма гидродинамических депрессий:

=1+1+1=3 °С. (18)

По температурам вторичных паров определим их давления и теплоту парообразования (табл. 56, [9]).

Таблица 2 - Параметры вторичных паров в корпусах

Температура, оС	Давление, МПа	Теплота парообразования, кДж/кг
tвп1=161,5	Рвп1=0,670	rвп1=2084,0
tвп2=135,9	Рвп2=0,335	rвп2 =2159,8
tвп3=27,1	Рвп3=0,0041	rвп3=2432,1

Гидростатическая депрессия обусловлена разностью давлений в среднем слое кипящего раствора и на его поверхности. Давление в среднем слое кипящего раствора Р_ср каждого корпуса определяется по уравнению:

(19)

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

с- плотность кипящего раствора, кг/м³;

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

Для выбора значения Н необходимо ориентировочно оценить поверхность теплопередачи выпарного аппарата F_оp. При кипении водных растворов можно принять удельную тепловую нагрузку аппаратов с естественной циркуляцией q = =20 000 - 50 000 Вт/м², аппаратов с принудительной циркуляцией q = 40 000 - 80 000 Вт/м². Примем q= 50000 Вт/м². Тогда поверхность теплопередачи 1-го корпуса ориентировочно равна:

м² (20)

м² (21)

м² (22)

где r₁ -- теплота парообразования вторичного пара, Дж/кг.

По ГОСТ 11987--81 [4], трубчатые аппараты с принудительной циркуляцией и вынесенной греющей камерой (тип 2, исполнение 2) состоят из кипятильных труб высотой 6 м при диаметре d_н = 38 мм и толщине стенки д_ст= 2 мм.

При пузырьковом (ядерном) режиме кипения паронаполнение составляет е=0,4--0,6. Примем е = 0,6.

Плотность раствора Na₂CO₃, находиться по табл. 181 [10].

с₁ = 1041,0 кг/м³; с₂= 1068,0 кг/м³; с₃=1388,0 кг/м³. (23)

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

(24)

(25)

(26)

Этим давлениям соответствуют следующие температуры кипения и теплоты испарения растворителя (табл. П57 ; [9])

Таблица 3 - Давление, температура и удельная теплота испарения

Р, МПа	t,оС	r, кДж/кг
Р1ср=0,682	t1ср=163,1	r1 =2078,6
Р2ср=0,367	t2ср=139,6	r2=2151,0
Р3ср=0,0204	t3ср=59,7	r3=2358,0

Находим гидростатическую депрессию по корпусам:

=t_1ср- t_вп1=163,1-161,5=1,6 ^оС (27)

=t_2ср- t_вп2=139,6-135,9=3,7 ^оС (28)

=t_3ср- t_вп3=59,7-27,1=32,6 ^оС (29)

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

?=++=1,6+3,7+32,6=37,9 ^оС (30)

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

= 1,62. (31)

где Т - температура паров в среднем слое кипятильных труб, К;

_атм - температурная депрессия при атмосферном давлении [5].

Находим значение по корпусам:

^єС (32)

^°С (33)

^єС (34)

Сумма температурных депрессий:

^°С (35)

В аппаратах с вынесенной зоной кипения как с принудительной, так и с естественной циркуляцией кипение раствора происходит в трубе вскипания, устанавливаемой над греющей камерой. Кипение в греющих трубках предотвращается за счет гидростатического давления столба жидкости в трубе вскипания. В греющих трубках происходит перегрев жидкости по сравнению с температурой кипения на верхнем уровне раздела фаз. Поэтому температуру кипения раствора в этих аппаратах определяют без учета гидростатических температурных потерь. Перегрев раствора t_пер может быть найден из внутреннего баланса тепла в каждом корпусе. Уравнение теплового баланса для j-го корпуса записывается в следующем виде:

G_нj•c_нj• (t_кj-1 - t_кj) + M•c_нj•t_перj = щ_j• (_{Iвп j} - c_в•tк_j) (36)

Для первого корпуса t_кj-1 - это температура раствора, поступающего в аппарат из теплообменника-подогревателя.

В аппаратах с естественной циркуляцией обычно достигаются скорости раствора u = 0,6 - 0,8 м/с. Примем u = 0,8 м/с. Для этих аппаратов масса циркулирующего раствора равна:

M = v • S • с, (37)

где S- сечение потока в аппарате (м²), рассчитываемая по формуле:

(38)

где d_вн - внутренний диаметр труб, d_вн =0,034 м;

Н - принятая высота труб, H= 6 м.

м² (39)

M₁ = 0,8*0,083*1041 = 69,1 кг/с (40)

м² (41)

M₂ = 0,8*0,103*1068 = 88,0 кг/с (42)

м² (43)

M₃ = 0,8*0,136*1388 = 151,0 кг/с (44)

Таким образом, температура перегрева в j-м аппарате равна

?t_{пер j} = (45)

Для расчета удельной теплоемкости двухкомпонентных (вода + растворенное вещество) разбавленных водных растворов (х < 0,2) пользуются приближенной формулой:

(46)

Здесь 4190 Дж/(кг·К) - удельная теплоемкость воды; Х - концентрация растворенного вещества, масс. доли.

С_н1=4190(1-0,0402)=4020 Дж/(кг·К)=4,02 кДж/(кг·К) (47)

С_н2=4190(1-0,068)=3905 Дж/(кг·К)=3,905 кДж/(кг·К) (48)

Для концентрированных двухкомпонентных водных растворов (Х > 0,2) расчет ведут по формуле [9]

, (49)

где c₁ - удельная теплоемкость безводного растворенного вещества, Дж/(кг·К).

Удельную теплоемкость химического соединения при отсутствии экспериментальных данных можно ориентировочно рассчитать по уравнению

(50)

где М - молекулярная масса химического соединения;

с - его массовая удельная теплоемкость, Дж/(кг · К);

n₁, n₂, n₃, ... - число атомов элементов, входящих в соединение;

c₁, с₂, c₃, ... - атомные теплоемкости, Дж/(кг•атом К).

Мс_н3=кДж/(кг·атом·К) (51)

С_н3=4190(1-0,35)+154^. 0,35=2777,4 Дж/(кг·К)=2,77 кДж/(кг·К) (52)

С_в - теплоемкость воды, С_в =4,19 кДж/(кг·К)

?t_пер1 = (53)

?t_пер1 =

?t_пер2 = (54)

?t_пер2 =

?t_пер3 = (55)

?t_пер3 =

Температурные потери (56)

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

t_к1 = t_r2 + ?₁? + ?₁???+ = 160,5 + 0,65 + 1,0 += 165,4 °С (57)

t_к2 = t_r3 + ?₂? + ?₂???+ = 134,9+0,96+1,0 +=140,0 °С (58)

t_к3 = t _{бк 3} + ?₃? + ?₃???+=26,1+4,03+1,0+=34,21°С (59)

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

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

^°С; (60)

^°С; (61)

^°С. (62)

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

?Дt_п =11,8+20,5+100,69= 132,99 ^оС. (63)

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

(64)

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

Расход греющего пара в первом корпусе, производительность каждого корпуса по выпариваемой воде и тепловые нагрузки по корпусам определяются путем совместного решения уравнений тепловых балансов по корпусам и уравнениями баланса по воде для всей установки:

Q₁ = D •(I_r1 - i₁) = 1,03 • [G_н • с_н •(t_k1 - t_н) + w₁ • (I_вп1 - c_в • t_k1) + Q_1конц]; (65)

Q₂ = w₁ •(I_r2-i₂)=1,03•[(G_н-w₁)•c₁ •(t _k2-t_k1)+w₂ • (I_вп2 - c_в • t_k2) + Q_2конц]; (66)

Q₂=w₂ • (I_r2 - i₂)=l,03 •[(G_н-w₁-w₂)•c₂ •(t_k3-t_k2)+w₃ •(I_вп3 -c_в•t_k3)+Q_3конц]; (67)

W = w₁ + w₂ + w₃. (68)

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

Q₁, Q₂, Q₃ - тепловые нагрузки по корпусам, кВт;

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

I₁, I₂, I₃ - энтальпии греющих паров по корпусам, кДж/кг;

I _вп1, I _вп2, I _вп3 - энтальпии вторичных паров по корпусам, кДж/кг; t_н - температура кипения исходного раствора при давлении в 1-м корпусе;

i₁, i₂, i₃ - энтальпии конденсата по корпусам, кДж/кг [9];

с_в - теплоёмкость воды, кДж/(кг•К);

с₁, с₂, с₃ - теплоёмкости раствора с начальной концентрацией в первом корпусе и втором корпусах соответственно, кДж/(кг•К) [10]; Q_1конц, Q_2конц, Q_3конц - теплота концентрирования по корпусам, кВт;

t_н - температура кипения исходного раствора при давлении в 1-м корпусе, °С, ?_н? - температурная депрессия для исходного раствора, равная

t_нl = t_вп1 + ?_н? = 161,5 + 0,33= 161,83 °С (69)

При решении уравнений (65) - (68) можно принять:

; ;

Анализ зависимостей теплоты концентрирования от концентрации и температуры показал, что она наибольшая для третьего корпуса. Поэтому рассчитаем теплоту концентрирования для 3-го корпуса:

Q_3конц = G_cyx • ?q = G_п • Х_н • ?q , (70)

где G_cyx - производительность аппарата по сухому веществу Na₂CO₃, кг/с; ?q - разность интегральных теплот растворения при концентрациях Х₂ и Х₃, кДж/кг [10]

Q_3конц = 5,55 *0,03 • (338,1-288,7) = 8,22 кВт.

Сравним Q_3конц с ориентировочной нагрузкой для 3-го корпуса Q_3ор:

Q_3ор = (G_н - w₁ - w₂) c₂ (t_k3 - t_k2) + w₃(I_в3 - c_в•t_k3) (71)

Q_3ор = (5,55 - 1,41-1,69) • 3,905 • (63,73-140,56) + 1,97 • (2543,5 - 4,19•63,73) = 3749,6 кВт.

(72)

Q_3конц составляет значительно меньше 3 % от Q_3ор, поэтому в уравнениях тепловых балансов по корпусам пренебрегаем величиной Q_3конц.

I₁=2782,4 кДж/кг ; i₁=667,9 кДж/кг

I₂=2768,4 кДж/кг ; i₂=572,7 кДж/кг

I₃=2732,8 кДж/кг ; i₃=106,5 кДж/кг

Получим систему уравнений:

Q₁ = D (2782,4-678,5) = 1,03[5,55•4,06•(165,4- 161.83) •W₁•(2768,4 -4,19•165,4)]

Q₂ = W₁ (2768,4-567,3) = 1,03 [(5,55 - W₁) • 4,02• (140,0-165,4)• •W₂ (2732,8 - 4,19•140,0)]

Q₃ = W₂ (2732,8-106,5) = 1,03 [(5,55 - W₁ - W₂) • 3,905 • (34,21-140,0) • • W₃ (2543,5 - 4,19 •34,21)];

5,07 = W₁ + W₂ + W₃

Решение системы уравнений даёт следующие результаты: D = 1,498 кг/с;

W₁ = 1,436 кг/с; Q₁ = 3151,64 кВт;

W₂ = 1,651 кг/с; Q₂ = 3160,78 кВт;

W₃ = 1,9885 кг/с; Q₃ = 4336,02 кВт.

Результаты сведем в таблицу 4

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

Наименование параметра	1-й корпус	2-й корпус	3-й корпус
Производительность по упариваемой воде W, кг/с	1,436	1,651	1,9885
Концентрация растворов Х, %	4,02	6,8	35,0
Давление греющих паров Рr, Па	0,980	0,65453	0,32906
Температура греющих паров tr, °C	177,2	160,5	134,9
Температурные потери ??, град.	4,9	5,1	8,11
Температура кипения раствора tк, °С	165,4	140,0	34,21
Полезная разность температур ?tn, °С	11,8	20,5	100,69

Наибольшее отклонение вычисленных нагрузок по испаряемой воде в каждом корпусе от предварительно принятых (W₁= 1,41 кг/с, W₂ = 1,69 кг/с, W₃ = 1,97 кг/с) не превышает 3 %, поэтому в дальнейших расчётах не производим пересчёт концентраций и температур кипения растворов по корпусам. В случае, если это расхождение составит больше 5 %, необходимо заново пересчитать концентрации, температурные депрессии и температуры кипения растворов, положив в основу расчёта новое, полученное из решения балансовых уравнений распределение по корпусам нагрузок по испаряемой воде.

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

Выбираем конструкционный материал, стойкий в среде кипящего раствора Na₂CO₃ в интервале изменения концентраций от 3 до 35 % [2]. В этих условиях химически стойкой является сталь марки Х17. Скорость коррозии ее не менее 0,1 мм/год, коэффициент теплопроводности л_ст= 25,1 .

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

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

^. (73)

Примем, что суммарное термическое сопротивление стенки и на- кипи равно , при этом не учитываем термическое сопротивление загрязнении со стороны пара.

, м² К/Вт. (74)

Таблица 5 - Физические свойства кипящих растворов Na₂CO₃

Параметр	Корпус	Литература
	1	2	3
Теплопроводность раствора л, Вт/(м К)	0,681	0,665	0,625	[10]
Плотность раствора с, кг/м3	1041	1068	1388	[10]
Теплоемкость раствора с, Дж/ (кг К)	4020	3905	2770	[9]
Вязкость раствора м, Па с	0,169 10-3	0,196 10-3	0,734 10-3	[9]
Поверхностное натяжение у, Н/м	0,0453	0,05072	0,0703	[9]
Теплота парообразования rв, Дж/кг	2071,2 10-3	2150 10-3	2416,08 10-3	[9]
Плотность пара сн, кг/м3	3,716	1,962	0,038	[9]

Коэффициент теплоотдачи от конденсирующегося пара к стенке б₁ находится по уравнению:

(75)

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

с_ж1, л_ж1, м_ж1 --соответственно плотность (кг/м³), теплопроводность (Вт/(м·К)), динамическая вязкость (Па с) конденсата при средней температуре пленки t_пл=t_г1-Дt₁/2; Дt₁-- разность температур конденсации пара и стенки, град.

t_пл=t_г1-Дt₁/2=177,2-2/2=177,2°С; (76)

Расчет б₁, ведут методом последовательных приближений. В первом приближении примем разность температур конденсации пара и стенки Д₁ =2,0 ^оС.

Тогда при t_пл= 177,2°С:

r₁ =2030,8·10³ Дж/кг;

с_ж1=889,8 кг/м³;

л_ж1=0,67612 Вт/(м·К);

м_ж1 =0,158·10^-3 Па с

(77)

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

Дt_ст= б₁ Д?д/л=8207,23•2•2,87•10^-4=4,71 °С. (78)

И разность температур конденсации пара и стенки со стороны раствора и его температурой кипения Д:

Д₂= Дt₁- Дt_ст - Д₁=11,8-4,71-2=5,09 ^оС. (78)

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

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

Коэффициент теплоотдачи от стенки к кипящему раствору б₂, для пузырькового кипения в вертикальных кипятильных трубках при условии естественной циркуляции раствора находиться по уравнению:

Вт/(м²·К) (79)

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

= б₁ Д₁=6=8207,23= 16414,46 Вт/м²; (80)

= б₂ Д₂=6391,4= 32532,2 Вт/м². (81)

Как видим, .

Для второго приближения примем Д₁ = 3,0 ^оС.

Рассчитаем б₁ по соотношению

Вт/(м²·К).

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

^°С;

Дt₂=11,8-6,38-3,0=2,42 єС;

^{Вт/(мІ·К);}

=7416,07= 22248,21 Вт/м²;

=7670,7= 18409,7 Вт/м².

Очевидно, что и во втором приближении остается неравенство удельных тепловых нагрузок q^, ? q^,,.Это приводит к следующему приближенному расчету.

В третьем приближении по графической зависимости удельной тепловой нагрузки q от разности температур между паром и стенкой в первом корпусе и определяем Д₁=2,78 ^оС. Получим:

^;

Дt₂=11,8-6,03-2,78=2,99;

^;

=7957•2,78=22120,0 Вт/м²;

=7558•2,99=22598 Вт/м².

Как видим, .

Если расхождение между тепловыми нагрузками не превышает 3 %, расчет коэффициентов б₁ и б₂ на этом заканчивают. Находим К₁, Вт/(м²·К).по уравнению (73):

.

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

t_пл=t_г2-Дt₁/2=160,5-5,9/2=157,55°С; (82)

Тогда при t_пл= 157,55°С:

r₂ =2095,2·10³ Дж/кг;

с_ж2=909,0 кг/м³;

л_ж2=0,6832 Вт/(м·К);

м_ж2 =0,176·10^-3 Па с

^;

^°С;

Дt₂=20,5-10,6-5,9=4,0^оС;

^;

= 6258,9= 36927,51 Вт/м²;

= 9406,02•4,0= 37624,1 Вт/м².

Как видим,

Определим К₂:

⁽⁸³⁾

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

t_пл=t_г3-Дt₁/2=134,9-42/2=113,9°С; (84)

Тогда при t_пл= 113,9°С:

r₃ =2223,6·10³ Дж/кг;

с_ж3=948,0 кг/м³;

л_ж3=0,6854 Вт/(м·К);

м_ж3 =0,246·10^-3 Па с

;

^°С;

Дt₂ = 100,69-42,0-44,14=14,55 ^оС;

;

=3661,6•42= 153787,8 Вт/м²;

= 10849,2•14,55= 157855,68 Вт/м².

Как видим,

Найдем К₃:

. (85)

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

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

 (86)

где , Q_j, K_j - полезная разность температур, тепловая нагрузка, коэффициент теплопередачи для j-го корпуса.

(87)

= 36,94 °С; (88)

= 59,81°С. (89)

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

?= ++= 36,24+36,64+59,81= 132,99 °С. (90)

Рассчитаем поверхность теплопередачи выпарных аппаратов по формуле:

; (91)

м²

м²

м²

Найденные значения поверхности теплопередачи выпарных аппаратов отличаются незначительно от ориентировочно определённой ранее F_op. Поэтому в последующих приближениях нет необходимости вносить коррективы на изменение конструктивных размеров аппаратов (высота, диаметр и число труб). Сравнение распределённых из условий равенства поверхностей теплопередачи и предварительно рассчитанных значений полезных разностей температур представлено в таблице 6.

Таблица 6 - Сравнение распределённых и рассчитанных значений полезных разностей температур

Корпуса	1	2	3
Распределённые в 1-м приближении	36,23	36,94	59,81
Рассчитанные	11,8	20,5	100,69

Второе приближение

Как видно из таблицы 6, рассчитанные полезные разности температур (из условия равного перепада давления в корпусах) и распределённые в 1-м приближении (из условия равенства поверхности теплопередачи в корпусах) существенно различаются. Поэтому необходимо заново перераспределять температуры (давления) между корпусами установки. В основу этого перераспределения температур (давлений) должны быть положены полезные разности температур, найденные по результатам распределения общей полезной разности из условия равенства поверхностей теплопередачи.

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

В связи с тем, что существенное изменение давлений, по сравнению с первым приближением, происходит только в 1 и 2 корпусах, где суммарные температурные потери незначительны, во втором приближении принимаем значения ?', ?'' и ?''' для каждого корпуса такие же, как и в первом приближении. Полученные после перераспределения температур (давлений) параметры растворов и паров по корпусам представлены в таблице 7.

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

Параметры	1 корпус	2 корпус	3 корпус
Производительность по выпаренной воде щ, кг/с	1,436	1,651	1,9885
Концентрация растворов X, %	4,02	6,8	35,0
Температура греющего пара в 1-м корпусе tг, °С	177,2	-	-
Полезная разность температур ?tn, °C	36,23	36,94	59,81
Температура кипения раствора tк = tг - ?tп, °C	140,96	99,15	35,24
Температура вторичного пара tв = tк - (?'+?''), °C	137,07	95,05	28,12
Давление вторичного пара Рв, Па	0,346	0,088	0,004
Температура греющего пара tг = tв - ?''', °C	-	136,07	94,05

Рассчитаем тепловые нагрузки (в кВт):

Q₁ = 1,03 •[5,55•4,06•(140,96-137,07) + 1,436•(2741,3 - 4,19•140,96)] =3270,4 кВт;

Q₂ = 1,03•[4,114•4,02•(99,15-140,96) + 1,651•(2677,6 - 4,19•99,15)] =3155,0 кВт;

Q₃ = 1,03•[2,436•3,905•(35,24-99,15)+1,9885•(2561,4 - 4,19•35,24)]=4493,4 кВт.

Расчет коэффициентов теплопередач выполняются тем же методом описанным выше. Из физических характеристик меняются те, которые зависят от Т_вп, они приведены в таблице 8.

Таблица 8 - Физические характеристики

Параметр	Корпус
	1	2	3
Вязкость раствора м, Па с	0,200 10-3	0,298 10-3	0,840 10-3
Поверхностное натяжение у, Н/м	0,0513	0,0597	0,0714
Теплота парообразования rв, Дж/кг	2158,8 10-3	2272,9 10-3	2429,8 10-3
Плотность пара сн, кг/м3	1,82	0,504	0,027

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

t_пл=t_г1-Дt₁/2=177,2-12,5/2=170,95°С;

Тогда при t_пл= 170,95°С:

r₁ =2052,5·10³ Дж/кг;

с_ж1=896,0 кг/м³;

л_ж1=0,679 Вт/(м·К);

м_ж1 =0,162·10^-3 Па с

^;

^°С;

Дt₂=36,23-12,5-18,66=5,08^оС;

^;

= 5201,55= 65019,4 Вт/м²;

= 13083,9•5,08= 66466,5 Вт/м².

Как видим, . Определим К₁:

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

t_пл=t_г2-Дt₁/2=136,07-12,3/2=129,92°С;

Тогда при t_пл= 129,92°С:

r₂ =2179,0·10³ Дж/кг;

с_ж2=935,0 кг/м³;

л_ж2=0,686 Вт/(м·К);

м_ж2 =0,212·10^-3 Па с

^;

^°С;

Дt₂=36,94-18,03-12,3=6,61^оС;

^;

= 5107,9= 62827,2 Вт/м²;

= 9374,25•6,61= 61967,1 Вт/м².

Как видим, . Определим К₂:

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

t_пл=t_г3-Дt₁/2=94,05-22,2/2=82,95°С

Тогда при t_пл= 82,95°С:

r₃ =2302,2·10³ Дж/кг;

с_ж3=969,9 кг/м³;

л_ж3=0,676 Вт/(м·К);

м_ж3 =0,343·10^-3 Па с

;

^°С;

Дt₂ = 59,81-25,42-22,2=12,19 ^оС;

;

=3990,6•22,2= 88591,6 Вт/м²;

= 7253,2•12,19= 88416,4 Вт/м².

Как видим, . Найдем К₃:

.

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

= 36,88 °С;

= 60,12°С.

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

? 35,99+36,88+60,12=132,99 °С.

Сравнение полезных разностей температур , полученных во втором приближении, приведены в таблице 9.

Таблица 9 - Сравнение полезных разностей температур

Корпуса	1	2	3
Распределённые во 2-м приближении	35,99	36,88	60,12
Распределённые в 1-м приближении	36,24	36,94	59,81

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

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

м²;

м²;

м²;

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

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

Действительная поверхность теплообмена F_д 50,43 м²

Диаметр труб d 38Ч2 мм

Высота труб Н 6000 мм

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

Диаметр сепараторов d_c 1600 мм

Общая высота аппарата Н_а 19500 мм

Масса аппарата М_а 9500 кг

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

Толщину тепловой изоляции находим из равенства удельных тепловых потоков через слой изоляции в окружающую среду:

, (92)

где - коэффициент теплоотдачи от внешней поверхности изоляции к воздуху, Вт/(м² К) ;

- температура изоляции со стороны воздуха, °С; Для аппаратов, работающих внутри помещения выбирают в пределах 35 ч 45 єС, а для аппаратов, работающих на открытом воздухе в зимнее время - в интервале 0 ч 10 єС.;

- температура изоляции со стороны аппарата, єС (температуру t_ст1 можно принимать равной температуре греющего пара, ввиду незначительного термического сопротивления стенки аппарата по сравнению с термическим сопротивлением слоя изоляции);

- температура окружающей среды (воздуха), єС;

- коэффициент теплопроводности изоляционного материала, Вт/(мК).

В качестве изоляционного материала выбираем совелит, который содержит 85% магнезии и 15 % асбеста. Коэффициент теплопроводности совелита

Вт/(м*К)

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

(93)

Такую же толщину тепловой изоляции принимаем для второго и третьего корпусов.

2. Расчет теплообменных аппаратов

Тепловой расчет греющего пара

Температурная схема

190>190

160<30

?t_б=30 ?t_м=160

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

(94)

Средняя температура KOH = 95

(95)

р=995 м³/с - плотность KOH при концентрации =2

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

(Вт) (97)

Расход пара:

(98)

Горячая жидкость при средней температуре имеет следующие физико-химические характеристики: м₁=0,03016 Па с, кг/м³, С₁ =4106,2 Дж/кгК.

Приняли ориентировочный теплообменник по ГОСТ 15118-79

Диаметр, D, мм 400

Поверхность теплообмена, F,м² 80

Число ходов 1

Общее число труб, шт 111

Коэффициент теплоотдачи для сырья.

Уточняем Re:

Режим движения сырья - турбулентный.

Вычислим критерий Прандтля для сырья

, (72)

.

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

б=(0,517/0,021)*0,021*23340^0,8*(2,17)^0,43=2252,4 Вт/(м²*К)

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

Принимаем тепловую проводимость греющего пара ?5800Вт/(м²*К)

Коэффициент теплопроводности стали Вт/м*К табл. XXXVI [Павлов]

Тогда

1/Уr_ст=1/(1/5800+0,002/46,5+1/5800)=2580Вт/(м²*К)

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

Требуемая площадь поверхности теплообмена:

м²

Запас площади поверхности теплообмена:

Запас площади поверхности теплообмена достаточен.

Приняли кожухотрубчатый теплообменник по ГОСТ 15118-79

Диаметр, D, мм 600

Поверхность теплообмена, F,м² 80

Число ходов 1

Длина труб, м 4,0

Общее число труб, шт 154

Диаметр труб, мм 25*2

2.1 Расчет пластинчатого теплообменника

Предназначен для нагрева раствора KOH

Приняли по ГОСТ 15518-78

Тип пластины 0,5

Размер пластины, м:

Длина 1370

Ширина 500

Толщина 1

Эквивалентный диаметр канала, м 0,0080

Площадь сечения канала, м² 0,0018

Приведенная длина канала, м 1,15

Количество пластин, шт 154

Число последовательно соединенных пакетов определили

=

=0,01*((20000*(0,0018)²*(154)²)^0,33/(0,01)²=2,35

Число ходов z=3

Число каналов в одном пакете

1)для пара

m₁=n/2z=154/2*3=25,67

2)для раствора

m₂=((n/2)-1)/z = (154/2-1)/3=25,33

Площадь сечения пакетов

f₁=m₁f=25,67·0,0018=0,046 м²

f₂=m₂f=25,33·0,0018=0,046 м²

w₂=11,67/(995·0,046)=0,25 м/с

Критерий Re для раствора

Re=(0,25*0,008*995)/(0,3016*10^-3)= 6598

Критерий Прандля для раствора

Pr=(0,3016*10^-3*4106,2)/0,58 = 2,12

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

б=( 0,58/0,008)*0,0978*(6598)^0,73*(2,12)^0,43=6013 вт/(м²*К)

Критерий Re для пара определили

Re=(3,14*1,15)/(80*0,3016*10^-3 )= 150

Критерий Pr для пара

Pr=((0,3016*10^-3)*4480)/0,67 = 2

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

б= (0,67/1,15)*322*2^0,4*150^0,7=8224,7

Уr =0,635*10^-3

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

K₂ = 1/(1/6013+0,635*10^-3+1/8224,7)=1111

Проверим правильность принятого допущения относительно ?t

?t=?t_ср*К/б₂=78,3*1111/8224,7=10,5 ⁰С

Поверхность теплообмена

F=6229516/78,3*1111 = 71,6 м²

Определяем площадь поверхности запаса:

Приняли по ГОСТ пластинчатый теплообменник со следующими характеристиками

Площадь теплообмена, м² 80

Число пластин, шт 154

Площадь одной пластины 0,5

Эквивалентный диаметр канала, м 0,0080

Поперечное сечение клапана, м² 0,0018

Диаметр условного прохода, мм 150

Приведенная длина канала, м 1,15

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

Определяем скорость жидкости в штуцерах:

Определяем коэффициент трения:

Для однопакетной компоновки пластин Z=5 гидравлическое сопротивление

Сопоставим заданный напор с расчётным гидравлическим сопротивлением: 20000/15175>1. Как видим, действительное гидравлическое сопротивление находится в пределах принятого значения.

Заключение

Выполнили расчеты нахождения коэффициентов теплопередачи процесса выпаривания. Они равны: К₁=1120, К₂=1111,

В результате выбрали по ГОСТ 11987--81 [1, с. 183] выпарной аппарат со следующими характеристиками:

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

Диаметр труб d 38*2 мм

Высота труб H 6000мм

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

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

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

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

Масса аппарата M_а 3000кг

Затем провели расчеты теплообменных аппаратов. Выбрали кожухотрубчатый теплообменник по ГОСТ 15118-79

Диаметр, D, мм 600

Поверхность теплообмена, F,м² 80

Число ходов 1

Длина труб, м 1,37

Общее число труб, шт 154

Размещено на Allbest.ru