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

книги из ГПНТБ / Эккер, Г. Теория полностью ионизированной плазмы

.pdf
Скачиваний:
4
Добавлен:
22.10.2023
Размер:
12.01 Mб
Скачать

§ 1. Точные уравнения для плотн. распр. однокомп. системы 123

зультате суммирования получим

dF<-hi

г

р

d F ^ )

~ t

 

d r

X F(h) 2 ^ ( г ' , p', * к } = 0 .

( 1. 10)

Здесь учтены слагаемые, первоначально исключенные [поэтому в (1.9) использовался знак суммы со штрихом]. Это оправдано, поскольку такие слагаемые дают пренеб­ режимо малый вклад по сравнению с другими N — 1 чле­ нами суммы. Уравнения (1.10) также представляют собой систему совместных дифференциальных уравнений. Одна­ ко порядок этой системы уравнений определяется теперь числом различных сортов частиц.

Уравнения (1.10) можно переписать в виде

др

( п и

где

Е“ - j*' j S Р'. 0 «.(тгц^гг) <&'• <Ы2)

Для величины Е выполняется равенство

дг Е = 4 п 2 J ^ (V)(r >Р'> * )М р',

(1.13)

если вспомнить, что

А ( т Г ^ Т ) = - 4лб(г - г')-

(1.14)

Уравнения (1.11) и (1.12) называются уравнениями Климонтовича [1] для истинной плотности распределения частиц рассматриваемой системы. Эти уравнения являют­ ся точными и свободными от приближений.

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

124 Гл. 2. Неравновесные состояния. Общее описание

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

§ 2. РАСПРЕДЕЛЕНИЯ, УСРЕДНЕННЫЕ ПО АНСАМБЛЮ ГИББСА

2.1. Вывод цепочки уравнений Боголюбова — Борна — Грина — Кирквуда — Ивона (ББГКИ ) из уравнений Кли­ монтовича

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

/ jv ( ri> • • •> r N , Pi, . . ., P jv-, t).

(2 . 1 )

Плотность вероятности i-й частицы в данной точке (*г, *р) р-пространства определяется следующим образом *):

 

 

U Or. fP;

t) =

j

fNFi dti ... dpN.

(2.2)

В более общем виде плотность вероятности найти набор

частиц

(1,

. . ., s)

в

точках (хг, . . ., 5г, *р, . .

., 5р)

дается

частной

функцией

 

распределения

порядка

s:

U (хг, . . . / г ,

Д>,

. . . ,

8р;

t) =

j

fNFt ... FsdTi

. .. dpN . (2.3)

Обратим внимание на разницу в написании индексов. Индекс внизу справа (г,-, рг, г;-, р^) характеризует положение частиц

(точек) в Г-пространстве. Индекс вверху слева (»г, {р, ir, ?'р) характеризует точку наблюдения в д-пространстве, в которой мы рассматриваем эффекты от частиц г и

§ 2. Распределения, усредненные по ансамблю Гиббса

125

Здесь мы использовали небольшое видоизменение соот­ ношения (1.1):

Fi (*г, *р) = 6(*г — гг (<)) 6(*р — рг (t)). (2.4)

Чтобы получить цепочку уравнений ББГКИ, которым удовлетворяет частная функция распределения /8, необ­ ходимо умножить уравнения (1.9) на плотность вероятно­ сти (2.1) в фазовом пространстве, а затем проинтегриро­ вать по Г-пространству. Тогда

 

X ^ S

^

Д * ! . . . dpN = 0.

(2.5)

Плотность

вероятности

в

фазовом

пространстве

f N (rt, . . ., rN,

р4, . . р jvJ

t) не зависит

от координат

{г,*р, и поэтому выполнить

интегрирование во

втором

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

ношение:

 

 

 

 

 

 

 

 

 

f

(

\

_(

d/jv^i

\

_р

I

dfif

\

__

 

lN \ dt

 

\

dt .

 

/ir> ip

4

dt

}it

ip

 

 

 

(

djj^Fi \

 

d f ^

 

/ dftfFi

\

cy

 

 

= v—^ — )ir> ip—jPi~dT == v—^ — ; ir> ip *

 

Используя соотношение (2.6) и равенство

 

 

/

SfNFj

\

_ /

dfNFj

\

 

 

 

 

 

 

'

dt

Ч , *p

'

dt

' V V r . n.

 

 

 

 

 

+ S ( ^ - 4 - + i / T s r ) v v ,ff^<). <2-7>

получаем

 

 

 

 

{ (

dF‘

\

I

d^Fj

\

,

, N \

dt

) i I<ip

{

dt

) i r

ip>r. , P,

+ 2

( 2.8)

126 Гл. й. Неравновесные состояния. Общее описание

Отсюда находим

 

 

1 М - ж - )<,<„**••■ ip » = (

d/i(*r. ip; О

(2.9)

 

d t

)*Г, *р

 

 

Здесь учтено, что второй член в (2.8) исчезает при интегри­ ровании по частям и применении уравнений движения.

Окончательно, исходя из определения частных функций распределения /8, находим уравнение для средних распре­ делений

дк

д / 1

 

 

 

d t

т

д к

 

 

 

 

 

— f dV f d jp'

-----г —т— )

 

 

J

J

Vd lr

I *r— ?r' I '

 

 

 

S ' / 2 (*r,

V , {P ,

V ; t ) e j = 0 . ( 2. 10)

 

 

d l p

].

 

 

Это и есть первое уравнение цепочки зацепляющихся урав­ нений для функций распределения. Следующее уравне­ ние выражает / 2 через / 3 и т. д. Прежде чем вывести такую цепочку уравнений, упростим задачу путем введения сле­ дующей частной модели.

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

2.Будем считать, что плотность вероятности в фазовом пространстве / N является симметричной. Поскольку урав­ нение Лиувилля симметрично относительно координат, такое предположение всегда справедливо, если начальное распределение для f N симметрично относительно коорди­ нат частиц. Однако следует заметить, что в целом ряде интересных задач это предположение не выполняется. Например, требование симметрии нарушается в проблеме пробной частицы. С другой стороны, пользуясь развитым ниже формализмом, нетрудно написать соотношения так­

же и для несимметричных случаев.

Используя эти упрощающие задачу допущения, вели­ чину еге} можно заменить на е2. Тогда уравнение (2.10)

§

2.

Распределения, усредненные по ансамблю Гиббса

127

запишется

в виде

 

 

 

 

 

 

dfi .

dfi

 

 

 

 

 

 

dt ~

т

g гг

 

 

 

 

 

 

 

е2 (N — 1)

[ d V

( ( —--------—

)•

 

 

 

v

'

J

J

v дч \ ir— It' I /

 

 

 

~

/

2(*r, V ,

*Р, V ; ^ V

= o.

(2.11)

 

 

о гр

 

 

 

 

 

 

Это уравнение является первым уравнением цепочки ББГКИ. Оно связывает частную одночастичную функцию

распределения

с двухчастичной функцией распреде­

ления / 2.

 

Чтобы решить уравнение (2.11), мы должны получить уравнение для двухчастичной функции распределения / 2. С этой целью рассмотрим два уравнения Климонтовича для частиц сорта г и А: в виде

dFt

dFj

dt

' т

QгГ

 

 

е2—^— ’ { f d}r'

f

( — .---:—

)

X

 

 

 

d*p

lJ

J

Vd *r I *r—it' I

/

 

 

 

X F t C r , *p) 2

 

V) } d V = 0

(2.12)

dFk

ftp

dFk

 

 

 

 

 

 

dt

m

ght

 

 

 

 

 

 

 

 

- « ■ J - H i v f ( 4 - _ _ v _ ) x

 

 

 

a ftp

l J

J \

aftr i fcr —Jr' |

*

 

 

 

X Fh (kr,

ftP) 2 Fj (V, У )}

d V = 0.

(2.13)

Умножив (2.12) на Fhf N, a (2.13) на Ftf N и сложив оба уравнения, проинтегрируем полученный результат по Г- пространству, используя при этом соотношение (2.3).

128 Гл. 2. Неравновесные состояния. Общее описание

Тогда получим

dft

df2

ftP

df2

at

т

gir

m

ght

- e2 i dV j

5

J - F ^ ^ F j C v ', V ) * , . . . d p w-

5 ip

. .

9 ftp

2 F) (V , V ) drt ... dPiV = 0.

