ФИЗИКА ПЛАЗМЫ, 2015, том 41, № 4, с. 374-380
ДИНАМИКА ПЛАЗМЫ
УДК 533.951
ДИНАМИКА ПОНДЕРОМОТОРНОГО УСКОРЕНИЯ ИОНОВ В ЛАЗЕРНО-ПЛАЗМЕННОМ КАНАЛЕ © 2015 г. В. Ф. Ковалев****, В. Ю. Быченков**
* Институт прикладной математики им. М.В. Келдыша, РАН, Москва
** Физический институт им. П.Н. Лебедева, РАН, Москва *** Всероссийский научно-исследовательский институт автоматики, Центр фундаментальных и прикладных исследований, Росатом, Москва e-mail: vfkvvfkv@gmail.com, bychenk@sci.lebedev.ru Поступила в редакцию 28.08.2014 г.
Построено аналитическое решение задачи Коши для кинетического уравнения для ионов, описывающее их радиальное ускорение под действием пондеромоторной силы лазерного пучка, распространяющегося каналированно в прозрачной плазме. Для осесимметричной геометрии лазерно-плазменного канала получены временные и пространственные зависимости функции распределения ионов и найдены их интегральные характеристики, такие как плотность, средняя скорость и энергетический спектр. Для квазистационарного лазерного пучка аналитически описаны формирование пика плотности у границы канала и эффект опрокидывания ионного потока.
DOI: 10.7868/80367292115040010
1. ВВЕДЕНИЕ
Среди различных схем генерации ионов высоких энергий под действием мощного лазерного импульса широко обсуждается их радиальное ускорение из лазерно-плазменного канала, формирующегося при распространении света в прозрачной плазме докритической плотности [1—5]. Использование лазерно-ускоренных ионов из низкоплотной плазмы интересно с точки зрения получения нейтронов и различных радионуклидов, особенно медицинских [1, 6—8].
Динамика радиального ускорения ионов лазерным излучением в низкоплотной плазме (например, газовая струя [1], испаренная фольга [9], нанопористый углерод [10]) в значительной степени обусловлена силовым (пондеромоторным) воздействием лазерного импульса на электроны плазмы, которые выталкиваются пондеромотор-ной силой из канала, в результате чего возникает сильное радиальное электрическое поле разделения зарядов, ускоряющее ионы. Характерным примером служит задача об ускорении ионов, рассмотренная в работах [2, 5] и показавшая, что аналитическое описание ускорения частиц лазерным импульсом из плазменного канала, формирующегося в результате самофокусировки, представляет собой непростую задачу даже для приближенных подходов. Этим обусловлено, что существующие теоретические модели [2, 3, 5] лишь весьма качественно характеризуют радиальное ускорение частиц из лазерно-плазменного
канала. Соответственно, изучение пространственно-временного распределения лазерно-ускоренных частиц проводится большей частью с использованием кинетического численного моделирования, как правило, методом "частица-в-ячейке" (PIC) [11—13]. Определенное упрощение вносит численная одномерная электростатическая модель пондеромоторного ускорения, описывающая динамику радиального разлета плазмы под действием заданного поля [2, 5], когда учитывается только медленная динамика электронов плазмы, что отвечает усреднению по их быстрым осцилляциям в лазерном поле. Поскольку основные результаты, полученные для электростатической пондеромоторной модели базируются на численном PIC-моделировании, это затрудняет предсказание зависимостей пространственно-временных и спектральных характеристик ускоренных частиц из лазерного канала от произвольных параметров лазера и плазмы. Проведенное численное моделирование эффекта ускорения ионов под действием электрического поля, обусловленного пондеромоторной силой, показало формирование сильного скачка плотности плазмы на цилиндрической поверхности на границе лазерного канала — ярко выраженной плотной стенки канала [2, 5]. Это сопровождается сильным локальным нагревом электронов, последующим опрокидыванием стенок канала и инверсией поля разделения заряда [5].
В результате численного моделирования [2, 5] было также установлено, что электроны адиаба-
тически откликаются на пондеромоторное воздействие лазера, так что электрическое поле в канале с хорошей степенью точности определяется условием баланса электростатической силы, возвращающей электроны в канал и выдавливающей их пондеромоторной силы. Поэтому можно надеяться, что пока такой баланс выполняется, динамика ускорения ионов может быть исследована строго аналитически с помощью кинетического уравнения для ионов с заданным электростатическим полем. Такое исследование и является основной целью нашей работы.
2. КИНЕТИЧЕСКАЯ МОДЕЛЬ С ЗАДАННЫМ УСКОРЯЮЩИМ ЭЛЕКТРИЧЕСКИМ ПОЛЕМ
Рассматривается распространение цилиндрически-симметричного лазерного пучка в прозрачной плазме. Неоднородность электрического поля лазера по радиусу приводит к вытеснению электронов плазмы из области действия сильного электрического поля, т.е. к возникновению неоднородности электронной плотности и полю разделения зарядов, что, в свою очередь, приводит к перераспределению ионной плотности и ускорению ионов в радиальном направлении до высоких энергий. Динамику этого процесса можно описать с помощью кинетических уравнений для функций распределения частиц, электронов и ионов плазмы, в которых влияние электрического поля лазера учитывается в кинетическом уравнении для электронов в виде дополнительной электростатической силы, задающей действие усредненной по высокой частоте лазера пондеромоторной силы [14, 15]. Однако как показывают результаты численного моделирования [5], до тех пор пока в плазме не возникают значительные градиенты температуры и не происходит опрокидывание радиального профиля ионного потока, электрическое поле в плазме в области распространения пучка определяется условием баланса электростатической силы, -еЕ, возвращающей электроны в канал, и выдавливающей их понде-
ромоторной силы, ^ = -те2У\11 + а2/2 (а — стандартный безразмерный вектор-потенциал лазерного излучения), т.е. еЕ = ¥р. Это вполне естественно в силу малости массы электронов, которые, быстро осциллируя на лазерной и плазменной частотах, стремятся к установлению такого равновесного состояния. При этом также считается, что электронное давление мало по сравнению со световым.
С учетом этого факта можно сильно упростить математическую модель для описания ускорения ионов, оставив лишь кинетическое уравнение для ионов, в котором электрическое поле полагается заданным. Учитывая симметрию лазерного пуч-
ка, распространяющегося вдоль оси г, мы будем рассматривать ионное кинетическое уравнение для проинтегрированной по продольной и аксиальной компонентам скоростей функции распределения ионов в цилиндрической системе координат {г, ф, г} с учетом ее зависимости только от времени t, радиальной координаты г и радиальной компоненты ионной скорости V г. В результате приходим к задаче Коши для кинетического уравнения для (проинтегрированной по V ф и ) ионной функции распределения частиц плазмы (тепловым разбросом по скоростям Vф и для
определенности пренебрегаем) /' (т, х, V)
(1)
дт/ + vд х/ + рду/ = 0, / = х/, р = р(т, х), /| т=0 = /о( х, V).
В уравнении (1) используются безразмерные переменные: безразмерное время т = , — ленгмюровская частота ионов с массой М' = АМ
и зарядом е = 1е, где А и Z — массовое и зарядовое число ионов, М — масса протона; безразмерная координата х = г/Ь, Ь — масштаб локализации лазерного пучка по радиусу; безразмерная скорость ионов V = V,. /(юЬ'Ь); безразмерное
электрическое поле р = ^еЕ/(М'ЬюЬ')); безразмерная функция распределения ионов
/' = (я°/(юЬ'Ь))/, п0 — невозмущенная плотность ионов. Выбор в (1) нормировки /' на характерную скорость ионов юЬ'Ь, а не на типично используемую тепловую скорость, связан с тем, что в обсуждаемой здесь постановке задачи тепловой разброс ионной функции распределения по скоростям несущественен и, соответственно, тепловая скорость ионов мала по сравнению с используемой для нормировки величиной юЬ'Ь.
Учитывая связь электрического поля в плазме с градиентом лазерной интенсивности, выражение для электрического поля р принимает вид
Р = -
,л2 г2 ®ЬеЬ
1 +
а (т, х)
1/ 2
(2)
Здесь юЬе — ленгмюровская частота электронов с
массой
т
и зарядом е = -е,
величина
а2(т, х) = А(т)а2(х) характеризует безразмерную интенсивность лазера, а2(х) = а^10(х, а0), а0 =
= 0.85х 10-9X [мкм]/100 [Вт/см2], где X — длина волны лазера в мкм, 100 — интенсивность лазерного импульса в максимуме в Вт/см2, а функция А(т) определяет форму лазерного импульса.
Функция а 2(х) характеризует распределение ин-
2
с
тенсивности лазера по радиусу; например, в численных расчетах в [5] использовался вариант с
а2(х) = а2ехр(-х2), в котором 10 = ехр(-х2) не
2
зависит от а0 и равна единице на оси пучка,
Iо(0, ао2) = 1.
При исследовании уравнений вида (1) стандартным является метод характеристик. Наиболее просто это уравнение исследуется при не зависящем от времени электрическом поле, р = р(х), что отвечает квазистационарному режиму, при котором длительность лазерного импульса превосходит характерное время ускорения ионов. В этом случае характеристики уравнения
(1) интегрируются и его решение / = Г(/1, /2) записывается через инварианты / и /2
Ji
v 2
ф, J2 =Т-J
d Ф
У(Ф)л/2( ji + Ф)'
(3)
a (x) = 2
f f 1 +
1 +
2 a0
cosh 2 x
(7)
Подобно использованному в численном моделировании [5] экспоненциальному закону спадания с увеличением х лазерной интенсивности,
10 = ехр(-х2), выбранная нами зависимость (7) также характеризуется монотонным убыванием
интенсивности 10(х, а^) при удалении на бесконечность и имеет максимум, локализованный в
центре пучка, 10(0, а^) = 1.
Использование (5) в (3) порождает решение кинетического уравнения (1) следующего вида:
/ = Г (х, V), т = т^р^у2 + сс^Н -2 х,
2
V2 = V-, р
где потенциал Ф и ¥(Ф) связаны неявными соотношениями
sinh x = sinh x cosh т-
v cosh x sinh т
Г~2 -2 ' Vv + cosh x
(8)
p( x) = Э xФ = ¥(Ф).
(4)
Для лучшего понимания процесса ускорения ионов наибольший интерес представляет случай, когда интегралы (3) вычисляются в явном виде. Соответственно, требуется специальное задание пространственного распределения электрического поля р, которое тем не менее должно соответст
Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.