Динамика пучка в синхротронах с цифровыми системами подавления когерентных колебаний заряженных частиц тема автореферата и диссертации по физике, 01.04.20 ВАК РФ

Жабицкий, Вячеслав Михайлович АВТОР
доктора физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Дубна МЕСТО ЗАЩИТЫ
2013 ГОД ЗАЩИТЫ
   
01.04.20 КОД ВАК РФ
Диссертация по физике на тему «Динамика пучка в синхротронах с цифровыми системами подавления когерентных колебаний заряженных частиц»
 
Автореферат диссертации на тему "Динамика пучка в синхротронах с цифровыми системами подавления когерентных колебаний заряженных частиц"

ОБЪЕДИНЕННЫЙ ИНСТИТУТ ЯДЕРНЫХ ИССЛЕДОВАНИЙ

ДИНАМИКА ПУЧКА В СИНХРОТРОНАХ С ЦИФРОВЫМИ СИСТЕМАМИ ПОДАВЛЕНИЯ КОГЕРЕНТНЫХ КОЛЕБАНИЙ ЗАРЯЖЕННЫХ ЧАСТИЦ

Специальность: 01.04.20 — физика пучков заряженных частиц и ускорительная техника

9-2013-63

На правах рукописи

ЖАБИЦКИЙ

Вячеслав Михайлович

Автореферат диссертации на соискание ученой степени доктора физико-математических наук

005531788

Дубна 2013

005531788

Работа выполнена в Лаборатории физики высоких энергий им. В.И. Векслера и A.M. Балдина Объединенного института ядерных исследований.

Официальные оппоненты: Агафонов Алексей Вениаминович, доктор физико-математических наук,

Физический институт им. П.Н. Лебедева Российской академии наук, г. Москва;

Зенкевич Павел Романович, доктор физико-математических наук,

Институт теоретической и экспериментальной физики, г. Москва;

Фещенко Александр Владимирович, доктор физико-математических наук,

Институт ядерных исследований Российской академии наук, г. Москва. Ведущая организация:

Институт физики высоких энергий, г. Протвино Московской области.

Защита состоится «-» О ^¿-^-Г._2013 г. в .//. часов на заседании диссертационного совета Д720.001.03 при Лаборатории ядерных проблем им. В.П. Дже-лепова Объединенного института ядерных исследований по адресу: 141980 г. Дубна Московской области, ул. Жолио-Кюри, д. 6.

С диссертацией можно ознакомиться в научной библиотеке Объединенного института ядерных исследований.

Автореферат разослан (J Wjf^ 2013 г.

Ученый секретарь диссертационного совета, доктор физико-математических наук

Ю.А. Батусов

Общая характеристика работы

Представлены результаты теоретического исследования динамики пучка в синхротронах с системами подавления когерентных колебаний заряженных частиц на основе оригинальных моделей с целью уточнения имеющихся закономерностей и выявления ограничений, обусловленных использованием цифровых методов обработки сигналов в цепи обратной связи системы подавления. Выявлены количественные и качественные закономерности для сдвигов частот и декрементов затухания когерентных колебаний сгустков в синхротроне при наличии системы подавления, согласующиеся с экспериментальными данными, полученными на Большом Адронном Коллайдере (LHC) на протонном пучке. Актуальность работы

Разработка, проектирование, модернизация и эксплуатация ускорительных комплексов протонов и ионов, таких как LHC (ЦЕРН), У-70 (ИФВЭ, Протвино), FAIR (GSI, Дармштадт), NICA (ОИЯИ), включающих в себя несколько синхротронов с высокими требованиями физического эксперимента к интенсивности и фазовой плотности пучка или коллайдер с предельно возможной светимостью, сопровождаются решением научно-технических задач ускорительной физики и техники, среди которых можно выделить пять проблем, относящихся к заявленной теме диссертации:

1. Достижение высокой интенсивности пучка заряженных частиц.

2. Сохранение фазовой плотности пучка при переводе из синхротрона-инжектора в принимающий синхротрон.

3. Обеспечение повторяемости расчетных параметров пучка в различных ускорительных циклах.

4. Контроль параметров синхротрона с требуемой точностью без заметного ухудшения качества пучка при длительных циклах ускорения и накопления.

5. Внедрение цифровых технологий в управляющей системе, обеспечивающей демпфирование когерентных колебаний сгустков.

Повышение интенсивности пучков заряженных частиц является одной из основных задач в ускорительной физике и технике. Одним из факторов ограничения интенсивности пучка в синхротронах и накопителях являются когерентные неустойчивости. Число обнаруженных к настоящему времени различных типов когерентных неустойчивостей весьма велико. Обусловленные ими ограничения на интенсивность ускоряемого пучка заставляют разрабатывать методы их изучения и подавления. С этой целью в современных синхротронах широко используются системы подавления (СП) когерентных колебаний пучка заряженных частиц [Ы]. В таких системах демпфирующий дефлектор осуществляет коррекцию импульса пролетающих заряженных частиц пропорционально их отклонению от положения равновесия, измеренному в предшествующие моменты времени с помощью датчика положения центра тяжести пучка. Подобные СП позволили существенным образом повысить интенсивность ускоряемого пучка. В связи с этим в составе современных протонных синхротронов функционируют системы подавления когерентных колебаний, обеспечивающие, как это отмечено в литературном обзоре диссертации, декремент затухания колебаний а <, 0.012. При этом актуален вопрос о требованиях к СП'в зависимости от интенсивности ускоряемых сгустков и структуры ускоряемого пучка.

В современных ускорительных комплексах, в составе которых функционируют несколько синхротронов, при переводе сгустков из синхротрона-инжектора в принимающий синхротрон неизбежно возникают отклонения инжектируемых сгустков от расчетных положений их устойчивого равновесия. При этом инжектируемые заряженные частицы, циркулируя в синхротроне, совершают малые колебания вблизи своего положения устойчивого равновесия: бетатронные колебания в поперечном направлении (вертикальной и горизонтальной плоскостях) относительно мгновенной замкнутой орбиты (деформированной или расчетной) в ведущем магнитном поле и синхротронные колебания в продольном направлении относительно синхронной фазы. Из-за ошибок инжекции эти малые колебания частиц, вследствие их электромагнитного взаимодействия с другими

частицами, могут инициировать совместные колебания выделенной группы частиц (сгустка) и становиться когерентными. Вследствие энергетического разброса частиц в инжектируемом пучке малые когерентные колебания переходят в некогерентные (процесс естественного затухания когерентных колебаний), что сопровождается увеличением линейных размеров сгустков. Происходит уменьшение фазовой плотности пучка. Так, на ускорительном комплексе ЦЕРН, в составе которого цепочка из четырех синхротронов, нормализованный поперечный эмиттанс пучка возрастает почти в 1.5 раза в основном из-за ошибок инжекции. В связи с этим необходимо подавлять когерентные колебания сгустков после инжекции с высоким темпом за промежутки времени, которые должны быть короче времени перехода когерентных колебаний в некогерентные. Поэтому актуальны оценки предельных темпов подавления в контексте допусков на рост эмиттанса пучка в условиях двух конкурирующих процессов (естественного затухания и подавления когерентных колебаний).

Система подавления когерентных колебаний, будучи установленной на синхротроне и настроенной в соответствии с выбранной рабочей точкой по числу колебаний за оборот £?о, должна обеспечивать стабильное подавление когерентных колебаний при вариации £)о■ В литературе отсутствуют сведения о допусках на величину изменения которые необходимы при проектировании СП.

При исследовании динамики пучка в синхротроне при наличии СП используются дифференциальные уравнения малых колебаний [Ь2]:

т0у

где й — циклическая частота колебаний заряженной частицы с массой то в фокусирующем электромагнитном поле, у — релятивистский фактор, г — величина задержки сигнала в цепи обратной связи. Рекомендации [Ь2], получающиеся в сглаженном приближении при эмпирической замене импульсной силы Рл(х) на диссипативную, пропорциональную скорости х, как минимум необходимо уточнить при использовании цифровой СП, поскольку результирующая сила воздей-

ствия дефлектора на пучок пропорциональна суперпозиции сигналов с датчика положения, зафиксированных на нескольких оборотах.

Основным направлением проектирования и развития СП является применение методов цифровой обработки сигналов в цепи обратной связи. Существующие предсказания по тактовой частоте цифровых систем подавления и оценки влияния шумов квантования на рост эмиттанса пучка не приводят к каким-либо ограничениям, обусловленным использованием цифровых технологий. Вместе с тем отсутствуют предсказания по возможным ограничениям на величины декрементов когерентных колебаний при использовании цифровых СП.

Основные параметры ускорительного цикла контролируются и регулируются с помощью цепей обратной связи в процессе ускорения. При измерении параметров пучка в ряде случаев используются методы импульсного гармонического воздействия на сгусток, например, с помощью электромагнитного толкателя. Эти воздействия могут быть систематическими (например, с целью контроля качества пучка после инжекции) или случайными (например, при измерении частот бета-тронных колебаний в произвольный момент времени в пределах ускорительного цикла). Амплитуда возбужденных когерентных колебаний должна быть достаточной по величине для детектирования сигналов и измерения параметров колебаний с требуемой точностью. Однако система подавления, обеспечивая затухание когерентных колебаний, приводит к неизбежному сокращению промежутка времени, в течение которого измеряемый сигнал может быть идентифицирован на уровне шумов. Поэтому необходимы оценки для величин и продолжительности действия сил возбуждения когерентных колебаний при наличии СП с учетом их возможного влияния на рост эмиттанса.

Актуальность перечисленных выше физических задач и поиск путей их решения явились мотивацией к выполнению данной работы.

Цель диссертационной работы: - построение физико-математической модели, описывающей динамику пучка заряженных частиц после инжекции и при ускорении под влиянием управляю-

щих воздействий цифровой системы обратной связи, приводящих к затуханию когерентных колебаний;

- получение аналитических решений, позволяющих сопоставить расчетные данные (например, декременты затухания когерентных колебаний) с экспериментально наблюдаемыми характеристиками колебательного процесса в максимально приближенных к эксперименту условиях, в том числе при темпах подавления, близких к предельно возможным;

- сопоставление результатов разработанной модели с получающимися в линейном приближении решениями уравнений динамики пучка при использовании метода Крылова—Боголюбова—Митропольского и методов решения дифференциальных уравнений с отклоняющимся аргументом;

