Научная статья на тему 'Оптимизация динамики пучков траекторий c использованием гладких и негладких функционалов. Часть 1'

Оптимизация динамики пучков траекторий c использованием гладких и негладких функционалов. Часть 1 Текст научной статьи по специальности «Физика»

CC BY
56
11
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
УПРАВЛЯЕМАЯ ДИНАМИЧЕСКАЯ СИСТЕМА / АНСАМБЛЬ ТРАЕКТОРИЙ / ГЛАДКИЙ ФУНКЦИОНАЛ / НЕГЛАДКИЙ ФУНКЦИОНАЛ / ПРИНЦИП МАКСИМУМА / ПУЧОК ЗАРЯЖЕННЫХ ЧАСТИЦ / УСКОРИТЕЛЬ / OPTIMAL CONTROL / CONTROLLED DYNAMIC SYSTEM / TRAJECTORY ENSEMBLE / SMOOTH FUNCTIONAL / NON-SMOOTH FUNCTIONAL / MAXIMUM PRINCIPLE / CHARGED PARTICLE BEAM / ACCELERATOR

Аннотация научной статьи по физике, автор научной работы — Овсянников Дмитрий Александрович, Мизинцева Мария Александровна, Балабанов Михаил Юрьевич, Дуркин Александр Павлович, Едаменко Николай Семенович

Проблемам управления и оптимизации в динамических системах посвящено много различных работ. Интерес к этим задачам со временем не уменьшается. Возникают новые задачи при разработке технологических процессов в разнообразных областях науки и техники, в частности при проектировании и создании современной электрофизической аппаратуры. В данной работе рассматриваются проблемы оптимизации и управления пучками траекторий. Задачу совместной оптимизации программного движения и пучка возмущенных движений предлагается решать с помощью комбинирования гладких и негладких функционалов. В ч. 1 описывается математическая постановка задачи, дается представление вариации комбинации функционалов и формулируются условия оптимальности в форме принципа максимума. Использование гладких и негладких функций позволяет определить функционалы, наиболее точно отражающие требования к динамике пучка заряженных частиц в ускорителях. В ч. 2 будут приведены результаты применения предложенной в ч. 1 методики для оптимизации динамики заряженных частиц в ускорителе с пространственно-однородной квадрупольной фокусировкой.Many different works are devoted to the problems of control and optimization in dynamic systems. Interest in these tasks does not decrease with time. New challenges arise in the development of technological processes in various fields of science and technology, in particular, in the design and creation of modern electrophysical equipment. In this paper, the problem of optimization and control of trajectory beams is considered. The problem of joint optimization of the program motion and the beam of perturbed motions using a combination of smooth and nonsmooth functionals is investigated. The first part deals with the mathematical formulation of the given representation of the variation of investigated functional and provides optimality conditions in the form of a maximum principle. In the second part, the problem of optimization of dynamics of charged particles in the accelerator with spatially homogeneous quadrupole focusing will be considered. Using a combination of smooth and non-smooth functions allows you to set the functionals that most accurately reflect the requirements for the dynamics of the charged particle beam in accelerators.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Оптимизация динамики пучков траекторий c использованием гладких и негладких функционалов. Часть 1»

2020 ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА Т. 16. Вып. 1 _ПРИКЛАДНАЯ МАТЕМАТИКА. ИНФОРМАТИКА. ПРОЦЕССЫ УПРАВЛЕНИЯ_

ПРОЦЕССЫ УПРАВЛЕНИЯ

УДК 517.977 МБС 49Л5

Оптимизация динамики пучков траекторий е использованием гладких и негладких функционалов. Часть 1

Д. А. Овсянников1, М. А. Мизинцева1, М. Ю. Балабанов1, А. П. Дуркин2, Н. С. Едаменко1, Е. Д. Котина1, А. Д. Овсянников1

1 Санкт-Петербургский государственный университет, Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9

2 Институт ядерных исследований РАН, Российская Федерация, 117312, Москва, пр. 60-летия Октября, 7а

Для цитирования: Овсянников Д. А., Мизинцева М. А., Балабанов М. Ю., Дуркин А. П., Едаменко Н. С., Котина Е. Д., Овсянников А. Д. Оптимизация динамики пучков траекторий с использованием гладких и негладких функционалов. Часть 1 // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2020. Т. 16. Вып. 1. С. 73-84. https://doi.org/10.21638/11702/spbu10.2020.107

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

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

1. Введение. Обе части работы посвящены проблемам управления и оптимизации динамики пучков (ансамблей) траекторий. Задачи управления и оценки состоя© Санкт-Петербургский государственный университет, 2020

