Популярное

Мифы о звукоизоляции



Как построить дом из пеноблоков



Как построить лестницы на садовом участке



Подбираем краску для ремонта



Каркасные дома из дерева


Главная » Об индукционном периоде

ОБ ИНДУКЦИОННОМ ПЕРИОДЕ РЕАКЦИЙ ЖИДКОФАЗНОГО ОКИСЛЕНИЯ УГЛЕВОДОРОДОВ

Галина Г.К., Спивак С.И. (spivak@bsu.bashedu.ru), Комиссаров В.Д. Башкирский государственный университет

Характерной чертой реакций цепного окисления углеводородов является автоускорение реакции в начальный период, обусловленное накоплением гидроперекисей. Классический график накопления гидроперекиси, приведенный в известной монографии [1], имеет s-образный характер. Начальный период протекания реакции характеризуется очень медленным накоплением гидроперекиси. Кривая накопления достаточно длительный период времени параллельна оси абсцисс. Этот временной интервал медленного накопления гидроперекиси определяется как индукционный период.

Задачей настоящей работы является математическое моделирование стадии индукционного периода:

- вывод соотношений на кинетические константы, гарантирующих индукционный период;

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

Механизм цепного окисления углеводородов в жидкой фазе может быть представлен следующей схемой [1 ]:

(0) RH+O2 - R+HO2;

(1) R+O2 - RO2;

(2) RO2+RH - ROOH+R;

(3) ROOH - RO+HO;

(31) RO+RH - ROH+R;

(32) HO+RH - H2O+R;

(4) R+R - P1;

(5) RO2+R - P2;

(6) RO2+RO2 - P3.

(Здесь P1 , P2 , P3 - молекулярные продукты.) В настоящей работе приняты следующие единицы измерений: концентрации - в моль/л, константы скоростей реакций первого порядка -в 1/с, константы скоростей реакций второго порядка - в л/моль с. Нумерация реакций выполнена в соответствии с [1 ].

Свободные радикалы R и RO2 весьма активны и быстро исчезают по реакциям (4) -(6). Поэтому в реакции окисления очень быстро устанавливается стационарная (кинетически-равновесная) концентрация этих радикалов, при которой скорости образования и гибели радикалов равны. Это позволяет в соответствии с методом стационарных концентраций Бой С d [ R ] d[ RO2 ]

денштейна-Семенова положить равными нулю производные - и - и получить

dt dt

выражение для скорости реакции окисления.

Выведем формулу для скорости цепного окисления в случае достаточно большой концентрации растворенного кислорода, когда [RO2] >> [R] и обрыв цепей происходит только по реакции (6). В этом случае d [ R ]

Wi -kj[O2][ R] + k2[RH][ RO/] = 0 (7)

