Том 50, № 12 (2024)
ТОКАМАКИ
О применимости теории ГАМ к реальным токамакам с гофрированным магнитным полем
Аннотация
Исследовано влияние тороидальной гофрировки магнитного поля на пространственную структуру и частоту геодезической акустической моды (ГАМ) в плазме токамака. Показано, что тороидальная асимметрия магнитной конфигурации приводит к зацеплению колебаний электрического потенциала ГАМ с тороидально- и полоидально-неоднородными возмущениями давления плазмы. Для токамаков с большим аспектным отношением выведен закон дисперсии ГАМ, учитывающий неоднородность гофрировки в сечении токамака. При увеличении числа катушек n тороидального поля эффект гофрировки перестает зависеть от n. Показана применимость стандартной теории для нахождения частоты и пространственной структуры ГАМ в крупных токамаках.



Оптимизация дивертора ТРТ на основе моделирования кодом SOLPS-ITER
Аннотация
На основе моделирования анализируются возможные режимы работы дивертора и параметры пристеночной плазмы проектируемого токамака ТРТ. Показано, что при мощности на сепаратрисе 18 МВт, соответствующей приблизительно вдвое большей вложенной мощности, можно обеспечить низкую плотность потока тепла на пластины дивертора 5МВт/м2 при концентрации электронов на сепаратрисе менее 7 × 1019 м–3 и эффективном заряде менее 2, что является реалистичным сценарием для этой установки. В случае повышения мощности на сепаратрисе возможен режим работы с более высокой тепловой нагрузкой на пластины дивертора, все еще остающейся в технологических пределах. Моделирование также показывает желательность увеличения расстояния между сепаратрисой и элементами камеры и преимущество работы с угловой конфигурацией внешнего дивертора по сравнению с “итероподобной”.



Исследование деталей профиля электронной температуры при сдвиге плазменного шнура и мощности ЭЦРН до 1.5 МВт в токамаке Т-10
Аннотация
Вопрос о существовании внутренних транспортных барьеров (ВТБ) около низкорациональных значений коэффициента запаса устойчивости q = 1, 2, ... в обычной L-моде токамаков с центральным дополнительным нагревом и положительным магнитным широм остается открытым. Для прояснения существования таких ВТБ были проведены эксперименты с запрограммированным линейном по времени сдвиге плазменного шнура, при котором позиции измерения электронной температуры Te по излучению на второй гармонике ЭЦР сдвигаются относительно плазменного шнура, что позволяет изучить детали профиля Te. На токамаке Т-10 были проведены серии экспериментов с перпендикулярным вводом ЭЦ-излучения мощностью 0.4 и 0.85 МВт при центральном ЭЦРН и быстрым сдвигом шнура на 0.13а (малого радиуса плазмы а) за 60 мс в плазме с углеродным лимитером. В обоих случаях снаружи поверхности q = 1 обнаружены узкие (шириной около 0.03а) и слабые (спад коэффициента электронной теплопроводности χе примерно вдвое) ВТБ, исчезающие при изменении параметров разряда. Эти ВТБ в 2–3 раза уже и на порядок слабее, чем ВТБ около поверхности q = 1, предложенные ранее для объяснения результатов токамака RTP. ВТБ не обнаружены в L-моде в экспериментах с одновременной генерацией кo-тока и контртока в центре плазмы двумя гиротронами общей мощностью 1.5 МВт в плазме с вольфрамовым лимитером (сдвиг шнура на 0.1а за 30 мс). По-видимому, данная статья является первой журнальной публикацией про исследование деталей профиля Te при быстром сдвиге плазменного шнура.



