Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Баев Теория колебаниы 2015

.pdf
Скачиваний:
25
Добавлен:
12.11.2022
Размер:
2.85 Mб
Скачать

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

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

Если в уравнениях бетатронных колебаний перейти от времени к

азимуту θ, то они принимают следующий вид:

3.3.26

 

 

 

 

 

]!! ν ] 0;

 

 

 

 

 

 

H!! ν H 0;

 

3.3.27

здесь

 

!

 

dz/dθ;

!

 

dx/dθ; и

относительные частоты

 

 

 

 

 

бетатронных колебаний,

выраженные

νв единицах циклотронной

 

]

 

 

 

H

 

ν

частоты:

 

 

 

ν √1 % ;

 

3.3.28

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

% .

 

 

3.3.29

Найденные выше условия устойчивости бетатронных колебаний (3.3.2) можно представить в графическом виде, т.е. в виде диаграммы

устойчивости. Из

ν

1,

3.3.30

выражений (3.3.28)

и (3.3.29) следует выражение

 

 

251

 

представляющее уравнение окружности единичного радиуса. Физический смысл имеет ее четверть, размещающаяся в первом квадранте, как показано на рис. 3.3.2

ν

1

0

1 ν

Рис. 3.3.2. Диаграмма устойчивости бетатронных колебаний.

