Теоретическое исследование массоотдачи в многокомпонентных газовых смесях Текст научной статьи по специальности «Математика»

Л. Э. Осипова, Э. Ш. Теляков, Д. А. Аминева,

Б. Н. Матюшко

ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ МАССООТДАЧИ

В МНОГОКОМПОНЕНТНЫХ ГАЗОВЫХ СМЕСЯХ

Ключевые слова: Многокомпонентная диффузия, массоотдача и массопередача, гетерофазная система, математическое

описание.

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

Keywords: Multicomponent diffusion, returns the masses and the transfer of the masses, heterophase system, the mathematical description of.

On the basis of the serial output and analysis of integral equations of diffusion boundary layer for the gas phase system of gas (steam) - liquid justifies-the structure of the equation returns the masses, which is based on the allocation of a process to a specific standard conditions (equimolar process) and differs from known phenomenological equations of the presence of additional members, taking into account the effect on the mass transfer of the resulting mass flow through the surface of the phase. Justice of the equation is confirmed by the analysis of the known experimentexperimental data.

Введение

В основе описания многокомпонентного массопереноса в системах газ (пар) - жидкость практически всегда лежат уравнения диффузии, а переход от уравнений диффузии к уравнениям массоотдачи реализуется с использованием каких-либо из теорий массопередачи [1,2,3]. Невзирая на многолетнюю историю исследований данной проблемы, она до сих пор вызывает неослабевающий интерес у исследователей. Механизм переноса массы в многокомпонентных системах достаточно сложен. Так для многокомпонентных смесей (МКС) наличие (или отсутствие) собственного градиента концентраций для рассматриваемого компонента далеко не полностью определяет механизм его переноса. При разделении МКС возникают такие специфические явления, как диффузионный барьер, осмотическая и реверсивная диффузия [4]. Было показано [5], что при описании диффузии в МКС существенными преимуществами обладают матричные методы, которые обеспечивают высокую точность расчетов процессов переноса и, одновременно, - компактность математического описания. Применительно к массопе-редаче вопросы описания процессов переноса в МКС менее развиты и до сих пор основываются на феноменологических подходах [6]. Поэтому представляется целесообразным использовать матричные методы и при описании процессов массопере-дачи.

Теоретический анализ процесса массоотдачи

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

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

N = N° + ЫК (1)

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

ю.

3 3у 3 / ч5у

3 +юг-*- = — (D + sd)- 3

(2)

дг г дг дг' дг

При записи (2) сделаны допущения: плотность жидкости в слое постоянна, разделяемая смесь

- бинарная, баро- и термодиффузия отсутствуют. Дополним (2) уравнением неразрывности:

а». а», = 0 (3)

5Z

»г Мс

^»Г = 0

дг

(4)

(5)

и граничными условиями:

г = 0 у = уг

Г = ш У = Уш

Тогда из (3) имеем:

Г дю,

аг

<3г

о

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

Тогда:

ду

С»7----------+

2 дг

(

Л

аг=(у+Бо)^

д

(6)

Правая часть (6) определяет внутренний диффузионный поток переносимого компонента от поверхности раздела в ядро потока:

.ду

у + Ер )•

д

= -N010 =|^

о, г

(7)

С учетом (7) проведем преобразование уравнения (6):

ш ш ~ (г „ Л ш

Гд»

= с!г-Г— Г

.) 2 & л д л дг

0 0 Vo

^с1г

(8)

Второй интеграл правой части уравнения (8) (ь ь А

может быть взят по частям

|ис1у = иу|^ - Г vdu

^ а

Тогда:

(г - А

д

0 Vo

= У,

(

Гд й 4г=1й с,гЧ1т2с|г

oVo

-М

с1г=

дю2

~дТ

= /У„-у))^г

дг

с1г =

(9)

Подставим (9) в (8). Поскольку в рассматриваемом случае пределы интегрирования не зависят от продольной координаты z, последовательность операций дифференцирования по z и интегрирования по г может быть изменена. Тогда, учитывая, что N.. является константой, окончательно

получим:

|с»§с1г-|(уш- у))гс|г -М с^-УшЬ

д»

дг

со

с1 Г

^•|®г(Уш -у)- Nс(уг -у*,)

(10)

Полный поток переносимого компонента должен учитывать также перенос массы за счет кон-

векции (1). Тогда:

ш

^-^ЮУ-у)-^ -Уш)+ЧУг (11)

0

Уравнение (11) называется интегральным уравнением диффузионного пограничного слоя. В (11) составляющие, входящие в интегральный член, зависят только от продольной координаты z. Функции »» {т) и у{т) для многих практически важных случаев достаточно хорошо изучены, что позволяет непосредственно использовать это уравнение для расчетов процессов массопереноса. При использовании феноменологического подхода это уравнение раскрывает физический смысл коэффициентов переноса (массоотдачи) и может применяться для определения соответствующих коэффициентов.