- формулирование перечня спецификаций для физико-технических обоснований проектов цифровых систем подавления в новых синхротронах с учетом опыта создания такой системы для ЬНС (ЦЕРН) в совместном проекте ЦЕРН и ОИЯИ.

Научная новизна

Приведены результаты теоретического исследования динамики сгустков заряженных частиц в синхротронах с цифровыми системами обратной связи с помощью методов дискретного анализа:

• уточнены зависимости и выявлены ограничения для декрементов затухания и частот когерентных бетатронных колебаний от величины баланса фаз, численно равной разности между набегом фазы бетатронных колебаний частицы при пролете от измерительного датчика до демпфирующего дефлектора и изменением фазы управляющего сигнала в цепи обратной связи;

• впервые продемонстрировано, что при использовании цифровых фильтров в цепи обратной связи систем подавления с задержкой сигнала более чем на один оборот возникают дополнительные к бетатронным моды колебаний сгустков, которые влияют на область устойчивости движения частиц;

• введено новое понятие: сепаратриса для декрементов затухания когерентных поперечных колебаний; приведен метод расчета, позволяющий установить кон-

фигурацию сепаратрисы и карты изолиний для декрементов затухания; показано, что высота сепаратрисы зависит от параметров цифрового фильтра и величины баланса фаз.

Практическая значимость

Разработанные в диссертации методы позволяют упростить качественный анализ цифровых систем подавления когерентных колебаний, а также усовершенствовать методы и оценки их характеристик с учетом найденных аналитических решений без трудоемких численных расчетов.

Достоверность и эффективность предложенных расчетных методов позволяют использовать их при проектировании современных синхротронов для выявления новых цифровых методов обработки детектируемых сигналов, не нарушающих условий устойчивости движения частиц, но повышающих помехоустойчивость управляющих процессов и позволяющих работать с меньшими уровнями входных сигналов и демпфирующих воздействий.

Разработанные методы применены автором диссертации для анализа режимов работы системы подавления когерентных поперечных колебаний пучка в ЬНС (ЦЕРН) и сопоставления проектных и фактических характеристик системы подавления, созданной в совместном проекте ЦЕРН и ОИЯИ.

Основные положения и результаты, выносимые на защиту

1. Результаты исследования динамики пучка в синхротронах с системами обратной связи, основанного на использовании методов дискретной математики.

2. Результаты расчетных и экспериментальных данных для сдвигов частот и декрементов затухания когерентных поперечных колебаний протонов в ЬНС.

3. Оригинальный подход для учета флуктуаций положения рабочей точки по числу бетатронных колебаний в синхротроне с цифровой системой обратной связи, приводящий к новому понятию - сепаратриса для декрементов затухания.

4. Доказано, что сдвиги частот и декременты затухания когерентных колебаний, вычисленные в линейном приближении по величине коэффициента обратной связи для сглаженных уравнений динамики пучка с помощью метода Крылова—

Боголюбова—Митропольского и методов решения дифференциальных уравнений с отклоняющимся аргументом, совпадают с соответствующими выражениями, вычисленными при использовании методов дискретного анализа.

5. Аналитический подход к решению задачи о подавлении остаточных колебаний сгустков после инжекции в синхротрон с цифровой системой обратной связи, имеющей нелинейную передаточную характеристику.

6. Основные результаты разработки, наладки и эксплуатации системы подавления когерентных поперечных колебаний пучка в ЬНС, обладающей существенными элементами новизны в принятых технических решениях.

Апробация работы

Результаты, приведенные в диссертации, докладывались на научных семинарах ОИЯИ, обсуждались в ускорительном подразделении ЦЕРН в связи с работами по созданию систем подавления когерентных колебаний пучка в ЬНС, докладывались на Международных конференциях по физике и технике ускорителей заряженных частиц, проходивших

- в России: 11иРАС'94 [16, 17], КиРАС'98 [21], 11иРАС'08 [23], ЯиРАС'Ю [25, 26], ЯиРАС'12 [15],

- в зарубежных странах: НЕАС'93 [18], РАС'93 [10], ЕРАС'94 [19], ЕРАС'96 [20], СООЬ'97 [8, 9], РАС'97 [11], ЕРАС'98 [12], РАС'99 [13], ЕРАС'ОО [22], 8агап1зеу'05 [1], ЕРАС'08 [14], ШАТ'09 [24], Загаг^су'ОЭ [3], Багатзеу'П [4].

Достоверность и обоснованность

Основные результаты физико-математического исследования цифровых систем подавления когерентных колебаний пучка заряженных частиц в синхротронах подтверждены экспериментальными работами, выполненными на Большом Адронном Коллайдере (ЬНС, ЦЕРН) при изучении поперечной динамики пучка при рекордных энергиях протонов от 450 ГэВ до 4 ТэВ.

Публикации

Материалы диссертации опубликованы в 36 печатных работах, из них 9 статей в рецензируемых журналах [1-9] и 17 статей в трудах конференций [10-26].

Личный вклад автора

Содержание диссертации и основные положения, выносимые на защиту, отражают персональный вклад автора в опубликованные работы. Полученные лично автором результаты опубликованы в виде 8 статей [1-8] в рецензируемых журналах, 10 докладов [11-13, 15, 19-22, 24, 26] в трудах конференций и 5 сообщений ОИЯИ [27, 29, 33, 34] и ЦЕРН [35].

В перечне устных докладов на международных конференциях включены 2 доклада [23, 25], которые диссертант представлял от имени коллаборации "LHC Damper" как первые экспериментальные результаты по наладке системы подавления когерентных поперечных колебаний протонного пучка в LHC. Подготовка других публикаций [9, 10, 14, 16-18, 28, 30-32, 36] проводилась совместно с соавторами на паритетных началах.

Все представленные в диссертации результаты получены лично автором.

Структура и объем диссертации

Диссертация состоит из введения, обзора литературы, 7 глав, заключения, списка сокращений и условных обозначений, библиографии и приложения. Общий объем диссертации 224 страницы, включая 33 рисунка и 3 таблицы. Библиография включает 138 наименований на 18 страницах.

Содержание работы

Во Введении обоснована актуальность диссертационной работы, сформулирована цель и аргументирована научная новизна исследований, показана практическая значимость полученных результатов, представлены выносимые на защиту основные положения.

В Литературном обзоре приведены научные результаты и экспериментальные данные, которые сопоставляются в диссертации с полученными автором результатами.

В Главе 1 изложено обоснование физико-математической модели цифровой широкополосной системы подавления когерентных поперечных колебаний

«жёсткого» сгустка. Динамика такого одиночного сгустка описывается с помощью линейных матричных уравнений движения заряженной частицы в синхротроне [ЬЗ] и линейных разностных уравнений для цифровых фильтров с конечной и бесконечной импульсной характеристикой [Ь4]. Состояние частицы задается с помощью матрицы-столбца X (л, л) на обороте п в точке 5 орбиты с периметром С0. Её первый элемент есть х(п, х) - смещение частицы относительно мгновенной замкнутой (равновесной или деформированной) орбиты, а второй элемент х'(п, х) - угол наклона траектории частицы. Для вычисления функци-

Рис. 1. Структурная схема цифровой системы подавления когерентных колебаний пучка [1.2]

ональных зависимостей (x(n,s)) и {x'(n,s)) для ансамбля частиц «жёсткого» сгустка используется процедура усреднения с помощью стационарных в широком смысле функций распределения в соответствии с общепринятыми моделями и предположениями [L5]. Для смещения центра тяжести сгустка х [л, .у] и угла наклона его траектории используется модельное приближение о равенстве дискретных и непрерывных усредненных величин смещений X[я,s] = (x(n,s)> и их производных. Такое же соглашение принимается для импульсного воздействия демпфирующего дефлектора (Ax'd(n)) = Дх^[и] и импульсного толкателя (Ах'е(п)) = Д*£[и] (см. рис. 1). Тогда [5, 27]:

Х[п + 1,í] = Х[п,s + С0] = M(í)£[и, j] + Ax'¿[л] МЛЕ + Дх'е[п] МеЕ, (1) где в матрице-столбце Ê имеем Ei = 0 и Е2 = 1. Для матриц перехода M (.v2|.fi) из точки si в позицию Í2 используются общепринятые обозначения [L3]: M (s) s M (s + C„|í), Md = M(s+ CoM, Me = M (s + C„|ïe). Дискретное смещение сгустка в датчике положения дг[л,5р] = отождест-

вляется с дискретным отсчетом на входе радиотехнической цепи обратной связи, пропорциональным входному напряжению ^¡„[л]. Поэтому для цепи обратной связи можно записать [5, 27]:

Ут[л] = Сф] + 5in Kin ф],

N, N,

Ы»] = Ь[т) Уф - т] и[п - т] + ^ с[т] УШг[п - т] и[п - т],

"1=0 т = 1

Ддг>] = X5out Fout[n] = лг50ц,KoUt Vm[n -q]u[n-q], (2)

