книги / Современный катализ и химическая кинетика
..pdf2.3. Реакции и термодинамическое равновесие |
51 |
2.3.2.Химическое равновесие для неидеальных газов
Впроведенном ранее примере мы считали газы идеальными, что обычно является достаточно хорошим приближением. Однако при высоких давлениях и низких температурах взаимодействие между молекулами может привести к заметным отклонениям от идеальности. В этих случаях мы можем отождествлять активности с парциальными давлениями или концентрациями. Для учета этих отклонений вводят [3] фугитивности (летучести) fx, связанные с парциальными давлениями соотношением:
fx = ФхРх-> |
(2.39) |
где фх — коэффициенты фугитивности. Эти коэффициенты безразмерны, по этому фугитивности fx имеют размерность давления. Коэффициенты фугитив ности могут быть рассчитаны через вириальные коэффициенты и определены по табличным данным для относительных (к критическим значениям) давле ния и температуры газа. Коэффициент фугитивности может иметь значения как больше, так и меньше единицы в зависимости от данных значений давле ния и температуры. Идеальным газам отвечает значение = 1. Ниже мы крат ко изложим способ расчета равновесной степени превращения с учетом отлич ных от единицы коэффициентов фугитивности.
Как уже упоминалось, активности компонентов совпадают с их фугитивно стями (а не с парциальными давлениями), поэтому корректная форма уравне ния (2.22) имеет вид
zz /^ \ ^C^P* |
fc fo ( o\•a + b-c-d ■ |
<2-40) |
a \ “B |
J A J В |
|
Для идеальной смеси газов (которые сами себя ведут неидеально) фугитив ность *У компонента X дается выражением
fx = ФхРх -^xtot^xAot5 |
(2.41) |
где (*$< — коэффициент фугитивности чистого газа X при давлении ptot и дан ной температуре. Заметим, что коэффициент фугитивности зависит как от дав ления, так и от температуры. С учетом определения (2.41) получаем
*(П = |
Ус |
Ур |
A Ptot Ж Plot |
Y V |
( п \ vc + vd “ va “ vd |
rp |
7C ZD |
Ftot |
|||
Cva |
Ув |
ф^ф^ те I / J |
|||
|
J\ |
|
|
|
(2.42)
Y vc у vd
7C 7P
K(T)=
Y Va Y vb
7a 7b
Применяя это выражение к случаю синтеза аммиака, находим
*(П = |
°NH3 |
/nh3 |
А |
(2.43) |
|
n2 н2 |
f1/2 |
f3/2 |
|||
|
У n2 |
Ун2 |
|
|
52 —Глава 2- Кинетика химических реакций
При тех же предположениях, что были использованы ранее (постоянное давление и стехиометрическая смесь N2 + ЗН2), получаем
Обратите внимание, что поправка свелась к ведению правильной констан
ты равновесия К (Т)/К?™. При низких давлениях наличие этой поправки прак
тически не сказывается на величине мольной доли аммиака, однако при давле нии ptot =100 бар и выше поправка начинает проявлять себя. Результаты расче тов представлены в последнем столбце табл. 2.4 и на рис. 2.1. Следует заметить, что существуют процедуры учета поправок и в том случае, когда смесь не явля ется идеальной, однако их рассмотрение выходит за пределы данной книги.
2.4.ЗАВИСИМОСТЬ СКОРОСТИ РЕАКЦИЙ ОТ ТЕМПЕРАТУРЫ
Сповышением температуры скорость реакции обычно увеличива ется, хотя для каталитических реакций, как будет показано далее в этой главе, такое увеличение наблюдается не всегда. Существует (очень приближенное) правило, согласно которому скорость реакции удваивается при увеличении тем пературы на каждые 10 К.
2.4. Зависимость скорости реакций от температуры |
53 |
Температурная зависимость константы скорости элементарных реакций хорошо описывается уравнением Аррениуса:
k(T) = ve~EJRT, |
(2.45) |
где итак называемый предэкспоненциальный множитель и £а — энергия акти вации (кДж/моль). Можно выражать энергию активации и в джоулях на одну молекулу (очень малая величина), тогда газовую постоянную R в формуле (2.45) необходимо заменить на постоянную Больцмана кв.
Аррениус получил уравнение (2.45), основываясь на имеющихся эксперимен тальных данных и результатах собственных исследований по гидролизу сахарозы, фруктозы и глюкозы. Заметим, что предложенная Аррениусом температурная зави симость является чисто экспоненциальной и в ней предэкспоненциальный мно житель является постоянным. Теории скоростей реакций, обсуждаемые в гл. 3, показывают, что аррениусовская зависимость с высокой степенью точности яв ляется корректной. Однако постоянство (независимость от температуры) предэкспоненциального множителя не подтверждено ни одной из существующих теорий. В теории переходных состояний предэкспоненциальный множитель имеет
степенную зависимость от температуры: |
Тх. Но эта зависимость незначи |
тельна в сравнении с экспоненциальной, поэтому обычно ею пренебрегают.
Рис. 2.2. Схематическое представление профи
ля потенциальной энергии для рекомбинации |
Координата реакции |
двух атомов или диссоциации молекулы |
или расстояние X—X |
Физическая интерпретация уравнения Аррениуса основывается на предпо ложении о существовании энергетического барьера, с преодолением которого связано протекание полной реакции. На рис. 2.2 эта ситуация проиллюстриро вана для реакции образования молекулы из двух атомов. Когда происходит сближение атомов по правой ветви потенциальной кривой, то сначала проис ходит отталкивание между ними, связанное с ростом потенциальной энергии, идущим до достижения расстояния, на котором собственно и происходит реак ция, сопровождающаяся, например, перестройкой электронной конфигурации и образованием связи. На этом расстоянии потенциальная энергия достигает максимума. Формирование связи вызывает уменьшение потенциальной энер гии при дальнейшем сближении атомов. Сейчас мы отождествим высоту по-
54 |
Глава 2. Кинетика химических реакций |
тенциального барьера с энергией активации (ниже мы увидим, что это отожде ствление не вполне корректно). Отметим, что на рис. 2.2 схематически показа на также энтальпия реакции.
Энергию активации можно определить по значениям скорости реакции, измеренным при разных температурах. Зависимость логарифма константы от обратной температуры, как видно из формулы (2.45), имеет вид уравнения пря мой линии:
lnfc(T) = lnv-4i- |
(2-46) |
К 1
Угол наклона этой прямой определяется энергией активации, а точка пере сечения с осью ординат — пердэкспоненциальным множителем. Нахождение энергии активации, таким образом, сводится к элегантной процедуре, которая постоянно будет использоваться ниже и выражается формулой:
Еа |
= -R ^4 = RT2 |
. |
(2.47) |
а |
Э(1/Т) |
эт |
v ' |
Читатель может теперь легко рассчитать величину энергии активации, от вечающую правилу: скорость реакции удваивается при увеличении температу ры на каждые 10 К.
Как уже отмечалось выше, химическая реакция достаточно редко протека ет в одну стадию, так что нам необходимо приспособить данное нами опреде ление энергии активации к многоступенчатым реакциям. Поскольку нас не будут интересовать термодинамические эффекты, мы определим кажущуюся энергию активации через соотношение
(2.48)
Это определение является следствием построения экспериментальных дан ных для скорости реакции в аррениусовских координатах lnr+ — 1/Т. Однако в общем случае экспериментальные данные могут и не укладываться на линейную зависимость. Таким образом, кажущаяся энергия активации пригодна для опи сания химических реакций только в ограниченном интервале температур. При интерпретации энергии активации нужно быть особенно внимательным в отно шении порядка реакции, поскольку он обычно зависит от условий проведения эксперимента. Ниже приведен пример, в котором скорость прямой реакции за висит как от активированного процесса, так и от равновесных стадий. Этот при мер отражает ситуацию, часто встречающуюся в каталитических реакциях.
Пример
Предположим, что мы имеем дело с многоступенчатой реакцией, скорость пря мого процесса в которой характеризуется следующей зависимостью от давления
+кК{К2р2
Г |
1 + К2р ‘ |
(2.49) |
2.5. Решения уравнений скорости реакций: зависимости от времени концентраций —J |
55 |
Мы примем, что параметры к и Кх могут быть представлены в виде |
|
k(T) = veE*> RT |
(2.50) |
И |
|
кх (Г) = e~^IRT = e-^x/KTe^x/R = Vxe-^HX/RT |
(2,51) |
Непосредственное применение уравнения (2.48) дает |
|
= Еа + ДЯ, + ДЯ2 - -5^2- щ. |
(2.52) |
Это соотношение показывает, что кажущаяся энергия активации представ ляет собой истинную энергию активации, «дополненную» энтальпиями равно весных стадий, то есть в данном случае кажущаяся энергия активации зависит как от давления, так и от температуры системы. Заметим также, что здесь мы пренебрегли возможной неэкспоненциальной зависимостью параметров. Как будет показано в гл. 3, величины v, ДЯи Д5 в определенной степени являются функциями температуры.
2.5.РЕШЕНИЯ УРАВНЕНИЙ СКОРОСТИ РЕАКЦИЙ: ЗАВИСИМОСТИ ОТ ВРЕМЕНИ КОНЦЕНТРАЦИЙ
ВРЕАКЦИЯХ РАЗЛИЧНОГО ПОРЯДКА
Вобщем случае скорость реакции определяется по изменению со става системы за определенный временной интервал. Для ее нахождения мы должны установить зависимости от времени концентраций компонентов, осно вываясь на приведенных в предыдущих разделах дифференциальных уравнени ях, связывающих скорость реакции с концентрациями реагирующих молекул.
Рассмотрим реакцию первого порядка (очевидно, реакцию изомеризации, хотя диссоциация соединения R два фрагмента Pj и Р2, описывается точно так же)
R-^P, |
(2.53) |
скорость которой, записанная в дифференциальной форме, равна
'• = -^ = ВД- |
(2.54) |
Для получения зависимости от времени концентраций компонентов R и Р необходимо проинтегрировать это уравнение, что легко делается методом раз деления переменных:
/'7^ = /(^)d/ = ln^ = ’^ |
(2’55) |
|R]0 LRJ о |
URJo J |
56 |
Глава 2. Кинетика химических реакций |
Пределы интегрирования были выбраны в соответствии с начальным усло вием [R] = [R]o при t = 0. Из решения (2.55) легко находятся зависимости от времени концентрации реагента и продукта реакции:
[R] = [R]oe^; [Р] = [R]o (1 - |
(2.56) |
Выражение (2.56) показывает, что, построив зависимость In [R] от времени, можно легко установить, имеет ли данная реакция первый порядок. Для реакций первого порядка зависимость In [R] от времени должна быть линейной, а тан генс угла наклона соответствующей прямой совпадает с константой скорости.
Если реакция первого порядка действительно является необратимой, то по времени Г1/2, необходимому для превращения 50 % продукта в реагент, которое для всех экспоненциальных зависимостей называется периодом полураспада или полупревращения, также можно определить константу скорости:
к = —. |
(2.57) |
Л/2
Интегрирование кинетических уравнений для реакций второго порядка про водится чуть сложнее. Рассмотрим реакцию второго порядка следующего вида
А + В—>АВ, |
(2.58) |
для которой уравнение скорости таково
= d[AB] = _d[A] = _d[B] = А_
dz |
dz |
dz |
i л j |
Примем, что за некоторое время t прореагировало некоторое количество х реагента А, тогда концентрации реагентов А и В станут равными [А]о — х и [В]о — х соответственно. Используя данные обозначения, можно после разделения пе ременных переписать уравнение (2.59) следующим образом
(2.60)
Ио -х)([В]0 -х)
Стандартный путь интегрирования таких дифференциальных уравнений состоит в разбиении левой части на два слагаемых. То есть нам необходимо найти коэффициенты а} и я2, обеспечивающие выполнение равенства
#1 |
^7 |
_ |
1 |
Г? 611 |
([А|0-х)([В|„ -х)'
Эти коэффициенты легко находятся с помощью алгебраических преобразова ний, и мы имеем
f_J_____ L_ |
1 |
dx = | kdt. |
(2.62) |
|
[В]о - * |
||||
|
0 |
|
А 57
2.5. Решения уравнений скорости реакций: зависимости от времени концентраций
Данное уравнение легко интегрируется, и, с учетом начального условия х = О при t = 0, получаем
1 |
1п гьч *[ч- |
(2.63) |
[В]0-[А]0 |
[[В][А]0-х. |
Представляет интерес рассмотреть два частных случая. Пусть А и В являют ся одним и тем же компонентом, например, это могут быть радикалы О, объе диняющиеся в молекулу О2. В этом случае уравнение (2.63) неприменимо, по скольку [А]о = [В]0, и имеет место деление на ноль. Это говорит о необходимо сти вернуться к уравнению (2.60), которое в данном случае принимает вид
(2.64)
*([А]) .- 2
и интегрируется непосредственно:
1_______1
(2.65)
[Л]0-х [А]о
Таким образом, если мы построим зависимость величины 1/([А]ох) от вре мени, то получим прямую линию, тангенс угла наклона которой равен кон станте скорости.
Вторая интересная ситуация возникает, когда один из компонентов присутству ет в значительном избытке, например, пусть [В] » [А]. Тогда [В] = [В]о — х = [В]0, и уравнение (2.63) принимает простую форму, отвечающую реакции псевдо первого порядка
[А] |
= |В],й=ЧЛ] = (Л11е-1"1< |
(2.66) |
|
[А]о - *
Для нахождения константы скорости в данном случае нужно использовать метод, описанный при анализе реакций первого порядка, то есть необходимо построить зависимость In [А] от времени. При этом произведение ЦВ]0 будет выступать в качестве константы скорости реакции псевдопервого порядка, из которой легко находится константа скорости реакции второго порядка.
Наконец, мы должны рассмотреть реакции нулевого порядка, хотя они встре чаются достаточно редко. Порядок реакции становится псевдонулевым, если реагент присутствует в значительном избытке и его концентрация практически не изменяется в ходе реакции. Дифференциальная и интегральная формы урав
нения скорости реакции нулевого порядка R |
Р таковы |
|
||
d[P]_ |
d[*] |
[R] = [R]o - kt; |
(2.66) |
|
dt |
= k |
|||
dt |
|
|
|
На рис. 2.3 для сравнения показаны зависимости от времени степени пре вращения реагентов для реакций разного порядка. В табл. 2.5 сведены диффе-
58 —Глава 2' Кинетика химических реакций
ренциальные и интегральные формы уравнений скорости для реакций разного порядка.
|
|
|
Рис. 2.3. Зависимости концен |
|||
|
|
|
трации |
продукта от времени |
||
|
|
|
для реакций разного порядка |
|||
Таблица 2.5. Уравнения скорости реакции в дифференциальной |
|
|||||
|
и интегральной формах |
|
|
|
|
|
Порядок |
Дифференциальная форма |
Интегральная форма |
Размерность |
|||
d[A]_ |
||||||
реакции |
kt = |
k |
|
|||
dr |
|
|||||
|
|
|
|
|
||
0 |
-k |
[A]o |
- [A] |
моль • л-1 |
• c-1 |
|
1 |
~k [A] |
|
[A] |
c-‘ |
|
|
|
|
|
|
|
||
2 |
-£[A]2 |
1 |
1 |
моль • л-1 |
• c-1 |
|
[A] ’ [A]o |
||||||
|
|
|
|
|||
2 |
-MA][B] |
i |
ln [a]q [b] |
моль • л-1 |
• c-1 |
|
[B]0-[A]0 |
[A]0[A] |
|||||
|
|
|
|
2.6.ВЗАИМОСВЯЗАННЫЕ РЕАКЦИИ В ПРОТОЧНЫХ РЕАКТОРАХ:
ПРИБЛИЖЕНИЕ СТАЦИОНАРНОГО СОСТОЯНИЯ
Как уже отмечалось выше, практически все представляющие прак тический интерес реакции являются многоступенчатыми. Возникает вопрос, как кинетика полной реакции складывается из элементарных ступеней. Рас смотрим двухступенчатую реакцию, записанную в общем виде
(2.68)
i |
i |
i |
2.6. Взаимосвязанные реакции в проточных реакторах Л 59
Анализ кинетики таких реакций следует проводить с учетом типа использу емого реактора. Мы рассмотрим реакторы двух типов — трубчатый (колонча тый) и емкостной (рис. 2.4).
Рис. 2.4. Схемы трубчатого проточного и емкостно |
|
Резервуар |
го реакторов. В идеале трубчатый реактор работает, |
Труба |
|
как реактор вытеснительного типа, в котором газ тол |
|
|
|
|
|
кается поршнем, движущимся по трубе. В идеаль |
|
|
ном емкостном реакторе компоненты в резервуаре |
|
|
постоянно перемешиваются, так что отсутствуют гра |
|
|
диенты концентраций компонентов. При открытых |
|
|
вентилях, показанных на рисунке, идет постоянная |
Трубчатый реактор |
Емкостной |
подача реагентов в резервуар и выгрузка продукта. В |
||
этом случае мы имеем дело с емкостным реактором |
вытеснительного типа |
реактор |
с постоянным перемешиванием |
|
|
Трубчатые реакторы обычно используются в крупномасштабных технологи ческих процессах. Реагенты подаются в него с постоянной скоростью, а продукт появляется на выходе также с постоянной скоростью. Говорят, что такие реакто ры работают в стационарных условиях, подразумевая, что как скорости реакции, так и концентрации компонентов в данной зоне реактора не зависят от времени (имеются, конечно, малые флуктуации параметров около средних значений).
Емкостные реакторы используются для получения продуктов тонкого син теза. Реагенты загружаются в них и постепенно превращаются в продукт. Оче видно, что в этом случае и скорость реакции, и концентрации компонентов являются функциями времени. Мы начнем наш анализ с кинетики взаимосвя занных реакций, протекающих в стационарном режиме.
Без ограничения общности анализа мы можем упростить его, заменив реа генты, интермедиаты и продукты одиночными молекулами:
R < |
> I |
> Р. |
(2.69) |
|
к\ |
|
|
Уравнения скоростей реакций для каждого изучастников процесса имеют вид:
^ = -^[^ +^[1], |
(2.70) |
= |
(2.71) |
= |
(2.72) |
Предположив, что процесс идет в стационарных условиях, как это имеет место в промышленных ректорах при постоянных скоростях подачи реагентов и отвода продуктов, мы можем концентрацию интермедиата считать постоян ной, что соответствует стационарному приближению:
dFll Л |
(2.73) |
-Н = 0- |
dt
60 |
Глава 2. Кинетика химических реакций |
Это равенство позволяет исключить из уравнений скоростей реакций кон центрацию компонента I, которая в общем случае практически не поддается измерению:
= |
<2-74) |
AVj |
I |
Скорость образования продукта теперь определяется выражением
d[P] |
(2.75) |
d/ k^ + k^
Возможны два предельных случая. Если скорость второй стадии реакции мала, мы можем пренебречь соответствующей константой в знаменателе фор мулы (2.75):
dlP] = ^^[Rl = ^i[R] ^<<k-y |
(2.76) |
В этом случае скорость полной реакции лимитируется стадией превраще ния интермедиата в продукт, которая и является определяющей. Первой ступе ни при этом отвечает виртуальное равновесие, так что концентрация интерме диата зависит только от константы равновесия Кг
Кажущаяся энергия активации находится путем вычисления логарифми ческой производной, как это выражает формула (2.48):
^каж = |
\Т 4 = Еа. 2 + ДЯ,. |
(2.77) |
|
о 1 |
|
Энергия активации для реакции как целого обычно называется кажущейся энергией активации, и в данном случае является комбинацией параметров двух ступеней реакции. Если энтальпия первой стадии отрицательна и велика по абсолютной величине (то есть на первой стадии реакция экзотермична), то кажущаяся энергия активации может стать отрицательной. Это означает, что с увеличением температуры концентрация интермедиата падает быстрее, чем растет константа скорости второй ступени.
Второму предельному случаю отвечает быстрая вторая стадия, когда опре деляющей для скорости реакции является первая стадия. В этом случае уравне ние скорости принимает вид
M = ^+[R] |
(2.78) |
Энергия активации реакции как целого совпадает с энергией активации первой стадии Еа г
Заметим, что быстрая стадия, следующая за лимитирующей, является ки нетически несущественной. Наоборот, если быстрая стадия предшествует ли