Соотношения, подобные (10), были получены впервые для теплового (Г.Н.Кружилин) и для гидродинамического (Т.Карман) пограничных слоев. Детальный вывод интегрального уравнения теплового пограничного слоя представлен в [7]. В отличие от этих работ при выводе (11) учтено наличие ненулевой составляющей поперечной скорости на внешней границе потока, что отражено в граничных условиях (4). Учет этого обстоятельства привел к появлению в интегральном уравнении диффузионного пограничного слоя дополнительного члена МУг-Уш). Как нам представляется, это обстоятельство должно учитываться уже при обосновании самой структуры феноменологических уравнений переноса, поскольку во многих случаях наличие ненулевой поперечной компоненты скорости на внешней границе пограничных слоев может достаточно серьезно сказываться на распределении функций ю2 {г) и Ууг).

Действительно, при описании процессов массоотдачи в бинарных смесях наиболее широкое распространение получили феноменологические уравнения вида:

^ = Р-У (12)

Физический смысл, вкладываемый в понятие коэффициента переноса, следует из сопоставления интегрального уравнения диффузионного пограничного слоя (11) с феноменологическим уравнением (12). Для эквимолярного процесса (Ыс = °) получим:

ш

^ |»уш

Р* =-----0--- -------- (13)

Ду

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

ш

ш

0

0

0

ш

0

0

0

ш

0

0

сообразным принять эквимолярный процесс:

ВД

сЫ

н=-

0

