Расчет абсорберов
При расчете абсорберов обычно заданы: расход газа, его начальная и конечная концентрации (иногда вместо концентраций задается степень извлечения ср), начальная концентрация абсорбента. Основными определяемыми величинами являются: расход Ь абсорбента, диаметр О и высота Н абсорбера, его гидравлическое сопротивление А р.
Расход абсорбента. Расход Ь поглотителя находится по уравнению (XI, 14а). При этом значение конечной концентрации Хк выбирают, исходя из оптимального удельного расхода абсорбента (см. стр. 438).
Гидравлическое сопротивление абсорберов. При движении через колонну газ преодолевает гидравлическое сопротивление, причем разность давлений газа на входе в абсорбер и выходе из него должна быть равна сопротивлению, оказываемому его движению. Гидравлическое сопротивление абсорбера зависит^от конструктивных особенностей аппарата' и гидродинамического режима его работы, связанного со скоростью газа. Основное влияние на величину Ар оказывает скорость газа. Оптимальную скорость газа в абсорбере можно определить только при помощи технико-экономического расчета с учетом всех величин, зависящих от да (гидравлическое сопротивление, диаметр и высота аппарата).
Если абсорбция проводится под повышенным давлением, то, как указывалось ранее, потеря напора на преодоление гидравлического.сопротивления абсорбера в данном случае составляет незначительную долю общего давления в системе и не оказывает существенного влияния на экономические показатели абсорбционной установки. При этом целесообразно использовать наибольшие возможные скорости газа в абсорбере, близкие к предельной, т. е. равной, например (0,8—0,9) и>а, где
— скорость, отвечающая точке захлебывания. В тех случаях, когда необходимо малое гидравлическое сопротивление (обычно, если абсорберы работают при атмосферном давлении или под давлением ниже атмосферного), чтобы снизить расходы на электроэнергию, для перемещения газа через абсорбер принимают более низкие скорости газа.
Однако при проектировании любого колонного аппарата следует иметь в виду, что более значительной экономии на капиталовложениях
6. Расчет абсорберов 459 удается достигнуть при уменьшении диаметра колонны, т. е. при увели- чении скорости газового-потока в аппарате. Диаметр абсорберов. Диаметр £) определяют по принятой фиктивной скорости пи0 (м/сек) газа, пользуясь общим уравнением (Х,75), где УТ — объемный расход проходящего через колонну газа, м3/сек. Высота абсорберов. Высоту Н абсорбера рассчитывают по общему уравнению массопередачи (Х,77) или (Х,77а). Например, если движу- щая сила выражена в концентрациях газовой фазы, то н, м К yaS Д Кср где М — количество поглощаемого газа; Кд — коэффициент массопередачи; а — удельная поверхность контакта фаз; 5 — площадь сечения колонны; Л Кср — средняя движущая сила процесса. ' Как указывалось в главе X, при отсутствии данных о поверхности контакта фаз высота абсорбера может быть найдена другими способами, например через объемные коэффициенты массопередачи или число единиц переноса [уравнение (X,60) Расчет пленочных абсорберов. Гидравлическое сопро- тивление рассчитывают по уравнению, аналогичному урав- нению (II,93а): . .И ^отиРг “Рпл
d9 2
Здесь Н — высота поверхности, по которой стекает пленка; d4 — эквивалентный диаметр канала, по которому движется газ; wora =• ш + шж- ср — относительная скорость газа (при противотоке); шж. ср — средняя скорость движения жидкой пленки, определяемая по уравнению (11,149) или (II, 149а); рг — плотность газа; А, — коэффициент трения.
Коэффициент трения X в общем случае зависит от критерия Рейнольдса для газа Rer и безразмерного комплекса величин --- ^р^ж., где [1Ж— вязкость жидкости и 0 — ее поверхностное натяжение. Величину X рассчитывают по уравнениям:
ReP Rer> кр
<Х1-31>
при
0,11 + 0,9 ;к’-^-Цж V
Х== (У1’32)
где Rer = —°™ критерий Рейнольдса для газа; Rer> кр — критическое значение
Мт
критерия Рейнольдса с учетом влияния на режимы движения газового потока скорости движения и физических свойств жидкой пленки.
Значение Rer_ кр определяется .зависимостью
Rer, кр —
86
0,11 + 0,9
1,19
(XI,33)
Диаметр абсорбера рассчитывают следующим образом. Для трубчатых аппаратов по принятому значению скорости гм газа в трубах сначала определяют суммарную площадь их поперечного сечения
S = ii- м?
460 Гл. XI. Абсорбция Затем, задаваясь внутренним диаметром (1 труб (обычно 0,02—0,05 м), находят их число: *-Тй1йГ <хи4> Принимая шаг = (1,25—1,5) йИ, где с?н—наружный диаметр труб и толщину бтр стенок труб, определяют диаметр абсорбера по уравнению (Х,75). Скорость газа т3, соответствующая началу захлебывания (см. стр. 446), может быть найдена по уравнению ,е (тяг •*“) -л~''п {&)' (-£•)* ,х,-35) где рж — плотность жидкости, кг/м3\ — вязкость жидкости, Мн'сек1м--, V и О' — расходы жидкости и газа, кг/сек. Уравнение (Х1,35) по структуре является общим для пленочных и насадочных абсорберов и отличается (см. стр. 438) только величиной А. В данном случае 4-0.47+1.51,^ Высота труб Н в трубчатоц абсорбере равна Я = _^1£_ где Ртр — внутренняя поверхность всех труб; сгвн — внутренний диаметр трубы. Пренебрегая толщиной стекающей пленки, можно принять, что Ртр = = ^ — поверхности контакта фаз. Тогда "— ЛЛ Подставляя вместо' Р ее выражение из уравнения (Х,46), получают Я = . Ми -г (XI,36) ПЛ &внКу р При расчете коэффициента массопередачи Ку (или Кх) коэффициент массоотдачи в газовой фазе Рг может быть определен по уравнению ^ ми;=~ Рег (рг;)¥ (XI,37) где К — коэффициент трения. Соответственно высота единицы переноса для газовой фазы 2 8<у?е°’16(Рг;)3 кт —— (XI,37а) [о,44 + 3,6 N11,, В уравнениях (XI,37) и (XI,37а): 1Ми' = —диффузионный критерий Нуссельта (0Р — коэффициент молекуляр- и ной диффузии в газовой фазе, м2/сек)', Ргг = диффузионный критерий Прандтля. Рг*Л* Коэффициент массоотдачи в жидкой фазе рж можно рассчитать по уравнению = 5Ке£(Р,;)"(-^)р где Ми’к = — диффузионный критерий Нуссельта для жидкой пленки; йэ = —^— =* = 46 — эквивалентный диаметр жидкой пленки; Яеж =■ а>ж‘ сР--'?——критерий Рей- И'Ж нольдса для жидкой пленки; Рг’ = —^ диффузионный критерий Прандтля для Рж£<ж жидкости; О ж—коэффициент молекулярной диффузии в жидкой фазе; (р*§)] 3 — приведенная толщина пленки (см. стр. 116).г
(хз,з8)
6. Расчет абсорберов 461 Значения В, т, пар для различных режимов движения пленки представлены ииже: Режим движения В т п р 1?еж<300 0,888 0,45 0,5 0,5 300 <Иеж< 1600 Реж>1600 7,7-10-5 1,0 0,5 О ^пр •' 1,21-Ю^О,909? і — 2,18 0,5 3'2- і15тКЄн о 1,47Соответственно этим режимам находят: при Иеж •< 300при 300 < Ие* С 1600
А* = °,20ебпр{?е^8-^8> (Р4)0’5 (-А.(ХГ.396)
при Иеж > 16С0
Аж = 3250бпр (Р4)0-5 (Х1,39в)
Расчет насадочных абсорберов. Для расчета гидравлического сопротивления аппарата, предварительно определяется сопро- тивление сухой насадки
Лрсух = ?.-£--^ (XI.40)
4?
Здесь Н — высота насадки; <£> — эквивалентный диаметр на- садки (е г— свободный объем насадки, м3!м3, а — удельная поверхность насадки, м2!м3)\ да = истинная скорость газа в слое насадки (да0 —
фиктивная скорость газа, т. е. скорость газа, отнесенная ко всему сечению колонны); А. — коэффициент сопротивления, учитывающий суммарные потери давления на трение и местные сопротивления.
Значения коэффициента сопротивления X для различных насадок определяются по эмпирическим уравнениям. Так, для беспорядочно засыпанных кольцевых насадок находят: при ламинарном режиме движения газа (Иег <4 40)
. (ХМ1)
при турбулентном режиме движения газа (1?егГ> 40)
'щ- (XI-41а]
*4
где Иег= ™ —’^г- — критерий Рейнольдса для газа (рг и рг — плотность и вязкость газа).
Мт
Гидравлическое сопротивление орошаемой насадкй Арор больше сопротивления сухой насадки. Это объясняется тем, что некоторое коли чество жидкости задерживается в насадке вследствие смачивания ее по верхности и скопления в узких криволинейных каналах, образуемых со прикасающимися насадочными телами. При этом уменьшаются свободно« сечение и свободный объем насадки и соответственно увеличивается истин ная скорость газа да, в результате чего возрастает гидравлическое сопро тивление насадки.
Точный расчет Арор затруднителен, так как при одних и тех же зна чениях скорости газа и плотности орошения эта величина будет зависеп от способа загрузки насадки, возможной неоднородности размеров на садочных тел и т. п.
462 Гл. XI. Абсорбция Для приближенного расчета Ар0р при работе колонны в режимах ниже точки инверсии можно использовать эмпирическое уравнение ЛРор = ЮШ Арсух (XI,42) где и — плотность орошения, м3/(м2-сек); Ь — опытный коэффициент, значения которого приводятся в специальной литературе * [например, для насадки 25Х 25Х 3 мм при и ==• ■= (0,5 — 36,5) 10~3 м3/(мг-сек) величина Ь= 51,2)]. Смоченная поверхность насадки. В абсорберах, работающих в режимах ниже точки инверсии, насадка может не полностью смачиваться жидкостью. В этом случае поверхность массопередачи будет меньше поверхности насадки. Отношение удельной смоченной поверхности асм (мг/м3) ко всей удельной поверхности насадки называется коэффициентом смачивания насадки и обозначается через о}з. Таким образом а Коэффициент смачивания в значительной мере зависит от величины и способа подачи орошения на насадку или от числа точек орошения пор. Величина гр возрастает с увеличением и и пор до определенных значений этих параметров, после чего остается практически постоянной. Коэффициент возрастает также с увеличением размеров насадочных тел. Коэффициент смачивания насадки может быть найден по уравнению <ф = 1 —Ае~т (Х1.43) плот- В этом уравнении показатель степени т — ОДе” = С (-‘^>жЛ, где рж и |ЛЖ аЦж / ность и вязкость жидкости. Значения постоянных А, С и п приводятся в литературе** (например, для колец размерами 15—35 мм: А — 1,02, С =■ 0,16, п= 0,4). С помощью уравнения (XI,43) по известной удельной поверхности а насадки можно рассчитать ее удельную смоченную поверхность асм. Скорость газа при работе абсорбера в режимах ниже точки подвисания не оказывает заметного влияния на величину о}з. Выше точки подвисания коэффициент смачивания возрастает с увеличением скорости газа. Однако не вся смоченная поверхность активна для массопередачи, так как в определенных точках слоя насадки (например, в точках контакта между насадочными телами) могут образовываться застойные зоны. Если обозначить через ва удельную активную поверхность насадки, то доля ее активной поверхности г|за = аа/а для правильно уложенных насадок (кольца, трубки и т. д.) может быть приближенно определена как функция от плотности орошения (при 1/ <3 <• 0,003 м3/(м2-сек): “ 0,0005 + 0,8
При [/> 0,003 м3/(мг-сек) для уложенных насадок смачивается практически вся их геометрическая поверхность и г|за 1. Для неупорядоченно засыпанных насадок удельную активную поверхность приближенно можно найти по формуле
85[/ .. .
0,00125 + и (XI ,44а)
Диаметр абсорбера £> определяют по общему уравнению (Х,75). Фиктивную скорость газа принимают с учетом рассчитанной зо уравнению (XI,25) предельной скорости ш>3. Полученное значение £> 1еобходимо скорректировать с учетом того, что насадка должна по воз
* См., например: X о б л е р Т. Массопередача и абсорбция. М., «Химия», 1964. Зм. с. 356.
** См. Р а м м В. М. Абсорбция газов. М.,; «Химия», 1966. См. с. 443.
6. Расчет абсорберов 463 можности полностью смачиваться жидкостью, а плотность орошения связана с О зависимостью [в м31{м2-сек): »-5» где Ьв — расход абсорбента, м31сек. Найденное значение I/ подставляют в уравнение (XI,43) и определяют величину коэффициента <ф. Если при данной плотности орошения и значение 1|? близко к единице, то рассчитанную величину £> можно считать удовлетворительной. Если же нужно улучшить смачивание насадки, т. е. увеличить гр, то необходимо либо повысить расход поглотителя (с последующим пересчетом И)о), либо заменить принятую насадку на насадку больших размеров. В последнем случае возрастает фиктивная скорость газа и соответственно уменьшится площадь поперечного сечения колонны. Высота а бсорбера. Рабочая высота абсорбера (высота слоя насадки) определяется на основе требуемого объема насадки Утс, который, в свою очередь, зависит для данной насадки от поверхности массопередачи Л которую находят по общему уравнению массопередачи (Х,46) или (Х,46а). Тогда объем насадки ^нас = (XI,46) где 5 — площадь поперечного сечения колонны, м2. Отсюда, подставляя вместо ее выражение из уравнения (Х,46), получают И = ^нзс _ В М — Ба^Ку &Уср При расчете коэффициента массопередачи Ку или Кх коэффициент массоотдачи рг в газовой фазе для неупорядоченно загруженных насадок может быть определен по уравнению Ми’ = 0,407Не°’655 (РГр)3 (XI,48) Соответственно высота единицы переноса для газовой фазы в этом случае _2 Аг = 0,615^е®-345 (Рг')3 (XI,49) Для регулярных насадок Ыи^ОЛб?!^'74^;;^ (X 1,50) ИЛИ / I \0.47 Аг= !,5^ре0'26(рг;)3 (XI,50а) где I — высота элемента насадки. ■ В уравнениях (XI,48)—-(XI,50) в критериях №1’ = |Зг^э/Ог и Ие = заопре- деляющии геометрическим размер принят эквивалентный диаметр насадки = 4е/а. При расчете высоты насадки через высоту единицы переноса для газовой фазы значения /гг, вычисленные по уравнениям (XI,49) и (XI,50а), следует делить на долю активной поверхности г|)а, найденную соответственно на основании формул (XI,44) и (X 1,44а). Для кольцевой насадки коэффициент массоотдачи в жидкой фазе, отнесенный к единице полной геометрической поверхности насадки, можно определить по уравнению 1^4 = 0,0021 КеУ3(р4)0-5 (XI,51; в котором критерий №1ж = (Зжбпр/Т>ж рассчитан по приведенной толщине пленки (см. стр. 116). Соответственно йж = П98прЯе^(Р4)<’.5 (Х1,51а
464 Гл. XI. Абсорбция Для более точного расчета Н движущую силу процесса в уравнении (Х,46) или (Х,46а) следует определять с учетом обратного перемешивания газа в абсорбере (см. стр. 419). Расчет тарельчатых абсорберов. Величину гидравлического сопротивления Дрт барботажных тарелок рассчитывают как сумму трех частных сопротивлений: ДРч = Дрс. т “Г ДРа + ДРгж (X1,52) где А/?с. т — сопротивление сухой тарелки; Дра — сопротивление, обусловленное силами поверхностного натяжения жидкости; Дргж— сопротивление газо-жидкостного слоя на тарелке. Величина Арс.х определяется по уравнению в котором w0!B = w/Fсв — скорость газа в отверстиях тарелки. Коэффициент сопротивления тарелок изменяется в широких пределах (0,5—4) и зависит от конструкций тарелок. Значения £ приводятся в специальной литературе *. Потере давления Ар0 на преодоление сил поверхностного натяжения жидкости о при входе в слой жидкости на тарелке; = (XI,53) “э
Для тарелок, работающих в струйном режиме, значением Ар0 можно пренебречь.
Гидравлическое сопротивление газо-жидкостного слоя на тарелке Аргж принимают равным статическому давлению слоя:
^Ргж = ^оРж? = ^гжРгж? (X 1,54)
где h0 и /ггж — высота жидкости я газо-жид;:остного слоя на тарелке; рж и ргж — плотность жидкости и газо-жидкостнои смеси на тарелке.
Значения Аргж могут быть определены по эмпирическим уравнениям **.
Для провальных, ситчатых и клапанных тарелок величину /ггж можно рассчитать по уравнению:
Ещ |/77в = 0,25Fr0~1,25 (XI,55)
Рж
где Ей о = Д/?гж/рга'“тв — критерий Эйлера; Fr0 = x>2OTB/ghrlli — критерий Фруда.
Брызгоунос. Как отмечалось в главе II, при определенных скоростях газ начинает увлекать с собой капли жидкости, которые образуются при разрыве пузырьков, выходящих на поверхность барботажного слоя; при этом капли попадают с потоком газа на вышерасположенную тарелку. Унос жидкости газовым потоком приводит к снижению движущей силы процесса массопередачи, увеличению жидкостной нагрузки сливных устройств, потере абсорбента с уходящим из абсорбера газом и является одной из основных причин, ограничивающих возможность интенсификации тарельчатых аппаратов.
Величина уноса е ориентировочно не должна превышать 5—10% общего количества подаваемой в абсорбер жидкости. Величина уноса возрастает с увеличением скорости w газа в колонне и уменьшением высоты Нсп сепарационного пространства, причем Нсп — Я* — кгж, где #т — расстояние между тарелками.
* См. Справочник химика, Т. V., М., «Химия», 1966. См. с. 375.
** См., например: Стабников В. Н. Ректификационные аппараты. М., «Машиностроение», 1965- 353 с.
6. Расчет абсорберов 465 Для ситчатых тарелок, например, величина уноса жидкости (в кг/кг газа) может быть найдена по уравнению где а — поверхностное натяжение жидкости, Мн/м. Для предотвращения потерь абсорбента вследствие уноса над верхней тарелкой часто устанавливают сепарационное устройство (например, слой насадки или металлической сетки и т. п.). Аналогичные устройства имеются и в колоннах других конструкций. Поверхность контакта фаз. Обычно поверхность контакта фаз определяют как поверхность находящихся в барботажном слое пузырьков. В этом случае удельную поверхность контакта фаз а (в м*/м3) можно вычислить из уравнения (11,141): 6в где в — газонаполнение барботажного (пенного) слоя, м?/м'\ йср — средний поверхностнообъемный диаметр пузыря, м (см. стр. 112). Поверхность % контакта, отнесенная к единице площади тарелки, выражается как В барботажных абсорберах процесс осуществляют обычно при скоростях газа, значительно превышающих скорость свободного всплывания пузырька. При этом поверхностью контакта фаз является как поверхность газовых струй, которые проходят через барботажный слой, не разбиваясь в нем на отдельные пузырьки, так и поверхность капель, образующихся над этим слоем при разрушении пузырьков. Определение поверхности струй и капель затруднительно. Кроме того, при измерении диаметра пузырька возникают трудности, связанные с усреднением замеренных значений и получением достаточно точных величин £?Ср- Некоторые данные о поверхности контакта фаз приводятся в специальной литературе *. Диаметр абсорберов. Диаметр О аппарата определяют по принятой фиктивной скорости газа в колонне, пользуясь общим уравнением (Х,75). Высота абсорберов. Рабочую высоту И (расстояние между крайними тарелками) барботажного абсорбера находят методами, указанными в главе X. При расчете Я по уравнению массопередачи коэффициент массопередачи определяется с помощью уравнения (Х.47) или (Х,48). Так как расчет поверхности контакта фаз на тарелке затруднителен, при обработке опытных данных по массопередаче в тарельчатых аппаратах коэффициенты массоотдачи относят чаще всего к сечению 5Т тарелки (точно определяемая величина), либо к объему пены = /гГЖЗт или жидкости на тарелке У0 = /г05х (где Нтж и Н.0 — высота пены и слоя жидкости на тарелке). Соответственно число единиц переноса на тарелку (пг или пж) выражается следующим образом: для газовой фазы (XI,56) (XI,58) для жидкой фазы (XI. 59) л * См. Ра мм В. М. Абсорбция газов. М., «Химия», 1966. См. с. 559.
466 Гл. XI. Абсорбция Коэффициенты массоотдачи рг,5т и рж, <;т, отнесенные к единице рабочей площади 5Т тарелки, связаны с поверхностными коэффициентами массоотдачи рг и рж соотношениями, аналогичными выражениям (Х,79): для газовой фазы где и рж у—коэффициенты массоотдачи в газовой и жидкой фазах, отнесенные к объему жидкости на тарелке. Коэффициенты массоотдачи или числа единиц переноса на тарелку обычно рассчитывают по уравнениям, применимым только для конкретной конструкции тарелки. Для колпачковых тарелок, например, число единиц переноса в газовой фазе пт можно определить по уравнению: где Ргг = vг/Dr — диффузионный критерий Прандтля для газа; vr — кинематическая вязкость газа, м^/сек; £>г — коэффициент молекулярной диффузии в газе, ж2/се/с; /гпер — высота переливной перегородки, м; д Соответствующее уравнение для расчета числа единиц переноса в жидкой фазе пж имеет следующий вид: Здесь £>ж — коэффициент диффузии в жидкой фазе, мЧсек\ тж — средняя продолжительность контакта фаз (в сек), которая находится по уравнению: где 2Т — длина пути жидкости, т. е. расстояние между сливными перегородками, м\ <?ср — линейная плотность орошения, отнесенная к средней ширине тарелки, м3/{м-сек). При применении ситчатых (переточных) и провальных тарелок имеем: для газовой фазы диффузионный критерий Пекле для жидкой фазы; /ггж — высота гаЗо-жидкостной смеси на тарелке, м. При расчете необходимого числа пд действительных тарелок (по любому из указанных выше способов) в первом приближении можно принять, что на всех тарелках со сливными устройствами фазы движутся по взаимно перпендикулярным направлениям и в этом случае движущую силу процесса вычисляют по схеме абсорбции с перекрестным током. На тарелках без сливных устройств движущую силу процесса рассчитывают по схеме полного перемешивания фаз (см. главу X, стр. 428). Рассчитав пя, определяют высоту Я абсорбера (в м) где #с. п — высота сепарационного пространства, м; кв — расстояние от верхней тарелки до крышки абсорбера (в м), которое принимается на основании конструктивных соображений. Рг. — Рг, УЛгж — Рг, 1'А0 — РгаЛгж (XI,60) для жидкой фазы Рж, — Рж, V' ^гж Рж. 1/^0 РжаЛ, (XI,61) и, Г (Рг')°-5 = 0,776 + 4,63/гпер — 0.238а.']/" рр + 0,0292? (X 1,62) »* = 30500^ (б8Апер+ 1)Тя (XI,63) (XI,64) (XI ,65) для жидкой фазы (XI, 66) (XI,67) — расход жидкости, отнесенный к периметру сливиой перегородки, м3/(м-сек). .
8. Схемы абсорбционных установок 467