(2.14)

}фк

 

Рассмотрим два последних члена в левой части этого уравнения. Поскольку функция f N в фазовом простран­ стве зависит только от переменных с нижними индексами, расположенными справа, f N можно вносить под знак дифференцирования. Оба члена содержат слагаемые двух типов. Одни содержат / ’-функции, индексы которых отли­ чаются друг от друга и довольно просто интегрируются по Г-пространству, давая трехчастичную функцию рас­ пределения / 3 (’г, fer, V , *р, ftp, V ; t). Другие слагаемые содержат две / ’-функции с- совпадающими индексами:

5

\ *v !*■••• ^ (i V t? - v i

)•

- j r - Fh (hi,

йр) Fh(V , У ) Fi (*г, *p) fN (n, ... ,

pw; t)

- * l * * l * V И-ИгргЬт)-

У Ft (V, У ) / * ( кг, ftP) x

a *p

X fN (ri, . . . , pN; t) drt . . . dp^.

(2.15)

В этом выражении интегрирование проводится сначала по координатам hr', йр' и гг', гр' соответственно, а затем по координатам Г-пространства. В результате получаем

—е2( ( —----

— ------

) —d f2 +

I \ 9 ij-

I i f — йг I

/ g in

+ ( - k i ^ ) - T r j A ■ <216>

 

§ 2. Распределения, усредненные по ансамблю Гиббса

129

Подставляя

эти

результаты

в

(2.14), мы

приходим

к уравнению

 

 

 

 

 

 

 

 

 

д/г

, У

df2 .

fep

df2

 

 

 

 

 

 

dt

‘ m

д 4

'

m

d hr

 

 

 

 

 

 

— e2

 

 

 

 

 

 

 

— — V

д/г

1

 

 

 

 

 

 

 

 

lr —hr|

!

d ftp

/

 

 

 

 

 

 

 

 

-1—

).

 

 

 

 

 

 

 

 

 

 