ния динамических систем в условиях неопределенности часто трактуются как задачи управления при неполной информации о начальных данных и внешних возмущениях [1]. Отметим также, что данные проблемы тесно связаны с задачами оптимизации динамики пучков заряженных частиц в ускорителях [2, 3]. При этом математическая теория управления пучками траекторий развивается в основном в двух направлениях: изучаются или гладкие, или негладкие функционалы. Многие исследования посвящены различным аспектам такой теории (см. [4-26]).

В настоящей работе рассматривается проблема совместной оптимизации программного и возмущенных движений с использованием комбинации гладких и негладких функционалов [27-30]. В ч. 1 описывается математическая постановка задачи, дается представление вариации комбинированного функционала и приводятся условия оптимальности в форме принципа максимума Понтрягина [31]. Функционал состоит из трех частей: программное движение оценивается функционалом Больца, а динамика пучка возмущенных движений — как гладкими, так и негладкими функциями. Гладкий функционал характеризует в среднем динамику всего пучка и его выходные параметры, негладкий функционал — наихудшие выходные параметра пучка и представляет в отдельности минимаксную задачу. Применение комбинации гладких и негладких функций дает возможность, в частности при проектировании, создании и эксплуатации ускорителей заряженных частиц, строить эффективные математические модели оптимизации динамики пучка с использованием функционалов, наиболее точно отражающих требования к разнообразным характеристикам ускоряемого пучка.

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

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

^ = f(t,x,u), ж(0) = ж0, (1)

= F{t, х, у, и), у(0)=уоеМо, (2)

at

где xo и уо — начальные условия; t G To = [0, T] С R1 — независимая переменная (как правило, время); x G Rn и у G Rm — фазовые векторы размерности n и m соответственно; u G Rr — r-мерная вектор-функция управления; T — фиксированное значение. Вектор-функции f (t,x,u) и F(t,x,y,u) считаются достаточно гладкими. Множество Mo — компакт ненулевой меры.

Совместно с системой (1), (2) для плотности распределения частиц [2, 3] д = g(t) = g(t,y(t)) на траекториях подсистемы (2) рассмотрим уравнение

^ = -д divyF(t, х, у, и) (3)

с законом распределения плотности до(уо), заданным на начальном множестве Mo:

д(0) = д(0,у(0)) = до (у о), уо G Mo. Здесь ди(ус!) — некоторая непрерывная неотрицательная функция.

Будем считать, что допустимые элементы управления и = и (г), г € То, составляют класс Б кусочно-непрерывных на [0, Т] вектор-функций со значениями в компактном множестве и С Нг. Также предположим, что система (1), (2) имеет единственное решение задачи Коши с фиксированным хо и произвольным уо € Мо на интервале То для всех допустимых управлений и € О.

Решение подсистемы (1) будем называть программным (или выбранным/рассчитанным) движением, решения подсистемы (2) с начальными условиями уо € Мо и фиксированным программным движением — возмущенными движениями.

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

Перейдем к математической постановке задачи одновременной оптимизации. На решениях системы (1), (2) и уравнения (3) с заданными начальными условиями и управляющей вектор-функцией и (г) введем следующие функционалы:

т

Ь (и) = I Р1(г,х(г)) Л + д1(х(Т)), (4)

о

т

12(и) = (г)) сИ + (Т)), (5)

о

1з(и)= тах {дз(хт,Ут) • ё(Т,ут)}, (6)

ут еМт,и

где

т1(г) = у ф2(г,х(г),уг)в(г,уг)¿г, Ы2(Т) = ! д2(ут)в(Т,Ут) ¿г.

Мт,и

Здесь множество М^и — поперечное сечение пучка траекторий подсистемы (2), исходящих из множества Мо при заданном управлении и(г) и соответствующем программном движении х(г) в момент г, функции Ф,С, у2, д\,д2, дз являются неотрицательными непрерывно дифференцируемыми функциями своих аргументов. Введем функционал

I (и)= ь(и)+ 12(и) + 1з (и), (7)

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

Задачу минимизации функционала (7) на элементах управления и € О будем называть задачей совместной оптимизации программного движения и возмущенных движений.

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

3. Вариация функционала. Предположим, что и и и + Ди — допустимые управления. Тогда вариация 61(и, Ди) функционала (7) имеет вид

61(и, Ди) = 61\(и, Ди) + 512(и, Ди) + 61з(и, Ди). Вариации 61\(и, Ди) и 612(и, Ди) функционалов (4), (5) известны [14, 16]:

