CC BY
14
Онлайн-доступ к журналу: http: / / mathizv.isu.ru
Серия «Математика»
2018. Т. 24. С. 68-81
УДК 519.6:621.384.6 MSG 65К10
DOI https://doi.org/10.26516/1997-7670.2018.24.68
Оптимизация динамики пучков заряженных частиц с использованием метода второго порядка
Д. А. Стариков
Санкт-Петербургский государственный университет, Санкт-Петербург, Российская Федерация
Аннотация. Задача оптимизации динамики пучка заряженных частиц рассматривается как задача программного управления ансамблем траекторий динамической системы с фиксированным конечным моментом времени. Состояние ансамбля динамических систем характеризуется плотностью распределения систем в фазовом пространстве, удовлетворяющей уравнению Лиувилля или уравнению Власова. Задача состоит в минимизации некоторого функционала, зависящего от конечного состояния ансамбля.
Предложено использовать алгоритм решения, основанный на вычислении первой и второй вариаций траектории динамической системы при вариации управления. В том случае, когда управление параметризовано, выражения для первой и второй вариации позволяют находить первые и вторые производные минимизируемого функционала по параметрам управления. Использование вторых производных может существенно ускорить процесс оптимизации по сравнению с алгоритмом, в котором используются только первые производные функционала.
Предложенный алгоритм реализован для пучка в канале с высокочастотной ква,-друпольной фокусировкой, который обычно используется в качестве начальной части ускорителей заряженных частиц. Рассмотрена простейшая задача оптимизации продольной динамики пучка в этом канале. Численное решение задачи проводится в рамках метода крупных частиц. Проведено сравнение эффективности методов первого и второго порядка. При этом метод второго порядка показывает существенное улучшение сходимости по сравнению с методом первого порядка.
Ключевые слова: оптимальное управление, ансамбль динамических систем, вторая вариация траектории, метод второго порядка, пучок заряженных частиц.
ОПТИМИЗАЦИЯ ДИНАМИКИ ПУЧКОВ 1. Введение
В работах Д. А. Овсянникова были разработаны методы численной оптимизации динамики пучков заряженных частиц, основанные на использовании первой вариации функционала, характеризующего качество управления, при вариации управления [6-8]. При этом задача оптимизации сформулирована как задача управления ансамблем динамических систем. Особенностью этой задачи является то, что ансамбль динамических систем характеризуется не только уравнениями динамики отдельной системы, но и плотностью распределения систем в фазовом пространстве, удовлетворяющей уравнению Лиувилля или уравнению Власова.
Численная реализация методов на основе первой вариации функционала может потребовать большого объема вычислений, в связи с чем и был предложен метод, основанный на использовании второй вариации функционала качества [4; 14; 15]. Актуальность разработки методов второго порядка для численного решения задач оптимального управления была отмечена также в работе [9].
Основную трудность при использовании метода первого порядка, представляет выбор направления и шага спуска, поскольку градиент функционала, который вычисляется в рамках метода первого порядка, определяет лишь гиперплоскость в многомерном пространстве параметров управления, при малом смещении вдоль которой значение функционала не меняется.
По сравнению с методом первого порядка, метод второго порядка позволяет находить не только первые, но и вторые производные функционала качества по параметрам, с помощью которых параметризовано управление. Направление и шаг спуска легко выбрать, если известны первые и вторые производные функционала.
Метод второго порядка, аналогичный представленному в настоящей работе, рассматривался также в работе [1], но там предлагалось находить вторые производные функционала качества как решения дифференциальных уравнений, число которых квадратично по числу параметров управления, то есть может быть достаточно велико.
2. Постановка задачи
Рассмотрим ансамбль управляемых динамических систем
f = f(t,x,u).
Здесь х - положение системы в фазовом пространстве Q, х € Q, dim Q = п, t — время, t € [to,T], а и = u(t) — управление, которое выбирается
из некоторого класса функций, например, кусочно-непрерывных. При этом предполагаем, что значения этих функций и(Ь) лежат в некотором компакте и, и(Ь) € и С Кт. Пусть системы распределены в некоторой области И фазового пространства или на некоторой движущейся поверхности Б в фазовом пространстве с плотностью, которую в обоих случаях будем описывать дифференциальной формой д(Ь,х) [3; 10-13], степень которой равна размерности фазового пространства п в первом случае и размерности поверхности Б во втором. Если в начальный момент ¿о задана начальная плотность распределения в области И или на некоторой начальной поверхности ¿>о, которую обозначим через до(х), д(Ьо,х) = до(%), то в обоих случаях плотность в последующие моменты времени может быть получена в виде [3; 10-13]
д(го + бг, Ffístx) = (2.1)
где Ff>st обозначает сдвиг Ли вдоль траекторий динамической системы на приращение параметра [2]. При этом, если взаимодействие между динамическими системами отсутствует, то уравнение (2.1) можно рассматривать как уравнение Лиувилля, а если системы взаимодействуют между собой, то уравнение (2.1) есть уравнение Власова.
Задача состоит в том, чтобы найти управление и из заданного класса, доставляющее минимальное значение функционалу
Ф (и) = I д(хт)д(Т,хт), (2.2)
к
характеризующему качество процесса управления. Здесь К обозначает И или Б в зависимости от размерности носителя распределения.
3. Численное решение
Уравнения для компонент вектора первой вариации х имеют вид
(15хг д/
(М с1хз
где
= -^-Ьх3 + 5ир, 6хг(г о) = 0, г = (3.1)
ди! дикди
(здесь и далее подразумевается суммирование по повторяющимся верхним и нижним индексам в соответствии с правилом Эйнштейна).
Решение системы неоднородных линейных уравнений (3.1) можно записать в виде
I
3
¿0
где 0(Ь, - функция Грина уравнения для первой вариации, представляющая собой тензор типа (|), компоненты которого удовлетворяют системе однородных уравнений
') ори) и , _
и условию = Е, где Е — единичная матрица. Кроме того,
компоненты матрицы Грина удовлетворяют также следующей системе уравнений
«КМ')
л- = (3'2)
Тогда компоненты вектора вариации траектории можно записать в виде
5х\г) = ! г = 1,Л, (3.3)
¿0
а первую вариацию функционала (2.2) в виде
т
5иФ = 11 ^С(Т,1')6иЦ1,х)е(1,х)сИ. (3.4)
П
Вводя сопряженную дифференциальную форму удовлетворяющую уравнению и условию
запишем первую вариацию функционала (3.4) в виде
т
6и Ф = -
П
Пусть и — кусочно-постоянная векторная функция £ :
1
■ J ! 1!)(1,,х)5иЛ1,,х)д(1,,х)<11. (3.5)
и = щ, г = 1, М, Ьм = Т.
Тогда функционал (2.2) можно рассматривать как функцию гМ параметров управления.
Из выражения (3.5) вытекает, что производные функционала по этим параметрам равны
дФ } [ ,д6и№,х)
дик
11ф^х)д6и^х)е(1,х)М. (3.6)
1 (о И
Далее применим метод крупных частиц, в рамках которого перейдем от непрерывной сплошной среды к дискретной (набор крупных частиц) и заменяем интегрирование по фазовым координатам суммированием по крупным частицам. Тогда минимизируемый функционал представим в виде
N
Ф(и) = 5>(ж0-)(Г)), (3.7)
3 = 1
а выражение для его производных по параметрам в виде
¿0
г Г ]=\
Здесь х^ — положение ¿-й крупной частицы в фазовом пространстве.
Для второй вариации траектории динамической системы при вариации управления нетрудно получить выражение [4; 14; 15]
г а
4о ¿0
I'
+
/ окт{И, г")5иГ(Г)) ей") + 5ЦЩ!)
t о ¿0
£й'. (3.8)
Здесь — тензор типа (2) , компоненты которого удовлетворяют
уравнениям
д£)г, (Ь Ь1) я;» я fт
и условиям
£>(*, £)}*; = 0, г^,к = 1~п.
Известия Иркутского государственного университета. 2018. Т. 24. Серия «Математика». С. 68-81
Выражая вариацию функционала через первую и вторую вариацию траектории в конечной точке, нетрудно найти вторые производные функционала (3.7) по параметрам:
д2Ф д(и^)2
<9Ф
дхз д{и:
д2Ф rdxj
d{xi)2ldv!y
(Г)]5
Q(T, ХТ),
(3.10)
где дхЗ/ди%(Т) и д2х^/д(и%)2(Т) определяются из выражений для первой и второй вариации (3.3) и (3.8) соответственно. Далее, в соответствии с методом крупных частиц, нетрудно перейти от интегрирования в (3.10) к суммированию по крупным частицам:
д2Ф
N г
9(Щ
к\2
Е
i=i
d2xL
—-—(Т) +
дхзд(и%)2К J
д2Ф д{хз)2
\—— (Т)]2 [ ди? { jj
(3.11)
4. Оптимизация ускорительной структуры с высокочастотной квадрупольной фокусировкой
В качестве примера применения представленного здесь метода рассмотрим задачу оптимизации ускорительной структуры с высокочастотной квадрупольной фокусировкой [5; 8; 16], которая требует значительного объема вычислений.
Ограничимся в этой работе только оптимизацией продольной динамики пучка. В простейшей модели канала с высокочастотной квадрупольной фокусировкой предполагается, что продольную компоненту электрического поля в некоторой ячейке можно представить в виде
z
4кТ f Ez = Uo-cos г? cos wi, rj(z) = / k(z') dz'. (4.1)
7Г J
zo
Здесь 2Uq — напряжение между электродами, со — частота колебаний поля, а — апертура этой ячейки (минимальное для этой ячейки расстояние от оси до электрода), к = тт/L, L — длина ячейки, которая меняется от ячейки к ячейке вдоль оси канала, rj(z) — фаза модуляции электродов, Т — эффективность ускорения. Таким образом, здесь мы рассматриваем k(z) как кусочно-постоянную функцию.
В рамках этой модели можно рассматривать продольное движение частиц независимо от поперечного.
В качестве координат в фазовом пространстве продольного движения рассмотрим приведенную энергию частицы 7 = 1 + W/(mc2) и ее фазу ip = cot, где t - тот момент времени, в который частица проходит
сечение с заданной продольной координатой z. Здесь W и т энергия и масса частицы, ас — скорость света. В качестве параметра траектории будем рассматривать безразмерную координату ( = z/X, где Л = 2ттс/ш.
Начальное распределение в фазовом пространстве зададим следующим образом. Пусть все частицы имеют одинаковую начальную энергию 7(Со) = 7о и равномерно распределены по начальной фазе <р((о) € [—27г,0]. Такое распределение описывается дифференциальной формой первой степени в фазовом пространстве поперечного движения, единственная компонента которой, заданная на моноэнергетической поверхности 7 = 7о, равна £>(0)^ = (27г)—1.
Уравнения продольной динамики для пучка малой интенсивности, учитывая выражение (4.1), запишем в виде
| = 2VT7(72 " 1)"1/2, (4-2)
d^y _
— = Ськи cos г] cos (р. (4.3)
ас,
Здесь ( = z/X, к = Хк, X = 2тгс/ш, Сь = 2еС/о/(7гтс2), а функция и(() = Т(Х() — единственная компонента вектора управления. При этом длины всех ячеек, а, следовательно, и к для каждой из ячеек мы считаем заданными.
В рамках метода крупных частиц каждой частице с номером г соответствуют две сопряженные функции — VV(¿) и Здесь г — номер крупной частицы, который записан в скобках в отличие от индексов компонент. Дифференциальные уравнения для этих функций имеют вид
dipip(i) — d(: = Сь ipj(i) kT cos r¡ sin ,
#7(¿) , 2тг
= %{%)■
¿С (7(2) ~~ I)3/2 '
Пусть управление — кусочно-постоянная функция: и(() = и^, ( € (0-1)0)) .1 = 1)^) гДе 0 границы ячеек, причем А(0 — 0-0 ~~ длина ячейки, а М — полное число ячеек. Тогда
= / ^>2'Ф1(г)Ськс08Г]С08<£^с1(.
3 /. г= 1 У-1
Для компонент функций Грина, которые используются при вычислении второй вариации и в уравнениях для компонент тензора И в
соответствии с уравнениями (3.2) имеем
= ощ.осьмо^^о^^с),
= 27Г(7(С02 - 1)-3/2С%((,(>),
^сс;с/)=27г(7(со2-1)-3/2^(С,С/)-
Здесь и далее мы опускаем номера крупных частиц при функциях Грина, фазовых переменных и аналогичных величинах для большей компактности приводимых выражений
Используя выражение (3.8), получим вторую вариацию фазы частицы при вариации управления в следующем виде
С
5М0 = I ^(С>С')МС,)СД(С,)со8(»7(С,))со8(^(С,))х
Со
с
X | СЩ1, (")5и((")Ськ((") сов^С")) соф«")) ¿С" ¿С'+
Со
С
+1 Щ<р((, ОЫОСьНС) (Мл«')) (-оФ(С)) х
Со
с
X I С^(С/,С//)^(С//)адС//)с08(г?(С//))с08(^(С//)Х<-
Со
с
-1 с^с,с/)^(с/)адс/)со8(г?(с/))8ш(^(с/))х
Со
с
X I С^(С/,С//)^(С//)адС//)с08(г?(С//))с08(^(С//)Х<-
Со
Аналогичный вид имеет и выражение для 527(С).
Запишем теперь дифференциальные уравнения для компонент тензора И, используя уравнения (3.9)
ЭД^,С,)=47Г(7(С,)2-1)-3/2Д^(С, С/)+6тг7(С/)(7(С/)2-1)-5/2^(С, СО-
дВ^,С) = СьЦ(') еМл«')) вш^О) ^(С, С0 + +2ТГ(7(С/)2-1)-3/2^(С,С/)-
гЮV (V -
= 2Ськ((') cos(í7(C,)) С0+
+Ськ((') со8(г?(С/)) со8(^(С'))С^(С> СО-=47г(7(С02-1)-3/2^(С, С0+67Г7(С/)(7(С/)2-1)-5/2^(С, СО-
гЮ1 (7 -
= ЗДСО сов(77(0) вшМС'ЯВЪК, С0+
+2тг(7(С/)2-1)-3/2^(С,С/)-
/9 Г)7 (V ^ _
+СД(С0 СО8(г?(С0) СО8(^(С0)Ц(С, СО-
Алгоритм численного решения состоит в том, что выбирается некоторое параметризованное управление, а затем оно улучшается в результате повторения следующих вычислений. Во-первых, находятся траектории крупных частиц и вычисляется функционал качества. После этого уравнения для динамических переменных, компонент функций Грина и компонент матриц £> интегрируются в обратном направлении, и вычисляются первые и вторые производные функционала качества по параметрам. Затем делается шаг, вектор которого к определяется этими производными. В настоящей работе шаг вычисляется следующим образом:
д2Ф <9Ф
при условии, что
д2Ф
аЩ > (45>
Если же условие (4.5) не выполнено, или же для вычисленного по формуле (4.4) значения оказывается невыполненным условие \к%\ < /¿¿о,
где Нго — некоторое максимальное значение \Нг\, то делается шаг абсолютной величиной Нго в направлении убывания функционала. Если же и дФ/дщ оказывается равным нулю в рамках некоторой точности, то на данном шаге полагается, что Нг = 0.
Схему (4.4) можно рассматривать как простейшую модификацию метода Ньютона
м д2Ф дФ
^ диьдщ " ,г° дип
3 = 1 '
При этом вместо матрицы вторых производных ||52Ф/с^с^|| используется матрица, в которой диагональные члены такие же, как и в матрице а все внедиагональные члены равны нулю. Применение такой упрощенной схемы обусловлено трудностями, связанными с обращением недиагональной матрицы высокого порядка.
Процесс останавливается при достижении требуемой точности.
Тестовые расчёты были произведены для ускорительной структуры с частотой 433 МГц, ускоряющей протоны от энергии 60 кэВ до 1850 КэВ.
На рис. 1 представлено распределение частиц по фазе и энергии на выходе из канала ускорения.
Рис. 1. Фазо-энергетический портрет пучка частиц на выходе из канала
ускорения.
На рис. 2 представлена зависимость значений функционала качества от номера шага оптимизационного процесса для методов первого и второго порядка с целью иллюстрации их сравнительной эффективности.
N. шаг
Рис. 2. Изменение функционала качества в процессе оптимизации.
Стартовые возмущения и, соответственно, оптимизация проводились вдоль всей структуры (точнее, от 22 ячейки и до конца структуры). Здесь линии А и В изображают процессы оптимизации методами первого и второго порядка соответственно.
5. Заключение
Видно, что метод второго порядка работает и дает существенное ускорение относительно метода первого порядка, даже в случае одной управляющей функции. Темп оптимизационного процесса для метода второго порядка очень высок на первых шагах, что компенсирует возросшее относительно метода первого порядка количество вычислений.
Таким образом, применение метода второго порядка для численного решения задач оптимального управления ансамблем динамических систем возможно и эффективно.
Список литературы
1. Бублик Б. Н., Гаращенко Ф. Г., Кириченко Н. Ф. Структурно-параметрическая оптимизация и устойчивость динамики пучков. Киев : Наукова думка, 1986. 304 с.
2. Дривотин О. И. Математические основы теории поля. СПб. : Изд-во С.-Петерб. гос. ун-та, 2010. 168 с.
3. Дривотин О. И, Овсянников Д. А. Решения уравнения Власова для пучка заряженных частиц в магнитном поле // Изв. Иркут. гос. ун-та. Сер. Математика. 2013. Т. 6, № 4. С. 2-22.
4. Дривотин О. И., Стариков Д. А. Применение методов второго порядка для численного решения задачи оптимального управления ансамблем динамических систем // Процессы управления и устойчивость. 2017. Т. 4, № 1. С. 113-117.
5. Капчинский И. M. Теория линейных резонансных ускорителей: Динамика частиц. M.: Энергоиздат, 1982. 240 с.
6. Овсянников Д. А.Математические методы управления пучками. Л. : Изд-во Ленингр. ун-та, 1980. 228 с.
7. Овсянников Д. А. Моделирование и оптимизация динамики пучков заряженных частиц. Л. : Изд-во Ленингр. ун-та, 1990. 312 с.
8. Овсянников Д. А., Дривотин О. И. Моделирование интенсивных пучков заряженных частиц. СПб.: Изд-во С.-Петерб. гос. ун-та, 2003. 176 с.
9. Федоренко Р. П. Приближенное решение задач оптимального управления. M. : Наука, 1978. 488 с.
10. Drivotin О. I. Covariant formulation of the Vlasov equation // Proc. of the 2011 Int. Particle Accelerators Conf., IPAC'2011. San-Sebastian, Spain, 2011. P. 2277-2279.
— accelconf.web.cern.ch/accelconf /ipac2011/papers/wepcll4.pdf
11. Drivotin O.I. Degenerate Solutions of the Vlasov Equation // Proc. of the 23rd Russian Particle Accelerators Conf., RUPAC'2012. St.-Petersburg, 2012. P. 376-378.
— accelconf.web.cern.ch/accelconf /rupac2012/papers/tuppb028.pdf
12. Drivotin О. I. Covariant Description of Phase Space Distributions// Vestnik of Saint-Petersburg University. Series 10. Applied mathematics. Computer science. Control processes. Issue 3. 2015. P. 39-52. https://doi.org/10.21638/11701/spbul0.2016.304
13. Drivotin О. I., Ovsyannikov D. A. Stationary Self-Consistent Distributions for a Charged Particle Beam in the Longitudinal Magnetic Field // Physics of Particles and Nuclei. 2016. Vol. 47, N 5. P. 884-913. https://doi.org/10.1134/S1063779616050038
14. Drivotin О. I, Starikov D. A. Second Order Method for Beam Dynamics Optimization // Proc. of the 24th Russian Particle Accelerators Conf., RUPAC'2014/ Obninsk, Oct. 2014. http://accelconf.web.cern.ch/AccelConf/ rupac2014/ papers/tupsal5.pdf
15. Drivotin О. I., Starikov D. A. Investigation of a Second Order Method of RFQ Channel Optimization // Proc. of the 25th Russian Particle Accelerators Conf., RUPAC'2016. St.-Petersburg, Nov. 2016. http://accelconf.web.cern.ch/AccelConf/ rupac2016/papers/tupsa013.pdf
16. Mathematical Models for Accelerating Structures of Safe Energetical Installation / О. I. Drivotin [et al.] // Proc. 6-th Europ. Part. Accel. Conf. EPAC'98. Stockholm, Sweden, 1998.
Стариков Дмитрий Александрович, аспирант, Санкт-Петербургский государственный университет, Российская Федерация, 198504, Санкт-Петербург, Петергоф, Университетский пр., 35, тел.: (812)4284726 (e-mail: radiumds@gmail.com)
Поступила в редакцию 30.04-2018
Charged Particle Beam Optimization with Use of a Second Order Method
D. A. Starikov
Saint-Petersburg State University, Saint-Petersburg, Russian Federation
Abstract. The problem of optimization of charged particle beam dynamics is considered. The problem is formulated as the control problem for a dynamical system ensemble with a fixed endpoint. A state of the dynamical system ensemble is described by a density of the systems in the phase space, which satisfies to the Liouville equation or to the Vlasov equation. The problem is to minimize a functional depending on terminal state of the ensemble.
It is proposed to use an algorithm based on calculation of the first and the second variations of trajectory of a dynamical system under the control function variation. If the control function is parametrized, expressions for the first and the second variations allow to find the first and the second derivatives of the functional being minimized over control parameters. Using of the second derivatives can make the optimization process sufficiently quicker, as compared with algorithm using only first derivatives.
The proposed algoritm is realized for a beam in the Radio Frequency Quadrupole (RFQ) channel, which is often used as initial part of a charged particles accelerator. The simplest problem of of optimization of longitudinal dynamics of the beam in this channel is considered. The numerical solution is finding on the base of the method of macroparticles. The comparison between the first order and the second order methods is conducted. The second order method shows sufficient increase of the rate of convergence as compared with the first order method.
Keywords: optimal control, dynamical system ensemble, second order variation, second order method, charged particle beam.
References
1. Bublik B.N., Garashchenko F.G., Kirichenko Н.Ф. Strukturno-parametricheskaya optimizatsiya i ustoychivost' dinamiki puchkov [Structural-Parametric Optimization and Stability of the Beam Dynamics]. Kiev, Naukova Dumka Publ., 1986, 304 p. (in Russian)
2. Drivotin O.I. Matematicheskie osnovy teorii polya [Mathematical Basics of the Field Theory] . Saint-Petersburg, Saint-Petersburg State University Publ., 2010, 168 p. (in Russian)
3. Drivotin O.I., Ovsyannikov D.A. Solutions of the Vlasov equation for a Charged Particles Beam in Magnetic Field. The Bulletin of Irkutsk State University. Series Mathematics [Izvestiya Irkutskogo Gosudarstvennogo Universiteta. Seriya Matematika], 2013, vol. 6, no. 4, pp. 2-22. (in Russian)
4. Drivotin O.I., Starikov D.A. Application of Second Order Methods for Numerical Solution of the Dynamical System Ensemble Optimal Control Problem. Protsessy upravleniya i ustoychivost', 2017, vol. 4, no. 1, pp. 113-117. (in Russian)
5. Kapchinskiy I.M. Teoriya lineynykh rezonansnykh uskoriteley: Dinamika chastits [Theory of Linear Resonant Accelerators: Particle Dynamics]. Moscow, Energoizdat Publ., 1982, 240 p. (in Russian)
6. Ovsyannikov D.A. Matematicheskie metody upravleniya puchkami [Mathematical Methods of Beam Control]. Leningrad, Leningrad University Publ., 1980, 228 p. (in Russian)
7. Ovsyannikov D.A. Modelirovanie i optimizatsiya dinamiki puchkov zaryazhennykh chastits. [Modeling and Optimization of Charged Particles Beam Dynamics. Leningrad, Leningrad University Publ., 1990, 312 p. (in Russian)
8. Ovsyannikov D.A., Drivotin O.I. Modeling of Intensive Charged Particle Beams. Saint-Petersburg, Saint-Petersburg University Publ., 2003, 176 p. (in Russian)
9. Fedorenko R.P. Priblizhennoe reshenie zadach optimal'nogo upravleniya [Approximate Solution of Optimal Control Problems]. Moscow, Nauka Publ., 1978, 488 p. (in Russian)
10. Drivotin O.I. Covariant formulation of the Vlasov equation. Proc. of the 2011 Int. Particle Accelerators Conf, IP AG'2011, San-Sebastian, Spain, 2011, pp. 22772279. accelconf.web.cern.ch/accelconf /ipac2011/papers/wepcll4.pdf
11. Drivotin O.I. Degenerate Solutions of the Vlasov Equation. Proc. of the 23rd Russian Particle Accelerators Conf, RUPAC'2012, St.-Petersburg, 2012, pp. 376378. accelconf.web.cern.ch/accelconf /rupac2012/papers/tuppb028.pdf
12. Drivotin O.I. Covariant Description of Phase Space Distributions. Vestnik of Saint-Petersburg University. Applied mathematics. Computer science. Control processes, 2015, Issue 3, pp. 39-52. https://doi.org/10.21638/11701/spbul0.2016.304
13. Drivotin O.I., Ovsyannikov D.A. Stationary Self-Consistent Distributions for a Charged Particle Beam in the Longitudinal Magnetic Field. Physics of Particles and Nuclei, 2016, vol. 47, no. 5, pp. 884-913. https://doi.org/10.1134/S1063779616050038
14. Drivotin O.I, Starikov D.A. Second Order Method for Beam Dynamics Optimization.Proc. of the 24th Russian Particle Accelerators Conf., RUPAC'2014, Obninsk, Oct. 2014.
http://accelconf.web.cern.ch/AccelConf/rupac2014/papers/tupsal5.pdf
15. Drivotin O.I., Starikov D.A. Investigation of a Second Order Method of RFQ Channel Optimization. Proc. of the 25th Russian Particle Accelerators Conf., RUPAC'2016, St.-Petersburg, Nov. 2016.
http://accelconf.web.cern.ch/AccelConf/rupac2016/papers/tupsa013.pdf
16. Drivotin O.I., et.al. Mathematical Models for Accelerating Structures of Safe Energetical Installation. Proc. of the 6-th Europ. Part. Accel. Conf. EPAC'98. Stockholm, Sweden, 1998.
Starikov Dmitriy Alexandrovich, Postgraduate, Saint-Petersburg State University, 35, Universitetskii pr., Petergof, St.-Petersburg, 198504, Russian Federation, tel.: (812)4284726 (e-mail: radiumds@gmail.com)
Received 30.04.2018
