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

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

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

n vx u0

us2

 

HyH̃y

 

 

 

u

 

 

 

 

 

 

4

 

 

 

и двух предыдущих уравнений находим:

 

2

1

 

 

1

 

 

Ht2

 

,

 

 

u02 us2

 

8

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

Ht2

8 u02

us2

 

 

(VII.2.9)

Поскольку

0,

то

 

в

 

волне,

распространяющейся со скоростью u0

u2

us,

магнитное поле

усиливается,

а

 

со

скоростью

u0 u3 us – ослабляется.

VII.3. Ударные волны в слабых полях

Рассмотрим случай, когда давление плазмы много больше магнитного

H2 us2.

(VII.3.1)

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

2

1

 

P2 P1

 

1

1

 

0,

(VII.3.2)

2

2

1

 

 

 

 

 

 

 

 

171

vt 0.

(VII.6.3)

Эти соотношения такие же, как в обычной газодинамике. Найдем скачок тангенциальной составляющей поля. Выражая vt из (VI.2.14) и подставляя в (VI.2.13), получаем:

2

 

Ht

 

 

Hn2

 

 

j

 

 

 

 

 

Ht .

VII.3.4)

 

4

 

 

 

 

 

 

Правая часть (VII.6.4) квадратична п о Hn, и ею можно пренебречь, т.е.

 

 

 

Ht

0 ,

(VII.3.5)

 

 

 

 

 

 

 

 

 

 

 

т.е.

 

 

 

 

 

 

 

 

 

 

Ht2

 

 

2

 

 

 

 

 

 

 

.

(VII.3.6)

 

 

Ht1

 

1

Так как 2

1,

 

наблюдается

усиление

магнитного поля

 

 

 

 

 

 

 

 

 

 

Ht2

Ht1 .

(VII.3.7)

Уравнения (VI.2.1 1 14) представляют и другую возможность: скачок плотности мал, а поток вещества через фронт – малая величина второго порядка

j2

Hn2

.

(VII.3.8)

 

 

4

 

С другой стороны, при 1 2 ударная адиабата превращается в том же приближении в

172

уравнение

P2 P1

 

1

2

 

2

 

 

 

Ht2

Ht1

.

(VII.3.9)

8

Условие 2 2 – это условие несжимаемости, т.е. уравнение (VII.6.9) - это ударная адиабата для несжимаемой жидкости. Тангенциальная составляющая магнитного поля не изменяется по направлению, но совершает скачок по величине. Скорость распространения разрыва, иначе фронта волны определяется так:

v2n v1n

 

j

 

Hn

,

(VII.3.10)

 

 

 

 

 

4

 

 

то есть равна альфвеновской скорости, вычисленной по нормальной составляющей поля. Так как в такой волне P2 P1, ослабляется тангенциальная составляющая магнитного поля:

Ht22 Ht12 .

(VII.3.11)

Рассмотрим теперь два частных случая волн произвольной магнитуды в произвольных полях.

Пусть перед фронтом волны магнитное поле перпендикулярно фронту, Ht1 0 . Тогда уравнение (VII.6.4) переходит в следующее:

j

2 Ht2

 

Hn2

Ht2.

(VII.3.12)

 

2

4

 

 

 

 

 

Пусть Ht2 0 и произвольно. Тогда

173

 

2

 

Hn2

 

2

Hn2

 

j

 

 

 

2,

т. е.

vn

 

.

(VII.3.13)

 

4

4 2

Пусть

 

теперь

Ht2

0,

т.е.

волна

распространяется параллельно полю, а ее фронт

полю перпендикулярен. В этом плоском случае

поле не влияет

на волну.

Второй

случай

– это

распространение волны поперек поля, Hn 0. Для этого случая u2 u3. Имеется только один тип волн.

Из уравнения

vt2 vt1

 

Hn

Ht2

Ht1

 

4 j

 

 

 

 

находим:

 

 

 

 

vt2

vt1 ,

(VII.3.14)

Выбором системы координат можно добиться vt 0. В этой системе вещество течет перпендикулярно фронту. Тогда из уравнения

(VII.6.4) находим:

 

 

 

 

 

 

 

 

Ht2

 

Ht1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

(VII.3.15)

 

 

 

 

 

 

 

 

 

2

1

Ввиду Hn

0 это уравнение дает:

 

 

 

 

 

 

 

 

 

 

 

H2

 

H1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

(VII.3.16)

 

 

 

 

 

 

 

 

 

2

1

Наряду с обобщенным давлением P P

H2

8

 

можно ввести обобщенную внутреннюю энергию

 

H2

,

т.к.

s P . В результате можно

 

1

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

174

написать:

2

1

 

P2 P1

 

1

2

,

(VII.3.17)

2

 

1 2

 

 

 

 

 

 

что и является ударной адиабатой для перпендикулярной волны.

VII.4. Генерация магнитного поля при турбулентном движении

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

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

Пусть вначале магнитное поле достаточно

