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

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

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

Литература

193

К § 1

Власов А. А., ЖЭТФ, 8, 291 (1938).

К § 2

Bohm D., Gross Е. Р., Phys. Rev., 75, 1851, 1864 (1949).

Fowler Т. К., Phys. Fluids, 4, 1393 (1961).

4, 391 (1962).

Fowler T. К., Journ. Nucl. Energy, Part C,

Harrison E. B., Proc. Phys. Soc. (London), 79, 317 (1962); 80,

432

(1962).

 

 

33, 196

(1929).

Tonks L., Langmuir / . , Phys. Rev.,

К § 3 и 4

 

 

 

 

Алътшуль Л. М., Карпман В. И., ЖЭТФ,

49, 515 (1965).

Backus G., Journ. Math. Phys., 1, 178 (1960).

Bernstein / . В., Greene J. M., Kruskal M. D., Phys. Rev., 108,

546

(1957).

 

 

 

 

Denavit J., Phys. Fluids, 9, 134 (1966).

 

Derfler H., Simonen T. 6\, Phys. Fluids, 12, 269 (1969).

Drummond

W. E., Pines D., Nucl.

Fusion

Suppl., Pt. 2, 1049

(1962).

Gould R. W.,

Phys. Fluids, 4, 139 (1961).

Fried B. D.,

Gary S. P.,

Phys. Fluids,

10, 570 (1967).

 

Kallen G., Intuitive Anaiyticity, в книге Preludes in Theoreti­ cal Physics in Honor of V. F. Weisskopf (ed. A. Shalit), North-Hol- land Publ., Amsterdam.

McGune J. E., Phys. Fluids, 9, 2082 (1966).

Taylor E. C., Phys. Fluids, 8, 2250 (1965).

Turski A. J., Ann. Phys., 35, 240 (1965).

