Численное моделирование динамики предварительно промодулированного электронного потока в пространстве дрейфа пролетного клистрона
П.О. Медакин, Р.Н. Никулин, О.А. Авдеюк, И.В. Приходькова, Н.В. Грецова, Л.В. Дружинина
Волгоградский государственный технический университет
Аннотация: В статье указано, что изучение процесса группировки электронов в пространстве дрейфа пролетного клистрона является актуальной задачей, позволяющей установить общие закономерности, применимые к более сложным моделям. В связи с этим исследовано поведение предварительно промодулированного электронного потока в пространстве дрейфа пролетного клистрона. Реализована численная модель, которая учитывает влияние полей пространственного заряда и взаимодействие заряженных частиц с элементами электродинамической системы. Проведена серия численных экспериментов с различными значениями тока и начальных скоростей электронов, а также их сравнение с теоретическими данными. В результате проведения численных экспериментов были получены данные, характеризующие динамику электронного потока в пространстве дрейфа пролетного клистрона, при различных значениях начальной скорости (0,5с, 0.9с) и катодного тока (10мА, 1А, 10А).
Ключевые слова: пролетный клистрон, математическая модель, численное моделирование, метод крупных частиц, метод «частица — частица», пространство дрейфа, распределение конвекционного тока, электронный поток, многопоточные вычисления, система дифференциальных уравнений.
Введение
Сфера применения мощных и сверхмощных пролетных клистронов достаточно обширна, особенное значение они имеют для обеспечения космической связи и для ускорителей элементарных частиц. Поэтому актуальны задачи проектирования и создания новых таких устройств, а также модернизация существующих приборов [1-3].
Однако, реальные мощные пролетные клистроны являются сложными устройствами, поэтому задача оптимизации большого количества параметров, влияющих на работу прибора, затратная и очень сложная. Вследствие этого, в данной работе была сделана попытка рассмотреть
М Инженерный вестник Дона, №1 (2021) ivdon.ru/ru/magazine/archive/n1y2021/6771
динамику электронного потока в пространстве дрейфа пролетного клистрона в ходе численного эксперимента.
Существует большое количество электровакуумных приборов СВЧ, различных по характеристикам и схемам, однако, подробное рассмотрение показывает, что все приборы имеют схожее взаимодействие электронного потока с электромагнитными полями, которое заключается в трех основных процессах:
- модуляция электронного потока по скорости,
- модуляция электронного потока по плотности,
- переход энергии электронов в энергию СВЧ сигнала.
Протекание данных процессов в разных приборах происходит по-разному, в одних независимо друг от друга, а в других одновременно.
Одним из приборов с четким разделением данных процессов является пролетный клистрон. В зависимости от назначения клистроны разделяют на усилительные, умножительные и генераторные. Так же существует классификация по количеству резонаторов.
Схема двухрезонаторного пролетного клистрона изображена на рис. 1.
Вход СВЧ Въюод СВЧ
Рис. 1. - Схема двухрезонаторного пролетного клистрона [4]
Формирование пучка электронов производится с помощью катода и ускоряющей сетки. Подогреватель вызывает термоэлектронную эмиссию электронов с катода, а разность потенциалов и0 между сеткой и катодом ускоряет все электроны до скорости, на которую рассчитан прибор.
Таким образом, все электроны, движущиеся к первой сетке входного резонатора, имеют одинаковую скорость. Между сетками входного
резонатора действует СВЧ поле: ш , следовательно,
электроны, попадающие в эту область на одном полупериоде, тормозятся, а на другом — ускоряются. Таким образом, первый резонатор служит для модуляции электронного потока по скорости и называется группирователем.
Далее, двигаясь по инерции, электроны попадают в пространство дрейфа, в котором «быстрые» электроны догоняют «медленные», образуя так называемые сгустки электронов. Данный процесс превращает модуляцию по скорости в модуляцию электронного потока по плотности, в связи с чем, пространство дрейфа называют еще областью группирования. Важно, чтобы электронный поток из-за взаимного отталкивания электронов, не касался стенок пролетной трубы, поэтому в пространство дрейфа вводят магнитную фокусирующую систему (МФС).
Очевидно, что частота следования сгустков совпадает с частотой модулирующего СВЧ поля во входном резонаторе, поэтому, пролетая между сетками выходного резонатора, сгустки вызывают наведенный ток подобной частоты. Таким образом, выходной резонатор служит для преобразования энергии промодулированного по плотности электронного потока в энергию выходного СВЧ сигнала [5,6].
Математическая модель и численная реализация
Реальные физические модели включают огромное количество взаимодействующих друг с другом частиц, поэтому одним из методов построения модели, который позволяет экономить ресурсы ЭВМ, является замена реальной системы системой с меньшим количеством частиц, но имеющей такие же свойства. Таким образом, математическая модель движения электронного потока построена на методе «крупных частиц», заключающегося в последовательном применении Лагранжевого и Эйлерового подходов. Начальный объем ансамбля частиц разбивается на некоторое количество не пересекающихся элементарных объемов, частицы каждого из которых заменяются одной крупной частицей с суммарными зарядом и массой частиц, входящих в начальный элементарный объем.
Для численного моделирования движения одной крупной частицы, необходимо записать уравнение движения в удобном для численных расчетов виде:
Так как частица движется с релятивистской скоростью, то:
(2)
С учетом полной производной по времени, имеем:
уЧ - г (3)
)'- + - !' + - 1'=- (3) <й с2 V йХ/ т
I I/ ~
Где у = /1 — —, ¥ — сила, действующая на частицу.
Для каждой компоненты декартовой системы координат выражение примет вид:
и
V
хУ I
У г\ , У У _рх
1 + — их I + иу — 1\иу + 1'г 1?х1'г ~ ~ }х
< + ууу |
, У* Ру V -V V =
г с2 у 2 т Уу
У
ч с
+ уу — уу + уг у
Г
с'
1 + г-у}
т
fx
Гу
л
(4)
Решая данную систему методом Крамера и преобразовав решения, получим систему из шести дифференциальных уравнений первого порядка:
[Тх{уу + 1?1) - [уУхУу - Ъ ихуг]
г
^х ' "У ' 1 ?
V
X
17,
V,
г3
[/у О! + V? ) -
У3
\fziyl + ¡-у) - Тх^х^-'г - ÍУVУVA
X = V
х
(5)
У = уу 2 =
Для численного решения системы дифференциальных уравнений (5) наиболее подходящим является классический для данного рода задач — метод Рунге - Кутты четвертого порядка [7-9].
Для качественного соответствия модели реальности, при расчетах сложных физических систем, необходимо учитывать взаимодействия между частицами. В данной модели для этого используется метод «частица-частица» и запаздывающие потенциалы Лиенара-Вихерта [3]:
В итоге данная модель позволяет учесть влияние поля пространственного заряда, в то время как стандартная — кинематическая теория клистрона этого не предусматривает.
В силу того, что расчет взаимодействий и интегрирование уравнений движения для отдельных частиц проходят независимо друг от друга, их можно провести в различных потоках или процессах, что позволяет использовать многоядерность современных процессоров. Для организации многопоточности вычислений используется технология Open MP. Генерация псевдослучайных чисел осуществляется с помощью алгоритма МЛ9937[10].
Результаты
Наиболее перспективным направлением в разработке пролетных клистронов является ka - диапазон [3], в связи с чем частота принята равной 37,5 ГГц, величина фокусирующего магнитного поля - 0.3Тл. Длина модулирующего резонатора вычисляется для каждого значения начальной скорости в соответствии с величиной коэффициента модуляции, принятого равным 0.98.
Рассмотрим динамику электронного потока в пространстве дрейфа пролетного клистрона при различных начальных скоростях электронов и значениях тока. На рисунках 1, 3 и 5 представлены траектории движения электронных потоков в пространстве дрейфа двухрезонаторного пролетного клистрона в момент группировки потока по плотности. На рис. 2, 4 и 6 представлены распределения конвекционного тока вдоль оси симметрии прибора. Вертикальными линиями на рисунках изображены координаты группировки электронного потока, рассчитанные теоретически. В результате
и
анализа полученных данных сформулирован вывод, что группировка нерелятивистских малоточных электронных потоков хорошо согласуется с кинематической теорией.
Рис. 2. - Электронный поток в пространстве дрейфа с начальной скоростью 0.5с, частотой 37,5 ГГц, ток 0.01А, глубиной модуляции 50%
Распределение тока вдоль оси х. Параметры:0.5с; 0.01А; 37,5 ГГц
Теоретачеагое расположение плоскостей группировки-
1
1 !
1
] |........................................1..........................................................
1 1 ...........
0 0.02 0.04 0.06 0.08 0.1
Рис. 3. - Распределение конвекционного тока вдоль оси симметрии прибора, (начальная скорость 0.5с, частота 37,5 ГГц, ток 0.01А, глубина модуляции
50%)
Эшароиный лоток в пространств* дрейфа. Парпметры:0.9с; 1А; 37,5 П"ц
0.0015 9.001 0.0005
-0.0005 -0.001 -0.0015
50%_
Рис. 4. - Электронный поток в пространстве дрейфа с начальной скоростью 0.9с, частотой 37,5 ГГц, ток 1А, глубиной модуляции 50%
Распределение той вдоль оси г. ПарометрыгО.Эс; 1А; 37,5 ГГц
12-1-1-1-1-1-:-1-:-1-
: Теоретически ржгеиюжение шюогостеи группировки-
С 0.02 0.04 0.06 0.06 0.1 0.12 0.14
хт
Теоретическое
¡¿¿•¡f., • 1
ЕвИВД
V*,' »■ ' V * ' Vf'
V
Е-'/*i-f • J '
Ws -I Л ■■ .
V-I .'A w'i». ■
Рис. 5. - Распределение конвекционного тока вдоль оси симметрии прибора (начальная скорость 0.9с, частота 37,5 ГГц, ток 1А, глубина
модуляции 50%)
Рис. 6- Электронный поток в пространстве дрейфа с начальной скоростью 0.9с, частотой 37,5 ГГц, ток 10А, глубиной модуляции 50%
Распределение ТОЮ вдоль оси X. Параметры:0.9с; 10А; 37,5 ГГц
Теоретическое расположение гъюсхостей груши ювки
0 0.02 0.04 0.06 0-08 0.1 0.12 0.14
х,т
Рис. 7. - Распределение конвекционного тока вдоль оси симметрии прибора (начальная скорость 0.9с, частота 37,5 ГГц, ток 10А, глубина модуляции 50%)
И Инженерный вестник Дона, №1 (2021) ivdon.ru/ru/magazine/archive/n1y2021/6771
Заключение
Исследование процесса группировки электронов в пространстве дрейфа пролетного клистрона является актуальной задачей, поскольку позволяет установить общие закономерности, применимые к более сложным моделям.
Следует отметить, что рассчитываемый прибор имеет идеализированные параметры (глубина модуляции 50%, М=0.98). Данные значения были выбраны для наиболее четкого наблюдения процесса группировки электронов. Однако, использование модулирующего напряжения, сравнимого по амплитуде с ускоряющим, является нецелесообразным для реального прибора. Так же использование коэффициента модуляции, близкого к единице, говорит о малой длине модулирующего резонатора, что невыполнимо для реальных приборов из-за технологических трудностей.
В результате проведения численных экспериментов были получены данные, характеризующие динамику электронного потока в пространстве дрейфа пролетного клистрона при различных значениях начальной скорости (0,5с, 0.9с) и катодного тока (10мА, 1А, 10А).
Литература
1. Lin Fu-Min. Analysis of the optimal gap width and gap-to-gap distance in %-mode double-gap cavities for broadband klystrons // Proceedings of Progress in electromagnetics research symposium. Hangzhou, China. 2008. pp. 1353-1356.
2. Мирошниченко А.Ю., Царев В.А., Губанов В.А., Акафьева Н.А. Моделирование электродинамических параметров многомодового двухзазорного резонатора для миниатюрных многолучевых приборов клистронного типа // Инженерный вестник Дона, 2013, №3. - URL: ivdon.ru/ru/magazine/archive/n3y2013/1902.
3. Акимов П.И. Проект многолучевого Клистрона ка-диапазона длин волн, 2013. URL: mwelectronics.ru/2013/Oral/3/10_AkimovPI_ proekt_ clistrona.pdf
4. Клюев Д. С., Осипов О.В. Приборы СВЧ и оптического диапазона. Конспект лекций. Самара: ФГОБУ ВПО ПГУТИ, 2014. 209 с.
5. Рошаль А. С. Моделирование заряженных пучков. Москва: Атомиздат,1979.224с.
6. Трубецков Д. И., Храмов А.Е. Лекции по сверхвысокочастотной электронике для физиков. Москва: Физматлит, 2003. - 496с.
7. Ковтун Д. Г. Трехмерный релятивистский электронный поток в скрещенных полях // Электромагнитные волны и электронные системы. 2004. Т. 9. № 2. С. 58-65.
8. Hockney R., Eastwood J. Numerical modeling by the particle method. Computer simulation using particles. McGraw Hill,1981. 640 p.
9. Медакин П.О., Никулин Р.Н., Авдеюк О.А., Поляков И.В., Грецова Н.В. Использование численного моделирования для расчета движения ионного пучка в лазерном масс - спектрометре ЭМАЛ - 2 // Инженерный вестник Дона, 2020, №7. URL: ivdon.ru/ru/magazine/archive/N7y2020/6552
10. Быков Д.В., Неретин А.Д. Прогнозирование производительности при реализации алгоритмов генерации случайных последовательностей больших размерностей на реконфигурируемых архитектурах с сопроцессорами // Инженерный вестник Дона, 2014, №2. URL: ivdon.ru/ru/magazine/archive/n2y2014/2414.
References
1. Lin Fu-Min. Proceedings of Progress in electromagnetics research symposium. Hangzhou, China. 2008. pp. 1353-1356.
2. Miroshnichenko A.YU., Carev V.A., Gubanov V.A., Akafeva N.A. Inzhenernyj vestnik Dona, 2013, № 3. URL: ivdon.ru/ru/magazine/archive/n3y2013/1902.
3. Akimov P.I. Proekt mnogoluchevogo Klistrona ka-diapazona dlin voln [Ka-band multibeam Klystron project]. URL: mwelectronics.ru/2013/Oral/3/10_AkimovPI_ proekt_ clistrona.pdf
4. Klyuev D. S., Osipov O.V. Pribory SVCH i opticheskogo diapazona . Konspekt lekcij [Microwave and optical devices. Lecture notes]. Samara: FGOBU VPO PGUTI, 2014. 209p.
5. Roshal' A. S. Modelirovanie zaryazhennyh puchkov [Charged beam modeling]. Moskva: Atomizdat,1979. 224 p.
6. Trubeckov D. I., Hramov A.E. Lekcii po sverhvysokochastotnoj elektronike dlya fizikov [Lectures on Microwave Electronics for Physicists]. Moskva: Fizmatlit, 2003. 496 p.
7. Kovtun D. G. Elektromagnitnye volny i elektronnye sistemy. 2004. T. 9. № 2. pp. 58-65.
8. Hockney R., Eastwood J. Numerical modeling by the particle method. Computer simulation using particles. McGraw Hill,1981. 640 p.
9. Medakin P.O., Nikulin R.N., Avdeyuk O.A., Polyakov I.V., Grecova N.V. Inzhenernyj vestnik Dona, 2020, № 7. URL: ivdon.ru/ru/magazine/archive/N7y2020/6552
10. Bykov D.V., Neretin A.D. Inzhenernyj vestnik Dona, 2014, № 2. URL: ivdon.ru/ru/magazine/archive/n2y2014/2414.