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

Морозов Введение в теорию горячей плазмы Ч.1 2011

.pdf
Скачиваний:
28
Добавлен:
12.11.2022
Размер:
1.82 Mб
Скачать

Im pe/2, т.е. сравним по порядку величины с электронной плазменной частотой.

Глава X. Столкновения в плазме

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

X.1. Интеграл столкновений в форме Ландау

Как уже указывалось, в левую часть кинетического уравнения (IV.2.6) входят электрические и магнитные поля, создаваемые как внешними источниками, так и собственными усреднёнными зарядами и токами в плазме. Но

211

помимо усреднённых полей в плазме существуют поля, которые существенно превышают усреднённые. Эти поля действуют на частицу, если она достаточно близко подходит к другой частице. Действие таких полей обычно сводят к мгновенным столкновениям. Интеграл столкновений в правой части кинетического уравнения учитывает действие таких полей. Оно сводится к мгновенному перебросу частицы из одного элемента импульсного пространства в другой. Так как столкновение происходит практически мгновенно, частица остаётся в том же элементе обычного пространства. Столкновения могут быть парными, тройными, и т.д. Однако в плазме обычно плотность частиц достаточно мала, и тройными, а также столкновениями более высокого порядка можно пренебречь. Итак, интеграл столкновений определяет баланс частиц, приходящих в результате столкновений в элемент фазового объёма drdp и уходящих из него.

Пусть в объём dr попадают две частицы с импульсами (это могут быть и обобщённые импульсы) p и p1. Столкновение "выбивает" частицу из своего импульсного объёма dp. При этом она может оказаться в произвольной точке импульсного пространства p. Вторая частица после столкновения будет иметь импульс p′′, который связан с импульсами p, p1, и pзаконом сохранения импульса (считается, что время столкновения настолько мало, что усреднённые силы не могут изменить суммарный импульс сколько-нибудь заметно):

212

p p1 pp1. (X.1.1)

Cуммарное число частиц, уходящих вследствие

столкновений из импульсного объёма dp в единицу

времени, определяется соотношением

dN dpf p fi p1 w p, p1

p, p1dp1dpdp1. (X.1.2)

i

 

В силу закона сохранения

импульса функция

w p, p1 p, p1 содержит функцию от p1. Здесь f p число частиц в единице фазового объёма с

импульсом p, fi p1 число частиц с импульсом p1, с которыми сталкивается первая частица. Эти частицы могут быть частицами разных сортов, например, электроны и ионы. Суммирование

ведётся

по

сортам

частиц.

Величина

w p, p1 p, p1

это вероятность того, что две

частицы с начальными импульсами p и p1 после

столкновения будут иметь импульсы pи p1.

С другой стороны, в результате столкновений

частицы могут приходить в объём dp. Их число

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

 

 

dN

dpf pfi p1wp, p1p, p1 dp1dpdp1. (X.1.3)

 

i

 

 

 

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

213

St f pfi p1wp, p1

p, p1 dp1dpdp1(X.1.4)

i

 

f p fi p1 w p, p1 p, p1 dp1dpdp1

i

Соотношение унитарности

wp, p1 p, p1 dp1dpw p, p1 p, p1 dp1dp

позволяет преобразовать выражение (X.1.4) к виду

St f pfi p1wp, p1

p, p1 dp1dpdp1. (X.1.5)

i

 

Кинетическое уравнение со столкновительным членом в виде (X.1.5) называется уравнением Больцмана.

Здесь vотн относительная скорость частиц

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

Выразим интеграл столкновений через переданный импульс q :

214

p

p q.

(X.1.6)

Из закона сохранения (X.1.1) имеем

 

p1

p1q.

(X.1.7)

Так как основную роль играют столкновения с малой передачей импульса, p q, p1 q, функции распределения в интеграле столкновений (X.1.5) можно разложить в ряд:

 

f pfi p1

f p fi p1

(X.1.8)

q

f p

 

fi p1 f p

fi p

.

p

p1

 

 

 

 

Процессы с малой передачей импульса можно приближённо описать диффузией в импульсном пространстве,

St ≈ −divpS

S

.

(X.1.9)

 

 

p

 

Вычислим плотность потока частиц в импульсном пространстве S . Рассмотрим единичную площадку, расположенную в некоторой точке p импульсного пространства, перпендикулярную вектору p. Пусть q _ проекция вектора q на направление вектора p. Число частиц, проходящих через эту площадку в единицу времени в направлении p, т.е. частиц, для которых q 0, можно выразить так:

215

S

d3q d3p1w p, p1, q f p fi p1 q .

(X.1.10)

 

i

q 0

 

Выбранную нами площадку пересекают только частицы, импульсы которых лежат в малом интервале от p q до p, и интегрирование по dp можно заменить просто умножением на длину интервала q .

В противоположном направлении поток частиц равен:

S d3q d3p1w p, p1, q f p q fi p1 q q . (X.1.11)

i

q 0

Используя условие

 

