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

Замараев - часть 1

.pdf
Скачиваний:
56
Добавлен:
28.03.2016
Размер:
974.78 Кб
Скачать

где ν средняя скорость молекул, l – длина свободного про- бега.

Поскольку l ~ P–1, то D ~ P–1 и kD ~ P–1, где P – дав- ление.

б) Радиоактивный распад атомных ядер

Рассматривая распад разных ядер одного и того же веще- ства как независимые события, получаем для числа ядер N(t)

уравнение

dN(t)dt = - b ×N(t)

ß

N = No × e−βt .

Это закон радиоактивного распада, известный из курса физики.

Константа b называется постоянной радиоактивного распада. По своему смыслу она является константой скоро- сти k1 для радиоактивного распада.

58

Глава 3 Формальная кинетика сложных реакций.

Определение кинетических констант из экспериментальных данных

Выше нами рассмотрены простые реакции, т. е. реакции, протекающие в одну стадию. Единственным исключением были обратимые реакции, рассмотренные в § 2.2 при анализе взаимосвязи между кинетикой и равновесием.

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

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

§ 3.1. Обратимые реакции

Кинетика обратимой реакции первого порядка

k1

А k–1 В.

рассмотрена в § 2.2.

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

В кинетическом отношении все возможные случаи обрати- мых реакций аналогичны друг другу. Поэтому подробно мы рассмотрим только один (наиболее сложный) случай, когда скорости обеих стадий пропорциональны произведению кон- центраций двух веществ.

Это случай реакции

А1 + A2

 

k+

В1

+ B2.

 

 

 

 

k

 

 

 

 

59

Обозначим через a1, a2, b1 и b2 начальные концентрации веществ A1, A2, B1 и B2 соответственно, а через xA1 , xA2 , xB1

и xB2 изменение концентрации этих веществ.

Очевидно, что

 

- xA

 

= - xA

 

= xB

= xB

 

= x .

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

Система кинетических уравнений имеет вид

 

 

 

 

 

 

 

 

 

-

 

dCA1

 

= -

 

dCA2

 

=

dCB1

=

 

dCB2

=

 

dx

;

(3.1)

 

 

 

 

 

 

 

dt

 

 

 

dt

 

 

dt

 

 

 

 

dt

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

= k

+

×(a

 

- x)×(a

2

- x) - k

×(b

 

+ x)×(b

2

+ x) .

(3.2)

 

 

 

 

 