г/1(и>Ди)= ]Щ!Ш5ф)л+ЩШ5х(т), (8)

,] дх ох

о

т

зЫиЛи) = 1'ф'Ы^)) I вы +

Мь.

+ + I 9тд-^6у(Т)3,ут. (9)

мт,

Вариации 6х и 6у для систем (1) и (2), входящие в равенства (8) и (9), удовлетворяют

следующим уравнениям [14, 31]:

^ = + (10)

¿¿у = дР(г,х,у,и)5х + дР(г,х,у,и) +

аЬ дх дх

и начальным условиям

6х(0) = 0, 6у(0) = 0.

В (10), (11) и далее символ Ди обозначает приращение функции при изменении одной только переменной и, например Ди/ = /(Ь, х,и + Ди) — /(Ь, х, и).

Можно показать, что вариация функционала (6), т. е. 61з(и, Ди), имеет вид

613(и,Аи)= шах \в(Т,уо)(^Щ^^6х(Т) +

уоеПт,и у \ дх

+ +д3(х(Т),у(Т))6в(Т,у0)]1 . (12)

При этом 6д в соответствии с уравнением (3) для плотности распределения частиц д удовлетворяет уравнению в вариациях

dSg ( д divy F д divy F

df ~e ^ Sx + —Q^Sy + A« divy F ) - Sg divy F (13)

с начальным условием Sg(0) = 0.

Принимая во внимание соотношение [3, 16]

d divy Sy д divy F д divy F _

можно переписать уравнение (13) следующим образом:

¿5д . . С divy 5у

—— = —од dlvu Р и, х.у.и)--г2-.

Л и у ,у, > л

Интегрируя уравнение для вариации функции плотности распределения с учетом условий divy 6у(0) = 0 и 6д(0) = 0, получаем равенство

6д = —д divy 6у.

Введем множество

Rt,u = |у0 : Уо € М0, (д3(х(Т), у(Т, х0, у0, и)) &(Т, у(Т, х0, у0, и))) = = ma„X {#з(х(Т),y(T,xo,yo,u)) g(T,y(T,xo,yo,u))}\.

yoEMo J

Теперь вариацию SI%(u, Au) (равенство (12)) функционала (6) можно представить так:

61з(и,Аи)= max ig(T,y0) ('6х{Т) +

уоея.т,и I \ дх

+ уо)-д3(хт, У(Т, уо)) divy 6у(Т, у0)

Введем сопряженные функции х(£,уь), и(^,уг), фт(г,уг), Ит^,уг), "т(г,уг), удовлетворяющие вдоль траекторий (1)-(3) следующим линейным интегро-дифферен-циальным уравнениям:

dx df* dip\ [OF* ■ =—

q~dx dyt'

dt дх дх J дх J дх

Mt,u Mt,u

d/л (дF , ^ Д , , д^*2

. —

йфт df* OF* д divy F*

'ym о Mm о ^m:

dt дх дх дх

d/j,m _ OF* д divy F*

7, о Mm Г; ^m:

dt дy дy

dt ~

и конечным условиям

х< т,Ы = -а«МХ))

^(T,yo) = -G'(w2) ■ дт Фт = -QT

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

дх

дд2 (ут )*

ду

ддз^(T), y(T, yo))*

fm(T, Уо) = -вТ

дх

ддз(хт ,Ут )*

Здесь

ду

Vm(T, Уо)=дз(хт,Ут)вт, Уо е Вт,u, yo е Mo.

дт = д(Т,ут), вт = д(гГ,ут), dG(w2)

G'(w2) =

) =

dw2 ' d^(w1)

сСм1

4. Вариация функционала и условие оптимальности. Предположим, что и(г) и и(г) — допустимые управления. Введем обозначение Аи(г) = и(г) — и(г). Тогда вариация функционала (7) может быть представлена в виде

SI(u, Au)= max J - f ((x + фт)* ■ Auf + f (f* ■ AuF)dyt +

увеПт,и 7 v J

K 0 Mt,u

+ fa ■ AuF + V*m ■ Au divy F) dt^ .

Введем функции Hi и H2:

Hl(t, х, y, X, Фт, fm, Vm, u) = (x + Фт)* ■ f (t, х, u)+tfn ■ F (t, х, y, u) + + Vm ■divy F(t, х, y, u), H2 (Ь,х,у, f,,v)= ц* ■ F(t,х,у,ь,). Тогда вариацию (14) функционала (7) можно переписать так:

SI(u, Au) = max < - [AuHi^,х^),у^,уо),х^,Уо),Фт(Ь,Уо),^т(Ь,Уо), yoERT,u I J v о

(14)

vm(t,yo),u(t))+ J AuH2(t,x(t),yt ,^(i,yt),u(t)),dy^j dt>. (15)

xu

Mt,u

Вариации (14) и (15) не являются классическими вариациями функционала (7), так как в общем случае они нелинейны по Au. Такие вариации используются в математической теории управления для получения условий оптимальности в форме принципа максимума Понтрягина.

Теорема. Предположим, что u° = uo(t) есть оптимальное управление, т. е. при этом управлении функционал (7) достигает минимума. Тогда для t G To = [0, T] выполняется равенство

max max H°(t,yo,u) = max H°(t,yo,u0(t)),

uEU yo€RTuo yo€RTuo

где

h o{t, yo,u)=hi (t, x°(t), y0t yo), x0t y o), Фт(t, y o), ^L(t, y o), (t, yo),u)+

+ JH2(t,xo(t),y?,^o(t,y?),u) dy0t.

При этом функции х0(Ь), у0(Ь,уо), х°(Ь,Уо), Фт(Ь,Уо), ^0т(Ь,Уо), ^(Ь,уо), У0,^О(Ь,Уг) соответствуют оптимальному управлению ио(Ь) и уо € ИТ,ио, у° € Мг,ио (Ь,у°).

5. Заключение. Вариация (15) функционала (7) и условия оптимальности позволяют строить направленные методы минимизации этого функционала. Отметим, что при совместной оптимизации программного движения и ансамбля возмущенных движений, т. е. при оптимизации функционала (7), динамика пучка траекторий влияет на программное движение. Математически это усматривается в наличии множителя (х + Фт) при приращении Ди/ под знаком первого из интегралов в правой части равенства (14), т. е. влияют не только интегральные характеристики пучка, но и «наихудшие» частицы, доставляющие максимум функционалу (6). При этом учитывается плотность распределения частиц в пучке.

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

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

Следует также отметить эффективность применения методики оптимизации динамики пучков траекторий при решении разнообразных задач электрофизики и ядерной медицины [33-35].

Mtu0

Литература

1. Куржанский А. Б. Управление и наблюдение в условиях неопределенности. М.: Физматлит, 1977. 394 с.

2. Овсянников Д. А. Математические методы управления пучками. Л.: Изд-во Ленингр. ун-та, 1980. 228 с.

3. Овсянников Д. А. Моделирование и оптимизация динамики пучков заряженных частиц. Л.: Изд-во Ленингр. ун-та, 1990. 312 с.

4. Ананьина Т. Ф. Задача управления по неполным данным // Дифференциальные уравнения. 1976. Т. 13. № 4. С. 612-620.

5. Пантелеев А., Семенов В. В. Синтез оптимальных систем управления при неполной информации. М.: Изд-во Моск. авиац. ин-та, 1992. 191 с.

6. Бортаковский А. С. Оптимальное и субоптимальное управления пучками траекторий детерминированных непрерывно-дискретных систем // Изв. Рос. акад. наук. Теория и системы управления. 2009. № 1. С. 18-33.

7. Kurzhanski A. B., Varaiya P. Dynamic optimization for reachability problems // Journal of Optimization Theory and Applications. 2001. Vol. 108. P. 227-251.

8. Kotina E. D. Discrete optimization problem in beam dynamics // Nucl. Instrum. Methods Phys. Res. A. 2006. Vol. 558. N 1. P. 292-294.

9. Altsybeyev V. V., Ovsyannikov D. A. Optimization of beam parameters in APF channel // 27th International Linear Accelerator Conference, LINAC 2014. Proceedings. Geneva, 2014. P. 722-725.

10. Ovsyannikov D. A. Mathematical methods of optimization of charged partical beams dynamics // Proceeding of the European Partical Accelerator Conference. Barselona, Spain, 1996. Vol. 2. P. 1382-1384.

11. Bondarev B., Durkin A., Ovsyannikov A. New mathematical optimization models for RFQ structures // Proceedings of the Particle Accelerator Conference. New York City. 1999. Vol. 4. P. 28082810.

12. Виноградова Т. К., Демьянов В. Ф. О принципе минимакса в задачах оптимального управления // Докл. АН СССР. 1973. Т. 213. № 3. С. 512-514.

13. Демьянов В. Ф., Виноградова Т. К., Никулина В. Н. и др. Негладкие задачи теории оптимизации и управления. Л.: Изд-во Ленингр. ун-та, 1982. 324 с.

