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

FTI_Lekcii_po_osnovam_fiziki_plazmy

.pdf
Скачиваний:
30
Добавлен:
29.05.2015
Размер:
1.25 Mб
Скачать

Формально радиус-вектор ведущего центра можно записать сле-

дующим образом:

 

 

 

 

 

 

 

 

 

 

 

R = r +

 

1

 

 

v h(r ) ,

(8.7)

 

 

 

 

 

 

 

 

 

 

ω

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

ρ = −

1

 

 

v h(r ) .

(8.8)

 

 

 

 

 

ω

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 8.1. Дрейфовое приближение: r - радиус – вектор частицы; R - ра-

диусвектор ведущего центра; ρ - радиус – вектор ларморовской спирали.

Рассмотрим, в общих чертах, метод получения дрейфовых уравнений. Продифференцируем выражение (8.7):

R = r +

1

 

 

+

1

 

 

 

 

 

 

 

 

ω

 

v h(r )

ω

 

v h(r )

 

 

H

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

ωH

 

v h(r ) .

 

 

 

 

 

(8.9)

ω2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

 

 

 

 

Подставляем в него уравнения движения (8.1), (8.2), а также учитываем, что

41

ωH =

ωH

+ (v)ωH ,

(8.10)

 

 

t

 

 

h =

h

+ (v )h.

(8.11)

t

 

Здесь h dh / dt - производная Лагранжа, а h / ∂t - производная Эйлера.

В полученном выражении величины h(r ) , ωH (r ) , E(r ) разла-

гаются в ряд Тейлора:

 

h(r ) = h(R) +(ρ )h(R),

(8.12)

и так далее. В итоге получается выражение для скорости ведущего центра

R = R(ρ,v, E(R), B(R),...),

(8.13)

где члены с ρ и v осциллируют с ωH , а члены, содержащие E(R), B(R) и т.д. изменяются медленно.

После усреднения получается следующее выражение для скорости движения ведущего центра:

R = v h

+

c

E B

+

v2

h α

B2

 

 

 

 

 

 

 

ωH

 

+

 

v2

h B .

 

 

(8.14)

 

 

 

 

 

 

 

 

 

 

 

 

2ωH B

 

 

 

 

 

Первый член этого уравнения описывает движение вдоль силовой линии, второй – электрический дрейф, третий – центробежный, а четвертый – градиентный дрейф.

42

На рисунке 8.2 представлена илюстрация вектора кривизны. Здесь α- вектор кривизны силовой линии (h / ∂s - производная вдоль

Рис. 8.2. Радиус кривизны силовой линии

 

 

силовой линии):

 

 

 

 

 

 

 

 

 

 

 

α = (h )h =

h .

 

 

(8.15)

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как

 

 

 

s ,

 

 

 

 

h =

 

 

 

 

 

(8.16)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

кр

 

 

 

α =

h

,

 

 

 

 

 

 

(8.17)

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

то

 

 

 

 

 

 

 

 

 

 

 

 

α

 

=

 

 

h

 

=

 

1

,

(8.18)

 

 

 

 

 

 

 

 

s

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

кр

 

 

где Rкр - радиус кривизны силовой линии.

Рассмотрим различные типы дрейфовых движений.

43

9. Электрический дрейф

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

На рисунке 9.1 магнитное поле направлено в сторону наблюдателя. Когда частица движется в направлении электрического поля, ее скорость увеличивается, и соответственно уменьшается, когда она движется в обратном направлении.

Рис. 9.1. Дрейф заряженной частицы в скрещенных электрическом и магнитном полях.

В результате, вдоль направления электрического поля частица не смещается, но в направлении, перпендикулярном ему возникает раз-

ность скоростей vD , так как составляющая скорости в одном направле-

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

vD

=

c

E B .

(9.1)

B2

 

 

 

 

 

 

 

 

 

 

44

 

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

10.Центробежный дрейф

При движении частицы со скоростью v вдоль искривлённой си-

ловой линии (рис. 10.1) с радиусом кривизны Rкр возникает дрейф,

обязанный своим происхождением центробежной силе инерции (так называемый центробежный дрейф)

F =

mv 2

.

(10.1)

 

 

ц

Rкр

 

 

 

 

 

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

Центробежную силу можно выразить через вектор кривизны силовой линии магнитного поля

Fц = mv 2α.

(10.2)

Она эквивалентна по своему действию электрическому полю

Eэкв = Fц / q.

(10.3)

45

Рис. 10.1. Центробежный дрейф: смещение заряженных частиц происходит в направлении, перпендикулярном плоскости рисунка

Отсюда, скорость центробежного дрейфа может быть записана следующим образом:

 

c

 

 

 

mv2c

 

 

(10.4)

v =

 

Eэкв

B

=

 

B α .

 

B2

qB2

 

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

11. Градиентный дрейф

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

Рис. 11.1. Градиентный дрейф

46

