Здесь {Xi}, {Yj} и {Zm} имеют то же значение, что и ранее. Соответствующая функция Ляпунова записывается в виде
|
|
|
|
|
|
|
|
|
k |
|
~ |
|
~ |
|
2 |
k,k |
~ |
~ |
|
|
2 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
∑εXi ,Xj |
|
|
|
|
|
||||||
|
|
|
|
|
Φ = ~ |
2∑εA1Xi (A1 |
− Xi ) + |
(Xi − Xj ) |
|
|
+ |
||||||||||||
|
|
|
|
|
|
|
1 |
|
i=1 |
|
|
|
|
|
|
i,j=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
P |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
1 |
|
k,l |
|
|
~ |
~ |
~ |
2 l,m |
~ |
~ |
~ |
2 |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
+ |
~ |
2 |
~ |
2∑εA2 Xi ,Yj |
(A2 |
Xi |
− Yj ) |
+ 2∑εYi ,P1Zj (Yi −P1 |
Zj ) |
|
|
+ |
(3.15) |
||||||||||
|
A |
1 |
|
i,j=1 |
|
|
|
|
|
|
|
i,j=1 |
|
|
|
|
|
|
|
||||
|
|
P |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
n,n |
|
~ |
~ |
|
2 |
|
n |
~ |
~ |
|
2 |
|
|||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
+ |
~ |
2 |
∑εZi ,Zj |
(Zi − Zj ) + 2∑εZi ,P2 |
(Zi |
−P2 ) |
|
. |
|
|||||||||
|
|
|
|
|
|
A |
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
i,j=1 |
|
|
|
|
|
|
|
|
|
|
|
|
||||
Очевидно, что сходным образом можно найти функции Ляпунова и для иных брутто-реакций с линейными по интермедиатам превращениям и произвольным числом исходных реагентов Ai и конечных продуктов Pj.
Существенно, что во всех описанных случаях функционалы Φ являются положительно определенными, что означает устойчивость стационарного состояния для соответствующих брутто-процессов. Для сильнонелинейных кинетических схем условия устойчивости надо рассматривать специально. Так, несложно показать, что положительная определенность функционала Φ действительно является необходимой для устойчивости стационарного состояния соответствующего химического процесса (естественно, в ситуации, если такие функционалы существуют).
Пример 5
Рассмотрим брутто-реакцию
R
P,
которая осуществляется через совокупность промежуточных превращений, имеющих автокаталитические стадии:
R
{Xi}, R + Xi
2Xi,
{Xj}
P.
171
Здесь, как и ранее, {Xi} обозначает произвольную совокупность мономолекулярных превращений интермедиатов Xi
(i = 1, …, k).
Для данной кинетической схемы функция, производная которой принимает нулевое значение в стационарном состоя-
нии относительно «внутренних» переменных ~ , может быть nXi
выражена в виде
|
k |
~ ~ |
2 |
k |
k |
~ |
|
~ |
|
2 |
|
||
Φ = 2∑εR,Xi |
(R − Xi ) |
|
+ ∑∑εXi ,Xj (Xi − Xj ) |
|
+ |
||||||||
|
i=1 |
|
|
|
i=1 j=1 |
~3 |
|
|
~ |
~2 |
(3.16) |
||
k |
~ |
~ |
2 |
|
k |
|
|
|
|||||
|
|
|
|
|
Xi |
|
R |
Xi |
|
||||
+ 2∑εXi ,P (Xi −P) + 4∑εRXi ,2Xi |
|
|
− |
|
|
|
. |
||||||
3 |
|
2 |
|
||||||||||
i=1 |
|
|
|
|
i=1 |
|
|
|
|
|
|||
Действительно, в стационарном состоянии по отношению к интермедиатам Xi мы имеем
|
|
d[Xi ] |
|
1 ∂Φ |
|
~ |
~ |
k |
~ ~ |
||
|
|
|
|
|
|||||||
JXi |
= |
|
= − |
|
~ |
= −εR,Xi |
(Xi −R)− ∑εXi ,Xj (Xi − Xj )− |
||||
dt |
2 |
||||||||||
|
|
|
∂Xi |
|
|
|
j=1 |
)= 0 . |
|||
|
|
|
|
~ |
~ |
|
~ ~ ~2 |
||||
|
|
− εXi ,P (Xi −P)+ εRXi ,2Xi |
(R Xi − Xi |
||||||||
Однако выражение (3.16) не всегда положительно. Это означает, что стационарное состояние рассматриваемой бруттореакции является уже не обязательно устойчивым. На самом деле, автокаталитические реакции при определенных соотношениях концентраций реагентов являются ярким примером процессов с неустойчивостью стационарного состояния.
3.3.3. Примеры анализа устойчивости стационарного состояния брутто-процессов, нелинейных относительно интермедиатов, кинетическими методами
Пример 1
Рассмотрим устойчивость простой нелинейной брутто-
реакции (1.27)–(1.28).
Рассмотрим малое отклонение δ от стационарного значе-
ния ~ ≡ ~ :
Y
Y
172
~ ≡ ~ +
Y Y y .
При этом
[Y] = [Y] + A y ,
где [Y] – стационарное значение концентрации интермедиата
Y, а А – некоторый положительный коэффициент, не зависящий от y.
Тогда, сохраняя члены, линейно зависящие от малого отклонения y, получаем
d[Y] |
|
dy |
|
~ |
|
~ |
|
|
|
|
~ |
2 |
~ |
|
|
|
= A |
|
= ε1δ − ε1 R |
− Y − y |
− 2ε2 |
Y − y |
−P |
≈ |
|||||||
dt |
dt |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
+ 4ε |
|
~ |
|
|
|
|
|
|
|
|
|
|
≈ − ε |
1 |
2 |
Y |
y. |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Таким образом, эволюция y во времени передается функ-
цией |
|
|
|
~ |
|
|
y = y0 |
|
|
ε1 |
|
, |
|
exp − |
+ 4ε2Y |
At |
||||
|
|
|
|
|
|
|
т. е. всегда экспоненциально убывает со временем.
~
Следовательно, стационарное решение Y для рассматриваемой схемы всегда является устойчивым.
Как видно будет далее, для эволюции сложных систем большое значение имеет устойчивость стационарного состояния интермедиатов при постоянстве внешних параметров. Особенно большое значение это имеет для схем превращений, нелинейных относительно интермедиатов.
Пример 2
Рассмотрим брутто-реакцию
R
P,
которая осуществляется по автокаталитической схеме с од-
ним интермедиатом-автокатализатором: R + X ε1 2 X
173
X ε2 P.
Найдем условия возникновения неустойчивости стационарного состояния.
Кинетика изменения концентрации интермедиата Х описывается уравнением
d[X] |
~ ~ ~ |
2 |
~ ~ |
dt |
= ε1(R X − X |
|
) − ε2 (X −P) . |
|
|
|
Найдем приращение диссипации энергии при возникновении малой флуктуации в концентрации (напоре) интермедиата:
~ = ~ +
X X x ,
= ~ − ~ Ar1 RT(lnR ln X) ,
|
|
~ |
|
|
~ |
|
|
|
|
|
|
|
|
|
|
x |
|
|||
|
|
|
|
|
|
|
|
|
|
+ |
|
|||||||||
δAr1 = −RT ln X + x |
−ln X = −RTln 1 |
~ ≈ |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
|||
|
|
|
|
δv2 |
= ε2x , |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
δAr2 = RT |
x |
, |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
~ |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
δP = δv1 δAr1 + δv2 δAr2 = |
|
|
|
|
||||||||||||||
|
~ |
~ |
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
x |
|||
= ε |
R x − 2Xx |
RT − |
|
|
|
|
+ ε |
x |
RT |
|
|
|||||||||
~ |
|
~ |
|
|||||||||||||||||
1 |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
X |
|||
|
|
|
|
|
~ |
~ |
|
|
x2 |
|
|
|
|
|
|
|||||
|
|
= RT |
ε1 2X |
−R |
+ ε2 |
|
|
|
|
. |
|
|
|
|
||||||
|
|
~ |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
X |
|
|
|
|
|
|
||||||||
−RT ~x , X
=
|
Неустойчивость |
|
возникает |
при |
δР < 0, т. е. при |
||||||||
ε |
~ |
~ |
+ ε |
|
> 0 . |
|
|
|
|
|
|
|
|
R − 2 X |
2 |
|
|
|
|
|
|
|
|
||||
1 |
|
|
|
|
|
|
|
|
~ |
|
|||
|
Найдем стационарное значение X : |
|
|||||||||||
|
|
|
|
|
~ |
|
|
~ |
~ |
|
ε2 |
~ |
|
|
|
|
|
|
2 |
|
(ε1R − ε2 )X |
|
|||||
|
|
|
|
|
X |
|
− |
|
|
− |
|
|
R = 0 , |
|
|
|
|
|
|
|
ε |
ε |
1 |
||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
174
|
~ |
ε |
2 |
|
|
ε |
2 |
~ |
|
|
|
|
2 , β = |
|
P . |
||
где α = R − |
ε1 |
|
ε1 |
|||||
|
|
|
|
|
||||
Отсюда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
= α ± |
α |
2 |
+ β . |
|
|
|
|
|
|
||||
|
|
|
|
X1,2 |
|
|
|
|
|
|
|
||||||||
Положительным является только решение |
|
|
|
|
|||||||||||||||
|
|
|
|
|
~ |
= |
α |
2 |
+ β + α . |
|
|
|
|
|
|
||||
|
|
|
|
|
X1 |
|
|
|
|
|
|
|
|||||||
Поэтому |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
~ |
|
+ ε |
|
= ε |
~ |
− 2α − 2 α |
2 |
+ β |
|
+ ε |
|
= |
|||||
ε R − |
2X |
|
|
R |
|
|
|
||||||||||||
1 |
|
|
|
|
2 |
|
1 |
|
|
|
|
|
|
|
|
|
2 |
|
|
~ ~ |
ε2 |
|
|
2 |
|
|
|
|
|
|
|
ε2 |
|
|
2 |
|
|
|
|
= ε1(R −R + |
|
− 2 α + β) + ε2 |
|
|
|
|
|
|
|
|
|
|
|||||||
|
= 2 ε1 |
ε1 |
− α + β . |
||||||||||||||||
|
ε1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Таким образом, неустойчивость возникает при
ε2 − α2 +β > 0 , ε1
или, что то же,
|
|
|
|
|
|
~ |
|
ε |
2 |
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
ε2 |
|
|
R − |
ε1 |
|
ε2 |
~ |
|||||
|
|
|
|
|
|
||||||||
|
|
|
|
> |
|
|
|
|
|
|
+ ε |
|
P . |
ε |
|
|
|
4 |
|
|
|
||||||
|
|
|
|
|
|
|
1 |
||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
Можно сделать чисто кинетическое рассмотрение наличия устойчивости путем рассмотрения кинетики релаксации флуктуации х:
dxdt ~ ε1 R~x − 2X~x − ε2x = ε1 R~ − 2X~ − ε2 x .
Очевидно, что неустойчивость возникает при
ε |
|
~ ~ |
|
− ε |
|
> 0 |
, |
R − 2X |
|
|
|||||
1 |
|
|
|
2 |
|
|
|
175