Физика плазмы, 2021, T. 47, № 9, стр. 829-855

Особенности радиального и аксиального распределений тороидального магнитного поля в осевом струйном выбросе на установке ПФ-3

В. И. Крауз a*, К. Н. Митрофанов b**, В. В. Мялтон a, И. В. Ильичев a, А. М. Харрасов a, Ю.В. Виноградова a

a НИЦ “Курчатовский институт”
Москва, Россия

b ГНЦ РФ “Троицкий институт инновационных и термоядерных исследований”
Москва, Россия

* E-mail: krauz_vi@nrcki.ru
** E-mail: mitrofan@triniti.ru

Поступила в редакцию 26.03.2021
После доработки 27.04.2021
Принята к публикации 29.04.2021

Полный текст (PDF)

Аннотация

Представлены результаты экспериментов по исследованию особенностей динамики и внутренней структуры осевого плазменного выброса на установке ПФ-3, выполненных в рамках программы лабораторного моделирования астрофизических джетов. Проведены измерения распределения тороидального магнитного поля в радиальном Bφ(r) и осевом Bφ(z) направлениях. Показано, что магнитное поле играет существенную роль в коллимации плазменного потока, что важно для верификации численных моделей распространения астрофизических джетов. Измерения магнитных полей с пространственным разрешением вдоль оси пролетной камеры позволили определить протяженность области с захваченным магнитным потоком. Исследования радиальной структуры центральной части плазменного выброса свидетельствует об увеличении в диаметре канала протекания центрального тока в направлении от головной к хвостовой части плазменного потока. Приближенно линейный характер распределения магнитного поля в виде функциональной зависимости Bφ(r) ~ r в области протекания центрального тока указывает на равномерное распределение плотности тока вдоль радиуса. Выявлены некоторые особенности протекания токов в головной и хвостовой части потока. Обсуждается взаимодействие плазменного потока с конструкционными элементами пролетной камеры.

Ключевые слова: осевой плазменный выброс, магнитное поле, динамика плазмы

1. ВВЕДЕНИЕ

Лабораторное моделирование астрофизических процессов является одним из интенсивно развиваемых направлений. Актуальность этого направления обусловлена естественными ограничениями, которые присущи наблюдательной астрономии и связаны как с удаленностью исследуемых объектов, так и с возможностями существующих средств наблюдения. В этих условиях особое значение приобретает развитие теоретических моделей и разрабатываемых на их основе численных кодов. При этом остро встает вопрос их верификации. Здесь особая роль отводится лабораторному моделированию. Этому способствует тот замечательный факт, что уравнения магнитной гидродинамики не имеют внутреннего масштаба. Это означает, что, несмотря на огромные различия в пространственных и временных масштабах, при соблюдении определенных соотношений подобия, процессы в астрофизической плазме и в лабораторной плазме описываются одними и теми же уравнениями. Эти соотношения достаточно четко сформулированы в ряде работ, например, [13].

Одним из объектов лабораторного моделирования являются выбросы (джеты), испускаемые молодыми звездными объектами. Астрофизические джеты представляют собой одно из ярчайших явлений, наблюдаемых во Вселенной и присуще широкому классу объектов – от активных ядер галактик до микроквазаров и молодых звезд [4]. Интерес к лабораторному моделированию именно выбросов молодых звезд обусловлен тем, что они обладают нерелятивистскими скоростями (~3 × 107 см/с), достижимыми в лабораторных условиях. Хорошо известны работы по моделированию джетов молодых звезд, выполненные на мощных лазерных установках [5, 6], Z-пинчах [7, 8], плазменных ускорителях [9, 10].

В последнее время все более широко используется термин “течение Хербиг-Аро” [11, 12], включающее в себя объекты, возникающие при взаимодействии выбросов с окружающей средой. Пример такого течения показан на рис. 1.

Рис. 1.

Оптическое изображение объекта HH 111 из работы [13].

Как видно, течение может быть достаточно сильно структурировано. Причина такой структурированности является в настоящее время объектом дискуссий. Естественно также возникает вопрос: как эти отдельные объекты связаны с “центральной машиной”.

Одним из устройств, позволяющих исследовать эти вопросы, является “плазменный фокус” (ПФ). ПФ является источником интенсивных плазменных потоков, возникающих в стадии пинчевания и распространяющихся в пролетной камере на расстояния ≥100 см, значительно превышающие их поперечные размеры [1417]. Отличительной особенностью экспериментов с ПФ является то, что поток распространяется в среде с конечной плотностью [18, 19], что делает возможным моделирование взаимодействия джета с фоновой плазмой. Оценки показали, что в этой схеме моделирования удовлетворяются основные требования к параметрам масштабирования, что позволяет распространять полученные результаты на астрофизические объекты [1921].

Среди основных результатов, полученных в схеме с ПФ установками, помимо компактности плазменного выброса при его распространении в пролетной камере, необходимо выделить такой важный результат, как обнаружение центрального тока, протекающего в керне плазменного выброса [22] и являющегося следствием захваченного магнитного потока, возникающего на стадии формирования потока в результате перезамыкания силовых линий. Понятно, что в этом случае топология магнитного поля в выбросе определяется топологией поля, создаваемой источником установки, что приводит к преобладанию тороидальной компоненты. Понятно, что в случае астрофизических джетов топология поля может быть иной. Тем не менее, и в этом случае наличие центрального продольного тока, а, следовательно, и тороидального магнитного поля, независимо от природы его возникновения, является одним из ключевых аспектов существующих теоретических моделей. Это тороидальное магнитное поле может играть существенную роль в обеспечении компактности выброса.

Таким образом, плазменный поток распространяется с собственным захваченным полем. Эти результаты были получены как на установке ПФ-3 в НИЦ “Курчатовский институт” [23], так и на установках КПФ-4 “Феникс” в Сухумском физико-техническом институте [25] и PF-1000 в Институте физики плазмы и лазерного микросинтеза в Варшаве [26]. Согласно разработанной в этих работах феноменологической модели на основании детального исследования радиального распределения тороидального магнитного поля, центральный ток замыкается на периферии выброса непосредственно на сам плазменный сгусток и не связан с основной цепью разряда. Однако детальных исследований распределения тороидального поля в продольном направлении выполнено не было. К этой задаче тесно примыкает и проблема распределения токов замыкания в головной и хвостовой областях выброса.

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

2. ПОСТАНОВКА ЭКСПЕРИМЕНТА

2.1. Схема эксперимента

Исследования осевого плазменного потока, возникающего в результате сжатия плазмы ТПО и образования пинча, проведены на установке ПФ-3 [26, 27], представляющей собой плазменный фокус с конфигурацией электродов типа Филиппова (см. рис. 2а). Разрядная система состоит из двух электродов (анод диаметром 92 см и штыревой катод типа “беличье колесо” диаметром 115 см, помещенные в вакуумную камеру диаметром 250 см), разделенных изолятором высотой 25 см. Эксперименты, представленные в настоящей работе, выполнены при стационарном напуске рабочего газа в камеру установки при следующих начальных условиях: рабочий газ Ne при давлении P0 = 2.0 Торр. Энергосодержание емкостного накопителя составляло ~370 кДж при зарядном напряжении 9 кВ, разрядный ток в системе – ≤2 МА. На этом же рис. 2а схематически показаны различные стадии сжатия токово-плазменной оболочки (ТПО) и образования пинча на оси установки в прианодной области (в анодной яме). Там же показана область распространения осевого струйного выброса из анодной области в пространство пролетной камеры.

Рис. 2.

Общий вид разрядной камеры установки ПФ-3 и ее электродной системы (а), пролетной камеры и расположения некоторых диагностических средств (б), (в): 1 – вакуумная камера; 2 – пролетная камера; 3 – магнитные зонды; 4 – оптические коллиматоры; 5 – смотровые окна оптических диагностик (кадровых ЭОП и камеры щелевой развертки K008). Во врезке к рис. (а) показан вид пролетной камеры сверху. На рис. (б) и (в) показаны варианты расположения магнитных зондов для измерения Bφ(r) и Bφ(z) распределений магнитного поля, соответственно.

Исследования динамики распространения осевого плазменного потока в окружающей среде из фонового газа проводились на значительных расстояниях от места генерации – на пролетной базе вдоль оси Z установки от 30 см до 50 см. Для этой цели используется специально сконструированная пролетная камера (см. рис. 2б, в) с набором диагностических патрубков, позволяющим проводить измерения большого количества параметров плазменного потока (пространственных размеров, скорости, магнитного поля, интенсивности собственного свечения и др.) на различных расстояниях от места генерации. Центры диагностических патрубков диаметром 4 см находятся на расстоянии 35 см от плоскости анода. Момент так называемой “особенности” на производной полного тока dI/dt условно принимается в качестве момента генерации осевого потока, а местом генерации в первом приближении принимается анодная плоскость (z = 0).