dt

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Константа равновесия в этом случае равна

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(b

~

×(b

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ x)

2

+ x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

KP

=

 

1

 

 

 

 

 

 

 

,

 

 

 

 

 

(3.3)

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

~

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

(a1 - x)

×(a2 - x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где x – значение x при равновесии.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом,

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x является корнем квадратного уравнения

 

KP × (a1 - x) ×(a2 - x) - (b1 + x) ×(b2 + x) = 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ß

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(K -1)× x2

{K (a + a ) + b + b }× x + (K × a ×a – b ×b ) = 0 .

14243

 

14444244443

 

 

 

144424443

 

P

 

 

 

P

 

1

 

 

 

 

1

2

 

 

 

 

 

 

P

 

1

 

 

 

 

1 2

 

a

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

Обозначим трехчлен в левой части уравнения через F(x). Примем для определенности, что a1 £ a2 и b1 £ b2. Очевидно, что прирост CB1 и CB2 не может быть больше a1. Следова-

тельно, интервал между –b1 и a1 есть интервал значений x, имеющих физический смысл.

Значения трехчлена F(x) на концах рассматриваемого ин- тервала:

F(a1) = - a12 - a1 ×b1 - a1 ×b2 - b1 ×b2 < 0 ,

F(-b1) = KP × (b12 + a1 × b1 + a2 × b1 + a1 × a2 ) > 0 .

60

Следовательно, в интервале (–b1, a1) трехчлен меняет знак. Значит, он имеет корень, лежащий в этом интервале. Этот корень имеет физический смысл (так как он лежит в ин- тервале (–b1, a1) и соответствует равновесному значению

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1 = x . Обозначим второй корень трехчлена F(x) через x2. То-

гда F(x) может быть записан в виде

 

 

 

 

 

 

 

 

F(x) =

 

 

 

 

~

- x)×(x2 - x) .

 

 

 

~

(KP -1)×(x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При этом x и x2 связаны соотношением

 

 

 

xˆ + x

2

= –

b

=

1

 

 

×{K

P

×(a

+ a

2

) + b

+ b

} .

(3.4)

 

 

 

 

 

 

a

 

KP -

1

 

 

1

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь обратим внимание на то, что правую часть уравне- ния (3.2) можно преобразовать в трехчлен F(x). Для этого на- до вынести за скобку k. Получим

dxdt = k

Учтем, что

I (x) =

× [KP × (a1 - x) × (a2 - x) - (b1 + x) × (b2 + x)],

 

 

 

 

 

 

 

dx

ß

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= k×F(x) ,

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

ß

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

= k×(KP

-1)×

 

~

- x)×(x2

- x) ,

 

 

 

 

 

 

(x

 

 

dt

 

 

 

 

 

 

ß

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

= (k+ - k)×(x - x)×(x2 - x),

 

 

 

dt

 

 

 

 

 

 

ß

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

dx

 

 

 

 

 

 

 

 

t

 

ò

 

 

 

 

 

 

 

= (k+ - k) ò dt .

 

 

 

 

~

 

- x2 )

 

o

 

(x - x)×(x

 

 

 

 

 

 

o

 

ò

 

 

 

 

dx

 

 

 

=

ò

A dx

+ ò

B dx

 

 

 

 

~

 

 

- x2 )

x

~

 

x - x2

 

 

 

 

(x - x)×(x

 

 

 

- x

 

61

(это разложение дробно-рациональной функции на простые

дроби).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тогда

 

 

 

 

 

 

P(x) = (x - x)×(x - x2 ) ,

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

~

 

 

 

 

~

 

 

A =

 

 

 

,

 

 

B =

 

 

 

 

 

 

 

 

 

 

,

 

 

 

P(x) = x

 

- (x

+ x2 )× x + x

× x2 ,

 

 

¢ ~

 

¢

 

 

 

 

 

 

 

 

 

 

 

 

 

P (x)

 

 

~

P (x2 )

 

 

 

 

~

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

~

 

¢

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¢

 

=

- x2 ),

 

 

¢

 

 

 

P (x) = 2x

- (x + x2 ) ,

 

P (x)

(x

 

P (x2 ) = x2

- x .

 

Итак,

 

1

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

dx

 

 

 

 

1

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I (x) =

 

 

 

× ò

 

 

 

 

 

 

+

 

 

 

 

 

 

 

×

ò

 

 

 

=

 

 

 

 

× ln

x - x

×

~

- x2

 

 

~

 

 

 

 

 

 

 

 

 

 

~

 

x

- x2

 

~

 

 

x - x2

 

x

 

 

 

 

x - x

 

 

 

 

 

x2 - x

 

 

 

 

 

x - x2

 

Кинетическое уравнение имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I (x) - I (0) = (k+ - k)× t ,

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

~

 

 

 

 

ß

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x - x

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

 

 

 

 

 

- ln

 

 

 

 

 

 

= (k

+ - k)×(x - x2 )× t ,

 

 

 

 

 

 

x - x2

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ß

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x - x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

 

 

 

 

 

 

×

 

 

 

 

 

 

=

 

(k

+

- k

)×(x

- x2 )× t ,

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

x - x2

 

 

 

x

 

 

 

 

 

 

 

ß

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x - x

 

 

 

 

= eαt ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x - x2

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

ß

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x × x2 -

 

× x2 = x

 

~

 

αt

 

 

 

 

× e

αt

,

 

 

 

 

 

 

 

 

 

x

× x × e

 

 

- x × x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

ß

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

× e

αt

 

 

 

 

~

 

 

 

×(1- e

αt

) ,

 

 

 

 

 

 

 

 

 

 

 

x ×(x2 - x

 

 

 

) = x × x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ß

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

× x2

 

×

 

 

1- eαt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x = x

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

~

 

 

αt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- x × e

 

 

 

 

 

 

 

 

 

 

 

ß

62

 

 

 

 

 

 

 

 

 

~

-x2 )×(k+ -k )×t

 

 

 

 

 

 

 

~

× x2

×

 

 

e

-( x

-1

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

x = x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

-x2 )×(k+ -k )×t

 

~

 

 

 

 

 

 

 

 

x2 × e

-(x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ß

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

× x2

×

 

 

e(x2 -x)×(k + -k )×t -1

 

 

 

 

 

x = x

 

 

 

 

 

 

 

 

 

 

 

.

(3.5)

 

 

 

 

 

~

 

 

 

~

 

 

 

 

 

 

 

 

x2 × e

(x2 -x)×(k+ -k )×t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- x

 

 

 

Рассмотрим метод определения констант скорости k+ и k

1.

~

 

 

Величина x определяется из опытных данных, если из-

вестен состав смеси при t ® ¥ (при t ® ¥

~

 

x ® x ).

 

2.

~

 

из

Зная x и состав смеси при t = 0, можно найти KP

уравнения (3.3), а затем и x2 по уравнению (3.4).

3. После этого по уравнению (3.5) можно найти (k+ k).

4. Зная (k+ k) и KP = k+ , можно найти k+ и kпорознь. k-

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

§ 3.2. Параллельные реакции

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

На практике чаще всего приходится иметь дело с тремя типами параллельных реакций.

А) Реакции типа

k1

Ak2 X , k3

соответствующие нескольким каналам одного и того же пре- вращения.

63

Учитывая независимость событий превращения по каждо- му из каналов, имеем

-

dCA

= (k1 + k2 + k3 ) × CA ,

 

 

dt

 

и если при t = 0, CA(0) = a и CX(0) = 0, то

CA = a × exp {- (k1

+ k2 + k3 )× t},

Cx = a × [1- exp {-

(3.6)

(k1 + k2 + k3 )× t}].

Очевидно, что превращение будет в основном идти по ка- налу с наибольшей константой скорости.

Б) Реакции типа

 

k1

X1

A

k2

X2 ,

k3

 

X3

 

 

в которых из одного вещества по параллельным каналам об- разуется несколько разных конечных продуктов. Здесь усло- вие материального баланса CA + Cxi = const отсутствует, но

выполняется условие материального баланса

3

CA + åCxi = const.

i=1

Примем, что при t = 0 CA = a и Cxi = 0 для всех веществ

Xi. Тогда для исходного вещества по-прежнему (как и в случае A) имеем

CA = a × exp {- (k1 + k2 + k3 )× t}.

Для продуктов получим

 

 

 

 

 

 

 

ki

× [a - CA ]=

ki

é

ì

üù

 

Cxi =

× a × ê1

- exp í- åki × týú

. (3.7)

 

 

 

 

k

 

 

k

ê

 

ú

 

 

å i

 

å i

î i

 

 

 

ë

þû

 

 

i

 

 

i

 

 

 

 

 

64

Как видим, доля каждого из веществ Xi в образовавшихся продуктах X1, X2, …Xiв любой момент времени составляет

åkiki .

i

Из анализа экспериментальной зависимости CA от t можно

найти åki , а из отношений

Cx

i1

=

k

i1

константы ki порознь.

Cxi2

ki2

i

 

 

В) Реакции типа

