слабое и не влияет на движение жидкости. Для простоты будем считать жидкость несжимаемой,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