
Быстрай Г.П. Термодинамика открытых систем Часть 2 (2006)
.pdf
Термодинамика открытых систем
ристики полученных хаотических решений.
|
|
10−2 с |
|
|
i* |
|
(а) |
|
|
2 |
|
|
|
|
0 |
|
|
|
|
0 |
|
0.075 |
t,с |
0.15 |
i эксп |
пА |
10−2 с |
(б) |
|
3.5 |
|
|
|
|
0 |
|
|
|
|
|
0 |
t, c |
0.1 |
|
Рис.5.11. а) Динамика тока одиночного ионного канала i*(t)
– решения уравнения (5.7); б) Экспериментальная запись активности Са2+-активируемого канала при [Cа2+]=10 мкмоль/л и потенциале V=20 мВ; данные [6].
При этом, в отличие от классических методов возможности описания проводимости каналов существенно возрастают.
Задачи к главе 5
Задача 5.1. Составьте алгоритм решения уравнения (5.7), которое следует представить в виде системы трех нелинейных дифференциальных уравнений – автономной системы уравнений (5.8). Получите результаты, приведенные на рис. 5.1, или близкие к ним. Проведите анализ статистический анализ гомофазных и гетерофазных флуктуаций внутренней термодинамической силы (а)(параметра порядка), фазовый портрет (б).
Задача 5.2. Используя решения уравнения (5.7) определите хаотическую динамику скорости изменения энтропии G*, т.е. то, что представлено на рис. 5.2 (а,б)). Как ведет при этом
105
Термодинамика открытых систем
параметр порядка η (в)?
Задача 5.3. Найдите хаотические решения, приведенные на рис. 5.3. и эволюцию “расстояния” между двумя расчетными траекториями уравнения (5.7) при заданных отличающихся начальных условиях. Расстояние между двумя соседними траекто-
риями η(t)/ и η(t)// задайте величиной δη(t)= η(t)/ − η(t)// , оп-
ределите tr – характерное время, за которое система забывает начальные условия.
Задача 5.4. Научитесь строить псевдофазовые портреты решений уравнения (5.7) для для любых произвольно заданных ∆; получите частные решения, приведеные на рис.5.4.
Задача 5.5. Постройте алгоритм получения хаотической динамики параметра порядка и показателя Ляпунова λ для термодинамической системы, описываемой отображением (5.14), который приводит к результату, приведенному на рис. 5.5.
Задача 5.6. Постройте бифуркационную диаграмму η= η( a* ) для отображения сборки (5.14) в интервале
−1.6 < a* < 0.1.
ЛИТЕРАТУРА
1.Пригожин И. Введение в термодинамику необратимых про-
цессов. М.: ИЛ, 1960. С.127
2.Николис Г., Пригожин И. Самоорганизация в неравновесных системах. М.: Мир, 1973. С.511.
3.Николис Г., Пригожин И. Познание сложного. М.: Мир, 1990. С.342.
4.Гиббс Дж. Термодинамика. Статистическая механика. М.:
Наука.1982. С.488.
5.Denbigh K.G. Note on Entropy, Disorder and Disorganization. Brit. J. Sci. T. 40 ,1989, p. 323.
6.Климонтович Ю.Л. Введение в физику открытых систем. Соровский образовательный журнал.N 8,1996. С. 109-116.
7.Леонтович М.Л. Введение в термодинамику. Статистическая физика. М.: Наука, 1983. C.416.
8.Ляпунов А.М. Общая задача об устойчивости движения. 2
106
Термодинамика открытых систем
издание Л.-М., 1935. С.
9.Базаров И.П. Термодинамика. М.: Высшая школа, 1976.С.447
10.Меркин Д.Р. Введение в теорию устойчивости движения.
М.: Наука,1987.
11.Быстрай Г.П. Метод функций Ляпунова в анализе открытых термодинамических систем// Вестник кибернетики. N 4.
ИПС СО РАН 2005. С.122-137.
12.Гленсдорф П., Пригожин И. Термодинамическая теория структуры, устойчивости и флуктуаций. М.: Мир, 1973.
С.279.
13.Семенченко В.К. Вступительная статья к книге Дьярмати И. Неравновесная термодинамика. М.: Мир, 1974. С.301.
14.Журавлев В.А. Термодинамика необратимых процессов в задачах и решениях. М.: Наука, 1979. С. 136.
15.Рубин А.Б. Биофизика. T.1. С.448. T.2. С. 467. М.: Книжный дом «Университет», 1999.
16.Gyarmati I. Оn the basic principles of the scattering processes and Its generalisation on the nonlinear promlems. Ann. D. Phys.1969, V.7. P.353.
17.Гилмор Р. Прикладная теория катастроф. М.: Мир.. 1984. T.1. С. 350. T.2. С. 285.
18.Быстрай Г.П., Пивоваров Д.В. Неравновесные системы. Свердловск: Изд-во Урал. госуни−та.1989. С.187.
19.Хакен Г. Информация и самоорганизация: макроскопический подход к сложным системам. М.: Мир, 1991. С.
20.Чернавский Д.С. Синергетика и информация. М.: Наука, 2001. С.245.
21.Кеплен С.Р., Эссиг Э. Биоэнергетика и линейная термодинамика необратимых процессов. М.: Мир, 1986. С. 382.
22.Быстрай Г.П. Методика оценки эффективности энергетических превращений в физических процессах, происходящих при воздействии на горные породы// Изв. вузов. Гор-
ный журнал.1988. N9. C.1.
23.Летников Ф.А. Синергетика геологических систем. Новоси-
бирск: Наука,1992. С.228.
24.Murray J.D. Mathematical Biology. Springer-Verlag. Berlin.
107
Термодинамика открытых систем
Heidelberg. New York. London. Paris. Tokyo. 1984.
25.Соболев С.Л. Локально-неравновесные модели процессов переноса // УФН. Т.167, N 10. C.1095
26.Лыков А.В. Тепломассообмен: (Справочник). М.: Энергия, 1978. С.480.
27.Шустер Г. Детерминированный хаос. Введение. М.: Мир. 1988. C.240.
28.Быстрай Г.П., Студенок С.И., Иванова С.И. Детерминированная модель гомофазных и гетерофазных флуктуаций в системе “жидкость−пар”// ТВТ. 2002. T.40. N 5. C. 779.
29.Быстрай Г.П., Студенок С.И., Иванова С.И. Детерминированный хаос при фазовых переходах первого рода в систе-
ме “жидкость−пар”// ТВТ. 2003. T.41. N 4. C.579.
30.Кольцова Э.М., Третьяков Ю.Д, Гордеев Л.С.и др. Нелинейная динамика и термодинамика необратимых процессов в химии и химической технологии. М.: Химия. 2001. С.407.
31.Лихтенберг А., Либерман М. Регулярная и стохастичесчкая динамика. М.: Мир. 1984. С.528.
32.Заславский Г.М. Стохастичность динамических систем. М.:
Наука, 1984. C.270
33.Быстрай Г.П. Детерминированный хаос при химических реакциях в межфазном слое при высоких температурах//
ТВТ. 2003. T.41. N 6. C.1.
34.Лифщиц E.M.,Питаевский Л.П. Физическая кинетика. Т. 10.
М.: Наука. 1979.С. 527.
35.Гельфер Я.М. История и методология термодинамики и статистической физики.Т.2. М.: Высшая школа. 1973.
36.Эткин В.А. Термокинетика (термодинамика неравновесных процессов переноса и преобразования энергии): Учебное пособие для вузов.- 2-е изд., - Тольятти, 1999.- 216 с.: илл. 17. Библиогр.: 180 назв.
37.Малинецкий Г.Г./ Хаос. Структуры. Вычислительный эксперимент. Введение в нелинейную динамику. М.: Эдиториал УРСС, 2000.
38.Пригожин И.Р. От существующего к возникающему. М.:
Наука, 1985.
108