175

слабое и не влияет на движение жидкости. Для простоты будем считать жидкость несжимаемой,const, divv 0. Движение тогда описывается уравнением Навье Стокса:

 

 

v

 

 

v, v

P

v.

(VII.4.1)

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь – коэффициент вязкости. Это уравнение

 

можно переписать так:

 

 

 

 

 

 

 

 

 

 

 

 

v

 

v2

v, , v

P

v.

(VII.4.2)

 

 

 

 

 

 

 

t

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

Возьмем rot от этого уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

rot v,

.

(VII.4.3)

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

 

мы

 

 

ввели

обозначение

 

1

rotv.

 

 

 

 

2

 

Cравним

 

 

это

 

 

 

уравнение

 

с

уравнением,

 

описывающим магнитное поле в проводящей

 

жидкости:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

rot v, H

 

c2

H.

(VII.4.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

4

 

 

 

 

Если величина v задана, то оно полностью описывает эволюцию магнитного поля.

Функции и H удовлетворяют одному и тому же уравнению, если 4c2 . В этом случае может существовать решение

H const .

(VII.4.5)

При стационарной (в среднем) турбулентности стационарно (в среднем) и магнитное поле H, т.е. в

176

среднем генерация и затухание уравновешены. Если4c2 , диссипативные потери будут меньше, чем

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

Эта аргументация не применима для двумерной турбулентности, например, в токамаках.

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

Глава VIII. Нелинейные волны в плазм

VIII.1. Простые волны

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

Для начала рассмотрим простейшую модель: пусть вдоль оси x движется неограниченный в пространстве пучок невзаимодействующих частиц с одинаковыми скоростями v. Уравнение, описывающее такой поток, имеет вид:

177

vt v vx 0.

.

Внесем теперь возмущение в теперь v v0 u, где u a sin kx. систему отсчета, движущуюся со Bместо уравнения (VIII.1.1) имеем:

(VIII.1.1)

поток. Пусть Перейдем в скоростью v0.

u

u

u

0.

(VIII.1.2)

t

x

 

 

 

Очевидно, что частицы, для которых u 0, будут догонять более медленные, движущиеся вправо. Вблизи точек, для которых sin kx 0, частицы будут скапливаться. Аналогичная картина будет и для частиц, у которых u 0. В результате профиль волн будет становиться все круче, а точки kx l l 0, 1, 2. . . будут точками сгущения плотности (см. рис.25). В дальнейшем это приводит к "опрокидыванию фронта".

Это называется "бунчировкой" частиц.

При изменении формы волны возникают высшие гармоники ряда Фурье. В нулевом приближении u 0 sin kx. Считая квадратичные члены малыми, в первом приближении получаем:

 

u 1

 

k 2 sin kx cos kx uiv

k 2

sin 2kx. (VIII.1.3)

 

t

 

 

 

 

 

 

2

 

Скорость u 1 растет линейно со временем.

Спустя время

1

первая

поправка u 1

k

становится масштаба нулевой, и метод теории возмущения перестает быть применимым.

178

Рис.25. Укручение фронта волны

VIII.2. Нелинейное распространение ионного звука

Рассмотрим теперь нелинейное распространение ионного звука. Для простоты будем считать Ti 0, а электронную теплопроводность устремим к бесконечности, т.е. положим Te const. Тогда в отсутствие магнитного поля уравнения, описывающие плазму, принимают вид:

n

 

nv 0 ,

(VIII.2.1)

t

x

 

 

 

179

 

v

v

v

 

 

Te

 

 

n

 

0

 

 

(VIII.2.2)

 

t

x

 

min

x

 

 

 

 

 

 

 

 

 

 

 

 

Здесь мы учли, что при Te

const

 

P

Te

n

.

 

x

x

Выделим из большого класса решений

интересующий нас класс – простые волны.

 

 

 

Если n и v – однозначные функции координат,

то n можно выразить через v :

v v n .

Тогда

уравнения (VIII.2.1) и (VIII.2.2) перепишем так:

 

dn

v

v

v

 

n

v

 

0 ,

t

x

 

x

 

v

 

 

 

 

 

 

 

 

 

v

 

 

v

2

1

dn v

 

 

 

 

v

 

cs

n

dv

 

 

0 .

 

t

x

x

Здесь мы обозначили cs2

 

 

Te/mi

 

Умножим второе из этих уравнений на вычтем из него первое:

(VIII.2.3)

(VIII.2.4)

dndv и

cs2

 

dn

 

 

2 dv

 

 

dv

 

 

n

 

dv

 

 

dx

 

n dx

0 ,

(VIII.2.5)

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

dn

 

2

n

2

.

(VIII.2.6)

 

 

cs

 

dv

 

 

 

Отсюда находим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cs dn

n.

 

 

(VIII.2.7)

 

 

 

 

 

 

dv

 

 

 

 

 

Подставляем это выражение в последний член уравнения (VIII.2.4):

180

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