d[R°2l= kj[O2][ R] -k2[RH][ RO2]-2k6[ RO2]2 = 0 , (8)



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

W =2k6[ RO2]2

[ RO2T2

2k 6

Скорость окисления углеводорода W=k2[RH][RO2] = k2[RH]

(10)

(11)

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

Длиной цепи v называется число полных звеньев цепи, приходящихся в среднем на каждый свободный радикал, образовавшийся в системе в результате реакций зарождения и вырожденного разветвления цепей.

Длина цепи равна отношению скоростей продолжения и обрыва или, что то же, инициирования цепей

1 k2

W V2k6

- [ RH ] -

(12)

(13)

- l j /-

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

Рассмотрим кинетику накопления гидроперекиси в начальный период окисления при достаточно большой концентрации растворенного кислорода, когда [RO2*] >> [R*] и длинных цепях (v>10). Расходованием гидроперекиси в начальный период можно пренебречь. Пусть механизм окисления включает в себя элементарные реакции (0), (1), (2), (3), (6).

Скорость образования гидроперекиси в такой реакции равна

d [ ROOH ]

k 2[ RH ]

W0 + k 3 [ ROOH ]

( 1 4)

Введем безразмерные переменные:

Тогда

dn dx

[ROOH ]

[RH]

dn dt dt dx

x =-k 2

22 где

- t , CO 0--

2k6 k3[ RH ]

(15)

t =-

dn dt

dt dx

[RH]

d dt

[ ROOH ] . [RH] .

k2 \

k2 \

k 3[ RH ]

x ;

d[ROOH] 2

k3[RH]

k2[RH]

W0 + k3[ ROOH]

[RH]

k2 К

k3[RH]

k3[RH]

= 2

W0 + k3[ROOH]

k 3[ RH ]

= 2

Wp + [ ROOH ]

k3[RH] [RH]

co0 +n ;



таким образом, . П^~+п , (16)

откуда в результате интегрирования получаем

т2 + 2тл/о>0 . (17)

При очень малых скоростях начального инициирования цепей (W0 очень мало): ~ т2 и [ROOH] ~ t2 , (18)

то есть накопление гидроперекиси должно протекать по очень простому математическому закону:

J[ROOH] k2L [RH].t. (19)

Таким образом, с математической точки зрения, индукционный период определяется как такой промежуток времени [2], на котором с некоторой точностью 8 выполняется соотношение (1 9):

irrooh]-klJE-[RH] t \< 8 (20)

2V k 6

В работе [2] индукционный период изучается на модельном механизме цепной реакции окисления углеводородов

RH + O2 - ROOH (21)

в жидкой фазе, который может быть представлен, например, следующей совокупностью элементарных стадий:

(a) RH-R + H

(b) R + RH + O2 -ROOH + R

(c ) ROOH + RH -R + ROH + OH (d) R + R -R-R ,

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

Дальнейший анализ проводится в предположении

k- £-а , kl £-в , JkL £-2r (а, в, у> 0) (22)

k1 k1 k1

где £ << 1. Это предположение физически соответствует факту наибольшего кинетического сопротивления стадии зарождения цепи.

Далее ставится задача об отыскании условий, при которых в пределе при £ - 0 концентрация продукта реакции (гидроперекиси ROOH) в зависимости от т = £8. t медленно изменяется на некотором начальном участке [0; Т]. Иначе говоря,

lim[ROOH ] ct2 , (23)

£-0



требуется найти с и 8 .

Предполагается также выполнение следующего условия

в < а < 2у

(24)

Физически условие (24) означает, что взаимодействие радикалов происходит наиболее легко, взаимодействие радикала с молекулой происходит труднее, а распад молекулы на радикалы -еще труднее.

Схеме реакций (a), (b), (c), (d) с учетом предположений (22) соответствует следующая система обыкновенных дифференциальных уравнений:

d [ RH ] dt

Wj - w2

d [ R ]

wj + w2 - 2w4

(25)

d [ ROOH ] dt

= w2 - w3

где [RH] - концентрация углеводорода, [R] - концентрация радикала и [ROOH] - концентрация продукта реакции (гидроперекиси). Здесь

wj = kj [RH]

w2 = е-аk2 [RH]-[ R]

w3 = e-p k3 [RH]-[ ROOH]

(26)

w4 = £

-27,

k4 [R]2

представляют собой скорости соответствующих элементарных стадий. Начальные условия :

[RH](t=0) = 1,

[ROOH](t=0) = [R](t=0) = 0.

(27)

В работе [2] доказано, что необходимым условием существования индукционного периода являются следующие условия на порядки констант скоростей элементарных стадий:

а > 7 , в 2(а - 7)

(28)

При выполнении условий (28) получаем формулу (29) для индукционного периода системы (25)

[ROOH]

k 3 k 2

т2

(29)



t = £

a-у+0,5/ т

(30)

Выполнение формулы (29) можно считать естественным определением индукционного периода. То есть можно определить индукционный период как такой промежуток времени [ T1; T2], на котором верно соотношение (29).

Получим формулу (29) из условий (28).

При t< T2 можно пренебречь изменениями [RH] и перейти к упрощенной системе дифференциальных уравнений для [ROOH] и [R*]:

d [ ROOH ] dt

= £-ak2[R] - £-ek3[ROOH]

d[R] = kj + £-%[ROOH] - £-2rk4[R]2 dt

(31)

Система (31) имеет единственное устойчивое положение равновесия: k3[ROOH] = k2£e - 2(a-Y)- z

[R] = £2у-a z , (32)

где z - положительный корень уравнения

k4 Z2- k2 Z - kj£2(a-Y) = 0

(33)

Решая (33) при условии a > у (первое из условий (28)), мы получаем из (32) условие накопления гидроперекиси [ROOH]: / < 2(a - у ) (второе из условий (28)). Выбирая масштабы измерения [ROOH] и [R*] в соответствии с (32)

[ROOH] = £

в -2(a- у)

У

[R] = £2у- a z получаем вместо (31 ) систему

(34)

(35)

dz dt

k2 Z - k3 y

= k3 y - k4 Z2 - kj£ 2(a-

(36)

При £ -> 0 (учитывая a > / (24) и a > у (28)) получаем систему

£ - = k2 z - k3 y

dt 2 3 y

k3 y = k4 z2

(37)

£



Отсюда следует формула (29).

При а < 7 аналогичная схема рассуждений приводит к линейной зависимости [ROOH] от t , то есть не происходит медленного накопления [ROOH] (а это значит, что индукционного периода при а < 7 нет).

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

Естественно для проведения вычислительного эксперимента выбрать такую реакцию, для которой достаточно надежно, в независимых условиях, определены константы скоростей элементарных стадий. Такой реакцией, например, является цепное окисление циклогексана в жидкой фазе. Экспериментально реакция исследовалась различными авторами в [1; 5-12]. Анализ этих работ позволил выделить следующие значения констант скоростей химических реакций, приведенные в таблице 1 .

Таблица 1 .

600 С

900 C

510-17

7-10-15

2,747106

3,5-106

0,13

0,72

5,6-10-10

4-10-8

5-106

6-106

2-106

4,4-109

6,6-109

2-109

2,7-109

2,8-106

3,33-106

Проверим выполнение предположений (22),(24) и необходимого условия существования индукционного периода (28) при £ = 10 .

При температуре 600С в случае £ = 10-20 выполнены предположения (22) и (24):

1) k2 /ko - £-а, k3/ko - £-Р , Wko - £-27 (а, р, 7 > 0);

2) р<а<2у