Термодинамика открытых систем
39.Малинецкий Г.Г., Потапов А.Б. Современные проблемы нелинейной динамики. М.: УРСС, 2002.
40.Быстрай Г.П., Ворох А.С., Андреев С.В. Детерминированный хаос в динамике тока одиночных ионных каналов био-
мембран. Биофизика. 2005. Т.50, вып.5. с. 851-861.
РЕКОМЕНДУЕМАЯ ЛИТЕРАТУРА Николис Г., Пригожин И. Самоорганизация в нерав-
новесных системах. М.: Мир, 1973. С.511.
Николис Г., Пригожин И. Познание сложного. М.:
Мир, 1990. С.342.
Климонтович Ю.Л. Введение в физику открытых систем. Соровский образовательный журнал.N 8,1996. С. 109116.
Базаров И.П. Термодинамика. М.: Высшая школа, 1976.С.447
Меркин Д.Р. Введение в теорию устойчивости движения. М.: Наука,1987.
Гленсдорф П., Пригожин И. Термодинамическая теория структуры, устойчивости и флуктуаций. М.: Мир, 1973.
С.279.
Журавлев В.А. Термодинамика необратимых процессов в задачах и решениях. М.: Наука, 1979. С. 136.
Рубин А.Б. Биофизика. T.1. С.448. T.2. С. 467. М.:
Книжный дом «Университет», 1999.
Гилмор Р. Прикладная теория катастроф. М.: Мир.. 1984. T.1. С. 350. T.2. С. 285.
Кеплен С.Р., Эссиг Э. Биоэнергетика и линейная термодинамика необратимых процессов. М.: Мир, 1986. С. 382.
Гельфер Я.М. История и методология термодинамики и статистической физики.Т.2. М.: Высшая школа. 1973.
Шустер Г. Детерминированный хаос. Введение. М.:
Мир. 1988. C.240.
Кольцова Э.М., Третьяков Ю.Д, Гордеев Л.С.и др. Нелинейная динамика и термодинамика необратимых процессов в химии и химической технологии. М.: Химия. 2001.
С.407.
109
Термодинамика открытых систем
Малинецкий Г.Г. Хаос. Структуры. Вычислительный эксперимент. Введение в нелинейную динамику. М.: Эдиториал УРСС, 2000.
Решения задач главы 1 1.1. Оба выражения в правой части (1.11) являются непол-
ными дифференциалами, так как производство энтропии представляет только часть прироста энтропии [12,15]. Однако производство энтропии при σe=0 можно преобразовать в полный дифференциал, следуя (1.12), так как выполняется:
|
di S |
= − |
1 |
|
dΛF |
|
|
|
(1.23) |
|||||
|
|
. |
||||||||||||
|
|
|
T dt |
|
||||||||||
|
|
|
dt |
|
|
|
σe =0 |
|
||||||
|
|
|
|
|
|
|
||||||||
Для устойчивых по Ляпунову термодинамических систем |
|
|||||||||||||
σ |
i |
= |
|
di S |
≥ 0 , |
dΛF |
|
|||||||
|
|
|
|
|
|
≤ 0 . |
|
|||||||
|
|
dt |
|
|||||||||||
|
|
|
dt |
|
Это и означает, что в реальных процессах свободная энергия уменьшается, т.к. ΛF = F − F0 > 0 .Этот результат является след-
ствием используемого принципа минимальности потенциала в состоянии равновесия.
Из (1.23) получаем для полного дифференциала Пригожина−Гленсдорфа соотношение di S = −(1 / T )dF T ,V [12]. Этот
дифференциал непосредственно связан с изменением свободной энергии Гельмогольца dF. Пригожин и Гленсдорф, не доказывая этого соотношения, объясняли последнее тем, что все процессы протекают в направлении уменьшения F до тех пор, пока свободная энергия не достигнет минимума в устойчивом равновесном состоянии. Наличие такой взаимосвязи между изменениями термодинамических потенциалов для неравновесных состояний и производством энтропии в литературе как правило не обсуждается, кроме единственного упоминания в [12] о взаимосвязи свободной энергии и производства энтропии в виде указанного выше дифференциала Пригожина−Гленсдорфа.
1.2. Производство энтропии в системе как это следует из
(1.11) равно:
110

Термодинамика открытых систем
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