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

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

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

3.2. Анализ Боголюбова ББКГИ-уравнений

127

Кроме того, / 8 зависит от времени неявно, через функцию/!. Примером такой зависи­ мости может служить функция f s для s ста­ тистически независимых частиц:

 

 

 

 

 

 

/. = ,П /і (*)•

 

 

(3.17)

 

 

 

 

 

 

 

 

 

 

h=l

 

 

 

 

 

 

 

Наконец, на гидродинамической стадии

 

 

 

 

 

 

 

Т о все функции распределения f s

 

 

 

зависят только

от

п, и и Т, т.

е. первых

 

 

 

пяти моментов функции /4.

изображена

 

 

 

 

На

рис.

3.3 схематически

 

 

 

последовательность только что рассмотрен­

 

 

 

ных случаев (когда при t = 0 только одна

 

 

 

частица из всего стационарного газа при­

 

 

 

ходит в движение).

 

 

 

 

 

 

 

 

 

При изменении состояния системы от

 

 

 

начального

хаотического до конечного хо­

 

 

 

рошо организованного равновесного со­

 

 

 

стояния нам требуются резко различаю­

 

 

 

щиеся уровни описания системы. Вначале

 

 

 

необходимо

знать

не меньше

чем N-ча-

 

 

 

стичную функцию распределения. В конеч­

 

 

 

ной, равновесной, стадии достаточно знать

 

 

 

термодинамические

 

переменные,

 

дающие

 

 

 

несравненно менее подробное описание си­

 

 

 

стемы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Анализ Боголюбова иерархии ББКГИ-

 

 

 

уравнений касается

кинетической стадии

 

 

 

и включает следующее ограничение:

 

 

 

 

 

 

lim

