Принцип максимальности производства энтропии в эволюции макросистем: некоторые новые результаты Текст научной статьи по специальности «Математика»

-□ □-

Запропоновано альтернативний споыб формалiза-ци принципу максимального виробництва ентропи для аналiзу широкого класу макросистем, зокрема термо-динамiчних. На його основi одержано рiвняння динамши таких систем. Найдено динамiчний iнварiант процессу еволюци. Доведено теорему про властивостi замкнутоI нерiвноважноi системи, що розвиваеться зггдно з принципом максимальностi виробництва ентропи. Надаеться визначення «внутрштго» часу системи

Ключовi слова: принцип максимальностi виробництва ентропи, принцип Циглера, еволющя складних

систем, макросистемы

□-□

Предложен альтернативный способ формализации принципа максимального производства энтропии для анализа широкого класса макросистем, включая термодинамические. На его основе получены уравнения динамики таких систем. Найден динамический инвариант процесса эволюции. Доказана теорема о свойствах замкнутой неравновесной системы, развивающейся в соответствии с принципом максимальности производства энтропии. Дано определение «внутреннего» времени системы

Ключевые слова: принцип максимальности производства энтропии, принцип Циглера, эволюция сложных систем, макросистемы -□ □-

УДК 517.956.3+519.246+519.218.7

|DOI: 10.15587/1729-4061.2014.31345|

ПРИНЦИП МАКСИМАЛЬНОСТИ ПРОИЗВОДСТВА ЭНТРОПИИ В ЭВОЛЮЦИИ МАКРОСИСТЕМ: НЕКОТОРЫЕ НОВЫЕ РЕЗУЛЬТАТЫ

Н . И . Дел ас

Кандидат технических наук Кафедра систем управлення летальних аппаратов Национальний авиационный университет пр. Комарова, 1, г. Киев, Украина, 03680 E-mail: [email protected]

1. Введение

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

Значительно реже используется другой вариационный принцип - принципа максимума производства энтропии Циглера. Он «заставляет» систему выбирать не только наиболее вероятное из своих возможных макросостояний, но и наиболее вероятную траекторию движения к этому макросостоянию. Сам Циглер рассматривал его, как более общую формулировку принципа максимума энтропии, «„более точную версию второго фундаментального закона термодинамики», называя ее фазовой версией [2].

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

2. Краткий анализ проблемы и постановка задачи

Возможно, одной из причин меньшей популярности принципа Циглера есть некоторая путаница, связанная с наличием его «антипода» - принципа минимума производства энтропии Пригожина [3]. Нужно заметить, что обстоятельства проявления этих двух принципов различны, и противоречий между двумя полярными формулировками не существует. Так, если принцип Пригожина сформулирован для открытых систем, находящихся в стационарном неравновесном состоянии (в состоянии «проточного равновесия»), то принцип Циглера представляет собой более общее требование, справедливое для нестационарных систем, эволюционирующих к своему стационарному состоянию. Когда система эволюционирует к своему стационарному состоянию, подчиняясь принципу максимума производства энтропии, величина максимума приращения энтропии с каждым последующим шагом уменьшается. Минимум производства энтропии Пригожина, по-существу, означает - минимакс в цепочке уменьшающихся максимумов ее приращений на каждом последующем шаге эволюции. Этот минимум характерен для стационарного неравновесного состояния, которое в термодинамике открытых систем играет ту же роль, что и равновесное состояние для изолированных систем.

Говоря об источниках, следует, прежде всего, отметить работу самого Циглера [2]. В русскоязычной литературе наиболее заметны работы Мартюшева Л. М.

©

и Селезнева В. Д., например, монография [4], где приводится хороший библиографический обзор. Однако в англоязычных источниках принципу максимальности производства энтропии (maximum entropy production principle - MEPP) уделено внимания больше. Среди этих публикаций , можно выделить такие, как [5-9].

В настоящее время формализация принципа Ци-глера в основном опирается на феноменологические представления [5], что производная термодинамической энтропии по времени может быть выражена через некоторые обобщенные величины - сопряженные термодинамические силы X; и потоки J :

dS/dt = £ XJi .

(1)

Решение задачи сводится к выбору значений Х1 и Ь , при которых dS/dt в необратимом процессе имеет положительный максимум. Задача упрощается, при использовании линейных соотношений Онзагера [10]:

где Цк = ьк1 (1,к = 1,2,...,п).

На практике данный подход иногда вызывает трудности из-за феноменологического характера вводимых величин, а также принятого допущения о линейном характере соотношений. Поиск альтернативного варианта формализации упомянутого принципа можно рассматривать как актуальную задачу, имеющую практическую направленность.

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

Стоит обратить внимание, что подавляющее большинство упомянутых макросистем представляют собой объекты, у которых на некотором ограниченном множестве «носителей» распределено ограниченное множество «ресурсов» [11, 12]. Между молекулами газа - распределена энергия, между галактиками - их масса, между людьми - материальные блага, между городами - их жители. Указанные термины удобно использовать для анализа некоторой абстрактной системы, которая могла бы относиться к самому широкому кругу задач.

3. Цель и задачи исследования

Цель настоящей статьи - получить уравнения эволюции некоторой универсальной неравновесной изолированной макросистемы, подчиненной принципу максимальности производства энтропии, выявить ее основные свойства.

В статье решаются такие задачи как: а) разработка альтернативного способа формализации данного принципа, исходя из общих комбинаторных представлений;

б) формализация определения «внутрисистемного» времени;

в) вывод уравнений эволюции для функции распределения числа элементов замкнутой неравновесной системы, подчиняющейся принципу максимальности производства энтропии;

г) доказательство теоремы о свойстве указанной функции распределения;

д) определение динамического инварианта процесса эволюции таких систем.

4. Формализация принципа максимальности производства энтропии

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

Схема описания процесса эволюции выбрана для наглядности дискретной и одномерной. Это не влияет на возможность ее обобщения, позволяя сосредоточиться, главным образом, на основных идеях подхода.

Вначале рассмотрим неравновесную замкнутую систему. На к-м шаге эволюции она состоит из N «носителей», среди которых неравномерно распределено количество Е «ресурсов». Каждый отдельный «носитель» располагает своей индивидуальной порцией е1 этого «ресурса». Можно построить дискретную шкалу разбиения всего диапазона е на М равных интервалов размером Де со значениями е1,... еМ (усреднёнными в пределах интервала). По этому признаку выделится М классов, в каждом из которых на к-м шаге содержится пк - «носителей», обладающих равным количеством индивидуальной доли «ресурса» е1, и соответственно содержится объем «ресурсов», равный Ек = пк е1.

Иными словами, можно выделить М ячеек «фазового» пространства с координатами е1,...еМ (рис. 1), для которых, с учетом принятых обозначений, очевидны равенства:

£ nk = Nk = N - баланс «носителей»,

£ Ek = £nkei = Ek = E - баланс, i=i i=i

«ресурсов», на k-м шаге эволюции,

(2)

(3)

где верхний индекс к - номер шага; нижний индекс 1 - номер ячейки.

= e / e - количество ячеек фазового пространства,

(4)

где Де - выбранный фиксированный объем ячейки.

Стоит заметить, что при решении подобных задач неопределенным часто оказывается верхний предел суммирования М, вычисляемый по формуле (4). Так, если размер ячейки фазового пространства Де может

быть по условию задан (например, в физике - постоянная Больцмана к), то максимальное значение «фазовой координаты» еМ, как правило, неизвестно. Иногда еМ устремляют в бесконечность (например, при расчете статистической суммы в распределении Максвелла [13]). Этот искусственный прием хоть и удобен в вычислительных целях, но допустим лишь для быстро убывающих экспоненциальных распределений: именно таким и является распределение Максвелла. В то же время, для распределений с «тяжелым хвостом» (например, степенных распределений с показателем q ) при расчете статистической суммы устремлять еМ в бесконечность уже нельзя - это приведет либо к слишком завышенному значению суммы, либо - к его неограниченному росту при q < 1.

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

Рис. 1. Распределение «носителей» и «ресурсов» в фазовом пространстве замкнутой макросистемы на ^м шаге эволюции

Рассмотрим динамику распределения элементов внутри замкнутой неравновесной системы (рис. 2), находящейся на к-м шаге своего пути к равновесному состоянию п" (е).

Рис. 2. Изменение количества элементов в нй ячейке на k+1-м шаге эволюции

На к + 1-м шаге, в 1-й ячейке содержится

£дпк е1 =£ДЕк = ДЕк = 0 .

(6)

В основе вариационного энтропийного принципа

лежит идея Больцмана, что система из множества доступных ей макросостояний чаще всего «выбирает» такое из них, которое может быть реализовано максимальным числом различных микросостояний. Образно говоря, реализуется такая внешняя конфигурация системы, которая обладает наибольшим количеством внутренних комбинаторных вариантов, способных ее воспроизвести.

Мощность (кардинальное число) множества микросостояний, благоприятных данному макросостоянию - есть статистический вес W, а его логарифм - есть энтропия Больцмана:

S = lnW.

(7)

(Обычно умножается на константу к=1,38х х10-23 Дж/град, которую выбором размерных величин можно превратить в единицу).

Используя некоторые положения теории разбиений [14] и принятые в начале статьи термины, статистический вес W можно определить, как число упорядоченных разбиений Я(п1,п2 ,...,пМ) на М непересекающихся подмножеств мощности п1, таких, что

М

Уп1 = N. Тут индекс < к > опущен. 1=1 Число упорядоченных разбиений совпадает с общим числом перестановок, за исключением количества тривиальных перестановок внутри каждого подмножества:

W = Я(п1,п2.,...,пМ) = -

№ М .

Пп!

(8)

Далее, применяя к (8) формулу Стирлинга: 1пт! = т ■ (1пт-1), можно энтропию Больцмана представить как:

М М

S = £ п1 ■ ^-Х(п1]дп1),

(9)

1=1

М

или Н = — = -У I —-1п—- . N £1N N '

(10)

В соответствии с определением вероятности, основанном на непосредственном ее подсчете [15], отношение:

Ек+1 = Ек + ДЕк = Ек + Дпк ■ е1 - «ресурсов».

Для замкнутой системы обмен с внешней средой «носителями» и «ресурсами» отсутствует, поэтому:

М

У Дпк = Д^ = 0, (5)

П;

N=р1

(11)

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

и

являются п, равновероятных элементарных событий, из их общего числа N.

Таким образом, соотношение (10) приобретает форму энтропии Шеннона:

Н = -> р, ■ 1пр, .

ДНк = Нк+1 - Нк =

пк + Дп

N

- х

х 1п-

+ > —1- ■ Ш—^

N 1=1 N N

1п

можно представить, как:

ДНк = -уДП^ ■ Ьп_±дп^-уДП^.

1=1 N N и N

Д°' ЧМ > (Дпк)'

(16)

(12)

Несмотря на то, что Н носит узкое название информационной энтропии, как видно, выражения (10) и (12) могут быть получены из анализа достаточно общей схемы.

Можно подытожить, что энтропия Больцмана (7) и тесно связанная с ней энтропия Шеннона (выражения (11), (12)) с точностью до константы определяются, как логарифм кардинального числа класса микросостояний, воспроизводящих данное м а кросостояние, и не зависят от природы «носителей» и «ресурсов» рассматриваемой макросистемы. Поэтому данные формы энтропии можно рассматривать как эквивалентные инструменты, пригодные для анализа не только термодинамических, или информационных, но и более широкого круга макросистем.

Приращение энтропии на ^м шаге эволюции замкнутой системы можно записать, используя (10) и принятые на рис. 2 обозначения. Для двух последовательных распределений пк(е) и пк+1(е) разность энтропии равна:

(13)

док д Дпк, пк+Дпк Дпк, пк+Дпк

ДНк = ->—■ 1п—1,-!—У—1-■ Ь-1—;—. (14)

1=1 N N 1=1 N пк

Выбирая Дпк достаточно малым, таким, что

пк + Дпк Дпк

— ~ —к", приращение энтропии на к - м шаге

Для замкнутой неравновесной системы второе слагаемое в выражении (15) равно нулю.

Контролируемый размер шага в процессе производства энтропии обычно рассматривают, как некоторый промежуток времени, за который происходит ее приращение (см. (1)). В настоящей статье, в качестве размера шага предлагается рассматривать величину параметра, который в неравновесной системе отражает порции внутренних нетривиальных перестановок. Таким параметром может быть среднее квадратичное приращений количества элементов в ячейках фазового пространства:

Формализм для принципа максимума производства энтропии излагаемый ниже, близок к идеям Джейнса, отличаясь видом анализируемого функционала и спецификой уравнений связи.

Он представляет собой процедуру определения такой функции Дпк = ^е) , при которой функционал ДНк {^е)} достигает условного максимума на каждом контролируемом шаге эволюции. В дискретном представлении задача эквивалентна нахождению условного максимума функции многих переменных ДНк(Дп1к,...ДпМ) в присутствии ограничений, задаваемых уравнениями связи. Для ее решения удобен метод множителей Лагранжа, где поиск условного экстремума Дн11(Дп11,...ЛпМ) сводится к определению безусловного экстремума новой функции Фк(Дпк,...ДпМ) , аддитивно включающей искомую функцию ДНк и ограничения, взвешенные множителями Лагранжа а, в, у..., которые впоследствии определяются из уравнений связи. В задаче, сформулированной выше, ограничениями выступают условия (5), (6) и (16), тогда;

г ^пк+дпк „ > ——!— 1 ^ N

Ф к(Дпк,...ДпМ) = ДНк(Дпк,...ДпМ) + а

... + Р

+7'

(17)

Экстремум Фк(Дпк,...ДпМ) определяется из условия:

Эфк-=о,...,^=0.

эдпк

эдпк

Тогда из (17), с учетом выражения (14), получим:

, пк + Дпк и у Дпк

1п—1-- = -1 + а + Ве, + — .

N 1 М Дпк

(18)

Это равенство обеспечивает максимум ДНк, так как его вторая производная - всегда отрицательна (множитель у, как видно из (24), всегда меньше нуля).

Следующей задачей является определение для (18) неизвестных значений а, в, у... Вначале определим сла-(15) гаемые: -1 + а + ве,.

Так как выражение (18) справедливо для произвольного размера шага Дпк, имеем право рассматривать крайний случай - когда этот размер максимален, и после исходного к-го состояния система сразу же переходит в равновесное состояние с распределением п". Соответственно в каждой ячейке приращение Дпк = 8пк, где

8пк = п"- пк.

Тогда из (18), получим:

-1+а + В е, = 1п—1----■^4"

^ 1 хт Л/Т

п" у 8пк

N - М ■Щ1

(19)

Среднее квадратичное 8п определяется аналогично, как в формуле (16):

^ М У н )2 =]ММ У (»"- пк )2.

(21)

Равновесное распределение п" считается известным, так как его можно заранее определить обычным способом, с использованием формализма Джейнса, через требование условного максимума энтропии для равновесного состояния. Исходное распределение пк также считается известным.

Подставляя (20) в исходное равенство (18), получим выражение для приращения количества элементов в ьй ячейке, на к-м шаге эволюции изолированной системы:

а накопленную сумму этих интервалов уДТк - как продолжительность процесса. Это дает определенные основания рассматривать параметр ДТк как некоторое безразмерное «внутрисистемное» время. Используя обозначения (22), а также формулы (23) и (24), выражения для этих величин можно представить как:

(Дпк/^2 ' ДНк/М ,

(26)

Дтк =

1 Д 8пк,

--У =-г1п

М1=1 8пк

N

Дп Дт

8пк хт, пк к -^гтг - N ■ Ь-1-к 8тк

-Дпк

Дтк =--

- к , Дп N

у М

(22)

где приняты обозначения:

-, 8тк =--

8nVN

у М

Аналогично определим выражение для безразмерного «внутрисистемного» времени прихода в равновесное состояние:

8тк =

8пкШ

1 м 8пк, п"

--У =-Чп —-

м£ 8пк N

(27)

Множитель у, можно определить, умножая выражение (20) на 8пк/8пк и суммируя по числу М ячеек. А, так как для замкнутых систем выполняются условия (5), (6), то:

^М 8пк/ . ч, . _ 8Ек Л

У (-1 +а + Р е| ) = (-1 + а)^—= 0,

м8пк

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

М

т=У

8пк

8п^ Ш =-Чп —-8пк N

(23)

Как видно из (26), интервал безразмерного внутрисистемного времени ДТк равен, образно говоря, отношению среднего квадрата от количества внутренних нетривиальных переходов (Дпк/^ к удельному объему происходящих при этом необратимых изменений, выраженных приращением энтропии ДНк/М .

Это в какой-то степени позволяет подойти к объяснению причины, почему в оптике «работает» принцип минимального времени Ферма. Выражение (26) подталкивает к мысли, что принцип Ферма - есть всего лишь проявление более общего принципа - принципа максимума производства энтропии.

Множитель у, можно представить и в иной форме, аналогично умножая (18) на Дпк/Дпк , и тоже, суммируя по М. Тогда, с учетом определения (14), получим:

5. Некоторые результаты, полученные на основе предложенного способа формализации

М Дпк

М ^ Дпк, пк+Дпк

—Чп—1--

N N

ДНк/М

(24)

Приращение энтропии ДНк на к - м шаге можно получить, решив совместно равенства (23) и (24):

Дпк М

ДНк = -Дп- у

8пк1=Г

8пк, п" —Чп—1-N N

(25)

Безразмерное «внутрисистемное» время

ДПк т

Введенный в (22) параметр ДТк = ---— отражает

У /м

определенный объем необратимых внутрисистемных изменений на к-м шаге эволюции. Контролируя ДТк, можно контролировать этот объем. Порции приращений ДТк можно рассматривать как интервалы, характеризующие протекание процесса эволюции,

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

Уравнение эволюции системы, выраженное через «внутрисистемное» время т

Эволюция изолированной системы может быть описана с помощью соотношения (22). Устремляя Дпк к нулю, промежуток Дтк может быть выбран сколь

угодно малым: dтk = 11т Дтк. Тогда в какой-то момент

Дпк

времени, который в дискретном представлении соответствует к-му шагу эволюции, выражение (22) можно

представить как: -Дпк- = ~пт - N ■ 1п-п"-. Это равенство

8пк _ , пк Дтк ' 8тк п1

справедливо для всех значений к, следовательно - в любой момент внутрисистемного времени т можно записать уравнение для изменения приращения числа элементов в 1-й ячейке изолированной системы:

п

= - ln^,

N dt N n

(28) dn|(t) = a(t)-[n"-n,(t)+ß(t)-(lnn"-lnn^t))] . (32)

N■ >(п"-п,(т)}^ Ь^

где а(т) =--!=М- (здесь учтены фор>> (п"- п,(т))2 1=1

мулы (21) и (27)).

Это и есть записанное через «внутрисистемное» время уравнение эволюции замкнутой неравновесной макросистемы при условии подчинения принципу максимальности производства энтропии.

Уравнение эволюции системы, выраженное через обычное время t дпк / N

Умножив выражение (22) на Дтк =--, и разде-

у/м

лив его на промежуток общепринятого времени М, с учетом (23) получим:

Anf = An

"AT = "At

k 'М.-х lnni±AnM

Snk M' n n"

Перейдя к терминам времени, уже можно отказаться от верхнего индекса, понимая, что Лпк = Дn1(t) и 8пк = 8n1(t). Тогда их среднеквадратичные значения (по аналогии с (16) и (21)) тоже соответственно равны:

- .1м 2 An(t)=ЛмЕКоо) и

- П м 2 П м 2

Sn(t) = Z(Sn,(t)) ^ммХ(иГ-n1(t)) , (29)

а, согласно (24):

Y == AH(t)/M M == An(t)/N .

Отсюда предельным переходом получим уравнение эволюции системы, выраженное через обычное время:

^ = a(t)■ (п"-па))-^Ш1 1пп£). (30) dt ^ 1 М ^ п" v у

Здесь введено обозначение:

a(t) = ^-Ит

Sn(t) At^o At

Ж - nj(t))

(31)

Входящее в (30) значение производства энтропии dH(t)

dt

определим, используя выражение (25), учиты-

dt

Здесь введено еще одно обозначение:

= N.j^nr - n](t) ] M 1=1 N N

(33)

Если макросистема близка к равновесному состоянию, уравнение (32) в пределе стремится к известному линейному виду:

dn1(t) / г

~drö c-(n'- ni(t)),

имеющему экспоненциальное решение.

Примечательно, что траектория эволюции замкнутой системы целиком определяется характером исходного распределения n1(t) и характером распределения в равновесном состоянии n". Оба эти распределения считаются известными (функция n" может быть предварительно найдена в соответствии с формализмом Джейнса, как решение на условный максимум энтропии для равновесного состояния системы).

Таким образом, вышеприведенный способ формализации принципа максимальности производства энтропии, основанный на достаточно общих посылках, без использования феноменологических конструкций, позволяет получить уравнение эволюции (32). Оно выявляет то обстоятельство, что вдали от равновесия зависимость становится существенно нелинейной. В этом случае значительную роль играет уже не разность величин n" -n1(t), а умноженная на производство энтропии разность логарифмов: lnn" -lnn1(t).

Возможно, этот факт и объясняет природу феномена «логарифмических очков» при субъективном восприятии человеком количественных величин. Ведь известен логарифмический характер некоторых законов психики (Госсена, Вебера-Фехнера). Пожалуй, вполне правдоподобной могла бы выглядеть гипотеза, что логарифмический масштаб субъективного восприятия у человека обусловлен существенной удаленностью его психики как сложной системы от ее, образно говоря, «равновесного» состояния.

Читателю, которому интересны эти вопросы, можно рекомендовать замечательную книгу Касьянова В. А. - «Субъективный анализ» [16].

Теорема о постоянстве суммы логарифмов

Можно показать, что справедлива следующая теорема. В изолированной (замкнутой) неравновесной системе, эволюционирующей в соответствии с принципом максимальности производства энтропии, в каждый момент времени t выполняется условие:

вая при этом обозначения (29). Тогда:

^lnn1(t) = const = ^lnn" .

(34)

AH(t) At "

1 An(t)

Sn(t) At ti

Sn(t)lnn"

N N

Устремляя Дt к нулю, получим удобное для практического использования уравнение эволюции для распределения в неравновесной замкнутой системе:

Доказательство вытекает непосредственно из соотношения (32). Суммируя обе части этого равенства по числу всех М ячеек, и учитывая, что для

м

изолированной неравновесной системы: >dn1(t) = 0,

1=1

М

У(п"- n1(t)) = 0, то для любого времени t из (32)

1=1 М "

следует: у(1пп"- 1пп^)) = 0. А, значит, справедливо 1=1

утверждение теоремы (29).

Таким образом, из всех возможных траекторий эволюции изолированной неравновесной системы, при условии максимальности производства энтропии, реализуется такое распределением элементов, что изме-

, nl(t) + dnl(t) „

нение суммы логарифмов у 1п—1-——1—в ячейках

с положительным приращением dn (t)( ), компенсируется изменением остальной части суммы в ячейках с отрицательным приращением dn1(t)(-).

Условие (29) можно трактовать также как постоянство суммы логарифмов вероятностей распределения в момент времени t :

равенство (37) справедливо, если выполнено требование:

df(t, e) df(t, e) de

at

Эе dt

= 0,

представляющее собой одномерную версию теоремы Лиувилля.

По условию сформулированной выше задачи рассматривалось распределение на множестве «носителей» лишь одного вида «ресурса», где 8 - фазовая координата этого распределения. В общем случае на одном и том же множестве «носителей» могут быть распределены различные виды «ресурсов», каждый из которых формирует свои фазовые координаты ц,V... . Тогда, при условии взаимной независимости этих распределений, вместо уравнения (37), нужно было бы записать:

£ lnp,(t) = const.

J

1 f 3f(t,e, ц,...) 3f(t,e, ^,...)de 3f(t,e, ц,...) d|i

f(t, e, ц...)

dt

de

dt

Эц

dt

+ ... Ide d^.. = 0 . (38)

Здесь Pj(t) =

n,(t) N

вероятность того, что в момент

f(t, e) = lim

P,(t)

де^о Ae

J lnf(t, e)de = const = I.

Множитель

f(t, e)

Учитывая, что множитель

1

f(t, е,ц...)

имеет положи-

времени t взятый наугад один из N элементов системы окажется «жителем» 1 - й ячейки. Если р1(т) разделить на ширину ячейки Де = еМ/М (см. формулу (4)), впоследствии устремляя ее к нулю, получим плотность вероятности с координатой е , в момент t:

тельный знак по определению вероятности, условие (38) выполняется, если:

9f(t, е, ц,...) 9f(t, е, |i,...)de 9f(t, е, |a,...)d|i

9t

Эе

dt

Эц

dt

+... = 0.

(35)

Тогда доказанная выше теорема (29) может быть сформулирована и для непрерывных распределений. Если неравновесная замкнутая система в соответствии с принципом максимума производства энтропии движется к своему равновесному состоянию, то в любой момент времени t выполняется условие:

(36)

Величину I можно рассматривать как динамический инвариант процесса эволюции системы. Его значение можно определить, имея распределение f"(е) для равновесного состояния:

eM

I = J lnf"(e)de .

0

Связь с теоремой Лиувилля о сохранении фазового объема.

Из (36) следует равенство нулю производной по времени:

fe) + Э£(М) del de = 0 . (37)

J0 f(t, e) I, at Эе dt) '

всегда положителен, поэтому

Последнее равенство, представляет собой многомерный вариант теоремы Лиувилля. Таким образом, если неравновесная замкнутая система в соответствии с принципом максимума производства энтропии эволюционирует к своему равновесному состоянию, то в любой момент времени t выполняется условие теоремы Лиувилля о сохранении фазового объема.

6. Выводы

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

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

1

п^), а также характером распределения в состоянии равновесия системы п". Показано, что вдали от равновесного состояния определяющую роль в уравнении эволюции играет не разность этих величин, а разность их логарифмов.

Дано математическое определение «внутрисистемного» времени.

Доказана теорема, утверждающая, что для замкнутой неравновесной системы, развивающейся в соответствии с принципом максимального производства энтропии, в каждый момент времени реализуется та-

кое распределение элементов по ячейкам фазового пространства, которое сохраняет постоянной сумму из логарифмов их количества. Приведена ее формулировка для случая непрерывных распределений. Показана связь с теоремой Лиувилля о сохранении фазового объема.

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

