Материал: Лекции по термодинамике неравновесных систем. Пармон

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

что полностью совпадает с термодинамическим рассмотрением.

Пример 3

Рассмотрим устойчивость стационарного состояния для брутто-процесса

R → P,

осуществляющегося по автокаталитической схеме Лотки– Вольтера с двумя интермедиатами Х и Y:

 

А + Х 2Х,

 

 

k1

 

 

k2

(3.17)

 

Х + Y 2Y,

 

k3

 

 

Y Р.

 

В обычной кинетической записи

 

d[X]

= k1 [A] [X] −k2 [X] [Y] ,

 

dt

 

 

d[Y] = k2 [X] [Y] −k3 [Y] . dt

Единственным ненулевыми стационарными решениями для схемы (3.17) являются

[X]o = k3 / k2, [Y]o = k1A / k2.

Пусть в момент пребывания системы в стационарном состоянии (Xo,Yo) возникла флуктуация концентраций сразу обоих интермедиатов:

[Х] = [Х]о + х,

[Y]= [Y]о + y.

Вслучае малости х и y кинетика развития этой флуктуации

вокруг точки ([X]o,[Y]o) описывается системой дифференциальных уравнений:

dxdt = (k1[A] −k2 [Y]o )x −k2 [X]o y = −k3 y ,

176

dydt = k2 [Y]o x + (k2 [X]o −k3 )y = k3 [Y]o x .

Для нахождения решения этой системы дифференциальных уравнений надо решить вековое уравнение

 

0 − λ

− k3

 

= 0,

 

 

 

 

 

k3 [Y]o

0 − λ

 

 

 

 

или, что то же,

 

 

 

 

 

 

λ2 + k32 [Y]o = 0 .

 

 

Отсюда

 

 

 

 

 

λ1,2 = ±ik3 [Y]o = ±ik3

k1[A]

.

k2

 

 

 

 

 

 

Таким образом, состояние системы после возникновения флуктуации из единственного ненулевого стационарного состояния

 

 

 

 

 

[X]o = k3 / k2,

 

 

 

 

 

 

[Y]o = k1[A] / k2

 

 

описывается выражением

 

 

 

 

[X]

 

[X]

 

 

 

 

 

 

 

 

 

o

 

 

 

 

 

 

=

 

 

+

exp(λ1t) +

exp(λ2t) .

 

[Y] t

 

[Y]o

 

 

 

 

Поскольку λ1 и λ2 – мнимые величины, это означает, что концентрации интермедиатов Х и Y осциллируют вокруг стацио-

нарного состояния ([Х]о, [Y]о) с частотой ω = k3

 

k1

[A] :

 

k2

 

 

 

 

 

 

 

 

 

 

 

 

[X]

 

[X]

 

 

 

 

 

 

 

 

 

 

 

 

o

 

 

 

 

ωt .

 

 

 

=

 

 

+

cos ωt +

 

sin

 

 

[Y] t

 

[Y]o

 

 

 

 

 

 

 

177

Это соответствует ситуации «центр» по Ляпунову (см. рис. 3.1 и подразд. 3.5.4).

Рис. 3.1. Фазовая диаграмма, описывающая незатухающие периодические решения вокруг стационарного состояния для схемы Лотки– Вольтера

Эволюцию рассматриваемой системы, но со схемой пре-

вращений, в которой все реакции обратимы

ε1

А + Х ε 2 Х,

Х+ Y ε3 2 2 Y, Y Р,

проще анализировать с использованием термодинамической форме записи кинетических уравнений

При этом

 

 

 

 

 

 

 

 

)= 0 ,

d[X]

 

~ ~ ~

2

 

 

~ ~ ~

2

dt

 

= ε1

(A X − X

 

)− ε2 (X Y − Y

 

 

d[Y]

~ ~

 

~

2

~ ~

 

 

 

dt

= ε2 (X Y − Y

 

)− ε3 (Y −P).

Стационарное решение находится из уравнений

~

~

 

~2

 

~

 

~

~2

) = 0

,

ε1(A

X − X

) − ε2 (X

Y − Y

~

~

 

~2

) − ε3

~

~

 

 

ε2 (X Y

− Y

(Y

−P) = 0 .

 

Отсюда

 

 

 

 

 

 

 

 

 

 

~2

 

ε3

 

~ ~ ε3 ~

 

 

 

 

 

 

 

 

 

P = 0

 

Y +

ε2

− X Y −

 

 

 

 

 

 

 

ε2

 

 

и, следовательно,

178

~

1

~ ε3

1

 

ε3

~ 2

 

ε3 ~

Y = −

 

− X +

±

 

 

 

 

+

P ,

 

 

 

2

4

 

ε2

− X

 

 

 

ε2

 

 

 

ε2

и т. д.

Развитие флуктуаций термодинамических напоров описывается системой дифференциальных уравнений:

 

dx

~

~

~ ~

~

λx

 

= ε1(A x

− 2X

x)− ε2 (Xy + Yx − 2Yy),

dt

 

~ ~

 

~

 

λy

dy

 

 

dt

= (Xy + Yx − 2Yy)− ε3y.

 

 

 

 

 

 

~

X

~

и Y

где коэффициенты λх и λy служат для перевода значений напоров в концентрации.

Отсюда

λx dxdt

λy dy

dt

~

~

~

~

~

= (ε1A − ε1

2X

− ε2Y)x + (− ε2 X + ε2

2Y)y;

~

~

~

 

 

= ε2Yx + (ε2 X

− 2ε2Y − ε3 )y.

 

 

Для описания эволюции флуктуаций надо найти показатели экспонент в решении этой системы уравнений:

 

~

~

 

~

− ε

~

+ ε

 

~

 

 

 

 

 

 

 

ε R − ε

2X − ε

Y − λ

X

 

2Y

= 0

,

1

1

~

2

 

~

2

 

~

2

 

 

 

ε2Y

 

 

ε2 X − 2ε2Y − ε3 − λ

 

 

~ ~

где X и Y – значения в стационарном состоянии.

При a priori заданной кинетической необратимости обеих нелинейных стадий схемы превращений стационарные значения

~

= ε3

ε2 ,

 

 

 

 

Xo

 

 

~

~

ε2 ,

 

 

 

 

Yo = ε1R

 

 

и

 

 

 

 

 

 

~ ~

 

~

 

 

 

 

 

 

ε1R − ε2Y − λ

 

− ε2 X

= 0.

 

 

~

~

− λ

 

 

ε2Y

ε2 X − ε3

 

179

Отсюда в точке ненулевого стационарного состояния

− λ

− ε

3

 

= 0.

 

~

 

 

ε1R

− λ

 

 

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

λ

2

~

= 0,

+ ε1ε3 R

λ1,2 = ±i

 

~

ε1ε3R ,

что соответствует решению, найденному традиционными кинетическими методами.

§ 3.4. Критерии эволюции в нелинейной термодинамике. Универсальный критерий эволюции Гленсдорфа–Пригожина

В открытых сильнонеравновесных системах линейные соотношения Онзагера не выполняются. Универсальным критерием эволюции таких систем является обобщение принципа минимальной скорости производства энтропии на нелинейные процессы. Ниже рассмотрено содержание этого критерия для систем как однородных, так и имеющих пространственную неоднородность.

3.4.1. Универсальный критерий эволюции для пространственно однородной системы

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

ренних необратимых процессов: dР ≤ 0 и ddtP ≤ 0 .

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

180

Источник: https://studfile.net/preview/16555341/