2.2. Магнитозондовый метод измерения распределения тороидального магнитного поля в плазменном потоке

В экспериментах использовались абсолютно калиброванные многоканальные магнитные зонды, позволяющие измерять тороидальную (азимутальную) компоненту магнитного поля в плазме при распространении осевого струйного выброса в дрейфовом пространстве пролетной камеры. Характеристики магнитных зондов для измерения тороидальной (Bφ) компоненты магнитного поля были следующие: пространственное разрешение ~4 мм; временное разрешение – не хуже 1.5 нс; чувствительность катушек зонда kp = (5.30 ± 10.01) × 10–8 В/(Гс/с), количество каналов измерения – 9, дискретность измерения пространственного распределения Bφ-поля – Δr или Δz = 8 мм. Общий вид и подробное описание конструкции магнитного зонда представлены на рис. 3 и в работах [23, 24]. Чувствительный элемент зонда калибровался как на величину вектора индукции магнитного поля, так и на его направление. Точность калибровки магнитных зондов (определения эффективной площади катушки) была не хуже 5%. Погрешность измерений магнитного поля в плазменном потоке с учетом обтекания оболочки зонда плазмой была не более 20% [2829].

Рис. 3.

Многоканальный магнитный зонд для измерения тороидальной (Bφ-) компоненты магнитного поля в плазменном потоке. Во врезке к рисунку показана конструкция чувствительного элемента магнитного зонда с указанием размеров катушек и расстояний между ними.

Набор диагностических патрубков, позволяет располагать магнитные зонды в диаметрально противоположных направлениях для измерения Bφ(r) распределения магнитного поля по обе стороны от оси (см. рис. 2б) в среднем сечении первой секции пролетной камеры (центр диагностических окон z1 = 35 см) и в осевом направлении с торца пролетной камеры (см. рис. 2в) для измерения Bφ(z) распределения на расстояниях 30–50 см от поверхности анода.

Использование двух таких зондов позволяло проводить одновременные измерения в 18-ти точках пространства пролетной камеры в зависимости от цели эксперимента. В случае измерения радиального распределения Bφ-поля измерения проводились в диапазоне радиусов r ∈ (–10, 10) см, в случае измерения Bφ(z) – в диапазоне z ∈ ∈ (30, 50) см от плоскости анода. Сигнал с i-й катушки зонда, расположенной в точке пространства (rizi, φi), был пропорционален производной индукции магнитного поля ($U_{p}^{i}\sim \partial B_{\varphi }^{i}{\text{/}}\partial t$). Для получения временной зависимости индукции магнитного поля $B_{\varphi }^{i}(t)$ в данной точке сигнал с катушки $U_{p}^{i}$ численно интегрировался с учетом полученных калибровочных коэффициентов $k_{p}^{i}$. Пространственные распределения магнитного поля Bφ(r) или Bφ(z) восстанавливались по семейству кривых $B_{\varphi }^{i}(t)$ в зависимости от расположения магнитных зондов в пролетной камере. Программным способом из зондовых сигналов восстанавливались радиальное и осевое распределения магнитного поля непрерывно во времени, и на основе полученных данных создавался видеоряд в виде avi-файла, состоящего из последовательности захваченных изображений размером 1120 × 700 пикселей, меняющихся со скоростью 4 кадра в секунду. Полученный видеоряд позволял оперативно проводить анализ полученных зондовых результатов. Таким образом, при помощи циклической обработки численных данных (по ~80 выстрелам) создана база зондовых измерений, анализ которой, например, позволил в дальнейшем определить положение протекания центрального тока и токов замыкания, протяженность области с магнитным полем и др.

2.3. Методики регистрации оптического излучения, кадровых изображений и щелевых разверток плазменного потока

Фотографирование плазмы в оптическом диапазоне проводилось с помощью электронно-оптических преобразователей (ЭОП) типа ЭП-16 [30] в (rz)-плоскости пролетной камеры. Каждый оптический регистратор представляет собой сборку из ЭОП, оптических фильтров, дистанционно управляемого цифрового фотоаппарата и системы объективов для переноса изображения. Регистрация плазменных потоков производилась в центральном сечении (z = 35 см) пролетной камеры. Экспозиция кадров составляла 30 нс. Поле зрения ЭОП-диагностики в приосевой области пролетной камеры захватывало область диаметром 13 см, полная разрешающая способность системы, приведенная к объекту, была не хуже 1 мм.

В качестве одного из диагностических инструментов в данной работе использована также скоростная камера с щелевой разверткой К008. Оптическая схема (см. рис. 4), включала в себя упорядоченный световод, представляющий собой систему из нескольких плоских жил, на одном конце объединенных в квадратную сборку и с разветвленным противоположным концом [31, 32]. Перед световодом располагалась оптическая система, представляющая собой телескоп Кеплера. Световод обращен соединенным концом к пролетной камере установки и расположен таким образом, чтобы регистрировать излучение с высоты z = 30.5 см, 35.0 см и 39.5 см над анодом. Несоединенные жилы на противоположном конце световода были собраны в линейку для выведения оптического излучения на камеру щелевой развертки. Для синхронизации был использован сигнал с оптического коллиматора, регистрирующего момент прохода ТПО определенного положения в камере. Ширина области зрения в радиальном направлении, перпендикулярном оси Z камеры составляла 9 см. Пространственное разрешение по объекту было не хуже 1 мм, временное разрешение – не хуже 400 нс для развертки 20 мкс и 120 нс для развертки длиной 6 мкс.

Рис. 4.

Схема регистрации потока в пролетной камере посредством щелевой камеры К-008. 1 – пролетная камера, 2 – оптическая ось системы, 3 – плазменный поток, 4 – собирающие линзы (телескоп Кеплера), 5 – фокальная плоскость, f1, f2 – фокусные расстояния линз, 6 – изображение плазменного потока, 7 – упорядоченный световод, 8 – входная щель камеры К-008, 9 – камера К-008, 10 – развертка потока.

В экспериментах также использовались оптические коллиматоры, собирающие световое излучение вдоль диаметра камеры, которое затем по световоду подается на фотокатод ФЭУ и регистрируется на осциллографе. Коллимирование обеспечивает регистрацию излучения из относительно небольшого плазменного объема: диаметр области, попадающей в поле зрения коллиматора, не превышает 5 мм в центре пролетной камеры. Это позволяет с достаточно высокой точностью регистрировать момент прихода плазменного потока в область наблюдения и определять среднюю на пролетном участке скорость. Использование двойного коллиматора, состоящего из двух таких каналов, как показано на рис. 2б, в, и разнесенных на небольшое расстояние (1.6 см), позволяет фактически определять мгновенную скорость в области наблюдения по сдвижке появления сигнала в различных каналах коллиматора.

3. ЭКСПЕРИМЕНТАЛЬНЫЕ РЕЗУЛЬТАТЫ

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

Радиальное распределение тороидального магнитного поля было объектом наших предыдущих исследований [23], основным результатом которых явилось подтверждение наличия центрального керна в осевом плазменном выбросе, в котором протекает осевой ток, а также обнаружение области протекания обратных токов замыкания. Увеличение количества измерительных катушек в магнитных зондах позволило уточнить динамику радиального распределения тороидального магнитного поля. Сопоставление данных оптического кадрового фотографирования (см. рис. 5) и магнитных измерений позволило обнаружить корреляцию между формой плазменного потока и временной зависимостью магнитного поля, регистрируемой зондами.

Рис. 5.

Кадровые оптические изображения плазменного потока (импульсы № 4581 и № 4587, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж) по данным ЭОП-регистрации при прохождении плазмой позиции z = 35 см в момент времени topt относительно момента особенности на производной полного тока: а) – на topt = 4.1 мкс; б) – на topt = 4.3 мкс. Масштаб сетки 1 см × 1 см.

Измерения были проведены вдоль центрального диаметра первой секции пролетной камеры (z = 35 см). Радиальные распределения тороидального магнитного поля Bφ(r, 35) восстановлены из семейства кривых $B_{\varphi }^{*}$(ri, 35, t) и $B_{\varphi }^{{**}}$(rj, 35, t) на некоторые моменты времени ti, указанные стрелками (↓) на рис. 6 и 8. Для примера на рис. 6 справа снизу в легенде показано соответствие представленных кривых радиусам расположения катушек зондов. На рис. 7 и 9 показаны соответствующие радиальные распределения магнитного поля. Так как два многоканальных зонда были расположены по обе стороны от оси Z в диаметрально противоположных направлениях, это позволяло исследовать как симметрию расположения области протекания центрального тока относительно оси пролетной камеры, так и радиальный размер этой области, а также область протекания обратных токов на периферии плазменного потока.