_с[

dz<

-р*Ду-1^(у-Увд)+1С (14)

—|ю*2(увд-у^г

0

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

ВД

А|0,2у„-У^г

—0-------------= Ф

^{“НУвд -у^г

Бс,

Р^ес

(15)

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

В качестве характеристической движущей силы процесса в уравнении (13), как правило, используется выражение Ду ^ уг ^ у . Это объясняется тем, что в большинстве практически важных случаев (например, при обобщении опытных данных с помощью аппарата теории подобия) исследователь не может оперировать фактической разницей концентраций на границах диффузионного слоя (Ду ^ уг ^ уш ). Можно утверждать, что практически все обобщающие уравнения по массоотдаче [2,8] получены с использованием среднерасходных концентраций. В то же время второй член в правой части интегрального уравнения диффузионного пограничного слоя (14) включает в себя вполне конкретную движущую силу (Ду ^ уг ^ уш ). Учитывая, что рассматриваемые движущие силы связаны между собой:

У =

Р =

Ю2ус1г

Ю,СІГ

(16)

(17)

Уг - У

уравнение (14) можно представить в окончательном виде:

N =с{5с;'^)рУ -у)-рН(уг -у)+Мс^ (18)

Следует отметить, что интегральные уравнения диффузионного пограничного слоя могут быть получены и другими приемами, помимо пока-

занного выше. Например, они следуют уже из балансовых соотношений для установившегося пограничного слоя. Поэтому можно утверждать, что структура уравнения неэквимолярной масоотдачи (18) для бинарной смеси принципиально не зависит от условий обтекания поверхности раздела фаз и должна быть универсальной. Анализ предельных вариантов формирования профилей концентраций в диффузионном пограничном слое показывает, что параметр р всегда лежит в интервале 1 < р < 2 и для большинства практически важных случаев может быть принят равным 1,5.

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

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

ІИІВМІЬІ

(19)

Непосредственное использование этого уравнения не позволяет получить структуру уравнения многокомпонентной массоотдачи. Соответствующее решение однако может быть получено при использовании приема линеаризации уравнений (19). Процедура линеаризации уравнений молекулярной диффузии достаточно подробно описана в работах [4,9]. При переходе к уравнениям конвективной диффузии задача принципиально не меняется. Тогда линеаризацией уравнения диффузии (19) с использованием диагонализирующей матрицы получим:

ду

дг

дг

+8г

дг

д_

дг

(20)

Как видим, система уравнений многокомпонентной конвективной диффузии распадается при линеаризации на ряд линейных дифференциальных уравнений, в которых вместо концентраций у выступают их линейные комбинации /, а вместо элементов матрицы многокомпонентной диффузии -элементы диагональной матрицы . Поскольку по

форме полученные дифференциальные уравнения совершенно аналогичны уравнению конвективной диффузии в бинарных смесях (2), можно использовать выше найденное решение (18):

(21)

Пренебрегая поправочной функцией (15) и проведя обратные преобразования, окончательно

ВД

0

ВД

0

и

для многокомпонентной смеси получим:

Н)=[в* - у)-рЧу - уМсЫ (22)

Элементы матрицы многокомпонентой мас-соотдачи могут быть рассчитаны разными способами [9]. Более предпочтительным представляется применение уравнения, предложенного специально для описания многокомпонентной эквимолярной массоотдачи [10]:

т г ,

К = £ру(уЫ -У]^)=[ ДУг -У) (23)

И

т

где вн=^рУ] (24)

И

и В = -ур, (25)

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

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

Общее математическое описание сопряженного процесса должно включать в себя:

• уравнения материального баланса (общего и покомпонентного) для каждой из фаз системы;

• уравнения теплового баланса для обеих фаз;

• уравнения массо - и теплоотдачи для обеих фаз;

• уравнения сопряжения потоков на границе раздела фаз (балансовые соотношения по значениям потоков и условия обеспечения термодинамического равновесия на границе раздела фаз).

Для записи уравнений массоотдачи в жидкой фазе использовать структуру (22) можно только в первом приближении. Действительно, при выводе этого уравнения для газовой фазы использован аппарат кинетической теории идеальных газов, и в частности - условие независимости бинарных коэффициентов диффузии от состава газовой смеси. Однако для жидкой фазы зависимость коэффициентов диффузии от состава носит более выраженный характер. Более того, в качестве справочных данных

по коэффициентам диффузии в литературе в основном используются коэффициенты диффузии для разбавленных растворов [1,2,8].

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

Следует также отметить, что неэквимоляр-ность процесса приводит к трансформации и феноменологического уравнения теплоотдачи, причем структура этого уравнения может быть получена также из интегрального уравнения теплового пограничного слоя [11].

Ещё один аспект разработанного описания связан с векторной формой записи уравнения мас-соотдачи (22). За положительное направление потока при записи принято направление от внешней границы потока (поверхность раздела фаз) к ядру потока. При описании процесса массопередачи, это обстоятельство необходимо учитывать, поскольку в этом случае необходимо закреплять направление уже межфазного потока: например, принимать за положительное направление поток массы от одной фазы (например, жидкой) - к другой (газовой).

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

Экспериментальная часть

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

Shy=ARemScn, (26)

причем при расчете критерия Shy используются

коэффициенты массоотдачи, рассчитанные по полному массовому потоку (12).

Отмечалось [12], что данная форма обобщения дает достаточно существенные средние отклонения для разных смесей (до ±30%). Поэтому были предприняты попытки повысить точность обобщения за счет использования более строгих обобщающих уравнений. Так в работе [13] была построена

специальная математическая модель внутрифазного эквимолярного массообмена в осесимметричном турбулентном газовом потоке и проведен специальный численный эксперимент. В модели уравнения (2), (3) были дополнены выражениями для расчета профилей скоростей и турбулентной вязкости в пограничном слое. Поэтому обобщения, полученные в результате моделирования, можно рассматривать как теоретические решения уравнения (13). При обобщении была использована форма критериального уравнения, предложенная ещё фон Карманом, и получено:

0,0345Ре075 Эс

(27)

Shy =

1+0,7 5Re-0,09{Sc -1 +ln [ +5Sc)/6]}

Сопоставление наиболее надежных критериальных уравнений по данным [2] с уравнением (27) представлено на рисунке. На этом графике данные (ряд 1 - ряд3) характеризуют процесс испарения чистого вещества в среду инертного газа (суммарный поток направлен от поверхности раздела в глубь газовой фазы); (ряд 4 - ряд 9) - процесс пленочной ректификации (примерно эквимолярные условия); (ряд 10 - ряд 11) - процесс абсорбции (суммарный поток направлен к поверхности раздела фаз). Данные, соответствующие ряду 12, получены по уравнению (27). Все расчеты сделаны для фиксированного значения критерия 8с=0,75.

5Ьт

-*-Ряд1 -¿^РядЗ -*-Ряд4 -Ж-Ряд5 —Рядб

—1—Ряд7 — Ряд8 РядЭ —♦-РядЮ —■— Ряд11 -*-Ряд12

Рис. 1 - Сопоставление литературных данных по обобщению массообменных характеристик для пленочных аппаратов: испарение чистых веществ в инерт (ряд1 - [14], ряды 2 и 3 - [2]); ректификация (ряд4 [2], ряды 5 и 6 - [15]; ряд7 - [16], ряды 8, и 9 - [2]); абсорбция (ряды 10 и 11 - [2], ряд12 - расчет по уравнению (27))

Как видим, данные по обобщению процессов ректификации лежат примерно в середине общего поля, обобщающего экспериментальные данные. При этом обобщения для процесса испарения расположены выше линий обобщения процесса рек-

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

Выводы

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

Обозначения