(

^

Ь

Ч

- г Ь

т

.

(3.18)

Р и с .

3.3. а — началь­

 

 

 

N-

 

 

 

 

 

Ѵ-*оо

 

 

 

 

 

 

 

 

ная стадия: 0<1 < т0, /_ѵ

где і/ѵ является конечным малым пара­

определяет состояние га­

метром.

(Соответствующим безразмерным

за; б кинетическая ста­

малым параметром будет г30/ѵ, где г0—ра­

дия:

т0 < t < t 0.

Части­

цы приобретают

некото­

диус взаимодействия.) В этом

предель­

рый

элемент сходства,.

ном случае (см. задачу 2.27) уравнение БИ г

fs = fs [/il; 6 — гидроди­

имеет

 

вид

 

 

 

 

 

 

 

 

намическая стадия: t0<

 

 

 

 

 

 

 

 

 

<1 t <_ Т0. Состояние га­

dFt

,

 

f

 

1

1

 

d(l + l)Qul+1Fi+ = 0.

за близко к равновесно­

 

 

 

му и задается через функ­

dt

 

■LiF I

 

 

■3

 

 

ции и, и,

Т.

 

 

 

 

 

 

 

i= l

 

 

=

 

(3.19)

 

 

 

Напомним,

 

что Fi =

Vlfi и і

(

Xi, Pi). Функциональная зави

 

=

х

 

 

 

симость между Fs и Fj имеет следующий вид:

 

 

(3.20)

 

 

 

 

 

 

 

 

 

Fs (1,

 

s: Fi).

 

 

128 Гл. II I . Формальное развитие теории кинетических уравнений

Уравнения (3.19) для I = 1, . . s запишутся так:

dt

— — - ! ^ + — ( d 2 ë 12F2,

(3.21)

т ds.i

V

J

и f

 

9FS

 

s

 

 

 

LSFs+ ~

2

 

®г. s + l - f 1S + I -

(3.22)

dt

 

i=l

Используя (3.20), производную по времени от Fs можно предста­ вить следующим образом:

д І \

/ б F4

dF,

(3.23)

I t ~

XöFi

ИГ

 

где скобками обозначено внутреннее произведение по всем функ­ циям, a 8FS/8F! — функциональная производная от Fs по F^ В простейшем случае, когда Fs является алгебраической функ­ цией от Flf эта производная запишется в виде

9FS _

dFg

dF1 (I)

dt

Гл

dFi (l)

dt

 

1=1

 

 

Исходными уравнениями анализа Боголюбова являются послед­ ние три уравнения (3.21) — (3.23). Формальное разложение, кото­ рое будет использовано для получения замкнутой системы уравне­ ний, задается равенством

Fs( 1,

JP1) = ^ 0) + 4 ^ 1,+ . . . ,

(3.24)

F! = Fi (Fi не разлагается, поскольку она является «независимой переменной»).

Подставляя (3.24) в (3.21), получим

dFj

»L ~ Ч

 

( d2ëiz

Г/1'

. P<i>

4-

dt

1

т

öxj

 

V

J

л

L 2

1 2

 

Н(0, + — А (

 

 

 

 

 

 

 

(3.25)

 

'

V

 

 

 

 

 

 

 

 

Это соотношение

служит

для определения коэффициентов ИФ;

 

 

Ат

 

 

 

Рі

d F i

 

■LiF1>

 

 

 

 

 

 

т

дхі

 

 

 

 

И(1) =

+

j

â2QaF ? \

 

(3.26)

A(l)= + j d2ëlzF t i)-

3.2. Анализ Боголюбова ББНГИ-уравнений

129

Если мы подставим разложение (3.25) для dFJdt в функциональ­ ное внутреннее произведение (3.23), то в результате будем иметь

9F7,

/ bFs

,

А 1' т А " ’+ - - ) =

dt

\6 F i

/ бö[Fj<»Р

+ (i/v) FSP+ . . . ]

А™ +±-А™ + . . . у . (3.27)

 

бFi

 

,

 

 

 

 

Путем перемножения членов, получим

* • > + * { < » ? . ~ > + ф . * • »

(3.28)

Это разложение можно записать в следующем виде:

dFs

: F>(0>F '0>+

[DWF ^ +

F>(1)F '0>]

(3.29)

dt

 

'

V

 

 

 

где операторы

и

Z)ll)

определяются равенствами:

>> =

/

 

5

f)(l) =

/_б

(3.30)

Д"

W

i ’

/ ’

_

W i

 

>■

Вернемся теперь к оставшимся уравнениям ББКГИ (3.22) и ис­ пользуем в них разложение (3.24). В результате получим

^ = - L , [ f i » + 4 - W + . . . ] +

2

j d ( s + i ) é {. e+1[/'i¥1+ 4 ^ i + - - - ] .

(3.31.)

г =

1

 

Приравнивая коэффициенты при одинаковых степенях к-1 в раз­ ложениях (3.29) и (3.31), придем к замкнутой системе функци­ ональных уравнений для переменных F® [1, . . ., s; FJ:

b imF(am= — LsF™,

(3.32а)

 

 

S

 

 

D™F? + D™F? = -

L,F?' +

2 j

d (s + 1) Ѳ*. S+1F'°>,

 

 

I

І = 1

 

 

n

s

 

 

2 b mF ? - l) = - LsF™ + 2 \ d (s + 1) Ѳ,. s + t F ^ .

(3.326)

J=0

i=l J

 

 

Последовательный ход

решения

этих

уравнений представлен

на рис. 3.4.

 

 

 

 

Резюмируя изложенное, отметим три этапа анализа Боголю­

бова: 1) предполагается,

что FS= FS [F^; 2) находится уравнение

для Fs [F4]; 3) построение завершается выводом двух разложений для dFJdt. Ими являются уравнения (3.29) и (3.31).

9-01243

130 Гл. II I . Формальное развитие теории кинетических уравнений

В данном формализме мы вывели замкнутую систему уравне­ ний для функций Fs, но две задачи остались нерешенными. Это, во-первых, получение конкретных выражений для операторов

D (t) и, во-вторых, метод вычисления Fi.

Именно на этих двух проблемах мы сосредоточим сейчас наше внимание. В частности, дадим стандартный метод, с помощью

Р и с . 3.4. Метод решения системы уравнений (3.32).

которого получается уравнение для Fi, справедливое до членов порядка ѵ~г. Для этого используем замкнутую систему уравнений, получаемую из (3.25) и (3.32а),

9 F I

I

P i

d F j

1

j d2è12*T,

(3.25')

dt

'

т

дХі

V

 

 

DmF ^ =

 

 

(3.32a')

Первый шаг состоит в том, чтобы решить второе из этих уравне­ ний для

F™ =F™( 1, 2; Fi)

при некоторых определенных «граничных условиях». Это дает явную функциональную зависимость F20>от Ft при соответствую­ щим образом выбранном начале отсчета времени. Последнее озна­ чает, что система частиц (со взаимодействиями, препятствующими существованию связанных состояний) становится некоррелиро­ ванной, когда протечет достаточно длительный интервал времени в обратном направлении. Уравнение (3.32а) решается именно при этом условии. Такое условие поставлено в связи с его дальней­ шим приложением к процессам столкновения и означает, что до столкновения частицы не коррелированъ!.

На втором шаге полученное таким образом решение для F20> подставляется в уравнение (3.25') и выводится соответствующее

выражение для Ѳ12ІД0>, что приводит к характерной форме боголюбовского обобщения уравнения Больцмана.

3.2. Анализ Боголюбова БВКГИ-уравнений

131

Для решения этих задач вводится оператор сдвига

 

= exp [tLi],

(3.33)

где Li — ^-частичный оператор Лиувилля (cp. с Â;,

который

задается уравнением (2.30)). Определение L t вытекает из уравне­ ния Лиувилля для изолированной системы из I частиц:

Ц г = \Ни Fü ^ - U F l

(3.34)

Решение этого уравнения имеет вид

Ft ( 1,

I; 0 = А(Ѵ г(1,

J; 0).

(3.35)

Оператор ДІ?Т действует на фазовые переменные Fi- Он «сдвигает» их обратно по времени на интервал т по отношению к значениям, определяемым Z-частичным гамильтонианом і / г. Для единичной свободной частицы мы получим

 - \ F 1 (х, р, 0 = ^ і ( х —-|;Т , р, tj .

(3.36)

В более общем случае (возвращаясь к обозначению фазовой переменной через z) для произвольной фазовой функции ф мы имеем

Â(t4>(zi > . . . , Z b г) = ф(Д(тг)2і, . . . , Д ^ 2 г; t).

(3.37)

Из (3.35) следует, что Ft постоянна вдоль траекторий системы в Г- пространстве І частиц. Ниже приведено несколько других свойств

оператора А:

 

До=1,

(3.38)

 

д 11 ^ т 2 = Ä ( T l + T 2)>

(3.39)

ад-т

— LA_t = — Д-XL.

(3.40)

дх

 

 

Последнее соотношение справедливо для операторов, не завися­ щих от времени. Решением неоднородного уравнения Лиувилля

дР (z, t)

LF (z, t) = P (z, t)

(3.41)

dt

 

 

является функция

t

 

 

(3.42)

F (z, t) = Â_tF (z, 0) + A_, j' dt'hfP (z, t')

(здесь использовалось свойство (3.40)).

Упомянутое выше граничное условие для «-частичной функции распределения как функционала от Fi легко задать через опера­

9*

132 Гл. III . Формальное развитие теории кинетических уравнений

торы Â. Как мы помним, это условие состоит в том, что система s взаимодействующих частиц в далеком прошлом была некор­

релированной, т. е.

,

 

lim Â%FS[1,

.. ., *; F, (1, t)] = lim Д«

ft Fi (h <)• (3-43)

Т-fc o

T—У оо

1 = 1

Поскольку данное условие справедливо для произвольной функ­

ции Ft, оно будет справедливо и для F't = lim At^ ^ i (1, t).

T-»- со

Подставляя

эту функцию в предшествующее равенство, получим

 

 

S

lim Д(1\F S[1,

.. ., s; А(Л Л (1, t)) = lim

П А(Л Ft (l, t). (3.44)

Т-Юо

T—>oo

1 = 1

Такая форма граничных условий используется при решении функ­ ционального уравнения (3.32а). Явный вид правой части (ПЧ) равенства (3.44) легко получить, используя (3.36). Это дает

ПЧ (3.44) = lim Â(_?T ft F A ( x t + ^ - x ) , рЛ =

T-*oo

г = 1

L '

H l

I

J

— lim Ц Fi [ Â ^ (хг+ - ^ - т ) , А ^рЛ =

T-»oo

L

'

i n

!

J

= II Л [ х Л , р Л ]. (3 .45)

г=і

Переменная pls> — это постоянный импульс, который имеет 1-я частица, взаимодействующая с совокупностью s частиц, если про­ следить за этой частицей из состояния р г в далекое прошлое. Переменная x';S) — это координата, которую имела бы l-я частица, если бы она двигалась с импульсом p)S) при начальной координате

Ali'ooX;. Такая интерпретация обусловлена той операцией, которая определена последним равенством:

X;S) = lim Â(_s)x ( Щ+ ~ т ) = Н т ( А(Д хг +

х ) .

Как мы увидим, эта интерпретация весьма полезна при нахож­ дении явного выражения для в уравнении (3.25').

Перейдем теперь к интегрированию уравнения (3.32а):

DmF ^ = —L»F?\

 

 

Согласно (3.26) и (3.30),

операция Dl )F f'

определяется

следую­

щим образом:

 

 

 

DtmF™ =

Л<°>(Fi) =

( - LiFi).

(3.46)

3.2. Анализ Боголюбова Б Б К ГИ-уравнений

133

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

любой функции F1: оно будет справедливо также и для AliiFt. Таким образом,

DWF?' [1, . . s; А Э Д ]= -

р — Г- U

• (3.47)

о

[Л- т^і] L

J

Вспоминая свойство (3.40), мы можем записать это уравнение так:

/)«»/?'<».

бFT dbll\Fi

(3.48)

 

' ÖIÂ-tFj]

что дает нам желаемый вид функциональной операции:

Z)(0)F ‘0,[1, . . . , 8; ДЗД] -

дх F T [ 1,

ATXF,

(3.49)

Подставляя это выражение в (3.32а), получим

 

 

F T [1, . . . ,

s; Д -’т-Fі] =

- L SFT [1, . ..,

s; Д З Д .

(3-50)

Рассмотрим пробное решение

 

 

 

F T [1

Д-т^і] =

xF ? (1.........

Л ).

(3.51)

Подстановка его в (3.50) дает

± № xF T (1, . .., в; Fd = - Ls№ xF T (1, - • •, s; F,). (3.51a)

Если мы вспомним (3.40) и заметим, что FT в таком виде, как она записана в уравнении (3.51а), не зависит от т, то увидим, что это уравнение является тождеством. Отсюда следует, что решением уравнения (3.50), а следовательно, и уравнения (3.32а) будет про­ извольный функционал от функционала, заданного соотношением

(3.51).

Другой формой решения может быть следующая:

F T (1, .. . , s; Fi) = Д( Ѵі® (1, . ■■, s; Дт’Еі).

(3.52)

Это решение пригодно для всех интервалов т, так что

F T ( 1, . .., s; ^ ) = И т A<Ѵі°>(1, ■■■, s; Д^ ) (3.53)

T->co

также является решением. Используя предшествующее граничное условие, а также (3.44) и (3.45), получим

F T (U . . . ,s ; E i ) = [I iM x(A Pis)).

(3.54)

i=i

 

134 Гл. II I . Формальное развитие теории кинетических уравнений

Это решение удовлетворяет уравнению (3.32а) и граничному усло­ вию (3.44). Таким образом, первая часть задачи о получении одно­ го уравнения для Ft из уравнений (3.25') и (3.32а') завершена. Из (3.54) видно, что решение низшего порядка Р$т в анализе Бого­ любова обнаруживает s-частичную корреляцию между частицами (через переменные x(S), p(S)).

Следующий шаг состоит в том, чтобы получить в уравнении

(3.25') соответствующее выражение для 6>i2F'0). С этой

целью

напомним сначала соотношение (2.155) и последующие:

 

Ь2 = К 2- & 12 )

 

к 2 = Рі

д

Р2

д

 

Зх-і

т

д %2

 

 

 

Отсюда вытекает, что

 

 

 

 

e 12F<°> = K2F‘°> -L 2F ‘0’.

(3.55)

Чтобы получить надлежащее выражение для L2F'0), вернемся

опять к уравнениям (3.32а') и (3.26):

 

 

- L 2F™ = D"'F™ = ( ^

- , Л«”) ,

(3.56)

/ісо) _____Рг_

dFi .

ч

 

Ат ' дхі

Выполняя функциональное суммирование, соответствующее внут­ реннему произведению ( , ), и вспоминая решение F '0>, которое

задается посредством (3.54), получим явный вид для L 2Fl2m:

F2F r

р'2>

Р22)] * д а Л К 2), p H -

 

(3.57)

+ -

Fi [х)2), р l^ — F, [x'2\ p H -

Наконец, чтобы получить нормальную форму обобщенного

уравнения Больцмана, произведем в операторе К 2 замену перемен­ ных Х[ и х2 на Хі и г = х2 — х4, что дает

К = {Р2~ Рі ) .JL4_.P1__

К2, Б ■ к 2.

(3.58)

2 \ т ) д г ' m öxi

 

 

Используя выражение (3.57) для L2F'0) и (3.58) для K 2F f \ под­ ставим (3.55) в (3.25') и придем к обобщенному уравнению Больц­ мана в форме Боголюбова:

dFj

Pi

djj

V j d2Kz, bF ‘0) + І - j d2 (K2- L 2) F ‘0). (3.59)

dt

m

Ö X j

В случае, когда F ‘0> однородна по всей области взаимодействия, второй интеграл обращается в нуль (К2 и Ь2 оба содержат про­

3.2. Анализ Боголюбова ББКГИ-уравнений

135

странственные производные), в то время как К 2>Б переходит

вобычный больцмановский оператор столкновений. Чтобы это показать, в соотношении

J К 2. бF T = j dx2dp2V - ^ F {F1[р)2>] F, [р«’]} (3.60)

(где аргумент х в F\ выпадает, поскольку F^ однородна по всей области столкновения) выполним интегрирование по х2 в цилин­ дрических координатах с осью цилиндра г), параллельной отно­ сительной скорости

 

 

 

 

V = -?2~ Pi

; т] = - ^ .

 

(3.61)

 

 

 

 

 

 

т

1

V

 

N

Радиус цилиндра Ъи азимут ф лежат в плоскости,

нормальной к У.

Отсюда следует,

что

 

 

 

 

 

 

 

 

 

 

dx2= Ьсіф db dp,

 

 

 

 

 

 

 

 

=

 

 

 

(3-62)

Интегрируя сначала по т], получим

 

 

 

 

+°°

 

 

 

 

 

 

 

 

 

 

j

drj

{F±[рП Ft [р™]} = Fi [pH Fi [pn + l -

(3.63)

 

—oo

 

 

 

 

 

 

 

 

Области после столкновения соответствует ц =

+оо (рис. 3.5).

Если из этой области (ц =

+ °°) проследить в бесконечно далекое

прошлое, то каковы будут зна­

 

х2

 

 

чения рі и р2

в случае парного

 

 

 

взаимодействия

частиц?

Ими

 

 

 

 

должны быть те значения им­

 

 

 

 

пульсов, которые перейдут

в р4

 

 

 

 

и р2 после столкновения, т. е.

 

 

 

 

импульсы р(

и

р' инверсивного

 

 

 

 

столкновения.

(Эти понятия бу­

 

 

 

 

дут

подробно

обсуждаться

в

 

 

 

 

гл. IV.)

до столкновения со­

 

 

 

 

Области

 

 

 

 

ответствует т] =

—оо (рис. 3.5).

 

 

 

 

Каковы будут значения рі и р2

 

 

 

 

при парном взаимодействии ча­

 

 

 

 

стиц,

если из этой области (ц =

 

 

 

 

= —оо) проследить в бесконеч­

 

 

 

 

ное прошлое, которое представ­

Р и с .

3.5. Связь между знаком р и

ляет

собой

состояние «до

до

 

областью столкновения.

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

стиц? В этой области две частицы движутся независимо одна от

136 Гл. I l l . Формальное развитие теории кинетических уравнений

другой как свободные, поэтому их импульсы равны рі и р2. Сле­ довательно,

Fi feO Fl [рП |t~ = F, (p() F, (pj) - F1(Pl) Fl (Pa). (3.64)

Подставляя это соотношение в (3.63), а затем в (3.59) (и удержи­

вая только однородный вклад К 2,в)> получим обычное уравнение Больцмана:

^ +

^ I

р Э -

 

- Л ( Р і) ^ і (Р2)]-

(3.65)

Уравнение Больцмана подробно рассматривается в гл. IV, где представлены три других его вывода. Более прямой из них не опирается на последовательность ББКГИ-уравнений. Но весь­ ма примечательно, что анализ Боголюбова дает независимый (т. е. опирающийся только на уравнение Лиувилля) метод, кото­ рый позволяет получить больцмановский интеграл столкновений

вявном виде. В данный момент важно, что такой метод существует.

Вглаве IV будет дана всесторонняя трактовка уравнения Больц­ мана, которая обеспечит более глубокое его понимание.

3.3. Групповые разложения

Совсем недавно были предложены другие методы разложения для исследования ББКГИ-уравнений. Сюда относятся работы Ростокера и Розенблюта (1960), Гернсея (1960), Ленарда (1960), Балеску (1960), Сэндри (1963) и Дюпри (1961). Однако эти работы

впротивоположность предыдущим больше связаны с выводом определенного кинетического уравнения, чем с получением реше­ ния для Ѵ-частичной функции распределения. (Кинетическим называется любое уравнение в переменных (t, |, х) для неизвест­ ной функции fi.) Тем не менее все эти методы объединяет то, что

вкаждом из них решение ищется в виде некоторого разложения. В частности, метод Боголюбова связан с разложением по обратным степеням параметра ѵ — удельного объема.

Существуют другие способы разложения, непосредственно опирающиеся на свойства Fz, F3, . . ., F N в предельном случае, когда k-1-> 0 . Что это за свойства? Очевидно, что, когда объем, приходящийся на одну частицу, становится большим, а следова­ тельно, и расстояние между частицами также возрастает, тогда фаза какой-либо частицы становится менее связанной с фазой любой другой частицы. Это утверждение о независимости состоя­ ний, записанное через функцию распределения Fs, выглядит так:

^ = П Л ( 0 -

(3.66)

!=Т

 

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