Файл: Быстрай Г.П. Термодинамика открытых систем Часть 2 (2006).pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 29.06.2024
Просмотров: 85
Скачиваний: 0
Термодинамика открытых систем
T |
di S |
= − |
d( ΛF + Λe ) |
; |
dΛe |
≡ T |
de S |
, |
(1.24) |
|
|
|
|||||||||
dt |
dt |
dt |
dt |
|||||||
|
|
|
|
|
|
где ΛF=F−F0 – термодинамический потенциал неравновесной системы, а Λe − термодинамический потенциал внешней среды. Последнее и является ответом на вопрос в данной задаче. Если перейти к дифференциалам, то из (1.24) получаем результат, приведенный А.Б. Рубиным [15]. Аналогичные соотношения могут быть получены для других потенциалов.
1.3. После совершения одного оборота цикла через τ система вновь вернется в первоначальное состояние, следовательно скорость продуцирования энтропии, или диссипации энергии согласно (1.24), в единицу времени равна определиться уравнением (1.25). В этом случае изменение значения термодинамического потенциала неравновесной системы через время τ0 будет равно нулю ∆Λ*F τ=0. В (1.25) τ0 − время совершения одного оборота цикла (считается достаточно малым). Для внешней среды ∆Λe τ≠0, так как именно за счет взаимодействия с внешней средой и совершается оборот цикла с производимой им за это время работой. Это и доказывает что протекание неравновесных процессов в цикле сопровождается остаточными изменениями в окружающей среде.
1.4. Уравнение (1.25) позволяет сравнивать между собой различные циклы в отношении их энергетической эффективности. Действительно, если имеются две системы, для которых
∆Λ*1 = ∆Λ*2 , то при τ1< τ2 следует что β1 f β2 . Иными словами,
скорость диссипации энергии в первом цикле больше, чем во втором, при том же значении совершенной работы. Этим самым доказывается результат, полученный впервые Т. Мицунойей.
1.5. Для открытой системы энтропия может как увеличиваться так и уменьшаться со временем, так как при стремлении F→F0 функция ΛF(t) в (1.12) уменьшается во времени dΛF(t)/dt<0, а при удалении/отклонении от состояния равновесия dΛF(t)/dt>0. Таким образом, уменьшение энтропии является неустойчивым по Ляпунову процессом, т.е. оно не выполняется на бесконечном интервале времени.
111
Термодинамика открытых систем
1.6. Для доказательства выделим в структуре обратимых потоков через границу, составляющую с теплом d0S/dt:
|
|
|
|
|
|
|
|
|
de S |
= |
d0 S |
+ |
de/ S |
|
, |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
dt |
|
|
dt |
|
|
dt |
|
|
|
|
|
|
|
|||||
de/ S/dt − |
|
|
все остальные потоки через границу. |
В результате с |
||||||||||||||||||||||||
учетом уравнения (1.13) и неравенства (1.14) получаем: |
|
|||||||||||||||||||||||||||
|
d |
i |
S |
|
dS |
d |
0 |
S |
|
d / S |
|
|
|
dS |
|
|
dU |
0 |
|
dV |
|
|||||||
T |
|
|
|
=T |
|
−T |
|
|
|
|
+ |
|
e |
|
|
=T |
|
|
|
− |
|
− P |
|
≥ 0 |
||||
dt |
|
dt |
|
dt |
|
dt |
dt |
|
dt |
|
dt |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
приde/ S / dt = 0 . Из этого неравенства следует неравенство
(1.26). Неравенства типа (1.26) обычно изучаются в учебных курсах по равновесной термодинамике и приводятся без строгого доказательства.
Решения задач главы 2
2.1. Представим коэффициент Lii в виде:
Lii ( X i ) = k1 − k2 X i + k3 X i2 .
Тогда из (2.11) следует (2.12), в котором
с = |
L0ii |
|
|
b1 |
= |
|
|
|
Lie |
|
, |
|
|
b2 |
= |
Lei |
; |
|||||||||
; |
|
|
|
|
L |
L |
|
|
|
L |
L |
|||||||||||||||
|
Lee |
|
|
|
|
|
|
|
|
|
ee |
|
|
|
|
|
|
ee |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
ii |
|
|
|
|
|
|
|
|
ii |
|
||||
|
с |
= |
|
k2 |
|
b1 |
|
X e |
|
2 |
; |
c |
2 |
= |
c1 |
. |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
1 |
|
|
|
|
|
|
Lie |
|
|
|
|
|
|
|
|
|
3x*0 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Здесь мы учли, что принцип симметрии коэффициентов для нелинейных процессов не выполняется: Lie ≠ Lei , для уп-
рощения предполагалось выполнение равенства X c ≈ X e . Здесь
также процесс с индексом “e” “приводит в движение” процесс “i” и 0 ≤ φ ≤1 при условии, что знаки у величин сy и b1 и b2 различ-
ны. Расчеты по (2.12) показали, что кривая энергетических превращений для нелинейных процессов может лежать как выше так и ниже кривой для линейных превращений (рис. 2.4б). Числен112
|
|
|
|
|
Термодинамика открытых систем |
|||||||||||||
ные расчеты также показали, |
что наибольшие отклонения φ |
|||||||||||||||||
имеют место |
в области максимума кривой в сторону превыше- |
|||||||||||||||||
ния |
эффективности |
линейных |
|
процессов при с1 ≈ с2 : |
||||||||||||||
x* = |
x* |
+ x* |
1 |
|
|
|
|
|
|
x* + x* |
|
2 |
|
|
|
|||
1 |
2 |
≈ |
|
|
, |
или |
|
≈ |
|
. |
|
|
||||||
|
|
|
|
|
|
|
|
|||||||||||
0 |
|
2 |
|
3 |
|
|
|
|
|
|
1 |
2 |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
Решения задач главы 3 |
|
||||||||
3.1. В уравнении (3.14) источник тепла равен |
||||||||||||||||||
|
|
|
|
|
|
|
W |
|
|
T 2 |
|
∂σe |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
= βT − αT |
3 |
|
||
|
|
|
|
|
|
|
C |
ρ = |
|
|
|
. |
||||||
|
|
|
|
|
|
|
C ρ |
∂T |
|
|||||||||
|
|
|
|
|
|
|
V |
|
|
V |
|
|
V |
|
|
|
|
|
В случае линейного источника имеем α=0. После интегрирования последнего уравнения получаем, что в данной задаче функ-
ция источников σe содержит слагаемые разных знаков, характеризующие стоки и источники тепла (рис.3.1):
σe = − |
C ρα |
1 |
T 4 |
|
β |
T 2 |
|
||
V |
|
|
− |
|
. |
||||
T02 |
4 |
2α |
|||||||
|
|
|
|
|
|
В результате скорость изменения энтропии (кинетический потенциал) в такой задаче является сложной функцией температуры:
dS |
= σe + J q X q = − |
αСV ρ |
T 4 |
+ |
βCV ρ |
T 2 |
+ |
λ |
( T )2 . |
||
dt |
|
2 |
2 |
||||||||
2 |
|
|
|
|
|
||||||
|
|
4T0 |
|
|
2T0 |
|
T0 |
|
|||
Производство энтропии в этой задаче равно |
|
||||||||||
|
σi = J q X q |
= |
λ |
( T )2 . |
|
|
|||||
|
|
|
|
||||||||
|
|
|
|
|
|
T 2 |
|
|
|
||
|
|
|
|
|
0 |
|
|
|
|
|
Можно доказать, что теорема Пригожина для описываемого уравнения выполняется, т.к. в такой нелинейной системе
имеется одно стационарное состояние при σe = const .
− σe
(αСV ρ/ T02 )
0
0.2 |
|
113 |
|
|
|
|
|
|
|
|
|
0.4 |
2 |
1 |
0 |
1 |
T |
Термодинамика открытых систем
Рис.3.1. Зависимость нелинейного источника тепла для нелинейного уравнения теплопроводности (3.14) при
β/α = 1.2 .
Решения задач главы 4 4.2. В этом уравнении функция внешних источников равна
W |
|
|
T 2 |
|
∂σe |
= βT − αT 3 . |
||
|
|
= |
0 |
|
|
|
|
|
C |
ρ |
C |
|
∂T |
||||
|
ρ |
|
|
|||||
V |
|
|
V |
|
|
|
V |
|
После интегрирования последнего уравнения получаем, что в
данной задаче функция источников σe содержит слагаемые разных знаков, характеризующие источники и стоки тепла соответственно:
σe = − |
C ρα |
1 |
T 4 |
|
β |
T 2 |
|
||
V |
|
|
− |
|
. |
||||
T02 |
4 |
2α |
|||||||
|
|
|
|
|
|
В результате в такой задаче скорость изменения энтропии (кинетический потенциал) является сложной функцией температуры и содержит слагаемые разных знаков:
dS |
= σe + J q X q = − αСV ρT 4 |
+ βCV ρT 2 + |
|
||
dt |
4T02 |
2T02 |
|
T Tτ• . |
|
|
+ |
λ |
( T )2 − |
λ |
|
|
|
|
|||
|
|
T 2 |
T 2 |
T |
|
|
|
0 |
|
0 |
|
Характерно, что данный тип нелинейности не приводит к нарушению теоремы Пригожина при постоянных граничных услови-
ях (σe=const).
Решения задач главы 5
114
|
Термодинамика открытых систем |
|
|
|
|
|
|
|
Термодинамика открытых систем |
|
||||||||||||||||
ТЕРМОДИНАМИКА ЛОКАЛЬНО-НЕРАВНОВЕСНЫХ СИСТЕМ |
Быстрай Г.П. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
Алгоритм решения уравнения (2.18) |
|
|
|
|
|
|
|
|
Потенциальная функция задается значениями переменной. Несмотря на ее |
|
|||||||||||||||
Задаются |
Задаются параметры уравнения |
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
хаотический характер эта функция прорисовывается хорошо |
|
|
||||||||||||||||
начальные |
амплитуда внешней силы, |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
условия |
частота внешней силы, |
|
время релаксации, |
|
|
|
|
|
1 |
(Zn,1) |
4 |
1 |
a (Zn,1) |
2 |
|
|
|
|
||||||||
|
время последействия |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
0.3 |
b := 1.7 |
|
ω := 2 |
|
θ := 1 |
τ := 0.24 |
ω θ = 2 |
a := −1.5 |
|
Fn,1 := 4 |
|
+ 2 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
x:= 0.02 |
|
Решение уравнения (системы автономных уравнений, их три) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
0.01 |
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
3 |
|
|
|
|
|
|
Fn,1 |
|
|
|
|
|
|
Fn,1 |
|
|
|
|
|
D(t,x) := |
− 1 − τ |
3 (x0) |
|
+ a x1 − (x0) |
|
+ a x0 + b (1 |
+ τ ω tan(x2)) cos(x2) |
0 |
|
|
|
|
|
0 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
θ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ω |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Проверка начальных условий |
|
|
|
Оператор решения |
|
Z := rkfixed(x,0,100,10000,D) |
|
|
1 4 |
|
2 |
0 |
2 |
4 |
1 0 |
500 |
1000 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
Z0,1 = 0.3 |
Z0,2 = 0.02 |
|
Z0,3 = 0.01 |
|
Количество |
|
|
|
|
n := 0..10000 |
|
|
|
|
|
Zn,1 |
|
|
|
|
t(n) |
|
||||
|
|
расчетных точек во |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Значение |
|
|
|
|
|
времени |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
переменной в |
Y0:= |
−a |
|
|
|
Переход к |
|
|
t(n) := |
2000n |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|||
аттракторе |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
реальному времени |
10000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Расчет: Динамика |
|
|
|
|
|
Фазовый портрет |
|
|
|
|
t(n) |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
(переменная,скорость ее изменения) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
переменной.Y0 - |
|
|
|
|
|
|
|
t(n) |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
Это странный аттрактор |
|
|
|
500 |
|
|
|
|
|
|
|
|
|
||||||||
значения в аттракторах |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t(n) |
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
1000 4 |
|
|
|
|
|
|
|
|
|
|
Z n,1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
0 |
2 |
4 |
|
|
|
|||
Y0 |
|
|
|
|
|
|
|
|
Z n,1 |
|
|
|
|
|
|
|
|
|
Zn,1,Y0,− Y0 |
|
|
|
|
|
||
0 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− Y0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
2 |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Псевдофазовый портрет хаотических |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Zn+15,1 |
|
|
||||||||
0 |
500 |
|
1000 |
1500 |
2000 |
|
|
|
|
0 |
|
|
пульсаций(зависимость последующих решений |
|
|
|||||||||||
|
|
|
|
|
|
|
от предыдущих), t=15 время задержки |
|
|
|
|
|||||||||||||||
|
|
|
t(n) |
|
|
|
|
|
|
|
|
Z n,2 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
2000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
Расчет: Динамика |
|
|
t(n) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Zn,1 |
|||
|
|
|
скорости изменения |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
переменной. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
0 |
5 |
0 |
|
5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Z n,2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
115 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
116 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Термодинамика открытых систем |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Алгоритм определения показателей Ляпунова |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
0.3000000008 |
Задаются начальные условия, слабо |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
y := |
|
0.02 |
|
|
|
|
отличающиеся от заданных на странице 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
0.01 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Решается второе уравнение (полностью аналогичное приводимому на первой странице) |
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
с начальными условиями, слабо отличающиеся от заданных на странице 1 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
1 − τ 3 (y )2 |
|
|
|
− (y |
|
|
|
|
|
+ b (1 + τ ω tan(y )) |
cos (y |
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
− |
+ a y |
1 |
)3 + a y |
|
2 |
) |
|
|
|
|
|
|
||||||||||||||||||||||||||||
G(t,y) := |
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
|
|
0 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
θ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ω |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Находится расстояние между двумя |
|
|
|
|
|
|
|
M := rkfixed(y ,0,100,10000,G) |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
траекториями |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n := 0.. 10000 |
|
|
|
|
t(n) := |
2000 n |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
δn,1 := |
|
Z |
|
|
− M |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s := 0.. 1164 |
|
|
|
|
|
|
|
|
10000 |
|||||||
|
n,1 |
n,1 |
|
|
Задается показатель Ляпунова |
|
|
|
λ |
Начальное расстояние между |
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
двумя траекториями |
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
λ := 0.018 |
D(s) := |
|
δ |
,1 |
|
|
eλ s |
|
|
|
|
|
− 10 |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
δ0,1 |
= 8 × 10 |
|
|
|
|
|
|
|
|
||||
Наклон пунктирной линии подбирается |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
значениями |
|
λ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
0.1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
0.01 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Расчет времени |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
1 |
.10 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
прогнозирования |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
δ |
n,1 |
|
1 |
.10 |
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
1 |
.10 |
5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10000 |
|
|
|
|
|
|
|
1 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
D(s) |
1 |
.10 |
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L := ∑ |
|
|
|
δn, |
1 |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
1 |
.10 |
7 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4180 |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
1 |
.10 |
8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n = 5820 |
|
|||||||||||
|
|
|
|
|
|
|
9 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
1 |
.10 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
1 |
.10 |
|
10 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
. |
|
|
11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t := |
1 |
ln |
L |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
1 |
.10 |
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L = 0.821 |
|
λ |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
0 |
|
|
500 |
1000 |
1500 |
2000 |
|
|
|
δ0,1 |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t(n) ,s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n := 2000 |
10000 |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
10000 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n = 1 × 10 |
|
|
|
|
|
n := 1164 |
2000 |
117
Термодинамика открытых систем
АЛГОРИТМ ПОСТРОЕНИЯ БИФУРКАЦИОННОЙ ДИАГРАММЫ
ДЛЯ ОТОБРАЖЕНИЯ СБОРКИ |
|
|
|
|
|
|
|
F(τ , N , x0, h , a) := for k 0 .. N |
|
|
|
|
|
||
|
x0 ← x0 |
|
|
|
|
|
|
|
|
|
(xk) |
3 |
+ |
a xk |
|
|
x |
|
|
|
|||
|
← x − h |
|
|
|
|
|
|
|
k+ 1 |
k |
|
|
|
|
|
|
|
1 − τ 3 (xk)2 + a |
|||||
|
x |
|
|
|
|
|
|
Устанавливается шаг |
Выбираются численные значения |
|
|
||||
изменения параметра |
|
|
|||||
других параметров и начальных условий |
|
||||||
|
|
||||||
a := −1.6 , −1.598 .. 0.1 |
|
G( a) := |
F( 0.14 , 100 , 0.01 , 0.6 , a) |
||||
|
|
F( a ) := F( 0.14 , 100 , −0.01 , 0.6 , a) |
|||||
|
Построение бифуркационной |
|
|
|
|
|
|
|
диаграммы |
|
|
|
|
|
|
|
1.5 |
|
|
|
|
|
|
G ( a) 39 |
|
|
|
|
|
|
|
G ( a) 40 |
1 |
|
|
|
|
|
|
F ( a) 39 |
|
|
|
|
|
|
|
F ( a) 40 |
0.5 |
|
|
|
|
|
|
G ( a) 21 |
|
|
|
|
|
|
|
G ( a) 22 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
G ( a) 23 |
0.5 |
|
|
|
|
|
|
F ( a) 21 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F ( a) 22 |
1 |
|
|
|
|
|
|
F ( a) 24 |
|
|
|
|
|
|
|
|
1.5 |
|
|
|
|
|
|
|
1.5 |
1 |
0.5 |
|
|
0 |
|
a
118