Доклады БГУИР
Doklady BGUIR
2017, № 7 (109) 2017, No. 7 (109)
УДК 621.385.6
РЕЛЯТИВИСТСКИЕ ЛАМПЫ БЕГУЩЕЙ ВОЛНЫ О-ТИПА С ЛЕНТОЧНЫМИ ЭЛЕКТРОННЫМИ ПУЧКАМИ
А.В. АКСЕНЧИК, И.Ф. КИРИНОВИЧ
Белорусский государственный университет информатики и радиоэлектроники, Республика Беларусь
Поступила в редакцию 3 октября 2017
Аннотация. Сформулирована математическая модель расчета релятивистских ламп бегущей волны (ЛБВ) О-типа на волнообразно изогнутых прямоугольных волноводах с ленточными электронными пучками с учетом пространственного заряда. Проведена оптимизация параметров и получены предельные значения коэффициента полезного действия (КПД) вариантов приборов для разных размеров электронного пучка, ускоряющих напряжений 500 кВ, токов электронного пучка 200 А. Показано, что значение КПД одно- и двухсекционных ЛБВ может достигать 56 %, коэффициент усиления по мощности - 30-50 дБ.
Ключевые слова: усилитель О-типа, ленточный пучок, лампа бегущей волны, прямоугольный волновод, оптимизация.
Abstract. The mathematical model of calculation of relativistic traveling wave tubes (TWT) of O-type on folded rectangular waveguides with sheet electronic beams taking into account a space charge was formulated. Parameters optimization was spent and the maximum efficiency of variants of devices for the different sizes of the electronic beam, accelerating voltages 500 kV, currents of electronic beam 200 A were obtained. It was shown, that the efficiency of one- and two-section TWTs can reach of 56 %, gain - of 30-50 dB.
Keywords: O-type amplifier, sheet beam, traveling wave tube, rectangular waveguide, optimization.
Doklady BGUIR. 2017, Vol. 109, ]Чо. 7, pp. 66-72
Relativistic traveling wave tubes of O-type with sheet electronic beams
A.V. Aksenchyk, I.F. Kirinovich
Введение
В опубликованных ранее работах [1, 2] проводились расчеты релятивистских ЛБВ О-типа на волнообразно изогнутых прямоугольных волноводах с электронными пучками цилиндрической формы. Однако в настоящее время для повышения эффективности и мощности вакуумных приборов СВЧ предлагается использование ленточных электронных пучков [3-5]. Ленточные электронные пучки с сечением прямоугольной формы позволяют увеличить ток пучка, уменьшить плотность пространственного заряда, обеспечить более высокую эффективность приборов в широкой рабочей полосе частот.
На рис. 1, а показана схема замедляющей системы прибора [1, 2, 6], в котором осуществляется дискретное взаимодействие прямолинейного электронного потока с полем ТЕ волны волнообразно изогнутого прямоугольного волновода. Назовем такой волновод - WB (wavy bending) волновод. Электронный поток (ЭП) проходит через прямоугольные отверстия в таком волноводе посередине широкой стенки (в максимуме поперечного электрического поля). На рис. 1, б показан ленточный пучок прямоугольной формы (сечение заштриховано) в канале дрейфа размером a 1*b1. В скобках указаны составляющие системы координат волновода, изображенного на рис. 1, а. При подаче на вход СВЧ сигнала волна, проходящая по волноводу, модулирует в зазорах ЭП по скорости. При оптимальных длинах каналов дрейфа
и фазах поля в зазорах модуляция ЭП усиливается за счет взаимодействия в последующих зазорах электромагнитной (ЭМ) волны и сгруппированного электронного потока. Модуляция по скорости переходит в модуляцию по плотности, и в последних зазорах происходит отбор энергии от сгруппированных сгустков. Таким образом, взаимодействие ЭМ волны в волноводе и электронов происходит дискретно в волноводных зазорах.
•с
*Ы\
г\
\J
\J
\J
beam
X(Z) al
б
Рис. 1. Схема прибора: а - схема замедляющей системы ЛБВ; б - ленточный пучок прямоугольной формы (сечение заштриховано) в канале дрейфа размером а1*Ь1
Для обеспечения синхронизма необходимо так подобрать длины отрезков волновода и длины каналов дрейфа, чтобы электрон при движении вдоль оси У попадал бы в поле Еу ТЕ
волны одной и той же фазы. В настоящей работе описана математическая модель и проведены оптимизационные расчеты ЛБВ О-типа на WB волноводе (TWT WB) в СМ-диапазоне длин волн.
Математическая модель процессов взаимодействия электронного потока с электромагнитным полем в ЛБВ
Описанная TWT WB (рис. 1, а) моделируется цепочкой эквивалентных четырехполюсников. Один четырехполюсник моделирует один изгиб волнообразно изогнутого прямоугольного волновода. Подробно цепочка эквивалентных четырехполюсников описана в [2, 6, 7].
На входе цепочки, слева, подключен генератор входного сигнала Е0 с внутренним сопротивлением 20, равным эквивалентному сопротивлению волновода ZW . Затем следует согласующий четырехполюсник М0, описывающий отрезок волновода до первого зазора. После последнего зазора, справа, подключен согласующий Мп+1 четырехполюсник для согласования волновода с нагрузкой 2п. Ввиду того, что четырехполюсники моделируют отрезки одного и того же волновода, они оказываются согласованными и при изменении частоты. Будем считать, что нагрузка с сопротивлением 2п согласована с волноводом, имеющим эквивалентное сопротивление ZW на опорной частоте.
Матрица передачи А четырехполюсника, моделирующего изогнутый отрезок волновода длинной £ , имеет следующий вид [6, 7]:
A =
ch( jti) Zw • sh( jti)
Z" sh( j ti)
ch( jti)
(1)
продольное
где Г = Г' — уГ" - постоянная распространения волны; Г = К-у/Т—Хт/К -волновое число; К = юф.0ц0 ; х2 = Х2тп = (тк/ а)2 +(пк/ Ь)2 - поперечное волновое число;
R
t =
i+2b (v* с )2
- коэффициент затухания для волны ТЕ10 в прямоугольном волноводе;
Zо • Ьл
ф-(х / х. )2
Xсг - критическая длина волны в волноводе, для волны ТЕ10 Xсг = 2а; X = с / / - длина волны
входного сигнала, / - его частота,
скорость света в вакууме;
7 л Ь л/^/80 ТЕ
=-----, = - эквивалентное сопротивление волновода для волны ТЕ10; в0, ц 0 -
2 а
п Л / Ц0
соответственно диэлектрическая и магнитная проницаемости вакуума; к8 = /
g
g = 5,7 106сим/м - проводимость меди; ¿0 = 120л ; а - размер широкой стенки волновода; Ь - размер узкой стенки волновода.
Задача возбуждения TWT WB решается методом наложения, последовательно определяются амплитуды напряжений на зазорах, вызванные входным сигналом Е0 и наведенными (конвекционными) токами. Конвекционный ток сгруппированного электронного потока можно определить совместным решением уравнения движения и уравнения возбуждения волновода [7]. Подробно алгоритм расчета возбуждения цепочки четырехполюсников с учетом прямой и обратной волн описан в [6, 7].
Уравнение возбуждения волновода. Используем обычную для волноводов систему координат: ось Z направлена вдоль волновода, ось У - перпендикулярна широкой стенке волновода и совпадает с направлением движения электронов (рис. 1, а), ось Х -перпендикулярна узкой стенке волновода.
Представим возбужденное поле в виде суперпозиции полей свободных волн:
Е=У с± Л р , н=У с± А р . (2)
р р
Здесь под индексом р понимаем два индекса, которые определяют тип поля в волноводе. Используя уравнения Максвелла, применяя лемму Лоренца, можно получить уравнение возбуждения волновода [7]. С учетом того, что будем рассматривать возбуждение и распространение в прямоугольном волноводе волны поля ТЕтп, уравнение возбуждения
конвекционным током ./ принимает вид [7]:
Ср = 11J Ер^ , (3)
V
где знак ( + ) соответствует волнам, движущимся в положительном направлении оси Z рис. 1, а, знак ( - ) - в противоположном направлении.
Задачу решаем в одномерном приближении, т.е. поперечное движение электронов учитывать не будем. Для моделирования электронного потока используем метод крупных частиц. Электронный поток представляем состоящим из N частиц прямоугольной формы (рис. 1, б), распределенных, при у = 0, равномерно на периоде 0 ^ 2л .
Рассчитываем интеграл (3) численно, используя метод средних и допуская, что по сечению плотность тока постоянна - 11 = Jm (0) • г • h . Получим:
21
N0
Jш =±11 у . (4)
В данной модели волновод возбуждается поперечным током Jш, коэффициенты С±тп вне этого объема от координаты г не зависят и постоянны. Обозначим амплитуду возбужденного поля как Е = — ]С±тп%х /%2 . Дальнейшие преобразования подробно описаны в [7]. Окончательно получим выражение для наведенного тока в ^м зазоре с учетом безразмерных параметров:
Т
^ = Ь+1 у Г8°(Т-Т0) • )ат . (5)
ик У 0 ^-=11
Входящие в (5) безразмерные параметры описаны ниже.
Уравнения движения. Электронный поток моделируется с применением метода крупных частиц. Разбиваем электронный поток вдоль направления движения (по координате У) на Ne крупных частиц. Форма частицы приведена на рис. 1, б. Такая форма легко получается из формы частицы, приведенной в [8], путем задания соответствующих размеров [8].
Уравнение движения электрона во внешнем электромагнитном поле без учета торможения излучением (т. е. до у /c < 0,99) в одномерном приближении имеет следующий вид (учтем, что движение электронов в данной математической модели происходит по координате Y (рис. 1, а)):
dv
~г _ —^о
1-
• Ey, (6)
у
где = е / т0, е, т0 - соответственно заряд и масса покоя электрона, V - продольная скорость электрона, у - начальная скорость электронного потока.
Представим напряженность поля как Еу = Е0 + Ер, где Е0 = Re(EmeJ)) -
напряженность ВЧ поля в зазоре волновода (6), Ер - напряженность продольного поля пространственного заряда. Поле пространственного заряда Ер рассчитываем исходя из моделирования электронного потока Ые крупными частицами.
Уравнения движения крупной частицы в нормированных переменных, алгоритм возбуждения цепочки четырехполюсников аналогичны тем, что описаны в [6, 7]. Для расчета полей пространственного заряда плоских электронных пучков использовались формулы, полученные в [8]. Введем следующие безразмерные переменные:
T _ у / L, З = at, у 0 =
1- W
c
— , с La
V =—, u = at, -аУ/Зо _—, у, _ v0 v0
-1/2
v
1- (^)2 c
-1/2
З = ™ k, (7)
(Уо-1) ^. ; U± м = E d a = eI0
U ± ~ -
_ _ -. ^ - _. ; _ \ Л m Uj± _ EJ d a _
a ~ r, , ' Sk / j TT ' m ±mn ' p 2 Ml
2лр0 d m U0 a л-0m080a b
где y - продольная координата электрона, L - длина ЛБВ, I k - длина отрезка волновода; a'b' - размеры канала дрейфа прямоугольного сечения, d-длина зазора в волноводе, равна ширине узкой стенки волновода b, a - опорная частота, ;к - относительное напряжение на зазоре, полученное с учетом прямых и обратных волн, возбуждаемых электронным потоком в других зазорах, и рассчитываемое по алгоритмам, приведенным в [6, 7].
После нормировки (6) получим релятивистское уравнение движения крупной частицы:
dV _ -^Ha Re&e'("'+v+9k))b0(T -Tok) + apFp], (8)
dT У ,V
£ _ ^-1). (9)
dT —
где
1 | u , - ui |
Fp _ N E 8 p ("-J~)sign(uj - ut), (10)
80(T - T0) - весовая функция пространственного распределения поля зазора, рассчитывается
методом сеток, 8 (у) - весовая функция распределения напряженности поля
пространственного заряда, определяется с использованием периодизированных полей крупных частиц прямоугольной формы (подробно описана в [6, 8]).
Уравнения (5), (8), (9) образуют самосогласованную нелинейную систему уравнений TWT WB. После определения наведенных в зазорах волновода токов Jk, возбуждающих соответствующие четырехполюсники, используется алгоритм расчета ЛБВ, описанный в [6, 7]. Заметим, что для выполнения условия синхронизма, т. е. постоянства фазы СВЧ поля при прохождении электронным сгустком зазоров в волноводе (рис. 1, а), использовалась оптимизация распределения коэффициента фазы четырехполюсников (эквивалентных отрезкам волновода). В четырех-шести точках задавалось первое приближение распределения коэффициента фазы на всех четырехполюсниках, а затем для аппроксимации использовался аппарат атомарных функций. Подробное описание этой методики приведено в [7].
2 \3/2
Результаты расчетов и их обсуждение
При выборе размеров прямоугольных волноводов, каналов дрейфа, ускоряющих напряжений, токов луча необходимо учитывать следующие факторы.
Угол пролета электрона в зазорах волноводов должен быть мал. В СМ-диапазоне при размерах узкой стенки 0,5-1,0 см лишь при ускоряющих напряжениях более 200 кВ можно получить приемлемые результаты. При выборе тока электронного луча необходимо учитывать следующее. В связи с небольшим эквивалентным сопротивлением волновода и низким сопротивлением связи при малом токе электронного луча не удается получить эффективного взаимодействия с ЭМ волной. Только при превышении некоторого порога тока луча получаются приемлемые коэффициент усиления и КПД более 10 %. Таким образом, TWT WB в СМ-диапазоне являются приборами большой мощности.
Проведены расчеты релятивистских TWT WB с ускоряющим напряжением 500 кВ, 10 = 200 А, длиной волны входного сигнала X = 6 см, размерами волновода - а = 4,8 см, Ь = 1 см. Оптимизировались длины каналов дрейфа и коэффициенты фаз четырехполюсников. Чтобы не было электромагнитной связи между щелями связи по каналам дрейфа, размеры канала дрейфа необходимо выбирать так, чтобы не было распространяющихся мод (Ню) по каналу [6]. Поэтому выбираем размеры канала дрейфа а1 х Ь1 = 2,5 х 0,6 [см]. Размер пучка электронов hx г = 2х 0,4 [см]. Тогда минимальная длина волны, на которой может работать усилитель без возбуждения, должна быть больше 5 см. Для числа зазоров N = 33 при оптимизации получен электронный КПД це = 0,566, коэффициент усиления Кр = 34,4 дБ - вариант А1. Частотная характеристика,
зависимость коэффициента усиления по мощности Кр от относительной частоты w = f / ^
приведена на рис. 2, кривая 1. Полоса частот усиления около 8-9%.
Кр [с1В]
50
2
1
Г/ у. д. , / \
^^ / / 3 —\ V Л 1 / 1
ЛЛ/ ^ / \уУ
0,91
0,94
0,97
1,03
1,06
1,09
Рис. 2. Зависимость коэффициента усиления по мощности Кр от относительной частоты w= f /
1 - вариант А1, 2 - вариант А2, 3 - вариант А3
Проведен расчет TWT WB, вариант А2, для размеров канала дрейфа а1х Ь1 = 2х 0,6 [см], размер пучка электронов hx г = 1,5 х 0,4 [см]. Для числа зазоров N = 33 получен электронный КПД це = 0,565, коэффициент усиления Кр =32,2 дБ - вариант В2. Частотная характеристика
этого варианта приведена на рис. 2, кривая 2. Полоса частот усиления - около 8-9 %.
Проведен расчет TWT WB, вариант А3, для размеров канала дрейфа а1х Ь1 = 1,2 х 1,2 [см], размер пучка электронов h х г = 0,8 х 0,8 [см]. Для числа зазоров N = 33 получен электронный КПД це =0,524, коэффициент усиления Кр =33,9 дБ. Частотная характеристика этого варианта
приведена на рис. 2, кривая 3. Этот вариант эквивалентен TWT WB с пучком цилиндрической формы с соответствующими размерами пучка и канала дрейфа. Полоса частот усиления почти в два раза меньше, чем в вариантах А1, А2 с ленточными пучками, рис. 2.
Проведены расчеты релятивистской двухсекционной [6] с последовательным соединением секций. Ускоряющее напряжение 500 кВ, !0 = 200 А. Первая секция - регулярная, число зазоров N = 17, вторая секция - нерегулярная, число зазоров N = 17. Размеры канала дрейфа а1х Ь1= 2,5 х 0,6 [см], размер пучка электронов h х г = 2х 0,4 [см]. Получен электронный КПД це = 0,598, коэффициент усиления Кр = 54,14 дБ - вариант А4.
Заключение
Предложены математическая модель и алгоритм расчета релятивистских ЛБВ на волнообразно изогнутом волноводе с ленточными электронными пучками. Использовалась методика оптимизации распределений коэффициентов фаз эквивалентных четырехполюсников. Для аппроксимации распределений коэффициентов фаз использовался аппарат атомарных функций. Показано, что в СМ-диапазоне расчетный электронный КПД достигает 56 %, коэффициент усиления по мощности односекционных ЛБВ - 34 дБ, двухсекционных - 54 дБ, полоса усиления - 8 %. Полоса частот усиления ЛБВ с пучком цилиндрической формы почти в два раза меньше, чем в ЛБВ с ленточными пучками. Распределение КПД и усиления в полосе частот - неравномерно.
Список литературы
1. Аксенчик А.В., Кураев А.А. Релятивистские ЛБВ на волнообразно изогнутом прямоугольном волноводе // Материалы 13-й Междунар. Крымской микроволновой конф. «КрыМиКо». Севастополь, 2003. С. 241-242.
2. Аксенчик А.В. Моделирование и оптимизация нерегулярных ЛБВ, ЛОВ типа «О» // Докл. БГУИР. 2004. № 2 (6). С. 78-90.
3. Reseach on X-band sheet beam electron optics system / Ding Zhao [et al.] // IVEC 2012. P. 173-174.
4. Theory and experimental investigation on the high performance transport of sheet electron beam for the XSBK and XSBK / Cunjun Ruan [et al.] // IVEC 2012. P. 147-148.
5. Multi- kW sheet beam amplifiers at Ka and W bands / Pascour J. [et al.] // IVEC 2014. P. 19-20.
6. Аксенчик А., Кураев А., Киринович И. Усилители и генераторы О-типа миллиметрового и субмиллиметрового диапазонов длин волн (теория и оптимизация). Германия: Palmarium Academic Publishing, 2014. 154 c.
7. Аксенчик А.В., Кравченко В.Ф. Оптимизация ЛОВ и ЛБВ О-типа на волнообразно изогнутых прямоугольных волноводах для миллиметрового и субмиллиметрового диапазона длин волн // Успехи современной радиоэлектроники. 2007. № 10. С. 46-56.
8. Аксенчик А.В., Киринович И.Ф. Поля пространственного заряда плоских пучков // Вес. Нац. акад. навук Беларуси Сер. фiз.-тэхн. навук. 2011. № 2. С. 97-107.
References
1. Aksenchik A.V., Kuraev A.A. Reljativistskie LBV na volnoobrazno izognutom prjamougol'nom volnovode // Materialy 13-j Mezhdunar. Krymskoj mikrovolnovoj konf. «KryMiKo». Sevastopol', 2003. S. 241-242. (in Russ.)
2. Aksenchik A.V. Modelirovanie i optimizacija nereguljarnyh LBV, LOV tipa «О» // Dokl. BGUIR. 2004. № 2 (6). S. 78-90. (in Russ.)
3. Reseach on X-band sheet beam electron optics system / Ding Zhao [et al.] // IVEC 2012. P. 173-174.
4. Theory and experimental investigation on the high performance transport of sheet electron beam for the XSBK and XSBK / Cunjun Ruan [et al.] // IVEC 2012. P. 147-148.
5. Multi- kW sheet beam amplifiers at Ka and W bands / Pascour J. [et al.] // IVEC 2014. P. 19-20.
6. Aksenchik A., Kuraev A., Kirinovich I. Usiliteli i generatory O-tipa millimetrovogo i submillimetrovogo diapazonov dlin voln (teorija i optimizacija). Germanija: Palmarium Academic Publishing, 2014. 154 s. (in Russ.)
7. Aksenchik A.V., Kravchenko V.F. Optimizacija LOV i LBV O-tipa na volnoobrazno izognutyh prjamougol'nyh volnovodah dlja millimetrovogo i submillimetrovogo diapazona dlin voln // Uspehi sovremennoj radiojelektroniki. 2007. № 10. S. 46-56. (in Russ.)
8. Aksenchik A.V., Kirinovich I.F. Polja prostranstvennogo zarjada ploskih puchkov // Ves. Nac. akad. navuk Belarusi. Ser. fiz.-tjehn. navuk. 2011. № 2. S. 97-107. (in Russ.)
Сведения об авторах
Аксенчик А.В., д.ф.-м.н., профессор, профессор кафедры вычислительных методов и программирования Белорусского государственного университета информатики и радиоэлектроники.
Information about the authors
Aksenchik A.V., D.Sci., professor, professor of Computational methods and programming department of Belarusian state university of informatics and radioelectronics.
Киринович И.Ф., к.ф.-м.н., доцент кафедры инженерной психологии и эргономики Белорусского государственного университета информатики и радиоэлектроники.
Адрес для корреспонденции
220013, Республика Беларусь, г. Минск, ул. П. Бровки, д. 6, Белорусский государственный университет информатики и радиоэлектроники тел. +375-44-702-00-95; e-mail: aksenchik@bsuir.by Аксенчик Анатолий Владимирович
Kirinovich I.F., PhD, associate professor of engineering psychology and ergonomics department of Belarusian state university of informatics and radioelectronics.
Address for correspondence
220013, Republic of Belarus, Minsk, P. Brovka st., 6, Belarusian state university of informatics and radioelectronics tel. + 375-44-702-00-95; e-mail: aksenchik@bsuir.by Aksenchyk Anatoly Vladimirovich