Ректификационная установка для разделения бинарной смеси "бензол-толуол"
Введение
Ректификация – разделение жидких однородных смесей на составляющие вещества или группы составляющих веществ в результате противоточного взаимодействия паровой смеси и жидкой смеси.
Возможность разделения жидкой смеси на составляющие её компоненты ректификацией обусловлена тем, что состав пара, образующегося над жидкой смесью, отличается от состава жидкой смеси в условиях равновесного состояния пара и жидкости.
Сущность процесса ректификации рассмотрим на простейшем примере разделения двухкомпонентной смеси, как и в случае нашего задания по курсовому проектированию, где требуется спроектировать ректификационную установку для разделения смеси «бензол-толуол». При ректификации исходная смесь делится на две части: часть, обогащенную легколетучим компонентом (ЛЛК), называемую дистиллятом, и часть, обедненную ЛЛК, называемую кубовым остатком.
Обычно ректификационный аппарат состоит из двух частей: верхней и нижней, каждая из котрых представляет собой организованную поверхность контакта фаз между паром и жидкостью.
В нижней части исходная смесь взаимодействует с паром, начальный состав которого равен составу кубового остатка. Вследствие этого из смеси извлекается легколетучий компонент.
В верхней
ступени пар начального состава
соответствующий составу
Пар для питания ректификационной колонны получается многократным испарением жидкости, имеющей тот же состав, что и кубовый остаток, а жидкость – многократной конденсацией пара, имеющего состав, одинаковый с составом дистиллята.
Задание
Рассчитать ректификационную установку для разделения бинарной смеси бензол – толуол.
Исходные данные для курсового проекта:
1. Производительность по сырью F = 10 т/ч
2. Концентрация низкокипящего
а) в сырье
б) в дистилляте
в) в кубовом остатке
3. Тип тарелок – ситчатые.
4. Давление в колонне
5. Смесь подается в колонну
при температуре начала
Содержание расчетно-
1. Расчет материального баланса.
2. Построение равновесной и
3. Определение скорости пара и диаметра колонны.
4. Расчет гидравлического
5. Определение числа тарелок и высоты колонны.
6. Расчет теплового баланса
7. Расчет и подбор
8. Расчет и подбор насоса.
9. Расчет и подбор штуцеров.
Заключение
Графическая часть
1. Общий вид колонны с узлами- 1 лист формата А1.
2. Общий вид кипятильника с узлами – 1 лист А1.
1. Расчет материального баланса
Температуры кипения и молекулярные массы разделяемых компонентов
tк, °С |
Мол. масса, кг/кмоль | |
бензол |
80,2 |
78 |
толуол |
110,8 |
92 |
Таким образом, бензол является легкокипящим компонентом, а толуол – тяжелокипящим компонентом.
Обозначим массовый расход дистиллята через GD кг/ч и расход кубового остатка через GW кг/ч.
Из уравнения материального баланса
по потокам:
GD + GW = GF
по бензолу:
0, 96· GD + 0,02 · GW = 0,50· GF
Откуда:
GD = 10000 – 4894 = 5106 кг/ч = 1,418 кг/с
GF = 10000 кг/ч = 2,777 кг/с
Все расчеты в данном случае ведутся для легкокипящего компонента, а значит х – концентрация бензола. Для дальнейших расчетов необходимо пересчитать составы фаз из массовых в мольные по соотношению:
где x - мольная доля ЛКК – бензола,
- массовая доля бензола, % (масс.)
Мб – мольная масса бензола, кг/кмоль
Мт – мольная масса толуола, кг/кмоль
Пересчитываем массовый состав смеси в мольный:
Молекулярные массы исходной смеси, дистиллята и кубового остатка определяем по формуле:
= 0,541 × 78 + (1 - 0,541) × 92 = 84,43 кг/моль
0,966 × 78 + (1 - 0,966) × 92 = 78,48 кг/моль
= 0,025 × 78 + (1 - 0,025) × 92 = 91,65 кг/моль
2. Построение равновесной и рабочей линии (у – х) и диаграммы t – х, у
Данные по равновесию приведены в таблице 1 (по данным табл. 2, стр. 153, [1]):
Таблица 1. Равновесные составы жидкости х и пара у, мол. %, температуры кипения t для смеси бензол – толуол при давлении 0,1013 МПа.
Смесь – бензол - толуол.
x% |
0 |
5 |
10 |
20 |
30 |
40 |
50 |
60 |
70 |
80 |
90 |
100 |
y% |
0 |
11,5 |
21,4 |
38 |
51,1 |
61,9 |
71,2 |
79 |
85,4 |
91 |
95,9 |
100 |
t0C |
110,6 |
108,3 |
106,1 |
102,2 |
98,6 |
95,2 |
92,1 |
89,4 |
86,8 |
84,4 |
82,3 |
80,2 |
По этим данным строим кривые в координатах x – y (рис. 1) и t – x, y (рис. 2).
Минимальное флегмовое число Rmin определяем по уравнению:
где xD – концентрация бензола (ЛКК) в дистилляте, мол. доли;
xF – концентрация бензола в сырье, мол. доли;
y*F – мольная доля бензола в паре, равновесном с жидкостью, состава хF; мол. доли.
- равновесная концентрация
Тогда минимальное флегмовое число равно:
Действительное флегмовое
2.3. Построение рабочих линий на диаграмме х – у.
Мольный расход сырья определяем по уравнению:
кмоль/час
Расход дистиллята находим следующим образом:
кмоль/час
Расход остатка определим по уравнению:
кмоль/час
Относительный мольный расход питания находим по уравнению:
Уравнение рабочей линии верхней (укрепляющей) части колонны:
Уравнение рабочей линии нижней (исчерпывающей) части колонны:
3. Определение скорости пара и диаметра колонны
Средние концентрации жидкости
а) в верхней части колонны:
б) в нижней части колонны:
Средние концентрации пара находим по уравнениям рабочих линий:
Средние температуры пара определяем по диаграмме t – x,y (рис. 2) в зависимости от средних концентраций пара.
Средняя температура пара в колонне:
Средние мольные массы паров в верхней и нижней частях колонны:
г/моль
г/моль
Определяем плотность пара в верхней и нижней частях колонны:
кг/м3
кг/м3
Средняя плотность пара в колонне:
кг/м3
Температура вверху колонны при уD = 0,966 равняется 82,0 оС, а в кубе-испарителе при хW = 0,025 она равна 109,4 оС.
Плотность жидкого бензола при 82,0 оС
кг/м3
Плотность жидкого толуола при 82 оС:
кг/м3 (табл. IV, стр. 512, [2])
Плотность дистиллята равна:
кг/м3
Плотность жидкого бензола при 109,4 оС
кг/м3
Плотность жидкого толуола при 109,4 оС:
кг/м3 (табл. IV, стр. 512, [2])
Плотность остатка равна:
кг/м3
Определим среднюю плотность жидкости в колонне:
кг/м3
Объемный расход проходящего через колонну пара V при его средней температуре в колонне, определяемой как среднее арифметическое, находится по уравнению:
где Т и Т0 – абсолютные температуры соответственно средняя пара в колонне и температура при нормальных условиях,
D – мольный расход дистиллята, равен 64,95 кмоль/ч (см. стр. 5, п. 2.3.)
Тогда объемный расход пара равен:
м3/с
Диаметр колонны может быть определен по формуле:
Принимаем расстояние между тарелками h = 500 мм. Для ситчатых тарелок по графику (рис. 7.2., стр. 323, [3]) находим С = 0,075. Скорость пара в колонне определяем по уравнению:
м/с
Тогда диаметр колонны равен:
м
Приняв стандартный размер обечайки равным 1,4 м уточним рабочую скорость пара.
м/с
4. Определение гидравлического сопротивления тарелок и давления в кубе колонны
Выбираем тарелку из ряда стандартных. Выбираем ситчатые тарелки ТС-Р с приведенными ниже характеристиками (стр. 216, [2])
(ОСТ 26-01-108-85)): Диаметр тарелки – 1400 мм; |
||||||
Свободное сечение колонны – 1,54 м2; |
||||||
Сечение перелива – 0,087 м2; Периметр слива – 0,86 м; |
||||||
Относительное свободное сечение тарелки – 10,7 %; | ||||||
Масса – 72 кг | ||||||
Рассчитаем гидравлическое сопротивление тарелки в нижней и верхней части колонны по уравнению
где - сопротивление сухой тарелки,
- сопротивление, вызываемое
силами поверхностного
пж - сопротивление парожидкостного слоя на тарелке.
а) Верхняя часть колонны.
Гидравлическое сопротивление сухой тарелки
=
где - скорость пара в прорезях;
- коэффициент сопротивления, равный для ситчатых тарелок со свободным сечением 7 – 10% (см. [2] стр. 354). Принимаем = 1,82.
- средняя плотность пара в верхней части колонны
Находим скорость пара в отверстиях тарелки . Скорость пара в колонне
м/с
Тогда
где Fсв = 10,7% - свободное сечение тарелки.
Гидравлическое сопротивление сухой тарелки равно
Сопротивление, обусловленное силами поверхностного натяжения
где s - поверхностное натяжение жидкости, dэ – эквивалентный диаметр отверстия, м. Для ситчатой тарелки равен диаметру отверстия.
При поверхностное натяжение бензола равно 20,36·10-3 Н/м, толуола – 20,71·10-3 Н/м (см. [2] стр.527). Тогда
Н/м
(см. [2] стр.527)).
=20,45 Па
Сопротивление парожидкостного слоя на тарелке:
Δрпж = 1,3 · hпж ·ρжср.в. · g · k
где g – ускорение силы тяжести, м/с2
k – относительная плотность пены (при расчетах принимают k = 0,5)
высота парожидкостного слоя на тарелке:
hпж = hп + Δh
величину Δh – высоту слоя над сливной перегородкой рассчитываем по формуле:
где Vж – объемный расход жидкости, м3/с, П – периметр сливной перегородки, м
Объемный расход жидкости в верхней части колонны:
м3/с
где г/моль
Периметр слива равен 0,86 м.
Тогда
м
Высота парожидкостного слоя на тарелке
hпж = 0,05 + 0,025 = 0,075 м
Сопротивление парожидкостного слоя на тарелке вверху колонны:
Δрпж = 1,3 · 0,075 ·812,5 · 9,81 · 0,5 = 388,57 Па
Общее гидравлическое сопротивление
тарелки в верхней части
б) Нижняя часть колонны.
Гидравлическое сопротивление сухой тарелки
При поверхностное натяжение бензола равно 18,46·10-3 Н/м, толуола – 19,11·10-3 Н/м (см. [2] стр.527). Тогда
Н/м
(см. [2] стр.527).
Сопротивление, обусловленное силами поверхностного натяжения.
=18,92 Па
Объемный расход жидкости в нижней части колонны:
м3/с
где г/моль
Тогда
м
Высота парожидкостного слоя на тарелке
hпж = 0,05 + 0,045 = 0,095 м
Сопротивление парожидкостного слоя на тарелке внизу колонны:
Δрпж = 1,3 · 0,095 ·782,9 · 9,81 · 0,5 = 474,26 Па
Общее гидравлическое сопротивление тарелки в нижней части колонны 205,02 + 18,92 + 474,26 = 698,2 Па
Проверим, соблюдается ли при расстоянии между тарелками h = 0,5 м необходимое для нормальной работы тарелок условие
Для тарелок нижней части колонны, у которых гидравлическое сопротивление больше
= =0,163 < 0,500 м
Следовательно, условие соблюдается.
Общее гидравлическое сопротивление тарелок в колонне
Р = 606,52 · 11 + 698,2 · 16 = 17842,92 Па
Давление в кубе колонны равно:
Робщ = 101325 + 17842,92 = 119 167,92 Па
.
5. Определение действительного числа тарелок и высоты колонны
Наносим на диаграмму х – у рабочие линии верхней и нижней части колонны (см. рис. 1) и находим число ступеней изменения концентрации пт. Их число соответствует количеству теоретических тарелок и равно для верхней части 6, для нижней части - 8 . Число тарелок рассчитываем по уравнению:
Для определения среднего к.п.д. тарелок
η находим коэффициент
При этой температуре Рб = 1204 мм рт. ст - давление насыщенного пара бензола, Рт = 492,5 мм рт. ст. - давление насыщенного пара толуола (стр. 565, [2]). Тогда
a = 1204/492,5 = 2,45.
Вязкость жидкости на питающей тарелке при температуре 95,15 оС определяем по формуле:
lnm = х × lnmб + (1 – x) × lnmт
где mб = 0,274 ·10-3 Па×с – вязкость бензола при температуре 95,15 оС (стр. 516, [2]), mт = 0,283 ·10-3 Па×с – вязкость толуола при температуре 95,15 оС.
Тогда ln mж = 0,541 ×ln 0,274 + (1 – 0,541) × ln 0,283
откуда mж = 0,278 ·10-3 Па×с
Произведение am = 2,45 × 0,278= 0,681
По произведению am = 0,681 находим к.п.д. тарелки h=0,53 (стр. 323, [2]).
Тогда число тарелок в верхней части колонны
n' = 6/0,53 ≈ 10
в нижней части колонны:
n'' = 8/0,53 ≈ 15
Общее число тарелок Nобщ = 10 + 15 = 25, с запасом 27, из них 11 в верхней части колонны, 16 – в нижней части.
Определим высоту колонны по формуле: , где
– расстояние между тарелками,
– число практических тарелок.
Тогда в верхней части колонны
в нижней части колонны
Но фактическая высота будет больше за счет следующих элементов:
1.Высота от верхнего днища до первой тарелки:
2.Ввод сырья, из расчета расстояния между 2 тарелками - 1,0 м.
3.Расстояние от нижней тарелки до куба колонны - 1,5 м.
4. - десятиминутный запас колонны, принимаем 3,0 м.
- Высоту юбки примем 2 м.
- Расстояние для установки люков – 0,8 м, количество люков – 2, общая высота 2 ∙ (0,8 – 0,5) = 0,6 м.
Общая высота колонны: 1,7 + 5,0 + 1,5 + 7,5 + 1,5 + 3,0 + 2 +0,6 = 22,8 м.
Принимаем высоту колонны 23 м.
6. Тепловой баланс колонны
Количество тепла, получаемой в кубе-испарителе от греющего пара, находим по уравнению:
Количество тепла, отнимаемого водой в дефлегматоре, определяем по уравнению:
Вт
где теплоту конденсации дистиллята определяем по правилу аддитивности:
Дж/кг
где rб = 392,4 · 103 Дж/кг - теплота конденсации бензола при температуре верха колонны, равной 82оС (табл. ХLV, стр. 542 [2]), rт = 377,8 · 103 Дж/кг - теплота конденсации толуола при температуре верха колонны, равной 82 оС
Значение теплоты, отводимой с дистиллятом, определяем по уравнению:
где GD – расход дистиллята, кг/с,
сD – удельная теплоемкость дистиллята, Дж/(кг·оС)
tD – температура верха колонны, оС.
Удельную теплоемкость дистиллята рассчитывается по правилу аддитивности:
Дж/(кг·оС)
где сб= 1906,45 Дж/(кг·оС) – удельная теплоемкость бензола при температуре верха колонны, равной 82 оС (рис. ХI, стр. 562 [2]), cт = 1851,98 Дж/(кг·оС) – удельная теплоемкость толуола при температуре верха колонны, равной 82 оС.
Вт
Значение теплоты, отводимой с кубовым остатком, определяем по уравнению:
где GW – расход кубового остатка, кг/с,
сW – удельная теплоемкость кубового остатка, Дж/(кг·оС)
tW – температура кубового остатка, оС.
Удельную теплоемкость кубового остатка также рассчитывается по правилу аддитивности:
Дж/(кг·оС)
где сб= 2053,1 Дж/(кг·оС)– удельная теплоемкость бензола при температуре низа колонны, равной 109,4 оС (рис. ХI, стр. 562 [2]), cт = 1885,5 Дж/(кг·оС) – удельная теплоемкость толуола при температуре низа колонны, равной 109,4 оС.
Вт
Значение теплоты, вносимой в колонну с сырьем, определяем по уравнению:
где GF – расход сырья, кг/с,
сF – удельная теплоемкость сырья, Дж/(кг·оС)
tF – температура сырья, оС.
Удельная теплоемкость сырья равна:
Дж/(кг·оС)
где сб = 1948,35 Дж/(кг·оС)– удельная теплоемкость бензола при температуре ввода сырья, равной 91 оС (рис. ХI, стр. 562 [2]), cт = 1822,65 Дж/(кг·оС) – удельная теплоемкость толуола при температуре ввода сырья колонны, равной 91 оС.
Вт
Принимаем потери тепла равными 3 % от полезного количества теплоты
Вт
Определим расход горячей
струи, подаваемой в куб
Теплоту парообразования кубового остатка рассчитываем по правилу аддитивности:
кДж/кг
где rб = 371,75 кДж/кг – теплота парообразования бензола при температуре кипения кубового остатка, равной 109,4 оС (табл. ХLV, стр. 541 [2]), rт = 363,14 кДж/кг – теплота парообразования толуола при температуре 109,4 оС.
Расход горячей струи равен:
кг/с
7. Расчет и подбор теплообменного оборудования.
а) Подогреватель исходной смеси.
Исходная смесь, расход которой равен GF = 2,777 кг/с нагревается от температуры t2н = 18 °C до t2к = 91°C.
Удельная теплоемкость сырья при средней температуре (91 + 18)/2 = 54,5 оС равна:
Дж/(кг·оС)
где сб = 1801,7 Дж/(кг·оС) – удельная теплоемкость бензола при температуре, равной 54,5 оС (рис. ХI, стр. 562 [2]), cт = 1759,8 Дж/(кг·оС) – удельная теплоемкость воды при температуре 54,5 оС.
Тепловая нагрузка аппарата:
Вт
Тогда расход водяного пара составит (абсолютное давление водяного пара принимаем равным 2 атм):
кг/с = 619,56 кг/ч
где r = 2208 кДж/кг·град – удельная теплота конденсации пара при давлении 2 атм. (по табл. LVII, стр. 549, [2]).
Температура конденсации греющего пара при давлении Р = 2 атм равна 119,6 оС (табл. LVII, стр. 550, [2]). Составим температурную схему для процесса:
119,6 оС → 119,6 оС
91 оС ← 18 оС
___________________
Δtм = 28,6 град Δtб = 101,6 град
Средняя разность температур равна:
Примем коэффициент
Рассчитаем площадь
По табл. 2.3., стр. 51, [3] выбираем одноходовой теплообменник с диаметром кожуха 0,4 м, длина труб 3,0 м, поверхность теплообмена 26 м2.
Запас площади поверхности
б) Расчёт кожухотрубчатого конденсатора (дефлегматора).
Количество тепла, отнимаемого водой в дефлегматоре (см. п. 6):
Вт
температура конденсации tk = 82 °С.
Тепло конденсации отводим водой с начальной температурой t2н= 17°С.
Составим температурную схему для процесса:
82 оС → 82 оС
37 оС ← 17 оС
___________________
Δtм = 45 град Δtб = 65 град
Следовательно, средняя разность температур равна:
Примем ориентировочно Kор= 400 Вт/м2·К(табл. 4.8, стр. 172, [2]).
Рассчитаем ориентировочное значение
требуемой поверхности
По табл. 2.9., стр. 57, [3] выбираем двухходовой теплообменник с диаметром кожуха 0,6 м, длина труб 4,0 м, поверхность теплообмена 75 м2.
Запас площади поверхности
в) Расчет кипятильника.
Тепловая нагрузка аппарата (см. п. 6):
QК= 1572225 Вт
Расход пара для подогрева горячей струи:
кг/с
где r = 2117 кДж/кг·град – удельная теплота конденсации пара при давлении 5 атм. (по табл. LVII, стр. 549, [2]).
Температура конденсации греющего пара при давлении Р = 5 атм равна 151,1 оС (табл. LVII, стр. 550, [2]).
Следовательно, средняя разность температур равна:
Примем коэффициент
Рассчитаем площадь
Задаваясь числом Re2= 15000, определим соотношение n /z для конденсатора из труб диаметром dн= 25 ´ 2 мм:
где n – общее число труб;
z – число ходов по трубному пространству:
d – внутренний диаметр труб, м,
μ - вязкость жидкости, определяем по формуле:
lnm = х × lnmб + (1 – x) × lnmт
где mб = 0,288 мПа×с – вязкость бензола (стр. 514, [2])
mт = 0,295 мПа×с – вязкость толуола при 90,4 оС (стр. 516, [2])
Тогда
ln mж = 0,025 ×ln 0,288 + (1 – 0,025) × ln 0,295
откуда mж = 0,293 мПа×с
В соответствии с табличными значениями (табл. 2.3., стр. 51, [3]) соотношение n /z принимает наиболее близкое к заданному значению у теплообменника с диаметром кожуха D = 400 мм, диаметром труб 20´ 2 мм, числом ходов z = 1 и общим числом труб n = 181.
По табл. 2.3., стр. 51, [4] выбираем одноходовой теплообменник с диаметром кожуха 0,40м, длина труб 6,0 м, поверхность теплообмена 68 м2.
Рассчитаем действительное число Re:
Для смеси при температуре 109,4 оС находим критерий Прандтля по формуле
где λ – коэффициент теплопроводности смеси, определяем по формуле 4.8., стр. 150, [2]:
где ε – температурный коэффициент, для бензола и толуола
ε = 1,8·10-3 оС-1
λ30 рассчитываем по формуле:
где с – удельная теплоемкость жидкости, Дж/кг·К;
ρ – плотность жидкости, кг/м3;
М – мольная масса жидкости, кг/кмоль;
А – коэффициент, зависящий от степени ассоциации жидкости,
А = 4,22 · 10-8.
Тогда при 109,4 оС коэффициент теплопроводности смеси равен:
По графику 4.1., стр. 154, [2] находим критерий Нуссельта.
Для Re = 6494 значение . Для нагревающихся жидкостей в расчете коэффициентов теплоотдачи можно принимать
(стр. 152, [2]).
Тогда получим Nu = 22,5 · 5,040.43 ·1,1 = 49,62
Тогда коэффициент теплоотдачи кубовой смеси равен:
Коэффициент теплоотдачи конденсирующегося водяного пара рассчитываем по уравнению (стр. 162, [2]):
Значения коэффициента Аt при температуре 151,1 оС определяем по табл. 4.6., стр. 162, [2].
Аt = 7455
Значение ε принимаем равным 1. Подставляя численные значения, получим:
Сопротивление стенки и загрязнений определяем по формуле:
где δст – толщина стенки трубы, равная 0,002 м,
λст – теплопроводность стали, равная 46,5
rзагр1 – сопротивление загрязнений со стороны кипящей жидкости,
rзагр2 – сопротивление загрязнений со стороны конденсирующегося пара.
По табл. ХХХI, стр. 531 [2] принимаем 1/rзагр1 = 5800 и 1/rзагр1 = 5800
Тогда
м2·град/Вт
Коэффициент теплопередачи определяем по уравнению:
Требуемую площадь теплообменника определяем по формуле:
Рассчитаем площадь
Проверим принятое допущение
Для этого определим температуру стенки трубы теплообменника со стороны органической жидкости. Удельная тепловая нагрузка равна:
Температура стенки со стороны смеси равна:
Для смеси при температуре 144,82 оС находим критерий Прандтля
Удельная теплоемкость кубового остатка при температуре 144,82 оС равна: