Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Либов Р. Введение в теорию кинетических уравнений

.pdf
Скачиваний:
36
Добавлен:
25.10.2023
Размер:
13.16 Mб
Скачать

5.3. Моментный метод Грэда

297

Самой

простой

 

является

аппроксимация первого

порядка

 

 

 

2г — сг

у J- - г ~2

а«гз\УЩГ) .

 

 

(5.123)

 

 

 

 

= JF0

( і

' 1

 

 

 

 

 

Опуская в (5.119) все члены

с а(3)

получим

 

 

 

datf

 

' +

4

 

 

 

 

 

 

 

 

dt

дхг

дхг 1

"

дхг

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ди;

 

 

„ .

1

DC2

Т<2)

(5.124)

 

 

 

 

■ Ш + МГ- ° І-ІІ

С2

Dt

 

 

 

 

j а .

Если при вычислении / <2) пренебречь всеми членами выше а(2,г то в последнее уравнение будут входить {р, и, Т , а)2’ } и л и , что то же, {р, и, Т, Рц). Поскольку аг)2) —симметричный тензор

с нулевым следом, то он

(а<2)) имеет только шесть компонент.

В сумме с переменными р,

и, Т это составит одиннадцать скаляр­

ных величин, которые определяют состояние системы в прибли­ жении первого порядка.

Тензорное уравнение (5.124) содержит шесть скалярных величин. Эти уравнения вместе с уравнением неразрывности

* + ѵ . р и = о ,

уравнением импульсов

р ( 4 f + u , V u ) + v - ^ = 0

и условием

Тг а<2> = О

составят одиннадцать скалярных уравнений. Система замкнута.

Задача 5.7. Вычислить вектор теплового потока Q, используя аппроксимацию первого порядка (уравнение (5.123)).

Задача 5.8. Используя приближение (5.123), получить явное выражение для / <2) в частном случае взаимодействия максвеллов­ ских молекул. Другими словами, показать, что

где

Ъх= л sin2 Ѳcos2 ѲdQ.

298 Гл. V. Решение и дальнейшие свойства уравнения Больцмана

5.4. Релаксация к равновесному состоянию

а) Дальнейшее толкование необратимости

Мы уже встречались с понятием необратимости в разд. 2 гл. I, в разд. 4 гл. II и в разд. 1 гл. IV. В главе I в основном обсужда­ лась динамическая обратимость. В главе II явление обратимости кратко рассматривалось с точки зрения Ѵ-частичной функции

Р и с . 5.4. Эволюция

числовой плотности, описываемая уравнением диф-

*

фузии.

распределения. В главе IV было замечено, что любой закон, где существенно направление изменения времени, является необрати­ мым законом.

Примером необратимого закона может служить уравнение диффузии

дп дгп

(5.125)

dt дх%

Решение этого уравнения таково, что начальный пик плот­ ности рассеивается со временем таким образом, что устраняется пространственная неоднородность функции п. Процесс прекра­ щается, когда п становится однородной (рис. 5.4).

5.4. Релаксация к равновесному состоянию

299

Явление, происходящее в жидкости согласно уравнению (5.125), можно объяснить геометрически. В области, где кривизна кривой

п

п

Р и с . 5.5. «Сглаживающий» .механизм уравнения диффузии.

п (X) положительна (пхх < 0), п уменьшается, а в области с отри­ цательной кривизной (пхх > 0) п растет. Это ведет к сглаживанию кривой п = п (х) (рис. 5.5).

Задача 5.9. Газ находится в одномерном контейнере (0 ^ х ^

^ L).

В

начальный момент задается

распределение плотности

п (х =

0, t — 0) =

п0, п (х = L,

t = 0) = rii,

п0-< Пі,

 

п (х) — п0+

(гіі щ) для

0 <.ö<Cx*CL — б,

где б

L. Это начальное распределение изображено на рис. 5.6, а.

а)

Доказать, что решение п = щ +

(д4 — тг0)

удовлетворя­

ет уравнению диффузии (5.125) и начальным условиям. б) Почему это решение физически некорректно?

Ответ, б) Оно не дает корректной скорости и (х, t). Скорость и определяется уравнением

71 дп пи — — и — .

ОХ

Так как жидкость ограничена, то в каждый момент времени

и (х — 0) = и (х = L) = 0,

следовательно, начальные условия для п должны также включать (граничные) условия:

дп

дп

=

0 .

дх ж=0

дх

x= L

 

Полные начальные условия для п показаны на рис. 5.6, б.

300 Гл. V. Решение и дальнейшие свойства уравнения Больцмана

Рассмотрим теперь процесс, который подчиняется уравнению (5.125) при обратном течении времени. То есть события развора­ чиваются так, что время t последовательно уменьшается. Этого легко достигнуть, заменяя t на —1' и считая время t' возрастаю­ щим. Тогда обратный процесс описывается уравнением

дп

„ 92ге

(5.126)

dt'

дх2

 

Это уравнение не совпадает с (5.125), и п (х , t) не обязана совпа­ дать с п (х, t'). Однако, очевидно, что положительные значения

п

Р и с . 5.6. Начальные

условия задачи 5.9.

а — заданные начальные данные;

б — полные согласующиеся началь­

ные данные.

пхх вызывают убыль п, а отрицательные значения пхх вызывают

рост п. Таким образом, поведение плотности п противоположно поведению плотности п из уравнения (5.125). Процесс изменения

п показан на рис. 5.7.

Кроме того, различие между уравнениями (5.125) и (5.126) хорошо выявляется при исследовании устойчивости их решений. Пусть мы хотим исследовать устойчивость равновесного состоя­

ния п = п — п0. Для этого внесем возмущения, так что при

5.4. Релаксация к равновесному состоянию

301

, t — 0 плотность будет такова:

п = п = п0 + 8 cos kx = {гс0 + А (X, t)} 1і=0, (5.127)

где е <С и к — действительная положительная постоянная. Рассмотрим решение вида

А (X, і) = 8 cos кх ехрсог.

(5.128)

Оно принимает при t = 0 заданное начальное значение. Чтобы

п

п

п

Р и с . 5.7. Изменение числовой плотности, описываемое уравнением диф­ фузии при обращенном течении времени.

найти со, надо подставить это решение в исходное уравнение

((5.125) либо (5.126)). Подстановка в (5.125) дает

 

и = — Dk2.

(5.129)

Следовательно, начальные

возмущения затухают, и

равновесие

п = п0 устойчиво.

 

 

Для исследования устойчивости решения п = п0 в случае

обратного направления времени подставим решение

 

А =

8 cos кх exp ©Р

(5.130)

302 Гл. V. Решение и дальнейшие свойства уравнения Больцмана

в уравнение (5.126). Получим

+ D k2.

 

со =

(5.131)

Начальные возмущения А (0)

= е cos кх

растут со временем,

и равновесие п = п0 является неустойчивым.

Если уравнение (5.125) выражает истинный закон природы, то в мире с обратным течением времени состояние однородной плотности должно быть неустойчивым. Из этого факта следует простой эксперимент по определению направления течения вре­ мени. Пусть газ однородно заполняет некоторую колбу х). Встрях­ нем ее. Если будут образовываться капельки, то время течет назад. Поскольку существует эксперимент, определяющий направ­ ление течения времени, то уравнение (5.125) описывает необрати­ мый закон.

Для обратимого закона нельзя указать эксперимента, опре­ деляющего направление течения времени. Если бы все фундамен­ тальные законы природы были обратимы, то невозможно было бы определить, куда течет время, назад или вперед. Точно так же, как законы природы не зависят от того, в какой момент времени их начали применять, они также не будут зависеть и от направления течения времени. Время будет однородным и изотропным.