Моделирование переноса тепла и частиц в токамаках COMPASS и Т-10 с помощью транспортной модели канонических профилей
Аннотация
Представлены результаты моделирования переноса тепла и частиц плазмы, получаемой в режиме омического нагрева в токамаке Т-10 с круглым лимитером и D-образной плазмы с дивертором в токамаке COMPASS. Проведено моделирование в режиме H-моды с омическим нагревом плазмы и при ее дополнительном нагреве в результате инжекции нейтрального пучка (NBI) в токамаке COMPASS. Моделирование осуществлялось с помощью Транспортной модели канонических профилей (ТМКП) с использованием кода ASTRA. Полученные радиальные профили электронной температуры и плотности плазменного шнура согласуются с измеренными со среднеквадратичными отклонениями в пределах точности эксперимента 10–15٪. Расчеты продемонстрировали весьма схожие профили плотности плазмы, как в омическом режиме, так и при дополнительном NBI-нагреве в H-моде. Профили электронной температуры в H-моде с дополнительным нагревом имеют более высокие пьедесталы, чем в омической H-моде, что согласуется с измерениями. Сравнение показало, что омические режимы в COMPASS и Т-10 можно описать одинаковыми коэффициентами жесткости в уравнениях переноса тепла и частиц в плазме.



Применение двухдлинноволновой цифровой голографической интерферометрии для оптического неразрушающего контроля износа защитных элементов сферического токамака ГЛОБУС-М2
Аннотация
Продемонстрирована возможность применения метода двухдлинноволновой цифровой голографической интерферометрии для оценки износа защитных элементов сферического токамака Глобус-М2 после рабочих плазменных разрядов. На данном этапе работ защитные элементы были извлечены из разрядной камеры токамака и использовались в качестве образцов в голографическом комплексе. Представлена схема голографического интерферометра для регистрации первичных голографических изображений, в котором управление системами регистрации и контроля длины волны излучения осуществляется через программно-аппаратный комплекс в реальном масштабе-времени. Представлены результаты измерений формы элементов токамака. Показано, что при изменении разности длин волн меняется чувствительность метода измерений, а в предложенной конфигурации оптической схемы возможно определить минимальное значение изменения формы на уровне 10–30 мкм. При этом ошибка определения разности фаз, по которой проводится оценка профиля поверхности, в цифровом методе может достигать порядка 2π/40.



МАГНИТНЫЕ ЛОВУШКИ
Полноволновое моделирование электронного циклотронного нагрева плазмы на первой и второй гармониках в установке ГДМЛ
Аннотация
Проанализированы возможные сценарии дополнительного электронного циклотронного нагрева плазмы в различных конфигурациях проектируемой в ИЯФ СО РАН открытой магнитной ловушке следующего поколения ГДМЛ (газодинамическая многопробочная ловушка). Для этого был применен полноволновой импедансный подход к моделированию взаимодействия электромагнитных волн с субтермоядерной плазмой, позволяющий учитывать взаимодействие электромагнитных и квазиэлектростатических мод в окрестности электронного циклотронного резонанса в осесимметричной магнитной конфигурации. Рассмотрены схемы нагрева на первой гармонике гирочастоты обыкновенной электромагнитной волной и на второй гармонике необыкновенной волной. Для каждого случая определены диапазоны параметров установки, в которых такая схема нагрева может быть эффективной; проанализированы оптимальные с точки зрения эффективности нагрева параметры фокусировки вводимого квазиоптического СВЧ-пучка.



ДИНАМИКА ПЛАЗМЫ
Плазмообразование на поверхности конденсированного вещества под действием мощного рентгеновского излучения
Аннотация
В ряде экспериментов поверхности конденсированного вещества, например, электроды электрофизических установок, подвергаются воздействию мощного импульсного рентгеновского излучения с плотностью потока энергии ~1 ТВт/см2. Источником этого излучения могут быть, например, Z-пинчи, образованные при токовом сжатии многопроволочных лайнеров. Под действием этого излучения на поверхности электродов могут идти процессы испарения и плазмообразования. В работе проведено теоретическое рассмотрение этих процессов. В случае когда толщина слоя плазмы мала по сравнению с характерными размерами электродов, плазмообразование можно описывать одномерными уравнениями магнитной гидродинамики с учетом переноса излучения. Одномерные расчеты, проведенные для условий экспериментов, проведенных на установке Ангара-5-1 (плотность потока энергии, приходящая от пинча, ~0.2 ТВт/см2, время действия излучения ~15 нс, материал электродов Fe), показали, что характерная температура плазмы составляет в этом случае ~40 эВ, плотность ~3 мг/см3, а ее скорость разлета ~60 км/с. Интересно, что магнитные поля в этих экспериментах, сами по себе небольшие (~0.8 МГс) и неспособные привести к плазмообразованию, своим давлением сдерживают разлет плазмы и влияют на ее характерные величины и скорость разлета. Полученная в расчете скорость оказалась несколько меньше измеренной экспериментально с помощью рентгеновского электронно-оптического преобразователя (~90 км/с), что, возможно, связано с неодномерным турбулентным разлетом плазмы или погрешностями эксперимента.



