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

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

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

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

Отсюда находим порядок величины корреляционного чле­ на в уравнении Ландау — Фоккера — Планка в области ЛФП:

=

1

• f1-33»

или

с.ЛФП: ® ( ^ ) '

 

 

<‘-м>

Оценка этой величины отличается от оценок соответ­ ствующих вкладов в областях Б и БЛБ на множитель

1пА.

Из соотношений (1.29), (1.31) и (1.34) очевидно, что если выполнено условие

1 пА > 1,

(1.35)

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

сти rw

г А,с. Ниже

мы будем

предполагать, что

условие

(1.35), которое

значительно

строже условия

А^> 1, выполняется. Таким образом, сделанное выше утверждение доказано.

Выводы

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

1п А > 1

при условии, что области, в которых данное приближение неприменимо, будут исключены соответствующим выбо­ ром границ таких областей, а именно:

область

Б

rMaKo = ^D)

область

ЛФП

Пиане

область БЛБ

^*мин==Па.

Важно отметить, что вышеуказанные методы можно скомбинировать таким образом, что для устранения рас-

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

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

= + IБЛБ ^ЛФП-

(1.36)

СТОЛКИ

Интерпретация этого метода очень простая.

В области г < гш члены столкновений Ландау — Фоккера — Планка (/дфп) и Боголюбова — Ленарда — Ба-

леску (/блб) взаимно

сокращаются, точное же описание

системы дает

член столкновений Больцмана (/б).

В области

rw ^

г ^

все три члена совпадают,

поэтому рассматривается вклад только от одного из них. Наконец, в области г > взаимно сокращаются

члены столкновений Больцмана и Ландау — Фоккера — Планка; остающийся член столкновений Боголюбова — Ленарда — Балеску точно описывает столкновения в этой области и не имеет расходимости.

Уравнения такого типа, как уравнение (1.36), были получены различными авторами [2—5].

1.2. Уравнение Больцмана

Вывод уравнения Больцмана из цепочки уравнений ББГКИ

Уравнение Больцмана для разреженной системы частиц с короткодействующими силами может быть полу­ чено различными путями. Для последовательности изло­ жения выведем это уравнение из цепочки уравнений ББГКИ.

Учитывая, что

ПСв — О (1), Ппл — О (е) 1

(1.37)

и пренебрегая снова внешними силами, запишем два нервых уравнения цепочки в виде

N_ _д_

т д \-• j ^ г / г ( г , гь vb v2; t)d rd \z (1.38)

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

И

Г д .

д

, .

. д

1 З ф 12