Действительно, k2 /ko - £-15/20, k3/ko - £-7/20 , k(/ko - £-22/20. Отсюда а = 3/4, р = 7/20, 7 = 11/20, то есть являются положительными и 7/20 < 3/4 < 22/20. Выполнено также необходимое условие существования индукционного периода (28): а > 7, Р < 2(а - 7). В самом деле, 3/4>11/20 и 7/20 < 2(3/4 - 11/20) = 8/20 = 0,4.

При температуре 900С в случае £ = 1 0-20 также выполнены предположения (22) и (24):

1) k2 /ko - £-7/10, k3/ko - £-3/10 , k(/k0 - £-1, т.е. а = 0,7, р = 0,3, 7 =0,5 (27 = 1);

2) 0,3 < 0,7 < 1 ;

выполнено и необходимое условие существования индукционного периода(28): а > 7 (0,7 > 0,5) и р <2(а - 7) (0,3 < 2(0,7 - 0,5) = 0,4).

Таким образом, имеем теоретическое подтверждение существования индукционного периода для реакций накопления гидроперекиси.

В соответствии со схемой реакций (0)-(6) изменение концентраций [RH], [R], [RO2] и [ROOH] во времени описывается следующей системой дифференциальных уравнений: dxo/dt= -wo -w2 -w3j -w32: dxj/dt= wo -wj +w2 + w3j + w32 -2w4 -w5,

dx2/dt= wj - w2 -w5 -2w 6, (38)

dx3/dt= w2 - w3



dx3j/dt= w3 - w3j dx32/dt= w3 - w32

с начальными концентрациями: x0 (0)=0.1; x1(0)= 10- ; x2(0)=0; x3(0)=0; x31(0)=0; x32(0)=0; [O2] = 10-3.

Здесь обозначены концентрации веществ x0=[RH], x1=[R], x2=[RO2], x3=[ROOH], x31=[RO], x32=[HO] и скорости элементарных стадий w0=k0x0[O2], wkxO, w2=k2x2x0,

w3=k3x3 w31= k31x31x0, w32=k32x32x0, w4=k4x12 w5=k5x1x2, w6= k6 x22.

В работе [3] был предложен метод редукции систем дифференциальных уравнений химической кинетики, который состоит в расчете в каждый момент времени значений всех скоростей реакций wi и отношений модулей концентраций ко времени xi(t)/t (z=1,...,w). Выбирается некоторое значение £, и реакция считается значимой в данный момент времени, если wi > £ . В противном случае соответствующая реакция считается незначимой, ее скорость в правых частях системы дифференциальных уравнений зануляется на промежутке времени, на котором выполняется условие wi < £ .

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