14. Овсянников А. Д. Об одном классе задач оптимизации в электростатическом поле // Докл. АН. 2013. Т. 453. № 4. С. 383-385.

15. Овсянников А. Д. Управление программным и возмущенными движениями // Вестн. С.-Петерб. ун-та. Сер. 10. Прикладная математика. Информатика. Процессы управления. 2006. Вып. 4. С. 111-124.

16. Овсянников А. Д. Управление пучком заряженных частиц с учетом их взаимодействия // Вестн. С.-Петерб. ун-та. Сер. 10. Прикладная математика. Информатика. Процессы управления. 2009. Вып. 2. С. 82-92.

17. Овсянников А. Д. Математические модели оптимизации динамики пучков. СПб.: Изд-во «ВВМ», 2014. 181 с.

18. Ovsyannikov D. A., Altsybeyev V. V. On the beam dynamics optimization problem for an alternating-phase focusing linac // Phys. Part. Nuclei Lett. 2016. Vol. 13. N 7. P. 775-779.

19. Ovsyannikov D. A., Ovsyannikov A. D., Vorogushin M. F., Svistunov Yu. A., Durkin A. P. Beam dynamics optimization: models, methods and applications // Nucl. Instrum. Methods Phys. Res. A. 2006. Vol. 558. N 1. P. 11-19.

20. Ovsyannikov A. D., Durkin A. P., Ovsyannikov D. A., Svistunov Y. A. Acceleration of different ion types in single RFQ structure // Problems of Atomic Science and Technology. 2016. Vol. 3(103). P. 54-56.

21. Овсянников Д. А., Едаменко Н. С. Моделирование динамики пучков заряженных частиц // Вестн. С.-Петерб. ун-та. Сер. 10. Прикладная математика. Информатика. Процессы управления.

2013. Вып. 2. С. 60-65.

22. Ovsyannikov A. D., Ovsyannikov D. A., Altsybeyev V. V., Durkin A. P., Papkovich V. G. Application of optimization techniques for RFQ design // Problems of Atomic Science and Technology.

2014. Vol. 3(91). P. 116-119.

23. Ovsyannikov A. D., Ovsyannikov D. A., Chung S.-L. Optimization of a radial matching section // International Journal of Modern Physics A. 2009. Vol. 24(5). P. 952-958.

24. Ovsyannikov A. D., Ovsyannikov D. A., Chung S.-L. Optimization of matching section of an accelerator with a spatially uniform quadrupole focusing // Technical Physics. The Russian Journal of Appl. Phys. 2009. Vol. 11. P. 1663-1666.

25. Ovsyannikov D. A., Ovsyannikov A. D., Antropov I., Kozynchenko V. BDO-RFQ code and optimization models // Proceedings of the International Conference of Physics and Control. 2005. P. 282288.

26. Ovsyannikov D. A., Mizintseva M. A., Ovsyannikov A. D. Joint optimization of smooth and nonsmooth functionals on beams of trajectories // Оптимальное управление и дифференциальные игры: Материалы Междунар. конференции, посвященной 110-летию со дня рождения Льва Семеновича Понтрягина. М.: Изд-во Матем. ин-та им. В. А. Стеклова Рос. акад. наук, 2018. С. 203-205.

27. Balabanov M. Yu., Mizintseva M. A., Ovsyannikov D. A. Beam dynamics optimization in a linear accelerator // Вестн. С.-Петерб. ун-та. Прикладная математика. Информатика. Процессы управления. 2018. Т. 14. Вып. 1. С. 4-13.

28. Мизинцева М. А. Об одной задаче минимизации гладких и негладких функционалов //

Вестн. С.-Петерб. гос. ун-та технологии и дизайна. Сер. 1. Естественные и технические науки. 2018. № 2. С. 37-41.

29. Mizintseva M., Ovsyannikov D. Minimax problem of simultaneous optimization of smooth and non-smooth functionals // Constructive non-smooth analysis and related topics: Proceedings of the Conference CNSA 2017. IEEE. 2017. P. 218-222.

30. Mizintseva M., Ovsyannikov D. On the minimax problem of beam dynamics optimization // 25th Russian Particle Accelerator Conference, RuPAC 2016. Saint Petersburg, 2016. P. 360-362.

31. Понтрягин Л. С., Болтянский В. Г., Гамкрелидзе Р. В., Мищенко Е. Ф. Математическая теория оптимальных процессов. М.: Наука, 1969. 384 с.

