CC BY
33
Новое в прикладной физике
Изв. вузов «ПНД», т. 20, № 6, 2012 УДК 621.385.69
ЗОНЫ УСТОЙЧИВОЙ ОДНОМОДОВОЙ ГЕНЕРАЦИИ В ГИРОТРОНЕ СО СВЕРХРАЗМЕРНЫМ РЕЗОНАТОРОМ
В. Л. Бакунин, Г. Г. Денисов, Н.А. Завольский, М.А. Моисеев
Путём численного исследования модели многомодового гиротрона с фиксированной структурой ВЧ-поля, на плоскости параметров «магнитное поле - ток» анализируются зоны устойчивой одномодовой генерации на частоте рабочей моды. Показано, что форма и размеры зон существенно зависят от плотности спектра мод резонатора и числа учитываемых мод-сателлитов (паразитных мод). С увеличением плотности спектра уменьшается площадь зоны, точка наибольшего значения КПД устойчивой одномодовой генерации сдвигается в область малых токов при одновременном уменьшении этого значения. При больших токах устойчивая одномодовая генерация достигается в мягком режиме возбуждения рабочей моды, а в жёстком режиме возникает примесь паразитных мод с мощностями, много меньшими мощности рабочей моды.
Ключевые слова: Гиротрон, устойчивая одномодовая генерация, плотность спектра собственных мод резонатора, фазовая связь мод, паразитные моды.
Введение
Для нагрева плазмы в существующих и проектируемых установках управляемого термоядерного синтеза требуются непрерывные гиротроны со всё большей мощностью и частотой излучения, которые в настоящее время достигают значений 1.5...2.0 МВт и 140...170 ГГц, соответственно. Для обеспечения приемлемого уровня омической тепловой нагрузки резонатора (порядка 2 кВт/см2) приходится увеличивать радиус резонатора и рабочие моды до ТЕ28д2...ТЕ34д9 и более высоких [1-4].
При увеличении индекса рабочей моды сгущается спектр собственных частот мод, попадающих в полосу циклотронного резонанса, что повышает опасность возбуждения паразитных мод и может ограничить возможности увеличения мощности и частоты гиротронов [4-7]. Представляет интерес определение условий существования устойчивой одномодовой генерации на частоте рабочей моды в зависимости от плотности спектра мод резонатора.
Для наиболее точного анализа взаимодействия мод гиротрона необходимо использовать уравнения с нефиксированной продольной структурой высокочастотного поля в резонаторе [8-10]. Однако в ряде случаев при достаточно высокой добротности резонатора используется приближение фиксированной структуры [5, 11-13]. Подобные модели качественно правильно описывают многие практически важные вопросы: устойчивость рабочей моды, конкуренцию и кооперацию мод.
В настоящей работе в пространстве параметров гиротрона с фиксированной структурой ВЧ-поля путём численных расчётов определяются области (зоны), соответствующие устойчивой одномодовой генерации (УОГ) ВЧ-излучения на частоте рабочей моды. Также определяется тенденция изменения оптимальных параметров, соответствующих максимально возможному значению КПД при УОГ, и выходных характеристик гиротрона при увеличении поперечных индексов рабочей моды и сгущении спектра собственных частот резонатора в полосе циклотронного резонанса.
1. Приближение фиксированной продольной структуры ВЧ-поля.
Понятие устойчивой одномодовой генерации. Уравнения гиротрона.
Эквидистантные мультиплеты. Сравнение моделей с фиксированной и нефиксированной продольной структурой ВЧ-поля
В случаях, когда диаметр резонатора гиротрона значительно превышает длину волны, в полосу циклотронного резонанса попадает большое число мод. Структура ВЧ-поля в резонаторе становится нефиксированной как по продольной координате, так и по азимутальной. Теория такого гиротрона была развита в [6, 8, 9].
При умеренной плотности спектра мод, когда их резонансные кривые не перекрываются, то есть — > ыs/Qs, где - частота, Qs - добротность в-й моды, можно использовать разложение поля по дискретному спектру мод (модовое представление). Поперечная структура поля каждой из мод считается фиксированной и совпадающей со структурой поля ТБтр-волны в пустом волноводе.
Обычно добротность резонаторов гиротронов достаточно велика - больше или сравнима с минимальной дифракционной добротностью [7, 14]
где Ь - длина резонатора, X - длина волны, соответствующая критической частоте в-й моды резонатора. В данных условиях продольная структура ВЧ-поля близка к фиксированной, соответствующей собственной моде «холодного» резонатора (в отсутствие электронного пучка). Приближение полностью фиксированной структуры ВЧ-поля позволяет получать некоторые практически важные результаты. Ниже мы ограничимся исследованием именно такой модели гиротрона.
Понятие «устойчивая одномодовая генерация», используемое нами здесь и далее, подразумевает следующее. Во-первых, расстояние между собственными частотами рассматриваемых мод резонатора превышает ширину их резонансных кривых (см. начало данного раздела), благодаря чему частоты мод являются различными и различимыми, что исключает возможность существования многомодовой одноча-стотной генерации. Во-вторых, предполагается, что режимы автомодуляционной и
(1)
стохастической генерации рабочей моды, соответствующие одномодовой многочастотной генерации и не являющиеся устойчивыми в соответствии со сложившейся терминологией, отсутствуют. Одномодовая многочастотная генерация наблюдалась при расчётах более сложных моделей гиротрона в [8] (низкодобротный резонатор, нефиксированная продольная структура ВЧ-поля) и в [15] (высокодобротный резонатор, наличие отражения от окна или нагрузки). В нашем случае система является высокодобротной без отражения, и данный тип генерации в расчётах не наблюдался. Эти ограничения делают понятие УОГ однозначно определённым.
Уравнения гиротрона с фиксированной продольной структурой ВЧ-поля, работающего на первой гармонике циклотронной частоты, как показано в [11-13], имеют следующий вид:
-Т- + Fs = и^ ¡: (д) ((V), ехр(-г^ о
^ + г (А + \р\2 - 1) р = г ^ Fsfs (д) ехр (г^), (2)
Fs (т = 0) = Г(0\
р (д = 0) = ехр (г,), 0 < ,< 2л.
Здесь р - безразмерный комплексный поперечный импульс электрона, Fs(т) -безразмерная комплексная амплитуда в-й моды, зависящие только от безразмерного времени
т = 200'
где £ - размерное время, ю0 - частота рабочей моды, являющаяся несущей частотой, 00 - добротность рабочей моды, с которой в рамках рассматриваемой модели можно отождествить добротности всех учитываемых мод Qs; , - фаза электрона при влёте в резонатор; ^ - азимутальная координата трубок тока электронного пучка. Угловые скобки означают усреднение по этим переменным. Продольное распределение в-й моды ВЧ-поля представлено в виде
fs (д) = fsо (д) ехр ^од), где fs0 (д) - продольное распределение ВЧ-поля в «холодном» резонаторе;
в± ®0<г
д
2в||
- безразмерная продольная координата, г - размерная продольная координата, |3ц;х -компоненты скорости электрона по отношению к скорости света с; д = ди - безразмерная длина пространства взаимодействия электронов с ВЧ-полем.
Считаем, что продольная структура ВЧ-поля описывается гауссовой функцией, что соответствует реальной структуре полей открытых резонаторов гиротро-нов [7, 17, 18]
(д) = ехр
-К! - о2
0 < д < ди, (3)
с
тогда - ширина распределения поля по уровню е 3. В качестве безразмерной длины резонатора в этом случае принимается значение ширины распределения ВЧ-поля по уровню е 1:
_ _ Ц _2|3ц с "
Далее,
2 — Юо
Д«0 _ --
в± Юо
- безразмерная отстройка частоты т., в-й моды от частоты рабочей моды Юо;
2 Юо — Юн
А
в± Юн
- безразмерная расстройка частоты рабочей моды Юо и циклотронной частоты Юн, определяемая величиной направляющего магнитного поля;
I. _ 2.35 • 10-41(А)
Чк
|ЗЦЗ||/ !/,о (д)!2 йд
о
- безразмерный токовый параметр, где I(А) - ток электронного пучка в амперах,
I I .¡^ 1 I V--Ки /К.. I \ 2
с _
стЕ —
1 I .т— 1 (ЧтрЯь/Ят )'
^ Ч2тр — т2\ .т (Утр) )
- фактор связи электронного пучка с модой ТБт,р, где Кт, Кь - радиусы резонатора и электронного пучка, соответственно, утр - р-й нуль производной функции Бесселя
.т (х);
_ (ю, — Юо) ¿о — (т. — то) ^
- добавка к фазе в-й моды, обусловленная различием времён ¿о вхождения групп электронов в пространство взаимодействия с ВЧ-полем и различием азимутальных структур мод. Для дублета мод время ¿о устраняется из первого уравнения (2) тривиально - путём переопределения азимутальной координаты Аналогично исключается ¿о в случае эквидистантного мультиплета мод с одинаковыми радиальными индексами, азимутальные индексы и частоты которых удовлетворяют соотношениям фазовой связи [5, 6, 11, 13, 18]:
т3 _ то + (—1)' ю, « юо + (—1).
в + 1 2
в + 1] 2
Ат,
АЮ, (4)
Юо
|Ю2Д:-1 + Ю2к — 2Юо| <7^,
Оо
где квадратные скобки означают взятие целой части.
Гиротрон, как правило, имеет два типа взаимодействия мод - конкуренцию и кооперацию [1]. Во взаимодействии первого типа участвуют моды с различными азимутальными индексами и неэквидистантным частотным спектром, не имеющие
фазовой связи (4). Характер связи мод в этом случае - энергетический, а генерация устанавливается на частоте одной моды - имевшей в начальный момент наибольший инкремент. Эквидистантные мультиплеты мод, удовлетворяющих условию (4), участвуют во втором типе взаимодействия. Самовозбуждение мультиплетов наиболее опасно, поскольку за счёт фазовой связи такие моды имеют наименьшие стартовые токи [6, 7, 13, 18-20]. При нахождении зон УОГ мы ограничимся рассмотрением именно таких мультиплетов.
Покажем, что модель с фиксированной структурой ВЧ-поля, описываемой гауссовой функцией (3), позволяет получать качественно правильные результаты для значений орбитального КПД и для точек перехода УОГ в мно-гомодовую генерацию. На рис. 1 приведены зависимости орбитального КПД от тока электронного пучка при оптимальной настройке магнитного поля (то есть при оптимальном значении А = = Aopt) для моделей с фиксированной и нефиксированной продольной структурой ВЧ-поля. Как следует из рис. 1, результаты расчётов показывают качественное соответствие этих двух моделей. Максимум орбитального КПД в модели с фиксированной структурой составляет 70%, что соответствует результату [7, 17] для ^ = 15. В модели с нефиксированной продольной структурой ВЧ-поля максимум орбитального КПД составляет 85%, что хорошо соответствует КПД гиротронов с оптимизированной структурой ВЧ-поля [3, 21]. Максимум КПД в обеих моделях достигается при близких токах пучка: I ~ 140 A. Переход от УОГ к многомодовой генерации с возбуждением эквидистантных мультиплетов в модели с фиксированной структурой ВЧ-поля происходит при токе I ~ 90 A, в модели с нефиксированной продольной структурой - при близком значении I ~ 70...80 A.
2. Влияние числа учитываемых мод на устойчивость рабочей моды
В вопросе исследования условий существования УОГ рабочей моды большое значение имеет число учитываемых паразитных мод.
На рис. 2 показана зона генерации стационарного одномодового приближения (в уравнениях (2) учитывается лишь одна мода) для фиксированной гауссовой
20
0 | | | | | | ' | ' | ' | ' |' | | | | 0 40 80 120 160/, А
Рис. 1. Зависимость орбитального КПД п±,орь оптимизированного по магнитному полю, от величины тока электронного пучка для моделей гиротро-на с фиксированной (кривая 1) и нефиксированной (кривая 2) структурами ВЧ-поля. Рабочая мода - ТЕ28 д2 с частотой 170 ГГц, ускоряющее напряжение V = 100 кВ, питч-фактор д = 1.2, безразмерная длина резонатора = 15. Точки с цифрами на кривой, соответствующей фиксированной структуре, обозначают моменты перехода устойчивой одномодовой генерации в многомодовую (указано число учитываемых мод: 2 моды - дублет ТЕ27д2 - ТЕ28д2, 3 моды - эквидистантный триплет ТЕ27Д2 - ТЕ28Д2 - ТЕ29Д2, 5 мод - эквидистантный пентаплет ТЕ26,12 - ТЕ27д2 - ТЕ28д2 -ТЕ29д2 - ТЕзод2). Пунктиром показаны значения П±,орл в одномодовом приближении. Точка «макс» соответствует максимуму орбитального КПД в од-номодовом приближении с фиксированной структурой ВЧ-поля, составляющему 70%
Рис. 2. Зона генерации (зависимость орбитального КПД одномодового приближения г|± от расстройки А). Безразмерная длина резонатора и токовый параметр пучка близки к оптимальным по КПД: = 15, I = 0.6. Рабочая мода ТЕ28д2 имеет расстройку, оптимальную по КПД: А = 0.53. Показаны положения мод эквидистантного мультиплета, находящихся в фазовой связи с рабочей модой и их безразмерной отстройки относительно друг друга
продольной структуры ВЧ-поля (3). Верхняя ветвь кривой п±(А) (сплошная линия) является устойчивой, нижняя (пунктир) - неустойчива. Рабочая мода ТЕ28,12, имеющая оптимальное значение расстройки А _ 0.53, при котором достигается максимум орбитального КПД п± _ 70% [7, 17], находится в режиме жёсткого возбуждения колебаний. Левые паразитные моды ТЕ27Д2 и ТЕ26Д2 также находятся в пределах зоны генерации и могут самовозбуждаться. Правая паразитная мода ТЕ2д12 находится вне зоны генерации, однако её возбуждение возможно за счёт фазовой связи мод [5, 6, 11, 13, 18].
Попарное взаимодействие рабочей моды ТЕ28,12 с ТЕ27.12 или с ТЕ26,12
является несинхронным и не приводит к нарушению устойчивости ее генерации.
Учёт паразитной моды ТЕ29Д2, находящейся вне зоны генерации, приводит к уменьшению стартового тока паразитной генерации вследствие фазовой связи мод (4) в триплете ТЕ2712 - ТЕ2812 - ТЕ2д12 [9]. При учёте только трёх данных мод возбуждение триплета срывает генерацию на рабочей моде (рис. 3, а). Вначале амплитуды паразитных мод остаются на уровне шумов, а амплитуда рабочей близка к стационарному значению одномодового приближения, орбитальный КПД _ 63%. Нарастание амплитуд паразитных мод приводит к одновременному срыву рабочей моды и установлению генерации на левой паразитной моде ТЕ27Д2 с орбитальным КПД _ 27%. Амплитуда правой паразитной моды ТЕ29Д2, как находящейся вне зоны генерации (см. рис. 2), в итоге снова спадает до нуля.
Рис. 3. Срыв генерации с рабочей моды ТЕ28д2 на паразитную ТЕ27,12 в триплете ТЕ27,12 - ТЕ28д2 -ТЕ29Д2 (а). Установление многомодовой генерации с наибольшей амплитудой у рабочей моды в квадруплете ТЕ26д2 - ТЕ27д2 - ТЕ28д2 - ТЕ29д2 (б). Значения параметров: = 15, А = 0.5, I = 0.06
Учёт ещё одной левой паразитной моды ТЕ26,12 стабилизирует рабочую моду (рис. 3, б), амплитуда которой лишь уменьшается на 15% по сравнению со стационарным одномодовым значением. Устанавливается 4-модовая генерация с орбитальным КПД п± = 56%, который незначительно меньше стационарного одно-модового. Аналогичный результат получается при добавлении симметричной моде ТЕ26,12 моды ТЕэ0,12, расположенной далеко правее зоны генерации: устанавливается 5-модовая генерация с наибольшей амплитудой у рабочей моды и тем же орбитальным КПД п± = 56%.
Полученные результаты показывают, что необходимо, как минимум, учитывать все моды, находящиеся в пределах зоны генерации одномодового приближения плюс одну высокочастотную моду вне этой зоны. Влиянием остальных паразитных мод можно пренебречь. Поскольку ширина зоны генерации одномодового приближения (см. рис. 3) зависит от величины токового параметра, при одновременном изменении параметров А и I число мод, попадающих в эту зону, может изменяться. Вследствие этого при исследовании вопросов существования УОГ рабочей моды в пространстве параметров А, I следует учитывать уже все паразитные моды, находящиеся в полосе циклотронного резонанса и имеющие фазовую связь (4) с рабочей модой.
3. Зоны устойчивой одномодовой одночастотной генерации в пространстве параметров «магнитное поле - ток»
Выводы, сделанные в предыдущих разделах, позволяют ограничиться рассмотрением динамической системы - модели гиротрона с гауссовой фиксированной продольной структурой ВЧ-поля (3), описываемой всего 5 безразмерными параметрами.
• I - токовый параметр рабочей моды. В настоящем исследовании считаем токовые параметры всех паразитных мод также равными I, по сути пренебрегая различием Ст>р-факторов связи мод с электронным пучком. Тем самым мы создаём более благоприятные условия для возбуждения паразитных мод (в реальности величина Ст,р тем меньше, чем дальше отстоит паразитная мода, имеющая фазовую связь с рабочей модой, от последней). Условия для существования УОГ рабочей моды, напротив, становятся более жёсткими, и будут иметь характер «оценки сверху».
• А - параметр, характеризующий величину магнитного поля.
• ц - безразмерная длина резонатора, которую мы зафиксируем на значении 15, близком к оптимальному.
-2 Ю1 - Ю0 „ г- г
• А = —2---отстройка частоты ближайшего сателлита от частоты рабо-
Р]_ «0
чей моды. Поскольку мы ограничиваемся рассмотрением мод, удовлетворяющих (4), нетрудно видеть, что параметр А- является характеристикой плотности спектра частот эквидистантных мод. Чем меньше А-, тем плотнее спектр, и обратно.
• N - число учитываемых мод.
Зафиксировав параметры ц, А, N на плоскости параметров «магнитное поле-ток» А, I, будем искать области (зоны), соответствующие УОГ. Обратим внимание,
что в этих зонах существует генерация только на рабочей моде, совпадающая со стационарным режимом одномодового приближения при тех же параметрах, а амплитуды всех паразитных мод равны нулю.
В начальный момент времени рабочая мода возбуждена до амплитуды, превышающей стационарное одномодовое значение, а паразитные моды имеют амплитуды на 2 порядка меньше - на уровне шумов. Проблемы сценария приведения гиротро-на в данное состояние мы не касаемся, это требует дополнительного исследования [22-24], которое будет произведено нами в последующем.
Ниже приведены результаты численных расчётов зон УОГ и их анализа по мере перехода от низших (с разреженным спектром) ко всё более высоким рабочим модам (с плотным спектром). При этом мы удваиваем поперечные индексы рабочих мод, сохраняя соотношение т/р ~ 2, поскольку в гиротронах мегаваттной мощности, как правило, Щ ~ /2 [2-4]. Чем плотнее спектр, тем большее число мод, имеющих фазовую связь (4) с рабочей модой, которые необходимо учитывать, попадает в полосу циклотронного резонанса.
3.1. Мода ТЕ12,5 (плотность спектра Д = 0.48). Спектр собственных частот резонатора вблизи рабочей моды ТЕ^^ является разреженным: Д = 0.48 (рис. 4, б). Зона УОГ для данной рабочей моды при учете 5 мод в мультиплете показана на рис. 4, а. Линия стартового тока 1старт [17] замыкает зону при малых значениях токового параметра I, как и должно быть. Область внутри кривой стартового тока соответствует мягкому, вне - жёсткому режиму возбуждения рабочей моды. Максимальное значение орбитального КПД УОГ (отмечено точкой) равно максимальному орбитальному КПД одномодового приближения: = 70% и достигается внутри зоны в жёстком режиме при значении токового параметра, равном оптимальному в одномодовом приближении: I = 0.06. Часть зоны, расположенная в области больших токов, соответствует второму и следующим максимумам КПД (имеющим меньшее значение, чем первый) и представляют интерес с точки зрения теории гиротронов [3, 16, 20].
8со/ю, %
Рис. 4. Зона УОГ при = 15, плотности спектра Д = 0.48 (рабочая мода ТЕ12,б), N = 5 (эквидистантный пентаплет ТЕ1о,б - ТЕи,б - ТЕ12,5 - ТЕ13,5 - ТЕм^) (а). Зависимость Ст,р-фактора от частоты в окрестности рабочей моды ТЕ12,5 (разреженный спектр) (б)
3.2. Мода ТЕ28,12 (плотность спектра Д = 0.2). В случае рабочей моды ТЕ28.12 спектр собственных частот резонатора является умеренно плотным: Д = 0.2 (рис. 5, в). Для данной рабочей моды были рассмотрены два значения числа учитываемых мод: N = 3 и N = 5 (рис. 5, а). Переход от триплета к пентаплету приводит к распаду большой зоны на две - «нижнюю», расположенную в области сравнительно небольших значений токового параметра (показана отдельно на рис. 5, б) и «верхнюю», соответствующую большим значениям I. Обе меньших зоны расположены внутри большей, касаясь её, что обеспечивает одинаковые максимальные значения орбитального КПД УОГ. Это значение достигается на границе «нижней» зоны (рис. 5, б) и составляет = 65%, что несколько меньше максимума в одномодовом приближении. Рис. 1 является иллюстрацией этого явления: после того, как кривая Д = Допт (I), доставляющая значения оптимизированного по магнитному полю орбитального КПД, пересечёт границу зоны УОГ, необходимо оставаться на правой границе зоны, чтобы по-прежнему обеспечивать «чистую» генерацию рабочей моды, что даёт меньшие значения по сравнению с одномо-довым приближением. Значение орбитального КПД УОГ при токе I = 0.06 (являющемся оптимальным в одномо-довом приближении) составляет лишь = 53%, что заметно меньше 70%.
Заметим, что полученный результат согласуется с результатом [18], где исследовалась устойчивость генерации рабочей моды относительно возбуждения триплета. При Д = 0.2 генерация рабочей моды с высоким орбитальным КПД устойчива (I и 0.04, ц± < 65%), при этом она находится в жёстком режиме возбуждения колебаний. В области малых расстроек (I и 0.04, Д < 0.25) режим возбуждения колебаний - мяг-
Рис. 5. Сравнение зон УОГ при = 15, плотности спектра Д = 0.2 (рабочая мода ТЕ28,12), N = 3 (эквидистантный триплет ТЕ27,12 - ТЕ28,12 - ТЕ29,12, одинарная штриховка) и N = 5 (эквидистантный пентаплет ТЕ26,12 - ТЕ27Д2 - ТЕ28Д2 - ТЕ29Д2 -ТЕзо,12, двойная штриховка) (а). Первая зона УОГ для N = 5 крупным планом (б). Зависимость От,Р-фактора от частоты в окрестности рабочей моды ТЕ28,12 (умеренно плотный спектр) (в)
кий, однако попасть отсюда в область высоких КПД мешает самовозбуждение паразитных мод, поскольку область малых расстроек находится за пределами зоны УОГ.
3.3. Мода ТЕ56,24 (плотность спектра Д = 0.1). В случае рабочей моды ТЕ56.24 спектр собственных частот является плотным: Д = 0.1 (рис. 6, б). Мульти-плет состоит из 11 мод, попадающих в полосу циклотронного резонанса. Зона УОГ приведена на рис. 6, а. Максимальное значение орбитального КПД = 58% достигается на границе зоны. При токе I = 0.06 максимальное значение = 39%, причём оно достигается уже не в жёстком режиме, как в случае рабочей моды ТЕ28Д2, а в мягком режиме.
3.4. Мода ТЕц2 48 (плотность спектра Д = 0.05). Для высокой рабочей моды ТЕц2,48 спектр собственных частот является очень плотным: Д = 0.05 (рис. 7, б). В полосу циклотронного резонанса попадают 19 мод, находящихся в фазовой связи с рабочей модой. Зона УОГ приведена на рис. 7, а. Максимальное
0.06-
0.04
0.02-
0-
О.
г|± = 39%
т,р
о.оооз -
0.0002
56,24
=58% 0
0001-
61,24
->-1-1-1-1-1-1-1
0 2.0 4.0 6.0 Д -4.0 -2.0 0 2.0 4.0
а б 5ю/со, %
Рис. 6. Зона УОГ при = 15, плотности спектра Д = 0.1 (рабочая мода ТЕбб,24), N =11 (эквидистантный мультиплет ТЕ51>24 - ... - ТЕ56,24 - ... - ТЕб1>24) (а). Зависимость Ст>р-фактора от частоты в окрестности рабочей моды ТЕб6.24 (плотный спектр) (б)
I-
0.06 -
0.04 -
0.02 -
0
О.
1^= 35%
т,р
0.00008
0.00004 -
г|1= 52%
121,48
5ю/со,
0 0.2 0.4 0.6 А -4.0 -2.0
а б
Рис. 7. Зона УОГ при = 15, плотности спектра Д = 0.05 (рабочая мода ТЕ112,48), N = 19 (эквидистантный мультиплет ТЕ1оз,48 - ... - ТЕ112,48 - ... - ТЕ121,48) (а). Зависимость Ст,р-фактора от частоты в окрестности рабочей моды ТЕ112,48 (очень плотный спектр) (б)
значение орбитального КПД = 52%достигается также на границе зоны. При токе I = 0.06 наибольшее значение орбитального КПД = 35% достигается также в мягком режиме.
3.5. Выводы из полученных результатов. Основные результаты, полученные в настоящей работе, проиллюстрируем рис. 8. Здесь друг на друга наложены «нижняя» зона УОГ, полученная для рабочей моды ТЕ28Д2 и мультиплета из 5 мод, зона УОГ для рабочей моды ТЕ56,24 и мультиплета из 11 мод и зона УОГ для рабочей моды ТЕ112>48 и мультиплета из 19 мод. Рис. 8 вкупе с рис. 5, б, рис. 6, а и рис. 7, а иллюстрирует общую тенденцию - с увеличением поперечных индексов рабочей моды и сгущением спектра собственных частот резонатора:
• уменьшается площадь зоны УОГ;
• уменьшается максимальное значение орбитального КПД УОГ, соответствующий этому значению токовый параметр также уменьшается;
• при значении токового параметра, являющемся оптимальным в одномодовом приближении, генерация рабочей моды смещается из жёсткого режима в мягкий.
Также обращают на себя внимание ещё два интересных факта:
• несмотря на быстрое увеличение (удвоение) поперечных индексов рабочей моды и уплотнение спектра площадь зоны уменьшается достаточно медленно -это говорит о том, что при обеспечении сценария включения, приводящего к начальному моменту времени к возбуждению рабочей моды и удерживанию паразитных мод на уровне шумов, возможна УОГ даже на очень высокой моде в области очень плотного спектра собственных частот (ТЕ112Д8) с достаточно высоким КПД (п± = 52%);
• вследствие сложной формы границ зон при фиксированной расстройке Д возможно чередование зон УОГ рабочей моды с областями многомодовой генерации (и/или с областями неустойчивости рабочей моды).
Необходимо отметить, что полученный нами результат близок к результату [23], где также исследуются зоны УОГ в пространстве параметров гиро-трона. Однако в [23] в качестве рабочей моды принимаются высокие моды «шепчущей галереи», такие как ТБ80,4 - с менее развитой пространственной структурой и более разреженным спектром - Д варьируется от 0.97 до 0.13. Границы зон УОГ в [23] делятся на два типа: соответствующие возбуждению двух ближайших к рабочей моде паразитных мод и многомодовой генерации с участием большего числа паразитных мод. В нашем случае рабочие моды ТЕ12,5, ТЕ28,12, ТЕ56,24, ТЕц2,48 имеют более развитую пространственную структуру, а начиная с моды ТЕ56,24
0 0.2 0.4 0.6 А
Рис. 8. Сравнение зон УОГ при =15 для Д = 0.2 (рабочая мода ТЕ28,12, эквидистантный мультиплет из 5 мод, одинарная штриховка), Д = 0.1 (рабочая мода ТЕБ6,24, эквидистантный мультиплет из 11 мод, двойная штриховка), Д = 0.05 (рабочая мода ТЕ112,48, эквидистантный мультиплет из 19 мод, тройная штриховка)
спектр собственных частот резонатора более плотный, чем в [23]. Кроме того, при определении границ зон УОГ мы не производили их разделение по типам неустойчивости, и рассматривали все имеющие фазовую связь с рабочей модой паразитные моды, попадающие в полосу циклотронного резонанса, не ограничивая их количество. При этом полученный нами вывод о медленном убывании площади зоны УОГ при сгущении спектра согласуется с аналогичным результатом [23].
Влияние густоты спектра на устойчивость рабочей моды исследовалось также в работах [25, 26]. Результат, полученный в [26] для модели лазера на свободных электронах, состоит в том, что при токе электронного пучка, не превышающем или равном оптимальному по КПД, устанавливается одномодовая генерация - на моде, имеющей наибольший инкремент, причем это утверждение справедливо при сколь угодно плотном эквидистантном спектре мод зеркального резонатора. При увеличении тока пучка в [25, 26] рабочая мода становится неустойчивой относительно возбуждения ближайших сателлитов, при больших токах возникает режим автомодуляционной, а затем стохастической генерации. В отличие от модели [25, 26], модель гиротрона с фиксированной структурой ВЧ-поля, рассмотренная нами, описывается иными уравнениями, благодаря чему на высоких рабочих модах точка оптимального по КПД набора параметров оказывается вне зоны УОГ. Отметим, что в рамках нашей модели при определённом наборе параметров, в частности Д = 0.2 (рабочая мода ТЕ28Д2), И- = 15, I = 0.14, Д = 0.5, возможна многомодовая автомодуляционная генерация. Стохастическая же генерация даже при токах, значительно превышающих стартовый, при всех рассмотренных в данной работе значениях параметров гиротрона не наблюдалась (аналогичное этому заключение сделано в [5]).
В настоящей работе не рассматривались моды, не имеющие фазовой связи (4) с рабочей модой. Таковыми являются, в частности, моды обратного вращения и моды с радиальными индексами, отличными от такового у рабочей моды. Влияние этих мод на форму зон УОГ будет рассмотрено в дальнейшем. Кроме того, в спектре собственных колебаний резонатора присутствует большое количество мод с малыми Ст>р-факторами связи с электронным пучком, в 2 и более раз меньшими От р-фактора рабочей моды. Такие моды опасности для УОГ рабочей моды не представляют. Более подробно это будет рассмотрено в дальнейшем.
Заключение
В настоящей работе путём численного исследования модели одномодового ги-ротрона с фиксированной структурой ВЧ-поля получены и проанализированы зоны устойчивой одномодовой генерации на частоте рабочей моды в пространстве рабочих параметров «ток электронного пучка - магнитное поле». Как следует из полученных результатов, размеры зон зависят от густоты спектра паразитных мод, находящихся в полосе циклотронного резонанса и имеющих фазовую связь с рабочей модой, однако площадь зоны с увеличением рабочей моды уменьшается достаточно медленно. Однако наибольшее достижимое значение КПД устойчивой одномодовой генерации уменьшается, а соответствующая ему точка в пространстве параметров смещается в сторону уменьшения тока. При больших токах и высокой рабочей моде
устойчивая одномодовая генерация возможна только при значениях рабочих параметров, соответствующих мягкому режиму возбуждения рабочей моды.
Авторы выражают благодарность В.Е. Запевалову и М.Ю. Глявину за обсуждение работы и ценные критические замечания.
Библиографический список
1. Litvak А^., Denisov G.G., Myasnikov V.E., Tai E.M., Azizov E.A., Ilin V.I. Development in russia of megawatt power gyrotrons // J. Infrared Milli Teraherz Waves. 2011. Vol. 332. P. 337.
2. Аликаев В.В., Денисов Г.Г., Запевалов В.Е., Курбатов В.И., Литвак А.Г., Мясников В.Е., Тай Е.М. Гиротроны для УТС // Вакуумная СВЧ электроника: Сборник обзоров. Нижний Новгород: ИПФ РАН, 2002. С. 71.
3. Методы нелинейной динамики и теории хаоса в задачах электроники сверхвысоких частот. Т.1. Стационарные процессы // Под редакцией А.А. Кураева, Д.И. Трубецкова. М.: Физматлит, 2009.
4. Thumm M. State-of-the-art of high power gyro-devices and free electron masers. Karlsruhe: KIT, 2010.
5. Nusinovich G.S. Mode interaction in gyrotrons // Int. J. Electronics. 1981. Vol. 51, № 4. P. 457.
6. Петелин М.И. Электронная селекция мод в гиротроне // Гиротрон. Сборник научных трудов. Горький: ИПФ АН СССР, 1981. С. 77.
7. Nusinovich G.S. Introduction To The Physics of Gyrotrons // The Johns Hopkins University Press, Baltimore, 2004.
8. Ginzburg N.S., Nusinovich G.S., Zavolsky N.A. Theory of non-stationary processes in gyrotrons with low Q resonators // Int. J. Electronics. 1986. Vol. 61, № 6. P. 881.
9. Завольский Н.А., Нусинович Г.С. Нестационарные процессы в гиротроне с нефиксированной структурой ВЧ поля // Радиотехника и электроника. 1991. Т. 36, № 1. С. 135.
10. Завольский Н.А., Нусинович Г.С., Павельев А.Б. Устойчивость одномодовых колебаний и нестационарные процессы в гиротронах со сверхразмерными низкодобротными резонаторами // Гиротрон. Сборник научных трудов. Горький. ИПФ АН СССР. 1989. С.84.
11. Моисеев М.А., Нусинович Г.С. К теории многомодовой генерации в гиромоно-троне // Изв. вузов. Радиофизика.1974. Т. XVII, № 11. С. 1711.
12. Бляхман Л.Г., Нусинович Г.С. Динамика многомодовых электронных мазеров // Радиотехника и электроника. 1982, № 5.
13. Запевалов В.Е., Нусинович Г.С. К теории амплитудно-фазового взаимодействия мод в электронных мазерах // Изв. вузов. Радиофизика. 1989. Т. 32, № 3. С. 347.
14. Власов С.Н., Загрядская Л.И., Петелин М.И. Резонаторы и волноводы с модами шепчущей галереи для мазеров на циклотронном резонансе // Изв. вузов. Радиофизика. 1973. Т. 16, № 11. C. 1743.
15. Dumbrajs O., Glyavin M.Yu., Zapevalov V.E., Zavolsky N.A. Influence of reflections on mode competition in gyrotrons // IEEE Transactions on Plasma Science. 2000. Vol. 28. P. 588.
16. Cai S.Y., Antonsen T.M., Saraph G.Jr., Levush B. // Int. J. Electron. 1992. Vol. 72, № 5-6. P. 759.
17. Нусинович Г.С., Эрм Р.Э. КПД МЦР-монотрона с гауссовым продольным распределением высокочастотного поля // Электронная техника. Серия 1. Электроника СВЧ. 1972, № 8. С. 55.
18. Глявин М.Ю., Нусинович Г.С. Устойчивость одномодовых автоколебаний в ги-ротроне с синхронным взаимодействием мод // Радиотехника и электроника, 1991.
19. Dumbrajs O., Anderer J., Illy S., Piosczyk B., Thumm M., Zavolsky N.A. Multi-frequency operation of a gyrotron // IEEE Transactions on Plasma Science. 1999. Vol. 27, № 2.
20. Carmel Y., Chu K.R., Read M.E., Klim K.J., Arfin B, Dialetis D, Fliflet A. Mode competition, suppression, and efficiency enhancement in overmoded gyrotron oscillators // Int. J. Infrared and MM Waves. 1982. Vol. 3. P. 645.
21. Кураев А.А., Ковалев И.С., Колосов С.В. Численные методы оптимизации в задачах электроники СВЧ. Минск: Наука и техника, 1975.
22. Нусинович Г.С. Способы подачи напряжений на импульсный гиромонотрон, обеспечивающие высокий КПД в режиме одномодовой генерации // Электронная техника. Серия 1. Электроника СВЧ. 1974. № 3. С. 44.
23. Levush B., Antonsen T.M. Mode competition and control in high-power gyrotron oscillators // IEEE Transactions on Plasma Science. 1990. Vol. 18, № 3. P. 260.
24. Grudiev A., Raguin J.-Y., Schunemann K. Numerical study of mode competition in coaxial cavity gyrotrons with corrugated insert // Int. J. of Infrared and Millimeter Waves. 2003. Vol. 24, No. 2. P. 173.
25. Bogomolov Ya.L., Bratman V.L., Ginzburg N.S., Petelin M.I., Yunakovsky A.D. Non-stationary generation in free electron lasers // Optics Communications. 1981. Vol. 36, № 3. P. 209.
26. Ginzburg N.S., Petelin M.I., Sergeev A.S. On the mechanism of self-modulatin onset in free electron lasers // Optics Communications. 1985. Vol. 55, № 4. P. 283.
ИПФ РАН, Н. Новгород Поступила в редакцию 24.04.2012
После доработки 21.06.2012
ZONES OF STABLE SINGLE-MODE GENERATION IN OVERMODED GYROTRONS
V.L. Bakunin, G. G. Denisov, N.A. Zavolsky, M.A. Moiseev
During numerical simulations of the model of gyrotron with fixed field-structure, zones of stable single-mode generation of the work mode are analyzed on the plane of parameters «magnetic field - beam current». It is shown that area and sizes of the zones are strongly dependent on the spectral density of eigen-frequencies of the resonator and on number of considered parasit modes. With increase of spectral density the point of highest attainable efficiency of stable single-mode generation shifts to the region of lower currents and simultaneously this value of efficiency decreases. In the region of higher
currents stable single-mode generation exists in soft self-excitation regime of the work mode, while in the region of hard self-excitation of the work mode multimode generation with power of parasit modes much less than power of the work mode arises.
Keywords: Gyrotron, stable single-mode generation, spectral density of the resonator, phase relationship of the modes, parasit modes.
Бакунин Владимир Лазаревич - родился в Горьком (1986). Окончил магистратуру физического факультета Нижегородского государственного университета им. Лобачевского (2009), с момента окончания по настоящее время работает в Институте прикладной физики РАН. Область научных интересов - численные методы, исследование электронно-волнового взаимодействия и динамических режимов генераторов СВЧ-излучения мегаваттного уровня мощности.
603950 Н. Новгород, ГСП-120, ул. Ульянова, 46 Институт прикладной физики РАН E-mail: knightrider7@yandex.ru
Денисов Григорий Геннадьевич - родился в Горьком (1956). Окончил радиофизический факультет Горьковского госуниверситета по специальности «Радиофизика и электроника» (1978). Д.ф.-м.н. (2002), профессор, чл.-корр. РАН, зав. Отделом электронных приборов ИПФ РАН. Автор более 200 научных работ. Профессиональные интересы - мощные источники электромагнитного излучения, многомодовые системы, методы измерения, преобразования структур полей квазиоптических волновых пучков, антенны и линии передачи мощного коротковолнового излучения
603950 Н. Новгород, ГСП-120, ул. Ульянова, 46 Институт прикладной физики РАН E-mail: den@appl.sci-nnov.ru
Завольский Николай Александрович - родился в 1956 году в Ворсме (Нижегородская область). Окончил Горьковский государственный университет по специальности математика (1978). После окончания и до настоящего момента работает в Институте прикладной физики Академии наук СССР (в настоящее время ИПФ РАН) Область научных интересов: численные методы и разработка программного обеспечения в решении задач акустики, распространения радиоволн и электронно-волнового взаимодействия в гиротронах. Автор более 80 научных работ.
603950 Н. Новгород, ГСП-120, ул. Ульянова, 46 Институт прикладной физики РАН E-mail: zavolsky@appl.sci-nnov.ru
Моисеев Марк Андреевич - родился в Горьком (1943). Окончил Горьковский государственный университет по специальности радиофизика (1966). После окончания работал в НИРФИ, а затем до настоящего момента в Институте прикладной физики Академии наук СССР (в настоящее время ИПФ РАН). Область научных интересов: теоретическое исследование электронно-волновых процессов в мощных гироприборах, распространения радиоволн в квазиоптических системах. Автор более 80 публикаций.
603950 Н. Новгород, ГСП-120, ул. Ульянова, 46 Институт прикладной физики РАН E-mail: moiseev@appl.sci-nnov.ru
