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

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

ских ситуациях разница между θ и T невелика. Однако с концептуальной точки зрения появление θ в уравнении для теплового потока означает существование локально-неравновесного стационарного состояния q = −λ θ, которое может существенно

отличаться от локально-равновесного стационарного состояния q = −λ T .

На рис.4.2.представлено отношение температуры θ к равновесной при высокоинтенсивных процессах нагрева и охлаж-

дения для фиксированного значения отношения параметра τT к

T0.

1.4

•

T > 0

1.2

θ/T

1

•

T < 0

0.8

Рис.4.2. Отношение локально-неравновесной температуры к равновесной при высокоинтенсивных процессах нагрева и охлаждения при τT / T0 =1.25 10−5 .

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

75

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

волновые свойства процессов теплопереноса при низких температурах – распространение тепловой волны с конечной скоростью [25], отражение тепловой волны от теплоизолированной границы, а при падении на границу раздела двух сред частичное отражение и частичное прохождение в другую среду; интерференция тепловой волны. Описание процессов горения, плавления, кристаллизации и др. также невозможно без привлечения этих представлений.

Задачи к главе 4

Задача 4.1. Считая уравнения (4.5) нестационарными уравнениям переноса в локально неравновесных системах сформулируйте и докажите для них теорему 1 для устойчивых по Ляпунову процессов.

Задача 4.2. Определить производство энтропии и скорость изменения энтропии для локально неравновесных процессов, описываемых гиперболическим уравнением с нелинейным тепловым источником

76

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

Глава 5. ТЕРМОДИНАМИКА ХАОТИЧЕСКИХ СИСТЕМ

Классическая механика лишила иллюзии, что можно построить вечный двигатель первого рода, термодинамика - второго, квантовая механика - что мы можем одновременно сколь угодно точно измерять координату микрочастицы и ее импульс, теория относительности - что удастся передавать информацию в вакууме со сверхсветовой скоростью [37]. Сегодня нелинейная динамика развеяла иллюзию глобальной предсказуемости: мы не можем предсказать, начиная с какого-то горизонта прогноза, поведение многих достаточно простых систем.

Всвое время работа Лоренца, в которой пожалуй впервые описано возникновение детерминированного хаоса, была опубликована в метеорологическом журнале, но в течение 10 лет она не была замечена. Метеорологи сегодня полагают, что горизонт прогноза для погоды не превышает трех недель. Другими словами, как бы точно не измерялиcь параметры атмосферы, предсказать погоду с помощью имеющихся приборов через три недели в данном месте, вообще говоря, невозможно. Горизонт прогноза для состояния океана эксперты оценивают в месяц.

Сейчас многие специалисты по физике Солнца предполагают, что аналогичная ситуация имеет место с Солнцем. Например, известно такое явление, как минимум Маундера [37], когда

втечение почти 70 лет всплесков солнечной активности не было. И возникает вопрос, можем ли мы предсказать следующий аналогичный минимум? Все это же относится к преблеме глобального потепления или похолодания.

Внелинейной динамике обычно выделяют следующие основные характеристики хаотических движений, которые могут рассматриваться как диагностические [27]:

1.Сложный непериодический характер временной эволюции динамических переменных.

2.Экспоненциальный характер разбегания близких по начальным данным траекторий на аттракторе.

3.Диссипативный характер протекающих процессов, который говорит о сжатии фазового объема.

77

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

4.Непрерывная зависимость спектральной плотности мощности пульсаций от частоты в конечном диапазоне частот.

5.Специфический характер реакции системы на малое изменение параметров и на внешнее воздействие сигналом малой амплитуды.

6.Конечная или бесконечная последовательность бифуркаций, наблюдаемая при вариации некоторого управляющего параметра системы, в результате которых топологическая структура фазовых траекторий претерпевает ряд (конечный

или бесконечный) изменений, завершающихся рождением странного аттрактора.

Анализ необратимых процессов будет неполным без рассмотрения мелко- и крупномасштабных внутренних флуктуаций, которые всегда присутствуют в реальных системах. При наличии локального неравновесия такие флуктуации могут усиливаться и переводить нелинейную систему в новое стационарное состояние, соответствующее новой структуре. Речь по сути идет о неравновесных фазовых переходах в условиях отсутствия локального равновесия. Это означает, что критерий эволюции системы, устанавливающий направление движения начинает зависеть от флуктуаций, а они не могут не отражать присущие системе нелинейные свойства [15]. Поэтому введение в правую часть уравнений типа (2.5), (2.6) источника случайного шума и изучение получаемых следствий при численных расчетах является малоперспективным, так как сразу возникают вопросы по характеристикам этого шума, которые неизвестны. Можно высказать утверждение, что флуктуации должны содержать сами нелинейные уравнения, которые получены выше, но решения указанных уравнений их не содержат. Поэтому рассмотрим некоторые полезные нелинейные модели, основанные на последних, которые приводят к флуктуационным (хаотическим) режимам. Однако вначале сделаем переход от релаксационных уравнений термодинамики к уравнениям второго порядка для параметра порядка, которые возникают в связи релаксацией термодинамических сил и потоков.

Цель этой главы – показать возможность возникновения детерминированного хаоса в термодинамических системах. Под

78

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

детерминированным хаосом будем понимать хаос, возникающий в системе нелинейных дифференциальных уравнений, число которых не менее трех. Хаос возникает при некоторых условиях, которые следует определить.

5.1. Переход от релаксационных уравнений локальнонеравновесных систем к уравнениям второго порядка

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

ний для параметра порядка η(t) η= x* − x*0 и приведенной ско-

рости изменения энтропии G*(t+τr). Последняя является знакопеременной потенциальной функцией и для нее справедливо градиентное уравнение

