Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Современный катализ и химическая кинетика

..pdf
Скачиваний:
46
Добавлен:
12.11.2023
Размер:
12.26 Mб
Скачать

2.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)

Энергия активации реакции как целого совпадает с энергией активации первой стадии Еа г

Заметим, что быстрая стадия, следующая за лимитирующей, является ки­ нетически несущественной. Наоборот, если быстрая стадия предшествует ли­