УДК 541.183
ОПИСАНИЕ И АПРИОРНЫЕ РАСЧЕТЫ РАВНОВЕСИИ АДСОРБЦИИ СМЕСЕЙ ФЛЮИДОВ НА АКТИВНЫХ УГЛЯХ Т.А. Кузнецова, Н.Г. Крюченкова, М.В. Бородулина, Е.Н. Егоров, А.М. Толмачев
(кафедра физической химии; e-mail: tolmach@physch.chem.msu.su)
На примере новых экспериментальных данных по адсорбции бинарных растворов неэлектролитов на активном микропористом угле показана возможность использования системы уравнений, полученных ранее на основе модели Оно-Кондо, для количественного описания изотерм избыточной адсорбции компонентов растворов. Подтверждена возможность использования «метода характеристических кривых» для априорных расчетов равновесий адсорбции смесей флюидов.
Полученные ранее в рамках модели Оно-Кондо-Арановича [1] системы уравнений:
(1)[1 - ^ (-)] - в - А + N1 (- )[1 - N1 (1)] кТ 1
+ 3N 1 (2) - т 1 (-)] = 0 , (1)
Г 1, рас. = а1,ш I [ 1 (п) - N 1 (- ) ], (2) п=1
Г 1, эксп. = Ш0^01 - N 1 (- . (3)
были использованы для количественного описания изотерм избыточной адсорбции компонентов бинарных растворов неэлектролитов на активном угле БАУ-1, полученных нами при 293 К по описанной в [2] методике для замкнутой группы растворов: изооктан - толуол, толуол - СС14, СС14 -изооктан.
Алгоритм численного решения системы уравнений (1-3), описанный в [1], позволяет находить параметры уравнений:
B _ _L(n) £ 02 £ 01 + £ 22 £ 12 )
kT
А _ 2 £ 12 _ £ и _ £ 22
ao=2ai,n
(4)
(5)
(6)
и сравнивать зависимости рассчитанной и экспериментальной избыточной адсорбции (Г1) от состава равновесного раствора )).
В уравнениях (1-6): ), ^(п) - мольные
доли первого компонента в исходном и равновесном растворах и в соответствующем адсорбционном слое; £у, £0. - энергии взаимодействия адсорбат-ад-сорбат и адсорбат-адсорбент; а1ш - емкость монослоя (моль кг ); а0 - предельная величина адсорбции по ТОЗМ [3]; ш0 - общее количество молей в исходном растворе; g - навеска адсорбента (кг); Ь(п) - коэффициент, зависящей от числа слоев ад-сорбата в щелевидной поре (для использованной нами двухслойной модели адсорбционного раствора Ь(п) = 1,12). Результаты расчетов приведены в таблице и на рис. 1.
Как видно из приведенных данных, система уравнений (1-3) количественно описывает все три избыточных изотермы адсорбции, несмотря на то, что они относятся к разным типам изотерм. При этом в соответствии с [1] значения энергетических параметров (4) и (5) имеют физически достоверные значения. Это дополнительно подтверждается сравнением значений (е02 - £01)'10 (Дж/молекула), рассчитанных нами из значений костант В (величины £у находили по значениям теплоты испарения чистых адсорбтивов). Очевидно, что сумма этих разностей для циклически замкнутого набора систем должна быть равной нулю. Для приведенной выше последовательности изученных нами систем эта сумма равна:
1,14 - 0,91- 0,25= -0,02,
т.е. действительно близка к нулевому значению.
Однако рассчитанные значения емкости монослоя (и соответственно предельной величины адсорбции
Значения параметров А/кТ, В, а01 для различных адсорбционных систем
№ Система Т, К Д/кТ В ао1, рас. моль/кг ао1, эксп. моль/кг
1 Толуол(1) - изооктан(2) 293 0,28 1,91 2,28 2,45
2 Толуол(1) - СС14(2) 293 0,54 1,54 0,87 2,45
3 СС14(1) - изооктан(2) 293 0,06 0,95 1,68 2,69
а0) заметно отличаются от найденных независимо по изотермам адсорбции соответствующих паров, что, как было отмечено в [1], является следствием различия модели (молекулы одинаковых размеров) и реальных систем с разными размерами молекул компонентов.
В работе [4] подробно рассмотрен «метод характеристических кривых» как универсальный метод априорных расчетов адсорбционных равновесий в бинарных и многокомпонентных системах. Метод основан на уравнениях подобия
Т
1 "V
'з 1пгхух Л
д N
0ху ,Тх
— ТЛВ Ьз
д 1пг А У,
д N
0 АВ ,ТАВ
(7)
Т 1 V
д|п/ XУз
д Nх
0 ху ,Тху
— Тав(Ь3
гд1пгАуА л д N л
1
| Ьх
0 АВ , ТАВ
(8)
Nх — NA
тем. Выбор уравнения подобия проводится по значениям параметра
м —
1х 1а
I,
(11)
Г _МаЕ 01Нху [N + Р АВ (1 ~ , (9)
х 1х1вЕо2 Р AB[NA + Р ху(1 - '
связывающих производные характеристических кривых первых компонентов, т.е. зависимостей вида
ЯТЬ^ = ад (10)
для рассчитываемой (X, У) и базовой (А, В) сис-
Если М < 0,3 (что имеет место для большинства систем), необходимо использовать уравнение (7), а в остальных случаях - уравнение (8) [4]. Если же значения Ьх при всех значениях N1 близки к единице, то оба уравнения позволяют получить близкие результаты. Важно подчеркнуть, что все расчеты проводят на основании данных по адсорбции индивидуальных компонентов смесей и данным для базовой системы.
В уравнениях (7-11): - температура для рссчитываемых и базовых систем; г,у, = Р|У|/( Р1У1)ст.(пары) = х,у, (растворы), р,, х, - равновесное давление и мольная доля, р, ст. - давление насыщенного пара соответствующего чистого компонента, у, - коэффициент активности компонентов объемной фазы; N1 = Рцс,/ (РуС, + с^ - объемные доли (с, - молярная концентрация в рамках метода полного содержания) компонентов адсорбционного раствора; Ру = а^ / а0, - стехиометрический коэффициент взаимного вытеснения компонентов адсорбционного раствора [4]; 0у - суммарная степень заполнения объема микропор (при адсорбции из растворов 0у = 1); I, - коэффициент подобия характеристических кривых адсорбции соответствующих паров по стандартному пару (бензол), Е01, Е02 - характеристическая энергия адсорбции стандартного пара на адсорбентах соответственно для базовой и рассчитываемых систем по ТОЗМ [3].
Подробный алгоритм пересчета избыточных величин адсорбции в «абсолютные» и последующего
Ь
у
N V — N
хА
Г1, ммоль/г
Г1, ммоль/г
1,5 -|
0,5
-0,5
1
0,8 -
0,6 0,4 0,2 0
0 0,2 0,4 0,6 0,8 1
Рис. 1. Зависимости от ^(ю) для адсорбции на БАУ-1 при 293 К первых компонентов растворов; а: • - (толуол -изооктан); Д - (толуол - СС14); б: • - (СС14 - изооктан)
0 0,2 0,4 0,6 0,8 1
-500 -
-1000
-1500
-2000
-2500
Кх
Т1п2хГх
-500 -1000 -1500 -2000 -2500 -3000
0 0
Т1шхГх
Рис. 2. Экспериментальные (•) и рассчитанные по (7) (°) и (8) (Д) характеристические кривые компонента Х для систем: а: толуол - СС14 - БАУ-11, 438 К, 0Ц = 0,85(Х,У) [5], базовая система - (толуол - СС14 - БАУ-1), 293 К(А, В); б: толуол - СС14 - БАУ-1, 293 К(Х, У), базовая система - (бензол - СС14 - АУ-6), 303 К(А, В) [6]
априорного расчета характеристических кривых компонентов (ХККС) Х и У приведен в [4].
На рис. 2, а результаты расчета ХККС сопоставлены с полученными ранее экспериментальными значениями для адсорбции смеси паров при 0.. = 0,85 (значения Тху1пгх при N = 1 и 0.. = 0,85 нахо-
ху XX
дили по изотерме адсорбции соответствующего пара). Поскольку в данном случае значения Ьх близки к единице, оба уравнения подобия позволяют получить близкие результаты. Удовлетворительное
согласие рассчитанной и экспериментальной ХККС подтверждает возможность расчета равновесий адсорбции из паровой фазы по данным для адсорбции из жидкого раствора, полученным при другой температуре на активном угле с другим значением Е01.
На рис. 2, б сопоставлены ХККС, найденные нами экспериментально для компонента X и рассчитанные по данным для базовой системы (получены в [6] для адсорбции из раствора на другом активном угле при другой температуре). В этом случае
б
а
1
0
0
Ьх заметно отличается от единицы, а параметр М = 0,15. Из рис. 2 видно, что использование уравнения (7) позволяет получить удовлетворительное согласие ХККС, что является дополнительным под-Работа выполнена при поддержке РФФИ, п]
тверждением возможности проведения априорных расчетов адсорбционных равновесий в бинарных системах в рамках «метода характеристических кривых».
ект № 02-03-32641 и проекта НШ-1275.2003.3.
СПИСОК ЛИТЕРАТУРЫ
1. Толмачев А.М., Бородулина М.В., Арзамасцева А.Б., Крючен-
коваН.Г., Рябухова Т.О. // Вестн. Моск. ун-та. Сер. 2. Химия. 2001. 42. С. 244.
2. Крюченкова Н.Г., Кузнецова Т.А., Бородулина М.В., Егоров
Е.Н., Толмачев А.М. // ЖФХ. 2001. 75. С. 1145.
3. Dubinin M.M. // Progress in surface and membrane sci. N.Y.,
1975. 9. P. 1.
4. Толмачев А.М., Пронина Е.А., Бородулина М.В., Егоров Е.Н.,
Крюченкова Н.Г. // Вестн. Моск. ун-та. Сер. 2. Химия. 2002. 43. С. 139.
5. ТрубниковИ.Б. // Дис. ... канд. хим. наук. М., 1982.
6. Ларионов О.Г. // Дис. ... докт. хим. наук. М., 1975.
Поступила в редакцию 02.03.03
DETERMINATION AND A PRIORI CALCULATIONS OF ADSORPTION EQUILIBRIA FOR FLUID MIXTURES ON ACTIVATED CARBONS
T.A. Kuznetsova, N.G. Kruchenkova, M.V. Borodulina, E.N. Egorov, A.M. Tolmachev
(Division of Physical Chemistry)
New experimental data of adsorption for binary non-electrolyte solutions on activated carbons is obtained. The opportunity of using the system of equations based on Ono-Kondo model for quantitative determination of excess adsorption isotherms for components of solutions is shown. The opportunity of using the method of characterized curves for apriori calculations of adsorption equilibria for fluid mixtures is proved.