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

Морозов Введение в теорию горячей плазмы Ч.1 2011

.pdf
Скачиваний:
28
Добавлен:
12.11.2022
Размер:
1.82 Mб
Скачать

магнитного поля в проводящей среде за счёт конечной проводимости или диффузию переменного поля внутрь проводящего слоя с

границы.

Величину

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, то вектор дополнительного элемента поверхности можно представить в виде dSv, 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

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