УДК 621.384.6 Вестник СПбГУ. Прикладная математика. Информатика... 2022. Т. 18. Вып. 4 MSC 97М50
Моделирование динамики пучка в ускорителе RFQ для проекта DARIA*
Г. Н. Кропачев1'2, Т. В. Кулевой1, А. Л. Ситников1, Е. Р. Хабибуллина1, С. В. Виноградов3,
1 Научно-исследовательский центр «Курчатовский институт», Российская Федерация, 117218, Москва, пл. Курчатова, 1
2 Объединенный институт ядерных исследований, Российская Федерация, 141980, Московская обл., Дубна, ул. Жолио-Кюри, 6
3 Московский физико-технический институт, Российская Федерация, 141701, Московская обл., Долгопрудный, Институтский пер., 9
Для цитирования: Кропачев Г. Н., Кулевой Т. В., Ситников А. Л., Хабибуллина Е, Р., Виноградов С. В. Моделирование динамики пучка в ускорителе RFQ для проекта DARIA // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2022. Т. 18. Вып. 4. С. 568-582. https://doi.org/10.21638/11701/spbul0.2022.411
Проект DARIA (источник нейтронов, предназначенный для прикладных исследований и промышленного применения) направлен на разработку компактных генераторов нейтронов для университетов, научных центров и промышленности. В рамках проекта DARIA в Научно-исследовательском центре «Курчатовский институт» разрабатывается импульсный протонный ускоритель с пространственно-однородной квадрупольной фокусировкой (RFQ) с энергией 3 МэВ, рабочей частотой 162.5 МГц и током 100 мА. В статье описаны основные аспекты проектирования динамики пучка в RFQ, а именно: выбор параметров RFQ, оптимизация транспорта и трансмиссии пучка. Приведены результаты моделирования динамики пучка в RFQ.
Ключевые слова: динамика пучка, RFQ, транспорт, трансмиссия.
1. Введение. В рамках проекта DARIA [1] в Научно-исследовательском центре «Курчатовский институт» (ККТЭФ) разрабатывается импульсный линейный ускоритель с энергией 13 МэВ, резонансной частотой 162.5 МГц и током пучка 100 мА, предназначенный для создания компактного нейтронного генератора. Предполагается, что линейный ускоритель будет состоять из двух секций: ускорителя с пространственно-однородной квадрупольной фокусировкой (RFQ) и ускорителя с трубками дрейфа (DTL), между которыми будет осуществляться 6-мерное согласование пучка. В работе представлены результаты выбора параметров RFQ, а также численного моделирования динамики частиц. Приведем основные требования к RFQ :
* Работа выполнена при финансовой поддержке Министерства науки и высшего образования Российской Федерации в рамках Соглашения № 075-15-2022-830 от 27 мая 2022 г. (продолжение Соглашения № 075-15-2021-1358 от 12 октября 2021 г.).
© Санкт-Петербургский государственный университет, 2022
Рабочая частота, МГц............................
А/2...............................................
Ток инжекции, мА................................
Входной нормализованный эмиттанс, ^-мм- мрад
Энергия пучка, МэВ..............................
Максимальная напряженность поля, кВ/см......
Отношение аксептанса к эмиттансу..............
100 3.0 3.0 250 3
162.5 1
Транспорт, % ... Трансмиссия, %
>95 >95
Рабочая резонансная частота 162.5 МГц определялась требованием компактности линейного ускорителя. Ток протонов 100 мА выбран исходя из условия создания заданного уровня нейтронного потока на мишени в реакции Ве (р, п). В качестве инжектора в данном проекте планируется использовать газодинамический источник ионов на основе электронно-циклотронного резонанса (ЭЦР) [2]. Этот источник может обеспечивать требуемый ток с нормализованным эмиттансом не более 2 эт-мм-мрад [3].
С учетом роста эмиттанса в канале согласования пучка с RFQ (low energy transport channel — LEBT) был выбран входной нормализованный эмиттанс £ninp = Зэт-мм-мрад. Выходная энергия протонов 3 МэВ соответствует согласованному пучку для перезахвата в последующей секции DTL без потерь. Величина максимальной напряженности электрического поля па поверхности электрода Esmax выбирается из соотношения
Esmax = Екр - Kp,
где EKp — критерий Килпатрика: напряженность поля на поверхности электрода,
Kp
фициент Килпатрика, учитывающий режим работы ускорителя. Для импульсного режима, как правило, Kp = 1.8. Тогда для проекта DARIA Esmax = 250 кВ/см.
Транспорт и трансмиссия — важные характеристики ускорителя. Транспорт — это отношение количества частиц, выходящих из ускорителя, к количеству инжектируемых частиц; трансмиссия — отношение количества ускоренных частиц на выходе линейного ускорителя к количеству инжектируемых частиц. Требования к RFQ: транспорт и трансмиссия — более 95 %, т. е. минимизация потерь частиц в RFQ. Цель состоит в том, чтобы свести к минимуму образование нейтронов и последующую активацию RFQ.
Как известно, нормализованный аксептанс канала — максимально возможный нормализованный эмиттанс согласованного пучка, транспортируемого в канале. Отношение нормализованного аксептанса Vk к нормализованному входному эмиттансу £ninP определяет возможную трансмиссию пучка в линейном ускорителе в зависимости от действия кулоновских сил. В соответствии с формулами И. М. Капчинского [5], отношение Vk/eninp = 3 является минимальным, при котором обеспечивается требуемое значение трансмиссии пучка. Поэтому нормализованный аксептанс определяется как Vk = 9эт-мм-мрад.
2. Выбор основных параметров RFQ. В секции RFQ непрерывный пучок от источника преобразуется в отдельные сгустки. Таким образом, RFQ подготавливает пучок к его дальнейшему захвату на следующем участке в секции DTL.
При выборе параметров RFQ применялся метод, основанный на законах масштабирования в RFQ, который был описан в [6]. В статьях [7, 8] также использовался этот метод. Кроме того, принимались во внимание результаты оптимизации поперечного и продольного движений заряженных частиц в RFQ [9, 10]. Соответствующие результаты представлены па рис. 1, где показаны зависимости параметров RFQ: m, ¡j>t, До от изменения значения нормализование го аксептанса Vk. Зде с ь m — отношение максимального расстояния от оси до электрода к минимальному расстоянию (модуляция электрода), fit — набег фазы поперечных колебаний, До — среднее расстояние от оси до электродов (средний радиус).
Рис. 1. Параметры КРС^: т (1), цг {2), Яо (3) в зависимости от величины нормализованного аксептанса
Согласно рис. 1, для обеспечения требуемого аксептанса канала У = 9 п-мм х мрад) при заданной рабочей частоте 162.5 МГц набег фазы поперечных колебаний равен ~ 43°, средний радиус Я0 = 9 мм и максимальное значение модуляции т — 1.95. Тогда радиус апертуры канала определяется из соотношения
а = --—Яо.
т + 1
В работе [11] было показано, что фокусирующее поле электродов ИРС^ имеет наилучшую линейность при р ~ 1.14Я0, а напряженность поля на поверхности электродов минимальна при р ~ 0.7Я0, где р — радиус кривизны электрода, Я0 — средний радиус канала. В качестве компромисса отношение радиуса кривизны электрода к среднему радиусу р/Я0 остается постоянным на уровне 0.8 для уменьшения напряженности поля на поверхности. Основные расчетные параметры ИРС^ следующие:
Напряжение инжекции, кВ.................................................60
Амплитуда напряжения между смежными электродами, кВ...............173
Максимальный коэффициент модуляции электродов.......................1.95
Минимальный радиус апертуры, мм........................................6.1
Среднее расстояние от оси до электрода, мм...............................9
Радиус кривизны электрода, мм............................................7.2
Равновесная фаза, град.....................................................-90^—40
Нормализованный аксептанс канала, п-мм- мрад...........................9
Параметры RFQ выбраны исходя из условия приблизительного равенства набегов фаз продольных ¡iz и поперечных i колебаний ионов на периоде фокусировки:
it = (1.05 ^ 1.10|•
Конструкция RFQ разделена на участки: 1) формирователя сгустков (Shaper), который начинается с раструба: поперечного согласователя (Radial Matcher Section — RMS); 2) группирующего участка с ускорением (Gentle Buncher); 3) регулярного ускорения (Accelerator). Секция RMS состоит из шести ячеек и позволяет уменьшить на-
клон огибающей пучка на входе RFQ. Таким образом, снижаются требования к фокусирующим элементам LEBT. На первом участке начинается слабая модуляция электродов для подготовки пучка к группировке. В этой части синхронная фаза находится в пределах —90° ^—85° и выполняется поперечное и продольное согласование пучка в RFQ. Затем второй участок группирует пучок с минимальной фазой около —40°. Особенностью такой схемы является то, что нормализованный аксептанс канала определяется на стыке второго и третьего участков.
На участке регулярного ускорения средний радиус апертуры, модуляция электрода и синхронная фаза поддерживаются постоянными: До = const, a = const, m = const, p = const. В этой части RFQ частицы получают максимальный прирост энергии.
Использование адиабатического изменения параметров RFQ вдоль группирующего участка приводит к постепенному уменьшению фазовой длины сепаратрисы. При этом с ростом скорости частиц сгустки раздвигаются, сохраняя неизменными геометрические размеры, а следовательно, и плотность пространственного заряда. Такой режим оптимален для получения большого значения аксептанса и высокого предельного тока [12, с. 130-143]. На рис. 2 показаны основные параметры RFQ в зависимости от длины: энергия W, равновесная фаза p, фактор пролетного времени Т, коэффициент модуляции m, средний радиус До, эффективность фокусировки ж и нормализованный аксептанс Vk вдоль RFQ. Здесь эффективность фокусировки ж
Т
идеальных электродов в RFQ [11]:
m п m2 — 1
Т =--,
4 m2Io(ka) + Io(mka)'
4
ж = 1 - -Т10(ка), п
где Io — модифицированная функция Бесселя нулевого порядка; k = 2n/ß\ — волновое число.
Энергия: прирост энергии на метр на участке регулярного ускорения максимален и составляет 0.94 МэВ/м; на формирователе выходная энергия составляет 0.070 МэВ, выходная энергия группирователя — 0.7 МэВ.
Равновесная фаза: в начале формирователя она сохраняется постоянной и равной —90° —85°
—40°
тическое группирование частиц. На участке регулярного ускорения фаза постоянная — 40°
m
теле он меняется незначительно из-за небольшого изменения равновесной фазы; на группирователе модуляция меняется с учетом адиабатичности параметров; на участке регулярного ускорения величина модуляции постоянна и равна 1.95 для получения максимальной энергии на выходе.
Средний радиус: на участке RMS средний радиус изменяется по синусоидальному закону, обеспечивая радиальное согласование пучка с RFQ, после раструба он постоянен и равен 9 мм.
Фактор пролетного времени (эффективность ускорения): на участке RMS при m Т Т
Т
3.5
3.0
о
^ 2.0 §
^ 1.5 £
Ч- 1.0
0.5
0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0
Длина, м
Рис. 2. Основные параметры RFQ вдоль всей структуры: W (МэВ), —^ (рад), T, m, Ro (см), ж, Vk (п-мм-мрад)
Эффективность фокусировки: на участке RMS при отсутствии ускорения эффективность фокусировки имеет максимальное значение ж = 1; с увеличением модуляции она уменьшается до 0.4 и сохраняется вдоль всего участка регулярного ускорения.
Нормализованный аксептанс: как было сказано выше, особенностью данной конструкции RFQ является то, что нормализованный аксептанс определяется на стыке участков группировки и регулярного ускорения частиц (рис. 2), и он равен Vk = 9п -мм-мрад.
3. Численное моделирование динамики пучка в RFQ. Численное моделирование динамики пучка в RFQ проводилось методом крупных частиц с использованием программ Trace Win [13] и LIDOS [14]. Графическое представление результатов моделирования осуществлялось в программе Plot Win [15]. Потенциал кулоновского поля рассчитывался на сетке 32 х 32 х 32 для одного периода модуляции в декартовой системе координат. Число шагов интегрирования уравнений движения по периоду модуляции было принято равным 100. Число частиц в полном пучке — 105. На входе в RFQ использовались усеченное (±3а) гауссово распределение в поперечном четырехмерном фазовом пространстве, а также равномерное распределение в продольном фазовом пространстве. Обычно в RFQ применяется четырехмерное согласование пучка (стандартное согласование). Опишем начальные согласованные условия на входе в RFQ для данного случая. Частицы равномерно распределены на продольной фазовой плоскости с фазовой длиной 360° и разбросом по энергии ±1 %. Эти условия типичны для большинства источников ионов. В радиальном направлении инжектируется аксиальносимметричный сходящийся пучок. На рис. 3 показан фазовый портрет согласованного пучка на входе RFQ.
Согласованный пучок способствует уменьшению роста эмиттанса. Согласование пучка с ускорителем характеризуется взаимным расположением эллипсов аксептан-са канала и эмиттанса пучка. На рис. 4 приведены типичные эллипсы аксептанса канала (темный) и согласованный эмиттанс пучка (светлый) на входе RFQ. На нем видно, что эмиттанс инжектируемого пучка и аксептанс канала подобны, т. е. согла-
Рис. 3. Фазовый портрет пучка на входе ИРС^ (условие согласования) а — плоскость ^ —Xб — плоскость У — У'; в — плоскость Фаза — Энергия;
г — плоскость X — У.
сованы. Численное моделирование показало, что отношение аксептанса к эмиттансу приблизительно равно 3, что соответствует требованиям (см. с. 567-568).
Согласно моделированию динамики пучка, длина всей структуры равна
5.446 м, формирователя сгустков — 1.224 м, группирователя — 1.454 м, участка регулярного ускорения — 2.766 м. Напряжение инжекции составляет 60 кВ, что позволяет обеспечить инжекцию требуемого тока при пониженном влиянии кулоновских сил.
На рис. 5 приведен профиль пучка КРС^ для тока инжекции 100 мА. На нем видно, что из 100 000 инжектированных макрочастиц 99 694 достигли выхода ускорителя, т. е. транспорт составляет 99.7 %, что соответствует требованиям, сформулированным на с. 567-568.
Следует отметить, что:
— максимальный размер огибающих в поперечных плоскостях практически не увеличивается, что свидетельствует о хорошем согласовании пучка с каналом в радиальном направлении;
— наиболее плотная область пучка (его ядро) занимает около 0.75 площади апер-
Рис. 4- Поперечное согласование нормализованного аксептанса Ук (темный) и нормализованного входного эмиттанса еПгПр (светлый) на входе ИРС^
Рис. 5. Профиль пучка в КРС^: ^ У и фаза по всей длине структуры
туры канала, где действие сил электрического поля, как правило, линейно. Это дает минимальный вклад в рост фазового объема;
— происходят постепенное уменьшение фазовой длины сгустка и его группировка вокруг равновесной частицы;
— когда на формирователе начинают образовываться сгустки (т. е. раздвигаются и ускоряются) под действием продольных кулоновских сил, возникают фазовые колебания центра пучка. По мере увеличения энергии частиц и расстояния между сгустками эти колебания затухают [12];
— из-за действия кулоиовских сил небольшое количество частиц в группирова-теле не достигает энергии, близкой к энергии равновесной частицы.
На рис. 6 показан фазовый портрет на выходе для всех частиц, которые
прошли всю структуру. Из него следует, что незначительная часть ионов (примерно 0.3 %) практически не ускоряется, т. е. имеет энергию, близкую к энергии инжек-ции. Позже эти частицы будут потеряны в магнитных линзах согласующего канала между ИРС} и БТЪ.
Рис. 6. Фазовый портрет пучка на выходе ИРС^ (общее количество частиц) а — плоскость X — Xб — плоскость У— Ув — плоскость Фаза — Энергия;
г — плоскость X — У.
На рис. 7 приведен фазовый портрет ускоренных частиц на выходе Ш^. Видно, что трансмиссия ускорителя составляет 99.5 %, что соответствует требованиям, сформулированным на с. 567-568. Отметим, что пучок хорошо сформирован по фазам: < 20° и по энергиям: | < 0.1 МэВ, что упрощает дальнейший захват пучка в секции БТЪ.
Рис. 1. Фазовые портреты пучка на выходе ИРС^ (только ускоренные частицы) а — плоскость X — Xб— плоскость У — У'; в — плоскость Фаза — Энергия;
г — плоскость X — У.
Распределение нормализованных среднеквадратичных поперечных эмиттансов вдоль длины 11Е(3 иллюстрирует рис. 8. Наблюдается значительный рост эмиттанса на формирователе и группирующем участке с ускорением. Это объясняется связью между продольными и поперечными колебаниями частиц. Поскольку на вход 11Е(3 поступает почти монохроматический пучок, возникает пульсация сгустка (действие продольной составляющей собственного поля пучка). Длина сгустка пульсирует на формирователе и группирователе. Пульсация длины сгустка приводит к пульсациям плотности объемного заряда и, как следствие, к рассогласованию пучка и росту эмиттанса. Из этого рисунка видно, что рост среднеквадратичного эмиттанса составляет 27 %. Для распределения Гаусса, когда эмиттанс содержит 99 % пучка, рост эмиттанса достигает 65 %.
0.35 -I--—¡-I--Uj------
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 Длина, м
Рис. 8. Нормализованный поперечный эмиттанс пучка вдоль ускорителя RFQ 1 — плоскость X — X2 — плоскость Y — Y
Отметим, что выбор основных параметров RFQ был направлен на получение максимальных значений транспорта и трансмиссии при минимальном увеличении эмиттанса. Таким образом, после выбора параметров и последующего численного моделирования параметры RFQ можно при необходимости корректировать до тех пор, пока не будут достигнуты требуемые результаты моделирования.
На группирующем участке с ускорением наблюдаются небольшие потери частиц. Это показано на рис. 9, который иллюстрирует распределение транспорта и трансмиссии вдоль RFQ. В ускоряющей части RFQ практически отсутствуют как потери (рис. 9), так и рост эмиттанса (см. рис. 8).
Выше были указаны величины транспорта и трансмиссии, которые характеризуют динамику пучка в RFQ. В этом случае eninv = 3п-мм-мрад было выбрано в качестве входного значения нормализованного поперечного эмиттанса для 99 % пучка с током 100 мА. Однако при формировании пучка, выходящего из источника ионов, увеличение радиального эмиттанса может происходить из-за кулоновских сил и нелинейных полей фокусирующих элементов канала LEBT. Оценим возможности разработанного канала RFQ, когда нормализованный эмиттанс входного пучка будет выше расчетного при фиксированном токе инжекции 100 мА.
На рис. 10 видно, что даже при значительном увеличении расчетного значения входного эмиттанса в 2.5 раза транспорт канала снижается до 98 %, а трансмиссия — до 97 %. Это характеризует разработанный RFQ с нормализованным аксептансом канала Vk = 9п-мм-мрад, как устойчивый к значению входного эмиттанса в широком диапазоне его изменения.
Рассмотрим предельные возможности разработанного ускорителя RFQ по току ускоренных частиц на выходе структуры при фиксированном значении фазовой плотности (отношение входного тока к входному эмиттансу сохраняется). На рис. 11 показаны зависимости тока ускоренных частиц от тока инжекции. Как видно на нем, максимально достижимый ток составляет 300 мА. Однако увеличение тока ускорен-
ных частиц сопровождается уменьшением коэффициента захвата (рис. 11). Поскольку требуемый выходной ток составляет 100 мА, разработанный имеет трехкратный запас по максимальному выходному току пучка.
Приведем результаты моделирования динамики пучка в Ш^:
Рабочая частота, МГц......................................................162.5
А^.........................................................................1
Ток инжекции, мА..........................................................100
Входная / выходная энергия, МэВ..........................................0.06/3.3
Амплитуда напряжения между смежными электродами, кВ...............173
Максимальная модуляция..................................................1.95
Минимальный радиус апертуры, мм........................................6.1
Средний радиус, мм........................................................9
Радиус кривизны электрода, мм............................................7.2
Максимальная напряженность поля на поверхности (Кр = 1.8), кВ/см ... .250
Равновесная фаза, град.....................................................—90ч—40
Входной нормализованный эмиттанс, п-мм- мрад..........................3.0
Нормализованный аксептанс канала, п-мм- мрад...........................9
Отношение аксептанса к эмиттансу.........................................3
Длина, м....................................................................5.446
Транспорт, %...............................................................99.7
Трансмиссия, %.............................................................99.7
Рост поперечного эмиттанса................................................1.65
100.00
99.95
^ 99.90
§ 99.85 и
| 99.80 и
§ 99.75
е-
99.70
Л
§ 99.65 ё
5 99.60
Е&
99.55 99.50
1
2
).0 0 .5 1 вНАР .0 : ЕЫ .5 2 вЕР В1Л> .0 2 *ТЪЕ *СНЕ1 .5 1 3.0 3.5 4 АССЕЬЕКА! .0 4 П<Ж .5 5. Дл! 0 5. ша, м
Рис. 9. Распределение величин транспорта (1) и трансмиссии (2) вдоль БРС^
£
о? я и о S
ее Си н
н* Си о Н и Я
es £
100.0 99.5 99.0 98.5 98.0 97.5 97.0 96.5 96.0 95.5 95.0
3 4 5 6 7 8 9
Нормализованный эмиттанс входного пучка, Эт'мм'мрад
Рис. 10. Транспорт (1) и трансмиссия ускорителя (2) в зависимости от величины нормализованного эмиттанса входного пучка (при постоянном входном токе 100 мА)
Рис. 11. Выходной ток (1) и трансмиссия ускорителя (2) в зависимости от входного тока пучка (при постоянной фазовой плотности)
4. Заключение. В рамках проекта DARIA был разработан импульсный протонный ускоритель с пространственно-однородной квадрупольной фокусировкой (RFQ) с энергией 3.3 МэВ, резонансной частотой 162.5 МГц и током 100 мА (см. с. 567-568). При токе инжекции 100 мА и нормализованном эмиттансе 3^-мм-мрад транспорт в канале RFQ составляет 99.7 %, а трансмиссия — 99.5 %. При этом поперечный рост эмиттанса — 65 % (эмиттанс содержит 99 % пучка). Отметим, что значения транспорта не менее 98 % и трансмиссии не менее 97 % находятся в широком диапазоне изменения входного эмиттанса. При этом RFQ имеет тройной резерв по максимальному достижимому току, равному 300 мА. Таким образом, разработанный канал RFQ удовлетворяет приведенным в п. 1 (см. с. 577-578) требованиям.
Литература
1. Grigoriev S., Iashina E., Pavlov K. Spin-echo small angle neutron scattering for a compact neutron source DARIA // Journal of Synch. Investigation. 2019. Vol. 13. P. 1132-1134. https://doi.org/10.1134/S1027451019060314
2. Skalyga V., Izotov /., Razin S. High current proton beams production at simple mirror ion source 37 11 Review of Scientific Instruments. 2014. Vol. 85 (2). P. 02A702-1-02A702-3. https://doi.Org/10.1063/l.4825074
3. Барабин С. В., Кропачев Г. П., Лукашин А. Ю., Кулевой Т. В., Выбин С. С., Голубее С. В., Изотов И. В., Киселева У. М., Скалыга В. АГригорьев С. В., Коваленко Н. А. Измерения эмиттанса газодинамического электронно-циклотронного резонансного источника ионов // Письма в Журн. техн. физики. 2021. Т. 47 (10). С. 7-10. https://doi.org/10.21883/PJTF.2021.10.50964.18628
4. Kapchinsky I. М. About approximations of Kilpatrick criterion // PTE. 1986. N 1. P. 33-35.
5. Капчинекий И. М. Динамика частиц в линейных резонансных ускорителях. М.: Атомиздат, 1966. 247 с.
6. Wadlinger Е. A. Scaling laws for RFQ design procedures // IEEE Transactions on Nuclear Science. 1985. Vol. 32 (5). P. 2596-2598. https://doi.org/10.1109/TNS.1985.4333992
7. Kropachev G., Kulevoy Т., Sitnikov A. Proton linac for compact neutron source DARIA // Journal of Surface Investigation: X-ray, Synchrotron and Neutron Techniques. 2019. Vol. 13 (6). P. 1126-1131. https://doi.org/10.1134/S1027451019060399
8. Kulevoy Т., Fatkullin R., Kozlov A., Kropachev G., Selesnev D., Semennikov ASitnikov A. Compact multipurpose facility — BELA // Proceedings of 29th Linear Accelerator Conference. 2018. P. 349-351. https://doi.org/10.18429/JACoW-LINAC2018-TUP0012
9. Balabanov M. Yu., Mizintseva M. AOvsyannikov D. A. Beam dynamics optimization in a linear accelerator // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2018. Т. 14. Вып. 1. С. 4-13. https://doi.org/10.21638/11701/spbul0.2018.101
10. Антропов И. В., Овсянников А. Д. Моделирование и оптимизация динамики частиц в ускорителе с пространственно-однородной квадрупольной фокусировкой // Вестник Санкт-Петербургского университета. Сер. 10. Прикладная математика. Информатика. Процессы управления. 2009. Вып. 4. С. 3-24.
11. Валабин А. И., Липкин И. М., Угаров С. Б. О выборе формы электродов в линейном ускорителе с пространственно-однородной квадрупольной фокусировкой // ВАНТ. Сер. Техника физического эксперимента. 1988. Т. 1 (36). С. 20-22.
12. Капчинский И. М. Теория линейных резонансных ускорителей. М.: Энергоиздат, 1982. 241 с.
13. Uriot D. Status of TraceWin Code 11 Proceeding of IPAC. 2015. P. 92-94. https://doi.org/10.18429/JACoW-IPAC2015-MOPWA008
14. Bondarcv B. /., Durkin A. P., Murin B. P., Nikolaishvili G. Т., Shlygin O. Yu. LIDOS — unconventional helper for linac beam designing // AIP Conference Proceedings. 1993. Vol. 297. P. 377384. https://doi.org/10.1063/1.45334
15. DACM Software. URL: http://irfu.cea.fr/Sacm/logiciels/ (дата обращения: 14.07.2022).
Статья поступила в редакцию 19 августа 2022 г.
Статья принята к печати 1 сентября 2022 г.
Контактная информация:
Кропачев Геннадий Николаевич — канд. физ.-мат. наук; kropachev@itep.ru
Кулевой Тимур Вячеславович — д-р техн. наук; kulevoy@itep.ru
Ситников Алексей Леонидович — канд. физ.-мат. наук; aleksey.sitnikov@itep.ru
Хабибуллина Екатерина Рамисовна — Ekaterina.Khabibullina@itep.ru
Виноградов Станислав Владиленович — канд. физ.-мат. наук; s.vinogradov@gmail.com
Beam dynamics of the project DARIA RFQ*
G. N. Kropachev1'2, T. V. Kulevoy1, A. L. Sitnikov1, E. R. Khabibullina1, S. V. Vinogradov3
1 National Research Centre "Kurchatov Institute", 1, pi. Kurchatova, Moscow, 117218, Russian Federation
2
141980, Russian Federation 3 Moscow Institute of Physics and Technology, 9, Institutskiy per., Dolgoprudny, Moscow Region, 141701, Russian Federation
For citation: Kropachev G. N., Kulevoy T. V., Sitnikov A. L., Khabibullina E. R., Vinogradov S. V. Beam dynamics of the project DARIA RFQ. Vestnik of Saint Petersburg University. Applied Mathematies. Computer Science. Control Processes, 2022, vol. 18, iss. 4, pp. 568-582. https://doi.org/10.21638/11701/spbul0.2022.411 (In Russian)
The project DARIA (neutron source dedicated to applied research and industrial applications) directed to the development of compact neutron generators for universities, scientific centers and industry. In framework of the DARIA project, the 3 MeV, 162.5 MHz and 100 mA pulse proton Radio Frequency Quadrupole (RFQ) is under development at National Research Centre "Kurchatov Institute". In this article the main aspects of the beam dynamics design of this RFQ are described, namely the choice of RFQ parameters, the optimization of the transport and the transmission. Results of beam dynamic simulation in RFQ are presented. Keywords: dynamic simulation, RFQ, transport, transmission.
References
1. Grigoriev S., Iashina E., Pavlov K. Spin-echo small angle neutron scattering for a compact neutron source DARIA. Journal of Synch. Investigation, 2019, vol. 13, pp. 1132-1134. https://doi.org/10.1134/S1027451019060314
2. Skalyga V., Izotov I., Razin S. High current proton beams production at simple mirror ion source 37. Review of Scientific Instrumental, 2014, vol. 85 (2), pp. 02A702-1-02A702-3.
https://doi.org/10.1063/1.4825074
3. Barabin S. V., Kropachev G. N., Lukashin A. Yu., Kulevoy T. V., Vybin S. S., Golubev S. V., Izotov I. V., Kiseleva E. M., Skalyga V. A., Grigoriev S. V., Kovalenko N. A. Izmereniya emittansa gazodinamicheskogo elektronno-tsiklotronnogo rezonansnogo istochnika ionov [Measurements of the emittance of a gas-dynamic electron-cyclotron resonant ion source]. Letters to the Journal of Technical Physics, 2021, vol. 47 (10), pp. 7-10. https://doi.org/10.21883/PJTF.2021.10.50964.18628 (In Russian)
4. Kapchinsky I. M. About approximations of Kilpatrick criterion. PTE, 1986, no. 1, pp. 33-35.
5. Kapchinsky I. M. Dinamika chastits v lineynykh rezonansnykh uskoritelyakh [Particle dynamics in linear resonant accelerators]. Moscow, Atomizdat Publ., 1966, 247 p. (In Russian)
6. Wadlinger E. A. Scaling laws for RFQ design procedures. IEEE Transactions on Nuclear Science, 1985, vol. 32 (5), pp. 2596-2598. https://doi.org/10.1109/TNS.1985.4333992
7. Kropachev G., Kulevoy T., Sitnikov A. Proton linac for compact neutron source Daria. Journal of Surface Investigation: X-ray, Synchrotron and Neutron Techniques, 2019, vol. 13 (6), pp. 1126-1131. https://doi.org/10.1134/S1027451019060399
8. Kulevoy T., Fatkullin R., Kozlov A., Kropachev G., Selesnev D., Semennikov A., Sitnikov A. Compact multipurpose facility — BELA. Proceedings of 29th Linear Accelerator Conference, 2018, pp. 349-351. https://doi.org/10.18429/JACoW-LINAC2018-TUP0012
9. Balabanov M. Yu., Mizintseva M. A., Ovsyannikov D. A. Beam dynamics optimization in a linear accelerator. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2018, vol. 14, iss. 1, pp. 4-13. https://doi.org/10.21638/11701/spbul0.2018.101
10. Antropov I. V., Ovsyannikov A. D. Modelirovaniye i optimizatsiya dinamiki chastits v uskori-tele s prostranstvenno-odnorodnoy kvadrupol'noy fokusirovkoy [Simulation and optimization of particle dynamics in accelerator with spatially uniform quadrupole focusing]. Vestnik of Saint Petersburg
* This work was financially supported by Ministry of Science and Higher Education of the Russian Federation in framework of Agreement N 075-15-2022-830 from May 27, 2022 (Prolongation of Agreement N 075-15-2021-1358 from October 12, 2021).
University. Series 10. Applied Mathematics. Computer Science. Control Processes, 2009, iss. 4, pp. 3-24. (In Russian)
11. Balabin A. I., Lipkin I. M., Ugarov S. B. O vybore formy elektrodov v lineynom uskoritele s prostranstvenno-odnorodnoy kvadrupol'noy fokusirovkoy [On choosing the electrodes from in linear accelerator with spatially uniform quadrupole focusing]. VANT. Series Technique of physical experiment, 1988, vol. 1 (36), pp. 20-22. (In Russian)
[
]
13. Uriot D. Status of TraceWin Code. Proceeding of JPAC, 2015, pp. 92-94. https://doi.org/10.18429/JACoW-IPAC2015-MOPWA008
14. Bondarev B. I., Durkin A. P., Murin B. P., Nikolaishvili G. T., Shlygin O. Yu. LIDOS — unconventional helper for linac beam designing. AIP Conference Proceedings, 1993, vol. 297, pp. 377384. https://doi.org/10.1063/1.45334
15. DACM Software. Available at: http://irfu.cea.fr/Sacm/logiciels/ (accessed: June 14, 2022).
Received: August 19, 2022.
Accepted: September 1, 2022.
Authors' information:
Gennady N. Kropachev — PhD in Physics and Mathematics; kropachev@itep.ru
Timur V. Kulevoy — Dr. Sci. in Technics; kulevoy@itep.ru
Alexei L. Sitnikov — PhD in Physics and Mathematics; aleksey.sitnikov@itep.ru
Ekaterina R. Khabibullina — Ekaterina.Khabibullina@itep.ru
Stanislav V. Vinogradov — PhD in Physics and Mathematics; s.vinogradov@gmail.com