Рис. 6.

Результаты измерений тороидальных магнитных полей в плазменном потоке на расстоянии z = 35 см от поверхности анода установки (импульс № 4553, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): временные зависимости индукции тороидального магнитного поля, измеренной магнитными зондами на радиусах ri ∈ [–2.5; –1.7; –0.9; –0.1; 0.7; 1.5; 2.4; 3.3; 4.2] см [кривые $B_{\varphi }^{*}$(ri, 35, t)] и rj ∈ [–3.5; –4.2; –4.9; –5.8; –6.6; –7.4; –8.2; –9.1; –10.0] см [кривые $B_{\varphi }^{{**}}$(rj, 35, t)]. Справа вверху схематично показано распределение магнитного поля в центральной части плазменного потока в различных его сечениях. Справа снизу в легенде показано соответствие представленных кривых радиусам расположения катушек зондов.

Рис. 7.

Результаты измерений тороидальных магнитных полей в плазменном потоке на расстоянии z = 35 см от поверхности анода установки (импульс № 4553, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): (а–е) – радиальные распределения магнитного поля Bφ(r), восстановленные из кривых $B_{\varphi }^{*}$(ri, 35, t) и $B_{\varphi }^{{**}}$(rj, 35, t) рис. 6 на моменты времени t1t6, указанные стрелками (↓).

Рис. 8.

Результаты измерений тороидальных магнитных полей в плазменном потоке на расстоянии z = 35 см от поверхности анода установки (импульс № 4562, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): временные зависимости индукции тороидального магнитного поля, измеренной магнитными зондами на радиусах ri ∈ [3.5; 4.3; 5.1; 5.9; 6.7; 7.5; 8.4; 9.3; 10.2] см [кривые $B_{\varphi }^{*}$(ri, 35, t)] и rj ∈ [–3.5; –4.2; –4.9; –5.8; –6.6; –7.4; –8.2; –9.1; –10.0] см [кривые $B_{\varphi }^{{**}}$(rj, 35, t)]. Справа в легенде показано соответствие представленных кривых радиусам расположения катушек зондов.

Рис. 9.

Результаты измерений тороидальных магнитных полей в плазменном потоке на расстоянии z = 35 см от поверхности анода установки (импульс № 4562, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): (а–е) – радиальные распределения магнитного поля Bφ(r), восстановленные из кривых $B_{\varphi }^{*}$(ri, 35, t) и $B_{\varphi }^{{**}}$(rj, 35, t) рис. 8 на моменты времени t1t6, указанные стрелками (↓).

В рассматриваемом на рис. 6, 7 случае такое размещение зондов позволяет проводить измерения в области радиусов от r = 4.2 см справа до r = –10 см слева от оси камеры, т.е. с достаточно высокой степенью точности промерить распределение магнитного поля в области оси камеры. В некоторых экспериментах наблюдалась несимметричная форма радиального распределения тороидального магнитного поля Bφ(r), измеренного двумя магнитными зондами при прохождении переднего фронта плазменного потока (см. на рис. 7а распределение Bφ(r) в момент времени t1). Как следует из распределения Bφ(r) наблюдается несимметричный пролет плазменного потока относительно магнитного зонда, расположенного слева от оси пролетной камеры. Скорее всего, практически все катушки данного зонда находятся вне области локализации магнитного поля осевого струйного выброса. В следующий момент времени (см. распределение Bφ(r) в момент t2) уже не менее 6‑ти катушек зонда оказывается внутри плазменного потока и регистрируют магнитное поле. Общий вид распределений Bφ(r) на рис. 7 и 9 указывает на существование двух областей в плазменном потоке: области центрального тока, где Bφ(r) ~ r и области, где Bφ(r) ~ 1/rn, n ≈ 0.9–1.0. Характер изменения со временем радиального размера области с Bφ(r) ~ r свидетельствует об увеличении диаметра канала протекания центрального тока в направлении от головной к хвостовой части плазменного потока, как схематически представлено справа вверху на рис. 6 или рис. 19е. Полученный результат находится в хорошем согласии с результатами, полученными ранее в экспериментах на установке КПФ-4 (г. Сухум) [24], а также с моделью, разработанной в работе [33]. В период времени t2t6 диаметр области центрального тока увеличивается с 4.4 см до величины более 7 см. Приближенно линейный характер распределения магнитного поля в виде функциональной зависимости Bφ(r) ~ r в области протекания центрального тока указывает на равномерное распределение плотности тока jz вдоль радиуса в данной области, согласно следующему выражению jz(r) ∝ ∂Bφ/∂r + Bφ/r. Также возрастает и радиус области протекания обратного тока. Увеличение внешнего размера плазменного потока в направлении от его головы к хвостовой части наблюдается и на кадровых оптических регистраторах (см. рис. 5).

Аналогичная картина наблюдается в другом выстреле (№ 4562), когда катушки магнитных зондов были разнесены по радиусу от оси пролетной камеры (см. рис. 8 и 9). В данном случае измерения непосредственно на оси не производились, однако такое расположение катушек позволяло исследовать поперечные размеры плазменного потока вплоть до радиуса стенок пролетной камеры (r ≈ 10.3 см). В период времени t1t5 диаметр области протекания центрального тока увеличивается с ~7 см вплоть до 10 см (см. рис. 9). Тороидальное магнитное поле на периферии плазменного потока регистрируется вплоть до радиуса r ~ 8 см на моменты времени t2, t3, а в последующие моменты времени t4t6 – вплоть до радиуса стенки пролетной камеры. Таким образом, в хвостовой части плазменного потока область, занимаемая захваченным Bφ-полем занимает весь объем, ограниченный стенками пролетной камеры. Вероятней всего, обратные токи, протекающие на периферии плазменного потока, замыкаются на металлические стенки камеры.

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

Если сохранение компактности головной части плазменного потока в радиальном направлении при его распространении вдоль оси камеры показано ранее достаточно убедительно (с помощью как магнитных, так и оптических измерений) [15, 17, 23], то вопрос о компактности осевого плазменного выброса в продольном направлении и, особенно, о месте замыкания обратных токов оставался дискуссионным. В предлагаемой нами модели на основании предыдущих измерений распределения тороидальных магнитных полей предполагается, что обратные токи замыкаются непосредственно на сам плазменный сгусток и не связаны с основной разрядной цепью [2325]. Однако прямых измерений продольного распределения тороидальных магнитных полей, и, соответственно, протяженности продольного осевого тока, проведено не было.

В данной серии экспериментов такие измерения были проведены с помощью двух магнитных зондов, вводимых с торца пролетной камеры. Зонды были расположены симметрично относительно оси камеры на радиусах 3 см и 5 см, как показано на рис. 2в. Каждый зонд мог перемещаться вертикально на различные высоты, что, в отличие от случая радиального расположения зондов, позволило исследовать динамику движения сгустка с магнитным полем вдоль оси Z в диапазоне z ∈ (30, 50) см и более точно определить протяженность области с захваченным магнитным потоком.

На рис. 10 и 12 представлены временные зависимости индукции тороидального магнитного поля [семейство кривых Bφ(3, zi, t) и Bφ(5, zj, t)], измеренные магнитными зондами на радиусах rp = 3 см и rp = 5 см (см. рис. 2в), соответственно. На рис. 11 и 13 показаны аксиальные распределения магнитного поля Bφ(z), восстановленные из кривых Bφ(3, zi, t) и Bφ(5, zj, t) на некоторые моменты времени ti, указанные стрелками (↓).

Рис. 10.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 42) см от поверхности анода установки (импульс № 4523, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): временные зависимости индукции тороидального магнитного поля, измеренной магнитными зондами на радиусах rp = 3 см [кривые Bφ(3, zi, t)] и rp = 5 см [кривые Bφ(5, zj, t)], где zi ∈ [30.4; 31.2; 32.0; 32.8; 33.6; 34.4; 35.3; 36.2; 37.1] см и zj ∈ [34.9; 35.6; 36.3; 37.2; 38.0; 38.8; 39.6; 40.5; 41.4] см. В легенде показано соответствие представленных кривых позициям расположения катушек зондов вдоль оси Z.

Рис. 11.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 42) см от поверхности анода установки (импульс № 4523, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): (а–л) – аксиальные распределения магнитного поля Bφ(z) на радиусах rp = 3 см и 5 см, восстановленные из кривых Bφ(3, zi, t) и Bφ(5, zj, t) рис. 10 на моменты времени t1t11, указанные стрелками (↓).