Веденов А. А., Велихов E. П., Сагдеев P. 3.,{УФН, 73, 701 (1961).

Weitzner H., Phys. Fluids, 6, 1123 (1963); 7, 476 (1964). Weitzner H., Comm. Pure Appl. Math., 18, 307 (1965).

К § 5

Case К. M., Ann. Phys. (New York), 7, 349 (1959).

McGune J. E., Phys. Fluids, 9,

1788 (1966).

Zelazny R. S., Ann. Phys., 19, 177 (1962).

Van Kampen

N . G., Theoretical

Methods in Plasma Physics,

W iley, New

York, 1967.

 

1 3 - 0 1 2 9 1

Глава 4

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

§1. ВЫВОД КИНЕТИЧЕСКИХ УРАВНЕНИЙ ИЗ ЦЕПОЧКИ УРАВНЕНИЙ ББГКИ

1.1.Основные положения

Начнем рассмотрение опять исходя из цепочки урав­ нений ББГКИ [см. (2.2.18)1 для частных функций рас­ пределения. Поскольку нас интересуют только внутрен­ ние состояния кулоновской системы, мы не будем учиты­ вать внешние силы. В докритической области, когда выполнено условие Л^>1, можно ограничиться случая­ ми s <^ N. Тогда цепочка уравнений представляет собой следующую систему:

дdt+ 2 1'дпд± .

г, k

 

 

 

г = 1

 

 

 

^ ^ S dr' S

r',

v';

t)dx'. (1.1)

 

i=l

 

 

 

 

Здесь

<j>(гг, rft)

— потенциальная энергия

взаимодей­

ствия,

зависящая

только от расстояния

| гг — rk |.

Вернемся к переменным Г-пространства (гг,

vf), ко­

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

поскольку

различие между

переменными ц-пространства

(*r, *v)

и Г-пространства (гг, v() было важно только в форма­ лизме Климонтовича.

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

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

Наряду с этим будут получены условия, при которых кулоновскую систему с корреляциями можно прибли­ женно описать функциями распределения, в которых

корреляции не учитываются, т. е.

парную

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

ную функцию можно

представить

в виде

(см. гл. 3)

/з (Г|, уи гj, v/, t)

= /* (гг, v;;

<) h {tj, \ }; t). (1.2)

Классификация по параметру связи и плотности

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

Пронормируем цепочку уравнений (1.1) относительно rc, vc и фс. В результате получим

* S о, 8

= ПсвПпл f

[ 2 4-

Ф(гь Г')- 4- h i d v '.

(1.3)

J

J i=1 dri

d\t

 

Здесь введены следующие безразмерные переменные:

Г

V

г,

1

fs = v lSVSfs

(1.4)

 

Vc.

 

 

 

Уравнение (1.3) в безразмерных переменных содержит только два характерных параметра

Псв= - ^ и Пnn = nrl

(1.5)

U IV с

Назовем Псв параметром связи, а Ппл — параметром плотности. Эти параметры удобны для характеристики метода разложения по малому параметру.

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

е1, то можно учесть все встречающиеся случаи различ-

13*

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

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

1)

Псв=

О (е),

1')

Ппл= 0 (г),

 

2)

Псв =

<9 (1),

2') Ппл= О (1),

(1.6)

3)

Псв= О (е-1),

3') Ппл= О (е-1).

 

Четыре комбинации Псв и Ппл, возникающие из групп (2), (3) и (2'), (3') соответственно, не удовлетворяют основ­ ным условиям разложения по малому параметру. Эти четыре комбинации соответствуют критическим и над­ критическим состояниям, в которых имеется сильная связь в одной или нескольких подсистемах. Такие состоя­ ния представляют интерес в случае жидкостей и твердых тел, но это выходит за рамки нашего рассмотрения, поскольку мы исследуем лишь системы, находящиеся в докритическом состоянии. Поэтому такие комбинации параметров Псв и Ппл ниже исключаются из рассмот­ рения.

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

3) Псв = О (г-1) и 1') Ппл = О (е),

(1.7)

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

Таким образом, остаются следующие четыре случая:

а) Псв= <9(1), Ппл= 0(е) разреженный газ,

б) Псв= 0 (е ), Ппл= 0(е) слабая связь,

(1.8)

в) Псв = 0(е), Ппл> 0 (1 ) дальнодействующие силы.

Эта терминология говорит сама за себя и соответствует той, которая была введена Фриманом [1].

Каждому случаю (1.8) соответствует свое кинетическое уравнение: разреженному газу — уравнение Больцмана;

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

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

уравн. Б Б Г К И

197

Применимость обычных кинетических

приближений

 

к описанию плазмы

 

 

Исследование, проведенное выше,

справедливо

для

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

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

Область О < г ^ О (rw). Здесь через rw обозначен классический радиус взаимодействия *)

Г ш = 4 -

(!.9)

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

rc= rw, фс = @, ус = ( - ^ ') /2>

(1.10)

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

Пе» = 1 , Ппл= ^ 2 .

(1.11)

Так как мы рассматриваем системы, находящиеся в докритическом состоянии, то

Л = 12яиЯв =

з е 3/а

(1. 12)

'У'Алп е3 > 1,

*) В последующем разложении используется терминология и соотношения, приведенные в приложении.

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

откуда получаем, что выражение для е в этой области имеет вид

6 4яЛ2 ^ ^ ’

(1.13)

 

а соответствующее кинетическое уравнение есть уравне­ ние Больцмана (Б).

Область О (rw) < г < О (г0). Чтобы установить пре­ делы изменения параметров Псв и Ппл в этой области, вычислим сначала соответствующие значения для г = =г О (го), где через г0 обозначено среднее расстояние меж­ ду частицами, равное

Характерными величинами в этом случае являются

Гс = Г0,

в

\1/2

(1.15)

 

 

которые приводят к следующим значениям параметров связи и плотности:

П с в = ( ^ ) 1/3, n M = J L = 0 (i).

(1.16)

Исходя из значений г0 и rw, указанных выше, получаем

9

( 1. 17)

1 > Ппл > 4ЯД2 •

Из этих соотношений следует, что при Л ^> 1 уравнения Ландау — Фоккера — Планка (ЛФП) дают адекватное опи­ сание лишь для центральной части области rw < г < г0, на границе же этой области они неприменимы.

Область О (г0) < г < оо. Для определения пределов изменения параметров Псв и Ппл отметим сначала, что при г - У оо

Псв-*0, ппл оо.

(1.18)

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

Используя предыдущие результаты, находим

 

( ^ ) 1/з> П св> 0 ,

А <

Ппл< о о .

(1.19)

Отсюда следует, что при Л

1 в этой области

применимо

уравнение Боголюбова — Ленарда — Балеску

(БЛБ).

На фиг. 19

сплошными линиями показаны области

применимости

различных

кинетических

уравнений.

0

 

r w

Г0

1

 

°о

1

------------ 1

I

 

I

_________Больцман

 

 

 

 

 

 

Ландау-Фокквр-Плат

 

 

 

 

Боголю бов-Ленард-Балеску

 

 

Фи г . 19.

Область применимости трех

кинетических

уравне­

 

 

ний (см. текст).

 

 

 

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

У р а в н е н и е Б. Основные предпосылки приближе­ ния Больцмана при многократных столкновениях в обла­ сти значений rw ^ г ^ XD не нарушаются, так как они приводят только к слабым отклонениям. Аналогично парным взаимодействиям многократные столкновения не­ зависимо можно учесть с помощью линейной теории. Поэтому больцмановское приближение справедливо

ВПЛОТЬ ДО Г =

У р а в н е н и е БЛБ. В методике Боголюбова — Ленарда — Балеску, которая справедлива до гмин = г0, проводится точный учет коллективных эффектов экрани­ рования. Существенным ограничением данной методики является тот факт, что в ней рассматриваются лишь взаимодействия, приводящие к слабым отклонениям, а это условие нарушается только для г < гw. Поэтому при­ ближение Боголюбова — Ленарда — Балеску дает пра­

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

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

^"мин = ^и>

У р а в н е н и е ЛФП. Ландау в своей работе [16] утверждал, что его рассмотрение без учета экранировки может и должно быть справедливым вплоть до г = О (XD). Ниже (см. п. 1.4), основываясь на уравнении Боголюбо­ ва — Ленарда — Балеску, мы покажем, что утверждение Ландау было правильным.

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

несколько

изменятся.

 

от областей 0 ^

г ^

гш

и A,D ^

^

Покажем,

что

вклад

г -< оо

 

пренебрежимо

мал

по

сравнению

с

вкладом

от

области

Гш <

г <

A,D.

 

 

 

 

 

 

Чтобы доказать это утверждение, перепишем первые

два уравнения

цепочки

(1.1)

 

 

 

 

 

д/1

|

 

 

 

 

 

d/i

 

 

 

 

 

 

 

dt

1

ЙГ)

 

m

^

dvi

 

 

 

 

 

 

 

 

( df l \

 

 

 

 

N

д

Гл

f<9<f>(ri,

rj) .

 

 

,

 

( ~dt )столк„

=

m

W , *

J

J - ^ F T

2- Ы

*>' Vl’

t] dY*

И

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( ж + ¥ ‘-

dri

y> -dF jSi2-

 

 

 

 

 

dt

1

 

 

 

 

 

 

 

 

 

 

 

1

дфiz

 

 

 

 

 

 

 

 

 

 

 

m

dr7*

 

 

 

t/i (Г1>Vb t)fi(r2, v2; t) + g12)-

 

N

dfi

 

 

 

I 4

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J *» 1

 

 

V*; ^dv3 -

 

 

 

 

£

Jd

^

l

A

^ d l3 y ^

gi3dy3_

 

 

 

 

“ £ ^7 ‘ J M ^

 

y *' *) d y * ~

 

 

 

- £ J 4

*.-<>•

( 1.21)

 

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

Здесь использовано представление /2 и / 3 через корре­ ляционные функции g:

/2 (^l! ^2’ ^i, v2; t) = /1 (rj, Vi, t) /1 (r2, v2, t) -\-gi%

и

(1.22)

/ 3 (ги • • •> v 3; t) = / 1/ 1/1 + / 1 ^ 2 3 + fiSi3 +

f i g n + gm -

В уравнении (1.20) (E) является самосогласованным полем в приближении Власова. Теперь используем тот факт, что во всех трех обычных методах кинетического приближения, рассматриваемых до сих пор, пренебрегалоСь тройной корреляционной функцией g123, поэтому в (1.21) можно опустить последний член. Далее, введем безразмерные переменные (1.4) в уравнение (1.21) и полу­ чим, что для всех комбинаций рассматриваемых парамет­ ров разложения выполняется соотношение

gi2 = О (Псв/j/t).

(1.23)

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

[ 4 j]

= 0 { ^ — fi {ri, y i\t) \ d r z X

 

\ bt /столки

\ m v с 11'

J

 

 

у

j ^ n c ^ r a , V 2 ; t)dva) .

(1.24)

Используя

условие нормировки

 

 

j /1 (r2, v2; t)d \2=r-^-

(1.25)

и переходя к относительным координатам г = | г4 — г 2 |, получаем

оо

(V±)

\ ^ n

CBr2d r) . (1.26)

\ 6t /столки

\ vcm

J0 dr

/

Разобьем область интегрирования в (1.26) на три части, соответствующие трем областям Б, ЛФП и БЛБ, опре­

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

деленным выше. Тогда можно записать

( ~1т ) = 0 [ ( С б + Сл ф п + Сб л б ) /1 ] . (1 -2 7 )

\ 0 г /столки

Оценим величину этих коэффициентов.

Область Б. г ^ О (rw). Основной вклад в интеграл (1.26) дает область вблизи rw, поэтому на основании (1.11) можно подставить в интеграл значение Псв = 1 . В ре­ зультате получим

rw

 

Cv = ° { ^ r c H r * d r ) = 0(nvcrl),

(1.28)

Сб= ° ( т ) ’

(1-29>

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

Область БЛБ. г ^ О (A,D). Вследствие того что дебаев­ ское экранирование уменьшается по экспоненциальному закону, главный вклад в интеграл дает область вблизи XD, поэтому можно считать, что Псв = А-1. Оценка вели­ чины (1.26) в данной области приводит к

Сб л б = О

l . r. * ) _ 0 ( £

i a ) l(i.30)

ИЛИ

 

 

 

С б л б = О ( -£■) .

(1.31)

Это означает, что в области БЛБ корреляционный член в уравнении Боголюбова — Ленарда — Балеску по поряд­ ку величины равен корреляционному члену в уравнении Больцмана для области Б.

Область ЛФП. О (rw) < г <. О (A,D). В этой области параметр связи имеет порядок

(1.32)

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