Особенности многоволнового взаимодействия электронного потока и электромагнитного поля в релятивистском дифракционном генераторе
А. И. Слепков1*2, О. В. Галлямова26
1 Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра общей физики. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2. 2МОУ Видновская СОШ № 7.
E-mail: 11 slepkov@phys.msu.ru, ьolgalglm@googlemail.com
Статья поступила 08.12.2008, подписана в печать 27.02.2009.
Исследованы механизмы самосогласованного взаимодействия трубчатого релятивистского электронного потока с полем сверхразмерной периодической структуры релятивистского дифракционного генератора в области частот 27г-вида колебаний основной моды. Численный анализ проводился с помощью матричного многомодового метода в его нестационарной вариации. Изучены структура полей, их модовый состав и спектр излучения. Выявлены условия стабилизации частоты генерации, связанные с продольными резонансами поверхностной пространственной гармоники и объемными резонансами на критической частоте моды Е§п цилиндрического волновода.
Ключевые слова: дифракционное излучение, релятивистский дифракционный генератор.
УДК: 621.385.69. PACS: 84.40.Fe.
Введение
Одним из решений задачи повышения импульсной мощности микроволнового излучения является использование релятивистских электронных потоков (РЭП) и электродинамической системы с поперечными размерами, значительно превышающими длину волны излучения. Наиболее перспективными в этом направлении оказались микроволновые источники, основанные на продольном взаимодействии трубчатого релятивистского электронного потока с полями сверхразмерных периодических волноводов — многоволновые черепковские генераторы (МВЧГ) и релятивистские дифракционные генераторы (РДГ) [1]. Последние, несколько уступая по мощности, обладают важным преимуществом: длительность импульса излучения Tgen в РДГ близка к длительности импульса тока rpuis, в то время как в МВЧГ наблюдается срыв генерации и ограничение длительности импульса Tgen ~ O.lTpuis [1]. В частности, в МВЧГ получена импульсная мощность излучения 15 ГВт в сантиметровом диапазоне длин волн [2], в РДГ — около 4.5 ГВт при Arad = 7.2 мм [3]. Для решения проблемы селекции мод и получения одночастотной генерации в этих устройствах было реализовано взаимодействие потока и поля на частотах соответственно ж- и 2тт-видов колебаний основной аксиально-симметричной моды периодического волновода. Механизм излучения в первом из указанных типов приборов основан на эффекте Вавилова-Черен-кова. Во втором типе устройств в основе физических процессов, определяющих их действие, лежит дифракционное излучение (ДИ) РЭП: электронный пучок, двигаясь со скоростью оо вблизи периодических неоднородностей с периодом d, излучает набор распространяющихся обобщенных плоских волн, характеризующихся направлением eos вп = + (или волновыми числами Ац„ = f + ^p,
= 0, ±1, ±2,..., /? = у и k±n = у¡k2 — /г|2), при переотражении которых от противоположной стенки в сверхразмерной структуре устанавливаются незатухающие колебания. Когда вп действителен, возникает излучение на п-й пространственной гармонике, а под ДИ понимают совокупность излучающихся пространственных гармоник. В случае структур бесконечной длины при \Щп\ < w/c
говорят о быстрых или объемных гармониках, а в случае \Щп\ > ш/с — о поверхностных (медленных) гармониках. Граница между ними — линии скорости света ш = ±к2с.
В структурах конечной длины используется понятие ДИ с учетом того, что на бриллюэновской диаграмме [4] точки аномалий Вуда, определяемые как пересечение соответствующих пространственных гармоник с линиями скорости света, вырождаются в размытые области резонансного отклика системы [5].
В настоящей работе особенности ДИ в РДГ исследованы с точки зрения традиционных для электроники подходов, основанных на анализе самосогласованного взаимодействия электронного потока и поля периодической структуры. При этом ДИ вблизи аномалий Вуда [4] может трактоваться как возбуждение собственных волн при условии синхронизма потока и поля.
1. Математическая модель и метод численного исследования
За основу теоретического метода анализа работы дифракционного генератора был выбран матричный много-модовый метод [6], использовавшийся ранее для анализа многоволнового черенковского генератора. Метод позволяет рассматривать процессы установления колебаний в изучаемых системах в слабонестационарном приближении, предполагающем малое изменение амплитуды поля за период Т = 2тт/ш.
В качестве электродинамической структуры рассматривается идеально проводящий цилиндрический волновод с периодическими неоднородностями в виде полуторов на пьедестале, й — период структуры, У?т1п и Лтах — внутренний и внешний радиусы волновода (рис. 1). Вдоль оси симметрии системы пропускается трубчатый электронный поток со скоростью электронов оо. определяемой ускоряющим напряжением £/, Яп — радиус потока, ток пучка /. Параметры структуры и РЭП в численном эксперименте подбирались близкими к характеристикам экспериментально реализованного в работе [3] РДГ. При исследовании возбуждения многих мод сверхразмерного периодического волновода аксиально-симметричный нерегулярный волновод (рис. 1) представляется в виде последовательности участков гладких волноводов
РАДИОФИЗИКА, ЭЛЕКТРОНИКА, АКУСТИКА
31
о| о| и <4
дг
■дд
25
Ох, дг
г) / Яг1
( хйЯ± ;
3±
25 дг
(2)
Здесь £>±,...,55* — матричные коэффициенты, определяемые из уравнений Максвелла, — энергетический интеграл взаимодействия потока и поля, имеющий в пренебрежении силами пространственного заряда следующий вид:
т+2тг
м.
VI
У —
т= 1
х (кп.3)*]с15 = [е„.5_1 хД1.,] (¡Б =
е„.5_1 х ¿Б.
(3)
Система уравнений (2), (3) решается на каждом временном шаге с помощью матричной прогонки. Для каж-
дой крупной частицы в пренебрежении силами пространственного заряда записывается релятивистское уравнение движения: ^ = ^ Ие(Еиге'г), где р = тщ'у —импульс
крупной частицы, 7 =
— релятивистский фактор,
Рис. 1. Модель электродинамической структуры РДГ
фиксированного диаметра. Процессы в электронном пучке описываются с помощью модели крупных частиц, представляющих собой бесконечно тонкие кольца с зарядом ц. Вихревое электрическое и магнитное поля на каждом выделенном участке с номером 5 разлагаются по системам собственных функций выделенных объемов:
N
Еил = У ' [С/г^^' + Сп 3(г, t)En_sj ,
п=1 (1
Л' ^
н Е + СпАг, ПН„-] ,
п=1
где C^:s(г,t) — комплексные коэффициенты разложения, а Е^3(г,г), Н^3(г,г) — векторы электрического и магнитного полей для моды £ол_™па (с учетом аксиальной симметрии структуры). После подстановки разложений (1) в уравнения Максвелла уравнения возбуждения записываются в матричном виде относительно векторов л:1*1 = (С^.С*,... , . В слабонестационарном приближении они имеют вид
М$ — число частиц, пролетающих через 5-й регулярный участок в интервал времени (т, т+2тт). Уравнения возбуждения дополняются условиями сшивания полей на скачках радиуса гладкого волновода, записанными в виде проекционных соотношений, обеспечивающих непрерывность потока вектора Умова-Пойнтинга:
т = Ш — нормированное время, Еиг — продольная составляющая напряженности вихревого поля. В слабонестационарном случае для анализа развития процессов во времени используется двухслойная разностная схема с опережением. Метод позволяет решать самосогласованную задачу взаимодействия потока и поля, исследовать эффективность, модовый состав и спектр излучения.
2. Продольные и объемные резонансы в электродинамической системе РДГ
Резонансные свойства электродинамической системы РДГ исследовались в области частот 27г-вида основной аксиально-симметричной моды в приближении заданного тока. Электронный поток предполагался модулированным на частоте ш, воздействием поля на поток пренебре-галось. Метод поиска резонансов соответствовал предложенному в [7]: при фиксированном значении продольного волнового числа йц, варьировалась величина скорости электронов возбуждающего потока; при этом частота модуляции выбиралась в соответствии с уравнением и = Путем перебора параметров электронного
потока и структуры в широком диапазоне определялись области резонансного дифракционного рассеяния. Построенные на основе такого исследования резонансно-частотные характеристики соответствуют дисперсионным характеристикам бесконечно длинных систем.
В результате проведенных исследований было показано, что в дифракционном генераторе, как и в черепковском (действующем в области частот 7г-вида), важную роль играют продольные резонансы поверхностных компонент волн структуры, обусловленные ее конечной длиной. При одновременном возбуждении нулевой и второй пространственных гармоник основной аксиально-симметричной моды структуры распределение радиальной компоненты поля вдоль оси системы в области потока имеет максимум в центре (возбуждается продольная мода с индексом т = Ы—\). В частности, положение резонансного максимума, соответствующего 2с1/\ = 1.94, не меняется по частоте при изменении продольного волнового числа (рис. 2,о), свидетельствуя о фиксации частоты резонансного возбуждения прямой и обратной волн. Как показало решение самосогласованной задачи, частота генерации также определяется этим механизмом и вблизи 27г-вида не меняется при изменении скорости потока. Подобный механизм взаимодействия потока и поля типичен также для частоты 7г-вида, как было показано при исследовании режимов работы МВЧГ [8]. Заметим, что при взаимодействии вблизи частоты 27г-вида одновременно возбуждаются как прямая, так и обратная волны, а полное поле имеет сложную пространственную структуру с максимумом продольной компоненты на оси периодического волновода.
Поскольку для выбранной геометрии системы в области частот 27г-вида открытыми являются значительное число мод гладкого волновода, резонансно-частотная характеристика системы имеет несколько максимумов. Типичные зависимости мощности прямого излучения от частоты приведены на рис. 2,6. Кроме резонанса в непосредственной близости к 27г-виду на зависимостях присутствует резонанс, соответствующий частоте открытия
2 d/X
2.00
1.94
Ez, огн. ед.
6.1 6.2 6.3 6.4 Р, норм. ед.
5 6.6 k,d
0.8
0.4
б
1.75
1.80
1.85 1.90 2d¡\ 1.96
R, см 6
Рис. 3. Поперечные распределения поля Е, в двух резонансных максимумах: 1 — 2d/\ = 1.794, 2 — 2¿/A= 1.94
Ez, отн. ед. 0.5 г
Рис. 2. Фрагмент дисперсионной характеристики в области частот 2л--вида (а) и мощность излучения вперед при двух случаях положения потока в односекционной системе при Rmax = 6 см, Д. = 6.4: 1 — R„ = 4.7 см, 2 - R„ = 5.1 см (б)
/ 2 2
моды Еоп структуры kzn= y'(f) - (fj) , где Rw -средний радиус периодического волновода, — п-й корень уравнения /о(.г) = 0, Jq(x) — функция Бесселя. Для выбранных параметров системы п = 8. Как иллюстрирует рис. 2, б, кривая 2, величина этого резонанса сильно зависит от радиуса потока и может быть заметной при значительном удалении потока от стенок волновода. Наличие такого резонансного взаимодействия потока и поля является типичным для РДГ на цилиндрическом волноводе, соответствует взаимодействию потока с пространственными гармониками — 1-й и +1-й моды ¿oп и отсутствует для черепковского взаимодействия вблизи частоты 7г-вида. При этом для резонанса, обусловленного 27г-видом, уменьшение радиуса потока (удаление от стенок) приводит к уменьшению мощности взаимодействия, что свидетельствует о важной роли поверхностной волны. Положение резонанса, соответствующего открытию моды £ол> сильно зависит от радиуса структуры и при определенных условиях в самосогласованной задаче может использоваться для повышения мощности генерации.
Для выявления роли объемных резонансов менялся внешний радиус волновода при неизменной высоте неоднородности и фиксированном расстоянии до пучка. Для значения k,d = 1.02• 2тт радиус волновода менялся в диапазоне 5.5 4- 6.5 см с шагом 1 мм. При изменении радиуса структуры положение резонанса, соответствующего продольному резонансу поверхностной волны (резонанс периодической системы), остается относительно фиксированным, в то время как положение резонансов, связанных с открытием моды Ецп, меняется. На графике поперечных распределений полей Ег(г) в двух резонансных пиках (рис. 3) кривая 1, представляющая
-0.5
-1.0
■SA
\л
V
V
\ <v / V l> J
I
о
10
20 L, см 30
Рис. 4. Сравнение продольных полей при двух значениях внешнего радиуса: 1 — Rnmx = 6 см, 2 — Rnmx = 5.55 см
поле в выходном сечении в низкочастотном резонансе, демонстрирует подавление поверхностной компоненты и большую роль объемных компонент поля в сравнении с кривой 2, соответствующей 27г-виду основной моды периодического волновода.
Типичные зависимости распределения продольной компоненты поля Ег вдоль оси симметрии структуры при самосогласованном взаимодействии потока и поля в резонансном пике для двух значений внешнего радиуса приведены на рис. 4. При Rmax = 6.0 см (кривая 1) критическая частота открытия моды £ов отстроена от частоты 27г-вида основной моды, при Rmax = 5.5 см (кривая 2) эти частоты близки. В том случае, когда частота открытия соответствующей объемной моды близка к частоте 27г-вида (RmaK = 5.5 см), продольная структура поля имеет более сложный вид — возникает медленное изменение поля вдоль оси г, что соответствует возбуждению объемной моды £ов на частоте, близкой к критической.
Структура полей излучения, устанавливающегося в сверхразмерных структурах РДГ, достаточно сложна для анализа: ранее было показано [9], что вблизи границ высших полос прозрачности при возбуждении объемных полей в сверхразмерных периодических системах важную роль играют резонансы гибридных волн. Но из сложной многоволновой картины можно таким образом вычленить два основных механизма, определяющих процессы взаимодействия потока и поля в дифракционном генераторе.
РАДИОФИЗИКА, ЭЛЕКТРОНИКА, АКУСТИКА
33
3. Самосогласованное взаимодействие потока и поля в РДГ
Анализ самосогласованного взаимодействия потока и поля проводился для нескольких одно- и двухсекционных вариантов дифракционного генератора. Задача дополнительной оптимизации не ставилась, поэтому эффективность взаимодействия потока и поля была относительно невысокой. В односекционном варианте наблюдалась устойчивая генерация с КПД около 1%. Учет отражающего элемента (ограничителя потока) на входе в систему, обычно используемого в экспериментальных макетах, привел к повышению КПД до 2%, а секционирование системы позволило повысить КПД до 4% (рис. 5,о, кривая 1). Во всех случаях генерация была одночастотной. В распределении продольной компоненты поля Ег вдоль оси структуры (рис. 5, б) обращает на себя внимание значительное возрастание амплитуды приосевого поля в центре структуры: наличие двух периодических секций приводит к возникновению максимума поля в трубе дрейфа, играющей роль высокодобротного резонатора для объемного поля.
О 100 200
Ez, отн. ед.
300 Т/271 400
8.0x10"
4.0x10"
-4.0x10-
цтшиштшииии б
•W1
тшшишшиишг
0
15
30
45 L, см
Рис. 5. Установление генерации в двухсекционной системе: при Рптх = 6 см (кривая 1) и при Ртах = 5.55 см (кривая 2) (а); продольное распределение поля Е, на оси двухсекционной структуры (изображена в верхней части рисунка) в режиме генерации при Рпшх = 6 см (б)
Мы предположили, что при тех значениях параметров, где совмещаются объемный резонанс моды £ов с резонансным откликом на частоте ближайшего к 27г-виду резонанса периодической структуры, достигаются наиболее эффективные условия установления генерации при рассмотрении самосогласованной задачи.
Такое предположение подтвердилось увеличением КПД до 4% в односекционной системе. В аналогичном варианте двухсекционной системы значение КПД достигло 8% (рис. 5,о, кривая 2), генерация устанавливается быстрее в сравнении с Я = 6 см. Изменение эффективности генерации со временем объясняется отклонением
частоты генерации от опорной частоты oj и конкуренцией двух механизмов: объемных резонансов и резонансов поверхностной волны. Генерация одночастотная, но с более широким спектром. На продольных распределениях поля Е,- вблизи потока в каждой из секций по отдельности устанавливается (iV—1)-я мода поверхностной волны. В продольном распределении поля Ег наблюдается медленное изменение поля вдоль оси г в каждой из секций по отдельности, проявляющееся в стремлении к бесконечности длины волны соответствующей объемной моды на частоте, близкой к критической.
Отметим, что наличие отражений поверхностной волны от краев структуры и поперечные резонансы объемного поля — не единственный механизм установления обратных связей в РДГ. Существенную роль, как и в черепковском генераторе, играют внутренние обратные связи, обусловленные взаимодействием потока с обратной волной (механизм ЛОВ). В частности, в односекционной структуре с близкими параметрами и преобладающим механизмом ЛБВ (k,d ^2тт) даже при значительных величинах амплитуды тока достигался только режим усиления без перехода к генерации. Соответствующая двухсекционная структура при этом давала неустойчивую генерацию.
Заключение
В результате численного моделирования многоволнового взаимодействия трубчатого электронного потока с полями осесимметричных периодических электродинамических замедляющих систем конечной длины в области частот 27г-вида исследованы резонансные режимы самосогласованного взаимодействия потока и поля в одно- и двухсекционных системах РДГ, выявлены структуры полей, изучен модовый состав и спектр излучения. В частности:
1) установлена важная роль продольных резонансов поверхностного поля при фиксации частоты излучения в релятивистском дифракционном генераторе;
2) показано влияние на эффективность и характер генерации в РДГ режима совмещения механизмов дифракционного излучения на критической частоте собственой объемной моды эквивалентного гладкого волновода и на частоте 27г-вида основной моды структуры;
3) в численном эксперименте подтверждена перспективность использования двухсекционной сверхразмерной системы РДГ для повышения эффективности генерации высокой мощности в микроволновом диапазоне.
Список литературы
1. Бугаев С.П., Канавец В.М, Кошелев В.И., Черепенин В.А. Релятивистские многоволновые СВЧ генераторы. Новосибирск, 1991.
2. Бугаев С.П., Канавец В.И., Климов A.M и др. // Радиотехника и электроника. 1987. 32, № 7. С. 1488.
3. Бугаев С.П., Власов A.M., Канавец В.И. и др. // Релятивистская дифракционная электроника. Вып. 6. Горький, 1990. С. 185.
4. Шестопалов В.П. Дифракционная электроника. Харьков, 1976.
5. Власов А.Н., Королёва О.В., Максимов A.C. // Изв. РАН. Сер. физ. 1997. 61, № 12. С. 2289.
6. Слепков A.M. // Изв. РАН. Сер. Физ. 2003. 67, № 12. С. 1678.
7. Канавец В.М., Слепков A.M., Федоров A.B. // Вести. Моск. ун-та. Физ. Астрон. 1990. № 4. С. 30.
8. Канавец В.И., Максимов А.С., Слепков А.И. Ц Изв. РАН. 9. Галлямова О.В., Канавец В.И., Нифанов А.С. Слепков А.И. Сер. физ. 2000. 64, № 12. С. 2490. // Изв. РАН. Сер. физ. 2001. 65, № 12. С. 1703.
Peculiarities of multiwave interaction between electron beam and electromagnetic field in relativis-tic diffractional generator
A.I. Slepkov 1 . O.V. Gallyaraova2i>
1 Department of General Physics, Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia.
2Municipal Educational Institution Vidnoe Secondary General School No. 7, Moscow Region, Russia. E-mail: 11 slepkov@phys.msu.ru, bolgalglm@googlemail.com.
The mechanisms of self-consistent interaction between tubular relativistic electron beam and fields of super-dimensional axisymmetric periodical structure of relativistic diffractive generators (RDG) in the frequency range of 27r-type oscillations of fundamental mode are investigated with the help of computer simulation based on time-dependent matrix multiwave method. The electromagnetic field profiles, mode structure, and radiation frequency spectrum are studied. The conditions of the RDG generating frequency stabilisation related to the longitudinal resonances of the surface spatial harmonics and bulk resonances at the cutoff frequency of the EOn mode of the cylindrical waveguide are abtained.
Keywords: Smith-Parcel radiation, relativistic diffractional generators.
PACS: 84.40.Fe.
Received 8 December 2008.
English version: Moscow University Physics Bulletin 4(2009).
Сведения об авторах
1. Слепков Александр Иванович — докт. физ.-мат. наук, доцент, профессор; тел.: 939-30-38, e-mail: slepkov@phys.msu.ru.
2. Галлямова Ольга Валерьевна — учитель; тел.: 541-26-37, e-mail: olgalglm@googlemail.com.