Рис. 12.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 46) см от поверхности анода установки (импульс № 4533, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): временные зависимости индукции тороидального магнитного поля, измеренной магнитными зондами на радиусах rp = 3 см [кривые Bφ(3, zi, t)] и rp = 5 см [кривые Bφ(5, zj, t)], где zi ∈ [30.4; 31.2; 32.0; 32.8; 33.6; 34.4; 35.3; 36.2; 37.1] см и zj ∈ [39.9; 40.6; 41.3; 42.2; 43.0; 43.8; 44.6; 45.5; 46.4] см. Справа – оптическое кадровое изображение плазменного потока, зарегистрированное в момент времени topt = 4.3 мкс. В легенде показано соответствие представленных кривых позициям расположения катушек зондов вдоль оси Z.

Рис. 13.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 46) см от поверхности анода установки (импульс № 4533, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): (а–и) – аксиальные распределения магнитного поля Bφ(z) на радиусах rp = 3 см и 5 см, восстановленные из кривых Bφ(3,zi,t) и Bφ(5,zj,t) рис. 12 на моменты времени t1t9, указанные стрелками (↓).

Рассмотрим результаты выстрела (№ 4523), в котором измерения магнитного поля проводились при прохождении плазменного потока на пролетной базе z ∈ (30, 42) см от поверхности анода установки (см. рис. 10 и 11). Особенностью этих измерений является то, что при таком расположении зондов они частично перекрывают друг друга по высоте: часть катушек находится на одинаковом удалении от анода, но на разном расстоянии от оси. Наблюдаемое на рис. 10 различие в профилях сигналов с магнитных зондов, расположенных на различных радиусах rp = 3 см и 5 см, может быть обусловлено не только более удаленным расположением по высоте катушек зонда на rp = 5 см, но и неплоской формой переднего фронта плазменного потока. На аксиальных распределениях магнитного поля Bφ(3, zi, t) и Bφ(5, zj, t) в моменты времени t1t3, видно (см. рис. 11), что вначале Bφ-поле начинает регистрироваться катушками зонда на радиусе rp = 3 см. При этом распределение магнитного поля Bφ(3, zi, t) в момент времени t3 в головной части плазменного потока плавно спадает от своего пикового значения $B_{\varphi }^{{{\text{max}}}}$ ≈ 1.2 кГс, т.н. локального максимума магнитного поля (см. т. А – обозначена стрелкой “←” на рис. 11) до нуля в диапазоне расстояний z ∈ (30, 34) см. В моменты времени t3t8 наблюдается перемещение области $B_{\varphi }^{{{\text{max}}}}$(3, zi, t) вдоль оси Z со скоростью $V_{z}^{A}$ ≈ (2.2–3.6) × 106 см/с. При этом следует заметить, что величина Bφ-поля в точке А уменьшается в несколько раз по мере пролета плазмы от z = 30 см до z = 36 см, что свидетельствует о затухании собственных токов, циркулирующих в плазме. На аксиальном распределении Bφ(5, zj, t), измеренном катушками зонда на радиусе rp = 5 см, наблюдается та же картина (см. перемещение т. B на рис. 11 в моменты t10, t11).

Характер продольного распределения тороидального поля (наличие локального максимума на распределении Bφ(z)) свидетельствует, скорее всего, о сгустке достаточно малых продольных размеров (несколько сантиметров) в области протекания центрального тока с последующим расширением токового канала, что соответствует результатам измерения радиального распределения тороидального поля, приведенным в п. 3.1 настоящей статьи, и зарегистрированным компактным объектам в экспериментах на установке PF-1000 [17, 25]. После прохождения основного сгустка в большинстве случаев магнитное поле стремится к нулю, что свидетельствует о замыкании основного осевого тока непосредственно на сгусток. Наличие небольшого, отличного в ряде случаев от нуля магнитного поля, может быть следствием утечки части тока по остаточной плазме непосредственно на электроды системы, либо обусловлено прохождением генерируемых в разряде запаздывающих потоков с собственным магнитным полем.

Из аксиальных распределений тороидального магнитного поля определена протяженность области (Δz) с захваченным магнитным потоком вдоль оси Z пролетной камеры установки. Так из распределения магнитного поля Bφ(3, zi, t) рис. 11 в момент времени t6 видно, что величина Δz составляет величину около 5.7 см, в то время как в других выстрелах величина Δz могла быть несколько больше – 6–7 см. В последнем случае большое значение Δz обусловлено в основном длительной хвостовой частью плазменного потока. В п. 3.3 настоящей статьи будут более подробно рассмотрены особенности распределения магнитного поля в головной и хвостовой частях плазменного потока.

Далее, рассмотрим поведение осевого струйного выброса на бóльшей пролетной базе вдоль оси пролетной камеры. В выстреле № 4533 измерения тороидального магнитного поля проводились при прохождении плазменного потока на пролетной базе z ∈ (30, 46) см от поверхности анода установки (см. рис. 12 и 13). В данном случае не было частичного перекрытия позиций расположения катушек двух многоканальных магнитных зондов вдоль оси Z, как в рассмотренном ранее случае на рис. 11. Как и следовало ожидать, увеличилась задержка появления сигналов с зондов на rp = 3 см и 5 см с 0.7 мкс до 1.5 мкс. На оптическом кадровом изображении (см. рис. 12 справа), зарегистрированном на topt = 4.3 мкс вблизи момента появления сигналов с зонда на rp = 3 см, видно, что указанный зонд уже погружен в плазму. В это же время в области радиуса rp = 5 см не наблюдается оптического свечения плазмы, т.е. передний фронт плазменного потока туда не дошел. В рассматриваемом выстреле скорость смещения $V_{z}^{A}$ области с максимальным уровнем магнитного поля на аксиальном распределении Bφ(3, zi, t) лежит в диапазоне (2.2–3) × 106 см/с (см. распределения магнитного поля в моменты t1t5 на рис. 13). Разница в увеличенной на 4 см пролетной базе по сравнению со случаем, представленным на рис. 11, объясняет увеличенную до 1.5 мкс задержку в появлении сигналов с зондов на радиусах rp = 3 см и 5 см с учетом полученной оценки скорости $V_{z}^{A}$. В представленном примере перемещение области с магнитным полем прослеживается более отчетливо. Протяженность области с захваченным магнитным потоком составила Δz ~ ~ 7 см.

Следует заметить, что на представленных аксиальных распределениях магнитного поля Bφ(z) отсутствует область с Bφ(z) = Const (т.н. “плато”), связанная с прохождением протяженной вдоль оси Z области протекания центрального тока. Действительно, все распределения Bφ(3, zi, t) на рис. 11 и 13 имеют в голове потока плавно нарастающий до максимума профиль и сразу же спадающий после максимума профиль в хвосте потока. Однако, в некоторых выстрелах зарегистрирована иная картина.

На рис. 14 представлена радиальная оптическая щелевая развертка на различных расстояниях zopt = [30.5; 35.0; 39.5] см вдоль оси пролетной камеры. На всех подобных изображениях щелевых разверток в данной статье штриховыми линиями обозначены радиусы rp расположения магнитных зондов, а штрихпунктирными линиями – ось Z пролетной камеры. В правой части каждого изображения развертки построена шкала соответствия между отсчетами интенсивности и цветовыми градациями. Число отсчетов интенсивности для 8-битного представления изображения щелевой развертки лежит соответственно в диапазоне от 0 до 255.

Рис. 14.

Результаты сравнения данных с оптической радиальной развертки (вверху) и магнитных зондов (внизу) в импульсе № 4540, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж: 13 – временные зависимости индукции тороидального магнитного поля, измеренной магнитными зондами на радиусах rp = 3 см на высотах zp = [30.4; 35.3; 39.9] см, соответственно. Пространственные щели оптической развертки располагались на высотах zopt = [30.5; 35.0; 39.5] см.

На представленной щелевой развертке видно длительное (до 2 мкс) прохождение области плазменного потока с интенсивным свечением, обозначенное на рис. 14 фигурной скобкой с подписью “плато”. Средняя скорость движения данной области на пролетной базе zopt = 30.5–39.5 см, определенная по задержке появления свечения на различных высотах, составила $V_{z}^{{opt}}$ ≈ 3 × 106 см/с. В этом же промежутке времени t1t6 на аксиальных распределениях тороидального магнитного поля Bφ(3, zi, t) (см. рис. 15) наблюдается прохождение области Bφ(z) = Const. Такое “плато” на распределении Bφ(z) образуется, когда некоторые катушки зонда, в данном примере 5 катушек, расположенные на одном и том же радиусе rp = 3 см, находятся снаружи области протекания центрального тока, где Bφ ~ 1/r. В этом случае катушки будут регистрировать одинаковый уровень тороидального магнитного поля. Как следует из распределения Bφ(z), представленного в момент времени t3 на рис. 15в, область плато имеет протяженность $\Delta z_{{pl}}^{{mag}}$ ~ 3.3 см вдоль оси Z пролетной камеры. Полученное значение $\Delta z_{{pl}}^{{mag}}$ хорошо совпадает с оценкой наиболее ярко светящейся части плазменного потока $\Delta z_{{pl}}^{{opt}}$ ~3.9 см, полученной из оптической развертки следующим образом: $\Delta z_{{pl}}^{{opt}} = V_{z}^{{opt}} \cdot \Delta t$, где Δt ≈ 1.3 мкс. Следует также заметить, что небольшое “плато” наблюдается и на распределениях Bφ(5, zi, t), зарегистрированных зондом на радиусе rp = 5 см (см. рис. 15д, е), т.е. область с распределением тороидального магнитного поля в виде Bφ ~ 1/r может простираться не менее чем до радиуса 5 см.