Литература

1. Jaynes, E. T. Where do we Stand on Maximum Entropy?' in The Maximum Entropy Formalism [Text] / E. T. Jaynes; R. D. Levine, M. Tribus (Ed.). - M. I. T. Press, Cambridge, 1979. - 105 p.

2. Циглер, Г. Экстремальные принципы термодинамики необратимых процессов и механика сплошной среды [Текст] / Г. Ци-глер. - М.: Мир, 1966. - 134 с.

3. Пригожин, И. Введение в термодинамику неравновесных процессов [Текст] / И. Пригожин. - М.: Изд-во иностр. лит., 1960. - 127 с.

4. Мартюшев, Л. М. Принцип максимальности производства энтропии в физике и смежных областях [Текст] / Л. М. Мартю-шев, В. Д. Селезнев. - Екатеринбург: ГОУ ВПО УГТУ-УПИ, 2006. - 83 с.

5. Martyushev, L. Maximum entropy production principle in physics, chemistry and biology [Text] / L. Martyushev, V. Seleznev // Physics Reports. - 2006. - Vol. 426, Issue 1. - P. 1-45. doi: 10.1016/j.physrep.2005.12.001

6. Ozawa, H. The second law of thermodynamicsand the global climate system: A review of the maximum entropy production principle [Text] / H. Ozawa, A. Ohmura, R. D. Lorenz, T. Pujol // Reviews of Geophysics. - 2003. - Vol. 41, Issue 4. - P. 1-24. doi: 10.1029/2002rg000113

7. Kleidon, A. Non-equilibrium Thermodynamics and the Production of Entropy: Life, Earth and Beyond, Springer Verlag, Heidelberg principle [Text] / A. Kleidon, R. D. Lorenz, (Eds.) // Reviews of Geophysics. - 2005. - Vol. 41. - P. 1018-1041.

8. Dewar, R. C. Maximum entropy production and the fluctuation theorem [Text] / R. C. Dewar // Journal of Physics A: Mathematical and General. - 2005. - Vol. 38, Issue 21. - P. L371-L381. doi: 10.1088/0305-4470/38/21/l01

9. Niven, R. K. Steady state of a dissipative flow-controlled system and the maximum entropy production principle [Text] / R. K. Niven // Physical Review E. - 2009. - Vol. 80, Issue 2. - P. 021113. doi: 10.1103/physreve.80.021113

10. Гроот, С. Неравновесная термодинашка [Текст] / С. Д. Грот, П. Мазур. - М.: Мир, 1964. - 456 с.

11. Делас, Н. И. Предельно гиперболический закон распределения в самоорганизованных системах [Текст] / Н. И. Делас, В. А. Касьянов // Восточно-Евроропейский журнал передовых технологий. - 2012. - Т. 4, № 4 (58). - С. 13-18. - Режим доступа: http://journals.uran.ua/eejet/article/view/4901/4543

12. Делас, Н. И. Эволюция сложных систем с гиперболическим распределением [Текст] / Н. И. Делас // Восточно-Европейский журнал передовых технологий. - 2013. - Т. 3, № 4 (63). - С. 67-73. - Режим доступа: http://journals.uran.ua/eejet/article/ view/14769/12571

13. Левич, В. Г. Курс теоретической физики. Т. 1 [Текст] / В. Г. Левич. - М.: «Наука», 1969. - 912 с.

14. Эндрюс, Г. Теория разбиений [Текст] / Г. Эндрюс; пер. с англ. - М.: Наука. Главная редакция физико-математической литературы, 1982. - 256 с.

15. Вентцель, Е. С. Теория вероятностей [Текст] / Е. С. Вентцель. - М.: Наука. Главная редакция физико-математической литературы, 1969. - 576 с.

16. Касьянов, В. А. Субъективный анализ [Текст] / В. А. Касьянов. - Киев: НАУ, 2007. - 512 с.