
Морозов Введение в теорию горячей плазмы Ч.1 2011
.pdf
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