В - элементы квадратной матрицы коэффициентов многокомпонентной эквимолярной массоотдачи; BL - элементы диагональной матрицы коэффициентов массоотдачи, D, D и Dl - соответственно коэффициенты бинарной диффузии, элементы квадратной и диагональной матриц многокомпонентной диффузии, м2/с; N - поток вещества, м3/(м2 с) y - мольная концентрация газа; Ay -движущая сила процесса массоотдачи; Г иг - поперечная и продольная координаты соответственно; Р - коэффициент массоотдачи бинарной смеси, м/с; SD - коэффициент турбулентной диффузии, м2/с; q и у - псевдопотоки и псевдоконцентрации линеаризированного уравнения многокомпонентной массоотдачи соответственно; О

- скорость, м/с; Re = О / v - критерий Рейнольдса; Sc = v/D - критерий Шмидта; Sh = pi/D - критерий Шервуда; [], (6 и( ) - соответственно квадратная, столбцовая и диагональная матрицы.

Индексы

D и K - диффузионная и конвективная составляющие массового потока соответственно; i, j -

компоненты; f - поверхность раздела фаз; с -суммарный поток; z, r - направления осей координат; да - ядро потока; * - эквимолярные условия;

- среднемассовые величины.

Литература

1. В.В.Кафаров, Основы массопередачи. М.: Высшая школа, 1972, С.494.

2. В.М.Рамм Абсорбция газов. М.:Химия, 1971, С.767.

3. Качалова Т.Н., Качалова Л.И. Модернизация узла гидрирования этан-этиленовой фракции // Вестник Казан. технол ун-та, - 2012, №10, С.230.

4. Toor H.L. Diffusion in three-component gas mixtures // AIChE Journal. 1957. V. 3. № 2. Р. 198.

5. Теляков Э.Ш., Низамов И.И., Таренко Б.И., Бикмурзин А.Р. Экспериментально-расчетное исследование многокомпонентной диффузии в газовых смесях // Бутлеров-ские сообщения. - 2011. - Т.27. - № 15. - С.71-79.

6. Лаптева Т.В., Зиятдинов Н.Н., Первухин Д.Д. Эффективность работы методов решения задачи проектирования работоспособных ХТС // Вестник Казан. технол унта, - 2012, №11, С.268-271.

7. В.П.Исаченко, В.А.Осипова, А.С.Сукомел. Теплопередача, Москва, Энергоиздат, 1981, С.416.

8. К. Ф. Павлов, П. Г. Романков, А.А. Носков Примеры и задачи по курсу процес-сов и аппаратов химической технологии. Изд. 10-е, пер. и доп. Л., Химия,1987. -С.576.

9. И.А. Александров Массопередача при ректификации и абсорбции многокомпонентных смесей. Л.: Химия, 1975, с.320.

10. Константинов Е.Н., Николаев А.М. Массопередача при ректификации многокомпонентных смесей // Известия ВУЗов, Нефть и газ. 1964. № 1. С. 53.

11. Осипова Л.Э., Теляков Э.Ш., Николаев Н.А. Моделирование неэквимолярных процессов массо- и теплопе-

реноса в системах газ (пар) - жидкость // Теорет. основы хим. технологии. 1997. Т. 31, № 3, С. 274.

12. Лотхов В.А., Малюсов В.А., Жаворонков Н.М. Определение фазовых сопротивлений в процессе пленочной ректификации по профилю концентраций. Ламинарный поток пара. Теорет. основы хим. технол., 1967, т.1, №3, с.319-322.

13. Константинов Е.Н., Кузнечиков В.А. Математическая модель турбулентного массообмена в многокомпонентной смеси. Теорет. основы хим. технол., 1975, т.9, №2, с.163-169.

14. Von Krevelen D.W., Hoftijzer P.J. Studies of gas absorption: mass transfer in a film reactors (wetted wall columns). Rec. Trav. Chem., 1949, v.68, №4, p. 221-242/

15. Боярчук П.Г. Исследование кинетики массообмена в ректификационных аппаратах пленочного типа. Дисс.... канд. технических наук, М.: МИХМ, 1963, С.164.

16. Теляков Э.Ш., Николаев А.М. Способ определения частных высот единиц переноса в процессе ректификации. Теорет. основы хим. технол., 1969, т.3, №6, с.826-830.

© Л. Э. Осипова - кнад. техн. наук, доц. КГАСУ; Э. Ш. Теляков - д-р техн. наук,проф. каф. машин и аппаратов химических производств КНИТУ; Д. А. Аминева - инж. Волжского научно исслед. института углеводородного сырья (ВНИИУС), [email protected]; Б. Н. Матюшко - канд. техн. наук, зав. лаб. ВНИИУС.