32. Капчинский И. M. Теория линейных резонансных ускорителей. М.: Энергоиздат, 1982. 239 c.

33. Golovkina A., Ovsyannikov D., Olaru S. Performance optimization of radioactive waste transmutation in accelerator driven system // Cybernetics and Physics. 2018. Vol. 7. N 4. P. 210-215.

34. Ovsyannikov D., Zavadskiy S. Pareto-optimal choice of controlled dimension for plasma stabilization system // IFAC-Papers OnLine. Elsevier Science Publishing Company Inc. elSSN: 24058963. 2018. Vol. 51. N 32. P. 175-178.

35. Bazhanov P., Kotina E., Ovsyannikov D., Ploskikh V. Optimization algorithm of the velocity field determining in image processing // Cybernetics and Physics. 2018. Vol. 7. N 4. P. 174-181.

Статья поступила в редакцию 2 марта 2019 г. ^атья принята к печати 13 февраля 2020 г.

Контактная информация:

Овсянников Дмитрий Александрович — д-р физ.-мат. наук, проф.; d.a.ovsyannikov@spbu.ru

Мизинцева Мария Александровна — магистр, ассистент; m.mizintseva@spbu.ru

Балабанов Михаил Юрьевич — канд. физ.-мат. наук, ст. науч. сотр.; m.balabanov@spbu.ru

Дуркин Александр Павлович — канд. техн. наук, доц.; durkinap@mail.ru

Едаменко Николай Семенович — канд. физ.-мат. наук, доц.; n.edamenko@spbu.ru

Котина Елена Дмитриевна — д-р физ.-мат. наук, проф.; e.kotina@spbu.ru

Овсянников Александр Дмитриевич — канд. физ.-мат. наук, доц.; a.ovsyannikov@spbu.ru

Optimization of dynamics of trajectory bundles using smooth and nonsmooth functionals. Part 1

D. A. Ovsyannikov1, M. A. Mizintseva1, M. Yu. Balabanov1, A. P. Durkin2, N. S. Edamenko1,

E. D. Kotina1, A. D. Ovsyannikov1

1 St. Petersburg State University, 7-9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation

2 Institute for Nuclear Research of the Russian Academy of Science, 7a, pr. 60-letia Oktiabria, Moscow, 117312, Russian Federation

For citation: Ovsyannikov D. A., Mizintseva M. A., Balabanov M. Yu., Durkin A. P., Edamenko N. S., Kotina E. D., Ovsyannikov A. D. Optimization of dynamics of trajectory bundles using smooth and nonsmooth functionals. Part 1. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2020, vol. 16, iss. 1, pp. 73-84. https://doi.org/10.21638/11702/spbu10.2020.107 (In Russian)

Many different works are devoted to the problems of control and optimization in dynamic systems. Interest in these tasks does not decrease with time. New challenges arise in the development of technological processes in various fields of science and technology, in particular, in the design and creation of modern electrophysical equipment. In this paper, the problem of optimization and control of trajectory beams is considered. The problem

of joint optimization of the program motion and the beam of perturbed motions using a combination of smooth and nonsmooth functionals is investigated. The first part deals with the mathematical formulation of the given representation of the variation of investigated functional and provides optimality conditions in the form of a maximum principle. In the second part, the problem of optimization of dynamics of charged particles in the accelerator with spatially homogeneous quadrupole focusing will be considered. Using a combination of smooth and non-smooth functions allows you to set the functionals that most accurately reflect the requirements for the dynamics of the charged particle beam in accelerators.

Keywords: optimal control, controlled dynamic system, trajectory ensemble, smooth functional, non-smooth functional, maximum principle, charged particle beam, accelerator.

References

1. Kurzhansky A. B. Upravlenie i nabludenie v usloviah neopredelennosti [Control in case of uncertainty ]. Moscow, Fizmatlit Publ., 1977, 394 p. (In Russian)

2. Ovsyannikov D. A. Matematicheskie metody upravlenia puchkami [Mathematical methods of beam control]. Leningrad, Leningrad State University Publ., 1980, 228 p. (In Russian)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

3. Ovsyannikov D. A. Modelirovanie i optimizatsia dinamiki puchkov zarjazhennih chastits [Modeling and optimization of charged particle beam dynamics]. Leningrad, Leningrad State University Publ., 1990, 312 p. (In Russian)

4. Ananyna T. F. Zadacha upravlenia pri nepolnyh dannyh [The problem of control in case of partial data]. Differential Equations, 1976, vol. 13, no. 4, pp. 612—620. (In Russian)

