Материал: LS-Sb89541

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Величина P = v + e A – обобщенный импульс. Уравнение для потока частиц m

имеет вид

v × rot P = 0.

Для «нормальных», или регулярных, потоков: rot P = 0 . Такие потоки могут создаваться эквивалентными эмиттерами в отсутствие потока магнитного поля через эмитирующую поверхность. Для стационарного или лами-

нарного потока, когда P = 0 , уравнение непрерывности примет следующий

t

вид: div j = div(r × v) = Ñ(r × v) = 0 .

4.4. Моделирование интенсивных электронных пучков

Гидродинамическая модель потока описывает движение пучков, однако имеет многие допущения. Такой подход нельзя использовать для описания интенсивных потоков заряженных частиц. Мера интенсивности электронных потоков называется «первеанс»:

P =

I

,

U 3/2

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

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

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

26

тие самосогласования полей построено на взаимной зависимости напряженности поля и плотности пространственного заряда.

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

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

4.5. Движение частиц в плазме

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

Основные параметры плазмы: 1) заряд электронов и ионов; 2) масса электронов (me) и масса ионов (mi); 3) концентрация электронов (ne) и ионов

(ni); 4) средняя кинетическая энергия (или температура) электронов (Te) и

ионов (Ti). Заряд и масса – фундаментальные физические константы, а тем-

пература

электронов и ионов

 

определяется формулами kT

=

1

m

v2

и

 

 

 

 

 

 

 

 

 

 

e

3

e

e

 

 

1

 

 

 

 

 

 

 

 

 

 

kT =

m

v2

. Как правило, Т

 

> T

. Существует несколько подходов для

 

е

i

3

i

i

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

описания движения частиц в плазме.

Описание плазмы с помощью кинетических уравнений. В кинетиче-

ской теории каждый сорт частиц характеризуется функцией распределения f = f (r, v, t ), зависящей от координат, скоростей и времени, т. е. функция распределения имеет смысл плотности частиц в пространстве.

Полное число частиц определенного сорта определяется как

+∞

+∞

 

 

N = d

3r f (r, t, v)d

3r ,

 

−∞

−∞

 

 

 

+∞

 

 

где d 3r = dvxdvydvz , и плотность частиц: n = f (r, t, v)d

3r.

−∞

27

Таким образом, N = n(r, t)d 3r. Функция распределения f – это плотность частиц в шестимерном пространстве (шестимерное, так как r зависит от координат и v тоже зависит от координат).

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

 

 

 

 

 

 

 

 

 

 

 

 

f

+ div

6

(V

f ) = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где введено обозначение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div

6

(V f ) =

x × f +

y × f

+

z × f

+

 

v

x

× f

+

 

v

y

× f +

v

z

× f .

 

 

 

 

 

 

 

 

 

6

 

 

x

y

 

z

 

 

 

 

 

 

 

vx

 

 

 

vy

 

vz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заряженные частицы при движении создают электрические и магнитные

поля,

поэтому

средняя

скорость

потока

частиц

 

определяется

как v =

=

e

E(r, t) +

1

[v, B(r, t)]

. Тогда div

6

(V

 

f )

= V Ñ

f = vÑf + v f .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

6

 

 

6

6

 

 

 

 

 

v

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

частиц можно записать в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vÑf +

f +

e

E(r, t) +

1

[v, B(r, t)] f = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

– уравнение Власова. Уравнения Больцмана описывают короткодействующее взаимодействие, уравнение Власова – дальнодействующие взаимодействия.

f

Смысл уравнения: если t оставить в левой части, а все остальное пере-

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

Общий вид f

= -r f

- v f .

 

t

r

v

 

Уравнение Власова предполагает условие, что поток частиц несжимаем.

Левая часть уравнения Власова – полная производная

f = 0. Это и означает,

 

 

 

t

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

28

В плазме все частицы коллективно движутся и взаимодействуют между собой. Описать траектории движения этих частиц невозможно, поэтому для решения плазменных задач используют статистические методы. Основная статистическая характеристика для описания движения частиц в плазме – функция распределения. Функцию распределения f = f (r, v, t ) можно пред-

ставить в виде произведения:

f (r , v, t ) = n(r, t) × fv (r, v, t ) ,

где fv (r, v, t ) – функция распределения частиц по скоростям, определяющая относительное число частиц определенной скорости в интервале vdv.

Такие перескоки описывают, добавляя в правую часть уравнения Власо-

ва дополнительную величину

 

dfa

= Stab , где суммирование идёт по всем

 

dt

 

 

b

 

 

 

сортам частиц, а смысл Stab

число частиц, которое появляется или исчезает

в единицу времени в единице шестимерного объема фазового пространства. Интеграл столкновений для лоренцевой плазмы (плазма, которая состо-

ит из электронов и бесконечно тяжелых холодных ионов с зарядом z >> 1) имеет вид

Stei = 2pe4 z2niL Dθ, ϕ × f , me2V 3

где ni – число электронов, рассеивающихся в единице объема; Λ – кулонов-

ский логарифм; V3 – элемент фазового пространства; Dθ, ϕ × f – угловая часть

лапласиана в сферической системе координат. В лоренцевой плазме энтропия возрастает. Это справедливо при столкновении одного сорта частиц.

4.6. Гидродинамическая модель плазмы. Плазма как жидкость

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

Описание движения частиц в плазме. Уравнение движения отдельной заряженной частицы имеет вид (4.2). Используют допущения, что столкновения и тепловое движение в плазме отсутствует. В таком случае предполагают, что все частицы элемента жидкости движутся вместе и их средняя скорость (u) совпадает с отдельной скоростью частицы.

29

Нужно получить уравнение для элементов объема жидкости, зафиксированных в пространстве. В таком случае рассматривают в одномерном пространстве любой параметр жидкости G(x, t). Это может быть плотность, вязкость, скорость течения, сжимаемость, температура и т. п. Запишем уравнение для параметра жидкости G(x, t) в системе отсчета, движущейся с жидко-

стью

dG(x, t)

=

G +

G ×

dx

=

G + u ´

G .

 

 

 

dt

t

x dx

 

t

x

Для трехмерного случая получим dG = G + (u ×Ñ)G . dt t

Первое слагаемое в правой части уравнения – изменение величины G в фиксированной точке пространства; второе слагаемое – изменение величины G при перемещении наблюдателя вместе с жидкостью в область пространства с другими значениями G. Если рассматривать уравнение для скорости жидкости (где G(x, t) – это u – скорость жидкости), то

u

 

 

 

 

n × m

t

+ (u

×Ñ)u

= nq (E + [u

´ B]).

 

 

 

 

 

При учете тепловых движений в правую часть уравнения добавляют силу давления:

u

 

 

 

 

n × m

t

+ (u

×Ñ)u

= nq (E + [u

´ B]) - Ñp,

 

 

 

 

 

где p = nkT.

30