Интернет-журнал «Науковедение» ISSN 2223-5167 http ://naukovedenie.ru/
Том 8, №2 (2016) http ://naukovedenie. ru/index.php?p=vol8-2
URL статьи: http://naukovedenie.ru/PDF/04TVN216.pdf
DOI: 10.15862/04TVN216 (http://dx.doi.org/10.15862/04TVN216)
Статья опубликована 15.03.2016.
Ссылка для цитирования этой статьи:
Биргер Е.В., Лебедев В.В., Зизюк А.П., Иванова Е.В. Проектирование абсорберов для очистки воздуха от газо -и парообразных примесей // Интернет-журнал «НАУКОВЕДЕНИЕ» Том 8, №2 (2016) http://naukovedenie.ru/PDF/04TVN216.pdf (доступ свободный). Загл. с экрана. Яз. рус., англ. DOI: 10.15862/04TVN216
УДК 628.5
Биргер Евгений Владимирович
НОУ ВПО «Институт мировой экономики и информатизации», Россия, Москва
Декан факультета «Информатизации и экологии» Кандидат экономических наук E-mail: [email protected]
Лебедев Владимир Владимирович
ФГБОУ ВПО «Российский государственный университет туризма и сервиса», Россия, Москва
Доцент кафедры «Сервисного инжиниринга» Кандидат технических наук E-mail: [email protected]
Зизюк Алексей Петрович
ФГБОУ ВПО «Российский государственный университет туризма и сервиса», Россия, Москва
Магистр
E-mail: [email protected]
Иванова Елена Владимировна
ФБУ «Научно-исследовательский центр по проблемам управления ресурсосбережением и отходами», Россия, Мытищи
Научный сотрудник E-mail: [email protected]
Проектирование абсорберов для очистки воздуха от газо- и парообразных примесей
Аннотация. Разделительные массообменные процессы широко применяются в промышленности, как правило, в аппаратах колонного типа для очистки воздуха и воздушных смесей. Наибольшее применение они находят в технических устройствах очистки воздуха в жилых, общественных и производственных помещениях. Используются подобные системы при абсорбционной очистке воздуха от газообразных примесей на предприятиях сервиса. Для проектирования таких установок на современном уровне нужны эффективные методики, доступные для применения техническими специалистами в процессе разработки новых установок и совершенствовании находящихся в эксплуатации. В настоящей работе представлена методика расчета абсорбционных аппаратов для очистки воздуха от газообразных и парообразных загрязнителей. Работа может помочь специалистам и обучающимся освоить основные принципы работы аппаратов абсорбционной очистки, методы расчета процессов и конструктивных параметров аппаратов.
Ключевые слова: абсорбция; граница раздела фаз; массообменные процессы; разделение газовых смесей; коэффициент массопередачи; площадь поверхности раздела фаз; рабочая линия процесса; число единиц переноса
Введение
Разделительные массообменные процессы широко применяются в промышленности и, как правило, осуществляют в аппаратах колонного типа. Для разделения газовых смесей широко применяют процесс абсорбции [1].
Абсорбцией называется процесс избирательного поглощения газов или паров из газовой фазы жидкими поглотителями. Жидкий поглотитель называют абсорбентом. Поглощаемый компонент газовой смеси называется абсорбтивом. Поглощенный компонент газовой смеси называют абсорбатом [2].
В процессах абсорбции газовая фаза (газовая смесь, парогазовая смесь) взаимодействует с жидкой фазой (поглотителем) путем непосредственного контакта. Через поверхность раздела фаз осуществляется перенос вещества из газовой фазы в жидкую. Как правило, граница раздела фаз представляет собой либо зеркало жидкости (поверхностные абсорберы) или поверхность плёнки жидкости, стекающей по поверхности (плёночные абсорберы) или насадке (насадочные абсорберы), либо поверхность проходящих через жидкость пузырьков газа (барботажные абсорберы) (например, тарельчатые и др.), либо поверхность капель жидкости (распыливающие абсорберы) [3]. На режимы образования межфазной поверхности может влиять взаимное направление и скорость относительного движения фаз. Например, в противоточных насадочных абсорберах колонного типа, в которых газовая смесь движется вверх навстречу стекающей по насадке под действием силы тяжести жидкости, при определённых скоростях газовой смеси возникает режим подвисания жидкости, при котором режим формирования границы раздела фаз из плёночного переходит в пузырьковый (барботажный) [4].
В процессе абсорбции происходит выделение тепла, которое связано с энергией растворения газов или конденсации паров и растворения конденсата в жидкости. Это тепло характеризуют удельной теплотой растворения или теплотой абсорбции.
Движущей силой процесса абсорбции является положительная разность концентраций целевого компонента в газовой смеси и равновесной концентрацией этого компонента в газовой смеси над жидким поглотителем в рабочих условиях [5].
Увеличение давления компонента над абсорбентом увеличивает движущую силу и повышает скорость поглощения газообразного компонента в жидкости. Повышает также скорость абсорбции охлаждение абсорбента, поскольку при повышении температуры жидкой фазы за счет теплоты абсорбции растет равновесная концентрация, что снижает движущую силу процесса.
Десорбция - обратный абсорбции процесс - это выделение поглощённого газа или отгонки поглощённого пара из поглотителя. Факторы, определяющие скорость десорбции противоположны по действию факторам, способствующим абсорбции. Это - нагревание абсорбента, понижение концентрации целевого компонента над ним ниже равновесного и добавление в абсорбент растворяемых вытесняющих веществ, которые снижают растворимость целевого компонента в абсорбенте, повышая его равновесное давление над жидкостью.
Материальный баланс, линия равновесия и рабочая линия
В процессе абсорбции взаимодействуют две фазы - газовая и жидкая. Газовая фаза является отдающей, а жидкая - принимающей. В процессе участвует, как минимум, один целевой компонент (распределяемое вещество), который переходит через поверхность раздела фаз, растворяясь в жидкости.
Содержание целевого компонента в фазах выражают при расчётах в относительных массовых единицах, относя массовое количество (моль или кг) вещества целевого компонента в фазах к массовому количеству инертного вещества фазы. Инертное вещество в газовой фазе - это газовая смесь за вычетом целевого компонента - носитель, а в жидкой фазе - чистый абсорбент. Количество газа - носителя в газовой смеси обозначают О (кг/сек), а количество чистого абсорбента обозначают Ь (кг/ск). Также обозначают для краткости и фазы. Относительное массовое содержание целевого компонента в фазах О и Ь определяют по
Масса целевого компонета „ ,
отношению к инертному компоненту: -. В газовой фазе
Масса инертного компонента
относительное массовое содержание обозначают У (кг/кг), а относительное содержание целевого компонента в жидкой фазе - X (кг/кг). Относительные массовые содержания могут быть иначе выражены отношением количества молей - моль/моль, что удобно, если количество взаимодействующих инертных компонентов фаз выражено в моль/с [6]. Массовые или мольные доли отличает от относительного массового содержания то, что они отнесены к полной массе или полному количеству молей рассматриваемой фазы, поэтому для них справедливы следующие нормирующие уравнения:
I х = 1
]
I У] = 1'
]
где xj - массовая или мольная доля целевого компонента в жидкой фазе (абсорбенте);
у^ - массовая или мольная доля целевого компонента в газовой фазе.
Исходное содержание вещества в фазах О и Ь может также быть задано в объёмных или массовых долях, или в единицах массовой концентрации вещества (например, кг/м3) [6]. В этом случае следует перевести эти величины в единицы относительного массового содержания.
В процессе абсорбции относительное массовое содержание распределяемого вещества в газовой фазе уменьшается от начального Уу до конечного У2, при этом относительное массовое содержание распределяемого вещества в жидкой фазе увеличивается от начального Х1 до конечного Х2. Закон сохранения массы в процессе абсорбции можно выразить уравнением материального баланса:
о -(у - у2 )=ь -x - x) (1)
Уравнение (1) выражает закон сохранения массы, т.е. равенство в процессе абсорбции количества распределяемого вещества, теряемого газовой фазой, количеству вещества, поглощаемого жидкой фазой. Это уравнение справедливо для любой стадии процесса (как непрерывного, так и периодического).
В состоянии фазового равновесия между количеством распределяемого вещества в газовой и жидкой фазах существует функциональная однозначная зависимость:
У * = /(X ), (2)
где Y* - равновесное относительное массовое содержание распределяемого вещества в газовой фазе, соответствующее в условиях фазового равновесия его содержанию в жидкой фазе Х.
Характер зависимости определяется природой взаимодействующих веществ, их количеством, давлением и температурой в системе. По правилу Гиббса при фазовом равновесии в двухфазных однокомпонентных системах имеется одна степень свободы: С=К-Ф+2=1-2+2=1. Если давление и температура системы постоянны, то при фазовом равновесии может изменяться количество вещества в фазах [7]. Линия равновесия (2) выражает связь при фазовом равновесии при постоянном давлении и температуре между количеством распределяемого вещества в фазах G и Ь.
В практических расчётах предпочтительно пользоваться опытными данными, которые получают в изучаемой системе при различных парциальных давлениях распределяемого компонента при фазовом равновесии. Опытные данные могут быть выражены табличной зависимостью между значениями равновесного парциального давления р* в газе и концентрации растворённого вещества с в жидкости.
Равновесное относительное массовое содержание Y* можно вычислить по формуле:
т * Мг р
У =—к---г, (3)
Мин Р - Р
где Мк - молекулярная масса абсорбируемого компонента;
Мин - молекулярная масса инертного компонента (смеси компонентов) газовой фазы (носителя);
Р - общее давление газовой смеси.
Соответствующее относительное массовое содержание абсорбированного компонента в жидкости можно определить по формуле:
X = —, (4)
Р-с
где с - концентрация поглощённого компонента в жидкости, кг/м3;
р - плотность жидкой фазы в условиях равновесия, кг/м3.
По полученным табличным значениям можно построить линию равновесия на диаграмме в координатах Y-X.
В частных случаях возможно построение линий равновесия с использованием отдельных законов фазового равновесия газов, паров и жидкости.
Рассмотрим метод построения линии равновесия для трех видов систем, заданных в вариантах расчёта, аналитическим способом.
1. Если распределяемое (поглощаемое) вещество - газ, который химически не взаимодействует с жидкостью, а только способен незначительно в ней растворяться, то при небольших концентрациях его в газовой фазе при фазовом равновесии справедлив закон Генри, согласно которому равновесное парциальное давление газа над жидкостью пропорционально мольной доле газа, растворённого в жидкости:
Р = E • х, (5)
где х -мольная доля (моль/моль смеси) газа в жидкости, кг/кг смеси ( х < 1); Е - коэффициент пропорциональности (коэффициент Генри).
Растворитель при фазовом равновесии также находится в равновесии со своим паром. При небольших давлениях и концентрациях равновесие жидкости и паров над ней подчиняется закону Рауля, согласно которому равновесное парциальное давление пара над жидким раствором равна произведению упругости насыщенных паров над чистым растворителем при рабочей температуре на мольную долю чистого растворителя в растворе:
Pl = Pls '(1 - х) > (6)
где Pls - упругость насыщенных паров над чистой жидкостью (не содержащей растворённых веществ).
В парогазовой смеси над поверхностью жидкости действует также закон Дальтона, согласно которому давление газовой смеси равно сумме парциальных давлений всех компонентов. В применении к данному частному случаю можно в соответствии с законом Дальтона записать уравнение:
П = Р + Pl + Pn , (7)
где П - общее давление в системе;
р - парциальное давление распределяемого компонента в газовой смеси; Pl - парциальное давление паров поглотителя в газовой смеси;
pin - парциальное давление инертной части газовой смеси (например, газа - носителя).
Для газовой смеси, которую можно считать смесью идеальных газов (при небольших давлениях), мольные доли компонентов равны объёмным долям и определяются отношением парциального давления компонента к общему давлению смеси. Из уравнения (7) следует, что сумма мольных долей газовой смеси равна 1:
y + yL + Уп =1 (8)
p E • х
где y = — = = m • х - мольная доля распределяемого компонента (с учётом (5) в
E
состоянии равновесия, причём m = — );
PL PLS • (l - х) (л \
yL = =-—-= y¿s ' (1 — х) - мольная доля паров жидкости (с учётом (6) в
PLS
состоянии равновесия, причём yLS = );
Р
yin = 1 — у — yL - мольная доля инертного газа (с учётом (8)).
Относительную массовую долю распределяемого компонента в газовой фазе можно определить по формуле (см. (3)):
Г = Мж_ .y = Мж___y__(9)
Mn ym Mn 1 - y - Уь
где Мк - молекулярный вес компонента, кг/кмоль;
Min - молекулярный вес инертного газа, кг/кмоль.
Относительную массовую долю распределяемого компонента в жидкой фазе можно определить по формуле:
т. Mr x
X = TT "л-, (Ш)
ML 1 - x
где Ml - молекулярный вес жидкого поглотителя.
Можно изменять значения x, вычислять соответствующие значения y, а затем по (9) и (10) Y и X. Результаты можно свести в таблицу - это будет представлять табличную форму выражения линии равновесия.
Можно получить аналитический вид функциональной зависимости линии равновесия:
т * Mr m ■ X
Y = MK ^М--(11)
Mn MK ■(! - y )+(1 - m)-X
2. Если распределяемое вещество пар, а конденсат пара образует с жидким абсорбентом идеальную смесь жидкостей. Это означает полную взаимную растворимость конденсата и абсорбента, а давление паров обеих жидкостей над их смесью в состоянии равновесия подчиняется закону Рауля:
*
Р = x - Ps (12)
PL =(1 - x)- Pls , (13)
где Ps - упругость насыщенных паров компонента в равновесии с чистой жидкой фазой компонента.
Согласно закону Дальтона в состоянии равновесия в газовой смеси над жидкой смесью поглотителя и распределяемого вещества:
П = Pin + Р* + pL = Pin + x ■ Ps + (1 - x)- Pls (14)
Согласно определению мольных долей компонентов в газовой смеси (см. 8) их сумма будет равна 1, а определение мольных долей паров жидкого поглотителя и поглощаемого компонента с учётом (12) и (13) будет следующим:
x - Ря
y = = x - ys, (15)
v Ps где ys~
П
Уь = (1 x)n PLS =(1 - x)-yLS , (16)
PLS
где yLS
LS П
Равновесные относительные массовые доли распределяемого компонента в газовой и жидкой фазах можно рассчитать по формулам (9) и (10). Так же, как в п.1, можно вычислить табличные значения Y-X линии равновесия в исследуемом интервале значений х. Аналитическая зависимость для определения функции равновесия:
* Mr Уо • X Y =—— —-—--(17)
Mn MtL ^ — yLS )+(1 — ys )• X
3. Система аналогична п.2, но жидкая смесь неидеальна. В силу особенностей взаимодействия компонентов в смеси равновесное давление паров над жидкостью будет
отличаться от того, что предсказывает закон Рауля:
*
Р =У х • Ps (18)
PL =Yl -(1 — х )• Pls (19)
где y, yL - активности поглощаемого компонента и поглотителя.
Активности отражают отклонение от закона Рауля, учитывая эффект взаимодействия компонентов в смеси жидкостей. Активности зависят содержания компонентов в смеси.
Равновесные мольные доли паров в газовой фазе определяют аналогично (15) и (16) с учётом (18) и (19):
У • х • Do
У = FS = У • х • ys (20)
yL = yL '(1 PLS = УL -(1 — х)^ yLS (21)
Аналогично пп. 2 и 3 можно рассчитать равновесные относительные массовые доли распределяемого компонента в газовой и жидкой фазе и построить табличную зависимость Y-X для линии равновесия. Также можно воспользоваться аналитической формой зависимости для линии равновесия:
Y * = M—__УyS • X__(22)
Mn MtL -(1 — yL • yLS ) + (1 — У ys )• X
Для вычислений надо знать функции зависимости активностей компонентов от мольных долей компонентов в жидкой бинарной смеси.
Формулы для определения активностей компонентов бинарной смеси жидкостей приводятся в [1, с. 293]:
A__A
х A
lg У =-7 =-7 (23)
f а\2 Г л\2
1 +
V XL B
х A
1 + 7-Г —
. (1 — х) B
^ Уь
В
В
г
X1
„.г В 1 + — ■ —
х А у
\2
1 ■ В
А у
(24)
V
х
где А, В - коэффициенты, определяемые для каждой бинарной смеси по опытным данным;
х, xь - мольные доли в жидкости поглощённого компонента и чистого поглотителя.
Для таких смесей в пособии заданы равновесные мольные доли компонентов в табличной форме. Для нескольких столбцов (3-4) таблицы выписываем пары значений х-_у в серединах строк, так как на концах строк при х=0 и х=1, что соответствует чистым компонентам жидкости, активности равны у=уь=1. Затем вычисляем значения активностей компонентов по формулам:
У
У
х ■ У*
Уь —
Уь
Уь
хь ■ уи
(1 - х )■ У и
(25)
(26)
В формулах (23), (24), (26) мольная доля в жидкости чистого поглотителя для бинарной смеси выражается через мольную долю поглощённого компонента:
X^ — 1 X
(28)
Для полученных по формулам (25), (26) нескольких значений активностей вычисляем соответствующие значения коэффициентов А и В:
А — ^ у
г 1 л2 1 + Ъ ■ ^ уь
V
X ^ у
— ^ У'
У
1 +
1 - х . МУ.
У
х
В — lg Уь
Г 1
1 ■ ъУ
V
X
ь
^ Уь
^ Уь
1 +
X
}%У_ 1 - x lg Уь
(29)
(30)
Полученные значения коэффициентов А и В усредняют:
Е А
А —
В —
т
ЕВ
т
где т— количество расчётных точек.
Подставив полученные средние значения коэффициентов в формулы (23) и (24), можно, используя операцию потенцирования, определить окончательно функции активностей от мольной доли поглощённого компонента:
2
2
2
А
(31)
у = 1 если х = 0 или х = 1
в
у = 10 если х ф 0 или х ф 1
у = 1 если х = 0 или х = 1
(32)
Теперь можно приступать к расчёту линии равновесия (22).
Построение рабочей линии выполняем на основании уравнения материального баланса с учетом относительного движения фаз в аппарате. Возможны две схемы относительного движения газовой и жидкой фаз: прямоточная и противоточная, см. рис. 1.
Жидкий поглотитель, как правило, поступает в верхнюю часть колоны абсорбера и движется вниз под действием силы тяжести, стекая по насадке, с тарелки на тарелку или падая в виде капель. Газ может поступать в аппарат в том же сечении и двигаться в том же направлении - сверху вниз, что и жидкость. Такой аппарат работает по прямоточной схеме. Газ и жидкость выходят из аппарата в его нижнем сечении. Если газ поступает в нижнюю часть колоны, из которой выводится из аппарата отработанный жидкий поглотитель, и далее движется вверх навстречу потоку жидкости, то такой аппарат работает по противоточной схеме. Газ отводится из верхней части колоны, куда поступает жидкость.
На рис. 1 для прямоточной и противоточной схем обозначены на схемах аппарата направления потоков взаимодействующих фаз, начальные и конечные относительные массовые содержания распределяемого компонента в фазах, а также на графиках показано расположение рабочей линии и линии равновесия в режиме абсорбции.
Рабочая линия должна проходить выше линии равновесия при абсорбции: при этом движущая сила процесса массопередачи - разность между рабочим относительным массосодержанием и равновесным относительным массосодержанием поглощаемого компонента в газовой смеси будет положительной.
ПРямоток
С, У?
С, У:
"У .• У2 У!
П Р и т м в о т о ь
а)
Б "тУ
б)
Рисунок 1. Прямоточная а) и противоточная б) схемы относительного движения фаз.
1 - Рабочая линия; 2 - Линия равновесия.
Уравнение рабочей линии можно получить, если задать входное сечение, например по газу, и записать уравнение материального баланса распределяемого вещества в фазах между выбранным входным сечением и другим произвольным сечением аппарата, которое расположено между входным и выходным сечением, например по газу, где относительные массовые содержания поглощаемого компонента соответствует точке У-Х рабочей линии в данном сечении.
Для прямоточной схемы материальный баланс между входным по газу и произвольным сечением, отстоящим от входного по направлению движения газа, будет выглядеть так:
С-(71 - 7) = L■ (X -X1),
откуда получаем уравнение рабочей линии:
7 = 71 - (X - X )= 71 -/-(X - X ), (33)
где / = — - удельный массовый расход жидкого чистого поглотителя, равный
С
отношению массовых расходов чистого поглотителя и инертного газа-носителя.
Для противоточной схемы материальный баланс между входным по газу и произвольным сечением, отстоящим от входного по направлению движения газа, будет выглядеть так:
С - (71 - 7) = L - (X2 - X), откуда получаем уравнение рабочей линии:
7 = 71 +/- (X - X2 ) (34)
Рабочие линии - прямые линии. При прямотоке наклон рабочей линии к оси абсцисс (ось Х) будет отрицательным, а при противотоке - положительным.
Поскольку расход газа и начальная и конечная концентрации компонента в газе заданы, а в жидкой фазе задана начальная концентрация компонента, то в расчёте сначала надо определить расход жидкой фазы. На графиках, рис. 1, показано положение рабочих линий, при котором удельный расход жидкости будет минимальным. Это соответствует случаю, когда рабочая линия в выходном, или ином сечении коснётся равновесной, при этом массосодержание компонента в этом сечении будет равно равновесному.
Для прямотока это соответствует точке линии равновесия:
72 = / (X 2)
Для противотока это соответствует точке линии равновесия:
71 = / (X 2)
Вид функциональной зависимости / (X*) для линии равновесия может
соответствовать одной из формул (11), (17), (22). Координата искомой точки Х* может быть определена по 71 или 72 из формул, либо по графику линии равновесия.
Минимальный удельный расход чистого жидкого абсорбента определим по формуле:
/ _ 71 - 72
(35)
Л 2 - Л 1
Обычно, если не заданы другие условия, принимают начальную относительную массовую долю поглощаемого компонента в абсорбенте равной нулю: Х1=0.
Удельный расход абсорбента обычно принимают выше минимального I > /т;п. Обычно его значение определяют в процессе технико-экономической оптимизации, так как
при изменении удельного расхода от до да имеют место конкурирующие влияния данного фактора на технико-экономические показатели установки. При увеличении расхода с одной стороны повышается средняя движущая сила абсорбции, что приводит к уменьшению габаритов абсорбера, а с другой стороны растут затраты на перекачку абсорбента, связанные с ростом гидравлических потерь. При снижении расхода абсорбента до средняя движущая сила уменьшается, что ведет в пределе к увеличению необходимой поверхности массопередачи до да, при этом требуется неограниченное увеличение габаритов абсорбера.
Массовый расход абсорбента будет равен ь — l ■ G. Назначив расход, определим выходное относительное массосодержание абсорбата в абсорбенте:
У - У
X 2 — х + ^ (36)
Теперь можно окончательно построить рабочую линию процесса абсорбции.
Полученные аналитические выражения также могут быть применены при расчёте нестационарных периодических процессов абсорбции.
Расчёт размеров абсорбционных аппаратов
Расчет абсорбционной колонны насадочного типа.
Необходимо в начале расчета задать исходные данные. Во-первых, к ним относятся, полученные выше функции линий равновесия фаз в системе и рабочей линии, расходы газовой и жидкой фазы. Во-вторых, необходимо определить физические свойства газа и жидкости при рабочих параметрах (температуре и давлении): плотности, коэффициенты вязкости и др. В-третьих, необходимо выбрать тип насадки и задать для нее параметры размеров элементов, удельной площади поверхности, пористости слоя насадки, коэффициент извилистости.
Режимы работы колонны зависят от скоростей движения газовой смеси и жидкости (абсорбента). Газ обычно движется через насадку снизу вверх навстречу стекающей по насадке жидкости. Для характеристики движения сред в колонне рассчитывают фиктивные средние скорости движения:
G
—-— (37)
Рг ■
ь (38)
Ж г, П
рж ■
где рг, Рж - плотности газа и жидкости, кг/м3; *к - площадь поперечного сечения колонны, м2.
Главным геометрическим параметром, влияющим на гидравлические характеристики насадок, является эквивалентный диаметр:
4 ■ е ■ щ
йэ —-^ , (40)
а
где е - пористость слоя насадки; у - коэффициент извилистости; а - удельная площадь поверхности насадки, м2/м3.
Удельная площадь поверхности насадки изменяется в широких пределах: для деревянных хордовых - от 50 до 100 м2/м3; для других - от 60 до 450 м2/м3. Пористость или доля свободного объема насадок меняется: для хордовых - от 0,55 до 0,77; для беспорядочно уложенных керамических колец - от 0,7 до 0,8; для правильно уложенных колец - от 0,7 до 0,74.
Обычно скорость газа принимают, учитывая то, что она не должна превышать предел захлёбывания. Режим подвисания - это режим торможения поднимающимся газом стекающей по насадке вниз жидкости при повышении скорости газа до начала подвисания. При дальнейшем повышении скорости до предела подвисания начинается режим захлёбывания. При скоростях газа, превышающих предел захлёбывания, течение жидкости обращается вспять и начинается унос ее газом. Работа насадочной колонны в режиме начала подвисания считается оптимальной по образованию максимальной поверхности межфазного контакта и эффективности массоотдачи. Оценить скорость газа в колонне, соответствующую началу режима подвисания, можно, например, по уравнению, рекомендованному в [1, с. 335]:
СсЛ0,43
Яе = 0,045- Аг0 57 -I , (39)
где Яе - критерий Рейнольдса;
Аг - критерий Архимеда.
Определяем критерий Явг для газа:
Яе = ——, (40)
/лг
где лг - коэффициент динамической вязкости газа, Па-сек. Определим критерий Архимеда:
Л 3 2
Аг = 8• ^ \рг Р-Рг , (41)
Иг рг
где §=9,81 м/с2 - ускорение свободного падения.
По формуле (39) сначала определяем Яе, а затем скорость газа, соответствующую началу режима подвисания:
^ = Яе —(42)
Рг Э
Задаем значение скорости равное или несколько меньшее величины, найденной по ф-ле (42).
Определяем, задавшись значением скорости газа, площадь поперечного сечения колонны:
Б к =--(43)
Рг
Затем скорость жидкости, см. ф-лу (38). Определяем критерий Явг для газа:
Яе г —-, (44)
а 'Уг
где \г - коэффициент кинематической вязкости, м2/сек.
Определяем критерий Явж для жидкости:
4 ■ w ■ ш
Яе ж —-—, (45)
а ■ уж
где Уж - коэффициент кинематической вязкости, м2/сек.
Для определения коэффициента массопередачи Ку, кг/(м2-с), который отнесен к разности (у - У *) относительного массодержания абсорбтива в газовой фазе, где У * -равновесное содержание соответствующее содержанию абсорбата X в жидкой фазе ( У * — / (х), т.е. на линии равновесия), по формуле:
КУ—т^т (46)
— + —
Ру Рх
необходимо определить коэффициенты массоотдачи /5у - со стороны газа, кг/(м2-с), вх -со стороны жидкости, кг/(м2-с), а также коэффициент т.
Коэффициенты массоотдачи определяют по критериальным уравнениям, полученным на основе экспериментальных данных в виде уравнений регрессии.
Коэффициент массоотдачи в газовой фазе для насадочных колонн с орошаемой насадкой определяем по уравнению [1, с. 345]:
Ыиог — 0,407 ■ Яе/6" Ргог^
V йЭ у
(47)
рг ■ йэ
где Ыиог —--диффузионный критерий Нуссельта;
Вг
РГдГ — - диффузионный критерий Прандтля;
^ Т
V
V
вг - коэффициент массоотдачи в газовой фазе, м/с;
Вг - коэффициент диффузии абсорбтива в газовой фазе, м2/с;
Н - высота рабочей части колонны (насадки).
'я
У И-0'47
При первом вычислении можно положить в (47)
V йэ у
— 1.
После определения диффузионного критерия Нуссельта по ф-ле (47), вычисляем коэффициент массоотдачи:
Рг — Мы0г ■ ^ (48)
й Э
Коэффициент массоотдачи в жидкой фазе для насадочных колонн определяем по уравнению [1, с. 345]:
Ми Вж = 0,00216-Яе ж °'77-РгДж ^
(49)
Где МиПж =
Рж
(
Б
N V
2 Л/3
диффузионный критерий Нуссельта;
' ж V Рж ■ § у
вж - коэффициент массоотдачи в жидкой фазе, м/с;
Лж - коэффициент динамической вязкости жидкости, Па-сек.
В качестве определяющего размера в (49) принимается приведенная толщина пленки
жидкости 8 =
( 2 V3
г* ж
2
V рж ■ § У
м.
После вычисления диффузионного критерия Нуссельта по (49) определяем коэффициент массоотдачи:
Рж = Ми£
Б
8
(50)
Значения коэффициентов диффузии определяют по справочным данным или расчётным путём.
Для диффузии газа А в газе В (или наоборот) коэффициент диффузии, например, вычисляют по формуле [3, с. 288]:
Б = 0,042183 -
3/ Т /2
Р-{»А + ^
11
+
МА мв
где Бг - коэффициент диффузии, м2/с; Т - температура, К; р - давление (абсолютное), Па; Ма, Мв - мольные массы газов А и В; »а, ив - мольные объёмы газов.
Для диффузии растворенного компонента (А) в жидкости (растворителе, В) коэффициент диффузии можно приближенно вычислить по формуле [3, с. 289]:
Б™ = -
10"
1 1
+-
МА Мв
А - В + )2 ^
где Бж - коэффициент диффузии, м2/с; Л - динамический коэффициент вязкости, мПа-с;
Ма, Мв - мольные массы растворенного вещества (А) и растворителя (В); »а, ив - мольные объёмы растворенного вещества и растворителя;
2
6
А, В - коэффициенты, зависящие от свойств растворенного вещества и растворителя. В разбавленных растворах коэффициент диффузии можно вычислить по формуле [3, с.
289]:
12 т• Шы
В ж = 7,4-10 42--,
ж /и • и
где Вж - коэффициент диффузии, м2/с; М - мольная масса растворителя, г/моль; и - мольный объём диффундирующего вещества, см3/моль; Т - температура, К;
Л - динамический коэффициент вязкости растворителя, мПа-с; в - параметр, учитывающий ассоциацию молекул растворителя.
Коэффициенты массоотдачи вг, вж отнесены к разности массовых концентраций, кг/м3. Чтобы вычислить коэффициенты массоотдачи вУ, вх , которые отнесены к разности массосодержаний, надо использовать формулы:
д. =А,- ^^ - , ' (51)
' Г к-т 1+(ы%„)г
Рх =рж- (1+X >(1+X •) • <52)
где ^=8310 Дж/(молыК) - универсальная газовая постоянная; Т - температура, К; П - давление, Па; - У + У
У = 1 2 - среднее массосодержание абсорбтива в газовой фазе, кг/кг инертного
газа;
X = Х1 + Х2 - среднее массосодержание абсорбата в жидкой фазе. Коэффициент т представляет отношение:
у*
т = —, X
которое определяется по линии равновесия. В формулу (46) подставляем среднее значение коэффициента:
У '(X,) У '(X,
т = 0,5- —+ —
V ^ X2 )
Рассчитываем по ф-ле (46) коэффициент массопередачи Ку.
(53)
Величину рабочей поверхности границы раздела фаз, через которую происходит
массопередача определим по формуле:
—
Fap = — • Щ , (54)
KY
где G - расход через колонну газа-носителя (инертного), кг/с;
Ky - коэффициент массопередачи, кг/(м2-с);
Fap - величина рабочей площади массопередачи, м2;
Y
J
dY
число единиц переноса.
Y - Y
Определенный интеграл при расчете пу вычисляется в пределах от У2 до У], так как при абсорбции содержание абсорбтива в газовой фазе уменьшается (У2<У]). Вычисление интеграла в общем случае возможно численным методом. Для этого надо построить
табличную зависимость подынтегральный функции —1—- от У, пользуясь уравнениями
У — У
рабочей и равновесных линий процесса, либо получая точки с графиков. Общий вид таблицы см. ниже.
Таблица
щ =
Y
2
Y Y2 Yn Yi Примечания
Имеется
зависимость по
X X2 Xn Xi рабочей линии процесса У = / (X)
Имеется
зависимость по
Y* Y'(X 2) Y'(Xn) Y'(X!) равновесной линии процесса У * = /(X)
Y - Y * Y - Y'(X2) Y - Y *(Xn) Y - Y *(X )
1 1 1 1
Y - Y * Y - Y'(X2) Y - Y *( Xn) Y - Y*(X )
п - промежуточные точки разбиения интервала [У2;У]].
На рис. 2 изображен процесс графического построения подынтегральной функции.
т'-т
Рисунок 2. Графическое изображение определения пу. 1 - Рабочая линия; 2 - Линия равновесия 3 - подынтегральная функция
Определенный интеграл графически эквивалентен площади под кривой 3 -заштрихованная криволинейная трапеция. Эту площадь можно найти численно методом трапеций по приближенной формуле:
п
£ 0,5 ■(/, + )-(у+1 - у.)
(55)
где К - число интервалов [У;У+1] разбиения интервала [У2У1].
Определив по (55) число единиц переноса, по формуле (54) определим рабочую площадь массопередачи. Она связана с геометрической площадью поверхности насадки по формуле:
^ = Н-8к-а-р, (56)
где произведение Н -Бк - а определяет геометрическую площадь поверхности насадки (Н - высота; Бк - площадь поперечного сечения; а - удельная площадь поверхности насадки, приходящаяся на единицу ее объема);
Ф - коэффициент активности насадки или смачивания. Он обычно меньше 1 (ф <1).
При оптимальной степени орошения насадки жидкостью ф ~1. Оптимальную степень орошения можно оценить по формуле:
и = Ь-а,
опт >
где а - удельная площадь поверхности насадки;
Ь - коэффициент, величина которого, например, при абсорбции аммиака водой составляет 0,158, а при абсорбции паров органических жидкостей водой принимается равным 0,093.
Надо сопоставить величину »ж с и опт:
при »ж < иопт величина ф <1;
при »ж ^ иопт величина ф -1.
к
В первом случае можно оценить ф по формулам [4, кн. 2, с. 489]:
р = 0,71 \рж- )13 " для колец размером 15-35 мм, загруженных
беспорядочно, где рж-^ =0,14-5,6 кг/(м2^с) - массовая скорость жидкости; о=30-70 мН/м -поверхностное натяжение жидкости; т = 1,42 • ё ~0'7, где ё - номинальный диаметр кольца, см.
р = А '{рж')0,435 •&" - для колец Рашига и седел Берля размером 12,5-25 мм,
загруженных беспорядочно, при плохой растворимости абсорбата, где п = Ь • ё ~р; для колец А=2,26, Ь=0,83, р=0,48; для седел А=0,767, Ь=0,495, р=0,98.
Из уравнения (56) можно определить рабочую высоту насадки:
Н = —. (57)
Бк • а • р
Для уточнения расчета следует вернуться к формуле (47), подставить полученное значение Н, и повторить расчет.
Диаметр колонны определим по формуле:
ёк =Л1 ^ . (58)
ж
Заключение
Методика расчета абсорбционных аппаратов для очистки воздуха от газообразных и парообразных загрязнителей может помочь освоить основные принципы работы аппаратов абсорбционной очистки, методы расчета процессов и конструктивных параметров аппаратов. Скорость газа принимают, чтобы она не превышала предел захлёбывания, когда режим подвисания определяется поднимающимся газом стекающей по насадке вниз жидкости при постепенном повышении скорости газа. При дальнейшем повышении скорости до предела подвисания начинается режим захлёбывания. При скоростях газа, превышающих предел захлёбывания, течение жидкости обращается вспять и начинается унос ее газом. Работа насадочной колонны в режиме начала подвисания считается оптимальной по образованию максимальной поверхности межфазного контакта и эффективности массоотдачи и оценить скорость газа в колонне, соответствующую началу режима подвисания.
ЛИТЕРАТУРА
1. 2.
3.
4.
5.
6. 7.
А.Н. Плановский, П.И. Николаев. Процессы и аппараты химической и нефтехимической технологии. - М.: Изд-во Химия, 2007, 493 с.
М.А. Адливанкина, М.Н. Борзых. Расчёт абсорбционных установок. Уч. пос., ГОУВПО «МГУС». - М, 2006, 65 с.
Павлов К.Ф., Романков П.Г., Носков А.А. Примеры и задачи по курсу процессов и аппаратов химической технологии. Уч. Пос. для вузов. 10-е изд. - Л: Химия, 1987, 576 с.
Гельперин Н.И. «Основные процессы и аппараты химической технологии». -М.: Химия, 1981. В двух книгах - 812 с., ил. - (серия «Процессы и аппараты химической и нефтехимической технологии»).
Лебедев В.В., Тыгер Л.М., Губанов Н.Н. Исследования нестационарных периодических процессов аэробной очистки сточных вод ЖКХ активным илом в циклотенках // Материалы Всероссийской заочной научно-практической конференции «АКТУАЛЬНЫЕ ПРОБЛЕМЫ ЭКОНОМИКИ И УПРАВЛЕНИЯ В СФЕРЕ ТУРИЗМА И СЕРВИСА» (проведенной в ФГБОУ ВПО «РГУТиС» 4 марта 2013 года). Сервис в России и за рубежом. 2014. Т. 8. № 4 (51). С. 117-129.
Губанов Н.Н., Чукалина Е.М., Иванов В.А. Совершенствование способов утилизации осадков сточных вод // Современные проблемы туризма и сервиса: Материалы Всероссийской научной конференции аспирантов и молодых ученых. - М.: ФГБОУ ВПО «РГУТИС», 2013 г. - 288 с.
Иванов В.А., Комаров Н.М., Крымская Е.Я., Панова М.В. Водные ресурсы России, модели метода их сохранения и вызовы проекта // [Электронный ресурс]: http://publ.naukovedenie.ru Интернет-журнал «НАУКОВЕДЕНИЕ» №6 (19). Выпуск 6, ноябрь - декабрь 2013 (дата обращения: 05.03.2014).
Birger Evgeniy Vladimirovich
Institute of World Economy and Informatization, Russia, Moscow
E-mail: [email protected]
Lebedev Vladimir Vladimirovich
Russian State University of Tourism and Service, Russia, Moscow
E-mail: [email protected]
Zizyuk Aleksey Petrovich
Russian State University of Tourism and Service, Russia, Moscow
E-mail: [email protected]
Ivanova Elena Vladimirovna
Centre for Managing Problems of Resources Sparing and Waste Management, Russia, Mytishchi
E-mail: [email protected]
Bases of designing of absorbers for clearing of air from gazo - and vaporous impurity
Abstract. Dividing массообменные processes are widely applied in the industry. These processes, as a rule, carry out in devices of columned type.
They find wide application and in the technician of clearing of air in inhabited, public and industrial premises. Application абсорбционной clearings of air of gaseous impurity at the service enterprises is actual.
For designing of such installations the effective techniques accessible to application by technical experts are necessary.
In the present work the design procedure абсорбционных devices for clearing of air of gaseous and vaporous pollutants is presented.
Work can help experts and trained to master main principles of work of devices абсорбционной clearings, methods of calculation of processes and design data of devices.
Keywords: quantity of substance exchange processes; division of gas mixes; absorption; absorbent; absorbtiv; absorbat; border of section of phases; absorbers; a working line of process; number of units of carrying over; a factor quantity of substance transfer; the area of an interface of phases
REFERENCES
1. 2.
3.
4.
5.
6. 7.
A.N. Planovsky, P.I. Nikolaev. Processes and devices for chemical and petrochemical technology. - M: Publishing house Chemistry, 2007, page 493.
M.A. Adlivankina, M.N. Borzyh. Calculation of absorbent installations. Textbook of SEIHPE «MSSU». - M, 2006, page 65.
Pavlov K.F., Romankov P.G., Noskov A.A. Samples and tasks for study of chemical technology's processes and apparatus. Textbook for high schools. 1'th edition. - l: Chemistry, 1987, page 576.
Gelperin N.I. «the Basic processes and devices of chemical technology». - M: chemistry, 1981. In two books - page 812, silt. - (Series «Processes and devices of chemical and petrochemical technology»).
Lebedev V.V., Tygler L.M., Gubanov M.M. Research of non-stationary periodic processes of aerobic wastewater purification by housing and public utilities by the use of active sludge cyclotenks. // Materials from All-Russian extramural scientific-practical conference "ACTIAL PROBLEMS OF ECONOMICS AND MANAGEMENT IN SPHERE OF TOURISM AND SERVICE". (Event took plave in FSBEI HPE "RGUTiS" from 04 March 2013) Service in Russia and abroad 2014 T. 8. № 4 (51). S. 117-129.
Gubanov N.N., Chukalina E.M., Ivanov V.A. Modernization of methods for utilization of wastewaters remainders // Modern problems of tourism and service: Materials from All-Russia scientific conference of post-graduates and young scientists - in FSBEI HPE "RGUTiS", year 2013 - page 288.
Ivanov V.A., Komarov N.M., Krymskaya E.Ya., Panova M.V. Russia's water resources, models for methods of their saving and the project's challenges // "Electron resource": http://publ.naukovedenie.ru, Internet-hosted magazine "Naukovedenie", №6 (19). Issue № 6, November - December 2013 (accessed date: 05.03.2014).