( д

д \ - |

LW + Vl'

+

^Vz

Vl- * ' д 7

~ ~ ^ ~ д Г

' U v i

e v J J X

X h (t,r u yu y2; t ) = ^ ^ 2 ( w

J ” )-£ -J > dT*dy>-

(1-39^

V = 1

 

 

Здесь введена относительная

координата г =

r2 — rt.

Заметим, что в левой части уравнения (1.39) содержатся

члены, совпадающие

с

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

в правой части (1.38).

члены, проинтегрируем (1.39)

Чтобы исключить

эти

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

+j /z(r’ Г1’ •••)rfv2dr +

+ j

(v2 —v4).

df2 (r, rlt ...) dr dv2 =

 

 

 

dr

 

 

~ m d \ i’ J ^ р /г ( г , 4 , .. .)d \zdr.

(1.40)

Введем теперь

область

взаимодействия первой частицы

с s-й частицей

 

 

 

 

D s = {rs ||rs-ri]<a}

(1.41)

и область, дополнительную к данной, в конфигурацион­ ном Ys-пространстве >й частицы

Ds = {rs ||rs-n|>a} = Vs Ds.

(1-42)

Через о здесь обозначена область взаимодействия а « гс.

Далее, определим укороченную функцию распределения

порядка s:

S+1 dv.s+ l

j

drNdv

N’

(1.43)

 

fNC

 

dn

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

Из (1.42) и (1.43) получаем

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

 

 

 

/£ = /а-

(1-45)

 

Укороченная функция распределения /® дает плотность

вероятности найти частицы 1, 2,

. . ., s в положениях

г15

. . .,

\ s при условии, что ни одна из частиц s +

1, . . .

. .

., N

не находится в области взаимодействия

первой

частицы.

Отметим, что fi? (ru vt; t)

представляет

собой

плотность вероятности найти частицу 1 с координатами rj и vj, которая не взаимодействует ни с одной из других частиц. Заметим также, что, поскольку мы выделили частицу 1 среди остальных частиц, симметрия функции

относительно координат нарушается.

Вводя укороченные функции распределения в урав­ нения (1.38) и (1.40), получаем

(1.46)

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

+

/? = 1V j dv2j> fl (г, rt, • •.) (v2—vO dS,

 

s

(1.47)

где S является поверхностью сферы D%. Это уравнение представляет собой второе уравнение Грэда, получен­

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

ное из цепочки уравнений ББГКИ. Как показал Трэд [6], это уравнение является точным и может быть полу­ чено даже, если не пренебрегать функцией / 3.

В приближении Больцмана, когда Ппл<^1, можно считать, что функция совпадает с /3. В результате полу­ чаем уравнение

+ = N \ • • о (v2 —Vi)dS. (1.48)

Это еще не уравнение Больцмана, но оно существенно отличается от исходного общего соотношения, прибли­ жаясь к уравнению боль-

цмановского типа.

 

В общем случае пред­

 

полагается, что корреляци­

 

онный

 

член

учитывает

 

плавное

 

изменение числа

 

частиц, находящихся в за­

 

данном состоянии в ре­

 

зультате

непрерывного

 

влияния

сил

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

 

вия. Однако в (1.48) изме­

 

нение

числа

свободных

 

частиц (не взаимодейст­

 

вующих

 

с другими) опре­

 

деляется

в виде

разности

 

числа пар частиц, входя­

 

щих в сферу взаимодей­

 

ствия (уменьшение свобод­

 

ных частиц),

и числа пар

Ф и г . 20. Сфера столкновений

частиц,

покидающих сферу

для заданных прицельных пара­

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

(увеличе­

метров.

 

ние свободных частиц). Это обусловлено применяемым здесь методом укороченных

функций распределения и разложением по параметру Ппл. Чтобы привести уравнение (1.48) к виду уравнения Больцмана, введем полярные координаты (Ъ, ср) в диамет­ ральной плоскости (фиг. 20) сферы S, перпендикулярной вектору относительной скорости g = v2 — vt. При таком описании сфера проецируется на диск С. Так как сфера

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

проецируется на диск С с двух сторон, необходимо разли­ чать две полусферы:

5+ = {г € 5 ^ о}={^},

 

 

5- = { r6 5 |g.ds<0} = {r-}.

1

Полусферы S ± относятся соответственно к удаляющимся

и приближающимся частицам (см. фиг. 20).

в

виде

Поэтому уравнение

(1.48) можно переписать

о

 

 

Ж + У1'Й г = ЛГ 5 dvz 1 bdb J Slfz(r+, ■■•) —

оо

 

 

— f2(T~,...)]d(p. (1.50)

Функция /

2 (г±, r4, Vi, v2; t)

пока еще не

известна.

До сих пор

мы пользовались

условием

 

 

Ппл< 1 .

(1-51)

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

rt,

v(,

v;; t') = / 2 (r+, rlf vt, v2; t).

(1.52)

Здесь штрихом

обозначены величины до столкновения,

а через

г+, vl5

v2,

t — соответствующие величины

после

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

Введем еще два ограничения. Первое из них состоит

в том, что

 

cr< L, или a< ^vcT,

(1.53)

т, е. характерная длина L и характерное время Т изме­ нения функции /2 велики по сравнению с длиной взаимо­ действия и временем взаимодействия. Если это условие выполнено, мы можем записать

/а (г-'> *i, vj, v;; t') = / 2 (rl5 г4, v(, v^; t),

(1.54)

h (r+, rl5 Vi, v2; t) = / 2 (rt, rlt Yj, v2; t).

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

Второе ограничение заключается в том, что на поверх­ ности S~ до начала столкновения частицы находятся в молекулярном беспорядке, т. е. выполняются равенства

A (rt,

rlt

vt,

v2;

t)

= A (rj,

ViJ

t) A (**i,

v2; t),

(1 55)

A (ч,

ч,

A,

v;;

t)

= a (**!,

v;-,

о A ( A .

v*'> 0-

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

J [Мч, ч; 9Мч, A; *)—

—А(ч> Vi; 0А(ч> V2; t)] gb db d<p dv2. (1.56)

Другой способ получения уравнения Больцмана из уравнения Лиувилля

Уравнение Больцмана, приведенное выше, было полу­ чено методом Трэда [6]. Можно вывести уравнение Больц­ мана из уравнения Лиувилля, причем в основном здесь используются различные приближения. По-видимому, впервые серьезные математические попытки в этом на­ правлении предпринял Кирквуд [7]. Основным в прибли­ жении Кирквуда является предположение о сглаживании всех процессов во времени, которое по существу приводит к тому же самому эффекту, что и введение Трэдом укоро­ ченных функций распределения. Кроме того, в резуль­ тате такого предположения получается уравнение, в кото­ ром учитываются только ближние столкновения. Конеч­ но, в качестве основного допущения в приближении Кирк­ вуда используется также и предположение о молекуляр­ ном беспорядке.

В методе, предложенном Боголюбовым [8], уравнение Больцмана получается как первый член общего разложе­ ния. Основная идея метода Боголюбова заключается в утверждении, что функции распределения более высо­ кого порядка А зависят от времени только через функцию распределения первого порядка. Чтобы это утверждение было справедливым, функция распределения s частиц, которые находятся в достаточно перемешанном состоя­ нии, должна удовлетворять следующему соотношению:

1 4 - 0 1 2 9 1

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

fa (Гц V i , . . rs, V s) = fi (rlt V i ) . . . fi ( r s , V s ) И Любая

функция распределения более высокого порядка при обра­ щении времени должна автоматически достигать такого состояния перемешивания.

Грин [9] использовал приближение, аналогичное мето­ ду Боголюбова. Имеются также другие выводы уравнения Больцмана, в которых вместо уравнения Лиувилля исход­ ным уравнением является так называемое основное урав­ нение. В классической области уравнение Больцмана, например, было получено Проутом [10] и Проутом и Пригожиным [11]. Поскольку обсуждение общих свойств урав­ нения Больцмана не является целью настоящего исследо­ вания, читателя, интересующегося этим вопросом, мы от­ сылаем к статье Трэда [6], в которой дано его полное рас­ смотрение.

Приближение, используемое Больцманом [12]

В противоположность всем приведенным выше выводам больцмановского уравнения Больцман в своем оригиналь­ ном исследовании рассматривал однокомпонентную систе­ му. Для того чтобы иметь возможность в последующем использовать идеи гл. 2, запишем плотность распределе­ ния в виде, данном Климонтовичем [см. (2.1.1)]:

F (г, v; 0 = S 6 (r —M 0 )s (v — Vj (*)).

(1.57)

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

(1.58)

где через К обозначаются все внутренние силы взаимо­ действия. Однако истинная плотность распределения F не рассматривалась Больцманом. Более того, он исполь­ зовал среднюю плотность числа частиц, получающуюся в результате усреднения по времени и р-пространству, т. е

(1.59)

После соответствующего разбиения всего р-пространства на элементы Ар, а времени на интервалы АТ в результате

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

усреднения (1.59) получаем / б в виде ступенчатой функ­ ции. Очевидно, чтобы это определение средней плотности

имело смысл, должно выполняться

следующее условие:

Л/б < / б,

(1-60)

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

 

 

/БЛ ц » 1.

 

 

(1.61)

Учитывая указанные требования,

 

уравнение (1.58) можно

записать в виде

 

 

 

 

 

 

dt

i

у — — f

)

столкн

(1.62)

'

дт

\ dt

)

 

Здесь член, стоящий в правой части, определяется корре­ ляциями.

В действительности основная заслуга Больцмана заключается в получении приближенной функциональной зависимости корреляционного члена от /в в разреженном газе на основании концепции ближних столкновений.

 

Болъцмановская

концепция ближних

столкновений

включает несколько

постулатов:

 

и

1. Взаимодействия являются короткими, сильными

редкими

[Псв = 0 (1 ), Ппл< 1 ] . Это

означает, что:

а)

можно

пренебречь корреляциями более высокого

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

2. Если Т и L — характерные время и длина измене­ ния /в в конвективном ц-пространстве, то выполняются

неравенства

 

Г » т с и £ » Х с,

(1.63)

где тс есть среднее время между столкновениями, а Кс — средняя длина свободного пробега.

14*

212 Гл. 4.

Неравное, состояния Кулон, сист.

с учетом Корреляций

3.

Если на основании первого постулата ввести коор­

динаты взаимодействующих пар частиц до и после столк­

новения,

то координаты

пар частиц

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

являются стохастически независимыми (молекулярный

беспорядок).

постулатов

и фиг. 20, легко

Исходя из указанных

видеть, что столкновительный член (bfJbt)cwaKa можно представить, согласно правой части (1.56), как разность членов, один из которых учитывает возникновение состоя­ ния с г4 и vb а другой — разрушение этого состояния.

Особенного внимания заслуживает тот факт, что в принципе результат Больцмана сильно отличается от уравнения (1.56), полученного из цепочки ББГКИ. Уже по своему смыслу функции /j отличаются в том и другом случае. Такое отличие находит отображение также в усло­ виях (1.60) и (1.61), которые, если они нарушаются, могут сделать уравнение Больцмана недействительным или, более точно, сделать его неприменимым к одночастичной системе. Эти интересные проблемы и следствия из них с точки зрения экспериментальной проверки изучались различными авторами [13, 14].

1.3. Уравнение Л андау — Фоккера — Планка

Уравнение Ландау со слабой связью

Согласно приведенной классификации (1.8), уравнение Ландау характеризуется следующей величиной парамет­ ров связи и плотности:

Пев = 0(e), Ппл = О (е).

(1.64)

Сравнивая эти соотношения с основными условиями (1.37) для уравнения Больцмана, мы напомним, что уравнение Ландау можно получить из разложения уравнения Больц­ мана для предельного случая слабого взаимодействия. Такое разложение хорошо известно (см., например, [15]). Однако, не нарушая общую линию настоящего изложе­ ния, выведем уравнение Ландау непосредственно из цепоч­ ки уравнений ББГКИ. Это также обусловлено тем, что в последующем выводе уравнения Ленарда — Балеску мы будем исходить из такого приближения.

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