k

A ¾¾1 ® X1 ,

k

A + B ¾¾2 ® X2

требуют для кинетического описания решения системы двух линейно независимых уравнений:

 

 

 

 

 

 

 

 

 

ì-

dCA

= k

1

×C

A

+ k

2

× C

A

×C ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

í

 

 

 

dCB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

-

= k

2

×C

A

× C .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

 

 

 

dt

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

î

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Третье

 

 

и

 

 

 

четвертое

 

 

уравнения

 

 

 

dCX

= k1 × CA

и

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

dt

 

dCX2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= k

2

× C

A

× C не являются линейно независимыми, так как

 

 

 

dt

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dCX

1

= -

 

dC

A

+

dC

B

 

 

 

и

 

 

dCX

2

= -

dC

B

.

 

 

 

 

 

 

dt

 

 

 

 

 

dt

dt

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

Дифференциальные уравнения для CA и CB интегрируются в квадратурах, но выражения при этом получаются довольно громоздкие. Однако ситуация существенно упрощается при CB(0) >> CA(0). В этом случае CB(t) можно считать константой, т. е. принять CB(t) = CB(0) = b. Тогда кинетика реакций типа «В» будет описываться теми же уравнениями, что и для реак-

65

ций типа «Б» с эффективной константой скорости для второго канала k2эф = k2 × b, т. е.

