Теплофизика высоких температур, 2021, T. 59, № 3, стр. 323-329
Подавление плазменных неустойчивостей в плазме трехэлектродного стабилизатора тока и напряжения
А. С. Мустафаев 1, А. Ю. Грабовский 1, *, В. С. Сухомлинов 2
1 Санкт-Петербургский горный университет
Санкт-Петербург, Россия
2 Санкт-Петербургский государственный университет
Санкт-Петербург, Россия
* E-mail: schwer@list.ru
Поступила в редакцию 02.07.2020
После доработки 18.11.2020
Принята к публикации 22.12.2020
Аннотация
Экспериментально исследованы условия возбуждения колебаний тока и напряжения в плазме трехэлектродного стабилизатора тока и напряжения. Установлено, что в рассматриваемых режимах плазма обладает отрицательной проводимостью, которая и обусловливает колебания. Предложен высокоэффективный метод подавления плазменных неустойчивостей, основанный на управлении знаком дифференциальной проводимости плазмы путем регулировки концентрации тепловых электронов плазмы с помощью внешнего электрода. Предлагаемый метод позволяет достигать высокого уровня стабильности энергетических характеристик стабилизатора при плотности разрядного тока до 5 А/см2 и мощности 500 Вт/см2.
ВВЕДЕНИЕ
Одним из ключевых направлений развития систем вооружений в Российской Федерации является разработка малогабаритных сверхмощных ядерных энергетических установок (ЯЭУ) [1]. Для решения этой задачи требуются электронные компоненты, позволяющие реализовать полное управление током в цепях ЯЭУ. К таким компонентам (приборам) выставляются следующие требования на допустимые уровни излучений реактора на стенке приборного отсека: флюенс быстрых нейтронов (Еn > 0.1 МэВ) должен быть не более 1012 н/см2; поглощенная доза фотонов (γ-квантов) не должна превышать 106 рад. Таким требованиям удовлетворяют твердотельные полупроводниковые приборы, но с ростом электрической мощности и рабочей температуры применение для них локальной радиационной защиты резко ухудшает массогабаритные характеристики ЯЭУ. В этих условиях единственным решением является использование газоразрядных электронных приборов на базе неравновесной плазмы [2, 3].
Эксплуатация приборов плазменной энергетики сопряжена с проблемой возбуждения неустойчивостей при попытках повышения их энергетических параметров [4, 5]. Колебания и неустойчивости такого типа можно использовать в практических целях, однако их генерация в условиях плазменных источников и стабилизаторов оказывает разрушительное воздействие на их энергетические и электрокинетические характеристики. Очевидно, что борьба с неустойчивостями подразумевает управление функцией распределения электронов по скоростям (ФРЭС) в рабочих режимах плазменных приборов.
В настоящей работе для решения проблемы подавления неустойчивостей разработана специальная конструкция триода на базе гелиевого низковольтного пучкового разряда (НПР) с управляющим электродом (УЭ), расположенным вне зазора катод–анод. Такая конструкция позволяет управлять функцией распределения заряженных частиц для эффективного подавления возникающих плазменных неустойчивостей и колебаний.
ЭКСПЕРИМЕНТАЛЬНЫЕ УСТАНОВКА, ПРИБОР И МЕТОД ИССЛЕДОВАНИЯ
Схема экспериментальной установки и конструкция плазменного триода представлены на рис. 1 и 2 соответственно. Их детальное описание можно найти в работе [2]. Рассмотрим основные элементы установки.
В блок I входят вакуумная камера 1, катодный 2 и анодный 3 узлы. Специальные окна из сапфира 4 использовались для наблюдений и оптических измерений. Вакуумная система (блок II) включает турбомолекулярный насос, обеспечивающий максимально достижимое разрежение 5 × 10–8 Тор. Высокая стабильность разрядных условий достигалась длительным циклом обезгаживания при температуре 700 К.
Схема экспериментального прибора представлена на рис. 2. Его катод 1 изготовлен из пористого вольфрама, импрегнированного барием, и нагревается тантал-ниобиевой проволокой. Анод 2 и управляющий электрод 3 выполнены из поликристаллического молибдена. Диаметр катода – 10 мм, диаметры анода и УЭ равны и составляют 30 мм. Компоновка триода нетрадиционна – анод имеет центральное отверстие диаметром 2 мм, при этом УЭ вынесен из межэлектродного зазора и располагается непосредственно за анодом. Зазор катод–анод составляет 8 мм, УЭ расположен на расстоянии 1 мм за анодом. Плазменный канал триода ограничен металлическим конусным экраном 9, находящимся под потенциалом катода.
Диагностика плазмы НПР осуществлялась методом плоского одностороннего зонда, развитым в работах [6–8]. Поскольку в рассматриваемых экспериментальных условиях плазма обладает осевой симметрией, ФРЭC не зависит от азимутального угла (рис. 3):
В формуле (1) $\upsilon = \left| {\mathbf{\upsilon }} \right|;$ θ − полярный угол.
Ток на плоский ленгмюровский зонд рассчитывается как
где ${{\upsilon }_{n}}$ − нормальная к поверхности зонда составляющая вектора скорости электрона ${\mathbf{\upsilon }},$ U − потенциал зонда, $\varepsilon = {{m{{\upsilon }^{2}}} \mathord{\left/ {\vphantom {{m{{\upsilon }^{2}}} 2}} \right. \kern-0em} 2},$ $\phi {\kern 1pt} '$ и $\theta {\kern 1pt} '$ − соответственно азимутальный и полярный углы вектора ${\mathbf{\upsilon }}.$ Вторая производная зондового тока (2) по потенциалу зонда U соответствуетВыражение (3) можно переписать в виде
Для нахождения ФРЭС представим f(ε, θ) и в виде рядов по полиномам Лежандра:
(5)
$f\left( {\varepsilon ,\theta } \right) = \sum\limits_{j = 0}^\infty {{{f}_{j}}\left( \varepsilon \right){{P}_{j}}\left( {\cos \theta } \right)} ,$После подстановки (5) и (6) в (4) получаем соотношение между компонентами ${{f}_{j}}$ и ${{F}_{j}}{\text{:}}$
(7)
${{f}_{j}}(qU) = {{F}_{j}}(qU) + \int\limits_{qU}^\infty {{{f}_{j}}(\varepsilon )\frac{\partial }{{\partial (qU)}}{{P}_{j}}\left( {\sqrt {\frac{{qU}}{\varepsilon }} } \right)d\varepsilon } .$Выражение (7) является интегральным уравнением Вольтерра II рода. Решая его с помощью метода резольвент [9], получим
(8)
${{f}_{j}}(qU) = {{F}_{j}}(qU) + \int\limits_{qU}^\infty {{{F}_{j}}(\varepsilon ){{R}_{j}}(qU,\varepsilon )d\varepsilon } .$Подставим в уравнение (8) соотношение
и получим основную формулу, позволяющую реконструировать компоненты ${{f}_{j}}{\text{:}}$Таким образом, для диагностики плазмы методом плоского зонда необходимо зарегистрировать значения (qU, α) при различных углах зонда относительно оси разряда, рассчитать по формуле (9) ряд лежандровых компонентов и реконструировать полную ФРЭС согласно формуле (5). Формула (9) демонстрирует, что метод не требует априорной информации об анизотропии распределения заряженных частиц, однако для корректного описания ФРЭС в сильнонеравновесной плазме необходимо большое число компонент, что не всегда возможно в реальном эксперименте. В этой связи метод плоского зонда усовершенствован авторами для диагностики моноэнергетичных пучков заряженных частиц [10].
В процессе зондовых измерений тщательно выдерживались все требования теории Ленгмюра. С этой целью изготавливались зонды толщиной 30 мкм и диаметром 0.5 мм. Также учитывались все факторы, способные искажать зондовые вольт-амперные характеристики (ВАХ) [11].
Для получения значений второй производной зондового тока использовался метод демодуляции, реализованный в измерительном комплексе на базе ПК. В качестве дифференцирующего сигнала для увеличения чувствительности метода использовалось 100% модулированное напряжение ΔU = U0(1 + cosω1t)cos ω2t [11].
ЭКСПЕРИМЕНТАЛЬНЫЕ РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ
Распределение параметров в плазме НПР. В эксперименте роль плазмообразующего компонента играл гелий, поскольку он обладает самыми высокими среди других инертных газов потенциалами ионизации и возбуждения. Это позволяет добиться наиболее яркого проявления нелокальных эффектов в плазме. Экспериментальные исследования проводились в диапазоне токов 0.1−2 А и давлений гелия 0.1−10 Тор. Выбор таких условий обусловлен возможностью перевода плазмы как в локальный, так и в нелокальный режим. В последнем случае длины релаксации на электрон-электронных Lee и электрон-атомных Lеа столкновениях превосходят длину межэлектродного зазора d [12, 13].
Вторые производные плоского зонда, зарегистрированные при его различных ориентациях относительно оси разряда при PHe = 2.5 Тор показаны на рис. 4а. Видно, что структура ФРЭС формируется двумя обособленными группами электронов: тепловыми ft и пучковыми f0 с концентрацией и средней энергией nt, $\overline {{{\varepsilon }_{t}}} $ и n0, $\overline {{{\varepsilon }_{0}}} $ соответственно.
Распределение тепловых электронов по скоростям близко к изотропному максвелловскому, тогда как группа быстрых электронов отличается значительной анизотропией распределения. Сравнение энергий тепловых и пучковых электронов выявляет сильную неравновесность ФРЭС, энергия тепловых электронов – примерно 2 эВ, быстрых – 25 эВ.
Формирование такой структуры ФРЭС связано с распределением потенциала по межэлектродному зазору (рис. 4б, кривая 1). У анода образуется скачок потенциала φа ≈ $1.5\overline {{{\varepsilon }_{t}}} ,$ являющийся практически непреодолимым барьером для тепловых электронов плазмы. Электроны, эмитированные катодом, ускоряются на прикатодном скачке потенциала и образуют пучок с небольшим разбросом по энергии, пронизывающий межэлектродный зазор. Рождение медленных электронов происходит в результате ионизации атомов гелия электронами пучка. Потенциалы возбуждения и ионизации для He составляют Uм ≈ 19.8 В и Uион ≈ 24.6 В соответственно, что обусловливает доминирующую роль пучка в неупругих процессах и токопереносе.
Если рассматривать промежуток между анодом и управляющим электродом, то отверстие в аноде является неким подобием плазменного катода и токоперенос здесь осуществляется в основном тепловыми электронами из плазмы в зазоре катод–анод, что становится возможным благодаря их ускорению в сильном электрическом поле двойного слоя в окрестностях анода (рис. 4б, кривые 2, 3).
С целью изучения процессов, протекающих в рассматриваемых плазменных условиях, по зарегистрированным значениям восстановлены лежандровы компоненты разложения ФРЭС f0–f6, определяющие набор важнейших базовых параметров плазмы: концентрацию, плотность электронного тока, анизотропию электронного давления и др. [11]. В частности, по компоненту f0 можно рассчитать функцию возбуждения атомов и генерации ионов:
(10)
${\text{Г}} = 4\pi {{N}_{a}}\int\limits_{{{v}_{{{\text{пор}}}}}}^\infty {{v}\sigma _{{eа}}^{i}} ({v}){{f}_{0}}({v}){{{v}}^{2}}d{v},$Подавление колебаний тока и напряжения в плазменных приборах с отрицательной проводимостью. Известно, что НПР может быть неустойчив к возбуждению различного типа колебаний [12–14], разрушающих рабочие режимы плазменных электронных приборов.
Вольт-амперные характеристики гелиевого НПР, зарегистрированные в условиях трехэлектродного стабилизатора тока и напряжения, представлены на рис. 5. Кривая 1 имеет участок отрицательной дифференциальной проводимости Gд, что является причиной развития колебаний разрядного тока и напряжения [15].
Типичная форма таких колебаний представлена на рис. 6. Колебания происходят при полной модуляции тока в разряде на частотах 50−150 кГц и при амплитудных значениях напряжения порядка 30 В. В условиях постоянного давления частота колебаний практически не изменяется при варьировании уровня разрядного тока. Если ток возрастает, то глубина модуляции падает вплоть до 10%. При увеличении давления наполнителя частота колебаний нарастает линейно (рис. 7).
Аксиальное распределение параметров плазмы перед возникновением колебаний приведено в таблице. Функция генерации Г (показывает число ионов, родившихся в единице объема за 1 с) рассчитывалась по формуле (10) с помощью измеренного лежандрова компонента ФРЭС f0. Величина D = n/τ определяет уход ионов из плазменного канала в радиальном направлении, где τ ~ ${r \mathord{\left/ {\vphantom {r {{{{\bar {\upsilon }}}_{{{\text{др}}}}}}}} \right. \kern-0em} {{{{\bar {\upsilon }}}_{{{\text{др}}}}}}}$ − время радиальной диффузии ионов, скорость движения ионов ${{\bar {\upsilon }}_{{{\text{др}}}}}$ в радиальном поле напряженностью Еr ~ 5 В/см составляла величину порядка 2 × 105 см/с [16].
Z, мм | φ, В | n0, 109 см−3 | nt, 1011 см−3 | τдиф, с | D ~ n/τ, 1017 см−3 c−1 | Г, 1017 см−3 c−1 | Г/D |
---|---|---|---|---|---|---|---|
1 | 20.1 | 15 | 0.5 | 4.0 | 0.1 | 6.0 | 60.0 |
2 | 20.7 | 12 | 0.75 | 5.5 | 0.13 | 4.8 | 37.0 |
3 | 21.2 | 7 | 1.23 | 7.0 | 0.17 | 2.8 | 16.4 |
4 | 21.4 | 4 | 1.6 | 8.5 | 0.2 | 1.6 | 8.0 |
5 | 22 | 2.74 | 2.5 | 10 | 0.24 | 1.1 | 4.58 |
6 | 22.4 | 1.75 | 3.73 | 11.5 | 0.27 | 0.7 | 2.59 |
7 | 22.8 | 1.11 | 5.57 | 13 | 0.3 | 0.44 | 1.47 |
8 | 23.3 | 0.71 | 8.29 | 14.5 | 0.34 | 0.28 | 0.82 |
Из таблицы видно, что по мере удаления от катода в аксиальной области плазмы несколько возрастает D из-за увеличения радиального градиента концентрации и поля. Вместе с тем вследствие влияния радиальной диффузии и неупругих процессов происходит обеднение быстрой части ФРЭС, влекущее снижение функции генерации ионов Г. Вблизи анода величины Г и D сопоставимы (рис. 8), что дает возможность записать условие стационарности концентрации ионов в виде Г ≈ D. В случае нарушения этого условия (при Г ≥ D) формируется избыточное количество ионов в прианодном регионе с последующим возрастанием тока медленных электронов на анод и возникновением релаксационных колебаний, характер которых зависит еще и от анодной нагрузки. При Г ≤ D дефицит ионов может приводить к обрыву тока.
Таким образом, физическую картину возбуждения колебаний можно сформировать, рассматривая соотношение процессов генерации ионов в плазме и их ухода из межэлектродного зазора. Из формулы (10) видно, что с увеличением давления гелия функция генерации возрастает и диффузионный поток ионов на анод увеличивается (время пробега ионов в зазоре катод–анод порядка миллисекунд, что хорошо согласуется с частотами колебаний). В результате происходит компенсация прианодного объемного заряда тепловых электронов, ранее запертых в потенциальной яме между катодом и анодом. В результате их хаотический ток на анод it возрастет. В предположении нулевого прианодного потенциального барьера величину it можно записать в виде
где e − заряд электрона, ne − концентрация электронов вблизи анода, υt − тепловая скорость электронов, S − площадь анода. Оценки при PHe = 5 Тор, ne ~ 1011 см−3, Te ~ 1 эВ дают it ~ 1.8 A, что приближенно соответствует зарегистрированному значению (рис. 6).На основании проведенных исследований предложен способ подавления колебаний тока и напряжения [17]. Из рис. 5 видно, что отбор тока тепловых электронов на управляющий электрод позволяет управлять знаком Gд. Возможность такого управления позволяет эффективно подавлять колебания тока и напряжения, поскольку во всех режимах с положительной дифференциальной проводимостью (рис. 5, кривые 2−5) колебания отсутствовали и более не возбуждались.
ЗАКЛЮЧЕНИЕ
В работе представлены результаты зондовых исследований ФРЭС и распределения параметров в пучковой плазме трехэлектродного стабилизатора тока и напряжения. Показано, что в плазме присутствуют две группы электронов разных энергий и концентраций – быстрые электроны пучка и тепловые электроны плазмы. Экспериментально обнаружены колебания разрядного тока и напряжения в диапазоне давлений гелия 1−10 Тор. Установлено, что в исследуемых режимах ВАХ прибора имеет участок отрицательной дифференциальной проводимости, которая является причиной неустойчивости. На основании данных зондовых измерений рассчитана функция генерации ионов Г и величина D, определяющая их уход из плазменного канала. Показано, что развитие колебаний связано с нарушением условия Г ≈ D. Разработан метод подавления колебаний, заключающийся в управлении знаком дифференциальной проводимости плазмы путем отбора тока тепловых электронов на управляющий электрод. Предлагаемый метод обеспечивает высокий уровень стабильности рабочих параметров плазменного стабилизатора вплоть до плотности разрядного тока 5 А/см2.
Авторы выражают благодарность проф. Ю.А. Лебедеву за ценные замечания и постоянное внимание к работе.
Работа выполнена при поддержке гранта Российского научного фонда (РНФ) № 21-19-00139.
Список литературы
Ярыгин В.И., Ружников В.А., Синявский В.В. Космические и наземные ядерные энергетические установки прямого преобразования энергии. М.: НИЯУ “МИФИ”, 2016. 364 с.
Мустафаев А.С., Грабовский А.Ю. Низковольтный пучковый разряд в легких инертных газах для решения проблем стабилизации напряжения // ТВТ. 2017. Т. 55. № 1. С. 24.
Mustafaev A., Klimenkov B., Grabovskiy A., Kuznetsov V. Grid Current Control in the Unstable Mode of Plasma Discharge // J. Phys.: Conf. Ser. 2019. V. 1400(7). P. 077 024.
Мустафаев А.С. Динамика электронных пучков в плазме // ЖТФ. 2001. Т. 46. № 4. С. 111.
Sukhomlinov V.S., Matveev R., Mustafaev A.S., Timofeev N.A. Kinetic Theory of Low-voltage Beam Discharge Instability in Rare Gases // Phys. Plasmas. 2020. V. 27. № 6. P. 062 106.
Mustafaev A.S. Probe Method for Investigation of Anisotropic EVDF. Electron Kinetics and Applications of Glow Discharges // NATO Int. Sci. Session / Ed. By Kortshagen U., Tsendin L. N.Y.–London: Plenum Press, 1998. V. 367. P. 531.
Мустафаев А.С., Страхова А.А. 3D-диагностика функции распределения электронов в плазме // ЗГИ. 2017. Т. 226. С. 462.
Мустафаев А.С., Сухомлинов В.С. Функция распределения ионов в плазме собственного газа // Зап. Горн. ин-та. 2016. Т. 222. С. 864.
Вольтерра В. Теория функционалов, интегральных и интегродифференциальных уравнений. M.: Наука, 1982. 304 с.
Sukhomlinov V.S., Mustafaev A.S., Timofeev N.A. Probe Measurements of the Electron Velocity Distribution Function in Beams: Low-voltage Beam Discharge in Helium // J. Appl. Phys. 2018. V. 123. P. 143301.
Мустафаев А.С., Грабовский А.Ю. Зондовая диагностика анизотропной функции распределения электронов в плазме // ТВТ. 2012. Т. 50. № 6. С. 841.
Ахиезер A.И., Файнбìерг Я.Б. О взаимодействии пучка заряженных частиц с электронной плазмой // ДАН СССР. 1949. Т. 69. № 4. С. 555.
Харченко И.Ф., Файнберг Я.Б., Корнилов E.A. Взаимодействие пучка электронов с плазмой в магнитном поле // ЖТФ. 1961. Т. 31. С. 761.
Корнилов E.A., Ковпик O.Ф., Файнберг Я.Б. Механизм образования плазмы при развитии пучковой неустойчивости // ЖТФ. 1965. Т. 35. С. 1378.
Гаряинов С.Л., Абезгауз И.Д. Полупроводниковые приборы с отрицательным сопротивлением. М.: Энергия, 1970. 320 с.
Захарова В.М., Каган Ю.М. О движении ионов и атомов в плазме // Спектроскопия газоразрядной плазмы. Л., 1970. С. 291.
Мустафаев А.С., Грабовский А.Ю. Способ стабилизации электрических параметров в газоразрядных приборах с отрицательным сопротивлением. Патент РФ № 2498441. Б.И. № 31. 2013.
Дополнительные материалы отсутствуют.
Инструменты
Теплофизика высоких температур