здесь τr≡τr/t0 − время релаксации скорости изменения энтро-

пии, или время релаксации потока одного из самых длительных неравновесных процессов. Параметр порядка ηсвязан с отклонением термодинамического потока x=Ji (или силы Xi) от среднего значения x0 в приведенном виде

η= x − x0 , xc xc

здесь xс=Jiс – некоторый масштаб потока. В системе уравнений

(5.1)-(5.2) уже два параметра порядка η(t) и G* ( t ) .

Дифференцируя (5.2) по η и подставляя полученное выражение в (5.1) получаем дифференциальное уравнение второго порядка для локально-неравновесных систем

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

где f − обобщенная сила двухямного потенциала G* [19]. Член с τr можно не учитывать когда время релаксации скорости изменения энтропии существенно меньше времени действия внешних сил τr /∆t<<1. Термодинамические уравнения (5.1)-(5.3) характеризуют локально-неравновесные процессы. Следует обратить внимание на то, что в нелинейном уравнении (5.3) сила f, пара-

ны в один и тот же момент времени t.

Можно рассмотреть частный случай, когда внешняя сила H* изменяется по гармоническому закону. Это означает, что управляющие параметры в (5.3) можно представить в виде

b* = −H * + H *s = −H *0 cosωt ,

H 0* = H * − H s* . Наличие в (5.3) времени релаксации и переход к

ДУ второго порядка является необходимым, но недостаточным условием возникновения хаоса.

5.2. Хаос при внешних гармонических воздействиях в условиях запаздывания

Если следовать [28, 29], то в реальных системах следует учитывать также последействие. Если внешнее воздействие представлено в момент времени t−τ, т.е. обобщенная сила

f = −∂G* / dη задана при последействии в виде f (t − τ) . При

делены в момент времени t. В результате приходим к уравнению для параметра порядка η с запаздыванием:

80

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

здесь τ − время последействия. Моделирование последействия в самом простом случае подразумевает разложение силы f(t−τ) в ряд; ограничиваясь двумя слагаемыми представим

В системе уравнений (5.1), (5.2) и (5.6) уже три параметра

порядка η(t), G* ( t ) и f(t−τ). В результате вместо уравнения

(5.3) получаем для локально-неравновесных систем с запаздыванием каноническое однородное уравнение второго порядка, содержащее время ретардации (запаздывания), в котором:

где декремент затухания и амплитуда внешней силы равны соответственно

Γ(t) =1− τ(3η2 + a* ) >0,

H 0`* = H0* (1+ τωtg(ωt)) ;

здесь τ = τ/ t0 − приведенное время ретардации (запаздывания); τr = τr / t0 − приведенное время релаксации скорости изменения энтропии (определяется наибольшим временем релаксации тер-

модинамического потока). Переменный параметр γ = Γ( t ) явля-

τr

ется параметром диссипации. Численные решения нелинейного уравнения для η (5.7), которое представлялось системой трех нелинейных дифференциальных уравнений, показывают на наличие в широкой области значений управляющих параметров не только регулярных, но и хаотических решений (рис.5.1, рис.5.2).

Как видно из рисунков при некоторых значениях пара-

81

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

метров в уравнении имеет место как гомо− так и гетерофазный хаос, реализующийся по типу странного аттрактора. В такой нелинейной термодинамической системе параметр порядка “мечется” между двумя симметричными стационарными состояниями (фазами)

η+ = − a* , η− = − − a* ,

оба из которых являются неустойчивыми. Подставляя хаотические решения в выражение для скорости изменения энтропии G* получаем для этой функции хаотические значения (при a*<0,

рис.5.2).

2

η

0

б

Рис.5.1. Моделирование гомофазных и гетерофазных флуктуаций внутренней термодинамической силы (а), фазовый портрет (б). а*=−1.5, ω = 2.6 , τ = 0.216 , η(0)=0.3, b0*=1.8. Фа-

зовый портрет соответствует двум аттракторам.

Хаотическую динамику решений уравнения (5.7) можно представить как хаотические колебания в одной из потенциаль-

82

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

ных ям с перебросом время от времени в другую яму благодаря периодической внешней силе.

Описываемый хаос в термодинамической системе является детерминированным, т.к. решаемое численными методами уравнение не содержит источников шума; он обусловлен нелинейными особенностями уравнения, проявляемых при периодическом воздействии на систему. Хаос означает наличие состояние беспорядка и нерегулярности. Для таких нелинейных систем динамические законы однозначно определяют эволюцию во времени состояния системы при известной предыстории.

0

в50

Рис.5.2. Хаотическая динамика скорости изменения энтропии G* (а,б) и параметра порядка η (в) (см. условия рис.2.2). Положительные значения G* соответствуют устойчивым состояниям по Ляпунову, отрицательные – структурной устойчивости.

83

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

5.3. Сжатие фазового объема. Диссипативность локально-неравновесной термодинамической системы

Напомним, что если система частиц описывается уравнениями Гамильтона, то при движении частиц фазовый объем остается неизменным. Диссипативные системы обладают той особенностью, что при их движении фазовый объем сжимается. Причем он сжимается к аттрактору более низкой размерности, чем исходное пространство [27, 30]. Покажем, что рассматриваемая нами термодинамическая система, описываемая уравнением (5.7), является диссипативной. Для этого представим ДУ для локально-неравновесной среды в виде автономной системы трех дифференциальных уравнений

H 0`* = H 0* (1+ τω tgZ ) .

Эта система уравнений имеет три степени свободы. Поверхность S, ограничивающая произвольно выбранный в фазовом пространстве {X,Y,Z} фазовый объем V, эволюционизирует так, что каждая ее точка движется по траектории, определяемой системой уравнений (5.8). Отсюда по теореме о дивергенции

или

84