w p, p1, q w p, p1, q

(X.1.12)

и разлагая в ряд Тейлора функции распределения, вычитаем из выражения (X.1.10) выражение (X.1.11):

S

S S

d3qd3p1w p, p1, q

(X.1.13)

 

 

 

 

 

i q 0

 

 

 

 

 

 

 

 

 

 

f p

fi p

1

 

fi p

1

 

f p

 

 

q q .

 

 

 

 

 

 

 

 

 

p1

 

p

 

 

 

Вследствие соотношения (X.1.12) интеграл по q (X.1.13) можно представить как половину интеграла от до . Вероятность w можно выразить через сечение рассеяния и относительную скорость

216

|v v1 |

 

 

 

 

 

 

 

 

 

 

 

 

 

 

wd3q |v v1 |d .

 

 

 

 

(X.1.14)

Тогда величину S можно представить так:

 

 

f p

fi p1

fi p1

f p

 

 

 

 

S d3p1

 

 

 

B , (X.1.15)

 

 

 

p

 

i q 0

 

 

 

 

 

 

p1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

1

q q |v v1 |d .

 

 

 

 

(X.1.16)

 

2

 

 

 

 

Вычислим тензор B для кулоновских столкновений. При малых углах рассеяния переданный импульс q перпендикулярен относительной скорости, т.е.

q q v v1 0.

(X.1.17)

Следовательно, тензор B является поперечным по отношению к вектору относительной скорости

v v1.

Из соотношений (X.1.15) и (X.1.17) видно, что для максвелловских функций распределения S обращается в ноль, так как производная по скорости от максвелловской функции пропорциональна скорости.

Сечение может зависеть только от относительной скорости. Единственный возможный поперечный тензор должен иметь вид:

217

B

 

1

B

v v1 v v1

,

(X.1.18)

2

v v1 2

 

 

 

 

 

 

где скаляр B имеет вид

B

1

q2|v v1 |d .

(X.1.19)

2

Рассмотрим две частицы в системе их центра инерции. В этой системе при малых углах рассеяния

q |v v1 | .

Здесь

m1m2

приведенная масса.

m1 m2

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

 

B

1

2|v v1 |3 2d .

 

 

 

 

(X.1.21)

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

Величину

t

 

1 cos d

 

1

2d

 

2

называют транспортным сечением.

 

 

 

 

 

 

 

 

Воспользуемся

формулой

Резерфорда

для

сечения для малых углов рассеяния

 

 

 

 

 

 

 

 

d

4 eei 2

 

dO

8 eei 2

 

 

d

.

(X.1.22)

2|v v1

|4 4

 

2|v

v1 |4

3

 

 

 

 

 

Здесь dO 2 sin d элемент телесного угла, e и ei заряды частиц. Транспортное сечение, таким образом, выражается так:

218

t

4 eei 2

L,

(X.1.23)

2|v v1

|4

 

 

 

Величина L d / называется кулоновским

o

логарифмом. Легко видеть, что кулоновский логарифм расходится как на нижнем, так и на верхнем пределе. Расходимость на верхнем пределе связана с тем, что мы проводили вычисления при малых углах рассеяния. Однако слабая логарифмическая расходимость позволяет ограничить углы рассеяния величиной max 1 c хорошей степенью точности. Расходимость на нижнем пределе связана с тем, что кулоновское сечение медленно падает с ростом прицельного расстояния, т.е. с уменьшением угла рассеяния. Однако в плазме поле рассеивающей частицы экранируется другими частицами на расстоянии порядка дебаевского радиуса rD. Поэтому интеграл можно обрезать на нижнем пределе величиной min.

Угол рассеяния при малых можно оценить так:

 

 

q

 

; q

F

 

|eei |

 

a

 

 

 

 

 

 

p

 

 

 

 

 

 

.

(X.1.24)

 

 

a2

|v v1 |

 

 

 

 

 

Здесь F кулоновская сила, a минимальное

расстояние между частицами. Для минимального угла рассеяния можно положить a rD. Таким образом,

219

min

|eei |

1

 

; L ln

rD |v v1

|2

.

(X.1.25)

rD

 

|v v1 |2

 

|eei |

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заметим,

 

что

 

при

электрон-ионном

 

столкновении me, |v v1 | vTe.

 

 

 

 

 

 

 

Окончательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S 2 eei 2L

 

 

 

 

 

 

(X.1.26)

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fi p1

fi p1

f p

 

 

 

 

 

 

d3p1

f p

 

 

 

 

 

 

 

p

 

 

 

 

q 0

 

 

 

 

p1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|v v1 |2 v v1 v v1

 

 

 

 

 

 

 

 

 

 

 

|v v1 |3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интеграл столкновений в форме

 

 

 

 

 

 

 

 

 

 

 

 

 

St

S

,

 

 

 

 

 

 

 

(X.1.27)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

где S имеет вид (X.1.26), называется интегралом столкновений в форме Ландау.

220

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]