5. Panteleev A., Semenov V. V. Sintez optimalnyh sistem upravlenia pri nepolnoy informatsii [Synthesis of optimal control systems with incomplete information ]. Moscow, Moscow Aviation Institute Publ. House, 1992, 191 p. (In Russian)

6. Bortakovsky A. S. Optimalnoe i suboptimalnoe upravlenia puchkami traektoriy determinirovannyh nepreryvno-diskretnyh sistem [Optimal and suboptimal control of beams of trajectories of deterministic continuous-discrete systems]. Proceedings of the Russian Academy of Sciences. Theory and control systems, 2009, no. 1, pp. 18—33. (In Russian)

7. Kurzhanski A. B., Varaiya P. Dynamic optimization for reachability problems. Journal of Optimization Theory and Applications, 2001, vol. 108, pp. 227—251.

8. Kotina E. D. Discrete optimization problem in beam dynamics. Nucl. Instrum. Methods Phys. Res. A, 2006, vol. 558, no. 1, pp. 292-294.

9. Altsybeyev V. V., Ovsyannikov D. A. Optimization of beam parameters in APF channel. Proceedings of the 27th International Linear Accelerator Conference, LIN AC 2014. Geneva, 2014, pp. 722725.

10. Ovsyannikov D. A. Mathematical methods of optimization of charged partical beams dynamics. Proceedings of the European Partical Accelerator Conference. Barselona, Spain, 1996, vol. 2, pp. 13821384.

11. Bondarev B., Durkin A., Ovsyannikov A. New mathematical optimization models for RFQ structures. Proceedings of the Particle Accelerator Conference. New York City, 1999, vol. 4, pp. 28082810.

12. Vinogradova T. K., Demyanov V. F. O printsipe minimaksa v zadachah optimalnogo upravlenia [On the minimax principle in optimal control problems]. Reports of the USSR Academy of Sciences, 1973, vol. 213, no. 3, pp. 512-514. (In Russian)

13. Demyanov V. F., Vinogradova T. K., Nikulina V. N. et al. Negladkie zadachi teorii optimizatsii i upravlenia [Non-Smooth problems of optimization and control theory]. Leningrad, Leningrad University Publ., 1982, 324 p. (In Russian)

14. Ovsyannikov A. D. Ob odnom klasse zadach optimizatsii v elektrostaticheskom pole [On a class of optimization problems in electrostatic field]. Reports of the Academy of Sciences, 2013, vol. 453, no. 4, pp. 383-385. (In Russian)

15. Ovsyannikov A. D. Upravlenie programmnym i vozmushennym dvizheniamy [Control of the program and perturbed movements]. Vestnik of Saint Petersburg University. Series 10. Applied Mathematics. Computer Science. Control Processes, 2006, iss. 4, pp. 111-124. (In Russian)

16. Ovsyannikov A. D. Upravlenie puchkom zarjazhennih chastist s uchetom ih vzaimodeystvia [Beam control of charged particles taking into account their interaction]. Vestnik of Saint Petersburg University. Series 10. Applied Mathematics. Computer Science. Control Processes, 2009, iss. 2, pp. 82-92. (In Russian)

17. Ovsyannikov A. D. Matematicheskie modeli optimizatsii dinamiki puchkov [Mathematical models of beam dynamics optimization ]. Saint Petersburg, Publishing House "BBM", 2014, 181 p. (In Russian)

18. Ovsyannikov D. A., Altsybeyev V. V. On the beam dynamics optimization problem for an alternating-phase focusing linac. Phys. Part. Nuclei Lett., 2016, vol. 13, no. 7, pp. 775—779.

19. Ovsyannikov D. A., Ovsyannikov A. D., Vorogushin M. F., Svistunov Yu. A., Durkin A. P. Beam dynamics optimization: models, methods and applications. Nucl. Instrum. Methods Phys. Res. A, 2006, vol. 558, no. 1, pp. 11-19.

20. Ovsyannikov A. D., Durkin A. P., Ovsyannikov D. A., Svistunov Y. A. Acceleration of different ion types in single RFQ structure. Problems of Atomic ¡Science and Technology, 2016, vol. 3(103), pp. 5456.

21. Ovsyannikov D. A., Edamenko N. S. Modelirovanie dinamiki puchkov zarjazhennih chastits [Modeling of charged particle beam dynamics]. Vestnik of Saint Petersburg University. Series 10. Applied Mathematics. Computer Science. Control Processes, 2013, iss. 2, pp. 60-65. (In Russian)

