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

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

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

§ 1. Вывод кинет, уравнений из цепочки уравн. Б БГК И 213

Используя разложение (1.22) и пренебрегая, как всегда, корреляционными функциями третьего порядка, можно первые два уравнения цепочки (1.20) и (1.21) пере­ писать в виде х)

dt дт1 V оt /столки

 

J ^ S t a ( r t, . . . , v 2;t)d r2dv2

(1.65)

и

 

 

 

 

 

д

Vz- i;)иг2 gi2'

 

 

 

 

\ dt ^ V i - 3ri

 

 

 

 

 

д

v*;

 

v2; O + ftsH -

 

г /9/, 9дфФ 1Яз .

 

+ -N

, dfi

дф2з

£13) dr3dv3.

(1.66)

 

аГ1 ё23

1 av2

аг5

 

 

Здесь по-прежнему не учитываются внешние силы. Кроме того, в этих уравнениях пренебрегается влиянием само­ согласованного поля на поток в пространстве скоростей, что дает возможность отбросить третий член в левой части уравнения (1.20) и четвертый и шестой члены в ле­ вой части уравнения (1.21). Очевидно, что для однородной системы такая процедура выполняется точно. Неоднород­ ность системы можно считать допустимой, если соответ­

ствующее самосогласованное

поле мало по сравнению

со средним корреляционным

полем.

Следует подчеркнуть, что малую пространственную зависимость ft в (1.65), (1.66) и последующих формулах следует понимать в указанном выше смысле.

Как и в случае уравнения Больцмана, мы стремимся исключить член gi2 из (1.65) с помощью уравнения (1.66). Теперь вспомним, что в безразмерном представлении (1.4)

первый член в правой части

(1.66)

пропорционален

Псв,

а второй — произведению

ПсвПпл.

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

вто­

рым членом в (1.66) можно пренебречь,

как членом более

высокого порядка малости [см. (1.64)].

 

 

1) Хотя используются частные функции

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

в силу основного предположения, данного на стр. 126, эффект от всех полевых частиц, представленный правой частью уравнения, записывается в виде произведения N на эффект только от одной Нолевой частицы с координатами (г2, v2).

214 Гл. 4. Неравное, состояния кулон, сист. с учетом корреляций

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

Учитывая (1.23), пренебрежем также членом gl2 в пра­ вой части (1.66). Тогда мы получим уравнение в частных производных относительно функции gl2:

Это уравнение можно переписать в виде

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

v1= v2 = 0, г12 (<') = r12 (t) + \i2(t'— t)

(1.69)

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

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

вправой части уравнения (1.68) через

А(гь г2, vt, v2;' t) =

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

 

t

°°)* (1.72)

g}2

t) — J -'1 ( • • • i t ) dt + gj2 (• • M

§ 1. Вывод кинет, уравнений из цепочки уравн.

Б Б Г К И 215

Коррелированное

 

состояние, соответствующее

време­

ни

t,

получается

 

из некоррелированного

состояния

#12 (■

• •>

—оо) =

0

при t = —оо.

на

длине

Предполагая,

что

изменением

функции f t

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

гс) можно

пренебречь,

получаем,

что подынтегральное выражение зависит только от отно­ сительного расстояния между частицами. Вводя новую

переменную

х =

t — t',

находим

 

 

 

gi2(|ri — г2|, vt, v2; t) =

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

= j