где x = ± 1 — коэффициент, учитывающий полярность подаваемого на дефлектор DK управляющего напряжения, Sm - чувствительность измерительного датчика ВРМ, К,п — коэффициент усиления по напряжению входных усилителей, 5хр - погрешность измерения бетатронных смещений колеблющейся частицы в месте расположения ВРМ, Sout - передаточная характеристика дефлектора, K<,ut - коэффициент усиления по напряжению выходных усилителей, и [и] - функция Хевисайда, N{ - порядок фильтра с конечной импульсной характеристикой (КИХ-фильтр), Ь[т] - коэффициенты КИХ-фильтра в управляющем процессоре DSP, Nr - порядок фильтра с бесконечной импульсной характеристикой (БИХ-фильтр), с[т] - коэффициенты БИХ-фильтра [L4], Kfltr - выходное напряжение после фильтра, Fout - напряжение на выходе цепи обратной связи. В (2) учтено, что синхронизация отсчётов Vm[n] и Kout[nj может осуществляться с задержкой на q оборотов. Иначе говоря, если г(0,.ур) есть момент пролета частицы сквозь ВРМ после инжекции, г = rpd — время пролета частицы от ВРМ до DK, то корректирующий сигнал должен возникнуть в DK в момент времени t(q, jd) такой, что [5]

r + ?7b = f(i,Jd)-r(0,Jp). (3)

Величина q далее называется кратностью задержки. Обычно q = 0 на ускорителях с большой длиной дорожки, однако в электронных синхротронах для выполнения условий синхронизации с учётом реальной задержки в цепи обратной связи приходится выбирать отличные от нуля значения q.

Для решения полученной системы линейных разностных уравнений (1) и (2) можно использовать преобразование Лорана (^-преобразование) [Ь4]:

оо

у(г.*) =2{у[п,5]} = ММ*-". = 0 Уп < 0, а е [О.С0]. (4)

п=0

= (г,*) гп~1 ¿2 = £Ке8[у(2г,*) г?"1; гг],

где Г - замкнутый контур, окружающий все особые точки 2, функции у (г, л) на комплексной 2-плоскости, причем \гГ\ < 1. После выполнения Л-преобразова-ния в уравнениях (1) и (2) получаются следующие выражения {5, 27]:

Х(г,.т) =-;-г— I Н--—- МйЕ

+ Z{Ax'e[n]} Mjy (5)

Здесь I - единичная матрица, X[0,î] - начальное состояние сгустка, а матрица M имеет вид

~ - ê(z) ---

Mhs M(2o,z,j) = M (s) + - : MdTM(sp\s), (6)

(pdpp)l/2

где fii — функции Твисса [L5] в точках расположения ВРМ и DK, Т - квадратная 2x2 матрица, в которой Т21 = 1, а все остальные элементы равны нулю. Следует подчеркнуть, что оптические функции Твисса и набеги бетатронных фаз в матрицах M (i), M ¿(s + C0|.vd) и M(jp|î) заданы на орбите, где число бетатронных колебаний за оборот равно go- Это отражено в определении матрицы M в (6), где зависимость от Qtl указана в явном виде как M((2n, z, л). Функция g(z), входящая в (5) и (6), определена следующим образом [5]:

g(z) = z'" xK(z) Sin Soui p • (7)

В (7) введено обозначение K(z) = Kin KoUt H(z) для передаточной функции радиотехнической части цепи обратной связи, зависящей от передаточной функции [L4] цифрового фильтра H(z) = Z{VmAn\Y-

H(z) = Z{VatI[n]} = ^[0] + £ b[m] - £ c[m] . (8)

Собственные значения zr есть решения характеристического уравнения [5, 27]:

det (zrT — М(2о. = zr2-zrtrM(öo,zr,i) + detM(öo,zr,i)

= z?-(2 cos ¡x + g(zr) sm(ß - r¡))zr + 1 - g(zr) sin r] = 0, (9)

где ¡i = 2л Q0 есть набег фазы бетатронных колебаний за оборот, а т) - от ВРМ до DK. Уравнение (9) при g = 0 есть квадратное уравнение, имеющее два комплексно-сопряжённых корня zli2 = exp(±j 2kQq), т.е. хорошо известное решение для бетатронных колебаний [L5]. В этом случае величина arg z 1,2 = 2тг{0о} есть приращение фазы бетатронных колебаний за оборот ({go} - дробная часть Q0: —0.5 < {Q0} < 0.5), которое фиксируется как наблюдаемая величина в сигнале с датчика положения [L2], Этому приращению можно сопоставить циклическую частоту Qß и абсолютную величину Qß для {Qo}'-

Qß = |{ßo}|, Slß = Qßü)a, &>o = 2п/ To, QßT0 = 2jtQß. Величины Qß и Qß далее будут использоваться в качестве характеристики минимальной моды бетатронных колебаний. С учетом принятого в радиотехнике правила соответствия К(<у) = К(z = exp(j соТ0)) для g(z) можно записать [5]:

g(z) = g z~í ЩО-, g = Sin SoutКо /¡Цр, Ко s где g = - коэффициент обратной связи, причем [Ax¿(í2jg)| (ßdßp)l/2 =

g При g ф 0 физический смысл имеют решения, для которых все моды

колебаний являются затухающими (|zr| < 1). Следовательно, для асимптотических решений [L4] можно записать для X(z, s) из (5):

lim = Iim(z — 1)X(z,j) = 0. (10)

П-УОО 2-И

Таким образом, все полюсы zr функции (z — l)X(z,.v) должны лежать внутри круга единичного радиуса, то есть влияние погрешностей 8хр в X(z, s) должно быть исключено. Это возможно, если [5, 24, 27]

N,

g(z = 1) = 0 => ¿[0] + £ b[m] = 0.

m = l

Последнее достигается благодаря применению в цепи обратной связи узкополосного заграждающего фильтра (notch filter, Hn(z) = (1 — z-1)) [24, 27].

Приближенные решения характеристического уравнения при g «1 1 есть

Zl, = ехр (=42* (е„ - = е-«*«». (11)

где введены следующие обозначения [5, 15]:

о) = (±Qp +к)ао, Ф = г] + 2irq Q0 - argK(ü>),

а = {Q} = {Qo}-^*cos*, (12)

2К« 4л-К<>

ВД=КЛ«1Я(г=е^), Ко=|К(^)|, *5шФ(%)>0.

Определяющим параметром в (11) является баланс фаз Ф, величина которого равна разности набега фазы бетатронных колебаний (r¡ + 2nq Q0) для частицы при пролете от датчика положения до дефлектора (с учетом кратности задержки q) и изменения фазы детектируемого сигнала arg К(&>) в цепи обратной связи. Декремент затухания максимален в резонансном случае, когда [5, 15]

| sin Ф| = 1 {Q} = {Qo}. (13)

Таким образом, оптимальные условия для затухания колебаний могут быть получены при q = 0 не только в случае | sin r¡\ = 1, что является обычным требованием при проектировании СП, но и при отклонении от этого условия для набега фазы бетатронных колебаний от датчика положения до дефлектора, если такое отклонение может быть скомпенсировано цифровым фильтром.

Поскольку g(z) ос z~4~n< при Nr < Nf, то характеристическое уравнение (9) является алгебраическим уравнением степени (2 + q + Nf) относительно z. Коэффициенты уравнения (9) — действительные числа. Следовательно, такое алгебраическое уравнение имеет парные комплексно-сопряжённые корни и/или различные вещественные корни (в частных случаях некоторые из вещественных корней могут совпадать). Решая алгебраическое уравнение (9), можно получить уточненные зависимости для декрементов затухания при росте величин коэффициента обратной связи g с учетом передаточной функции //(z), зависящей от структурной схемы цифрового фильтра согласно (8). Таким образом, при использовании цифровых фильтров в цепи обратной связи систем подавления с

15

задержкой сигнала более чем на один оборот возникают дополнительные моды колебаний сгустков, которые ограничивают область устойчивости движения частиц при высоком темпе подавления колебаний [3, 5, 27].

Рис. 2. Расчетные декременты затухания и результаты измерений горизонтальных (а) и вертикальных (б) когерентных поперечных колебаний пучка протонов с энергией 3.5 ТэВ в ЬНС [5, 15]

Так, в системе подавления когерентных поперечных колебаний пучка в ЬНС используются последовательно соединенные по1:сЬ-фильтр и фильтр Гильберта [5]. На рис. 2 продемонстрировано совпадение экспериментальных данных с расчетными зависимостями а (#) для декрементов затухания колебаний при рекордной энергии 3.5 ТэВ протонного пучка в ЬИС в широком диапазоне величин коэффициента обратной связи < 0.08) вплоть до значений, для которых декремент затухания колебаний (а < 0.05) по крайней мере в четыре раза превосходит достигнутый на других протонных синхротронах уровень (а ~ 0.012) [5, 15].

Результаты первой главы опубликованы в работах [3, 5, 8, 10, 15, 20, 24, 27].

В Главе 2 приведены результаты исследования проблемы о конфигурации области устойчивого затухания когерентных поперечных колебаний сгустков при вариации основных параметров синхротрона, влияющих на бетатронное движение частицы [4, 8, 19]. Динамика взаимодействующих сгустков в синхротроне для отклонений хт[п] на каждом «-м обороте (т е [1, АГь], Къ - число сгустков) описывается с помощью линейных дифференциальных уравнений, в которых электромагнитное взаимодействие сгустков задано с помощью хвостовых полей, и линейных разностных уравнений для напряжений Ут[п], возникающих

в цифровых фильтрах. Для решения полученной системы линейных уравнений использован метод Z-преобразования. Возникающая после ^-преобразования система линеиных уравнении для \т (z) = Z{xm[n]} с постоянными коэффициентами с помощью общеизвестной процедуры [L6] приводится к нормальным координатам v,„ (z). В итоге получается характеристическое уравнение для собственных значений zr, по форме совпадающее с характеристическим уравнением для «жёсткого» сгустка (9), в котором, однако, фазовые набеги для поперечных колебаний должны быть вычислены с учетом эффективного числа бетатронных колебаний за оборот Qm(zr) = Qt + AQm(zr) для конкретной моды колебательного процесса, описываемого с помощью нормальных координат vm(zr) (здесь Qc = <20 -(- Ago, где А<2о - квазистатическая поправка к go)- Каждому собственному значению zr можно сопоставить циклическую частоту <Л(к), так что )о>(к)То = j 2пк + In zr(Qm), где к - целое число. Поскольку в первом приближении по силе взаимодействия сгустков при выключенной системе подавления имеем zr = exp(±j 2kQc), то для комплексных поправок AQт(0(к)) в этом же приближении получаются выражения, совпадающие с известным результатом для неустойчивостей сгустков [32] и сплошного пучка [35] в зависимости от частоты a>(fc) = (±{Qc}+k)a>a. Наряду с периодичностью частот собственных волн, равной кы{для каждого к происходит расщепление ш^) на Къ сателлитов вследствие того, что AQm есть функция от К^1 lnz,. Для большого числа сгустков Кь максимальный разброс AQimaxl для сателлитных мод возникает в синхротронах с большой протяженностью орбиты в случае проникновения электромагнитного поля сгустков сквозь стенки вакуумной камеры с апертурой 2b [32]:

1+j2яу0 b2

где Т\ - время диффузии поля в проводящие стенки камеры; уь - параметр, характеризующий взаимодействие сгустков; Nb - число частиц в сгустке; Гц -классический радиус частицы [L5]; - среднее значение -функции по длине орбиты. Разброс частот сателлитных мод пропорционален инкременту неустойчивости Vb&b. Набег фазы бетатронных колебаний г) от датчика до дефлектора

зависит от (<2с + Яс А0<шах)(&>№)))- Поэтому обеспечить строгое выполнение оптимальных условий (13) для подавления всех мод когерентных поперечных колебаний взаимодействующих сгустков одновременно невозможно [32].

Таким образом, наряду с зависимостью декремента затухания а от коэффициента обратной связи g возникает зависимость декремента от эффективной величины Q или от мгновенной величины <2 при случайном смещении ¿>20 положения рабочей точки относительно расчетной величины (?0- Трехмерную

Рис. 3. Карты: (а) изолиний a{g,SQo) = щ [4] и (б) изоамплитуц <2е - Qp) - щ [5]

зависимость как функции двух переменных можно представить на

плоскости (я, Я<2о) в виде карты изолиний [4] — семейства линий, для каждой из которых величина декремента имеет постоянное значение а/, например:

а&Яво ■)=(*/. «0 = Т0/таес, а, = ¿/80 V / е [1.8],

где гйес - постоянная времени фазового перемешивания частиц. На рис. За приведен пример такой карты изолиний для СП с фильтром Гильберта [4]. Поскольку декремент подавления когерентных колебаний Г0/ га должен превосходить величину а0 = Т„/ тЛсс, то замкнутую кривую = а0 можно отождествлять с сепаратрисой для декрементов затухания когерентных поперечных колебаний. Высота сепаратрисы зависит от параметров цифрового фильтра и баланса фаз [4]. Это означает, что требования к стабильности положения рабочей точки в синхротроне должны формулироваться с учетом высоты сепаратрисы.

В [5] выявлены закономерности для резонансной зависимости амплитуд а. вынужденных колебаний сгустков в ответ на ограниченное во времени гармоническое воздействие с циклической частотой сое = (<2е — 0,р)аМ) от параметров цепи обратной связи с цифровым фильтром. Амплитуда вынужденного колебания возрастает при приближении £)е к <2р, но не обращается в бесконечность, как это происходит при резонансе в отсутствие СП. Установлено, что коэффициент резонансного усиления амплитуды вынужденных колебаний есть /?етах) = [ /2а. Трехмерную зависимость ае(п = г/Г(), ()е - как функцию двух переменных можно представить на плоскости («, ()е — (?//) в виде карты изоамплитуд ае(п, <2е — ()р) = а/ (см. рис. 36). Виден эффект смещения максимума резонансной кривой относительно <2р. Установлено, что этот эффект зависит от баланса

(О) 0.281 0.280 0.279 0.278 0.277 і 0.276 0.275

1/2« 10!

-160 -140 -120 -100 -80 -60 -40

-160 -140 -120 -100 -80 -60 -40

Рис. 4. Расчетные зависимости величины {()} (а) и коэффициента резонансного усиления (б) от фазового параметра <р фильтра Гильберта (сплошные линии) в сопоставлении с результатами экспериментальных измерений [7]

фаз [5]. Данный эффект был применен для тонкой настройки баланса фаз системы подавления когерентных поперечных колебаний протонного пучка с энергией 450 ГэВ в ЬНС. На рис. 4 приведены экспериментальные данные, подтверждающие наличие линейной поправки к когерентной частоте бетатронных колебаний при малом отклонении баланса фаз от оптимального значения [7].

Результаты второй главы опубликованы в работах [3-5, 7, 8, 16, 19, 28, 29, 31, 32, 34, 35].

В Главе 3 приведены результаты исследования многоканальных систем подавления когерентных поперечных колебаний «жёсткого» сгустка.

Основанием для изучения многоканальных систем подавления стали экспериментальные данные, полученные на синхротроне SPS (ЦЕРН) при исследовании устойчивости высокоинтенсивного пучка заряженных частиц при приближении числа бетатронных колебаний к полуцелому числу. Было обнаруже-

103

Il О-►

Li о1

tune-0.582 tune=0.536 tune=0.522

unstabil

Сч

P

too low gain

good damping 1

10

10"

10" 10"' feedback gain (a.u.)

Рис. 5. (а) Результаты измерений постоянной затухания rd и (б) расчетные величины декрементов затухания а = Т0/rd в зависимости от g для значений Q0\ 26.582 (сплошная линия), 26.536 (пунктирная), 26.522 (точечная), 26.500 (штрих-пунктирная) [36]

но, что для величин коэффициента обратной связи, превышающих инкремент неустойчивости, пучок устойчив при Qq = 2fr582, причем постоянная времени затухания уменьшается с ростом g и стабильность пучка в целом улучшается (см. рис. 5а) [36]. Однако при Qn — 26.536 диапазон величин g, при которых захваченный пучок был устойчив, резко сокращался. Ещё уже данный диапазон был при Qо = 26.522, причем, начиная с некоторых значений g, не удавалось сохранить пучок даже после инжекции. На рис. 5б представлены расчетные зависимости [36], полученные с учетом результатов предыдущих глав диссертации для области устойчивости пучка по диапазону допустимых величин коэффициента обратной связи при приближении числа бетатронных колебаний к полуцелому числу в синхротронах с одноканальной системой подавления. Видно существенное её сокращение, что соответствует экспериментальным измерениям на SPS.

С учетом этих результатов для классической (одноканальной) системы подавления приведено физико-математическое обоснование многоканальной и бы-

строй (двухканальной) [9, 18] систем подавления когерентных поперечных колебаний пучка заряженных частиц. В результате выполненного исследования [30, 31, 33] динамики пучка в синхротроне с многоканальной системой подавления когерентных поперечных колебаний, основанного на использовании метода дискретных преобразований для решения уравнений бетатронных колебаний, установлено, что при различных расстановках датчиков и дефлекторов, а также связях между ними, включая перекрестные, собственные значения характеристического уравнения зависят от набегов фаз бетатронных колебаний между датчиками и дефлекторами, но не зависят от параметров матрицы Твисса ß и а в местах их расположения. Получены аналитические выражения для декрементов затухания колебаний для быстрой СП при оптимальной расстановке датчиков и дефлекторов.

Результаты третьей главы опубликованы в работах [9, 17, 18, 30, 31, 33, 36].

В Главе 4 рассмотрены основные физические понятия, математический аппарат и методы решения системы линейных разностных уравнений [6, 26], описывающих продольную динамику заряженных сгустков в синхротроне при наличии системы подавления их когерентных продольных (синхротронных) колебаний. Продемонстрировано, что на плато магнитного поля и при ускорении демпфирование когерентных синхротронных колебаний возможно только при фазовой модуляции ускоряющего ВЧ-напряжения. Показано совпадение экспериментальных данных, полученных на синхротроне PETRA II (бустер для HERA в DESY), и расчетных величин для декрементов затухания когерентных продольных колебаний в апериодическом режиме.

Результаты четвертой главы опубликованы в работах [6, 26].

В Главе 5 показано, что дифференциальные уравнения с отклоняющимся аргументом [L7] описывают динамику пучка в синхротронах с системами обратной связи в рамках сглаженного приближения, традиционно используемого в теоретических исследованиях [L2] и заключающегося в замене локального демпфирующего воздействия непрерывно действующей корректирующей силой.

В [7] продемонстрировано, что для решения таких дифференциальных уравнений с отклоняющимся аргументом, описывающих динамику пучка в сглаженном приближении в синхротроне с системой обратной связи, можно использовать метод интегральных преобразований (одностороннее преобразование Фурье) и метод усреднения Крылова—Боголюбова—Митропольского (КБМ) [Ь8]. Например, в синхротроне с системой обратной связи с простейшим КИХ-фильтром сглаженное уравнение поперечной динамики пучка имеет вид [7, 27]

х(1) + П~х(1) = Тп

г^- £ Ь[т)(х{1 -х-(т+ д)Т0) + 8х0) и (Г - г - (т + $)Т0), (14)

^ ^ т—0

где Но — нормировочная постоянная и £2 = Qo соо- Если использовать метод одностороннего преобразования Фурье для решения уравнения (14), то возникает следующее характеристическое уравнение [7]:

—оР~ + Я2 - £^Н(а>) = 0.

То Н0

Данное уравнение, будучи трансцендентным, имеет бесконечное множество решений. Однако при слабой зависимости передаточной характеристики фильтра от частоты, когда

О ¿Н(со)

«1,

|Я(£2)|

01= П

¿ш

можно воспользоваться методом Эйлера для решения характеристического уравнения, что в первом приближении по малому параметру £ для собственных значений а> г = £2 + <Уц А (? + j а/То приведет к следующим выражениям [7]:

а =--лгвтФ, АО =---——лгсовч', (15)

2 Но 4тг Н0

где величина баланса фаз есть

Ф(£2) = ^ + 2л-<72о-агвЯ(£2), т) = 0. г = 2лгбог/Г0. Здесь //(Я) — передаточная характеристика фильтра:

Н(П) = >я£2Го ЕЕ 2ят£о.

т= 0

В [7] показано, что согласно методу усреднения КБМ для приближенных уравнений а = —аа и ^ = ((20 + АбМь относящихся к основной гармонике, для а и Д<2 получаются выражения, совпадающие с (15).

В синхротроне сглаженное уравнение продольной динамики пучка при наличии системы обратной связи с простейшим КИХ-фильтром имеет вид [7]:

I (г) + =

Ь[т]($(1 - х - (т + «)7о) + и(г-т-(т + д)Т0)

gx d

TqHo dt

(16)

Lm=0

где f = &4>s - фазовое отклонение частицы от синхронной фазы, £2 = Í2S = vscuo - циклическая частота синхротронных колебаний, vs - число таких колебаний за оборот, <5£о - статическая погрешность измерения фазового отклонения частицы от синхронной фазы. Выполняя одностороннее преобразование Фурье, получаем следующее характеристическое уравнение [7]:

+ Q2s +j = 0.

lono

Если обеспечена слабая зависимость передаточной характеристики фильтра от частоты, то, воспользовавшись методом Эйлера, получаем следующие выражения для собственных значений [7]:

«V = Í2s + j = ^s + «oAvs + j £,

T0H0 1 o

где

Ф(Й5) -t] + 2лqvs - arg#(£2s), r¡ = £2sr = 2jtvsx/T0 .

as = | ШШ x eos ф(0.), A„ = f HÍM ж Sin Ф(П§>.

2 #o 4jr "o

Эти же выражения для ors и Avs нетрудно получить с помощью укороченных уравнений для основной гармоники при использовании метода усреднения КБМ. Вновь следует подчеркнуть, что определяющим параметром асимптотических свойств когерентных синхротронных колебательных мод является величина баланса фаз Ф(£28). Подбирая параметры фильтра b[m], можно обеспечить оптимальные условия (sin Ф (S2s) = 0) дая подавления когерентных синхротронных

23

колебаний, когда декремент затухания колебаний максимален и колебательный процесс имеет резонансный характер (Дvs = 0).

Таким образом, в линейном приближении по величине коэффициента обратной связи для сдвигов частот и декрементов затухания когерентных поперечных и продольных колебаний получаются аналитические выражения, которые совпадают с соответствующими выражениями оригинальных моделей, где использовались методы дискретной математики для учета локальности корректирующего воздействия, пропорционального локальному когерентному смещению сгустка.

Результаты пятой главы опубликованы в работе [7].

В Главе 6 приведены некоторые результаты исследования динамики пучка в синхротроне с системой подавления когерентных поперечных колебаний с нелинейной передаточной характеристикой цепи обратной связи:

У^Д*>]=<Г/(*[«]), €=gx,

где f{x[n\) есть нелинейная функция отх[и]. Используя (1) для смещений «жёсткого» сгустка в датчике положения на трех последовательных оборотах при наличии корректирующего демпфирующего воздействия A.x'¿[n] в отсутствие внешнего воздействия толкателя (Дх£[и] = 0), в [12]'получено следующее уравнение:

х[п + 2] - 2 cosO) х[п -I- 1] + дс[я] =

= y/Md (sin(ц - т/) Ax'd[n + 1] + sm(n) Ax'Jn]) • 0 7)

В [13] для уравнения (17) с малым параметром е, по аналогии с методом усреднения КБМ для нелинейных дифференциальных уравнений, получены следующие укороченные уравнения для амплитуды ап и фазы фп основной гармоники х[п] = ansmi¡fn:

d 2П

о

2л-

= - J f{Qn sin fn) sin(^„ + n) .

o

Эти выражения по форме совпадают с укороченными уравнениями для нелинейного дифференциального уравнения, в которых фазовый сдвиг г] = фрд, связанный с задержкой демпфирующего сигнала, присутствует в явном виде. Отсюда, например, для кубической нелинейности /(*[«]) = х\п\ + (х[я])3£з имеем [12]:

Это выражение для зависимости амплитуды колебаний от времени совпадает с хорошо известным результатом для нелинейного дифференциального уравнения для осциллятора со слабым трением [Ь8]. Для величины отношения амплитуд линейных и нелинейных колебаний получаем [12, 21]:

Таким образом, «жесткая» нелинейность > 0) приводит к небольшому увеличению скорости убывания амплитуды когерентных поперечных колебаний, причем этот эффект квадратичным образом зависит от начальной амплитуды колебаний До. В случае «мягкой» нелинейности (#з < 0) возникает небольшое уменьшение скорости убывания амплитуды колебаний.

В [13, 21, 22] представлены основные закономерности затухания когерентных поперечных колебаний сгустка для нелинейных характеристик: кулоновское трение, ограниченное по величине корректирующее воздействие. Показано, что такие нелинейные характеристики приводят к появлению дополнительных мод колебаний и к ограничениям областей устойчивого движения частиц.

В [2, 6] установлены количественные закономерности увеличения эмиттан-са пучка Де, возникающего вследствие ошибок инжекции Аа0, в зависимости от постоянных времен затухания когерентных колебаний га, нарастания когерентной НеуСТОЙЧИВОСТИ Гщ5( И Перехода КОГереНТНЫХ КОЛебаНИЙ В НеКОГереНТНЫе Г(1сс в синхротронах с системами обратной связи с линейной и нелинейной передаточной характеристикой цепи обратной связи. Так, для линейной передаточной

Доехр(—па)

8* . а = — БН!//.

а„ =

характеристики получено [2, 6]:

Д£ = £-£0 = (Д«^ (18)

£() е0 2 \ тл г instJ

где Сто - среднеквадратичный размер инжектируемого пучка. Обозначим:

-L--L > к

**inst ^dcc

Если Дао — Сту, то относительный рост эмиттанса Ае/е при к = 1 составляет 12.5 %. Такую величину прироста эмиттанса можно оценивать как предельно допустимое значение. Если предполагается, что относительный рост эмиттанса Ае/е должен быть не более 5 %, то есть < 0.1, то к > 2.

Результаты шестой главы опубликованы в работах [2, 6, 12, 13, 21, 22]. В Главе 7 приведено описание системы подавления когерентных поперечных колебаний пучка в LHC ("LHC Damper"), которая предназначена [23, 25]:

- для демпфирования когерентных поперечных колебаний, возникающих из-за ошибок инжекции,

- для подавления когерентных поперечных неустойчивостей,

- для очистки межсгусткового пространства от фоновых частиц в сепаратрисах, которые не используются в режиме столкновения сгустков,

- для возбуждения когерентных поперечных колебаний при проведении измерений, связанных с исследованиями поперечной динамики пучка.

Созданная система подавления поперечных колебаний протонного пучка с рекордно высокой энергией 0.45 - 4 ТэВ, работающая в составе крупномасштабного сверхпроводящего синхротрона-накопителя, обладает существенными элементами новизны в принятых технических решениях [1, 14, 23]:

• большая величина интеграла отклоняющего поля в демпфирующих дефлекторах (около 900 кВ на оборот, на пластинах дефлектора номинальное напряжение ±7.5 кВ на частоте 1 МГц);

• цифровая обработка информационных сигналов вследствие длительных задержек сигналов обратной связи, кратных периоду обращения пучка 7» = 88.93 мкс по орбите с периметром Со = 26658.883 м;

• малая нижняя граничная частота системы поперечной обратной связи (большой периметр орбиты) и учет особенностей изменения импедансов устройств и пучка с частотой (полоса частот 3 кГц - 1 МГц для уровня —3 дБ от Ктах);

• отсутствие резонансов вплоть до проектной частоты (40 МГц) следования сгустков для усилителей и до ускоряющих частот (400 МГц), отвечающих структуре пучка, для дефлекторов;

• адаптация конструкции высоковольтных дефлекторов пучка к криогенно-ва-куумному окружению и к образованию электронных облаков.

В ОИЯИ и на российских предприятиях создано уникальное оборудование (мощные широкополосные усилители и высоковакуумные дефлекторы) для системы подавления "LHC Damper" в совместном проекте ЦЕРН - ОИЯИ, предназначенное для работы в новой области энергий. Проектные параметры этого оборудования достигнуты, а для амплитуды корректирующего воздействия — превышены: получена величина отклоняющего поля более 11 кВ в полосе частот до 100 кГц, что превышает проектный уровень 10.5 кВ [25].

На рис. 6 представлены измеренные зависимости величин когерентных бе-татронных смещений одного сгустка с числом протонов ~ 1 х 1011 от номера оборота после инжекции [25]. Получена величина постоянной времени затухания колебаний га = 40 Г(), что соответствует проектному значению. Минимальная величина rd составляет около 10 оборотов при энергии инжекции [L1],

Рис. 6. Подавление горизонтальных ошибок при инжекции (beam 1) при выключенной (а) и включенной (б) СП; сигналы с датчиков Q1 (верхний луч) и Q9 (нижний луч) [25]

Для настройки системы подавления когерентных бетатронных колебаний протонного пучка в LHC были использованы методики, основанием для которых стали результаты, полученные автором настоящей диссертации (см. рис. 2 и рис. 4). Качественные и количественные закономерности для допусков и ограничений для систем подавления когерентных поперечных колебаний взаимодействующих сгустков подтверждены экспериментальными данными, полученными на LHC при формировании протонного пучка высокой интенсивности при ин-жекции и максимальной энергии протонного пучка 4 ТэВ.

Система подавления когерентных поперечных колебаний протонного пучка и пучка ионов свинца в LHC успешно использовалась и безотказно функционировала в период с мая 2010 г. до окончания Run I в феврале 2013 г. для подавления ошибок инжекции с максимально возможным декрементом затухания колебаний, для стабилизации пучка при ускорении и при соударениях [L1],

Результаты седьмой главы опубликованы в работах [1, 5, 14, 15, 23, 25].

В Заключении сформулированы основные результаты диссертации.

1. Разработана оригинальная концепция исследования физики пучков в синхротронах с системами подавления когерентных колебаний заряженных частиц, что позволило выявить качественно новые ограничения и закономерности устойчивого затухания колебаний частиц, обусловленные цифровыми методами обработки сигналов в цепи обратной связи.

2. Предложены новые подходы для решения фундаментальных задач динамики сгустков в синхротронах с системами обратной связи, основанные на методах дискретного анализа, позволяющие проводить аналитические расчеты устойчивых режимов затухания когерентных колебаний сгустков в широком диапазоне величин коэффициента обратной связи:

- установлено, что определяющим параметром для собственных частот когерентных поперечных (бетатронных) и продольных (синхротронных) колебаний является величина баланса фаз, численно равная разности между набегом фазы

когерентных колебаний частицы при пролете от детектора до демпфирующего дефлектора и изменением фазы детектируемого сигнала в цепи обратной связи;

- доказано, что при использовании цифровых фильтров в цепи обратной связи систем подавления с задержкой сигнала более чем на один оборот возникают дополнительные к бетатронным моды колебаний сгустков, которые ограничивают область устойчивости движения частиц.

3. Продемонстрировано совпадение расчетных зависимостей и экспериментальных данных для декрементов затухания когерентных поперечных колебаний протонного пучка в LHC при энергии 3.5 ТэВ в широком диапазоне величин коэффициента обратной связи.

4. Доказана перспективность разработанной физико-математической концепции для решения задач устойчивости взаимодействующих сгустков в синхротронах с системами обратной связи:

- получен новый результат о расщеплении когерентной частоты взаимодействующих сгустков на сателлитные моды, что обусловлено структурой пучка, и сформулировано заключение о невозможности обеспечения строгого выполнения оптимальных условий для подавления когерентных колебаний сгустков вследствие разброса частот сателлитных мод, который пропорционален инкременту неустойчивости;

- продемонстрировано, что область устойчивости пучка по диапазону допустимых величин коэффициента обратной связи при приближении числа бетатрон-ных колебаний к полуцелому числу в синхротронах с классической (одноканаль-ной) системой подавления существенно сокращается, что соответствует экспериментальным измерениям на синхротроне SPS (ЦЕРН);

- представлено физико-математическое обоснование новой быстрой системы подавления когерентных поперечных колебаний с двумя цепями обратной связи, которая может обеспечить подавление бетатронных колебаний независимо от их фазы в местах расположения измерительных датчиков.

5. Введено новое понятие «сепаратриса для декрементов затухания когерентных

поперечных колебаний», ставшее основой подхода для учета влияния флуктуа-цнй положения рабочей точки по числу бетатронных колебаний на декременты затухания когерентных колебаний сгустка с учетом структурной схемы цепи обратной связи. Установлено, что требуемая стабильность положения рабочей точки в синхротроне должна обеспечиваться с учетом высоты сепаратрисы.

6. Выявлены закономерности для резонансной зависимости амплитуд вынужденных колебаний сгустков в синхротроне от параметров цифровой цепи обратной связи. Уточнена зависимость смещения максимума резонансной кривой от величины баланса фаз. Данный эффект был применен для тонкой настройки баланса фаз системы подавления когерентных поперечных колебаний протонного пучка в ЬНС. Приведены экспериментальные данные, подтверждающие наличие линейной поправки к когерентной частоте бетатронных колебаний при малом отклонении баланса фаз от оптимального значения.

7. Доказано, что сдвиги частот и декременты затухания когерентных колебаний, вычисленные в линейном приближении по величине коэффициента обратной связи для сглаженных уравнений динамики пучка с помощью метода Крылова-Боголюбова—Митропольского и методов решения дифференциальных уравнений с отклоняющимся аргументом, совпадают с соответствующими выражениями, вычисленными при использовании методов дискретного анализа.

8. Изложен аналитический подход к решению задачи о подавлении остаточных колебаний сгустков после инжекции в синхротрон с цифровой системой обратной связи, имеющей нелинейную передаточную характеристику. Для решения разностных нелинейных уравнений динамики пучка использован приближенный подход, аналогичный методу усреднения Крылова—Боголюбова—Митропольского. Выявлены основные закономерности затухания когерентных колебаний сгустка для ряда нелинейных передаточных характеристик.

9. Установлены закономерности увеличения эмиттанса пучка, возникающего вследствие ошибок инжекции, в зависимости от декремента затухания когерентных колебаний, инкремента когерентной неустойчивости и постоянной перехода

когерентных колебаний в некогерентные в синхротронах с системами обратной связи с линейной и нелинейной передаточной характеристикой. 10. Приведены основные результаты разработки, наладки и эксплуатации системы подавления когерентных поперечных колебаний пучка в LHC, обладающей существенными элементами новизны в принятых технических решениях.

Список публикаций

1. Жабицкий В.М. Системы подавления поперечных колебаний пучка в современных синхротронах // Письма в ЭЧАЯ. 2006. Т. 3, вып. 7(136). С. 96-102.

2. Жабицкий В. М. О росте эмиттанса пучка вследствие ошибок инжекции в синхротронах с нелинейными системами подавления когерентных поперечных колебаний // Письма в ЭЧАЯ. 2008. Т. 5, вып. 1(143). С. 86-93.

3. Жабицкий В.М. Устойчивость ионного пучка в синхротронах с цифровыми системами подавления когерентных поперечных колебаний II Письма в ЭЧАЯ. 2010. Т. 7, вып. 7(163). С. 742-747.

4. Жабицкий В. М. Устойчивость пучка в синхротронах с цифровой системой подавления когерентных поперечных колебаний в зависимости от частоты бетатронных колебаний // Письма в ЭЧАЯ. 2012. Т. 9, N0.4-5(174-175). С. 563-568.

5. Жабицкий В. М. Отклик пучка на переходный процесс в синхротронах с цифровой системой подавления когерентных поперечных колебаний // Письма в ЭЧАЯ. 2013. Т. 10, No.l(178). С. 120-129.

6. Жабицкий В. М. О росте продольного эмиттанса сгустков из-за ошибок инжекции в синхротронах с системой демпфирования когерентных продольных колебаний // Письма в ЭЧАЯ. 2013. Т. 10, No.2(179). С. 236-244.

7. Жабицкий В. М. Динамика пучка в синхротронах с системами обратной связи // Письма в ЭЧАЯ. 2013. Т. 10, No.6(183). 11 с. (в печати).

8. Zhabitsky V. M. Tranzverse Feedback System with a Digital Filter // Nucl. Inst. & Meth. in Phys. Res. A. 1997. Vol. 391. Pp. 96-99.

9. Gleibman E.M., Zhabitsky V.M., Ivanov I.N., Kononov G.I., Lebedev N.I., Likhachev A. G., Melnikov V. A., Malakhov N. A., Pilyar N. V., Rukojatkina T. V., Scheulin A. S., Merker S. E., Mikheev M. S. System for initial betatron oscillation damping for UNK I stage beam // Nucl. Inst. & Meth. in Phys. Res. A. 1997. Vol. 391. Pp. 100-102.

10. Zhabitsky V.M., Korenev I.L., Yudin L.A. Transverse Feedback System with Digital Filter // Proc. of the 1993 Particle Accelerator Conference, 17-20 May 1993, Washington, USA. IEEE, 1993. Pp. 2543-2545. ISBN: 0-7803-1203-1.

11. Zhabitsky V.M. Damping of Injection Oscillations // Proc. of the 1997 Particle Accelerator Conference, 12-16 May 1997, Vancouver, Canada. IEEE, 1998. Pp. 2329-2331. ISBN: 0-7803-4376-X.

12. Zhabitsky V. M. Theory of a Transverse Feedback System with a Nonlinear Transfer Function // Proc. of the Sixth European Particle Accelerator Conference, 22-26 June 1996, Stockholm, Sweden. Institute of Physics, 1998. Pp. 1380-1382. ISBN: 0-7503-0579-7.

13. Zhabitsky V. M. Nonlinear Damping of Injection Oscillations // Proc. of the 1999 Particle Accelerator Conference, 29 March - 2 April 1999, New York, USA. IEEE, 1999. Pp. 1222-1224. ISBN: 0-7803-5573-3.

14. Baudrenghien P., Hofle W., Killing F., Kojevnikov I., Kotzian G., Louwerse R., Montesinos E., Rossi V., Schokker M., Thepenier E., Valuch D., Gorbachev E. V., Lebedev N. I., Makarov A. A., Rabtsun S. V., Zhabitsky V. M. LHC Transverse Feedback System and its Hardware Commissioning // Proc. of the 11th European Particle Accelerator Conference, 23-27 June 2008, Genoa. The European Physical Society Accelerator Group (EPS-AG), 2008. Pp. 3266-3268. ISBN: 978-92-9083-315-4.

15. Zhabitsky V.M. Transient Response of a Digital Transverse Feedback System in Synchrotrons // Proc. of XXIII Russian Particle Accelerators Conference Ru-PAC-2012, 24-28 September 2012, St.Petersburg, Russia. Saint Petersburg State University, 2012. Pp. 73-75. ISBN: 3-95450-125-0.

16. Жабицкий B.M., Коренев И. Л., Юдин Л. И. Подавление резистивной неустойчивости сгруппированного пучка dиспользованием цепей обратной связи // XIII Всесоюзное Совещание по ускорителям заряженных частиц, 13-15 октября 1992, Дубна. Т. I. ОИЯИ D9-92-455, Дубна, 1993. С. 78-81.

17. Жабицкий В.М., Чернов П. А. Влияние ошибок в расстановке датчиков и толкателей на параметры системы подавления резистивной неустойчивости в I ступени УНК // XIII Всесоюзное Совещание по ускорителям заряженных частиц, 13-15 октября 1992, Дубна. Т. I. ОИЯИ D9-92-455, Дубна, 1993. С. 82-86.

18. Ivanov I. N., Zhabitsky V. М., Korenev I. L., Yudin L. A., Kopylov L. I. Requirements to the Feedback System to Damp Resistive Instability (I Stage UNK) // Proc. of the XVth International Conference on High Energy Accelerators, 20-24 July 1992, Hamburg, Germany. Vol. II. World Scientific, 1993. Pp. 1034-1036.

19. Zhabitsky V. M. Theory of Multi-Bunch Resistive Wall Instability Damping using Feedback System with a Digital Filter // Proc. of the Fourth European Particle Accelerator Conference, 27 June - 1 July 1994, London, England. World Scientific, 1994. Pp. 1090-1092.

20. Zhabitsky V. M. Transverse Feedback System with a Digital Filter and Additional

Delay // Proc. of the Fifth European Particle Accelerator Conference, 10-14 June 1996, Sitges, Spain. Institute of Physics, 1996. Pp. 1833-1835.

21. Жабицкий B.M. Нелинейный режим подавления остаточных колебаний пучка после инжекции в синхротрон // XVI Совещание по ускорителям заряженных частиц, 20-22 октября 1998, Протвино. Т. 2. ИФВЭ, 1999. С. 134-139.

22. Zhabitsky V. М. Nonlinear Damping of Injection Oscillations // Proc. of the Seventh European Particle Accelerator Conference, 26-30 June 2000, Vienna. European Physical Society Interdivisional Group on Accelerators (EPS-IGA), 2000. Pp. 1933-1935.

23. Gorbachev E.V., Lebedev N.I., Makarov A. A., Pilyar N.V., Rabtsun S.V., Smolkov R.A., Zhabitsky V.M., Baudrenghien P., Hofle W., Killing F., Ko-jevnikov I., Kotzian G., Louwerse R., Montesinos E., Rossi V., Schokker M., Thepenier E., Valuch D. LHC Transverse Feedback System: First Results of Commissioning // Proc. of the XXI Russian Particle Accelerators Conference Ru-PAC-2008, 28 September-3 October 2008, Zvenigorod, Russia. Moscow, 2008. Pp. 97-100.

24. Zhabitsky V. M. Beam Stability in Synchrotrons with Digital Filters in the Feedback Loop of a Transverse Damper // Proc. of the 11th International Conference on Heavy Ion Accelerator Technology, 8-12 June 2009, Venizia (Italy). Istituto Nazionale di Fisica Nucleare (INFN), 2009. Pp. 306-310.

25. Zhabitsky V.M., Hofle W., Kotzian G., Montesinos E„ Schokker M., Valuch D. Beam Tests of the LHC Transverse Feedback System // Proc. of XXII Russian Particle Accelerators Conference RuPAC-2010, 27 September-1 October 2010, Protvino, Russia. IHEP, Protvino, 2010. Pp. 275-279.

26. Zhabitsky V. M. Digital Longitudinal Feedback Systems in Synchrotrons // Proc. of XXII Russian Particle Accelerators Conference RuPAC-2010, 27 September-1 October 2010, Protvino, Russia. IHEP, Protvino, 2010. Pp. 224-226.

27. Жабицкий В. M. Применение Z-преобразования для решения задачи о подавлении свободных колебаний пучка в ускорителе. Сообщение ОИЯИ Р9-91-91. Дубна. 1991.—14 февраля. 18 с.

28. Жабицкий В. М., Иванов И. Н., Коренев И. Л., Юдин JI. И. Подавление рези-стивной неустойчивости в I ступени УНК при цифровой обработке сигнала в цепи обратной связи. Сообщение ОИЯИ Р9-91-99. Дубна. 1991, —22 февраля. 12 с.

29. Zhabitsky V. М. Theoretical Treatment of a Classical Transverse Feedback System Using Z-transform. JINR Commun. E9-91-156. Dubna. 1991.-9 April. 10 p.

30. Жабицкий B.M., Коренев И. JI., Юдин Л.И. Помехоустойчивость системы подавления когерентных поперечных колебаний пучка в синхротроне. Сообщение ОИЯИ Р9-91-246. Дубна. 1991.-31 мая. 14 с.

31. Жабицкий В. М., Коренев И. Л., Юдин Л. И. Система подавления резистивной

неустойчивости в I ступени УНК с БИХ-фильтром в цепи обратной связи. Сообщение ОИЯИ Р9-91-494. Дубна. 1991.-13 ноября. 11с.

32. Жабицкий В. М., Коренев И. Л., Юдин Л. И. Подавление резистивной неустойчивости сгруппированного пучка с использованием цепей обратной связи. Сообщение ОИЯИ Р9-92-309. Дубна. 1992, — 17 июля. 20 с.

33. Жабицкий В. М. Влияние ошибок в расстановке датчиков и толкателей на параметры системы подавления резистивной неустойчивости в I ступени УНК. Сообщение ОИЯИ Р9-92-310. Дубна. 1992.-17 июля. 14 с.

34. Zhabitsky V. М. Beam Stability in Synchrotrons with Digital Filters in the Feedback Loop of a Transverse Damper. JINR Commun. E9-2009-82. Dubna. 2009. — 2 June. lip.

35. Zhabitsky V. M. Transverse Damper System for LHC. SL-RFS Note 91-14, LHC Note 175. CERN, Geneva. 1991.-13 December. 16 p.

36. Zhabitsky V. M., Hofle W., Kotzian G. Beam Stability in Synchrotrons with Notch and All-Pass Filters in the Feedback Loop of a Transverse Damper. CERN-BE-2009-013. Geneva. 2009.-29 April. 12 p.

Цитированная литература

Ll. Hofle W. Progress in Transverse Feedbacks and Related Diagnostics for Hadron Machines // Proc. of 4th International Particle Accelerators Conference IPAC-2013, 12-17 May 2013, Shanghai, China. Physical Review Special Topics - Accelerators and Beams (PRST-AB), 2013. Pp. 3990-3994.

L2. Lonza M. Multi-bunch feedback system // CAS - CERN Accelerator School: Digital Signal Processing, 31 May - 9 June 2007, Sigtuna, Sweden. CERN-2008-003, Geneva, 2008. Pp. 285 - 330.

L3. Courant E.D., Snyder H. S. Theory of the Alternating-Gradient Synchrotron // Annals of Physics. 2000. Vol. 281. Pp. 360^180.

L4. Гоноровский И. С. Радиотехнические цепи и сигналы. М.: Советское радио, 1977. 608 с.

L5. Лебедев А.Н., Шальнов А. В. Основы физики и техники ускорителей. М.: Энергоатомиздат, 1991. 528 с.

L6. Ландау Л. Д., Лифшиц Е. М. Механика. М.: Наука, 1973. Т. I. 208 с.

L7. Эльсгольц Л. Э., Норкин С. Б. Введение в теорию дифференциальных уравнений с отклоняющимся аргументом. М.: Наука, 1971. 296 с.

L8. Боголюбов Н. Н. Собрание научных трудов в двенадцати томах. М.: Наука, 2005. Т. III. 605 с.

Получено 24 июня 2013 г.

Отпечатано методом прямого репродуцирования с оригинала, предоставленного автором.

Подписано в печать 24.06.2013. Формат 60 х 90/16. Бумага офсетная. Печать офсетная. Усл. печ. л. 2,1. Уч.-изд. л. 1,75. Тираж 100 экз. Заказ № 58014.

Издательский отдел Объединенного института ядерных исследований 141980, г. Дубна, Московская обл., ул. Жолио-Юори, 6. E-mail: publish@jinr.ru www.jinr.ru/publish/

 
Текст научной работы диссертации и автореферата по физике, доктора физико-математических наук, Жабицкий, Вячеслав Михайлович, Дубна

Объединенный институт ядерных исследований

На правах рукописи УДК 621.384.634.5

05201 450141

Жабицкий Вячеслав Михайлович

Динамика пучка в синхротронах с цифровыми системами подавления когерентных колебаний заряженных частиц

01.04.20 - Физика пучков заряженных частиц и ускорительная техника

ДИССЕРТАЦИЯ на соискание ученой степени доктора физико-математических наук

Дубна-2013

Содержание

Введение............................................................................6

Обзор литературы................................20

Глава 1. Физико-математическая модель цифровой широкополосной системы подавления когерентных поперечных колебаний......26

1.1. Вводные замечания и упрощения...................26

1.2. Дискретная модель системы подавления колебаний сгустка .... 30

1.3. ^-преобразование как метод решения................39

1.3.1. Основные определения....................39

1.3.2. Общее решение........................41

1.3.3. Характеристическое уравнение ...............42

1.3.4. Асимптотические решения..................43

1.3.5. Приближенные решения...................44

1.3.6. Идеальная система подавления................45

1.4. Реализации цепей обратной связи...................48

1.5. Отклик пучка на дельта-импульс...................51

1.6. Сравнение с экспериментальными результатами ..........53

1.7. Выводы и результаты..........................55

Глава 2. Допуски и ограничения для систем подавления когерентных поперечных колебаний сгустков......................58

2.1. Дискретная модель системы подавления колебаний сгустков ... 59

2.1.1. Базовые дифференциальные уравнения...........60

2.1.2. ^-преобразование как метод решения ...........62

2.2. Допуск на отклонение рабочей точки.................64

2.2.1. Основные соотношения и допущения............64

2.2.2. Карта изолиний для декрементов затухания........65

2.3. Отклик пучка на гармоническое воздействие............69

2.3.1. Резонансная кривая......................69

2.3.2. Коэффициент резонансного усиления............71

2.4. Подавление неустойчивости сгруппированного пучка.......73

2.4.1. Основные уравнения.....................73

2.4.2. Общее решение........................76

2.4.3. Одиночный сгусток......................78

2.4.4. Два идентичных сгустка...................82

2.4.5. Эквидистантные идентичные сгустки............84

2.4.6. Разброс частот когерентных колебаний сгустков......89

2.4.7. Подавление когерентных колебаний сгустков .......91

2.5. Сравнение с экспериментальными результатами ..........92

2.6. Выводы и результаты..........................95

Глава 3. Многоканальные системы подавления когерентных поперечных колебаний................................98

3.1. Особенность решений характеристического уравнения для классической системы вблизи полуцелого числа бетатронных колебаний Q0..................................99

3.1.1. Приближенное решение....................99

3.1.2. Сравнение с экспериментальными результатами......101

3.2. Дискретная модель многоканальной системы............102

3.2.1. Базовые уравнения ......................102

3.2.2. Характеристическое уравнение ...............105

3.2.3. Быстрая система подавления.................106

3.2.4. Сравнение с экспериментальными результатами......112

3.3. Выводы и результаты..........................113

Глава 4. Физико-математическая модель системы подавления коге-

рентных продольных колебаний......................114

4.1. Дискретная модель системы подавления...............114

4.1.1. Базовые уравнения ......................114

4.1.2. Система обратной связи с амплитудной модуляцией . . .120

4.1.3. Система обратной связи с фазовой модуляцией......120

4.1.4. Система обратной связи с фазовой модуляцией и дополнительной задержкой.....................122

4.1.5. Пример оценки параметров системы подавления.....123

4.2. Сравнение с экспериментальными результатами..........125

4.3. Выводы и результаты..........................126

Глава 5. Сглаженное приближение для описания динамики пучка в синхротроне с системой обратной связи.................128

5.1. Основные уравнения..........................128

5.2. МетодьТрешения............................131

5.2.1. Метод интегральных преобразований............131

5.2.2. Метод Крылова—Боголюбова—Митропольского......133

5.3. Системы обратной связи с цифровым фильтром ..........136

5.3.1. Демпфирование когерентных бетатронных колебаний . . .136

5.3.2. Демпфирование когерентных синхротронных колебаний . 139

5.4. Выводы и результаты..........................140

Глава 6. Системы подавления когерентных поперечных колебаний с нелинейной передаточной характеристикой цепи обратной связи . .141 6.1. Дискретная модель системы подавления с нелинейной передаточной характеристикой цепи обратной связи............142

6.1.1. Базовое разностное уравнение................142

6.1.2. Модификация метода Крылова—Боголюбова—Митропольского для решения разностных уравнений.......143

6.2. Примеры систем подавления с нелинейной передаточной характеристикой цепи обратной связи ...................146

6.2.1. Кубическая нелинейность...................146

6.2.2. Кулоновское трение......................148

6.2.3. Нелинейный режим с ограничением силы коррекции . . .149

6.3. О росте эмиттанса пучка в синхротроне вследствие ошибок ин-жекции при наличии системы подавления..............150

6.3.1. Физико-математическая модель...............151

6.3.2. Система подавления с линейной передаточной характеристикой ............................154

6.3.3. Система подавления с нелинейной передаточной характеристикой ...........................156

6.4. Выводы и результаты..........................159

Глава 7. Система подавления когерентных поперечных колебаний пучка в коллайдере LHC ............................161

7.1. Физико-технические спецификации системы подавления .....161

7.2. Основные параметры функционирующей системы подавления . . 164

7.3. Некоторые экспериментальные результаты..............174

7.4. Выводы и результаты..........................179

Заключение....................................180

Список сокращений и условных обозначений ...............185

Литература ....................................191

Приложение А. Некоторые математические соотношения........209

Введение

Представлены результаты теоретического исследования динамики пучка в синхротронах с системами подавления когерентных колебаний заряженных частиц на основе оригинальных моделей с целью уточнения имеющихся закономерностей и выявления ограничений, обусловленных использованием цифровых методов обработки сигналов в цепи обратной связи системы подавления. Выявлены количественные и качественные закономерности для сдвигов частот и декрементов затухания когерентных колебаний сгустков в синхротроне при наличии системы подавления, согласующиеся с экспериментальными данными, полученными на Большом Адронном Коллайдере (ЬНС) на протонном пучке.

Объектом исследования является сгруппированный пучок заряженных частиц, подверженный действию детерминированных возмущающих сил в синхротроне с системами обратной связи. Колебания частиц внутри сгустка, диффузионные явления, возбуждение колебаний случайными силами и радиационные эффекты, характерные для электронных синхротронов, не рассматриваются.

Диссертация выполнена в Объединенном институте ядерных исследований (ОИЯИ). В нее вошел цикл теоретических и экспериментальных исследований, проведенных автором в связи с созданием «Системы подавления когерентных поперечных колебаний пучка в ЬНС». Исследования осуществлялись в соответствии с Проблемно-тематическим планом научно-исследовательских работ и международного сотрудничества ОИЯИ, Соглашением о сотрудничестве между ОИЯИ и Европейской организацией ядерных исследований (ЦЕРН) по их участию в проекте по созданию Большого Адронного Коллайдера (ЬНС).

В диссертации осуществлено решение научно-технической проблемы, важной для физики пучков заряженных частиц и актуальной для нового поколения протонных и ионных синхротронов, которые создаются или совершенствуются с целью проведения эффективных исследований в области физики элементарных частиц и релятивистской ядерной физики.

Актуальность работы

Разработка, проектирование, модернизация и эксплуатация ускорительных комплексов протонов и ионов, таких как LHC (ЦЕРН), У-70 (ИФВЭ, Протвино), FAIR (GSI, Дармштадт), NICA (ОИЯИ), включающих в себя несколько синхротронов с высокими требованиями физического эксперимента к интенсивности и фазовой плотности пучка или коллайдер с предельно возможной светимостью, сопровождаются решением научно-технических задач ускорительной физики и техники, среди которых можно выделить пять проблем, относящихся к заявленной теме диссертации:

1. Достижение высокой интенсивности пучка заряженных частиц.

2. Сохранение фазовой плотности пучка при переводе из синхротрона-инжектора в принимающий синхротрон.

3. Обеспечение повторяемости расчетных параметров пучка в различных ускорительных циклах.

4. Контроль параметров синхротрона с требуемой точностью без заметного ухудшения качества пучка при длительных циклах ускорения и накопления.

5. Внедрение цифровых технологий в управляющей системе, обеспечивающей демпфирование когерентных колебаний сгустков.

Повышение интенсивности пучков заряженных частиц является одной из основных задач в ускорительной физике и технике. Одним из факторов ограничения интенсивности пучка в синхротронах и накопителях являются когерентные неустойчивости. Число обнаруженных к настоящему времени различных типов когерентных неустойчивостей весьма велико. Обусловленные ими ограничения на интенсивность ускоряемого пучка заставляют разрабатывать методы их изучения и подавления. С этой целью в современных синхротронах широко используются системы подавления (СП) когерентных колебаний пучка заряженных частиц [37]. В таких системах демпфирующий дефлектор осуществляет

коррекцию импульса пролетающих заряженных частиц пропорционально их отклонению от положения равновесия, измеренному в предшествующие моменты времени с помощью датчика положения центра тяжести пучка. Подобные СП позволили существенным образом повысить интенсивность ускоряемого пучка. В связи с этим в составе современных протонных синхротронов функционируют системы подавления когерентных колебаний, обеспечивающие, как это отмечено в литературном обзоре диссертации, декремент затухания колебаний а < 0.012. При этом актуален вопрос о требованиях к СП в зависимости от интенсивности ускоряемых сгустков и структуры ускоряемого пучка.

В современных ускорительных комплексах, в составе которых функционируют несколько синхротронов, при переводе сгустков из синхротрона-инжектора в принимающий синхротрон неизбежно возникают отклонения инжектируемых сгустков от расчетных положений их устойчивого равновесия. При этом инжектируемые заряженные частицы, циркулируя в синхротроне, совершают малые колебания вблизи своего положения устойчивого равновесия: бетатронные колебания в поперечном направлении (вертикальной и горизонтальной плоскостях) относительно мгновенной замкнутой орбиты (деформированной или расчетной) в ведущем магнитном поле и синхротронные колебания в продольном направлении относительно синхронной фазы. Из-за ошибок инжекции эти малые колебания частиц, вследствие их электромагнитного взаимодействия с другими частицами, могут инициировать совместные колебания выделенной группы частиц (сгустка) и становиться когерентными. Вследствие энергетического разброса частиц в инжектируемом пучке малые когерентные колебания переходят в некогерентные (процесс естественного затухания когерентных колебаний), что сопровождается увеличением линейных размеров сгустков. Происходит уменьшение фазовой плотности пучка. Так, на ускорительном комплексе ЦЕРН, в составе которого цепочка из четырех синхротронов, нормализованный поперечный эмиттанс пучка возрастает почти в 1.5 раза в основном из-за ошибок инжекции. В связи с этим необходимо подавлять когерентные колебания сгустков после ин-

жекции с высоким темпом за промежутки времени, которые должны быть короче времени перехода когерентных колебаний в некогерентные. Поэтому актуальны оценки предельных темпов подавления в контексте допусков на рост эмиттанса пучка в условиях двух конкурирующих процессов (естественного затухания и подавления когерентных колебаний).

Система подавления когерентных колебаний, будучи установленной на синхротроне и настроенной в соответствии с выбранной рабочей точкой по числу колебаний за оборот <2 о, должна обеспечивать стабильное подавление когерентных колебаний при вариации <2о- В литературе отсутствуют сведения о допусках на величину изменения <2о> которые необходимы при проектировании СП.

При исследовании динамики пучка в синхротроне при наличии СП используются дифференциальные уравнения малых колебаний [38]:

т0у

где £2 — циклическая частота колебаний заряженной частицы с массой то в фокусирующем электромагнитном поле, у — релятивистский фактор, т — величина задержки сигнала в цепи обратной связи. Рекомендации [38], получающиеся в сглаженном приближении при эмпирической замене импульсной силы /^(х) на диссипативную, пропорциональную скорости х, как минимум необходимо уточнить при использовании цифровой СП, поскольку результирующая сила воздействия дефлектора на пучок пропорциональна суперпозиции сигналов с датчика положения, зафиксированных на нескольких оборотах.

Основным направлением проектирования и развития СП является применение методов цифровой обработки сигналов в цепи обратной связи. Существующие предсказания по тактовой частоте цифровых систем подавления и оценки влияния шумов квантования на рост эмиттанса пучка не приводят к каким-либо ограничениям, обусловленным использованием цифровых технологий. Вместе с тем отсутствуют предсказания по возможным ограничениям на величины декрементов когерентных колебаний при использовании цифровых СП.

Основные параметры ускорительного цикла контролируются и регулируются с помощью цепей обратной связи в процессе ускорения. При измерении параметров пучка в ряде случаев используются методы импульсного гармонического воздействия на сгусток, например, с помощью электромагнитного толкателя. Эти воздействия могут быть систематическими (например, с целью контроля качества пучка после инжекции) или случайными (например, при измерении частот бета-тронных колебаний в произвольный момент времени в пределах ускорительного цикла). Амплитуда возбужденных когерентных колебаний должна быть достаточной по величине для детектирования сигналов и измерения параметров колебаний с требуемой точностью. Однако система подавления, обеспечивая затухание когерентных колебаний, приводит к неизбежному сокращению промежутка времени, в течение которого измеряемый сигнал может быть идентифицирован на уровне шумов. Поэтому необходимы оценки для величин и продолжительности действия сил возбуждения когерентных колебаний при наличии СП с учетом их возможного влияния на рост эмиттанса.

Актуальность перечисленных выше физических задач и поиск путей их решения явились мотивацией к выполнению данной работы.

Цель диссертационной работы:

- построение физико-математической модели, описывающей динамику пучка заряженных частиц после инжекции и при ускорении под влиянием управляющих воздействий цифровой системы обратной связи, приводящих к затуханию когерентных колебаний;

- получение аналитических решений, позволяющих сопоставить расчетные данные (например, декременты затухания когерентных колебаний) с экспериментально наблюдаемыми характеристиками колебательного процесса в максимально приближенных к эксперименту условиях, в том числе при темпах подавления, близких к предельно возможным;

- сопоставление результатов разработанной модели с получающимися в линейном приближении решениями уравнений динамики пучка при использовании метода Крылова—Боголюбова—Митропольского и методов решения дифференциальных уравнений с отклоняющимся аргументом;

- формулирование перечня спецификаций для физико-технических обоснований проектов цифровых систем подавления в новых синхротронах с учетом опыта создания такой системы для ЬНС (ЦЕРН) в совместном проекте ЦЕРН и ОИЯИ.

Научная новизна

Приведены результаты теоретического исследования динамики сгустков заряженных частиц в синхротронах с цифровыми системами обратной связи с помощью методов дискретного анализа:

• уточнены зависимости и выявлены ограничения для декрементов затухания и частот когерентных бетатронных колебаний от величины баланса фаз, численно равной разности между набегом фазы бетатронных колебаний частицы при пролете от измерительного датчика до демпфирующего дефлектора и изменением фазы управляющего сигнала в цепи обратной связи;

• впервые продемонстрировано, что при использовании цифровых фильтров в цепи обратной связи систем подавления с задержкой сигнала более чем на один оборот возникают дополнительные к бетатронным моды колебаний сгустков, которые влияют на область устойчивости движения частиц;

• введено новое понятие: сепаратриса для декрементов затухания когерентных поперечных колебаний; приведен метод расчета, позволяющий установить конфигурацию сепаратрисы и карты изолиний для декрементов затухания; показано, что выс