первого члена берём по частям:
dpxdpydpzpx ∂∂pfx vyHz − vzHy −n VyHz − VzHy .
Другие компоненты уравнения (V.1.1) вычисляются аналогично. Окончательный результат, переписанный в векторной форме, даёт:
v
v, H ∂∂pf
dp n V, H .
V.4. Одножидкостная гидродинамика
Уравнения (V.1,7), (V.2.18), |
и (V.3.1), |
написанные по отдельности для электронов и разных сортов ионов, представляют собой уравнения многожидкостной гидродинамики. Рассмотрим пока более простой случай двухжидкостной гидродинамики, когда плазма состоит из электронов и одного сорта ионов. Для простоты будем считать, что ионы однозарядные, Z 1. В этом случае можно перейти к более простой системе одножидкостной гидродинамики.
Для электронов уравнение (V.2.18) перепишется
так:
mene |
∂ |
Vje |
∂ |
Vie − |
∂Pe |
− |
(V.4.1) |
∂t |
∂xj |
|
|||||
|
|
|
∂xi |
|
|||
136
|
|
1 |
Ve, H i |
|
ji |
|
∂ ije |
|
− ene |
Ei |
ene |
|
− |
|
− RTi. |
||
c |
|
|
||||||
|
|
|
|
|
∂xj |
|||
Мы заменили разность скоростей ионов и
электронов через плотность тока j ne Vi − Ve и ввели проводимость плазмы
|
|
|
|
|
|
|
|
|
|
e2ne |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
me ei |
. |
|
|
|
|
|
|
|
(V.4.2) |
|||
Здесь ei |
|
|
|
− частота электрон-ионных |
|||||||||||||||
столкновений. Для ионов аналоничное уравнение |
|||||||||||||||||||
принимает вид: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
mini |
|
∂ |
|
|
Vji |
∂ |
|
Vii − |
∂Pi |
|
|
(V.4.3) |
|||||||
|
∂t |
∂xj |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
∂xi |
|
|
|||||||
eni Ei |
1 |
|
i |
, H i |
|
− eni |
ji |
|
|
∂ iji |
RTi. |
||||||||
|
c |
|
V |
|
|
|
− |
|
|||||||||||
|
|
|
|
∂xj |
|||||||||||||||
Введём массовую скорость V |
|
meVe miVi |
. Если |
||||||||||||||||
|
me mi |
|
|||||||||||||||||
упорядоченная скорость электронов не слишком превосходит упорядоченную скорость ионов, то массовая скорость практически совпадает с последней. Примем во внимание также, что mi mi, условие квазинейтральности ne ni, пренебрежём вязкостью и сложим уравнения (V.4.1)
и (V.4.3):
137
min |
dV |
−P |
1 |
j, H |
(V.4.4) |
|
dt |
|
c |
||||
Здесь мы ввели |
полное |
|
давление |
плазмы |
||
P Pi Pe.
Уравнение (V.4.4) можно переписать несколько
иначе, используя уравнение rotH 4c j.
min ddtV −P 41 , H , H .
Раскрывая двойное векторное произведение в правой части, получаем:
min |
dV |
|
|
− P |
H2 |
|
1 |
|
H, H. |
(V.4.5) |
|||
dt |
8 |
4 |
|||||||||||
|
|
|
|
|
|
|
|
|
|||||
Величину |
H2 |
называют магнитным давлением, а |
|
||||||||||
8 |
|
||||||||||||
величину P |
|
|
H2 |
|
полным давлением. |
|
|
||||||
|
|
8 |
|
|
|||||||||
Это уравнение должно быть дополнено уравнением для тока. Воспользуемся для этого уравнением (V.4.1). Оценим члены в левой и правой частях. Пусть характерное время изучаемого процесса . Тогда первый член в левой части имеет величину порядка nmeVe/ . Первый член в правой части имеет порядок nTe/a nmev2Te/a, где a -
характерный размер задачи. Второй член в левой части имеет порядок men Ve 2/a. Очевидно, что
если Ve vTe и Ve v2Te/a, левой частью в уравнении (V.4.1), т.е. инерцией электронов можно
пренебречь. В результате получаем обобщённый
138
закон Ома:
1 |
∂Pe |
|
− |
Ei |
|
1 |
Ve, H i |
|
(V.4.6) |
||||
ene |
∂xi |
|
c |
||||||||||
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
ji |
|
|
∂ ije |
− RTi |
|
|
||||
|
|
|
|
|
∂xj |
|
|
|
|||||
|
|
|
− |
|
ene |
. |
|
|
|||||
|
|
|
|
|
|||||||||
Электронная вязкость во многих задачах
пренебрежимо мала. В металлических проводниках градиент электронного давления пренебрежимо мал (чего нельзя сказать о горячей плазме), а относительная скорость электронов и ионов практически равна электронной скорости. Поэтому в металлическом проводнике уравнение (V.4.6) переходит в обычный закон Ома с учётом эффекта Холла.
Из уравнения (V.4.1) получаем одно из важнейших уравнений в теории магнитного удержания − условие равновесия. В стационарном случае производная по времени обращается в ноль. Для неподвижной плазмы уравнение равновесия принимает вид:
P |
1 |
j, H . |
(V.4.7) |
c |
Во многих случаях плазма в термоядерных установках вращается. Инерционный член min V, V оценивается так же, как и для электронов. Поэтому если скорость вращения плазмы много меньше тепловой скорости ионов, то уравнение (V.4.6) пригодно для расчёта равновесия
139
и в этом случае. Следует, однако, заметить, что ионы примеси часто значительно массивнее ионов основной плазмы, и при расчётах следует учитывать их инерцию.
V.5. Теорема вмороженности
Рассмотрим бесстолкновительную плазму, скорость движения которой много меньше тепловой скорости ионов. Тогда в уравнении (V.4.3) можно пренебречь инерцией ионов. Пренебрежём также вязкостью и термосилой. Тогда плотность тока в плазме выражается в виде:
j |
Pi |
− |
Ei |
1 |
V, H i . |
(V.5.1) |
ne |
c |
В металлическом проводнике Pi 0. Возьмём rot от уравнения Максвелла rotH 4c j
и подставим туда (V.5.1): |
|
|
|
|
|
|
|
||||||||
|
c2 |
|
H −c , E |
|
|
c |
, E . |
||||||||
|
4 |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
В металлических |
проводниках |
обычно |
|||||||||||||
сonst. Используя |
ещё |
одно |
уравнение |
||||||||||||
Максвелла, rotE − |
1 |
|
∂H |
, в результате получаем: |
|||||||||||
c |
∂t |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
∂H |
− rot V, H |
c2 |
|
|
H. |
(V.5.2) |
|||||||
|
|
∂t |
4 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||
Правая часть уравнения (V.5.2) описывает диффузионное расплывание неоднородностей
140