На рисунке 11.1. показано движение заряженных частиц в условиях градиента магнитного поля. Силовые линии магнитного поля выходят из плоскости рисунка в сторону наблюдателя. Напряженность магнитного поля больше в верхней части рисунка и меньше в нижней. Это приводит к вращению частиц вокруг силовых линий магнитного поля по сложной траектории. Ларморовский радиус больше в области слабого поля и меньше в области более сильного. В результате происходит дрейф заряженных частиц со скоростью:

v =

v2

h B .

(11.1)

 

 

 

2ωH B

 

 

 

 

 

 

 

Электроны и ионы дрейфуют в противоположные стороны.

12. Кинетическое уравнение с самосогласованным полем

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

 

f = f (r , v,t).

 

(12.1)

Функция f определяется

таким образом,

чтобы

величина

f (r ,v,t)d3rd3v имела

смысл числа частиц в элементарном объеме

d 3rd 3v шестимерного

фазового

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

Другими

словами,

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

Шестимерный вектор R6 ={x, y, z, vx , vy , vz } - это радиус -вектор частицы в фазовом пространстве, а его производная по времени имеет смысл шестимерной скорости V6 ={x, y, z, vx , vy , vz }.

Для того, чтобы найти полное число частиц сорта N , необходимо проинтегрировать функцию распределения по пространству координат и скоростей:

47

 

N = d 3r d3v f (r, v,t).

(12.2)

−∞

−∞

 

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

 

n = d 3v f (r , v,t),

(12.3)

−∞

 

 

j = d 3v f (r , v,t)v.

(12.4)

−∞

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

f

+ div (V

f ) = 0.

(12.5)

 

t

6

6

 

 

 

 

 

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

div (V f ) =

 

 

xf +

 

yf +

zf +

v

 

f

 

 

 

 

 

x

6 6

x

 

 

y

 

 

z

vx

 

 

 

 

 

 

 

 

+

 

 

vy f +

vz f .

 

 

 

(12.6)

 

 

 

 

 

 

 

 

 

 

vy

vz

 

 

 

 

Если в него подставить уравнения движения

 

 

 

 

 

 

 

 

r = v,

 

 

 

 

 

 

 

(12.7)

48

 

e

 

1

 

 

 

v =

 

E(r ,t) +

,

(12.8)

 

c

v, B(r ,t)

 

m

 

 

 

 

то можно показать, что

 

 

 

 

 

 

 

 

 

 

div6V6 = 0.

 

 

(12.9)

Отсюда следует, что

 

 

 

 

 

 

 

div (V

f ) = V

6

f

= v f + v

f .

 

(12.10)

 

6

6

6

 

 

v

 

 

 

 

 

 

 

 

 

 

 

Окончательно, вынося V6

из под оператора дифференцирования

div6 в выражении (12.5), запишем кинетическое уравнение в виде:

f

 

e

 

1

 

 

f

 

 

t

+ v f +

 

E(r ,t) +

 

v, B(r ,t)

v

= 0.

(12.11)

 

 

 

m

 

c

 

 

 

В этом уравнении электрическое и магнитное поля содержат компоненты, медленно меняющиеся в пространстве и во времени, они соот-

ветствуют

плавным изменениям скорости заряженных частиц

( E , B

), и микрополя (Eсл, Bсл ), соответствующие столкновениям

частиц и, соответственно, дискретным изменениям скоростей:

E =

E + Eсл,

(12.12)

B =

B + Bсл.

(12.13)

Микрополя отличны от нуля на расстояниях порядка средних расстояний между частицами. Поскольку работа с функцией распределения

предполагает, что в фазовом пространстве объемом d 3rd 3v содержится достаточно большое количество частиц, уравнение (12.11) следует усреднить по большому ансамблю частиц. В первом приближении это

49

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

B , но теперь уравнение (12.11) хотя формально и сохраняет свой прежний вид, но уже не учитывает столкновения частиц. Поля, входящие в него являются самосогласованными, так как определяются усредненной плотностью заряда и усредненным током в плазме. Эти поля описываются уравнениями Максвелла.

Кинетическое уравнение (12.11) с самосогласованным полем назы-

вается уравнением Власова. Если в нем f / ∂t оставить с левой стороны, а все остальные члены перенести в правую часть, то уравнение можно представить в виде:

f

= −r

f

v

f .

(12.14)

t

 

r

 

v

 

Отсюда следует его физический смысл: изменение функции распределения в данной точке пространства обусловлено приходом в нее частиц из других точек и с другими скоростями.

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

одной точки {r1,v1} в другую {r2 ,v2} , причем r1 = r2 .Такие пе-

рескоки описываются путем добавления в правую часть уравнения Власова дополнительного члена:

 

dfa

= fa +v f +

 

 

 

 

dt

 

 

 

 

t

 

 

 

 

qa

 

E(r

,t) + 1 v,

B(r,t) fa

= Stab .

(12.15)

 

m

 

 

c

v

b

 

a

 

 

 

 

 

Выражение Stab

называется интегралом столкновений частиц сор-

та a с частицами сора b . Суммирование проводится по всем сортам

50

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