Рис. 15.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 46) см от поверхности анода установки (импульс № 4540, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): (а–е) – аксиальные распределения магнитного поля Bφ(z) на радиусах rp = 3 см и 5 см, восстановленные на моменты времени t1t9, указанные стрелками (↓) на рис. 14.

Динамика распространения осевого струйного выброса вдоль оси Z пролетной камеры установки анализировалась по результатам обработки данных с оптических коллиматоров (о) и магнитных зондов () и представлена на рис. 16. На верхней части рисунка по оси абсцисс отложена задержка появления сигнала с оптического коллиматора или магнитного зонда на различной высоте z от поверхности анода относительно “особенности” на производной полного тока dI/dt, по оси ординат отложена позиция коллиматора (или зонда) на высоте z. Иными словами, представлена временная зависимость z(t) появления свечения плазмы (или магнитного поля) на расстоянии z от анода установки. Там же кривыми 1 и 2 представлена интерполяция экспериментальных точек (о и ) полиномом n-степени (n = 5), соответственно.

Рис. 16.

Динамика распространения осевого струйного выброса вдоль оси Z пролетной камеры установки по результатам обработки данных с оптических коллиматоров () и магнитных зондов (). Временные зависимости: сверхуz(t) появления свечения плазмы (или магнитного поля) на расстоянии z от анода установки. Кривые 1, 2 – интерполяция экспериментальных точек ( и ) полиномом n-степени (n = 5), соответственно; снизу – кривые 3, 4 аксиальной скорости движения плазмы, полученной путем дифференцирования кривых 1 и 2, соответственно.

Интересным с точки зрения анализа кривых 1 и 2 является значение z0 = z(0) в момент особенности на dI/dt (t = 0). Видно, что кривые 1 и 2 пересекают ось ординат в диапазоне значений z0 = = 7–11 см. Полученное таким образом значение z0 в первом приближении можно считать оценкой сверху высоты пинча в момент генерации осевого струйного выброса вблизи “особенности” на производной полного тока dI/dt. Информация о значении z0 необходима для надежной оценки скорости осевого струйного выброса в момент его генерации и движения в начале пролетной базы.

На нижней части рис. 16 представлены временные зависимости скорости движения плазменного потока Vz = dZ/dt (кривые 3 и 4), рассчитанные путем численного дифференцирования кривых 1 и 2. Следует заметить, что кривые 3 и 4, полученные по данным оптических коллиматоров и магнитных зондов, оказались близки и подобны друг другу. Характер полученных кривых 3 и 4 указывает на торможение плазменного потока по мере его движения на пролетной базе z ∈ ∈ (10, 65) см. Из представленных данных рис. 16 следует, что на стадии пролета плазменным потоком центральной области первой секции пролетной камеры (z = 35 см) скорость его движения вдоль оси установки лежит в диапазоне Vz = (3–5) × 106 см/с. Близкое значение скорости Vz получено из аксиальных распределений магнитного поля Bφ(z), измеренных магнитными зондами при пролете плазмой центрального сечения первой секции камеры (см. рис. 11 и 13), а также оптических разверток (см. напр. рис. 14).

По совокупности многих серий экспериментов на установке ПФ-3 построена зависимость максимального уровня азимутального магнитного поля $B_{\varphi }^{{{\text{max}}}}$, зарегистрированного зондом в позиции rp = 2 см, от расстояния z над поверхностью анода (z = 0). На рис. 17а представлена указанная зависимость $B_{\varphi }^{{{\text{max}}}}(z)$.

Рис. 17.

Результаты измерений магнитных полей и оптического свечения плазмы осевого струйного выброса на различной дистанции z вдоль оси установки, отсчитанной от поверхности анода установки (z = 0): а) – диссипация магнитного потока Bφ-поля в плазменном потоке. Максимальная величина индукции тороидального магнитного поля $B_{\varphi }^{*}$, регистрируемая магнитным зондом в позиции rp = 2 см; б) – мгновенная скорость плазменного потока Vz, восстановленная по данным оптических коллиматоров. Аппроксимация экспериментальных точек (⚫ и ⚪) выполнена в виде ƒ(z) = ƒ0 ⋅ exp(–Kz).

Так, в случае, когда измерения проводились на расстоянии zz0 = 10 см от поверхности анода максимальный уровень тороидального магнитного поля в месте расположения зонда достигает значения $B_{\varphi }^{{{\text{max}}}}\sim 30$ кГс. В предположении, что осевой (центральный) ток в плазменном потоке совпадает с осью Z камеры установки, то пересчет Bφ-поля в точке измерения (rp = 2 см) на величину тока в предположении функциональной зависимости Bφ ~ 1/r дает оценку I = $5B_{\varphi }^{{{\text{max}}}}{{r}_{p}}$ = $10B_{\varphi }^{{{\text{max}}}}$ ≈ ≈ 300 кА. Эта оценка в 35 раз больше величины тока (~8.5 кА), измеренного зондом в первой секции пролетной камеры на высоте z = 35 см от поверхности анода в наших предыдущих экспериментах [22]. Таким образом, происходит диссипация магнитного потока Bφ-поля в осевом плазменном выбросе по мере его продвижения вдоль оси установки. Согласно следующей аппроксимации экспериментальных точек (⚫) рис. 17а в виде $B_{\varphi }^{{{\text{max}}}}(z)$ = Bφ0 ⋅ exp(–KBz) сделана оценка дистанции zmag пролетной базы, на которой происходит уменьшение величины $B_{\varphi }^{{{\text{max}}}}$ в e‑раз.

Подобным образом на рис. 17б представлена мгновенная скорость плазменного потока Vz, восстановленная по результатам обработки сигналов с оптических коллиматоров на различной высоте z от поверхности анода установки. Следует заметить, что полученное значение zmag ≈ 3.3 см оказалось на порядок меньшим, чем дистанция zopt ≈ ≈ 30 см, на которой происходит уменьшение скорости движения плазменного потока Vz в e-раз, определенной по данным с оптических коллиматоров [34]. При этом на расстоянии zz0 = 10 см от поверхности анода значение аксиальной скорости Vz составляет величину около 9 × 106 см/с. Вполне возможно, что именно с такой скоростью происходит отрыв плазменного потока на стадии его генерации – стадии, когда происходит разрушение пинча вблизи особенности на производной полного тока dI/dt. Следует заметить, что несколько иная оценка скорости Vz ≈ 7 × 106 см/с получена по данным рис. 16. Скорее всего, разница в значениях скоростей объясняется как различным способом интерполирования данных – экспоненциальной функцией ƒ(z) = ƒ0 ⋅ exp(–Kz) и полиномом 5-й степени, так и отсутствием надежных данных с оптических коллиматоров в области z ∈ (10, 30) см (см. рис. 17б).

3.3. Особенности распределения магнитного поля в головной и хвостовой частях плазменного потока

В этой части статьи представлены особенности распределения магнитного поля в головной и хвостовой частях плазменного потока, которые были зарегистрированы в некоторых экспериментах (см. рис. 18–23). Следует выделить наиболее ясно прослеживающиеся детали на профиле аксиальных распределений Bφ(z) рис. 19, 21 и 23: обратное направление тороидального магнитного поля в голове и в хвосте потока, обозначенное на рисунках в виде Bφ“–”, а также отдельные пики (локальные максимумы) на профиле Bφ(z), обозначенные буквами “A1” и “A2”.

Рис. 18.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 42) см от поверхности анода установки (импульс № 4519, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): временные зависимости индукции тороидального магнитного поля, измеренной магнитными зондами на радиусах rp = 3 см [кривые Bφ(3, zi, t)] и rp = 5 см [кривые Bφ(5, zj, t)], где zi ∈ [30.4; 31.2; 32.0; 32.8; 33.6; 34.4; 35.3; 36.2; 37.1] см и zj ∈ [34.9; 35.6; 36.3; 37.2; 38.0; 38.8; 39.6; 40.5; 41.4] см. В легенде показано соответствие представленных кривых позициям расположения катушек зондов вдоль оси Z.

