Разделение смеси жидкостей на составляющие компоненты
Содержание
Введение…………………………………………………………
1. Технологическая схема
и ее описание……………………………………………
2. Технологический расчет………………………………………………………………
2.1. Материальный баланс
колонны…….…………………………………………….……
2.2.Определение минимального
флегмового числа…..……………………………
2.3.Определение давления
в колонне………………………………………………………
2.4.Определение
2.5.Средние массовые расходы по жидкости и пару…………………………………...…... 12
2.6.Тепловой баланс колонны ………………………………………………………………… 13
3.Гидравлический расчет…………………
4. Расчет штуцеров……………………………………………….………
Заключение …………………………………………………
Список использованной литературы ……………………………….………………………..…22
Введение
Ректификация - массообменный процесс разделения однородной смеси летучих компонентов, осуществляемый путем противоточного многократного взаимодействия паров, образующихся при перегонке, с жидкостью, образующейся при конденсации этих паров.
Разделение жидкой смеси основано на различной летучести веществ. При ректификации исходная смесь делится на две части: дистиллят - смесь, обогащенную низкокипящим компонентом, и кубовый остаток - смесь, обогащенную высококипящим компонентом.
Процесс ректификации осуществляется в ректификационной установке основным аппаратом которой является ректификационная колонна, в которой пары перегоняемой жидкости поднимаются снизу, а навстречу парам стекает жидкость, подаваемая в виде флегмы в верхнею часть аппарата.
Процесс ректификации может
протекать при атмосферном
Степень разделения смеси
жидкостей на составляющие компоненты
и чистота получаемого
В промышленности применяют тарельчатые, насадочные, пленочные трубчатые и центробежные пленочные аппараты. Они различаются в основном конструкцией внутреннего устройства аппарата.
Тарельчатые контактные устройства, в ректификационной колонне, подразделяются на: колпачковые, ситчатые, клапанные, решетчатые, провальные и т.п.
Преимущественное
Технологическая схема и ее описание
Рис. 1 - Принципиальная схема ректификационной установки
1 – емкость для исходной смеси; 2 – подогреватель; 3 – колонна;
4 – кипятильник; 5 – дефлегматор; 6 – делитель флегмы; 7 – холодильник; 8 – сборник дистиллята; 9 – сборник кубового остатка.
Исходную смесь из емкости 1 центробежным насосом подают в теплообменник 2, где она подогревается до температуры кипения. Нагретая смесь поступает на разделение в ректификационную колонну 3, где состав жидкости равен составу исходной смеси xF . Стекая вниз по колонне, жидкость взаимодействует с поднимающимся вверх паром, образующимся при кипении кубовой жидкости в кипятильнике 4. В результате массообмена с жидкостью пар обогащается легколетучим компонентом. Для более полного обогащения верхнюю часть колонны орошают в соответствии с заданным флегмовым числом жидкостью (флегмой) состава хР, получаемой в дефлегматоре 5путем конденсации пара, выходящего из колонны. Часть конденсата выводится из дефлегматора в виде готового продукта разделения - дистиллята, который охлаждается в теплообменнике 7и направляется в емкость 8. Из кубовой части колонны насосом непрерывно выводится кубовая жидкость – продукт, обогащенный труднолетучим компонентом, который охлаждается в теплообменнике 7 и направляется в емкость 9.
Таким образом, в ректификационной
колонне осуществляется непрерывный
неравновесный процесс
Технологический расчет
Материальный баланс процесса
Составляем материальный баланс для определения количеств и состава веществ, участвующих в процессах ректификации.
Материальный баланс колонны, обогреваемой паром:
,
где L—производительность установки по исходной смеси, D –производительность установки по дистилляту, R— производительность установки по кубовому остатку.
Материальный баланс для НК:
,
где хF, xД, хW— массовая доля легколетучего компонента в исходной смеси, дистилляте, кубовом остатке соответственно. Откуда
D=L(xL-xR)/(xD-xR)=15000(0,5-
R=L(xL-xD)/(xR-xD)=15000(0,95-
Для дальнейших расчётов выразим концентрации исходной смеси, дистиллята и кубового остатка в мольных долях:
1),
где - молекулярная масса метанола кг/кмоль), - молекулярная масса воды (Мводы=18 кг/кмоль).
2)
3)
Определение минимального флегмового числа
Для определения минимального флегмового числа строим кривую равновесия, предварительно выполнив расчет равновесного состава жидкости и пара смеси метанол—вода.
Таблица 1. Равновесный состав жидкости и пара смеси метиловый спирт—вода.
По данным табл. 1 построим диаграмму равновесия x,y:
Рис. 2 Диаграмма равновесия между паром (у) и жидкостью (х) при постоянном давлении.
Минимальное флегмовое число
Рабочее флегмовое число R
,
где β – коэффициент избытка флегмы (берем произвольно)
; ;
; .
Рассчитаем число B:
; ;
; .
Построим диаграммы равновесия x,y. На диаграммах отложим значения В, затем построим рабочие линии укрепляющей и исчерпывающей части колонны и нанесем линии обозначающие теоретические тарелки. По количеству пиков, определим число теоретических тарелок (Nт).
Рис. 3 - Диаграмма равновесия между паром (х) и жидкостью (у) при флегмовом числе R1=0,9.
Рис. 4 - Диаграмма равновесия между паром (х) и жидкостью (у) при флегмовом числе R2=1,07.
Рис. 5 - Диаграмма равновесия между паром (х) и жидкостью (у) при флегмовом числеR3=1,25.
Рис. 6 - Диаграмма равновесия между паром (х) и жидкостью (у) при флегмовом числеR4=1,29.
Рис. 7 - Диаграмма равновесия между паром (х) и жидкостью (у) при флегмовом числе R5=1,36.
Рис. 8 - Диаграмма равновесия между паром (х) и жидкостью (у) при флегмовом числе R6=1,44.
Полученное количество теоретических тарелок занесем в табл. 2:
Таб. 2 – Данные для нахождения Rопт
|
|||||||||||||||||||||||||||||||||||||||||
|
| |||||||||||||||||||||||||||||||||||||||||
где – минимальное произведение, пропорциональное объему ректификационной колонны (– число ступеней изменения концентраций или теоретических тарелок, а – расход паров).
По данным табл. 2 построим график зависимости:
Рис. 9 – График зависимости Nтот Nт·(R-1)
Из рис. 7 находим, что Rопт=1,25. Итак, число теоретических тарелок в верху колонны n=2, внизу m=6.
Определение давления в колонне
В данной работе в верхнем
продукте содержится 97,4% (мольных)
метанола, поэтому температура кипения
такой смеси будет
Поскольку при давлении в
Принимаем давление верха
πз.п= 112 + 9∙0,35 = 115,1 кПа;
πн = 115,1 + 9∙0,52 = 119,78 кПа.
Средние давленияв верхней и нижней частях колонны:
πср.в=(πз.п+πД)/2=(115,1+112)/
πср.н=(πз.п+πн)/2=(115,1+119,
Определение температурного режима в колонне
Давления вверху и внизу
Рис.8. Изобарные температурные кривые кипения и конденсации смеси метанол-вода
Так как сырье поступает при
температуре кипения (е=0), для
нахождения температуры tLдоста
Средние температуры низа и верха колонны составляют:
tв= (tD+tL)/2 = (70,9+81,7)/2=76,30С;
tн= (tR+tL)/2 = (99,2+81,7)/2=90,450С.
Средние массовые расходы по жидкости и пару
Т.к доля отгона е=0,паровой поток G0,вносимый с сырьем равен 0.
Жидкий поток, вносимый в колонну с сырьем g0=GF=15000(кг/ч)
Массовый расход пара Gв в верху колонны:
Т.к сырье поступает в жидкой фазе,Gm=G=Gв=17419(кг/ч)
Массовый расход флегмы, поступающей в зону питания:
Массовый расход флегмы, поступающей в нижнюю часть колонны:
Массовая концентрация нкк, поступающего на нижнюю тарелку, определяется соотношением:
где ,
Массовая концентрация нкк в паровом потоке, поднимающемся с последней тарелки нижней части колонны в зоне питания:
Массовая концентрация нкк в паровом потоке, покидающем зону питания:
Тепловой баланс колонны
Принимаем температуру
Тепловой поток, отводимый водой в дефлегматоре, рассчитывается по уравнению
QД=GB[rD+ cD(tD– tx)],
где rD - теплота конденсации паров, поступающих в дефлегматор;cD – теплоемкость жидкого дистиллята (флегмы);tD,tx – температуры верха колонны и флегмы, поступающей на орошение.
При этом средние значения удельной теплоты испарения rDи удельной теплоемкости cDнаходятся по правилу аддитивности:
rD=0,05∙2345,2+0,95∙1097,9 =1159,73кДж/кг;
cD=0,05∙4,18+0,95∙2,719 = 2,79 кДж/(кг∙К);
QД=17419,28[1159,73+ 2,79(70,9– 40)]/3600 = 6029,03 кВт.
Энтальпии сырья iL, дистиллята iD, кубовой жидкости iR определяются по правилу аддитивности при соответствующей температуре:
iL=(0,5∙4,19+0,5∙2,81)∙81,7=
iD=(0,05∙4,18+0,95∙2,76)∙70,9=
iR=(0,98∙4,18+0,02∙2,94)∙99,2=
С учетом тепловых потерь, принятых
равными 5% от полезно используемого
расхода теплоты, тепловой
Qкуб = 1,05∙( QД+ GDiD+GRiR–
GLiL)=1,05∙(6029+7741,9∙200,
В качестве теплоносителя в
кипятильнике колонны
GТ =Qкуб/((Iт.н. - Iт.к.)∙φ )= 6415,56/((2730-558,9)∙0,95) = 3,11 кг/с.
Принимаем, что вода в
Gх=QД/(iх.н.– iх.к.) = 6029,03/(4,18∙(38-25))=110,95 кг/с =0,173м3/с=622,8м3/час.
Массовый расход холодного
gх= (gп∙ rD)/( ItD – itx),
где ItD– энтальпия пара, поступающего в дефлегматор при температуре tD; itx– энтальпия жидкости при температуре tx.
ItD=0,05∙(2345,9+4,18∙70.9)+0,
itx= (0,05∙4,18+0,95∙2,6)∙40=107.2 кДж/кг.
gх= (9677.37∙ 1159.73)/( 1360.4 – 107.2) = 8954.9 кг/час.
Массовый расход горячего орошения g1 , стекающего с 1-й тарелки верхней части колонны
g1 ≈ gх+ gх∙
сх(tD - tx)/rD=8954.9
+ 8954.9∙(0,05∙4,18+0,95∙2,7)∙(
Проверка: g1≈gn; 9616.8≈9677.3
Табл.3. Теплофизические свойства воды и метилового спирта
свойство |
Вода |
Спирт | ||||||||||
Температура, 0С |
Температура, 0С | |||||||||||
40 |
65,4 |
71,5 |
73,3 |
89,6 |
98,8 |
40 |
65,4 |
71,5 |
73,3 |
89,6 |
98,8 | |
Удельная теплоемкость с, кДж/(кгК),[5, с.520,544] |
4,18 |
4,18 |
4,18 |
4,19 |
4,19 |
4,23 |
2,6 |
2,72 |
2,77 |
2,77 |
2,85 |
2,93 |
Давление насыщенного пара Р, мм. рт. ст.[5, с.519,547]. |
249 |
525 |
800 |
1600 |
||||||||
Удельная теплота |
2345,2 |
1097,3 |
||||||||||
Плотность жидкости ρ, кг/м3,[5, с.520,7,т.1,с.555]. |
992 |
981,7 |
977,1 |
976 |
965,3 |
958,8 |
745 |
750 |
744,4 |
742,5 |
725,4 |
715,3 |
Поверхностное натяжение σ, дин/см,[7,т.1,с.1010,1015]. |
69,65 |
65,22 |
64,14 |
60,81 |
59,07 |
20 |
~14,6 |
|||||
Вязкость пара µп, сП,[5, с.539,7,т.1,с.1003]. |
0,011 |
0,012 |
0,013 |
~0,007 |
~0,009 |
|||||||
Вязкость жидкости µж, сП, [5, с.538]. |
0,66 |
0,41 |
0,32 |
0,47 |
0,31 |
0,25 |
||||||
Гидравлический расчет
При средней температуре верхней части 76,3°С, плотности жидких метанола и воды равны соответственно: и .
При средней температуре нижней части 90,45°С, плотности жидких метанола и воды равны соответственно: и .
Средние мольные и массовые концентрации нкк в верху колонны:
В нижней части колонны:
Средние мольные концентрации нкк в паре-находим по уравнениям рабочих линий
для верха
для низа ,где
yв=0,55∙0,635+0,4=0,75
yн= 1,69∙0,185-0,0114=0,3
Средние молярные массы пара в верхней и нижней частях колонны:
Mв=32∙0,75+18∙0,25=28,5(кг/
Мн=32∙0,3+18∙0,7=22,2(кг/
Средние давления вверхней и нижней частях колонны:
Пв=(112+115)/2=113,5(кПа)
Средние плотности пара в верхней и нижней частях колонны:
Cредние плотности флегмы:
Cредние обьемные расходы пара:
Максимально допустимая линейная скорость пара:
Диаметр верхней части колонны
Диаметр нижней части колонны
Выбор диаметра колонны
Из двух найденных диаметров и выбираем максимальный (), и в соответствии с этим выбираем из принятого в химической промышленности нормального ряда диаметров колонн [2, с. 122].
Действительные рабочие скорости паров
Для нашего процесса оптимальными будут ситчатые тарелки-т.к. они отличаются простотой устройства,легкостью монтажа, осмотра и ремонта. Гидравлическое сопротивление этих тарелок невелико. Ситчатые тарелки устойчиво работают в широком диапазоне скоростей газа.
По ГОСТ 26-805-73 принимаем следующие размеры ситчатой тарелки: диаметр отверстий dо= 4 мм, высота сливной перегородки hп = 40 мм; свободное сечение тарелки (суммарная площадь отверстий) – 10% от общей площади тарелки; площадь, занимаемая двумя сегментными переливными устройствами, - 20% от общей площади тарелки.
Расчет гидравлического
сопротивления тарелок
Верхняя часть колонны
Скорость пара в отверстиях тарелок ωотв=ωв/0,12 = 1,4/0,12=11,6 м/с.
Коэффициент сопротивления неорошаемых ситчатых тарелок I=1,75[5, с.354]. Тогда гидравлическое сопротивление сухой тарелки
==129,5 Па.
Сопротивление, обусловленное силами поверхностного натяжения,
ΔРσ=4σ/d0==4∙14,6∙0,001/0,004 =14,6 Па.
Периметр сливной перегородки П=1,32 м [5, с.354].
Принимаем отношение плотности парожидкостного слоя(пены) на тарелке к плотности жидкости к′ =0,5.
Средний объемный расход флегмы в верхней части колонны
Vф,в= gn/ρф,в=9677,37/(3600∙797,4)=
Высота слоя над сливной перегородкой
Δh=( Vх,в/(1,85∙П∙к′))2/3=(
0,0034/(1,85∙1,32∙0,5))2/3=0,
Высота парожидкостного слоя на тарелке
hп.ж.=hп.+Δh=0,04+0,02=0,06 м.
Сопротивление парожидкостного слоя
ΔРп.ж.=1,3∙hп.ж.∙ρф,в∙к′∙g=1,
Общее гидравлическое сопротивление
тарелки в верхней части
ΔРв =ΔРсух +ΔРσ +ΔРп.ж.= 129,5+14,6+254,2= 398,3 Па.
Нижняя часть колонны
Скорость пара в отверстиях тарелок ω′отв=ωн/0,12 = 1,88/0,12=15,6м/с.
Гидравлическое сопротивление сухой тарелки
==123 Па.
Сопротивление, обусловленное силами поверхностного натяжения,
ΔР′σ=4σ/d0==4∙60,77∙0,001/0,
Средний объемный расход жидкости в нижней части колонны
V′ф,н= g/ρф,н==24677,38/(3600∙885)=0,
Высота слоя над сливной перегородкой
Δhп.ж′=( V′х,н/(1,85∙П∙к′))2/3=( 0,077/(1,85∙1,32 ∙0,5))2/3=0,034м.
Высота парожидкостного слоя на тарелке
h′п.ж.=hп.+Δh′=0,04+0,034=0,
Сопротивление парожидкостного слоя
ΔР′п.ж.=1,3∙h′п.ж.∙ρх,в∙к′∙g=
Общее гидравлическое сопротивление тарелки в нижней части колонны
ΔР′н=ΔР′с+ΔР′σ+ΔР′п.ж.= 123+60,8+417,6= 601,4 Па.
Проверим, соблюдается ли при h=0,4 м необходимое условие для нормальной работы тарелок h>1,8∙ΔР/(ρжg)
Для тарелок верхней части колонны: 1,8∙398,3/(797,4∙9,81)=0,09.
Для тарелок нижней части колонны: 1,8∙601,4/(885∙9,81)=0,12.
Следовательно, вышеуказанное
условие соблюдается и
Минимальная скорость пара в отверстиях ωотв.min, достаточная для того, чтобы ситчатая тарелка работала всеми отверстиями:
Для верхней части колонны
ωотв.min=0,67∙(g∙ ρх.в.∙hп.ж/ξρв)1/2=0,67∙(9,81∙ 797,4.∙0,06/1,75∙1,1)1/2 =10,46 м/с;
Для нижней части колонны
ω′отв.min=0,67∙(g∙ ρх.н.∙h′п.ж/ξρн)1/2=0,67∙(9,
Фактические скорости в отверстиях тарелок ωотви ω′отв больше минимальных ωотв.minи ω′отв.min, следовательно, все отверстия тарелок будут работать равномерно.
Величина межтарельчатого уноса жидкости еуопределяется по уравнению
еу=А(0,052∙hσ-1,72)/(hv∙ψ2) ∙ (ω/εэф∙ S)3,7 ,
где А,v, ψ-коэффициенты; для h≥ 0,4 м , А=0,159, v=0,95, ψ=0,9[2, с.194]; hσ- глубина барботажа; для ситчатых тарелок hσ≈ hп=0,04 м; ω- скорость пара в колонне; εэф– эффективная рабочая площадь тарелки, вычисляемая по графику[2, с. рис.111-7]; для отношения П/D=1,76/2,4=0,73, εэф=0,65.
Величина S определяется из уравнения
S=1,15∙10-3 ∙(σ/ρп)0,295 ∙ ((ρж-ρп)/µп)0,425,
где σ,ρж – поверхностное натяжение и плотность жидкости; ρп,µп – плотность и вязкость пара.
Рассчитаем величину межтарельчатого уноса жидкости еуотдельно для верхней и нижней части колонны.
Верхняя часть колонны
Средняя вязкость паров µп.в.:
Мв/µп.в.=увМм/µм + (1-ув)Мвода/µвода ,
где µми µвода – вязкости метанола и воды.
Откуда µп.в=29,34/(0,81∙32/0,0086+(1-
S=1,15∙10-3 ∙(14,6∙10-3 /1,1)0,295 ∙ ((797-1,1)/0,009∙10-3 )0,425=0,77;
еу=0,159∙(0,052∙40-1,72)/(4000
Нижняя часть колонны
Средняя вязкость паров µп.н.:
Мв/µп.н.=унМм/µм + (1-ун)Мвода/µвода ,
где µми µвода – вязкости метанола и воды.
Откуда µп.в=23,11/(0,365∙32/0,009+(1-
S=1,15∙10-3 ∙(60,77∙10-3 /0,82)0,295 ∙ ((885-0,893)/0,0106∙10-3 )0,425=1,24;
еу=0,159∙(0,052∙40-1,72)/(4000
Приближенно определим средний коэффициент полезного действия тарелок.
Верхняя часть колонны
Коэффициент относительной
летучести для средней
Вязкость жидкости µж.в.:
lgµф.в.=0,635∙lg0,298 + (1-0,635)∙lg 0,392 , откуда µж.в.=0,33 мПа∙с;
α∙µф.в = 3,9∙0,33 =1,27.
По графику[5, с.314; 6, с.63] находим ηср.в=0,45.
Нижняя часть колонны
Коэффициент относительной
летучести для средней
Вязкость жидкости µж.н.:
lgµф.н.=0,185∙lg0,25 + (1-0,185)∙lg 0,32 , откуда µж.н.=0,31 мПа∙с;
α∙µж.н = 3,6∙0,31 =1,107.
По графику находим ηср.н=0,48.
Число реальных тарелок в верхнейn и нижнейm частях колонны составляет:
n=nт/ηср.в =2/0,45 =4;
m=mт/ ηср.н= 6/0,48 =12,5 ≈13.
Рассчитанные значения гидравлического
сопротивления тарелок ΔРв=398,
Высота колонны равна:
Н =ho+hн
+hт.н. + hз.п.
+hт.в. +hв
=1,5∙2+2,2+0,4∙(4-2)+0,4∙(13-
Рис.9. Определение высоты колонны
Назначение штуцеров:
А – выход паров с верха колонны в шлемовую трубу; В – ввод флегмы; С – ввод сырья; Е – ввод паров из кипятильника; К – вывод кубовой жидкости в кипятильник колонны.
Расчет штуцеров
Диаметры штуцеров (см. рис.9) рассчитываются по уравнению
dшт=(Vп/(0,785∙ωдоп))1/2,
где Vп – объемный расход потока в трубопроводе, при этом допустимые скорости потоков ωдоп принимаются в соответствии [1,с.19табл.3].
Плотность жидких продуктов холодного орошения ρх, сырья ρсырья, кубового остатка ρкуб.ост. рассчитываются в зависимости от температуры и состава:
1/ρх =0,95/745+0,05/992, откуда ρх =754,4 кг/м3,
1/ρсырья= 0,5/734,3+0,5/971,7, откуда ρсырья = 836,5 кг/м3,
1/ρкуб.ост. = 0,02/725+0,98/965 , откуда ρкуб.ост. =958,7 кг/м3.
Плотность паров, поступающих из кипятильника колонны ρп,R, и паров, уходящих с верха колонны ρп,D, рассчитываются при соответствующих температурах и давлениях:
ρп,R= (0,01∙32+(1-0,01)∙18)∙273∙117,
ρп,D= (0,95∙32+(1-0,95)∙18)∙273∙113,
Диаметр штуцера А для вывода паров из колонны в дефлегматор
dА=(11816,5/(0,785∙3600∙1,24∙
Принимаем по ГОСТ 12830-67 dА=500 мм[3, с.218].