Рассмотрим опять эксперимент со свободным расширением. Предположим, что в начальный момент газ заполняет весь сосуд. Но потом (что совершенно неправдоподобно) он внезапно сжимает­ ся. Или, что более впечатляюще, пусть все, что мы наблюдаем, внезапно обратит свое направление движения, подобно обратному движению киноленты. Люди будут молодеть, земля станет вра­ щаться в противоположную сторону. Тогда, конечно, время поте­ чет в обратном направлении. Но этого нет!

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

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

Мы показали, что необратимое уравнение диффузии дает дисперсионное соотношение со = —Dk2 для прямого направления)*

*) При температуре и давлении, когда газовая фаза термодинамически устойчива.

5.4. Релаксация к равновесному состоянию

времени и со = A-Dk? — для обратного. В более общем случае линеаризация обратимого закона около равновесного состояния должна включать как возрастающие, так и затухающие решения, так что оба движения возможны для прямого течения времени. Если существуют только возрастающие решения, то, несомненно, должен быть возможен эксперимент, устанавливающий различие между двумя направлениями течения времени.

Задача 5.10. а) Рассмотреть дисперсионное соотношение для уравнений Максвелла в вакууме и доказать обратимость этих уравнений.

б) Подобным же образом исследовать обратимость уравнений

Н_

dJ

=0,

dt

дх

 

/ = аЕ,

дЕ _ q дх е0

В чем заключается основное различие между этими уравнениями и электромагнитным уравнением в вакууме?

б) Равновесные распределения. Микросостояния и макросостояния. Априорный подход

Уравнение диффузии

дп д2п

(5.125)

dt дх2

и ему подобные необратимы и, следовательно, находятся в про­ тиворечии с фундаментальными законами механики, которые являются обратимыми. Однако эксперимент подтверждает спра­ ведливость уравнения (5.125). А раз оно всегда удовлетворяется (при условии, что имеется достаточное число частиц, чтобы можно было ввести плотность п), то, по-видимому, закон диффузии (на­ зываемый иногда законом Фика) является верным. Но если этот закон справедлив, то существует различие между временем, теку­ щим в прямом и обратном направлениях (стрелка времени).

Экспериментальное подтверждение уравнения (5.125), каза­ лось бы, ставит под сомнение справедливость законов механики (которые, по-видимому, также всегда верны). Как согласовать эти два факта?

Хорошего «примирения» можно достичь, если ввести в тео­ рию вероятностные законы. Это так называемый априорный под­ ход — вывод следствий из определений или принципов, рассматри­ ваемых как самоочевидные. Введем прежде всего постулат равных априорных вероятностей, согласно которому все состояния на энергетической поверхности равновероятны.

304 Гл. V. Решение и дальнейшие свойства уравнения Больцмана

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

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

Этим идеям можно дать количественную основу, формулируя их через микросостояния и макросостояния. Микросостояние изо­ бражается вектором в Г-пространстве. Макросостояние опреде­ ляется через ц-пространство. Это фазовое пространство для одной частицы (ц —для молекулы; Г — для газа). Хотя координаты

и импульсы для N частиц системы изображаются одной точкой

в67Ѵ-мерном Г-пространстве, в 6-мерном ц-пространстве они будут представлены множеством из N точек (отдельной точкой изображаются координаты и импульсы каждой частицы).

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

Когда определено микросостояние, то известны координаты и импульсы всех частиц. Каждой частице присваивается индекс. То есть мы должны считать, что частица 1 находится в (qt, pt), частица 2 —в (q2, р2) и т. д. для того, чтобы определить микро­ состояние.

С другой стороны, чтобы определить макросостояние системы,

мы просто

говорим,

что

пх частиц находятся в ячейке 1,

п2

в ячейке 2 и т. д. Какие частицы из всей совокупности N частиц

находятся

в ячейке

1,

не существенно. Макросостояние

несет

в себе много меньше информации об объекте, чем микросостояние. Чтобы задать микросостояние системы, нужно знать всю ее струк­ туру. Чтобы знать макросостояние, необходимо знать лишь сово­ купность чисел заполнения.

Можно проиллюстрировать различие между этими двумя уров­ нями описания на примере системы из двух частиц. Система тако­ ва, что каждая частица может быть только в одном из четырех возможных состояний. Обозначая частицы через а и Ъ, можно