На основании программы работы [4] произведены расчеты длины индукционного периода. В таблице 2 представлены результаты этих расчетов.

Таблица 2.

t0 = 600 C

t0 = 900 C

5-10-4

1739 час. 72,5 сут.

13,3 мин.

6-10-4

2474 час. 103 сут.

13,94 час

7-10-4

3211 час. 134 сут.

31,2 час

810-4

3950 час. ~ 165 сут.

48, 5 час.

9-10-4

4691 час. 195 сут.

65,9 час.

10-3

5435 час. 226 сут

83,3 час

Здесь 8 есть точность, с которой выполняется формула (19) (8 =yj[ROOH ] -[rh ] t).

М k 6

Графики зависимости длины индукционного периода нд от точности 8 при температуре 600 C представлены на рис. 1, при температуре 900 C - на рис. 2.

Анализ рис. 1 и 2 показывает, что длины индукционных периодов существенно зависят от точности 8 , с которой выполняется условие (20), являющееся математическим определением индукционного периода. С ростом 8 длина индукционного периода увеличивается. Это естественно, ибо рост 8 означает уменьшение точности определения индукционного периода.

Длина индукционного периода, как это видно из таблицы 2, существенно зависит от температуры. С ростом температуры развитие реакции происходит более активно, что приводит к уменьшению индукционного периода.

Существенное значение для расчетов имеет величина порога значимости £ . Он может быть выбран из разных соображений, В наших расчетах £ был выбран из условия, что с нулевого момента времени должна быть значимой реакция (0) в исходной системе (0) - (6). В противном случае весь сложный механизм (0) - (6) не может быть реализован. Значимость реакции (0) означает выполнение следующего условия для £ :

w0 > £ (39)



Из(39)следует

£ < 10-18 . (40)

В наших расчетах было зафиксировано £ = 10- .

Таким образом, использование математического определения индукционного периода как удовлетворение соотношения (20) в совокупности с системой условий (22), (24), (28) позволяет рассчитать длины индукционных периодов в зависимости от характеристик реакции. Рассчитанные значения могут иметь важное значение для определения оптимальных технологических режимов соответствующих процессов в промышленности, что будет предметом дальнейших исследований.

120000

80000

40000 Н

2.2E-4

2.3E-4

2.4E-4 2.5E-4

Рис. 1.

2.6E-4

2.7E-4 5



10000 -л

8000 -

6000 -

4000 -

2000 -

4.6E-4 4.8E-4 5.0E-4 5.2E-4 5.4E-4 8




ЛИТЕРАТУРА

1. Эммануэль Н.М., Денисов Е.Т., Майзус З.К. Цепные реакции окисления углеводородов в жидкой фазе. М.:Наука, 1965. 375 с.

2. Спивак С.И., Шабат А.Б., Шмелев А.С. В сб.: Нестационарные процессы в катализе. Материалы всесоюзной конференции. Часть 1. Новосибирск. 1979. С.118-121.

3. Тропин А.В., Масленников С.И., Спивак С.И. Кинетика и катализ. 1995. Т.36. №5.

С.658-664.

4. Вайман А.М., Спивак С.И. В сб.: Вопросы математического моделирования и механики сплошных сред. Бирск, 1 997. С.73-81 .

5. Денисов Е. Т. Кинетика гомогенных химических реакций. М. Высшая школа . 1 988.

6. Семенченко А.Е., Соляников В.М., Денисов Е.Т. Нефтехимия. 1971. Т. 11. № 4. С.555-

562.

7. Николаев А.И., Еникеева Л.Р., Сафиуллин Р.Л. Хим. Физика. 1984. Т.3. №5. С.711-714.

8. Мкокеу A.I., Safiullin R.L., Komissarov V.D. React. Kinet. Catal. Lett. 1986. V.31. № 2. P.355-359.

9. Галимова Л.Г., Масленников С.И., Николаев А.И. Известия АН СССР. Серия химическая. 1980. № 1. С.2464-2469.

10. AtHnson R. Chem. Rev. 1986. V.86. № 1. P.69-187.

11. Droege A., Tully F. Journ. Phys. Chem. 1987. V.91. № 5. P. 1222-1225.

12. Druliner T.D., Krusic P.T., Lehr G.F., Tolman C.A. Journ. Org. Chem. 1985. 50. P.5838-

5845.



© 2017 РубинГудс.
Копирование запрещено.