Файл: Быстрай Г.П. Термодинамика открытых систем Часть 2 (2006).pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 29.06.2024

Просмотров: 71

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Термодинамика открытых систем

T

di S

= −

d( ΛF + Λe )

;

dΛe

T

de S

,

(1.24)

 

 

dt

dt

dt

dt

 

 

 

 

 

 

где ΛF=FF0 – термодинамический потенциал неравновесной системы, а Λ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. Для открытой системы энтропия может как увеличиваться так и уменьшаться со временем, так как при стремлении FF0 функция Λ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