22. Ovsyannikov A. D., Ovsyannikov D. A., Altsybeyev V. V., Durkin A. P., Papkovich V. G. Application of optimization techniques for RFQ design. Problems of Atomic Science and Technology, 2014, vol. 3(91), pp. 116-119.

23. Ovsyannikov A. D., Ovsyannikov D. A., Chung S.-L. Optimization of a radial matching section. International Journal of Modern Physics A, 2009, vol. 24(5), pp. 952-958.

24. Ovsyannikov A. D., Ovsyannikov D. A., Chung S.-L. Optimization of matching section of an accelerator with a spatially uniform quadrupole focusing. Technical Physics. The Russian Journal of Appl. Phys., 2009, vol. 11, pp. 1663-1666.

25. Ovsyannikov D. A., Ovsyannikov A. D., Antropov I., Kozynchenko V. BDO-RFQ code and optimization models. Proceedings of the International Conference of Physics and Control, 2005, pp. 282288.

26. Ovsyannikov D. A., Mizintseva M. A., Ovsyannikov A. D. Joint optimization of smooth and nonsmooth functionals on beams of trajectories. Optimal control and differential games. Proceedings of the International Conference dedicated to the 110th anniversary of the birth of Lev Semenovich Pontryagin. Moscow, Steklov Mathematical Institute of RAS Publ., 2018, pp. 203-205.

27. Balabanov M. Yu., Mizintseva M. A., Ovsyannikov D. A. Beam dynamics optimization in a linear accelerator. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2018, vol. 14, iss. 1, pp. 4-13.

28. Mizintseva M. А. Ob odnoi zadache minimizatsii gladkih i neglagkih functsionalov [On one problem of minimization of smooth and non-smooth functionals]. Vestnik of Saint Petersburg State University of Technology and Design. Series 1. Natural and Technical Sciences, 2018, no. 2, pp. 37-41. (In Russian)

29. Mizintseva M., Ovsyannikov D. Minimax problem of simultaneous optimization of smooth and non-smooth functionals. Constructive non-smooth analysis and related, topics. Proceedings of Conference CNSA 2017. IEEE, 2017, pp. 218-222.

30. Mizintseva M., Ovsyannikov D. On the minimax problem of beam dynamics optimization. Proceedings of the 25th Russian Particle Accelerator Conference, RuPAC 2016. Saint Petersburg, 2016, pp. 360-362.

31. Pontryagin L. S., Boltyansky V. G., Gamkrelidze R. V., Mishchenko E. F. Matematicheskaia teoria optimalnyh protsessov [Mathematical theory of optimal processes]. Moscow, Nauka Publ., 1969, 384 p. (In Russian)

32. Kapchinsky I. M. Teoria lineinyh resonansnyh uskoritelei [Theory of linear resonance accelerators]. Moscow, Energoizdat Publ., 1982, 239 p. (In Russian)

33. Golovkina A., Ovsyannikov D., Olaru S. Performance optimization of radioactive waste transmutation in accelerator driven system. Cybernetics and Physics, 2018, vol. 7, no. 4, pp. 210-215.

34. Ovsyannikov D., Zavadskiy S. Pareto-optimal choice of controlled dimension for plasma stabilization system. IFAC-Papers OnLine. Elsevier Science Publishing Company Inc. eISSN: 2405-8963, 2018, vol. 51, no. 32, pp. 175-178.

35. Bazhanov P., Kotina E., Ovsyannikov D., Ploskikh V. Optimization algorithm of the velocity field determining in image processing. Cybernetics and Physics, 2018, vol. 7, no. 4, pp. 174-181.

Received: March 02, 2019.

Accepted: February 13, 2020.

Authors' information:

Dmitri A. Ovsyannikov — Dr. Sci. in Physics and Mathematics, Professor; d.a.ovsyannikov@spbu.ru

Maria A. Mizintseva — Magister; m.mizintseva@spbu.ru

Mihail Yu. Balabanov — PhD in Physics and Mathematics, Senior Researcher Worker; m.balabanov@spbu.ru

Alexander P. Durkin — PhD in Technics, Associate Professor; durkinap@mail.ru

Nikolai S. Edamenko — PhD in Physics and Mathematics, Associate Professor; n.edamenko@spbu.ru

Elena D. Kotina — Dr. Sci. in Physics and Mathematics, Professor; e.kotina@spbu.ru

Alexander D. Ovsyannikov — PhD in Physics and Mathematics, Associate Professor; a.ovsyannikov@spbu.ru

i Надоели баннеры? Вы всегда можете отключить рекламу.