5.4. Релаксация к равновесному состоянию

 

305

представить микросостояние в виде диаграммы I д

I

I Ь

\

Так как таблица L J — I— 1— I представляет

состояния

одной

частицы, то она изображает ц-пространство. Одним из возможных

делений ц-пространства является 1 , 1 , 1 , т. е. разбиение его

на две ячейки. Тогда типичные макросостояния задаются парами чисел заполнения, например (1, 1); (2, 0); (0, 2).

Макросостояние (1, 1) представлено графически на рис. 5.8. Как много микросостояний соответствует этому одному макро-

f

1 • I а

I

Р и с . 5.8. Макросостояние,

соответствующее последовательности чисел

 

заполнения

(1,1).

состоянию? Мы изобразили эти микросостояния на рис. 5.9. Вось­ ми микросостояниям соответствует одно и то же макросостояние. Таким образом, имеется восемь путей реализации одного макро­ состояния.

С точки зрения этих понятий постулат равных априорных вероятностей означает, что все микросостояния изолированной системы, находящейся в равновесии, равновероятны х).

I

а I

I b I

I

1

I

I 0

I

 

1

I

а I____і____I

ь I

I

ь '____ I___ I

а

I

I

1 Д I

Ъ I

I

1

I

ь I а

1

 

I

I

I _ £ j ____ L-A..J

I

I

Ь 1

1

д I

Ри с . 5.9. Микросостояния, соответствующие макросостоянию

| • | е |

Каким бы, однако, всеобъемлющим ни было знание микро­ состояния системы, оно гораздо детальнее того, которое тре-

1) Результирующая сила для изолированной системы равна нулю. Следовательно, все состояния изолированной системы равновесные. Сила, вызывающая свободное расширение, не относится к типу F — та.

2 0 - 0 1 2 4 3

306 Гл. V. Решение и дальнейшие свойства уравнения Больцмана

буется в макроскопической физике, где мы чаще имеем дело с макросостояниями. Например, распределение энергии п (Е ) dE, задающее число частиц с энергией dE около Е , характеризует макросостояние. Одночастичное распределение по скоростям I d3^, задающее число частиц со скоростями в элементе d3^ около

I, также определяет макросостояние.

В рассмотренной элементарной задаче с двумя частицами, каждая из которых может находиться в одном из четырех упомя­ нутых состояний, мы нашли, что существует восемь микросостоя­

ний, соответствующих макросостоянию 1 * 1 * 1 . Анало­

гично можно установить, что существуют два микросостояния,

дающие одно макросостояние 1 ** 1 і , и два микросостоя­

ния, реализующие одно макросостояние I

I ** I . Ддярас-

смотренной системы эти три макросостояния

являются единст­

венно возможными.

 

Итак, в целом имеется двенадцать микросостояний, и все они равновероятны. Это означает, что в любой заданной (достаточно многочисленной) совокупности производимых измерений систе­ мы все микросостояния встречаются одинаково часто. При доста­ точно большом числе экспериментов 4/6 измерений дадут макро­

состояние

) •

1 •

I ,

1/6 измерений придется на макросостоя­

ние 1 * * 1

I

и

1/6

— на (

I I. . Этот простой пример

иллюстрирует два момента. Во-впервых, как видно из рис. 5.9, одному макросостоянию соответствует много микросостояний. Во-вторых, вероятность макросостояния (называемая иногда тер­ модинамической вероятностью) пропорциональна полному числу различных микросостояний, реализующих данное макросостояние.

Так как концепция макросостояния основана на разбиении ц-пространства на ячейки с некоторыми шестимерными объемами, то уместно вывести правило для получения вероятностей макро­ состояний через объемы Г-пространства, которое является есте­ ственной системой координат для микросостояния.

Поскольку любое заданное макросостояние определяется через объем в р-пространстве, относящиеся к нему микросостояния

Соседние файлы в папке книги из ГПНТБ