Рис. 19.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 42) см от поверхности анода установки (импульс № 4519, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): (а–з) – аксиальные распределения магнитного поля Bφ(z) на радиусах rp = 3 см и 5 см, восстановленные из кривых Bφ(3, zi, t) и Bφ(5, zj, t) рис. 18 на моменты времени t1t8, указанные стрелками (↓); е) – структура осевого плазменного потока. Показаны: стрелками – схема циркуляции центрального (Ic) и обратного токов (Ir); пунктирными линиями – радиальное распределение азимутального магнитного поля в плазменном потоке Bφ(r) в его центральной части (Bφ ~ r) и на периферии (Bφ ~ 1/rn) в сечениях 1 и 2.

Рис. 20.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 42) см от поверхности анода установки (импульс № 4529, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): временные зависимости индукции тороидального магнитного поля, измеренной магнитными зондами на радиусах rp = 3 см [кривые Bφ(3, zi, t)] и rp = 5 см [кривые Bφ(5, zj, t)], где zi ∈ [30.4; 31.2; 32.0; 32.8; 33.6; 34.4; 35.3; 36.2; 37.1] см и zj ∈ [34.9; 35.6; 36.3; 37.2; 38.0; 38.8; 39.6; 40.5; 41.4] см. В легенде показано соответствие представленных кривых позициям расположения катушек зондов вдоль оси Z.

Рис. 21.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 42) см от поверхности анода установки (импульс № 4529, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): (а–и) – аксиальные распределения магнитного поля Bφ(z) на радиусах rp = 3 см и 5 см, восстановленные из кривых Bφ(3, zi, t) и Bφ(5, zj, t) рис. 20 на моменты времени t1t9, указанные стрелками (↓).

Рис. 22.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 42) см от поверхности анода установки (импульс № 4522, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): временные зависимости индукции тороидального магнитного поля, измеренной магнитными зондами на радиусах rp = 3 см [кривые Bφ(3,zi,t)] и rp = 5 см [кривые Bφ(5,zj,t)], где zi ∈ [30.4; 31.2; 32.0; 32.8; 33.6; 34.4; 35.3; 36.2; 37.1] см и zj ∈ [34.9; 35.6; 36.3; 37.2; 38.0; 38.8; 39.6; 40.5; 41.4] см. В легенде показано соответствие представленных кривых позициям расположения катушек зондов вдоль оси Z.

Рис. 23.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 42) см от поверхности анода установки (импульс № 4522, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): (а–и) – аксиальные распределения магнитного поля Bφ(z) на радиусах rp = 3 см и 5 см, восстановленные из кривых Bφ(3, zi, t) и Bφ(5, zj, t) рис. 22 на моменты времени t1t9, указанные стрелками (↓).

В некоторых выстрелах на переднем фронте временных зависимостей индукции тороидального магнитного поля (см. напр. семейство кривых Bφ(3, zi, t) на рис. 18) зарегистрирован знакопеременный характер изменения Bφ-поля. В данном случае магнитный зонд регистрирует прохождение переднего фронта плазменного потока, внутри которого существуют области с противоположным направлением силовых линий Bφ‑поля. В моменты времени t1t4 на аксиальных распределениях Bφ-поля наблюдается прохождение области с Bφ “–” (см. рис. 19) вдоль оси в диапазоне z ∈ (32, 36) см. Такая ситуация связана с азимутальной неоднородностью растекания тока в головной части плазменного потока, как схематически показано на рис. 19е. Аналогичная ситуация может реализовываться и в хвостовой части потока (см. напр. семейство кривых Bφ(3, zi, t) на рис. 20 и распределения Bφ(z) на рис. 21 в моменты t5t8). Также в некоторых случаях, в приосевой области плазменного потока (см. распределения Bφ(z) на рис. 19 в моменты t4, t5 или на рис. 23 в моменты t3t5) наблюдается прохождение нескольких областей с повышенным уровнем Bφ-поля (см. т. А1 и А2). Это может быть связано как с прохождением двух близко идущих друг за другом плазменных сгустков (см. т. А1 и А2 на рис. 19д), так и с прохождением основного плазменного потока (см. т. А2 на рис. 23в–д) с впередиидущим фор-сгустком (т. А1). Так, в первом случае расстояние между точками А1 и А2 составляет z(A1)–z(A2) ≈ 1.8 см, а во втором фор-сгусток отстоит от основного потока на z(A1)–z(A2) ≈ 4 см. В пользу предположения о нескольких идущих друг за другом сгустков свидетельствуют также данные с кадровых ЭОП-регистраторов (см. рис. 24) и радиальных оптических разверток (см. рис. 25).

Рис. 24.

Кадровые оптические изображения плазменного потока (импульсы № 4548 и № 4554, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж) по данным ЭОП-регистрации при прохождении плазмой позиции z = 35 см в момент времени topt относительно момента особенности на производной полного тока: а) – на topt = 4.7 мкс; б) – на topt = 4.4 мкс. Масштаб сетки 1 см × 1 см.

Рис. 25.

Результаты сравнения данных с оптической радиальной развертки (вверху) и магнитных зондов (внизу) в импульсе № 4535, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж: 13 – временные зависимости индукции тороидального магнитного поля, измеренной магнитными зондами на радиусах rp = 3 см на высотах zp = [30.4; 35.3; 39.9] см, соответственно. Пространственные щели оптической развертки располагались на высотах zopt = [30.5; 35.0; 39.5] см.

Например, из оптической щелевой развертки рис. 25 следует, что впереди основного тела плазменного потока могут лететь несколько фор-сгустков. В представленном случае в явном виде наблюдается до 3-х фор-сгустков. Один из них, обозначенный на всех развертках надписью “фор-сгусток”, расположен на расстоянии около 3 см от оси пролетной камеры. В том же месте пролетной камеры был расположен магнитный зонд, сигналы с которого синхронизованы с изображением оптической развертки и представлены ниже на этом же рисунке (кривые 13). Выбраны сигналы с 3-х из 9-ти катушек магнитного зонда, позиции которых zp = [30.4; 35.3; 39.9] см вдоль оси пролетной камеры соответствовали позициям пространственных щелей камеры оптической развертки zopt = [30.5; 35.0; 39.5] см. Из сравнения данных магнитных зондов и оптических разверток следует, что в момент пролета указанного на рисунке фор-сгустка положения зонда на радиусе rp = 3 см на кривых 13 присутствует пьедестал. Далее, магнитные зонды регистрируют сильное увеличение Bφ-поля до пикового значения ~0.6 кГс и последующий его спад практически до нуля. Пиковые значения магнитного поля регистрируются зондом во время пролета положения зонда наиболее ярко светящейся части плазменного потока. Следует заметить, что временной спад оптического излучения происходит быстрее, чем спад магнитного поля. В рассматриваемом случае через ~15 мкс после пролета первого плазменного потока регистрируется рост уровня Bφ‑поля выше 0.3 кГс, регистрируемого зондом в позиции zp = 30.4 см. В это же время на оптической развертке в позиции zopt = 30.5 см регистрируется увеличение оптического свечения, связанного с приходом второго плазменного потока.

Прохождение двух близко идущих друг за другом сгустков, как наблюдалось на аксиальных распределениях тороидального магнитного поля (см. напр. рис. 19), также было зарегистрировано на оптических щелевых развертках (см. рис. 26). Из сравнения зондового сигнала в позиции zp = = 30.4 см (кривая 1) и данных оптической развертки в позиции zopt = 30.5 см, прослеживается некоторая корреляция 2-х пиковой формы зондового сигнала с прохождением двух явно выраженных областей интенсивного свечения в плазменном потоке.

Рис. 26.

Результаты сравнения данных с оптической радиальной развертки (вверху) и магнитных зондов (внизу) в импульсе № 4536, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж: 13 – временные зависимости индукции тороидального магнитного поля, измеренной магнитными зондами на радиусах rp = 3 см на высотах zp = [30.4; 35.3; 39.9] см, соответственно. Пространственные щели оптической развертки располагались на высотах zopt = [30.5; 35.0; 39.5] см.

Попутно следует заметить, что на оптических развертках, зарегистрированных в других позициях вдоль оси пролетной камеры zopt = 35.0 см и 39.5 см не наблюдается прохождения двух областей повышенного свечения плазмы. При этом катушки магнитного зонда, расположенные в близких позициях zp = 35.3 см (кривая 2) и zp = = 39.9 см (кривая 3), регистрируют тоже только один пик роста Bφ-поля.

3.4. Взаимодействие плазменного потока с конструкционными особенностями пролетной камеры

