Морозов Введение в теорию горячей плазмы Ч.1 2011
.pdfмагнитного поля в проводящей среде за счёт конечной проводимости или диффузию переменного поля внутрь проводящего слоя с
границы. |
Величину |
Dm |
c2 |
называют |
4 |
коэффициентом диффузии магнитного поля в проводнике. При достаточно высокой проводимости можно положить
nePi E ≡ Ei 1c V, H i .
В металлических проводниках с большой проводимостью это даёт:
E 0. |
(V.5.3) |
В плазме градиент давления не обязан быть
равным нулю, и некритическое использование выражения (V.5.3) может привести к ошибкам.
В металлическом проводнике при высокой проводимости диффузионный член в (V.5.2) перестаёт играть роль, и можно пользоваться уравнением идеальной гидродинамики
∂H |
rot V, H . |
(V.5.4) |
|
∂t |
|||
|
|
Вплазме в этом уравнении появится
дополнительный член: |
|
|
|
||
|
∂H |
rot V, H − |
1 |
n, Pi . |
(V.5.5) |
|
∂t |
n2e |
|||
|
|
|
|
141
В этом разделе ограничимся тем случаем, когда
достаточно ограничиться уравнением (V.5.4). Это условие приводит к "вмороженности" плазмы в магнитное поле. Рассмотрим произвольный замкнутый контур L , движущийся вместе с плазмой, и введём магнитный поток H dS
S
через этот контур. Пусть за малый отрезок времени контур L превратился в контур L′. Тогда изменение потока в первом приближении можно представить как изменение потока через старый контур за счёт изменения магнитного поля плюс изменение потока за счёт изменения площади контура S при неизменном магнитном поле:
t |
∂H |
dS |
H dS′. |
(V.5.6) |
∂t |
||||
S |
|
S |
|
|
Так как элемент контура движется со скоростью плазмы v, то вектор дополнительного элемента поверхности можно представить в виде dS′v, dL t, а интеграл (V.5.6) в виде (см. рис. 23)
t |
|
∂H |
dS H, v dL . |
|
∂t |
||||
|
S |
L |
|
|
Используя теорему |
Стокса |
и уравнение |
(V.1.40), получаем:
142
t |
|
∂H |
rot v, H dS |
0, |
(V.5.7) |
∂t |
|||||
|
S |
|
|
|
|
Рис. 23. Изменение контура интегрирования
т.е. магнитный поток через любой замкнутый контур, движущийся вместе с идеально проводящей плазмой (с учётом сделанных выше предположений) не меняется. Магнитные силовые линии как бы приклеены к плазме и "отлипают" только в меру нарушения сделанных предположений. Это обстоятельство сильно облегчает качественное рассмотрение многих процессов.
143
Глава VI. Линейные МГД-волны и течения в идеально проводящей плазме.
VI.1. МГД-волны малой амплитуды в холодной плазме
В предыдущих главах были рассмотрены волны малой амплитуды в незамагниченной плазме. В замагниченной плазме появляется выделенное направление, связанное с внешним магнитным полем, поэтому спектр возмущений становится значительно более сложным. Рассмотрим малые возмущения в однородной неподвижной холодной плазме. С одной стороны, пусть частота этих возмущений достаточно велика по сравнению с частотой столкновений, т. е. плазму можно считать идеально проводящей. С другой стороны, она должна быть достаточно мала, чтобы можно было пренебречь током смещения в уравнениях Максвелла. Пусть невозмущенная плазма помещена в однородное магнитное поле H0. Невозмущенные величины плотности, давления и магнитного поля будем обозначать индексом "0", а возмущенные волнистой чертой наверху, например H H0 H̃. При этом будем считать возмущенные величины малыми первого порядка величины. Пренебрегая величинами второго порядка малости,
144
линеаризованную систему уравнений, описывающих плазму в МГД − приближении, можно представить в виде:
|
|
|
|
|
|
|
̃ |
|
|
|
|
|
(VI.1.1) |
|
|
|
|
divH 0, |
|
|
|
|
|||||
|
|
|
|
̃ |
|
|
|
|
|
|
|
|
|
|
|
|
|
∂H |
rot v. H0 , |
|
|
|
(VI.1.2) |
||||
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
∂t |
|
|
|
|
|
|
|
|
|
|
|
|
∂ñ |
|
div n0ṽ 0, |
|
|
(VI.1.3) |
|||||
|
|
|
∂t |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|||
̃ |
|
1 |
|
|
|
1 |
|
|
|
|
|
||
∂v |
|
|
|
̃ |
|
̃ |
|
. |
(VI.1.4) |
||||
∂t |
− n0mi |
|
|
||||||||||
p |
4 min0 |
rotH, H0 |
|
Будем искать решение системы уравнений (VI.1.1-4) в фурье-представлении
|
|
̃ |
|
|
(VI.1.5) |
|
k, H 0, |
|
|||
|
̃ |
|
|
|
(VI.1.6) |
H H0 k, ṽ , |
|||||
ñ |
|
k, ṽ |
, |
(VI.1.7) |
|
|
n0 |
|
min0 ṽ kmiu2s ñ − 41 H̃ k, H0 − k H0, H̃ . (VI.1.8)
Здесь мы выразили возмущенное давление p̃ через возмущение плотности ñ; считая движение изэнтропическим, s const:
145
̃ |
|
∂p |
|
2 |
|
|
p |
|
|
|
sñ mus |
ñ , |
(VI.1.9) |
|
∂n |
|
где us – скорость звука в незамагниченной среде. Исключаем ñ с помощью уравнения (VI.1.7)
2 |
k, ṽ |
|
̃ |
− k H0 |
̃ |
|
|
H k, H0 |
, H |
|
|||
v kus |
|
− |
|
|
, |
|
|
4 min0 |
|
H̃ −ṽ k, H0 H0 k, ṽ ,
|
|
|
̃ |
|
|
|
|
|
k, H 0. |
|
|
Выберем |
ось |
x |
вдоль |
направления |
|
распространения волны k, k k, 0, 0 |
̃ |
||||
Из |
уравнения (VI.1.12) |
видно, что H k т.е. |
|||
̃ |
|
|
|
|
|
0, Hy , |
Hz . |
|
|
|
|
H |
|
|
|
Выберем плоскость xy так, чтобы она проходила
через k и H0: H0 H0x, H0y, 0
При этом легко видеть, что z-компоненты уравнений (VI.1.10-11) содержат только z-компоненты всех возмущенных векторов, а уравнение (VI.1.12) удовлетворяется тождественно. Тогда z-компонента системы (VI.1.10-12) отщепляется и превращается в систему двух линейных алгебраических уравнений:
(VI.1.10)
(VI.1.11)
(VI.1.12)
146
|
H0x |
̃ |
̃ |
|
uvz |
4 mn0 |
Hz 0, |
H0xṽz uHz 0. |
(VI.1.13) |
Здесь мы ввели фазовую скорость u k .
Эта система уравнений имеет нетривиальное решение, если ее определитель обращается в ноль.
u |
H0x |
|
|
|
4 mn0 |
|
0, |
(VI.1.14) |
|
|
|
|||
H0x |
u |
|
|
|
|
|
|
|
|
или
u1 |
|
H02x |
vA |
H0x |
. |
(VI.1.15) |
|
4 min0 |
H0 |
||||||
|
|
|
|
|
Здесь мы ввели альфвеновскую скорость
vA |
H0 |
. |
(VI.1.16) |
|
|||
|
4 min0 |
|
Знак означает, что волна может распространяться как вдоль направления оси x, так и против него. Угол между осью x и полем H0 может быть произвольным. При увеличении угла фазовая скорость уменьшается от альфвеновской до нуля при распространении перпендикулярно магнитному полю. Возмущение магнитного поля H̃ перпендикулярно как невозмущенному магнитному полю, так и направлению распространения волны, т.е. направлено вдоль оси z. Возмущение плотности пропорционально x-составляющей скорости:
147
|
ñ |
n0 |
ṽx |
, |
|
|
(VI.1.17) |
||
|
u |
|
|
|
|||||
а скорость смещения |
|
плазмы ̃ |
вдоль оси |
z |
|||||
пропорциональна возмущению поля Hz: |
|
|
|||||||
vz − |
|
|
̃ |
. |
|
|
(VI.1.18) |
||
|
Hz |
|
|
||||||
|
|
|
4 min0 |
|
|
|
|||
Уравнение (VI.1.15) можно переписать в |
|||||||||
произвольной системе координат: |
|
|
|
||||||
|
H0xkx |
|
|
H0, k |
, |
(VI.1.19) |
|||
|
|
|
4 min0 |
||||||
4 min0 |
|
|
|
|
|
||||
откуда находим групповую скорость волны |
|
||||||||
ugr |
|
∂ |
|
|
H0 |
. |
|
(VI.1.20) |
|
∂k |
|
|
|||||||
|
|
|
|
|
4 min0 |
|
|
|
.
Очевидно, что групповая скорость направлена строго вдоль невозмущенного магнитного поля.
Рассмотрим теперь остальные компоненты системы (VI.1.10-12). Y–компонента уравнения
(VI.1.10) дает:
|
H0y |
̃ |
|
̃ |
|
Hy. |
(VI.1.21) |
|
|||
uvy − 4 min0 |
Y– компонента уравнения (VI.1.11) дает:
148
̃ |
̃ |
̃ |
(VI.1.22) |
Hyu vxH0y − vyH0x. |
X – компонента уравнения (VI.1.10) дает:
|
us2 |
H0y |
̃ |
|
|
vx u − |
u |
|
4 min0 |
Hy. |
(VI.1.23) |
X – компонента уравнения (VI.1.11) дает тождественный ноль.
Для получения дисперсионного уравнения определитель системы (VI.1.21-23) следует приравнять к нулю. В результате получаем биквадратное уравнение
|
2 |
2 |
|
2 |
2 H02x |
|
2 |
2 H02y |
|
|
||
u |
|
− us |
u |
|
− vA |
|
u |
|
vA |
|
, |
(VI.1.24) |
|
|
H02 |
|
H02 |
корни которого имеют вид:
|
|
1 |
|
|
2 |
2 |
|
|
H0x |
|
|
||||
u2,3 |
|
|
|
|
us |
vA |
2vAus |
|
|
|
|
|
(VI.1.25) |
||
2 |
H0 |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
2 |
|
2 |
|
H0x |
|
|
|
|
||
|
|
|
|
us |
vA − 2vAus |
|
|
. |
|
|
|||||
|
|
|
H0 |
|
|
Колеблются x и y-компоненты скорости и y– компонента магнитного поля. При этом векторы H̃ и ṽ лежат в плоскости H0k.
Проанализируем полученные ветви колебаний в пределе большого и малого отношений магнитного
149
давления и давления плазмы. Рассмотрим сначала случай
H02 |
|
min0us2 |
. |
(VI.1.26) |
|
8 |
2 |
||||
|
|
|
Вэтом случае волна u2 ≈ us, т.е. превращается
вобычный звук. Волна u3 превращается в волну с фазовой скоростью.
u3 |
vA |
H0x |
. |
(VI.1.27) |
|
||||
|
|
H0 |
|
В несжимаемом случае us → , и остается только один тип волны, который и исследовался Альфвеном. В этой волне
|
̃ |
|
|
̃ |
H |
|
|
4 min0 . |
(VI.1.28) |
||
v − |
Очевидно, что в этом случае v H0, т.е. в гидродинамическом смысле это строго поперечная волна. Рассмотрим теперь противоположный случай, т.е.
v |
2 |
2 |
|
H02 |
|
2T |
|
|
u |
, |
|
|
|
. |
|
|
4 min0 |
mi |
|||||
|
A |
s |
|
|
|
Это равнозначно малости параметра
|
2n0T |
|
us |
|
2 |
1. |
||
|
|
vA |
|
|||||
|
H02 |
|
||||||
|
|
8 |
|
|
|
|
|
|
150