^ (ri(0 —ViT, r2(<) — v2T, v1? v2; t

t) dx.

(1.73)

Определяя поток J

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

 

 

 

 

 

(Ч г)

 

= —

 

 

(1.74)

 

 

 

 

\

o t /столки

 

 

 

 

из формул

(1.65),

(1.71)

и (1.73) находим

 

 

 

 

 

 

 

о

 

 

 

 

 

X R

r ( ^

I - ^

)

/l(ri’ Vi; 0 /i(rb V 2;

 

dx.

(1.75)

L

1

'

1

2 '

 

 

 

J ri2—V12T

 

Используя

разложение

 

 

 

 

 

l/i (Гь V t ;

0/i(*i,

V 2 ; / ) ] t _ T =

 

 

 

 

 

 

 

 

=

( /i/i) f - ( i- x ) [ - ^ - ( /1/i)]t + . . .

(1-76)

и подставляя его в выражение (1.75), получаем уравнение Ландау с запаздыванием.

Если функция распределения заметно не изменяется в течение времени взаимодействия {Т^> tc), можно огра­ ничиться первым членом разложения (1.76), т. е. пре­ небречь запаздыванием в системе. В результате мы полу­ чим

•>= - N Ш 2 J <Е><Е>.. {/• ('*• ™ <)a/i(x i; 0 -

- / , ( r . . v , ; l ) l),l(y

' ,>}<iv2. (1.77)

216 Гл. 4. Неравное, состояния кулон, сист. с учетом корреляций

Здесь коэффициент ( Е)(Е )Г2 определяется

выраже­

нием

 

<Е)(Е)Г2= j dr2 j Е [г12])(Е [г12 —v12t] dx,

(1.78)

о

 

аЕ является кулоновским полем.

Тензор ( Е)(Е )Г2 представляет собой автокорреля­

ционный тензор второго порядка, зависящий от vt и v2 и вычисленный в приближении прямолинейных траекто­ рий. Индекс г2 указывает, что усреднение проводится только по г2.

Уравнения (1.74) и (1.77) вместе представляют собой уравнение Ландау, хотя в них не приводятся еще точные значения автокорреляционных функций поля. Уравнение Ландау в виде, приведенном здесь, соответствует урав­ нению Больцмана для слабых взаимодействий [см. (1.8)]. Поэтому интересно отметить, что при его выводе мы вынуждены были сделать те же самые основные допу­ щения. Действительно, условие gl2 —оо) = 0 соответ­ ствует предположению Больцмана о начальном моле­ кулярном беспорядке. Условия (1.53) для пространствен­ ной и временной зависимости функции распределения отражаются здесь в требованиях L гс и tc.

Вычисление корреляционных функций поля

Это вычисление представляет собой геометрическую задачу. Для 6е решения удобно выбрать цилиндрическую систему координат с осью z, направленной вдоль вектора относительной скорости g = v2 — vt. Радиальную коор­ динату мы вновь обозначим через Ь (см. стр. 207).

В этой системе координат поле выражается следующим образом:

 

Ъcos ф'

 

Б(г» —г* + gr)

е

(1.79)

[b2+ (Z- gT)2]3/2 b sin ф

,2 —ГР

§ 1. Вывод кинет,

уравнений из

цепочки уравн.

Б Б Г К И

217

При этом корреляционный тензор записывается как

 

Е |0)(Е |, = е2 {[Ь2 + z2] [Ь2 + (z- gx)2]}- 3/2 х

 

 

Ъ2 cos2 ф

Ь2 sin ф cos ф

(z— gx) Ъcos ф\

 

Ь ^ тф совф

Ь2з т 2ф

(z — §х)Ь8т ф | ,

(1.80)

(bz cos ф

bz sin ф

z (z — gx)

/

 

а элемент объема как

 

 

 

 

dr2 = dzbdbdq.

(1.81)

Что касается пределов интегрирования по прицельному параметру Ъ, вспомним обсуждение области применимо­ сти метода Ландау — Фоккера — Планка, которое при­ вело к значениям Ьмин = rw и Ьмакс = X,D. (Точное дока­ зательство пределов применимости этого метода будет дано в п. 1.4 в связи с выводом уравнения Боголюбова — Ленарда — Балеску.) Интегрирование по углу ф в (1.80) тривиально и дает для недиагональных членов значение, равное нулю. Интегрирование по z дает для члена A zz тоже нулевое значение вследствие соотношения симметрии

j

Azz (

z, gx)dx =

J

Azz (z, gx) dx.

(1.82)

о

 

 

о .

 

 

 

Члены A xx

и A jn, интегрируются

элементарно,

так что

корреляционные

функции поля

принимают вид

 

 

 

 

 

/1

0

0\

 

 

<Е)(Е>Г2 = 2яе2- ^ Ц о

1

0 ,

(1.83)

или

 

 

 

\0

0

0 /

 

 

 

 

 

 

 

 

 

(Е)(Е>Г2= 2яе2 ( ^

А - )

(йг1—g ) ( g ) ,

(1.84)

где через I обозначается единичный тензор. Подставляя

этот результат и выражение (1.77) в (1.74), получаем окончательный вид уравнения Ландау [16]

(М±.\ \ б*

х |/ i (ть

О ЛТ I «2 \ 2 1

:2я г ( - )

v2; t) 9fi (ft, vi; d\i

А

9

f

g l ^^

-^ х

 

t)

7i (4,

Vj! t)

9fi (r«, v2; t)

} d \2.

 

 

 

 

d\2

(1.85)

 

 

 

 

 

218 Гл. 4. Неравное, состояния кулон, сист. с учетом корреляций

Уравнение Ландау является кинетическим уравне­ нием для случая слабого взаимодействия, не учитываю­ щего эффекты экранирования. Интересно из цепочки уравнений ББГКИ получить эквивалентное уравнение Фоккера — Планка. Это можно сделать, исходя из (1.77) и связывая корреляционные функции поля с коэффи­ циентами переноса в пространстве скоростей.

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

"В последующем изложении мы используем два первых момента изменения скорости в единицу времени: коэффи­ циент динамического трения

 

о

(1-8б>

 

 

и коэффициент диффузии в пространстве скоростей

 

Тс

Тс

(Av1)(Avi) = ±

\ dt'

J dt"E |4<)(Ег") . (1.87)

 

Ь

о

Здесь Е — полное поле, действующее на рассматривае­ мую частицу 1; скобки () означают усреднение по фазо­ вым координатам всех остальных частиц; тс — среднее время свободного пробега, которое пропорционально вели­

чине

TpA/lnA [см. (3.1.7)].

 

 

в

Удобно рассмотреть сначала коэффициенты диффузии

пространстве скоростей.

Простая замена переменной

t"

=

t' + т в (1.87) дает

 

 

 

 

Тс

т с - t '

 

 

(Av1)(Av1) = ^ - ( ^ - ) 2 j dt'

j <E|r )(E|r+ t>dt,

(1-88)

где подынтегральное выражение представляет собой авто" корреляционную функцию поля. Оценим (1.88), исполь" зуя два приближения, которые уже применялись ранее-

§ 1. Вывод кинет, уравнений из цепочки уравн. БВГКИ 219

а)

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

на

расстояниях,

больших

гс, не

дают вклада в интеграл,

так что существует максимальное

время

Д = r j v c,

в

течение

которого

могут взаимодей­

ствовать частицы,

 

причем

 

 

 

 

 

( Е |Г)(Е \t.+x ) =

0 для

т > Д.

(1.89)

б)

Изменения

Д

малы

в

смысле,

указанном

выше,

т. е. Д

(t — х) ж Д

 

(t),

когда использовалось «приближе­

ние в отсутствие запаздывания». Точнее говоря, если Т есть характерное время изменения Д , то условие Т^> tc

приводит

к

тому, что ( Е

| г)(Е

| t'+x)\ не зависит

от t'

для t' <

Т.

Отсюда

следует,

что

 

 

 

 

 

То

 

 

 

 

(Av1)(Av1) =

( ^ ) 2^

\d t'

j

(Е |0)(Е |t) dx =

 

 

 

с

о

 

1 Ц

< ( С

 

 

 

 

 

 

 

оо

 

 

 

 

= 2 ( ^ ) 2

j< E |0)(E|t)dr.

(1.90)

 

 

 

 

 

 

о

 

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

диффузии в

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

скоростей

(1.90) может

быть

записан как сумма N средних полей, возникающих при

парных взаимодействиях, т. е.

 

 

 

<Av1)(Av1>= 2N ( - ^ ) 2 j

(Е)(Е>Г2Д (г15 v 2; t)dvz.

(1.91)

Здесь мы использовали выражение

(1.78).

 

Рассмотрим теперь коэффициент динамического трения

 

 

Тс

 

 

 

 

 

( A v ^ - ^ - i -

[ <E[r12(f')]>d*'.

(1-92)

 

п г Т с

J0

 

 

 

 

в котором подынтегральное

выражение

можно записать

в виде

 

 

v

г

 

 

 

 

 

 

 

Е [г12 (*')] =

Е [г12 |о + g |о

-

2 - М

j

Е [г12 (Г")] dt”] .

 

L

 

0

0

 

(1.93)

 

 

 

 

 

 

220 Гл. 4. Неравное, состояния кулон, сист. с учетом корреляций

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

Е [г12 (£')] Е [г12 -f- gt']

(1.94)

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

 

E[r)2(0 ] = E [ri2 + g r ] - ^ | ( / 2 -£_ j (*'_ **)Е ||.Л \

(1.95)

о

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

<Лт*> ■= - (= •)1i

} л '

1

(»' ■-п < Ц | , ■:Е |-> ■**•

(1.9!))

Вводя новую переменную x = t' — t", находим

 

 

То

Г

 

 

- ( - S - ) 4

J

\

 

(1.97)

 

о

о

 

 

или, усредняя по значениям т — t',

<ЛЧ> = 2 ( £ ) 2-dvik - тv, ]

*' 1 <Е|.)(Е|»>*=

 

 

0

0

 

е \2

д

tc

 

 

(1.98) .

= 2 Ш

^ Г - 1 ( 1 - ^ ) < Е |« ) ( В |.) Л .

Так как, согласно

общему предположению, отношение

xc/tc л; Л/In Л 1, окончательно получаем

 

 

оо

 

 

<Ду,) = 2 ( ^ . ) г ^ - .

|< Е |,)(Е |,> Л =

 

 

о

 

 

2 N ( ~k)

 

j <Е)(Е)Гя/1 (rj, v„; t)dv9,

(1.99)

§ 1 . Вывод кинет, уравнений из цепочки уравн. В В Г К И 221

что приводит к важному соотношению между коэффи­ циентом динамического трения и коэффициентом диффу­ зии в пространстве скоростей:

<Avi>=±' ^ r ‘<Avi^ Avi^

(1.100)

Уравнение Фоккера Планка для слабого взаимодействия

С помощью соотношений (1.91) и (1.100) легко преобра­ зовать уравнение Ландау для слабого взаимодействия к уравнению Фоккера — Планка. С этой целью в урав­ нения (1.74) и (1.77) подставим выражение (1.91) и в ре­ зультате получим

( i r

L

[ т <Av.)<Av, >- ! £ -

 

- j v ( - i -)2 J

<E)(E)r,/i(r„ v,; t).

?>riv,] . (1.101)

Интегрирование по частям и использование выражения

(1.91) дает

 

6 t

 

= •/- •f4 (AviXAvj)' '

d \ i

{-гг)

 

 

\

 

/столки d \ i \ 2 >

 

 

 

 

 

( 1- 102)

Отсюда, учитывая (1.100), находим

( ж - ) ™

„ „ н - - ^ г '« Дт‘>« +

 

Это уравнение является уравнением ФоккераПланка

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

Приведенный выше вывод уравнения Фоккера — Планка не совпадает с выводом, данным самими автора­ ми. Следуя общей линии нашего изложения, мы получили

222 Гл. 4. Неравное, состояния кулон, сист. с учетом Корреляций

это уравнение из цепочки уравнений ББГКИ для част­ ных значений параметров Псв = О (е) и Ппл = О (е). Кроме того, мы использовали приближение парных столк­ новений, интегрирование по прямолинейным траекто­ риям и приближение в отсутствие запаздывания.

Для сравнения ниже приведем другой вывод урав­ нения, данный Чандрасекаром [19], который является более-коротким и не требует приближения парных столк­ новений, интегрирования по прямолинейным траекто­ риям и отсутствия запаздывания. Но при этом, как мы увидим из дальнейшего, уравнение не будет содержать всей той информации, которая была получена выше. Будем исходить из основного уравнения Чепмена — Колмого­ рова для марковских процессов

/ (v; £-f-A<)= j /(v — Av; t) PAt (v — Av, Av)dAv. (1.104)

Здесь функция / есть одночастичная функция распределе­ ния однородной системы в отсутствие внешних сил. Через PAt (v, Av) обозначена вероятность того, что частица, которая в момент времени t имела скорость v, за интервал времени At изменит свою скорость на величину Av. Рас­ кладывая подынтегральное выражение (1.104) в следую­ щий ряд:

/ (v — Av; t) PAt (v — Av, Av) =

 

 

= / (v; t) PAt (v, Av) -

Av

{fPAt) +

 

 

 

 

+ 4 ( Av- | t)2( / ^ ) -

+

---

,

(1-Ю5)

из (1.104)

получаем

 

 

 

 

 

 

/(v;

■+ A«)->(v; . ) = _ ^ , , [(Av)/1 +

 

 

 

 

 

 

 

+ T ^ - { - i r - 1<AvXAv>fl } -

+

•■•

(‘ -106)

При

этом

были использованы

следующие

соотношения:

 

 

P a, cZAv =

|

At х )

\ .

 

(1.107)

 

 

 

[ Ai(Av)(Av) J

 

 

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