В серии экспериментов, обсуждаемой в этом разделе статьи, дистанция распространения осевого струйного выброса была ограничена высотой пролетной камеры установки. При этом нижний край торцевого фланца пролетной камеры располагался на высоте z ≈ 51 см над поверхностью анода установки (см. рис. 2в). В предыдущих экспериментах было показано, что плазменный поток может распространяться на расстояния не менее z = 95 см [35], поэтому в данном случае происходило взаимодействие плазменного потока с торцевым фланцем сразу же в этой секции пролетной камеры.

На рис. 27–30 представлены результаты измерений магнитных полей как вблизи центрального сечения пролетной камеры (z ≈ 30–37 см) на радиусе rp = 3 см, так и вблизи торцевого фланца (z ≈ 44–50 см) на радиусе rp = 5 см. Как видно из распределений магнитного поля Bφ(3, z) рис. 28 и 30 в период времени t1t4 зондом, расположенным на радиусе rp = 3 см, регистрируется прохождение области плазменного потока с максимальным уровнем Bφ-поля (обозначен т. A) со скоростью Vz ≈ (2.1–3.5) × 106 см/с. В некоторых случаях, как было рассмотрено в п. 3.3, наблюдается прохождение двух областей с различным уровнем тороидального магнитного поля (обозначены как А1 и A2). В следующий период времени t5t9 на аксиальном распределении магнитного поля, измеренном зондом на rp = 5 см, также наблюдается прохождение области с пиковым уровнем магнитного поля ~0.3–0.4 кГс (обозначен т. B). Следует заметить, что на распределении Bφ(5, z) в области значений z = 46–50 см со временем развивается повышенный примерно в 2 раза уровень магнитного поля по сравнению с т. B. Данная область отмечена как Bφ“↑”. Такое кратное увеличение уровня Bφ-поля, захваченного плазменным потоком, связано со взаимодействием головной части потока с торцевым фланцем пролетной камеры вблизи позиции z ≈ 50 см. При таком взаимодействии происходит торможение плазмы на фронте потока (см. данные z(t) на рис. 16) и деформация токового контура, связанного с магнитным потоком Bφ-поля. При этом, согласно хотя бы частичному сохранению магнитного потока через данный контур, в нем увеличивается индукция тороидального магнитного поля. Кинетическая энергия направленного вдоль оси Z движения плазменного потока частично преобразуется в энергию магнитного поля. Такое неупругое взаимодействие плазменного потока с торцевым фланцем пролетной камеры приводит к его отражению в обратную сторону за счет давления сжатого потока Bφ-поля и внутреннего давления плазмы. Расширение плазменного потока приводит к уменьшению величины тороидального магнитного поля, которое наблюдается начиная с t > 20 мкс (см. семейство кривых Bφ(5, zj, t) на рис. 29). Следует заметить, что такое отражение плазменного потока от торцевого фланца пролетной камеры было обнаружено ранее в экспериментах на установке ПФ-3 [3135].

Рис. 27.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 50) см от поверхности анода установки (импульс № 4542, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): временные зависимости индукции тороидального магнитного поля, измеренной магнитными зондами на радиусах rp = 3 см [кривые Bφ(3, zi, t)] и rp = 5 см [кривые Bφ(5, zj, t)], где zi ∈ [30.4; 31.2; 32.0; 32.8; 33.6; 34.4; 35.3; 36.2; 37.1] см и zj ∈ [43.9; 44.6; 45.3; 46.2; 47.0; 47.8; 48.6; 49.5; 50.4] см. В легенде показано соответствие представленных кривых позициям расположения катушек зондов вдоль оси Z.

Рис. 28.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 50) см от поверхности анода установки (импульс № 4542, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): (а–и) – аксиальные распределения магнитного поля Bφ(z) на радиусах rp = 3 см и 5 см, восстановленные из кривых Bφ(3, zi, t) и Bφ(5, zj, t) рис. 27 на моменты времени t1t9, указанные стрелками (↓).

Рис. 29.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 50) см от поверхности анода установки (импульс № 4543, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): временные зависимости индукции тороидального магнитного поля, измеренной магнитными зондами на радиусах rp = 3 см [кривые Bφ(3, zi, t)] и rp = 5 см [кривые Bφ(5, zj, t)], где zi ∈ [30.4; 31.2; 32.0; 32.8; 33.6; 34.4; 35.3; 36.2; 37.1] см и zj ∈ [43.9; 44.6; 45.3; 46.2; 47.0; 47.8; 48.6; 49.5; 50.4] см. В легенде показано соответствие представленных кривых позициям расположения катушек зондов вдоль оси Z.

Рис. 30.

Результаты измерений тороидальных магнитных полей в плазменном потоке на пролетной базе z ∈ (30, 50) см от поверхности анода установки (импульс № 4543, газ Ne, P0 = 2.0 Торр, U0 = 9 кВ, W0 = 370 кДж): (а–и) – аксиальные распределения магнитного поля Bφ(z) на радиусах rp = 3 см и 5 см, восстановленные из кривых Bφ(3, zi, t) и Bφ(5, zj, t) рис. 29 на моменты времени t1t9, указанные стрелками (↓).

4. ВЫВОДЫ

Проведенные на установке ПФ-3 эксперименты по исследованию радиального и осевого распределений тороидального магнитного поля, захваченного осевым плазменным выбросом, позволяют сделать следующие выводы:

– Общий вид радиальных распределений тороидального магнитного поля Bφ(r), зарегистрированных в эксперименте, указывает на существование двух областей в плазменном потоке: области центрального тока, где Bφ(r) ~ r и области, где Bφ(r) ~ 1/rn, n ≈ 0.9–1.0;

– Характер изменения со временем радиального размера области с Bφ(r) ~ r свидетельствует об увеличении в диаметре канала протекания центрального тока в направлении от головной к хвостовой части плазменного потока. При этом зарегистрировано, что диаметр области центрального тока увеличивается с 4 см до величины более 7 см. Полученный результат находится в хорошем согласии с результатами, полученными ранее в экспериментах [24], а также с моделью, разработанной в работе [33]. Приближенно линейный характер распределения магнитного поля в виде функциональной зависимости Bφ(r) ~ r в области протекания центрального тока указывает на равномерное распределение плотности тока вдоль радиуса в данной области;

– Тороидальное магнитное поле на периферии плазменного потока регистрируется вплоть до радиуса r ~ 8 см. В хвостовой части плазменного потока область, занимаемая захваченным Bφ-полем занимает весь объем, ограниченный стенками пролетной камеры. Вероятней всего, обратные токи, протекающие на периферии плазменного потока, замыкаются на металлические стенки камеры (r ≈ 10.3 см);

– Измерения магнитных полей с пространственным разрешением вдоль оси Z пролетной камеры позволили определить протяженность области (Δz) с захваченным магнитным потоком. Минимальная величина Δz, зарегистрированная в этой серии экспериментов, составила величину около 5 см. В некоторых случаях, в приосевой области плазменного потока наблюдалось прохождение нескольких областей с повышенным уровнем Bφ-поля. Это могло быть связано как с прохождением двух близко идущих друг за другом плазменных сгустков, так и с прохождением основного плазменного потока с впередиидущим фор-сгустком;

– Оценка скорости перемещения Vz области с максимальным уровнем тороидального магнитного поля вдоль оси Z на участке 30–36 см составила величину (2.2–3.6) × 106 см/с, что близко к оценке скорости движения плазменного потока, определенной из обработки результатов оптических диагностик (коллиматоров и разверток);

– Во время движения осевого плазменного выброса происходит его торможение и диссипация магнитного потока, захваченного плазмой. Согласно аппроксимации экспериментальных измерений магнитного поля в виде $B_{\varphi }^{{{\text{max}}}}(z)$ = Bφ0 ⋅ · exp(–KBz) сделана оценка дистанции zmag пролетной базы, на которой происходит уменьшение величины $B_{\varphi }^{{{\text{max}}}}$ в e-раз. Полученное значение zmag ≈ 3.3 см оказалось на порядок меньшим, чем дистанция zopt ≈ 30 см, на которой происходит уменьшение скорости движения плазменного потока Vz в e-раз, определенной по данным с оптических коллиматоров;

– Неупругое взаимодействие плазменного потока с верхним торцевым фланцем пролетной камеры, приводящее к его частичному отражению, сопровождается сжатием магнитного потока Bφ‑поля. При этом происходит кратное увеличение уровня Bφ-поля, захваченного плазменным потоком. Происходящее при отражении плазмы от торцевого фланца расширение плазменного потока приводит к уменьшению величины тороидального магнитного поля.