CX1

=

 

 

k1

 

 

×a ×[1

- exp {- (k1

+ k2b)× t}] ,

k1

+ k

2 ×b

 

 

 

 

(3.8)

 

 

 

 

k2 ×b

 

 

CX2

=

 

 

×a ×[1- exp {- (k1

+ k2b)× t}] ,

 

k1

+ k2 ×b

 

 

 

 

 

 

где, как и раньше, a = CA(0) – начальная концентрация веще- ства A.

§ 3.3. Последовательные реакции

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

В качестве характерных примеров рассмотрим два случая:

Случай А. Цепочка последовательных необратимых мо-

номолекулярных реакций

A1 ¾¾®k1 A2 ¾¾k2 ® A3 ¾¾k3 ® A4 ¾¾k4 ® ... ¾¾kn¾1 ® An ¾¾kn ® P .

Для рассматриваемого случая система кинетических урав-

нений имеет вид

ì

 

dCA

1

 

= - k1 × CA1 ,

ï

 

 

 

 

 

 

dt

 

 

ï

 

 

 

 

 

 

ï

 

 

dCA

2

 

= k1 × CA1 - k2 × CA2 ,

ï

 

 

 

 

 

 

 

dt

 

 

 

ï

 

 

 

 

 

 

ï

 

 

dCA

 

3

 

= k2 ×CA2 - k3 ×CA3 ,

í

 

 

 

 

 

 

 

dt

 

 

 

ï

 

 

 

 

 

 

ï

 

 

 

 

 

 

 

 

ï ........................................

ï

 

 

 

 

 

 

 

 

ïdCA

n

 

= kn1 ×CAn1 - kn × CAn .

ï

 

 

 

 

 

 

dt

 

 

î

 

 

 

 

 

 

 

Это система из линейных уравнений с постоянными коэф- фициентами. Из курса дифференциальных уравнений из- вестно, что такая система интегрируется. Решение представ-

66

ляет собой линейную комбинацию экспонент eλit , где li кор- ни характеристического уравнения.

В рассматриваемом случае характеристическое уравнение имеет вид

(-k1 - l)

0

0

.....

0

0

 

 

 

 

k1

(-k2 - l)

0

.....

0

0

 

= 0 .

0

k2

(-k3 - l)

.....

0

0

 

.....................................................................................

 

0

0

0

.....

kn1

(-kn

- l)

 

Решив это характеристическое уравнение, можно найти все значения li. Нетрудно увидеть, что l1 = –k1, l2 = –k2, l3 = –k3, … ln = –kn . Далее можно выразить все концентрации

Ci в виде линейных комбинаций экспонент eλit . Мы не будем проделывать все эти громоздкие вычисления. Отметим лишь три обстоятельства.

1. Очевидно, что зависимость от t для CA1 будет иметь вид

CA1 = a × ek1t , т. е. l1 = - k1 . (3.9) 2. Выражение для CA2 будет зависеть только от k1 и k2, но не от ki с i > 2. Следовательно, CA2 можно представить в виде

CA2 = A × ek1t + B × ek2t .

Из условия CA2 = 0 при t =0 получаем следующее соот- ношение между коэффициентами A и B:

0 = A + B Þ B = - A .

Следовательно,

CA2 = A × (ek1t - ek2t ).

Значение коэффициента A можно получить из условия

dCA2

 

= k1

× CA

 

= k1 ×a ,

 

 

dt

 

 

 

t=o

ß

1

t=o

 

 

 

 

 

 

 

 

 

67