Газодинамическая и кинетическая стадии газового Z-пинча
Аннотация
Рассматривается феноменологическая газодинамическая модель сжатия перетяжки газового Z-пинча, через торцы которой плазма вытекает с большой скоростью. Проведенные расчеты показали, что создаются условия, при которых релаксация ионного компонента плазмы запаздывает на фоне макроскопической динамики сжатия. Поэтому описание Z-пинчей на стадии максимального сжатия должно учитывать кинетику ионов. Такой подход может раскрыть механизм ускорения ионов до высоких энергий и объяснить высокую интенсивность нейтронного излучения на финальной стадии сжатия перетяжки.



ИЗЛУЧЕНИЯ ПЛАЗМЫ
Усиление излучения слоем фотоионизованной плазмы инертного газа в магнитном поле
Аннотация
Изучено взаимодействие циркулярно поляризованной электромагнитной волны со слоем фотоионизованной плазмы инертного газа, находящегося в магнитном поле. Дан детальный анализ коэффициентов отражения и прохождения волны в условиях, когда частота волны совпадает с циклотронной частотой фотоэлектронов. Выявлена возможность сильного увеличения коэффициентов отражения и прохождения при отрицательных малых по абсолютной величине значениях мнимой части диэлектрической проницаемости фотоионизованной плазмы. Показано, что в слое фотоионизованной плазмы, полученной в процессе быстрой многофотонной ионизации атомов ксенона при атмосферном давлении, существует возможность усиления напряженности поля терагерцового излучения более чем на два порядка.



ДИАГНОСТИКА ПЛАЗМЫ
Взаимодействие мощного потока водородной плазмы со сверхзвуковой газовой струей и вольфрамовой мишенью
Аннотация
Приведены результаты исследования взаимодействия мощного потока водородной плазмы со сверхзвуковой газовой струей перед вольфрамовой мишенью. Азот или неон, инжектированные перед поверхностью мишени, служат надежным способом защиты от прямого воздействия водородной плазмы на вольфрам. Экспериментально показано, что образующаяся плазма газовой струи является мощным источником коротковолнового линейчатого излучения. Плотность поглощаемой вольфрамовой мишенью энергии ≈ 25 Дж/см2 оказывается вдвое меньше энергии, поглощаемой вольфрамом при импульсном воздействии потока водородной плазмы без газовой струи ≈ 50 Дж/см2. Максимально достигаемая поверхностью вольфрама температура составляет ≈ 3700 К с применением газовой струи и ≈ 5800 К – без газовой струи. Наличие газовой струи-завесы перед вольфрамом приводит к локализации испаренного вольфрама вблизи мишени на расстояниях до 1 см от поверхности.



НИЗКОТЕМПЕРАТУРНАЯ ПЛАЗМА
Упрощенный подход к описанию кинетических свойств примесных ионов в слабоионизованной плазме гелия
Аннотация
С помощью приближенного подхода, справедливого при большом отношении массы иона к массе нейтральной частицы, вычислена подвижность ионов , , и в гелии и получено хорошее согласие с имеющимися результатами расчетов методом Монте-Карло и экспериментальными данными при высоких значениях приведенного электрического поля. Этот упрощенный подход использован для определения средней энергии ионов и константы скорости диссоциации этих ионов в гелии. Проведено сравнение полученных характеристик ионов с результатами более точных расчетов методом Монте-Карло. Получено хорошее согласие между этими двумя подходами для средней энергии ионов, но для константы скорости диссоциации различие оказывается достаточно серьезным, достигающим порядка величины и более. Показано, что это различие связано с особенностями энергетического распределения ионов в случае большого отношения массы иона к массе нейтральной частицы.



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