Если точка, изображающая частицу на плоскости (ν , ν , лежит на окружности единичного радиуса, то ее бетатронные колебания устойчивы, если же нет, то бетатронные колебания такой частицы

неустойчивы.

 

В современных ускорителях магнитное

поле, за редким

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

этого вводится обобщенный азимут a как величина следующего вида:

¾ 2π Π ;

3.3.31

здесь s– координата, отсчитываемая вдоль траектории частицы;

 

периметр этой траектории; очевидно, для круговой

траектории

 

Π

обобщенным азимутом становится обычный азимут цилиндрической системы координат.

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

252

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

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

 

I!! ν I 0,

ν

 

 

3.3.32

однако в

уравнениях (3.3.32) коэффициент

уже

не

постоянная

величина, а функция обобщенного азимута Š

 

функции

 

 

 

 

свойство

 

 

 

Далее

следует отметить очень важное ¾ .

 

 

 

 

Š

 

После прохождения периода пространственной структуры

магнитного

 

ν

 

¾ .

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

 

Таким образом, уравнения бетатронных колебаний (3.3.32) можно

переписать в следующем более общем виде:

3.3.33

где

I!! ¾ I 0,

¾ (

¾ ;

3.3.34

 

здесь T – период пространственной структуры магнитного поля. Уравнение (3.3.33) с условием (3.3.34) является уравнением Хилла,

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

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

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

253

θ

T

r

 

g

Рис. 3.3.3. Кольцевой магнит с прямолинейными промежутками

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

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

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

полюсами магнита постоянно, а в прямолинейных промежутках равно

нулю. При такой аппроксимации функция

a будет иметь вид,

254

 

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

символ

 

m

 

длину участка равновесной орбиты под полюсом магнита.

В данном случае матрица периода MT

будет равна произведению

 

a

 

 

 

 

 

 

 

 

 

двух матриц:

<[ <~ · <,

 

 

 

3.3.35

где ~

 

 

 

 

 

 

 

и

матрица

магнита

и

матрица

прямолинейного

промежутка соответственно.

 

 

 

 

 

 

<

 

 

 

<

 

 

 

 

 

 

 

 

Построим сначала матрицу магнита. В пределах магнита

уравнение

Хилла

представляет

 

 

уравнение

гармонического

осциллятора, решение которого известно:

 

 

 

 

 

 

 

 

 

 

!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

3.3.36

 

 

 

 

I θ I

cos ν

θ ν

 

sin

θ ;

 

 

 

 

I! θ I νŠsin î I! cos ν θ ;

3.3.37

здесь I и I! начальные условия движения заряженной частицы.

ν, D

g

Sm

s

T

a

Рис. 3.3.4. Прямоугольная аппроксимация коэффициента

вуравнении Хилла

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

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

255

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

&

& &¾ &a

3.3.38

&θ &¾ &a &θ.

Объединим это выражение с выражением (3.3.31),

представив

периметр траектории в виде:

 

3.3.39

Π 2π 1 Λ ;

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

Тогда производная (3.3.37) преобразуется к виду

 

&

1

&

3.3.40

&θ 1 Λ &¾ ,

и, ведя ее в решения (3.3.36), (3.3.37) уравнения движения заряженных частиц в области магнита, будем иметь в них производные, вычисленные по обобщенному азимуту:

I θ I

cos ν θ

ν 1 Λ

sin ν θ ;

3.3.41

 

 

I!

 

 

I! θ I ν 1 Λ sin ν θ I! cos ν θ .

Теперь можно записать матрицу магнитного сектора

~

 

 

cos ν θ ,

 

sin ν θ

< ì

 

ν

 

1 Λ í .

 

 

 

 

 

 

 

 

 

ν

 

1 Λ sin ν

θ ,

cos ν θ

 

 

 

 

 

 

 

 

3.3.42

3.3.43

256

Матрица прямолинейного промежутка, длинной g , очевидно, имеет вид

 

 

 

1,

g

 

 

 

3.3.44

 

 

<

0,

1! ;

 

 

 

 

или, перейдя к обобщенному азимуту,

g

 

 

 

 

 

 

 

 

 

 

@1,

 

 

 

 

 

3.3.45

 

 

<

 

1 Λ A.

 

 

 

 

 

 

0,

 

 

1

 

 

 

 

 

 

 

<[

Перемножив матрицы (3.3.43) и (3.3.45), найдем матрицу

пространственного периода магнитного поля

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θ

 

 

 

¯

cos ν θ ,

 

 

 

g

 

sin ν

 

 

 

[

 

 

1 Λ cos ν θ ν

 

1 Λ ²

,

<

 

 

 

 

 

 

g

 

 

 

 

 

 

- ν

1 Λ sin ν

θ ,

 

 

cos ν θ ν

 

 

 

 

(3.3.46)

 

 

 

sin ν θ

 

°

 

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

 

1

ν

g

 

3.3.47

1 N cos ν

θ 2

sin ν θ N 1.

Распишем эти условия для каждого вида бетатронных колебаний отдельно. Условия устойчивости вертикальных бетатронных колебаний:

1

g

3.3.48

1 N cos √%θ 2

% a$% √%θ N 1;

условия устойчивости радиальных бетатронных колебаний:

257

1 N cos √1 %θ 12 1 % g sin √1 %θ N 1. 3.3.49

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

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

3.4. Уравнение огибающей пучка

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

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

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

258

& H 5 H 0,
&

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

Ниже речь пойдет об одном из таких уравнений, которых сегодня известно несколько [19]. Это уравнение огибающей Владимирского Капчинского или ВК-уравнение, описывающего поведение границы пучка.

Запишем уравнение движения частицы относительно равновесной орбиты

3.4.1

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

Представим его общее решение в виде (2.2.43), тогда с помощью

выражения (2.2.44) можно получить следующее выражение:

 

&H

3

&|χ|

3

3.4.2

&

&

cos ψ θ |χ| sin ψ θ .

Освобождаясь в выражениях (2.2.43) и (3.4.2) от тригонометрических функций, найдем один из интегралов уравнения

движения, который

представляет уравнение эллипса на

фазовой

плоскости &H/& в фиксированный

момент времени

 

&|χ|

 

H

 

 

&H

-

3.4.3

3

,|χ| & H

&

|χ|! .

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

подберем значение константы

 

, а также начальные

условия для

?|χ|/?

 

 

 

259

 

 

и

модуля

фундаментальных

решений уравнения движения

 

 

3

 

 

 

 

 

 

 

 

 

начальный

 

 

таким образом, чтобы эллипс (3.4.3) в

 

|χ 0 |

 

момент времени совпал с границей фазового объема пучка, изображенного на рис. 3.4.1.

При этом оказывается [13], что интеграл уравнения движения (3.4.3) определяется длинами N и M полуосей эллипса, изображенного на рис. 3.4.1, которым аппроксимируется эмиттанс пучка

3 3 √< .

 

3.4.4

&H

 

 

 

&

<

 

 

 

H

 

 

 

 

Рис. 3.4.1. Эмиттанс пучка

 

Поскольку всем частицам,

находящимся на границе

фазового

объема пучка, должно соответствовать полученное из начальных

условий значение интеграла уравнения

 

движения

 

 

 

,

то и в

любые последующие моменты времени, а не

только

в начальный

 

3 3

 

 

момент, эллипс (3.4.3) при условии, что

 

3 3 , должен оставаться

границей фазового объеме пучка.

 

 

cos ψ θ 1,

3 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Частица с

 

 

оказывается на

границе

пучка,

когда

 

 

 

поэтому

 

ρ 3 |Ê |,

 

 

 

 

 

3.4.5

 

ρ

 

 

 

 

 

 

 

 

 

где

– расстояние от равновесной орбиты до границы пучка.

 

 

Перейдя к уравнению (2.2.44) для

модуля

|χ |

,

получим ВК-

уравнение огибающей:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

ρ

5 ï

3

 

0

 

 

 

 

 

3.4.6

 

 

 

 

&

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

260