Таким образом, детальное исследование распределения магнитного поля позволяет сделать вывод о возможности существования автономных плазменных сгустков с вмороженным магнитным полем, не связанных с “центральной машиной” – плазменным фокусом в нашем случае. Магнитный поток тороидального магнитного поля (Bφ-поля) оказывается отсоединен от основного контура разрядного тока установки. Скорее всего, это происходит в момент генерации осевого плазменного выброса. Для плазмы высокой проводимости, выбрасываемой вдоль оси установки, выполняется отношение Bφ/(ρ ⋅ r) ≈ const (условие вмороженности магнитного поля в плазму), где ρ – плотность плазмы, r – характерный размер объема, содержащего захваченный магнитный поток. По мере движения плазмы происходит затухание магнитного потока, вызванного конечной электропроводностью плазмы, за характерное время td ~ (2πσ/c2) ⋅ r2. Следует заметить, что в лабораторных условиях сохранение магнитного потока происходит в средах с высокой проводимостью σ. Однако в условиях космической плазмы более существенны большие характерные размеры r рассматриваемых контуров и, соответственно, большие времена затухания магнитного поля по сравнению с временем изучаемого процесса. Тем не менее, полученные данные могут быть использованы для объяснения характера структур, наблюдаемых в астрофизических джетах и построения соответствующих физических моделей.

Работа выполнена при финансовой поддержке РФФИ (грант № 18-29-21006).

Список литературы

  1. Ryutov D., Drake R.P., Kane J., Liang E., Remington B.A., Wood-Vasey W.M. // The Astrophysical J. 1999. V. 518. № 2. P. 821.

  2. Remington B.A., Drake R.P., Ryutov D.D. // Rev. Mod. Phys. 2006. V. 78. P. 755.

  3. Ryutov D.D. and Remington B.A. // Plasma Phys. Control. Fusion. 2002. V. 44. № 12B. P. 407.

  4. Бескин В.С. Осесимметричные стационарные течения в астрофизике. М.: ФИЗМАТЛИТ, 2005. 384 с.

  5. Albertazzi B., Ciardi A., Nakatsutsumi M., Vinci T., Beard J., Bonito R., Billette J., Borghesi M., Burkley Z., Chen S.N., Cowan T.E., Herrmannsdorfer T., Higgin-son D.P., Kroll F., Pikuz S.A., Naughton K., Romagna-ni L., Riconda C., Revet G., Riquier R., Schlenvoigt H.P., Skobelev I.Y., Faenov A.Y., Soloviev A., Huarte-Espino-sa M., Frank A., Portugall O., Pepin H., Fuchs J. // Science. 2014. V. 346. I. 6207. P. 325.

  6. Belyaev V.S., Zagreev B.V., Lobanov A.V., Matafonov A.P., Bisnovatyi-Kogan G.S., Moiseenko S.G., Toropina O.D., Gromov A.I. // Astronomy Reports. 2018. V. 62. P. 162.

  7. Lebedev S.V., Frank A., Ryutov D.D. // Rev. Mod. Phys. 2019. V. 91. P. 025002.

  8. Byvank T., Banasek J.T., Potter W.M., Greenly J.B., Seyler C.E., Kusse B.R. // Phys. Plasmas. 2017. V. 24. P. 122701.

  9. Bellan P.M. // J. Plasma Phys. 2018. V. 84. P. 755840501.

  10. Lavine E.S., You S. // Phys. Rev. Lett. 2019. V. 123. P. 145002.

  11. Haro G. // Astron. J. 1950. V. 55. P. 72.

  12. Herbig G.H. // Astrophys. J. 1950. V. 111. P. 11.

  13. Reipurth Bo, Bally J. // Annual Rev. Astron. Astrophys. 2001. V. 39. P. 403.

  14. Krauz V., Myalton V., Vinogradov V., Velikhov E., Ananyev S., Dan’ko S., Kalinin Yu., Kharrasov A., Mitrofanov K., Vinogradova Yu. // Proc. 42nd EPS Conf. Plasma Physics, Lisbon, Portugal, 2015. V. 39E. ISBN 2-914771-98-3. P. 4.401.: http://ocs.ciemat.es/EPS2015PAP/pdf/P4.401.pdf

  15. Krauz V.I., Beskin V.S., Velikhov E.P. // International Journal of Modern Physics D. 2018. V. 27. P. 1844009.

  16. Полухин С.Н., Джаманкулов А.М., Гурей А.Е., Никулин В.Я., Перегудова Е.Н., Силин П.В. // Физика плазмы. 2016. Т. 42. С. 1080.

  17. Krauz V.I., Paduch M., Tomaszewski K., Mitrofanov K.N., Kharrasov A.M., Szymaszek A., Zielinska E. // Europhysics Letters (EPL). 2020. V. 129. P. 15003. https://doi.org/10.1209/0295-5075/129/15003

  18. Ананьев С.С., Данько С.А., Мялтон В.В., Жужунашвили А.И., Калинин Ю.Г., Крауз В.И., Ладыги-на М.С., Марченко А.К. // Физика плазмы. 2016. Т. 42. С. 282.

  19. Skladnik-Sadowska E., Dan’ko S.A., Kwiatkowski R., Sadowski M.J., Zaloga D.R., Paduch M., Zielinska E., Kharrasov A.M., Krauz V.I. // Phys. Plasmas. 2016. V. 23. P. 122902.

  20. Бескин В.С., Истомин Я.Н., Киселев А.М., Крауз В.И., Митрофанов К.Н., Мялтон В.В., Нохрина Е.Е., Собьянин Д.Н., Харрасов А.М. // Изв. вузов. Радиофизика. 2016. Т. 59. С. 1004.

  21. Kalashnikov I., Chardonnet P., Chechetkin V., Dodin A., Krauz V. // Physics of Plasmas. 2018. V. 25. P. 062901.

  22. Митрофанов К.Н., Крауз В.И., Мялтон В.В., Велихов Е.П., Виноградов В.П., Виноградова Ю.В. // ЖЭТФ. 2014. Т. 146. С. 1035.

  23. Митрофанов К.Н., Крауз В.И., Мялтон В.В., Виноградов В.П., Харрасов А.М., Виноградова Ю.В. // Астрономический журнал. 2017. Т. 94. С. 152.

  24. Крауз В.И., Митрофанов К.Н., Войтенко Д.А., Астапенко Г.И., Марколия А.И., Тимошенко А.П. // Астрономический журнал. 2019. Т. 96. С. 156.

  25. Krauz V.I., Mitrofanov K.N., Paduch M., Tomaszewski K., Szymaszek A., Zielinska E., Pariev V.I., Beskin V.S., Istomin Ya.N. // J. Plasma Phys. 2020. V. 86. P. 905860607.

  26. Filippov N.V., Filippova T.I., Khutoretskaia I.V., Mial-ton V.V., Vinogradov V.P. // Phys. Lett. A. 1966. V. 211. № 3. P. 168.

  27. Митрофанов К.Н., Крауз В.И., Грабовский Е.В., Мял-тон В.В., Виноградов В.П., Падух М., Шольц М., Карпински Л. // Физика плазмы. 2015. Т. 41. С. 413.

  28. Глазырин И.В., Грабовский Е.В., Зукакишвили Г.Г., Карпеев А.В., Митрофанов К.Н., Олейник Г.М., Самохин А.А. // Вопросы атомной науки и техники (ВАНТ). Сер. Термоядерный синтез. 2009. Вып. 2. С. 67.

  29. Митрофанов К.Н., Крауз В.И., Грабовский Е.В., Мял-тон В.В., Падух М., Грицук А.Н. // Приборы и техника эксперимента (ПТЭ). 2018. Т. 61. № 2. С. 78.

  30. Виноградов В.П., Крауз В.И., Мокеев А.Н., Мял-тон В.В., Харрасов А.М. // Физика плазмы. 2016. Т. 42. С. 1033.

  31. Ананьев С.С., Данько С.А., Мялтон В.В., Калинин Ю.Г., Крауз В.И., Виноградов В.П., Виноградова Ю.В. // Вопросы атомной науки и техники (ВАНТ). Сер. Термоядерный синтез. 2013. Т. 36. № 4. С. 102.

  32. Крауз В.И., Митрофанов К.Н., Харрасов А.М., Ильичев И.В., Мялтон В.В., Ананьев С.С., Бескин В.С. // Астрономический журнал. 2021. Т. 98. С. 29.

  33. Бескин В.С., Калашников И.Ю. // Письма в Астрономический журнал (ПАЖ). 2020. Т. 46. № 7. С. 494.

  34. Крауз В.И., Войтенко Д.А., Митрофанов К.Н., Мялтон В.В., Аршба Р.М., Астапенко Г.И., Марколия А.И., Тимошенко А.П. // ВАНТ .Сер. Термоядерный синтез. 2015. Т. 38. Вып. 2. С. 19.

  35. Крауз В.И., Виноградов В.П., Мялтон В.В., Виноградова Ю.В., Харрасов А.М. // Вопросы атомной науки и техники. ВАНТ. Сер. Термоядерный синтез. 2018. Т. 41. Вып. 3. С. 48.

Дополнительные материалы отсутствуют.