Морозов Введение в теорию горячей плазмы Ч.1 2011
.pdf≡ Im pe/2, т.е. сравним по порядку величины с электронной плазменной частотой.
Глава X. Столкновения в плазме
Как уже указывалось, в горячей плазме целый ряд процессов происходит достаточно быстро, и столкновениями частиц при их изучении можно пренебречь. С другой стороны, многие процессы, такие, например, как установление равновесия, перенос тепла и частиц, распространение медленных возмущений существенным образом зависят от столкновений. В главе V мы уже рассматривали некоторые свойства интеграла столкновений, стоящего в правой части кинетического уравнения. В этой главе мы изучим интеграл столкновений подробнее.
X.1. Интеграл столкновений в форме Ландау
Как уже указывалось, в левую часть кинетического уравнения (IV.2.6) входят электрические и магнитные поля, создаваемые как внешними источниками, так и собственными усреднёнными зарядами и токами в плазме. Но
211
помимо усреднённых полей в плазме существуют поля, которые существенно превышают усреднённые. Эти поля действуют на частицу, если она достаточно близко подходит к другой частице. Действие таких полей обычно сводят к мгновенным столкновениям. Интеграл столкновений в правой части кинетического уравнения учитывает действие таких полей. Оно сводится к мгновенному перебросу частицы из одного элемента импульсного пространства в другой. Так как столкновение происходит практически мгновенно, частица остаётся в том же элементе обычного пространства. Столкновения могут быть парными, тройными, и т.д. Однако в плазме обычно плотность частиц достаточно мала, и тройными, а также столкновениями более высокого порядка можно пренебречь. Итак, интеграл столкновений определяет баланс частиц, приходящих в результате столкновений в элемент фазового объёма drdp и уходящих из него.
Пусть в объём dr попадают две частицы с импульсами (это могут быть и обобщённые импульсы) p и p1. Столкновение "выбивает" частицу из своего импульсного объёма dp. При этом она может оказаться в произвольной точке импульсного пространства p′. Вторая частица после столкновения будет иметь импульс p′′, который связан с импульсами p, p1, и p′ законом сохранения импульса (считается, что время столкновения настолько мало, что усреднённые силы не могут изменить суммарный импульс сколько-нибудь заметно):
212
p p1 p′ p′1. (X.1.1)
Cуммарное число частиц, уходящих вследствие |
|
столкновений из импульсного объёма dp в единицу |
|
времени, определяется соотношением |
|
dN ∑ dpf p fi p1 w p, p1 |
→ p′, p1′ dp1dp′dp1′ . (X.1.2) |
i |
|
В силу закона сохранения |
импульса функция |
w p, p1 → p′, p′1 содержит −функцию от p′1. Здесь f p − число частиц в единице фазового объёма с
импульсом p, fi p1 − число частиц с импульсом p1, с которыми сталкивается первая частица. Эти частицы могут быть частицами разных сортов, например, электроны и ионы. Суммирование
ведётся |
по |
сортам |
частиц. |
Величина |
w p, p1 → p′, p1′ − |
это вероятность того, что две |
|||
частицы с начальными импульсами p и p1 после |
||||
столкновения будут иметь импульсы p′ и p1′ . |
||||
С другой стороны, в результате столкновений |
||||
частицы могут приходить в объём dp. Их число |
||||
определяется выражением: |
|
|
||
dN′ |
∑ dpf p′ fi p1′ w′ p′, p1′ → p′, p1 dp1dp′dp1′ . (X.1.3) |
|||
|
i |
|
|
|
Тогда столкновительный член, описывающий баланс между двумя этими процессами, принимает вид
213
St ∑ f p′ fi p1′ w′ p′, p1′ |
→ p, p1 dp1dp′dp1′ − (X.1.4) |
i |
|
−∑ f p fi p1 w p, p1 → p′, p′1 dp1dp′dp′1
i
Соотношение унитарности
w′ p′, p′1 → p, p1 dp1dp′ w p, p1 → p′, p′1 dp1dp′
позволяет преобразовать выражение (X.1.4) к виду
St ∑ f p′ fi p1′ w′ p′, p1′ |
→ p, p1 dp1dp′dp1′ . (X.1.5) |
i |
|
Кинетическое уравнение со столкновительным членом в виде (X.1.5) называется уравнением Больцмана.
Здесь vотн − относительная скорость частиц
случае столкновений нейтральных частиц сечение быстро спадает с ростом прицельного параметра. В случае же кулоновских столкновений сечение падает с ростом прицельного параметра, а значит, и с падением переданного импульса, медленно. Таким образом, в плазме основной вклад в интеграл столкновений дают столкновения с малой передачей импульса.
Выразим интеграл столкновений через переданный импульс q :
214
p′ |
p q. |
(X.1.6) |
Из закона сохранения (X.1.1) имеем |
|
|
p1′ |
p1−q. |
(X.1.7) |
Так как основную роль играют столкновения с малой передачей импульса, p q, p1 q, функции распределения в интеграле столкновений (X.1.5) можно разложить в ряд:
|
f p′ fi 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