МИНИСТЕРСТВО НАУКИ И ВЫСШЕГО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ
ВЫСШЕГО ОБРАЗОВАНИЯ
«НОВОСИБИРСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»
Электрофизических установок и ускорителей
Кафедра _______________________________________________________________________
(полное название кафедры)
Утверждаю
ЭФУиУ
Зав. кафедрой _______________
А.В. Бурдаков
_____________________________
(подпись, инициалы, фамилия)
«___» _______________ 2021г.
ВЫПУСКНАЯ КВАЛИФИКАЦИОННАЯ РАБОТА БАКАЛАВРА
Бердников Константин Николаевич
_______________________________________________________________________________
(фамилия, имя, отчество студента – автора работы)
Расчет динамики электронного пучка в линейном ускорителе ВЭПП-5
_______________________________________________________________________________
(тема работы)
_______________________________________________________________________________
Физико-технический факультет
_______________________________________________________________________________
(полное название факультета)
03.03.02 Физика
Направление подготовки _________________________________________________________
(код и наименование направления подготовки бакалавра)
_______________________________________________________________________________
Руководитель
от НГТУ
Руководитель
от организации
Автор выпускной
квалификационной работы
______________________________________
(фамилия, имя, отчество)
______________________________________
(фамилия, имя, отчество)
Арсентьева Мария Васильевна
Бердников Константин Николаевич
д.ф.-м.н.
______________________________________
(ученая степень, ученое звание)
______________________________________
(ученая степень, ученое звание)
ФТФ, ФТ-71
______________________________________
(факультет, группа)
______________________________________
(подпись, дата)
______________________________________
(подпись, дата)
______________________________________
(подпись, дата)
Аракчеев Алексей Сергеевич
Новосибирск 2021
______________________________________
(фамилия, имя, отчество)
МИНИСТЕРСТВО НАУКИ И ВЫСШЕГО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ
ВЫСШЕГО ОБРАЗОВАНИЯ
«НОВОСИБИРСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»
Электрофизических установок и ускорителей
Кафедра _______________________________________________________________________
(полное название кафедры)
УТВЕРЖДАЮ
Введите
здесь
свой
текст
Зав. кафедрой __________________
ЭФУиУ
А.В.
Бурдаков
(фамилия, имя, отчество)
__________________
(подпись, дата)
ЗАДАНИЕ
НА ВЫПУСКНУЮ КВАЛИФИКАЦИОННУЮ РАБОТУ БАКАЛАВРА
Бердникову Константину Николаевичу
студенту _________________________________________________________________
(фамилия, имя, отчество)
03.03.02 Физика
Направление подготовки ____________________________________________________
(код и наименование направления подготовки бакалавра)
__________________________________________________________________________
Физико-технический факультет
__________________________________________________________________________
(полное название факультета)
Расчет динамики электронного пучка в линейном ускорителе ВЭПП-5
Тема _____________________________________________________________________
(полное название темы выпускной квалификационной работы бакалавра)
__________________________________________________________________________
Анализ продольной динамики электронного
Исходные данные (или цель работы) __________________________________________
пучка
от
момента
инжекции
до
режима
стабильного ускорения.
__________________________________________________________________________
Измерение мощностей клистронов ВЭПП-5.
__________________________________________________________________________
Моделирование электромагнитных полей от электронной пушки до конца первой
__________________________________________________________________________
ускоряющей секции. Моделирование продольной динамики пучка в программе
__________________________________________________________________________
ASTRA, определение оптимального режима работы
__________________________________________________________________________
Структурные части работы __________________________________________________
1. Введение
2. Схема элементов начальной части линейного ускорителя инжекционного
__________________________________________________________________________
__________________________________________________________________________
комплекса
__________________________________________________________________________
2.1. Элементы ускоряющей и магнитной системы
__________________________________________________________________________
2.2. Система питания
__________________________________________________________________________
3. Анализ продольной динамики
__________________________________________________________________________
3.1. Уравнения продольной динамики в поле бегущей волны
3.2. Условие захвата в режим ускорения
___________________________________________________________________________
3.3. Теоретические оценки захвата частиц в ускорение полем бегущей волны
___________________________________________________________________________
4. Измерение мощностей клистронов
___________________________________________________________________________
4.1. Измерение с помощью диодного детектора мощности
___________________________________________________________________________
5. Моделирование продольной динамики пучка
___________________________________________________________________________
5.1. Параметры пучка
___________________________________________________________________________
5.2. Расчет энергии частицы при пролете ускоряющей секции с учетом SLED
___________________________________________________________________________
5.3. Варьирование мощности группирователя
___________________________________________________________________________
6. Заключение
___________________________________________________________________________
7. Список литературы
___________________________________________________________________________
___________________________________________________________________________
8. Приложение 1. Поля магнитных элементов
___________________________________________________________________________
Задание согласовано и принято к исполнению.
Руководитель
от НГТУ
Руководитель
от организации
Аракчеев Алексей Сергеевич
Арсентьева Мария Васильевна
Студент
Бердников Константин Николаевич
______________________________________
(фамилия, имя, отчество)
______________________________________
(фамилия, имя, отчество)
______________________________________
(фамилия, имя, отчество)
д.ф.-м.н.
______________________________________
(ученая степень, ученое звание)
______________________________________
(ученая степень, ученое звание)
ФТФ, ФТ-71
______________________________________
(факультет, группа)
______________________________________
(подпись, дата)
______________________________________
(подпись, дата)
______________________________________
(подпись, дата)
Тема утверждена приказом по НГТУ № _________ от «____» ___________ 2021г.
изменена приказом по НГТУ № _________ от «____» ___________ 2021г.
ВКР сдана в ГЭК № _______, тема сверена с данными приказа
7
7
___________________________________________________
(подпись секретаря государственной экзаменационной комиссии по защите ВКР, дата)
_________________________________________________
(фамилия, имя, отчество секретаря государственной
экзаменационной комиссии по защите ВКР)
Аннотация
Инжекционный комплекс ВЭПП-5 представляет собой интенсивный источник электронных и позитронных сгустков на энергии 430 МэВ, обеспечивающий все потребности работающих в Институте ядерной физики СО РАН
установок на встречных электрон-позитронных пучках. В данный момент
наблюдаются значительные потери пучка в линейном ускорителе между электронной пушкой и спектрометром, установленным после второй ускоряющей
секции. С целью возможного повышения эффективности ВЭПП-5 проведен
анализ продольной динамики электронного пучка от момента инжекции с
начальной энергией 200 кэВ до режима стабильного ускорения в первой ускоряющей структуре. Получены теоретические оценки захвата пучка в ускорение. Проведены оценки набираемой энергии пучком в ускоряющей структуре
с учетом системы умножения мощности типа SLED. Представлены измерения
мощности и спектров клистронов при помощи диодного детектора мощности,
построены распределения электрических и магнитных полей линейного ускорителя ВЭПП-5 от пушки до конца первой ускоряющей секции. В полученных
полях проведено моделирование динамики пучка при различных амплитудах
электрического поля группирователя. Показано влияние поля группирователя
на полное токопрохождение, а также на захват частиц с учетом необходимого
энергетического разброса.
Annotation
The VEPP-5 injection complex is an intense source of electron and positron
bunches at an energy of 430 MeV, which provides all the needs of the installations
operating at the Institute of Nuclear Physics of the SB RAS on colliding electronpositron beams. At the moment, significant beam losses are observed in the linear
accelerator between the electron gun and the spectrometer installed after the second
accelerating section. In order to possibly increase the efficiency of VEPP-5, the analysis of the longitudinal dynamics of the electron beam from the moment of injection
with an initial energy of 200 keV to the regime of stable acceleration in the first
accelerating structure is carried out. Theoretical estimates of beam capture into acceleration are obtained. The energy gain by the beam in the accelerating structure is
estimated with allowance for the SLED-type power multiplication system. Measurements of the power and spectra of klystrons using a diode power detector are presented, the distributions of electric and magnetic fields of the VEPP-5 linear accelerator from the gun to the end of the first accelerating section are constructed. In the
fields obtained, the beam dynamics was simulated at various amplitudes of the electric field of the buncher. The effect of the bunching field on the total current passage,
as well as on the capture of particles, taking into account the required energy spread
is shown.
Оглавление
1
Введение ................................................................................................... 7
2 Схема элементов начальной части линейного ускорителя
инжекционного комплекса ..................................................................................... 9
3
2.1
Элементы ускоряющей и магнитной системы ............................... 9
2.2
Система питания ............................................................................. 11
Анализ продольной динамики ............................................................. 13
3.1
Уравнения продольной динамики в поле бегущей волны .......... 14
3.2
Условие захвата в режим ускорения ............................................. 15
3.3
волны
Теоретические оценки захвата частиц в ускорение полем бегущей
19
4
Измерение мощностей клистронов...................................................... 21
4.1
5
Измерение с помощью диодного детектора мощности .............. 21
Моделирование продольной динамики пучка.................................... 24
5.1
Параметры пучка ............................................................................. 24
5.2 Расчет энергии частицы при пролете ускоряющей секции с
учетом SLED ...................................................................................................... 25
5.3
Варьирование мощности группирователя .................................... 27
6
Заключение............................................................................................. 31
7
Список литературы................................................................................ 32
8
Приложение 1. Поля магнитных элементов. ...................................... 33
1 Введение
Инжекционный комплекс ВЭПП-5 (Рис. 1) — это интенсивный источник
электронных и позитронных сгустков на энергии 430 МэВ, обеспечивающий
все потребности работающих в Институте ядерной физики СО РАН установок
на встречных электрон-позитронных пучках.
Рис. 1. Схема инжекционного комплекса ВЭПП-5.
1 – ускоритель электронов; 2 – ускоритель позитронов; 3 – накопитель-охладитель.
Комплекс состоит из линейного ускорителя электронов на энергию
285 МэВ, линейного ускорителя позитронов на 430 МэВ, а также накопителяохладителя с каналами впуска и выпуска пучков в коллайдеры института ядерной физики.
Линейный ускоритель также имеет участок поворота частиц, поэтому
энергетический разброс пучка ограничен акцептансом поворотного магнита,
что в данном случае накладывает ограничение на разброс энергии около 1.5%.
Параметры пучка в линейном ускорителе представлены в Таблица 1.
Таблица 1. Параметры пучка в линейном ускорителе.
Энергия инжекции
Длительность пучка (ширина на полувысоте)
Энергия на выходе линейного ускорителя
Число электронов
Энергетический разброс
200 кэВ
4.6 нс
285 МэВ
2 × 10
1.5%
7
Рис. 2. Осциллограмма тока пучка на выходе пушки (одна клетка по оси ординат соответствует 0.85 А)
Эффективность работы коллайдеров повышается с увеличением числа частиц в пучке, поэтому важной задачей является поддержание высокого уровня
числа частиц в пучке.
Пучок после прохождения второй ускоряющей секции попадает в спектрометр, который показывает, что на данный момент на него опадает около
50% от числа частиц в пушке. В рамках данной работы ставится задача определить, какое влияние на продольную динамику оказывают элементы ускорителя и показать, какого токопрохождения можно добиться в предположении
отсутствия потерь частиц за счет поперечной динамики. Результаты работы в
будущем могут быть использованы для повышения эффективности работы инжекционного комплекса.
Таблица 2. Параметры пучка после прохождения второй ускоряющей секции.
Энергия
Число электронов
100 МэВ
1 × 10
Наибольшие потери пучка в линейном ускорителе ВЭПП-5 происходят на
том этапе, пока скорость пучка равна 𝛽𝑐 ≈ 0.7𝑐. Поэтому в данной работе рассматривается участок ускорителя от электронной пушки до конца первой
ускоряющей секции. На участке от конца первой ускоряющей секции до конца
второй ускоряющей секции располагается канал с квадруполями, в котором
пучок уже сформирован, поэтому на данном участке происходят небольшие
потери пучка. Основные потери пучка происходят на участке, где происходит
8
формирование пучка: от электронной пушки до конца первой ускоряющей
секции. Поэтому дальнейшее рассмотрение продольной динамики пучка ведется до конца первой ускоряющей структуры.
2 Схема элементов начальной части линейного ускорителя инжекционного комплекса
2.1 Элементы ускоряющей и магнитной системы
Схема элементов ускорительного тракта от пушки до первой ускоряющей
секции (УС1) представлена на (Рис. 3).
Рис. 3. Схема элементов ускорительного тракта до начала первой ускоряющей секции.
1 – катод; 2 – магнитные линзы; 3 – магнитные экраны; 4 – набор токовых колец
(тонких соленоидов); 5 – первая ускоряющая секция; 6 – соленоид первой ускоряющей
секции; 7 – группирователь.
Группирователь (Таблица 4, Рис. 4) [1] выполнен на основе резонансных
ячеек со стоячей электромагнитной волной. Он состоит из четырех связанных
цилиндрических резонаторов (3 чашки + 1 трансформатор типа волны). В таком группирователе непрерывный начальный пучок преобразуется в последовательность сгустков с частотой регулярных ускоряющих структур для дальнейшего захвата в режим ускорения.
9
Таблица 3. Соответствие входной мощности
группирователя амплитуде электрического поля.
Амплитуда электрического поля, МВ/м
Мощность группирователя, кВт
0.5
1
1.5
2
2.5
3
3.5
4
4.5
0.9
3.6
8.1
14.4
22.4
32.3
44.0
57.4
72.7
Таблица 4. Параметры резонаторов группирователя.
Собственная добротность
Коэффициент связи
Шунтовое сопротивление
3349
1.3
0.4 МОм
Рис. 4. Распределение продольного электрического поля на оси группирователя при
входной мощности 8.1 кВт.
В инжекционном комплексе используются ускоряющие структуры (Таблица 5) с постоянным импедансом и набегом фазы на ячейку 2π/3, работающие
на бегущей волне. Структуры изготовлены на основе диафрагмированного
волновода, и работают на бегущей волне с продольной составляющей электрического поля.
10
Таблица 5. Параметры ускоряющих структур линейного ускорителя.
Параметры
Рабочая частота
Рабочий вид колебаний
Собственная добротность
Значения
2855.5 МГц
1.3 104
Период
34.99 мм
Внутренний диаметр резонатора
83.75 мм
Диаметр отверстия диафрагмы
25.9 мм
Толщина диафрагмы
Длина
6 мм
2.93 м
Количество ускоряющих ячеек
(включая два трансформатора типа
волны)
85
Фазовая скорость (с – скорость света
в вакууме)
1𝑐
Групповая скорость (с – скорость
света в вакууме)
Шунтовое сопротивление
0.021c
51 МОм/м
Коэффициент затухания по полю
0.108 м-1
Время заполнения структуры
0.456 мкс
2.2 Система питания
Четыре клистрона SLAC 5045 являются источниками питания 14 ускоряющих секций линейного ускорителя. Каждый клистрон обеспечивает СВЧ импульсы длительностью 3.5 мкс. В течение 3 мкс мощность запасается в резонаторах системы SLED (SLAC Energy Development) [2], после чего происходит
переворот фазы волны генератора. В последующие 0.5 мкс мощность, накопленная резонаторами системы SLED, поступает в ускоряющие структуры.
Мощность от клистрона 3 после накопления в системе SLED делится двумя 3
дБ мостами между первыми тремя ускоряющими структурами в соотношении
2:1:1, причем малая часть мощности поступает в группирователь. Таким обра11
зом, в первую структуру направляется около половины СВЧ мощности, отраженной от системы SLED, а оставшаяся половина делится поровну между второй и третьей ускоряющими структурами.
Система умножения мощности типа SLED состоит из волноводного трехдецибельного щелевого моста и двух высокодобротных цилиндрических резонаторов (Рис. 5). Диаметр резонатора D = 196 мм, высота H = 346.9 мм. В качестве рабочей моды колебания выбрана мода Н015, поскольку она имеет на
рабочей частоте 2855 МГц наибольшую собственную добротность резонаторов порядка 𝑄 ≈ 100000 [3]. При этом на данной частоте существует также
дополнительная мода колебания Е115 с добротностью порядка 𝑄 ≈ 50000.
Для снятия вырождения резонатор сделан не полностью цилиндрическим (см.
Рис. 5). За счет такого изменения формы изменяется емкость резонатора на
моде колебания Е115, что приводит к изменению ее резонансной частоты.
Кроме этого, ввод мощности осуществляется через отверстие связи, смещенное относительно центра резонатора, что также способствует снятию вырождения.
Рис. 5. Система умножения мощности SLED.
12
Рис. 6. График отраженной мощности от времени.
Волноводный трехдецибельный мост представляет собой направленный
ответвитель первого рода. В стационарном режиме амплитуды отраженных от
него волн 𝑏 = (𝑏 , 𝑏 , 𝑏 , 𝑏 ) связаны с амплитудами падающих волн 𝑏 =
(𝑏 , 𝑏 , 𝑏 , 𝑏 ) матрицей рассеяния S [4]:
𝑏⃗ =
𝑏
𝑏
𝑏
𝑏
0
= 𝑆 ∙ 𝑎⃗ = 0
1
−𝑖
𝑎
1 −𝑖
−𝑖 1 ∙ 𝑎
𝑎
0 0
𝑎
0 0
0
0
−𝑖
1
Зависимость отраженной мощности от системы умножения мощности
SLED от времени выражается формулой (1).
𝑃(𝑡) = 𝑃 ∙
( )
− 𝑎(𝑡)
∙ ( )
∙ ∙
∙
+
∙
∙ ( )
∙
, (1)
где 𝑄 – добротность резонаторов, a(t) – нормированный сигнал от клистрона,
k – коэффициент связи, 𝛼 – коэффициент затухания ускоряющей структуры.
3 Анализ продольной динамики
На выходе пушки полная длительность пучка составляет около 7 нс, поэтому для более эффективной работы ускорителя пучок должен быть по большей части сгруппирован перед влетом в ускоряющую секцию. С учетом того,
что фазовая скорость пучка меньше скорости света, при влете в ускоряющую
структуру он начинает скользить по фазе. Этот процесс будет продолжаться
до того момента, пока скорость пучка не станет равной фазовой скорости
волны, и сопровождаться группировкой пучка.
Рассмотрим более подробно, что происходит с пучком с начальной энергией 200 кэВ при инжекции в ускоряющую структуру без предварительной
13
группировки [5]. Данный анализ позволит в дальнейшем судить о том,
насколько эффективна работа группирователя в линейном ускорителе ВЭПП5.
3.1 Уравнения продольной динамики в поле бегущей волны
Уравнение изменения энергии заряженной частицы на оси ускоряющей
структуры на основе диафрагмированного волновода с бегущим типом волны
можно представить в виде:
(2)
= 𝑒𝐸 cos 𝜑,
где 𝐸 – амплитуда основной гармоники, 𝜑 – фаза частицы относительно
волны, 𝑒 – заряд частицы. Если 𝑣 = 𝑐𝛽 (с – скорость света в вакууме)– фазовая скорость волны, то фазу 𝜑 можно найти как
𝑑𝑧
−𝑡
𝑣ф
𝜑=𝜔
Так как
= 𝑣 = 𝑐𝛽 – скорость частицы, то после преобразований полу-
чим:
=
. (3)
−
Произведя замену переменных: 𝜉 = , 𝛾 =
и𝐺=
, где 𝜆 – длина
волны ускоряющей структуры, а 𝑚 – масса покоя электрона, преобразуем
уравнения (2) и (3) к виду:
= 𝐺 cos 𝜑,
= 2𝜋
−
.
В случае, если фазовая скорость волны равна скорости света, уравнения
приобретают окончательный вид
(4)
= 𝐺 cos 𝜑.
= 2𝜋 1 −
,
(5)
Выражения (4)-(5) представляют собой уравнения движения для продольной фазовой плоскости. Из (4) и (5) можно получить гамильтониан:
14
𝐻 = −𝐺 sin 𝜑 + 2𝜋 ∫ 1 −
𝑑𝛾 = −𝐺 sin 𝜑 + 2𝜋𝛾(1 − 𝛽) (6)
где 𝛾 – релятивистский фактор частицы. Решения системы (4)-(5) представляют собой кривые на фазовой плоскости (𝛾, 𝜑), вдоль которых может двигаться частица, причем ее гамильтониан сохраняется и в любой момент времени равен
𝐻 = −𝐺 sin 𝜑 + 2𝜋𝛾 (1 − 𝛽 ), (7)
где 𝜑 – начальная фаза волны, в которую попадает частица при инжекции, 𝛾
и 𝛽 – соответствующие кинетические параметры в момент инжекции.
3.2 Условие захвата в режим ускорения
Определим, при каких условиях частица может быть захвачена волной в
ускорение и непрерывно увеличивать свою энергию, пока находится в поле
этой волны. При 𝛾 → ∞, 𝛽 = 1 − (2𝛾 ) ⇒ 𝛾(1 − 𝛽) → 0. Отсюда, подставляя в выражение для гамильтониана (6), находим:
sin 𝜑
→
→ sin 𝜑 = −
(8)
где 𝜑 - фаза, которая будет являться асимптотой для фазовой траектории в
случае непрерывного ускорения частицы вдоль этой траектории. Так как
|sin 𝜑 | ≤ 1, то достигнуть сколь угодно большой энергии могут только те частицы, у которых 𝐻 < 𝐺. Частицы, для которых 𝐻 > 𝐺, всегда имеют конечную энергию, а соответствующая фазовая траектория будет осциллировать.
Требование 𝐻 < 𝐺 равносильно
2𝜋𝛾
1−
> 2𝜋𝛾 − 𝐺(1 + sin 𝜑 ).
(9)
Это неравенство является иррациональным, при его решении может быть
два случая.
Первый случай заключается в том, что выражение в правой части неравенства (9) больше нуля, тогда имеем право возводить неравенство в квадрат:
4𝜋 𝛾
2𝜋𝛾 − 𝐺(1 + sin 𝜑 ) > 0
1
1−
> 2𝜋𝛾 − 𝐺(1 + sin 𝜑 )
𝛾
После некоторых преобразований получаем
15
𝛾 >
+
)
(
𝐺<(
(
)
.
(10)
)
Как видно, можно считать, что данный сценарий можно применить к
условно малым полям. В зависимости от параметра 𝐺 при заданных фазах инжекции можно определить минимальную начальную энергию, которая необходима, чтобы частицы захватились в режим стационарного ускорения. При
этом параметр 𝐺 не может быть любым, поскольку для него также действует
ограничение, связанное с энергией и фазами инжекции. То есть возникает множество параметров [𝛾 , 𝐺], которое может удовлетворять (10) и зависит от фаз
инжекции.
Второй случай возникает, если в (9) в правой части неравенства выражение меньше нуля, тогда решением неравенства является вся область допустимых значений для правой части, то есть
2𝜋𝛾 − 𝐺(1 + sin 𝜑 ) < 0
1
1−
≥0
𝛾
или
𝐺>(
𝛾 ≥1
).
(11)
Как видно, это можно отнести к случаю условно больших полей. Данное
условие говорит о существовании множества значений амплитуд ускоряющих
полей и начальных фаз, для которых возможно ускорение частиц с любыми
начальными энергиями. Само множество амплитуд полей ускоряющей структуры зависит от начальных параметров инжектируемых частиц и не может
быть произвольным.
Захват в режим стационарного ускорения говорит о том, что через неопределенный промежуток времени частица попадет в такие условия, что ее энергия будет только увеличиваться. Заметим, что критерии (10) и (11) не говорят
ничего о том, как будет происходить сам захват в режим ускорения. Другими
словами, нельзя сказать, что будет с частицами в промежуток времени между
инжекцией и стационарным ускорением. Очевидно, что одним из возможных
и неблагоприятных сценариев является изменение направления движения частицами. При этом, данное обстоятельство не противоречит попаданию части-
16
цей в режим стационарного ускорения. Она может совершить несколько продольных осцилляций и в итоге попасть в нужную фазу для последующего увеличения энергии. То есть, с точки зрения продольной динамики, критерии захвата будут выполнены, но в реальном ускорителе, скорее всего, такие частицы будут потеряны, особенно, при наличии внешнего магнитного поля.
С помощью фазовой траектории – зависимости гамма-фактора частицы
от фазы – рассмотрим, как изменяется направление движения частицы. Уравнение фазовой траектории можно получить из выражения для гамильтониана
(6):
𝛾(𝜑) =
(
)
(
)
. (12)
Найдем условие, при котором частица может остановиться в некоторой фазе
𝜑 , потеряв всю энергию. В этом случае в (12) необходимо, чтобы 𝛾(𝜑 ) обращалось в 1:
1=
(𝐻 + 𝐺 sin 𝜑 ) + 4𝜋
4𝜋(𝐻 + 𝐺 sin 𝜑 )
В итоге, имеем
. (13)
sin 𝜑 =
Из (13) следует, что фазы остановки 𝜑 не существуют тогда и только тогда,
когда выполняется одно из неравенств совокупности
2𝜋 − 𝐻
< −1
𝐺
.
2𝜋 − 𝐻
>1
𝐺
Подставив выражение для гамильтониана (7) через параметры инжекции 𝛾 и
𝛽 , можно переписать предыдущую совокупность неравенств как
. (14)
Заметим, что 2𝜋 1 − 𝛾 (1 − 𝛽 ) > 0 при любых 𝛾 и 𝛽 , следовательно, второе
неравенство системы (14) описывает область отрицательных амплитуд полей,
а потому может быть опущено. Таким образом, неравенство
𝐺<
(
)
(15)
17
определяет область, в которой частица будет двигаться без остановки.
В итоге можно заключить, что существует несколько сценариев инжекции частиц в ускоряющую структуру с бегущей волной. Первый сценарий описывается системой (10). В этом случае частица попадает в режим ускорения,
но ее начальная энергия должна быть выше определенного значения. Второй
случай связан с условием (11). Здесь частица также попадает в режим ускорения, ее начальная энергия может быть любой, но поля структуры должны быть
выбраны с учетом данной энергии. Стоит заметить, что в обоих случаях
начальная энергия частиц определяет ширину области захвата, то есть диапазон начальных допустимых фаз и амплитуд полей. Кроме данных двух случаев, существует критерий (15), которые не противоречит захвату частиц в режим ускорения, но говорит о том, что перед этим частица может потерять полностью свою начальную энергию, после чего она начнет ускоряться в обратном направлении. Если при этом выполняются условия захвата (10) и (11), то
через какое-то время частица снова начнет терять свою энергию, изменит
направление движения во второй раз и попадет в область стационарного ускорения. Таким образом, возможны две точки изменения направления движения
частицы.
Вышеперечисленные случаи (10), (11) и (15) отображены на Рис. 7 при
энергии инжекции частиц 𝑊 = 0.2 МэВ с учетом связи параметра 𝐺 с амплитудой ускоряющего поля. Закрашенные области пронумерованы соответственно четырем случаям. Видно, что при увеличении амплитуды электрического поля область фаз, когда частица может захватиться в режим ускорения
и не изменить направление своего продольного импульса, уменьшается. То
есть большие амплитуды ускоряющих полей не всегда выгодны для ускорения
частиц с нерелятивистскими начальными энергиями.
18
Рис. 7. Закрашенные области соответствуют различным случаям 1-4 захвата в
ускорение: область 1 – захват в ускорение при энергии инжекции 0.2 МэВ, нет точек полной остановки, область 2 – захват в ускорение при энергии инжекции 0.2 МэВ, существуют точки полной остановки, область 3 – захват в ускорение при любой энергии инжекции, нет точек полной остановки, область 4 – захват в ускорение при любой энергии
инжекции, существуют точки полной остановки.
3.3 Теоретические оценки захвата частиц в ускорение полем бегущей
волны
С помощью упомянутых выше условий захвата определим, какой процент
частиц попадает в область захвата – с учетом смены направления движения и
без него.
Используя условия захвата и отражения, разобьём всю область фаз [−𝜋/
2; 3𝜋/2] на несколько областей в зависимости от того, как ведет себя частица
при инжекции.
На границах области захвата выполняется равенство
𝐺 = −𝐺 ∙ sin 𝜑 + 2𝜋𝛾 (1 − 𝛽),
откуда находим значение синуса фаз инжекции, соответствующих границам
области захвата:
sin 𝜑 =
2𝜋𝛾 (1 − 𝛽)
− 1.
𝐺
Таким образом, границами области захвата будут фазы инжекции, равные (левая и правая, соответственно)
asin
(
)
− 1 и 𝜋 − asin
(
)
−1 .
19
Отсюда может быть найдена доля частиц, которая может быть захвачена в
ускоряющее поле бегущей волны бесконечной протяженности:
2𝜋𝛾 (1 − 𝛽)
𝐺
2𝜋
𝜋 − 2 asin 1 −
Рис. 8. Теоретическая зависимость доли захваченных частиц (%) от амплитуды
поля.
Далее определим границы области фаз инжекции, в которой частица
сразу начинает ускоряться (то есть не останавливается и не меняет направление движения). На границах выполнено равенство
𝐺=
2𝜋 1 − 𝛾 (1 − 𝛽 )
.
1 − sin 𝜑
Отсюда находим значение синуса фаз, лежащих на границе области захвата
без остановки:
sin 𝜑 = 1 −
2𝜋 1 − 𝛾 (1 − 𝛽 )
.
𝐺
Отсюда же видно, что в полях с
𝐺 < 𝜋(1 − 𝛾 (1 − 𝛽 ))
отражения нет при любой фазе инжекции (для энергии инжекции 200 кэВ это
примерно 6 МВ/м)
20
Рис. 9. Теоретическая зависимость доли частиц (%), захваченных и двигающихся
без остановки от амплитуды поля.
Таким образом, границами области захвата с последующим движением без
остановки будут фазы инжекции, равные (левая и правая, соответственно)
asin 1 −
(
)
и 𝜋 − asin 1 −
(
)
Представленные зависимости характерны для инжекции непрерывного
пучка. В нашем же случае пучок является недогруппированным, поэтому, зная
амплитуду ускоряющего поля, можно качественно оценить долю захваченных
частиц в ускорение. Для первой ускоряющей структуры характерная амплитуда ускоряющего поля составляет 35 МВ⁄м. Поэтому, учитывая все вышеизложенные рассуждения, можно ожидать, что для непрерывного пучка захват
частиц в ускорение, учитывая вторичные сгустки, будет больше 70% (см. Рис.
8) , а захват в ускорение для первого сгустка будет не хуже, чем 53% (см. Рис.
9).
4 Измерение мощностей клистронов
4.1 Измерение с помощью диодного детектора мощности
Детекторы мощности СВЧ-сигналов (Рис. 10) предназначены для выделения огибающей высокочастотных сигналов и применяются для измерения
уровня мощности в заданном диапазоне частот или на определенной рабочей
частоте.
21
Рис. 10. Упрощенная схема диодного детектора мощности.
При подключении внешней нагрузки конденсатор детектора совместно с
импедансом нагрузки образуют фильтр низких частот, необходимый для выделения низкочастотного продетектированного напряжения.
Фильтр — устройство для выделения желательных компонентов спектра
электрического сигнала и/или подавления нежелательных. Частотные спектры
всех клистронов были сняты при помощи анализатора спектра. На Рис. 11
представлен спектр клистрона 3, откуда видно, что присутствует вторая гармоника на удвоенной рабочей частоте, амплитуда которой отличается от основной на 10 дБ. Для подавления этой гармоники при проведении измерений
был использован фильтр на 3 ГГц.
Рис. 11. Спектр клистрона 3.
Так как целью работы является моделирование продольной динамики
пучка ВЭПП-5 на этапе его инжекции в ускоряющую структуру, то, для того
чтобы применить теоретические выкладки, необходимо получить значения
мощностей, которые ответвляются в ускоряющие секции. Схема измерений
при помощи диодных детекторов мощности приведена на Рис. 12. Для ослабления сигнала падающей волны клистрона с целью его измерения использу-
22
ется ответвитель на 70 дБ, а также дополнительные аттенюаторы с ослаблением порядка 30 дБ. Суммарная аттенюация для каждого клистрона указана в
подписях к Рис. 13-Рис. 16.
Клистрон
Ответвитель
+
аттенюаторы
Фильтр
Диодный
детектор
мощности
Рис. 12. Схема измерения мощности.
Измеренные сигналы с диодных детекторов мощности приведены на
Рис. 13–Рис. 16. Сигнал осциллографа был откалиброван при помощи СВЧ
генератора мощности Agilent N583A. Указанная мощность клистронов получена по среднему значению сигнала в области «полочки».
Рис. 13. Сигнал с первого клистрона. Суммарная аттенюация 99.02 дБ.
𝑃 = 43.45МВт.
Рис. 14. Сигнал со второго клистрона.
Суммарная аттенюация 101.31 дБ.
𝑃 = 52.84 МВт.
Рис. 15. Сигнал с третьего клистрона.
Суммарная аттенюация 101.55 дБ.
𝑃 = 44.22 МВт.
Рис. 16. Сигнал с четвертого клистрона.
Суммарная аттенюация 101.94 дБ.
𝑃 = 42.46 МВт.
Так как пучок является нерелятивистским только в первой ускоряющей
секции, то в рамках изучения продольной динамики электронного пучка особый интерес представляет мощность, которая ответвляется в первую ускоряющую секцию. Значение этой мощности может быть определено на основе измерения мощности клистрона 3.
23
С учетом параметров системы SLED, а также выходной мощности клистрона 3, можно определить, какая мощность поступает в первую ускоряющую секцию. Для моделирования динамики используется оценка мощности
клистрона 3, полученная на основе измерений при помощи диодного детектора мощности.
5 Моделирование продольной динамики пучка
Два аспекта динамики пучка – продольная и поперечная динамика – влияют на его прохождение в ускорительном тракте. В поперечной динамике на
токопрохождение влияют магнитные элементы ускорителя (распределение
полей магнитных элементов и ошибки поля, поперечная перефокусировка и
недофокусировка пучка). В продольной динамике на токопрохождение влияют группировка пучка, распределение поля в ускоряющей структуре. В свою
очередь, распределение поля в ускоряющей структуре зависит от мощности
клистрона, а также от синхронизации инжекции пучка с поступлением мощности от системы умножения мощности SLED. В рамках данной дипломной
работы рассмотрена только продольная динамика пучка в предположении оптимизированной поперечной динамики, то есть параметры магнитной системы
подобраны так, чтобы токопрохождение было максимально возможным.
5.1 Параметры пучка
Для моделирования динамики пучка использовался программный пакет
ASTRA [6]. В качестве параметров пучка использовались измеренные параметры пучка инжекционного комплекса (Таблица 6, Рис. 2). Рабочая частота
ускоряющего поля в ускоряющей секции 𝑓 = 2856 МГц. Следовательно, пучок полной измеренной длительности около 7 нс является для такой структуры
непрерывным. Поэтому для моделирования можно принять длительность
пучка равной периоду ускоряющей секции, то есть 350 пс. Соответствующее
продольное фазовое пространство приведено на Рис. 17.
Таблица 6. Параметры пучка для моделирования.
Ток пучка, А
Длительность пучка, пс
5.25
350
24
Рис. 17. Продольное фазовое пространство пучка в начале моделирования.
5.2 Расчет энергии частицы при пролете ускоряющей секции с учетом
SLED
В инжекционном комплексе мощность в ускоряющие секции поступает в
виде импульса длительностью 0.5 мкс. Ускорение пучка производится в тот
момент, когда СВЧ мощность достигает нагрузки, ускоряющая структура заполнена, и набор энергии пучком максимален. Поэтому важно оценить набор
энергии пучком.
Распределение ускоряющего поля в начале ускоряющей структуры относительно момента переворота фазы SLED представлено на Рис. 18. С помощью
этой зависимости можно оценить, какое максимальное ускоряющее поле видит пучок в ускоряющей структуре. Для измеренной мощности клистрона 3,
равной 44.22 МВт, максимальная амплитуда ускоряющего поля в первой ускоряющей секции составляет около 35 МВ/м.
25
Рис. 18. Амплитуда ускоряющего поля в начале ускоряющей секции относительно
момента переворота фазы.
Пучок при влете в ускоряющую структуру является недогруппированным, поэтому до момента, пока 𝛽 пучка не станет равной единице, пучок будет
догруппировываться. Далее, когда 𝛽 = 1, пучок «вмораживается» в ускоряющее поле и видит одну амплитуду до конца ускоряющей структуры. Данный
процесс происходит в первых ячейках ускоряющей структуры, поэтому, считая момент инжекции пучка оптимальным, можно принять, что пучок пролетает от начала до конца ускоряющей структуры при максимальной амплитуде
ускоряющего поля, которую выдает SLED.
Учитывая распределение поля в ускоряющей структуры относительно
момента переворота фазы, можно оценить, какую энергию наберет частица
при пролете ускоряющей структуры (Рис. 19). Максимальный набор энергии,
который соответствует оптимальному моменту инжекции пролету частицы
при максимальной амплитуде ускоряющего поля, равен 70 МэВ.
26
Рис. 19. Набор энергии частицы при пролете ускоряющей секции относительно момента переворота фазы.
Следовательно, в качестве оценочных параметров при проведении моделирования можно принять амплитуду ускоряющего поля равной 35 МВ⁄м,
при этом набираемая энергия в первой ускоряющей секции составит около 70
МэВ.
5.3 Варьирование мощности группирователя
Для наилучшего захвата в ускорение необходимо определить оптимальную амплитуду электрического поля группирователя. Она должна быть такой,
чтобы пучок после пролета группирователя был недогруппирован, так как
процесс группировки продолжается после группирователя, а также в нескольких ячейках ускоряющей структуры, пока энергия пучка не станет достаточно
релятивистской.
Для подбора оптимальной амплитуды проводились моделирования для
различных мощностей группирователя. Процесс группировки для амплитуды
поля 1.5 МВ/м представлен на Рис. 20. После окончательного завершения
группировки пучка в ускоряющей секции получается четко различимый сгусток, в котором содержится ≈ 90% частиц, а также «хвост» пучка, состоящий
из тех частиц, которые сначала попали в замедляющую фазу, а потом были
ускорены. Динамика таких частиц подробно описана в разделе 3. Здесь и далее
под полным токопрохождением пучка будем понимать процент частиц, попавших при ускорении в первый ускоренный сгусток, а оценивать его будем на
выходе первой ускоряющей секции.
27
Рис. 20. Процесс группировки пучка для группирователя с амплитудой электрического
поля 1.5 МВ/м. Координата опорной частицы указана вверху диаграммы, за начало отсчета принято положение катода.
Чтобы определить оптимальную амплитуду, необходимо сравнить токопрохождение и количество частиц (Рис. 21-Рис. 22), которое попадает в энергетический разброс 1.5%, что соответствует акцептансу поворотного магнита
линейного ускорителя.
28
Рис. 21. Зависимость токопрохождения от амплитуды группирователя.
Рис. 22. Зависимость количества частиц, которые попали в 1.5% энергетический
разброс, от амплитуды группирователя
Из Рис. 21-Рис. 22 видно, что с увеличением амплитуды полное токопрохождение уменьшается. Также уменьшается с ростом амплитуды группирователя количество частиц, которое попадает в энергетический разброс 1.5%.
В инжекционном комплексе начальная энергия пучка после пушки составляет 200 кэВ. Поэтому данный вид графиков (Рис. 21-Рис. 22) объясняется
тем, что с увеличением мощности группирователя частицы, которые попали в
замедляющую фазу, получают обратный импульс и возвращаются на катод
или теряются. При дальнейшем увеличении мощности группирователя, его
действие на частицы пучка становится схожим с действием ускоряющей
29
структуры, а при ускорении без предварительной группировки доля прошедших частиц будет стремиться к 50%. Исходя из проведенных моделирований
можно установить, что оптимальная мощность группирователя составляет
около 10 кВт. Токопрохождение падает почти линейно с повышением мощности –около 5% на 10 кВт мощности.
30
6 Заключение
В рамках дипломной работы проведен анализ продольной динамики непрерывного пучка с энергией инжекции 200 кэВ в ускоряющей структуре на
основе диафрагмированного волновода с постоянным импедансом. На основе
анализа получены условия захвата непрерывного пучка в ускорение структурой на бегущей волне. Получены теоретические оценки захвата пучка в ускорение как с учетом отраженных частиц, так и без них. Исходя из теоретических
оценок захвата, установлено, что для наилучшего захвата в ускорение непрерывного пучка необходимо инжектировать его в поля с небольшой амплитудой ускоряющего поля.
Проведены измерения мощности и спектров клистронов при помощи диодного детектора мощности.
С учетом измеренной мощности клистрона 3 и параметров системы умножения мощности SLED получены оценки набираемой энергии.
Построены распределения электрических и магнитных полей линейного
ускорителя ВЭПП-5 от пушки до конца первой ускоряющей секции. В данных
поля проведено моделирование динамики пучка при различных амплитудах
электрического поля группирователя. Показано влияние поля группирователя
на полное токопрохождение, а также на захват частиц с учетом необходимого
энергетического разброса 1.5%.
По теоретическим оценкам, без предварительной группировки захват в
ускорение должен составлять около 53%. В моделировании показано, что максимальное токопрохождение составляет 93% при мощности группирователя
около 10 кВт и при повышении мощности линейно спадает около 5% на 10 кВт
мощности. Однако стоит заметить, что такое высокое токопрохождение соответствует произвольному энергетическому разбросу, в то время как для прохождения поворота необходим разброс 1.5%. Поэтому в моделировании была
также определена зависимость количества частиц, которое попадает в данный
энергетический разброс от мощности группирователя. Показано, что максимальная доля частиц, которая удовлетворяет данному энергетическому разбросу, составляет около 73%.
При наличии системы калибровки системы мощностей на ВЭПП-5 данные результаты могут быть использовано для увеличения его эффективности.
31
7 Список литературы
1) В. М. Павлов. СВЧ-группирователь форинжектора ВЭПП-5. Препринт. Институт ядерной физики им. Г.И. Будкера СО РАН. 2002.
2) Farkas Z. D., Hogg H. A., Loew G. A., Wilson P. B. SLED: A Method of doubling
SLAC’s Energy // Proc. IX International Conference on High Energy Accelerators.
Stanford, California, 1974.
3) М. В. Мальцева, В. М. Павлов, С. Л. Самойлов, А. В. Андрианов, И. Л. Пивоваров. Измерения резонаторов системы умножения мощности типа SLED.
Институт ядерной физики им. Г.И. Будкера СО РАН. 2017.
4) Вольман В. И., Пименов Ю. В. Техническая электродинамика. М.: Радио и
связь, 2000.
5) M. V. Arsentyeva, K. N. Berdnikov, and A. E. Levichev, Analysis of the particle
longitudinal motion during the injection process into the disk loaded traveling wave
accelerating structure, J. Inst. 15, (2020) P10015.
6) https://www.desy.de/~mpyflo/
32
8 Приложение 1. Поля магнитных элементов.
Рис. 23. Амплитуда магнитного поля линзы на оси.
Рис. 24. Суммарная амплитуда магнитного поля токовых колец на оси.
33
Рис. 25. Амплитуда магнитного поля соленоида УС1 на оси.
34
Отзывы:
Авторизуйтесь, чтобы оставить отзыв