Адсорбция газов, паров, растворов
.pdf
|
|
|
|
111 |
|
|
|
A = − |
p |
s |
|
f |
s |
|
(5.24) |
G = RT ln |
|
= RT ln |
|
||||
|
|
|
p |
|
|
f |
|
где p – равновесное давление, а f – летучесть пара при температуре Т. Введение в формулу летучести вместо давления позволяет учесть
неидеальность газовой фазы.
Если адсорбция выражена в безразмерных единицах, то и дифференциальную мольную работу адсорбции целесообразно представить
в форме безразмерного отношения А/Е0, где Е0 |
— характеристическая |
энергия адсорбции. |
|
Тогда термическое уравнение адсорбции можно представить в общей |
|
форме: |
|
θ = f[(A / E0 ), n] |
(5.25) |
Уравнение (5.25) выражает функцию распределения микропор по дифференциальной мольной работе адсорбции, причем Е0 является одним из параметров этой функции. Большинство функций распределения в нормированной форме характеризуется двумя параметрами. В формулу также введен второй параметр, который условно обозначен через n.
Согласно уравнению (5.25), мы получаем выражение для так называемой характеристической кривой:
A = E0ϕ(θ, n) |
(5.26) |
Если для различных паров функция φ и параметр n остаются неизменными, то
A1 A2 = E01 E02 =β (5.27)
т. е. характеристические кривые в координатах А - θ являются «подобными». Другими словами, при одинаковых значениях θ отношения ординат характеристических кривых для двух паров постоянны и равны коэффициенту подобия β в том интервале изменения заполнения θ, в котором сохраняются исходные допущения неизменности функции φ и постоянства ее параметра n.
112
Из уравнения (5.26) следует, что Е = A для некоторого заполнения θ0 или характеристической точки, определяемых в общем случае из условия
ϕ(θ0 , n) = 1 |
(5.28), |
причем при неизменности функции φ заполнение θ0 будет одинаковым для различных паров. В связи с этим возникает возможность определения характеристической свободной энергии адсорбции по одной точке изотермы адсорбции, соответствующей заполнению θ0, согласно уравнению (5.28). Естественно, что абсолютное значение θ0 зависит от вида функции φ.
Рассматривая в термическом уравнении адсорбции (5.25) f как функцию распределения, мы, по существу, приняли допущение о температурной инвариантности этой функции, полагая, что ее параметры Е0 и n являются постоянными величинами для рассматриваемой адсорбционной системы. Так как Е0 = А для заполнения θ0, т. е. Е0 является одной из точек характеристической кривой, то допущение о температурной инвариантности автоматически приводит к независимости характеристической энергии адсорбции Е0 и, как следствие, параметра n от температуры.
Основное уравнение теории (5.25) в самом общем виде выражает распределение заполнения микропор θ по дифференциальной мольной работе адсорбции А, причем параметры распределения Е0 и n не зависят от температуры, если соблюдается температурная инвариантность характеристических кривых. Принимая температурную инвариантность уравнения (5.25) и учитывая известное в математической статистике распределение Вейбула, Дубинин и Астахов получили термическое уравнение адсорбции в аналитической форме:
|
|
n |
(5.29) |
|
θ = exp |
− A |
|
|
|
|
|
E0 |
|
|
Для большого числа адсорбционных систем в случае активных углей параметр n=2, а в случае цеолитов 3. Уравнение (5.29) позволяет перейти к уравнению характеристической кривой:
|
113 |
|
|
A = E0 [ln(1 |
θ |
)]1/ n |
(5.30) |
|
|
|
Далее из уравнения (5.30) следует, что при заполнении θ0 = 1/e = 0,368 характеристическая энергия адсорбции Е0 равна дифференциальной мольной работе адсорбции А.
Выражая в уравнении (5.29) степень заполнения через величины адсорбции из θ = aa0 , получим термическое уравнение адсорбции:
a = a0 |
exp |
− (A |
E |
)n |
|
|
|
|
|
(5.31) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
или в линейной форме: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
n |
= lna |
|
RT n |
ln |
n P |
(5.32), |
|||||
lna = lna − |
|
|
|
A |
|
|
− |
|
|
|
s |
|||||
0 |
E |
n |
|
|
|
|
|
0 |
|
E0 |
|
|
|
P |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
причем отрезок, отсекаемый прямой на оси ординат, равен lnа0, а угловой
коэффициент прямой равен (RT/E0)n. |
Если |
показатель |
степени |
известен (2 для активных углей и 3 для цеолитов), |
то из графика легко |
определяются предельная величина адсорбции а0 и характеристическая энергия адсорбции Е0 на основании одной экспериментальной изотермы адсорбции (см. рис. 5.4.). Отклонения от линейной зависимости в области малых величин адсорбции связано с тем, что (5.32) не выполняется при a / a0 < 0.1, а отклонения в области больших величин адсорбции проявляются в случае активных микропористых углей с заметно развитой поверхностью мезопор за счет дополнительной адсорбции на этой поверхности при относительных давлениях близких к единице.
|
|
114 |
|
|
lna |
|
|
|
|
1.5 |
|
|
|
|
0 |
|
|
|
|
–1.5 |
|
|
|
|
1 |
2 |
3 |
4 |
–ln2(PS/P) |
Рис 5.4. Схема определения параметров уравнения Дубинина-Радушкевича Уравнение (5.32) является эмпирическим, поскольку основано на
допущении о возможности использования распределения Вейбула для описания распределения заполнения микропор θ по дифференциальной мольной работе адсорбции А. Более того оно термодинамически некорректно, так как при Р, стремящемся к нулю, не дает уравнения Генри.
Однако, в отличие от уравнения БЭТ, основное уравнение ТОЗМ количественно описывает изотермы адсорбции паров на микропористых углях (в интервале 0.15< θ<1) и цеолитах (в интервале 0.4< θ<1) (см. таблицы
5.12-5.16).
Кроме того, в рамках ТОЗМ дополнительно показано, что: В случае микропористых активных углей можно рассчитывать:
Изотермы адсорбции большого числа адсорбтивов на данном адсорбенте по изотерме адсорбции одного (стандартного, как правило бензола) адсорбтива, поскольку коэффициенты подобия β с хорошей точностью могут быть рассчитаны по соотношениям:
β = |
Пi |
,βi E0i = βC6H6 E0,C6H6 |
(5.33), |
|
|||
|
ПС6H6 |
|
115
Где Пi, ПC6H6 – парахоры соответствующих адсорбтивов, βC6H6 = 1
При этом значения a0 c хорошей точностью находятся по соотношению
a0 |
= |
W0 |
(5.34), |
|
V0 |
||||
|
|
|
Поскольку отношения плотностей адсорбатов и адсорбтивов для больших групп веществ практически одинаковы и в (5.34) можно использовать табличные значения мольных объемов жидких адсорбтивов (V0). При этом,
значения W0 находятся с ошибкой, зависящей от соотношения мольных
объемов адсорбата и адсорбтива, но для последующих расчетов в рамках метода полного содержания эта ошибка не имеет значения, поскольку важен лишь факт постоянства W0.
Кроме того, на основании сравнения значений E0 с шириной микропор (x0), оцененной по данным углового рассеяния электронов, было установлено эмпирическое соотношение:
x0 = |
B |
= |
26.056 -1.53.10-5 E03.5 |
(5.35) |
βE0 |
βE0 |
Недавно было показано21, что в случае цеолитов коэффициенты подобия при использовании азота в качестве стандартного пара с хорошей точностью
могут быть рассчитаны по соотношению
β = 0,359 10 |
30 α |
|
EN |
2 |
|
П |
|
N2 |
|
|
|
(5.36), |
|||
|
|
|
|
||||
|
|
E ПN2 |
где: EN2 , E -потенциалы ионизации молекулы азота и адсорбтива (эВ), αN2 -
поляризуемость молекул азота.
Интересно отметить, что подобие характеристических кривых в случае адсорбции на цеолитах выполняется в существенно большем интервале величин адсорбции (0.2<θ<0.9), чем уравнение (5.32).
Отметим в заключение этого раздела, что в случае невыполнения постулата о температурной инвариантности характеристических кривых
116
температурная зависимость адсорбции на микропористых углях и цеолитах может быть рассчитана на основании установленной экспериментально в работах А. А. Фомкина и сотр. «линейности изостер адсорбции», выполняющейся в широких интервалах температур и, в том числе, при переходе через критические температуры адсорбтивов. Соответствующие линейные зависимости для любой постоянной величины адсорбции могут быть представлены в виде:
ln P = A − |
B |
(5.37), |
|
T |
|||
|
|
где А и В – константы, для определения которых необходимы экспериментальные данные как минимум при двух температурах.
В приведенных ниже таблицах 5.12-5.16 сопоставлены ошибки при описании экспериментальных данных по адсорбции различных адсорбтивов на микропористых активных углях и цеолитах, а в таблицах 5.17 -5.18- приведены параметры соответствующих уравнений17.
Таблица 5.12. Экспериментальная изотерма для системы пропан-АУ, 333К и относительные ошибки ( ,%) ее описания различными уравнениями.
p/ps |
а |
5.32 |
5.14 |
5.16а |
5.18а |
1,64*10-4 |
0,316 |
7,81 |
36,8 |
-32,7 |
16,5 |
2,74*10-3 |
1,222 |
-4,79 |
-6,98 |
2,01 |
1,17 |
4,83*10-3 |
1,560 |
-3,22 |
-6,93 |
3,09 |
0,24 |
1,39*10-2 |
2,320 |
-0,10 |
-3,37 |
1,74 |
-0,10 |
2,09*10-2 |
2,626 |
0,41 |
-2,19 |
0,74 |
-0,25 |
2,64*10-2 |
2,793 |
0,20 |
-1,98 |
0,07 |
-0,66 |
4,26*10-2 |
3,151 |
0,59 |
-0,71 |
-0,49 |
-0,30 |
5,16*10-2 |
3,270 |
0,10 |
-0,88 |
-0,85 |
-0,70 |
6,59*10-2 |
3,445 |
0,28 |
-0,33 |
-0,80 |
-0,33 |
9,69*10-2 |
3,690 |
0,16 |
0,04 |
-0,65 |
-0,03 |
1,07*10-1 |
3,756 |
0,22 |
0,20 |
-0,53 |
0,14 |
1,23*10-1 |
3,842 |
0,41 |
0,51 |
-0,28 |
0,48 |
0,31 |
4,325 |
1,47 |
1,83 |
1,32 |
2,03 |
|
|
|
|
|
|
0,40 |
4,348 |
-0,05 |
0,21 |
0,35 |
0,35 |
|
|
|
|
|
|
0,60 |
4,438 |
-0,06 |
-0,09 |
0,00 |
-0,36 |
|
|
|
|
|
|
0,70 |
4,459 |
-0,06 |
-0,24 |
-0,28 |
-0,78 |
|
|
|
|
|
|
0,85 |
4,476 |
-0,07 |
-0,44 |
-0,76 |
-1,46 |
|
|
|
|
|
|
117
0,90 |
4,478 |
-0,07 |
-0,50 |
-0,93 |
-1,69 |
|
|
|
|
|
|
Таблица 5.13. Экспериментальная изотерма для системы бензол-АУ, 303К и
относительные ошибки ( ,%) ее описания различными уравнениями.
p/ps |
а |
5.32 |
5.14 |
5.16а |
5.18а |
3,05*10-6 |
0,591 |
7,45 |
- |
-4,84 |
17,1 |
5,31*10-6 |
0,673 |
3,03 |
19,6 |
-3,41 |
13,2 |
8,00*10-6 |
0,792 |
6,62 |
17,9 |
2,30 |
16,2 |
1,54*10-5 |
0,909 |
1,54 |
7,51 |
2,19 |
11,0 |
2,16*10-5 |
0,996 |
1,22 |
4,98 |
3,32 |
10,3 |
2,75*10-5 |
1,054 |
0,26 |
2,78 |
3,48 |
9,06 |
1,35*10-4 |
1,454 |
-7,28 |
-8,15 |
1,19 |
-1,05 |
2,23*10-4 |
1,641 |
-6,31 |
-7,00 |
1,36 |
-1,20 |
2,28*10-4 |
3,562 |
-0,81 |
2,87 |
-1,23 |
-0,05 |
5,11*10-2 |
3,832 |
-0,48 |
2,61 |
-0,79 |
0,21 |
8,99*10-2 |
3,985 |
-0,46 |
1,83 |
-0,63 |
0,03 |
0,145 |
4,105 |
-0,14 |
1,24 |
-0,38 |
0,04 |
|
|
|
|
|
|
0,182 |
4,151 |
-0,10 |
0,76 |
-0,39 |
-0,14 |
|
|
|
|
|
|
0,231 |
4,196 |
0,02 |
0,28 |
-0,39 |
-0,30 |
|
|
|
|
|
|
0,276 |
4,226 |
0,11 |
-0,09 |
-0,42 |
-0,45 |
|
|
|
|
|
|
0,335 |
4,265 |
0,44 |
-0,31 |
-0,32 |
-0,42 |
|
|
|
|
|
|
0,432 |
4,323 |
1,15 |
-0,36 |
-0,05 |
-0,17 |
|
|
|
|
|
|
0,501 |
4,343 |
1,32 |
-0,67 |
-0,11 |
-0,32 |
|
|
|
|
|
|
0,704 |
4,420 |
2,60 |
-0,51 |
0,36 |
0,17 |
|
|
|
|
|
|
0,784 |
4,445 |
3,07 |
-0,41 |
0,53 |
0,36 |
|
|
|
|
|
|
0,865 |
4,471 |
3,58 |
-0,23 |
0,74 |
0,60 |
|
|
|
|
|
|
0,914 |
4,462 |
3,37 |
-0,65 |
0,47 |
0,23 |
|
|
|
|
|
|
118
Таблица 5.14. Экспериментальная изотерма для системы азот-цеолит NaX,
90К и относительные ошибки ( ,%) ее описания различными уравнениями.
p/ps |
а |
5.32 |
5.18а |
5.16а |
5.14 |
1,1*10-8 |
0,291 |
-245 |
-224 |
-231 |
-174 |
2,6*10-8 |
0,693 |
-113 |
-86,1 |
-71,7 |
-76,2 |
6,4*10-8 |
1,29 |
-63,6 |
-37,1 |
-19,5 |
-36,5 |
8,1*10-8 |
1,51 |
-52,0 |
-26,9 |
-10,9 |
-26,8 |
1,3*10-7 |
2,18 |
-22,8 |
-2,31 |
4,43 |
-2,15 |
1,4*10-7 |
2,33 |
-14,2 |
4,82 |
7,40 |
4,95 |
2,7*10-7 |
2,92 |
-14,5 |
3,73 |
8,33 |
5,56 |
2,9*10-6 |
5,51 |
-1,63 |
7,81 |
3,58 |
17,2 |
4,9*10-6 |
5,88 |
-3,54 |
4,40 |
1,59 |
15,3 |
1,7*10-5 |
6,72 |
-6,08 |
-1,31 |
-1,62 |
11,2 |
6,7*10-5 |
7,50 |
-6,97 |
-4,51 |
-3,07 |
6,69 |
1,3*10-4 |
7,87 |
-6,49 |
-4,80 |
-3,07 |
4,77 |
2,6*10-4 |
8,25 |
-5,05 |
-3,92 |
-2,44 |
3,63 |
5,8*10-4 |
8,63 |
-3,66 |
-3,05 |
-1,72 |
1,87 |
4,3*10-3 |
9,28 |
-0,71 |
-1,08 |
-0,32 |
-1,61 |
1,2*10-2 |
9,53 |
0,91 |
0,05 |
0,30 |
-1,92 |
2,0*10-2 |
9,62 |
1,56 |
0,48 |
0,49 |
-1,78 |
2,8*10-2 |
9,68 |
2,04 |
0,81 |
0,65 |
-1,54 |
5,7*10-2 |
9,80 |
3,05 |
1,53 |
0,97 |
-0,83 |
0,119 |
9,88 |
3,75 |
1,94 |
1,03 |
-0,21 |
|
|
|
|
|
|
119
Таблица 5.15. Экспериментальная изотерма для системы бензол-NaX, 493К
и относительные ошибки ( ,%) ее описания различными уравнениями.
p/ps |
а |
5.32 |
5.18а |
5.16а |
5.14 |
4,9*10-5 |
0,07 |
73,9 |
-17,9 |
-101 |
7,08 |
6,1*10-5 |
0,12 |
76,9 |
17,7 |
-34,3 |
33,1 |
1,3*10-4 |
0,18 |
49,9 |
2,36 |
-17,7 |
13,1 |
1,8*10-4 |
0,24 |
40,4 |
4,92 |
-4,76 |
12,3 |
2,7*10-4 |
0,32 |
25,9 |
1,88 |
2,45 |
6,20 |
3,3*10-4 |
0,32 |
9,51 |
-12,5 |
-1,88 |
-8,89 |
8,6*10-4 |
0,69 |
3,27 |
-0,48 |
6,46 |
-0,68 |
8,9*10-4 |
0,74 |
6,93 |
3,68 |
6,92 |
3,44 |
1,2*10-3 |
0,87 |
5,34 |
3,54 |
4,65 |
3,17 |
1,9*10-3 |
1,10 |
0,78 |
-0,45 |
0,01 |
-0,67 |
2,0*10-3 |
1,10 |
-0,52 |
-1,78 |
-0,34 |
-1,99 |
2,2*10-3 |
1,14 |
-0,39 |
-1,75 |
-0,77 |
-1,91 |
2,4*10-3 |
1,20 |
0,50 |
-0,99 |
-1,09 |
-1,11 |
2,6*10-3 |
1,26 |
1,36 |
-0,31 |
-1,27 |
-0,38 |
3,3*10-3 |
1,39 |
2,43 |
0,21 |
-1,30 |
0,22 |
3,4*10-3 |
1,40 |
2,63 |
0,38 |
-1,23 |
0,39 |
3,8*10-3 |
1,47 |
2,64 |
-0,02 |
-1,02 |
0,00 |
4,5*10-3 |
1,55 |
2,37 |
-0,82 |
-0,55 |
-0,83 |
6,4*10-3 |
1,74 |
3,70 |
-0,59 |
2,07 |
-0,81 |
120
Таблица 5.16. Экспериментальная изотерма для системы циклогексан-.СaY,
413К и относительные ошибки ( ,%) ее описания различными уравнениями.
.
p/ps |
а |
5.32 |
5.18а |
5.16а |
5.14 |
6,1*10-6 |
0,215 |
79,4 |
-3,30 |
-43,1 |
19,8 |
1,1*10-5 |
0,300 |
67,6 |
-1,32 |
-23,0 |
9,35 |
2,1*10-5 |
0,414 |
52,9 |
-0,68 |
-9,44 |
0,51 |
4,1*10-5 |
0,561 |
36,9 |
-1,31 |
-1,31 |
-5,73 |
5,6*10-5 |
0,664 |
32,2 |
1,75 |
2,76 |
-3,73 |
8,1*10-5 |
0,772 |
24,2 |
0,82 |
4,08 |
-5,22 |
1,6*10-4 |
1,01 |
12,8 |
0,13 |
4,54 |
-4,98 |
2,2*10-4 |
1,13 |
7,18 |
-1,40 |
3,50 |
-5,50 |
3,1*10-4 |
1,28 |
6,01 |
0,49 |
3,17 |
-2,47 |
4,1*10-4 |
1,43 |
5,22 |
1,90 |
2,52 |
-0,01 |
5,6*10-4 |
1,53 |
1,06 |
-0,64 |
0,79 |
-1,57 |
7,6*10-4 |
1,69 |
1,51 |
0,98 |
0,22 |
0,94 |
1,2*10-3 |
1,86 |
-0,66 |
-0,31 |
-1,14 |
0,56 |
1,5*01-3 |
1,96 |
-1,61 |
-1,00 |
-1,77 |
0,27 |
2,1*10-3 |
2,10 |
-1,32 |
-0,68 |
-1,90 |
0,85 |
3,2*10-3 |
2,28 |
-0,26 |
0,07 |
-1,50 |
1,59 |
6,7*10-3 |
2,54 |
1,41 |
0,63 |
-0,28 |
1,46 |
9,5*10-3 |
2,63 |
1,63 |
0,16 |
0,07 |
0,45 |
1,2*10-2 |
2,68 |
1,53 |
-0,53 |
0,14 |
-0,69 |
1,9*10-2 |
2,81 |
3,62 |
0,60 |
1,52 |
-0,27 |