ФИЗИКА ПЛАЗМЫ, 2015, том 41, № 3, с. 286-293
НЕЛИНЕЙНЫЕ ЯВЛЕНИЯ
УДК 533.9
ПАРАМЕТРИЧЕСКИЙ РАСПАД НЕОБЫКНОВЕННОЙ ЭЛЕКТРОМАГНИТНОЙ ВОЛНЫ В РЕЛЯТИВИСТСКОЙ ПЛАЗМЕ © 2015 г. В. Г. Дорофеенко, В. Б. Красовицкий*, В. А. Туриков**
Институт передовых исследований, Вена, Австрия * Институт прикладной математики им. М.В. Келдыша РАН, Москва, Россия ** Российский университет дружбы народов, Москва e-mail: krasovit@mail.ru Поступила в редакцию 22.04.2014 г. Окончательный вариант получен 03.07.2014 г.
Рассмотрена параметрическая неустойчивость необыкновенной электромагнитной волны в плазме, предварительно нагретой до релятивистской температуры. Получена и исследована автомодельная система нелинейных уравнений в полных производных, учитывающая "тепловую" массу электронов. Малые возмущения параметров нагретой плазмы проанализированы на основе дисперсионного уравнения, определяющего фазовые скорости быстрой и медленной необыкновенных волн в линейном приближении. В отличие от холодной плазмы зона непрозрачности в области частот, превышающих верхнегибридную частоту электронов, исчезает, и асимптоты обеих ветвей сближаются. Теоретический анализ системы нелинейных уравнений показал, что инкремент распадной неустойчивости растет с увеличением начальной температуры электронов плазмы. Этот результат на качественном уровне подтверждается выполненным в работе численным моделированием процесса нагрева плазмы инжектированным из вакуума лазерным импульсом.
DOI: 10.7868/80367292115030026
1. ВВЕДЕНИЕ
В последние годы процесс взаимодействия мощных лазерных импульсов с плазмой интенсивно исследуется экспериментально и теоретически (см., например, обзор [1]). Такое взаимодействие имеет множество практических приложений от ядерной физики до медицины.
При взаимодействии лазерных импульсов ультрарелятивистской интенсивности с плазмой генерируется квазистатическое магнитное поле с индукцией свыше 100 МГс. Вблизи границы плотной плазмы образуется более разреженная плазма, к которой приложено сверхсильное магнитное поле. В этих условиях лазерная волна может распространяться в окрестностях электронного циклотронного и верхнегибридного резонансов, где присутствие сильного магнитного поля существенно изменяет картину взаимодействия электромагнитного излучения с плазмой [2].
Неустойчивость холодной плазмы в поле внешней электромагнитной волны с постоянной амплитудой, распространяющейся под произвольным углом к магнитному полю, аналитически и численно исследована в работе [3]. Установлено, что оптимальные условия для нагрева электронов плазмы возникают при распространении лазерного импульса поперек магнитного поля. Обратное воздействие возбуждаемых коле-
баний на волну учтено в работах [4, 5], где получена и исследована самосогласованная система нелинейных уравнений, состоящая из уравнения огибающей электромагнитной волны и уравнения для потенциала в плазме. Показано, что в первоначально холодной плазме имеет место распад необыкновенной электромагнитной волны на две волны половинной частоты, причем эффективность процесса зависит от величины инкремента неустойчивости, определяемого самосогласованной системой линеаризованных уравнений и пропорционального амплитуде волны накачки [4]. На нелинейной стадии неустойчивости экспоненциальный рост амплитуды вторичной волны с половинной частотой сменяется обратным процессом возвращения энергии в первичную волну и возникновением нелинейных колебаний в плазме поперек внешнего магнитного поля [5].
Одним из перспективных направлений использования предварительно нагретой плазмы для реализации УТС является Z-пинч и сопряженный с ним лазер. Роль последнего предполагает дополнительный нагрев уже сформировавшегося плазменного шнура с последующей транспортировкой плазмы на термоядерную мишень [6]. Для эффективного преобразования энергии лазерной волны в энергию электронов
плазмы может быть использована распадная неустойчивость [4, 5]. В настоящей работе параметрический распад необыкновенной электромагнитной волны на две волны половинной частоты рассмотрен в плазме с релятивистской температурой электронов [7].
В разд. 2 система уравнений Максвелла и уравнений гидродинамики с градиентом давления сведена к самосогласованной системе автомодельных уравнений в полных производных, описывающей распространение нелинейной волны в нагретой плазме.
Дисперсионное уравнение, определяющее фазовую скорость волны в линейном приближении, получено в разд. 3. Показано, что в отличие от холодной плазмы [5] непрозрачность плазмы на частотах, превышающих верхнегибридную частоту, исчезает.
Вывод формулы для инкремента параметрической неустойчивости необыкновенной электромагнитной волны в нагретой плазме представлен в разд. 4.
Эффективность нагрева плазмы инжектированным из вакуума лазерным импульсом подтверждается выполненным в разд. 5 численным моделированием.
2. СИСТЕМА НЕЛИНЕЙНЫХ УРАВНЕНИЙ
Система гидродинамических уравнений плазмы с релятивистской температурой электронов имеет вид [7]
I+V. 11 ю <«,, ]=.
Ex -
(Ву + Bo) V .
| - V. |) ю (5) Р. ] =
= е
Е. +
(Ву + Во) х
Т А
п д.
г \
п
IУ У
(1)
дп + . ) = о, дt д/
где
Рх,. = т^ х., 2
о (5)
5 =
тс Т
с
У =
V
1 - V-
2
с
= £з(5)
К2(5)
п-1/2
электронов Т приводит к появлению в уравнениях (1) "тепловой" массы [7]
те/1 = тО (5) =
5 Т
т |1 + I, Т < тс2,
2 тс ) (2)
С 2 А (2)
4Т
2
1 +
т с
8Т2
Т > тс2.
При этом компоненты самосогласованного электромагнитного поля в формулах (1)
Е =-1 дА*
х с д ? '
В = дАх
Ву =17,
(3)
(Ах — векторный потенциал и Е. — продольное электростатическое поле) определяются системой уравнений
' 1 ^ - Л
^с2 дt2 д.2
Ах = 4пе nv х,
дЕ. . , ч
—- = 4пе (п - п0).
д.
(4)
Из уравнений (1) и (4), следуют интегралы
Е. = 4пепо (— о - -), (5)
О(5)Рх + -Ах + тюв ( - .о) = о, (6)
с
где . о — координата электрона с рх = о в отсутствие поля, Е. = Ах = о.
Будем искать решение в виде бегущей волны, считая все величины зависящими от безразмерной фазы
¥ = ю р
? ■
V
рИ J
В безразмерных переменных еА
А =■
тс
е = ,
ти рс
а = у I 1
в V.
и р
( о - .)
с
исходная система уравнений (1) и (4) упрощается к виду
а- = -£- £, с=-2Чва - 4
в2 -в2 -\\4Я I
где
т, V, р — масса, скорость и импульс электрона, п — плотность, с — скорость света, Во — внешнее поперечное магнитное поле, Кп (5) — функция Макдональда. Учет релятивистской температуры
а = 4Я ^ р-С'Ш'
Я = а2 + (р2 -1)(1 + 02),
0 = А - ОвС, Р = Vри /с,
Ов = юв/Ир,
рт = (Т / т^с 2)1/2 =[5О (5)]-1/2,
(7)
Рис. 1. Фазовая скорость в(ш) для быстрой (БН) и медленной (МН) необыкновенных волн. Пунктирными линиями показаны условия параметрического резонанса в г = в(ш г). Параметр 0В = 0.67. Величина вт = 0.1.
причем, в силу формулы (2), величина вт изменяется в пределах от 0 при Т = 0 до 1/2 при Т = ^.
3. ДИСПЕРСИОННОЕ УРАВНЕНИЕ
Для малых возмущений А(/, ¿) = = А0 ехр I (шt - kz) дисперсионное уравнение линеаризованной системы уравнений (7) имеет вид
± = 1 _ 1 1 -ш2(1 -в2/в2) (8)
в2 ш21 + оВ -ш 2(1 _в2/в2)'
где ш — частота волны (в единицах ш р), в = ш/ ск — безразмерная фазовая скорость.
Зависимость в(ш) для двух значений вт при
фиксированном 0В представлена на рис. 1 и рис. 2. Условия параметрического распада быст-
Рис. 2. То же, что на рис. 1, для вт = 0.5.
рой необыкновенной (БН) волны с частотой шг на две медленные необыкновенные (МН) волны с половинной частотой шг /2 показаны пунктирными линиями. Границы области прозрачности холодной плазмы ш£ и шя для БН- и МН-волн (где в ^ ^) определяются формулой
= 2(а/)В + 4 + )в) (см. [5]).
ш
Условие параметрического резонанса ш = шг и ш' = ш г /2, когда частота возбуждаемой волны ш' равна половине частоты волны накачки ш, требует выполнения равенства (8) как для частоты накачки шг, так и для частоты вторичной волны шг /2. Можно показать, что для неизменной фазовой скорости волны в это условие реализуется только в ограниченной области значений внешнего магнитного поля 0В = 0Вг (шг, в т) и в = в г (ш г, вт \ где
)вг = ^ (5 - ш2(1 - в2)))(1 - в2) + д/ш4 (1 - в2)2 + ЩГ) -1 - в
2
т,
(9)
в г ш 2(1 + в2) + Уш44 (1 -в2)2 + 16 в2. (ю)
2(ш4 - 4)
Графики этих функций представлены на рис. 3. Как следует из этого рисунка, неустойчивость имеет место только в области резонансных частот
42
< ш г < шс
4 -в2
¡1 -вт
закрашенной на рисунке серым цветом. При этом
0 < 0
Вг
< 3,
4 -вт
<в г <
4. ИНКРЕМЕНТ НЕУСТОЙЧИВОСТИ
Метод расчета инкрементов распадных не-устойчивостей волн в плазме в общем виде изложен в обзоре [8]. Поэтому ниже мы приведем лишь краткое описание наших вычислений, аналогичных выполненным в [4] для холодной плазмы, необходимое для воспроизводимости полученных результатов. Линеаризация системы (7) вблизи линейной волны накачки, представленной вектором
У = (А, А', С, а) = 2 (ВД^ + У,*(ю)е -йв*),
дает:
х' =
^(0) + гу(0) • У (V)
х(^),
(11)
где ю — частота волны накачки, вектор-столбец Г (У) есть
Л в3
в2 -1-Тя Г (у) = - 1
в2 - 1\л/д
. р(= р-С'Ш
а х(у) представляет собой малое возмущение. Мы
не приводим явный вид матриц Г/(0) и ГУ'(0) из-за их громоздкости.
Ищем решение системы (11) методом вариации постоянных с последующим усреднением по периоду колебаний, полагая х = Z • а, где Z(у) — 4 х 4-матрица линейно независимых решений
уравнения Z' (у) = ГУ (0) • Z(у). В результате получаем:
где
а' = М • а,
М = ЦZ-1 (Ое/юу + О*е"/юу IZ
(12)
2п/ю
(/) = ю / / Шу,
z =
ю-А
У1 (ю)е/юу, У1 (-ю)е",юу, У1
ю-А
у
ю-А
У1
ю
у
где малая величина А = ю - 2ю' определяет расстройку параметрического резонанса, а вид функции У1 (ю)
1
1/У2
1 1 1 1 — Qвr
* \ --Рг
*
%
\ч
Л
^2 юг
Рис. 3. Внешнее магнитное поле Qвr и фазовая скорость волны Рг в зависимости от резонансной частоты юг. Величина Рт = 0.3.
1
/ю
(Р2 - 1)ю^
У1 (ю) =
ю2 (р2-Р2 )-р2
/Р(Р2 - 1)ю^д ю2 (Р2-Р2 )-р2
следует из линейной теории. Удерживая в М лишь резонансные слагаемые, имеем:
(13)
о = гу (0)У1.
Так как матрица О пропорциональна малой амплитуде У1, то зависимость а от свободной переменной является медленн
Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.