| ’r—3r' I

/

 

 

 

 

~

/

3Cr, ftr, v,

у

у

V; 0+

 

 

 

 

 

д гр

 

 

 

 

 

 

 

 

 

 

 

 

+ {— -------- )•

 

'

 

 

 

 

 

 

 

\

dhr | hr—V I

 

 

 

 

. - i -

U ({r, У V, У У V; 0} dV = 0.

 

(2.17)

 

д йр

 

 

 

 

 

 

•*

 

 

 

Это и есть второе уравнение цепочки ББГКИ, которое связывает двухчастичную функцию распределения / 2

счастной трехчастичной функцией распределения / 3. Легко видеть, что такая процедура может быть продол­

жена до получения общего вида цепочки уравнений:

d/s dt

у . *Р d/s

—1 "г а»Г i=l

с 2 у '

/

д______1___

\ J / .

^

V(Иг Иг— fer|

1

> dip

i, ft

 

1

 

 

 

 

i=1

 

 

 

/s+1( . . . ,

v ,

. . . , V ;

t) d у - о . (2 .18)

d »p

Для реальных физических задач обычно частные моле­ кулярные функции распределения /8 не представляют инте­ реса, так как в принципе невозможно различить между собой частицы одного и того же сорта. Значительно боль­ ший интерес представляет изучение общих функций распределения /<8>, которые дают среднюю плотность любых s частиц в наперед заданных положениях *). Общие§*

х) Для сравнения см. аналогичное рассуждение в начале

§ 4 гл. 1.

9 - 0 1 2 9 1

130 Гл. 2. Неравновесные состояния. Общее описание

функции распределения определяются следующим обра­ зом:

(2.19)

Вводя эти общие функции распределения /<S) в це­ почку уравнений (2.18), приходим к окончательному виду цепочки уравнений

а/<*>

у гр d f< V

____ g2

у

' I

_________ 1_______ \ , 5 / (

 

 

 

 

 

 

 

 

dfw

 

9t

Z i т ’ g i T

 

ZA \ g i t

| i r _ fe r | / . g i р

 

 

i= l

 

 

г, к

 

 

 

 

 

- e 2

[ d V

[ У

( - 4 ---------— )

 

 

 

J

 

J.

\

д гr

I г Г —Jr' I '

 

 

 

 

 

i= l

 

 

 

 

h

r

v

,

 

.... V; 0<*V= o.

(2. 20)

 

9 *p

 

 

 

 

 

 

 

В литературе принято использовать вместо частных функций распределения f s нормированные функции рас­

пределения /8, определеннее следующим образом:

/. = 77.- (2-21)

Для этих функций цепочка уравнений (2.18) принимает вид

d / s . у

d j s

„2 у

' /

__________1

 

\ .

J L _

^—1

го

д if

Zj

\ д ir lip — &r I

/

9 ip

i = l

 

 

i,

ft

l

I

 

F

•, iV —s f <2 v

f ^ v

2

--------— )

 

J

 

J

^

l S'г Иг—Jr'I /

 

 

 

 

1=1

 

 

5

/.+i ( . . . ,

v ,

V ; 0 = 0.

(2.22)

 

 

 

 

 

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

Литература

131

(2.20) и соответственно (2.22) неявно подразумевается пре­ дельный переход, при котором F - > o o h . / V - > - oo, h o N / V =

= const.

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

ЛИТЕРАТУРА

1. Климентович Ю. Л ., Статистическая теория неравновесных про­ цессов в плазме, изд-во МГУ, 1964.

Дополнительная литература

Chappell W. R., Journ. Math. Phys., 8, 298 (1967).

Chappell W. R., Microscopic kinetic theory, в книге Lectures in Theoretical Physics, Yol. IX, C, Kinetic Theory (ed. W. E. Brittin), Univ. of Colorado Press, Denver, Colorado, 1967.

Wu C. S., Plasma kinetic theory in the Klimontovich formalism, в книге Lectures in Theoretical Physics, Vol. IX, C, Kinetic Theory (ed. W. E. Brittin), Univ. of Colorado Press, Denver, Colorado, 1967.

9*

Глава 3

Неравновесные состояния кулоновской системы. Приближенное описание в отсутствие корреляций между частицами

§ 1. ПРИБЛИЖЕНИЕ ВЛАСОВА

1.1. Уравнение Власова

Рассмотрим первое уравнение цепочки (2.20), получен­ ное в гл. 2, и обобщим его на случай нескольких сортов частиц (р, v, . . .) *)

d/U) (Мт, М.р; t )

Вр 3/0) (^г, Мр; t )

d t

gV-t

( <Рт’ [ dyр' ( —------

^ ------

)•

J

J Р \ д»т | i * r - V | /

V

(ИГ) Vr'5Bp) vp';i)= 0. ( 1. 1)

д'

Определяя

корреляционную функцию

следующим

образом 2*):

t) =/<!> (иг , ^р ; t) / (1>(vr , Ч>; t) +

 

/ (2>(**г, vr, i*p, р;

 

+ ^ v ('4r , vr, |Xp , vp ;0 ,

(1.2)

находим из (1.1) уравнение для общей функции распре­

деления / а>

частиц сорта р:

g/Я)

dfО)

-вц<Е)-

а / < * >

dt

d

d

 

g^vC^r, v ,

*Р, V ;< ),

(1.3)

 

о

 

 

 

х)

Латинскими буквами обозначаются

различные

частицы,

а греческими буквами — разные

сорта частиц.

 

2)

Заметим, что корреляционные функции

отличаются в дан­

ном рассмотрении от корреляционных функций, применяемых для описания равновесных систем (см. примечание на стр. 73).

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