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

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

11-2012-72 На правах рукописи УДК 519.6:621.384.633

КИЯН Игорь Николаевич

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РЕЖИМОВ РАБОТЫ МНОГОЦЕЛЕВОГО ИЗОХРОННОГО ЦИКЛОТРОНА АИЦ-144

Специальность: 05.13.18 — математическое моделирование, численные методы и комплексы программ

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

О 4 0НТ2012

Дубна 2012

005052466

Работа выполнена в Лаборатории ядерных проблем Объединенного института ядерных исследований.

Научные руководители: кандидат физико-математических наук

Амирханов Илькизар Валиевич, Лаборатория информационных технологий ОИЯИ

кандидат физико-математических наук Карамышева Галина Анатольевна, Лаборатория ядерных проблем ОИЯИ

доктор физико-математических наук, профессор Севастьянов Леонид Антонович, Российский университет дружбы народов, профессор

кандидат физико-математических наук, с.н.с. Казаринов Николай Юрьевич, Лаборатория ядерных реакций ОИЯИ, начальник сектора

Ведущая организация: НИИ электрофизической аппаратуры им. Д.В.Ефремова, г. Санкт-Петербург

Защита состоится « О » Ю/^Р^А-иАЛ. 2012 г. в ^-У^часов на заседании диссертационного совета ^Д 720.001.04 в Лаборатории информационных технологий Объединенного института ядерных исследований, г. Дубна Московской области.

С диссертацией можно ознакомиться в библиотеке ОИЯИ.

Автореферат разослан «З 0» &-2012 г.

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

Официальные оппоненты:

профессор

^ - Иванченко Иосиф Моисеевич

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

В диссертации изложены результаты разработок и исследований, выполненных соискателем в Лаборатории ядерных проблем Объединенного института ядерных исследований (ЛЯП ОИЯИ) в период с 2004 по 2012 год в рамках проблемно-тематических планов научно-исследовательских работ ОИЯИ (№ 03-2-1102-2010/2012 "Совершенствование Фазотрона ЛЯП (ОИЯИ) и разработка циклотронов для физических и прикладных исследований") и протоколов о международном сотрудничестве, осуществляющемся между ОИЯИ и Институтом ядерной физики Польской академии наук (ИЯФ ПАН) (№ 3970-210/12 "Протокол о выполнении совместной научно-исследовательской работы"). Работа выполнена при финансовой поддержке РФФИ, грант № 10-01-00467-а.

Актуальность. В настоящее время одним из главных направлений в области развития ускорительной техники является разработка, производство и эксплуатация универсальных изохронных циклотронов. Они находят применение как в области научных исследований (производство радиоизотопов для исследований в области ядерной спектроскопии, электроядерная энергетика), так и в области медицины (производство радиоизотопов для различных видов томографии, протонная терапия раковых опухолей): DC-601, AIC-1442, VINCY3. Циклотроны данного вида являются многоцелевыми и могут использоваться для ускорения частиц с различным отношением числа стандартных единиц атомной массы к числу элементарных зарядов (A/Z). Энергия вывода в таких машинах может плавно варьироваться для каждого типа ускоряемых частиц.

Диссертационная работа посвящена созданию и применению на практике новой методики математического моделирования режимов работы многоцелевого изохронного циклотрона АИЦ-144, расположенного в ИЯФ ПАН и

1 http://fls2.iinr.rii/flnr/dc60.html

2 http://www.ifj.edu.pl/str/dc/index.htnil

3 http://www.vinca-at.org/tesla.html

предназначенного как для проведения протонной терапии меланомы глаза у пациентов (при работе в основном режиме), так и для производства радиоизотопов, используемых в научных исследованиях в области ядерной спектроскопии. АИЦ-144 был разработан группой сотрудников ИЯФ ПАН под руководством д-ра Е.Схвабэ4. Система вывода АИЦ-144 была разработана группой сотрудников ОИЯИ при участии сотрудников ИЯФ ПАН5. Рассматриваемая методика позволяет моделировать новые режимы работы АИЦ-144 и переходить на работу с ними без остановки и разборки циклотрона для проведения дополнительных магнитных измерений. Новая методика содержит в себе набор алгоритмически связанных методов для расчета токов в главной катушке, двадцати концентрических катушках, в двух парах гармонических катушек, а также для численного контроля точности расчета изохронного магнитного поля, рассчитываемого по аналитическим формулам М.М.Гордона6 (Прил., п.1). Общий алгоритм моделирования новых режимов работы АИЦ-144 включает в себя расчет токов во всех катушках циклотрона [К6МК13], численный контроль точности расчета изохронного магнитного поля [К14], аналитический расчет фазового движения7 и частот свободных колебаний [КЗ] ускоряемых ионов. В общий алгоритм входят рабочие алгоритмы, позволяющие по наибольшему значению тока выведенного пучка проводить: • Улучшение фазового движения ускоряемых частиц, которое проводится за счет регулировки амплитуды центрального бампа результирующего магнитного поля циклотрона (от токов в главной катушке и в двадцати концентрических катушках) [К 13].

4 Е.Схвабэ и др. Магнитная система автоматического изохронного циклотрона AIC-144. // Proceedings of the International Seminar on Isochronous Cyclotron Technique. Institute of Nuclear Physics, Krakow, 13th—18th November, 1978. Pp. 237-256. Krakow, 1978.

5 Nikolaj A. Morozov and others. A new extraction system for the upgraded AIC-144 cyclotron. // NUKLEONIKA, 2001. International Journal of Nuclear Research, V.46 (2). Pp. 51-57.

6 M.M.Gordon. Calculation of Isochronous Fields for Sector-Focused Cyclotrons // Particle Accelerators, 1983. Vol. 13. Pp. 67-84.

7 T.Stammbach. Cyclotrons. // Joint Universities Accelerator School. Archamps-France. 1997.

• Оптимизацию параметров магнитной системы циклотрона (регулировку амплитуды и фазы первой гармоники рабочего магнитного поля ускорителя (от токов в главной катушке, в двадцати концентрических катушках и в двух парах гармонических катушек), проводимую на радиусе заброса пучка ускоренных ионов в электростатический дефлектор) [К13].

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

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

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

дополнительных магнитных измерений. Моделирование в этом случае осуществляется на основе уже имеющихся измеренных и пересчитанных карт исходных магнитных полей. Такой подход требует разработки и развития, как ряда новых математических моделей, так и ряда программных и рабочих алгоритмов. При использовании нового подхода, время перехода на моделируемый режим работы значительно сокращается. Например, переход на основной режим работы АИЦ-144, осуществленный в июне 2009 года, занял всего шесть часов. Сюда вошло время на аналитические и численные расчеты, измерение тока внутреннего и внешнего пучка ускоренных ионов, а также время на улучшение фазового движения ускоряемых частиц. Еще два часа ушло на оптимизацию параметров магнитной системы циклотрона, проведенную на ускорителе в апреле 2011 года.

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

Цель диссертационной работы заключается:

1: В разработке новой методики математического моделирования режимов работы многоцелевого изохронного циклотрона АИЦ-144, позволяющей переходить на работу в моделируемых режимах без остановки и разборки ускорителя для проведения дополнительных магнитных измерений.

2. В создании на языке программирования С++ в среде MS Visual Studio .NET комплекса программ Cyclotron Operator HELP Program Complex 2004-2012 и динамически связанной библиотеки GaussDLL, реализующих разработанные математические модели и алгоритмы.

3. В практическом применении разработанного программного обеспечения при компьютерном моделировании новых режимов работы АИЦ-144 (В том числе основного режима работы циклотрона: протоны, кинетическая энергия вывода Ек -60/60,5 МэВ, частота ВЧ-генератора Frf =26,155/26,25 МГц).

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

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

1. Разработана новая методика математического моделирования режимов работы многоцелевого изохронного циклотрона АИЦ-144, позволяющая переходить на работу в моделируемых режимах без остановки и разборки ускорителя для проведения дополнительных магнитных измерений [К1]+[К14].

2. Впервые созданы математическая модель и алгоритм, которые позволяют проводить расчет тока в главной катушке АИЦ-144 на основе карт исходных магнитных полей, измеренных для трех уровней тока в катушке возбуждения главного магнита [К9].

3. Разработан новый метод расчета токов в двадцати концентрических катушках для рассчитанного уровня тока в главной катушке АИЦ-144, основанный на минимизации квадратичного функционала, в состав которого введена штрафная функция. Данная функция предназначена для направленного покомпонентного

ввода вектора решения (вектора нормированных значений оптимизированных токов в концентрических катушках) в границы области допустимых значений. Ввод найденного вектора решения в рамки заданных граничных значений осуществляется в направлении, определяемом в соответствии с расположением зон влияния концентрических катушек на основное магнитное поле АИЦ-144 (от центра к краю циклотрона) [К6]-г[К10].

4. Впервые созданы математическая модель и алгоритм, которые позволяют проводить расчет токов в двух парах гармонических катушек для рассчитанного уровня тока в главной катушке АИЦ-144. Расчет проводится путем решения неоднородной системы нелинейных алгебраических уравнений, решаемой итерационным способом [К11]+[К13].

5. Усовершенствован метод контроля точности расчета изохронного магнитного поля, рассчитываемого по аналитическим формулам М.М.Гордона. Показано, что для основного режима работы АИЦ-144 значение модуля относительной ошибки частоты обращения протонов во всем диапазоне изохронизации формируемого магнитного поля составляет | crF0 \<2,5Е-4, что сопоставимо по своей величине со значением модуля относительной погрешности измерения карт исходных магнитных полей, измеряемых с помощью датчика Холла, | ams \~1+2Е-4 [К14]. Практическая значимость

1. На языке программирования С++ в среде MS Visual Studio .NET в стандартах DBA (Dialog-based Application) и SDI (Single Document Interface) создано программное обеспечение, реализующее разработанные математические модели и алгоритмы:

• Комплекс программ Cyclotron Operator HELP Program Complex 2004-2012, установленный на сервере АИЦ-144 и используемый для моделирования новых режимов работы циклотрона.

• Динамически связанная библиотека GaussDLL, включающая в свой состав комплекс программ матрично-векторных операций и решения СЛАУ с

вещественными коэффициентами, размещенная на сайте ОИЯИ в библиотеке готовых программ JINRLIB [К 15].

2. С помощью комплекса программ Cyclotron Operator HELP Program Complex 2009 для АИЦ-144 было выполнено компьютерное моделирование основного режима его работы, предназначенного для ускорения пучка протонов с необходимыми для протонной терапии меланомы глаза параметрами: кинетической энергией вывода Ек~60 МэВ, значением и стабильностью тока выведенного пучка IЬсат ш >20 нА, <xlheamtxl ~±5 % соответственно [К7]+[К10].

3. Смоделированный с помощью созданного программного обеспечения основной режим работы АИЦ-144 был установлен на циклотроне в июне 2009 года без эмпирической подстройки токов в концентрических катушках и частоты ВЧ -генератора. Благодаря эффективной работе сотрудников Циклотронного отдела ИЯФ ПАН и высокой стабильности установленного режима работы, подразделение протонной терапии ИЯФ ПАН впервые на протяжении ряда лет получило возможность бесперебойной работы и использования пучка выведенных протонов:

• Для настройки медицинских приборов и аппаратуры, которая проводилась в течение около двух лет с момента установки смоделированного режима работы на АИЦ-144.

• Для проведения успешной протонной терапии меланомы глаза у первой партии пациентов, впервые проведенной в Польше на АИЦ-144 в период с февраля по апрель 2011 года [К12]-н[К14].

4. С помощью комплекса программ Cyclotron Operator HELP Program Complex 2011 для основного режима работы АИЦ-144 в апреле 2011 года была проведена оптимизация параметров магнитной системы циклотрона [К 13]. Выведенный пучок ускоренных протонов был успешно использован в ИЯФ ПАН для производства радиоизотопов, нарабатываемых с целью проведения научных исследований в области ядерной спектроскопии.

Практическая значимость разработанных методов и комплексов программ документально подтверждена: отчетами о проведенных на АИЦ-144 в период с 2004 по 2011 год физических экспериментах, подписанными начальником ЦО ИЯФ ПАН (Прил. п.2); письмом от дирекции ИЯФ ПАН к дирекции ЛИТ ОИЯИ BD/47/2011 от 27 апреля 2011 года об успешных результатах лечения раковой опухоли глаза у первых девяти пациентов (Прил. п.З). Новая методика и разработанное программное обеспечение, используемые в настоящее время для моделирования режимов работы АИЦ-144, могут найти широкое практическое применение в процессе эксплуатации аналогичных ускорителей, расположенных в различных российских и зарубежных научно-экспериментальных центрах. Апробация работы. Основные положения и результаты диссертации представлены и доложены на следующих российских и международных конференциях:

• VI Международный семинар по проблемам ускорителей заряженных частиц памяти В.П.Саранцева, 8-10 сентября, 2005. Алушта, Украина.

• XXI Russian Particle Accelerator Conference 2008 (RuPAC 2008), September 28 -October 3,2008. Zvenigorod, Russia.

• The International Conference Mathematical Modeling and Computational Physics 2009 (MMCP 2009), July 7-11,2009. Dubna, Russia.

• Всероссийская конференция (с международным участием): Информационно-телекоммуникационные технологии и математическое моделирование высокотехнологичных систем. 18-22 апреля, 2011. Москва, Россия.

• IX Международный семинар по проблемам ускорителей заряженных частиц памяти В.П.Саранцева 17-21 сентября, 2011. Алушта, Украина.

Результаты, включенные в диссертацию, докладывались на семинарах ИЯФ ПАН, Лаборатории информационных технологий и Лаборатории ядерных проблем Объединенного института ядерных исследований (ЛИТ ОИЯИ, ЛЯП ОИЯИ).

Публикации. Материалы диссертации опубликованы в 14 печатных работах, из которых 7 статей в рецензируемых научных журналах, рекомендуемых ВАК (в одном зарубежном и трех российских) [К2], [К5], [Кб], [К8], [К9], [К13], [К14], 3 статьи в научных сборниках [К10]-г[К12], 2 статьи в сообщениях ОИЯИ [КЗ], [К4] и 2 аннотации докладов в материалах конференций [Kl], [К7], а также в виде 1 распространяемого программного продукта [К 15].

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

Достоверность представленных в диссертации результатов моделирования новых режимов работы АИЦ-144 подтверждается:

• Результатами тестирования разработанного программного обеспечения: комплекса программ Cyclotron Operator HELP Program Complex 2004-2012 и динамически связанной библиотеки GaussDLL.

• Сопоставлением результатов компьютерного моделирования основного режима работы АИЦ-144, выполненного с помощью различных версий разработанного программного обеспечения.

• Методическими результатами контроля точности расчета изохронного магнитного поля, полученными путем сопоставления аналитических и численных

оценок, выполненных на основе использования расчетной карты результирующего магнитного поля АИЦ-144.

• Успешными результатами научно-методических натурных экспериментов, проведенных на АИЦ-144 в период с 2004 по 2011 год. (Отчеты, удостоверенные подписью начальника ЦО ИЯФ ПАН (Прил. п.2)).

• Успешными результатами проведения протонной терапии меланомы глаза у первой партии пациентов, впервые проведенной в Польше на АИЦ-144 в период с февраля по апрель 2011 года. (Письмо от дирекции ИЯФ ПАН к дирекции ЛИТ ОИЯИ ВО/47/2011 от 27 апреля 2011 года (Прил. п.З)).

Структура и объем диссертации. Диссертация состоит из введения, трех глав, заключения и приложения. Список работ, опубликованных по теме диссертации, содержит 15 наименований. Полный объем диссертации составляет ^1$5страниц машинописного текста, включая таблиц и ^У^рисунков.

СОДЕРЖАНИЕ РАБОТЫ

Во введении приводятся общая характеристика и краткое содержание диссертационной работы, ставится задача на разработку новой методики математического моделирования режимов работы многоцелевого изохронного циклотрона АИЦ-144, позволяющей переходить на работу в моделируемых режимах без остановки и разборки ускорителя для проведения дополнительных магнитных измерений, дается ретроспектива решения поставленной задачи, приводится описание основных параметров и магнитной системы АИЦ-144. Первая глава диссертации посвящена постановке задачи по моделированию основного режима работы АИЦ-144, предназначенного для ускорения пучка протонов с необходимыми для протонной терапии меланомы глаза параметрами. В этой главе рассматривается порядок проведения магнитных измерений, технология подготовки исходных данных, приводится описание разработанных математических моделей для расчета токов во всех катушках АИЦ-144 [К6]-г[К13]. В первой главе рассматриваются уравнения движения, записанные в

приближении замкнутых равновесных орбит (без учета напряженности электрического поля), численное решение которых позволяет проконтролировать точность расчета изохронного магнитного поля в области изохронизации формируемого магнитного поля циклотрона [К 14]. Все расчеты проводятся в безразмерных величинах (безразмерность параметра обозначается *). Под средними магнитными полями понимаются поля, усредненные по азимуту. В первом параграфе ставится задача на моделирование основного режима работы АИЦ-144: протоны, кинетическая энергия вывода Ек-60/60,5 МэВ, частота ВЧ-генератора F^ =26,155/26,25 МГц. В этом параграфе приводится содержание

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

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

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

Математическая модель для расчета тока в главной катушке АИЦ-144 1пшт шпк

строится с помощью интерполяционного полинома Лагранжа второй степени. Расчет проводится исходя из заданного значения радиуса рабочей точки (места пересечения среднего основного и изохронного магнитных полей в области конечных радиусов ускорения ионов) Rwork. Значение радиуса рабочей точки К,,гк определяет значение рабочего тока в главной катушке циклотрона Imam л.

Указанная математическая модель и расчетный алгоритм подробно рассмотрены в статье [К9].

Математическая модель для расчета токов в двадцати концентрических катушках АИЦ-144 /,, у =1,2..л (2=20) для рассчитанного уровня тока в главной катушке I I т,гк строится на основе использования метода наименьших

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

2

ер!ит к ( \ _т / \ ^ /

о

1

X =

7=1

J=1

\В0 ■ -^Lq

, _ Ij,max г' _ 11 / > 0 (2)

Kj) ш~ f ' 'j - т > j,max W

* i lim i ma*

V.lim У,шах

где Бу,тах(г*) - это среднее магнитное поле максимального вклада j-ой концентрической катушки в основное магнитное поле; Ir IjM, /,тах-расчетный, предельный и максимально допустимый токи в j-ой концентрической катушке; г - средний радиус замкнутой равновесной орбиты; R'seplum - радиус

расположения первого ножа электростатического дефлектора; X - переменный параметр штрафной функции; р - постоянный параметр штрафной функции; к -число свободных компонент вектора решения.

(3) 5L(/)=sL(r')+BLt(r') (4)

8 C.M.Davisson, R.B.Theus. Linear Programming Vs Least Square Method in Determing the Best Trim Coil Currents for the NRL cyclotron. // Particle Accelerators, 1976. Vol. 7. Pp. 103110.

Втшк{г ) = бВьитр (г )+ <5 ВеЛце (г ) (5) В иап\Г )— Втат\Г )+ У, В),сопи (О (6)

7=1

%,соп„{г')=В].т (7)

У,шах

где Д5 (г*) - это разность между средним требуемым В,тк(г) и средним исходным магнитным полем Вмп(г')-, В,¡{г') - изохронное магнитное поле;

Втшк(г') - маска изохронного магнитного поля; бВьшР(г') - поправка к изохронному магнитному полю, формирующая среднее магнитное поле в центре циклотрона; бВе^е(г') - поправка к изохронному магнитному полю, формирующая среднее краевое магнитное поле; Втып{г ) - среднее основное

магнитное поле; - среднее магнитное поле вклада }-ой

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

ой концентрической катушке; т - это число фиксированных компонент вектора решения. Для }-ой концентрической катушки нижнее и верхнее

граничные значения ее тока задаются в рамках соответствующих минимально и максимально допустимых значений. То есть, предельный ток в у'-ой концентрической катушке 1)Мт определяется следующим образом:

^ = / (0 < / <1 < I ) ' ^ > 0'7 =/'2-2' 2 = к + т (8)

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

нелинейных алгебраических уравнений, которая решается итерационным путем. На каждом шаге итерационного цикла формируется неоднородная система линейных алгебраических уравнений, которая решается с помощью метода Гаусса с выбором главного элемента по матрице коэффициентов. Устойчивость решения обеспечивается за счет исключения из расчета части задействованных концентрических катушек. Исключение проводится на основе использования нового критерия (минимального значения произведения функционального минимума на число обусловленности СЛАУ). Процесс построения указанной математической модели, описание расчетного алгоритма, результаты численных и физических экспериментов, проведенных на АИЦ-144, подробно рассмотрены в статьях [К6МК10].

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

Математическая модель для расчета токов в двух парах гармонических катушек АИЦ-144 1к, к=г+1,г+2 (г=20) для рассчитанного уровня тока в главной катушке 1„шЫ у10гк строится на основе использования гармонического

анализа рабочего магнитного поля циклотрона. Рабочее магнитное поле АИЦ-144 записывается в следующем виде:

К* (>■'■ <Р)= В1{г, <р)+ Ъв'Ыт (г*, <р)- ¡1), ¡1 = -А- (9)

*=--+! 'к, тах

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

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

1^*,тах — Д./о» < Ь < 0)1 А, = •< Г, и >0, к =г+1,г+2 (10)

*'1ш1 I/ (0< I <1 </ •П'-1*.™* ' 4 '

В результате ряда преобразований получается неоднородная система нелинейных алгебраических уравнений. При фиксированных значениях амплитуды и фазы первой гармоники рабочего магнитного поля, заданных на радиусе заброса пучка ускоренных ионов в электростатический дефлектор, неоднородная система нелинейных алгебраических уравнений превращается в неоднородную систему линейных алгебраических уравнений, которая решается стандартным путем. Процесс построения указанной математической модели, описание рабочего алгоритма оптимизации параметров магнитной системы АИЦ-144, результаты численных и физических экспериментов, проведенных на циклотроне, подробно рассмотрены в статьях [К11]+[К13].

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

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

В качестве исходного уравнения для расчета динамики пучка ускоряемых ионов берется уравнение Лоренца-Ньютона, описывающее движение одиночной заряженной частицы в электромагнитном поле циклотрона. ]5 = дЁ + д-[?хв] (11)

где р - это релятивистский импульс заряженной частицы; - заряд частицы; г - радиус-вектор ускоряемого иона; Е - вектор напряженности электрического поля; В - вектор индукции магнитного поля. В развернутом представлении уравнение Лоренца-Ньютона выглядит следующим образом:

Ё + [?хв]--^-(?-Ё)\ (12)

mr = q-

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

индукции В на ось OZ имеет отрицательные значения. В результате данного разложения получаются уравнения движения (система нелинейных дифференциальных уравнений второго порядка). Уравнения движения, записанные в приближении замкнутых равновесных орбит, решаются методом Рунге-Кутта. Математическая модель, выведенные аналитические выражения, использующиеся для расчета начальных условий, а также результаты аналитических и численных расчетов, выполненных для основного режима работы АИЦ-144, подробно приведены в статье [К14].

В шестом параграфе приводятся основные выводы, сделанные на основе анализа материалов, изложенных в первой главе.

Вторая глава диссертации посвящена описанию созданного языке программирования С++ в среде MS Visual Studio .NET в стандартах DBA (Dialog-

based Application) и SDI (Single Document Interface) программного обеспечения, реализующего разработанные математические модели и алгоритмы. В первом параграфе приводится описание комплекса программ Cyclotron Operator HELP Program Complex 2004-2012, предназначенного для компьютерного моделирования новых режимов работы АИЦ-144. Комплекс состоит из трех программ:

1) Magnetic Field Calculation Program 2012 - программа, предназначенная для расчета карт дополнительных магнитных полей.

2) Cyclotron Mode Calculation Program 2012 - программа, предназначенная для расчета токов во всех катушках циклотрона, а также для численного расчета частоты обращения и фазового движения ускоряемых частиц.

3) Cyclotron Analytic Model Program 2012 - программа для аналитического расчета изохронного магнитного поля, частот свободных колебаний, а также частоты обращения и фазового движения ускоряемых частиц.

Во втором параграфе приводится описание динамически связанной библиотеки GaussDLL, включающей в свой состав комплекс программ матрично-векторных операций и решения СЛАУ с вещественньми коэффициентами [К 15]. В третьем параграфе приводятся основные выводы, сделанные на основе анализа материалов, изложенных во второй главе.

Третья глава диссертации посвящена описанию технологии компьютерного моделирования новых режимов работы АИЦ-144, моделирование которых проводится с помощью комплекса программ Cyclotron Operator HELP Program Complex 2004—2012. В этой главе приводятся результаты физических экспериментов и эксплуатации АИЦ-144 в основном режиме его работы. В первом параграфе приводятся результаты компьютерного моделирования основного режима работы АИЦ-144: протоны, кинетическая энергия вывода Ек ~60/60,5 МэВ, частота ВЧ-генератора Frf =26,155/26,25 МГц. Во втором параграфе приводятся результаты запуска и эксплуатации АИЦ-144 в основном режиме его работы [К7]-г[К14].

В период с 1998 по 2011 год на АИЦ-144 был проведен ряд магнитных измерений. Карты основных магнитных полей были измерены для трех уровней тока в главной катушке 1тат ] =535/570 А; 1та1„ 2 =570/585 А; 1тат ¡=600/600 А в

1998+1999, 2006/2011 годах соответственно. Карты дополнительных магнитных полей были измерены при токе в главной катушке 1тат =570 А для концентрических катушек №1+20 и для двух пар гармонических катушек в 1998+1999 годах; для концентрических катушек №17+20 в 2011 году. Карты основных магнитных полей NOWE010.TXT, 2006R045.TXT и 2011R109.TXT, 2011R.152.TXT были измерены при токе в главной катушке 1таШ =570 А в 1999, 2006 годах и в 2011 году до и после установки вакуумной камеры в правильное положение соответственно. Содержащиеся в них основные магнитные поля значительно отличаются друг от друга по своей форме и величине амплитуды первой гармоники В, (г). На Рис. 1 показан уход среднего основного магнитного поля, оцененный на основе измерений, сделанных в 2011 году, по сравнению с измерениями, сделанными в 2006 и 1999 годах. Значения средних основных магнитных полей для 2006R045.TXT и 2011R152.TXT на радиусе начала краевого магнитного поля Ие^е=0,б м приведены к единой величине, полученной для NOWE010.TXT на том же радиусе. На Рис. 2 представлена амплитуда первой Гармоники Вх (г) основного магнитного поля, измеренного при токе в главной катушке 1тШ =570 А в 2006 году и в 2011 году до и после установки вакуумной камеры в правильное положение. Коррекция положения вакуумной камеры [К14] понадобилась из-за ее самопроизвольного сдвига, определенного из расчетов, сделанных во время оптимизации параметров магнитной системы АИЦ-144, проведенной в апреле 2011 года [К13]. Для обеспечения возможности оптимизации амплитуды и фазы первой гармоники рабочего магнитного поля амплитуда первой гармоники основного магнитного поля должна быть В{ =3±1 Гс на радиусе заброса пучка ускоренных ионов в электростатический дефлектор Ка,=0,б2м.

-1 Bmean(2011R152) - Bmean(2006R045) -2 Bmean(201TR152) -;Bmean(ftowe010)

--1 B1(2006R045)

--2B1{2011R109)

--ЗВ1Й011В15?!

R(m)

Рис. 1. Разность средних основных магнитных полей.

Рис. 2. Амплитуда первой гармоники основного магнитного поля.

Из Рис. 1 видно, что наибольший уход среднего основного магнитного поля наблюдается в краевой области. С 1999 по 2011 год он составил \&Втеа„ | «55 Гс на

радиусе R =0,7 м. Из Рис. 2 видно, что после установки вакуумной камеры в правильное положение в сентябре 2011 года величина амплитуды первой гармоники на радиусе заброса пучка ускоренных ионов в электростатический дефлектор Rexl =0,62 м была уменьшена от 5, -»73 Гс до В, я4 Гс. На Рис. 3-гРис. 5 показаны токи внутренних пучков ускоряемых протонов в зависимости от радиуса, измеренные на АИЦ-144 для основного режима его работы в 2007, 2009 и 2011 годах (Прил., п.2).

-Ibeam internal Cyclotron Operator HELP 2007 -Ibeam internal Cyclotron Operator HELP 2009

Рис. 3. Измеренный ток внутреннего пучка.

Рис. 4. Измеренный ток внутреннего пучка.

На Рис. 6 показан ток пучка выведенных протонов в зависимости от тока в главной катушке, измеренный на АИЦ-144 для основного режима его работы в 2011 году (Прил., п.2). Значение тока пучка выведенных протонов 1Ьеат =2 нА

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

- Ibeam internal Cyclotron Operator HELP 2011

Experiment 22.04.11;

.............Frfi 26.155 MHz!..........................'.........................!........................

.....I.............Udee = 60 kV ............J.................................................!.....................

Rdeflector = 62,8 cm I

Rbump = 16 cm -i- Redge = 59,5 cm j Imain = 577,23 A

егаЬгНЕЬР 2011

Experiment 22.04.11 Frf = 26,155 MHz

60 kV Rdellector = 62,8 on = 16 cm + s ■ 59,5 cm

R(m)

576.4 576.8 577.2 577.6 578.0 578.4 578.8 579.2 Imain coil (A)

Рис. 5. Измеренный ток внутреннего пучка.

Рис. 6. Измеренный ток выведенного пучка.

Из Рис. 3 видно, что в 2007 году внутренний пучок был ускорен без каких-либо потерь протонов в области изохронизации сформированного магнитного поля (от радиуса конца центрального бампа ЯЬитр до радиуса начала краевого магнитного поля Яы,:1,). Из рисунков Рис. 4, Рис. 5 видно, что в указанной области рабочих

радиусов к 2009 году появились потери ускоряемых протонов, которые несколько возросли к 2011 году. Ток внутреннего пучка, измеренный в эксперименте от 22.04.2011 года, по сравнению с измеренным в эксперименте от 04.06.2009 года, начинает спадать в краевом магнитном поле несколько раньше. Оптимальный ток в главной катушке АИЦ-144, установленный в эксперименте от 22.04.2011 года, по сравнению с установленным в эксперименте от 04.06.2009 года, несколько больше по своей величине. Появление и рост потерь ускоряемых протонов объясняются как уходом основного магнитного поля в краевой области (Рис. 1), так и увеличением амплитуды первой гармоники во всем диапазоне рабочих радиусов ускорения частиц от источника ионов до системы вывода (Рис. 2). Из

Рис. 6 видно, что при изменении тока пучка выведенных протонов от 90 % до 90 % от его установленного значения 1Ьат^я=2 нА, диапазон изменения тока в

главной катушке составил Л1татсо:,=0,44 А, что говорит о хорошем качестве изохронизации сформированного магнитного поля и гарантирует стабильную работу АИЦ-144 в установленном режиме.

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

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

1. Разработана новая методика математического моделирования режимов работы многоцелевого изохронного циклотрона АИЦ-144, позволяющая переходить на работу в моделируемых режимах без остановки и разборки ускорителя для проведения дополнительных магнитных измерений [К1]-=-[К14].

2. Впервые созданы математическая модель и алгоритм, которые позволяют проводить расчет тока в главной катушке АИЦ-144 на основе карт исходных магнитных полей, измеренных для трех уровней тока в катушке возбуждения главного магнита [К9].

3. Разработан новый метод расчета токов в двадцати концентрических катушках для рассчитанного уровня тока в главной катушке АИЦ-144, основанный на минимизации квадратичного функционала, в состав которого введена штрафная функция, предназначенная для направленного покомпонентного ввода вектора решения в границы области допустимых значений [К6]+[К10].

4. Впервые созданы математическая модель и алгоритм, которые позволяют проводить расчет токов в двух парах гармонических катушек для рассчитанного уровня тока в главной катушке АИЦ-144. Расчет проводится путем решения неоднородной системы нелинейных алгебраических уравнений, решаемой итерационным способом [К11]+[К13].

5. Усовершенствован метод контроля точности расчета изохронного магнитного поля, рассчитываемого по аналитическим формулам М.М.Гордона. Показано, что для основного режима работы АИЦ-144 значение модуля относительной ошибки

частоты обращения протонов во всем диапазоне изохронизации формируемого магнитного поля составляет |.оу0 \<2,5Е-4, что сопоставимо по своей величине со значением модуля относительной погрешности измерения карт исходных магнитных полей, измеряемых с помощью датчика Холла, | \~l-r2E-4 [К14].

6. На языке программирования С++ в среде MS Visual Studio .NET в стандартах DBA (Dialog-based Application) и SDI (Single Document Interface) создан комплекс программ Cyclotron Operator HELP Program Complex 2004-2012, реализующий разработанные математические модели и алгоритмы.

7. С помощью комплекса программ Cyclotron Operator HELP Program Complex 2009 для АИЦ-144 было выполнено компьютерное моделирование основного режима его работы, предназначенного для ускорения пучка протонов с необходимыми для протонной терапии меланомы глаза параметрами [К7]-ь[К10].

8. Смоделированный с помощью созданного программного обеспечения основной режим работы АИЦ-144 был установлен на циклотроне в июне 2009 года без эмпирической подстройки токов в концентрических катушках и частоты ВЧ-генератора. Установленный режим работы был использован при проведении успешной протонной терапии меланомы глаза у первой партии пациентов, впервые проведенной в Польше на АИЦ-144 в период с февраля по апрель 2011 года [К12МК14].

9. С помощью комплекса программ Cyclotron Operator HELP Program Complex 2011 для основного режима работы АИЦ-144 в апреле 2011 года была проведена оптимизация параметров магнитной системы циклотрона [К13]. Выведенный пучок ускоренных протонов был успешно использован в ИЯФ ПАН для производства радиоизотопов, нарабатываемых с целью проведения научных исследований в области ядерной спектроскопии.

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

Список работ, опубликованных по теме диссертации

К1. И.Н.Киян. Новый метод моделирования режимов работы многоцелевого изохронного циклотрона. // 6-ой Международный семинар памяти профессора В.П.Саранцева, 2005. Алушта, Украина. Сборник аннотаций D9-2005-123. Дубна 2005. С. 16.

К2. И.Н.Киян. Новый метод моделирования режимов работы многоцелевого изохронного циклотрона. // ПИСЬМА В ЭЧАЯ, 2006. Т.З, N7(136). Стр. 49-54. КЗ. И.Н.Киян, Р.Тарашкевич. Анализ области применения аналитических формул для расчета частот свободных колебаний в изохронных циклотронах. // Сообщение ОИЯИ Р9-2005-114, Дубна 2005.

К4. И.Н.Киян, Р.Тарашкевич. Методы расчета режимов работы и параметров изохронного циклотрона. // Сообщение ОИЯИ Р9-2006-40, Дубна 2006. К5. Igor N. Kiyan and others. Simulation of operation modes of isochronous cyclotron by a new iterative method. // NUKLEONIKA, 2007. International Journal of Nuclear Research, V.52 (1). Pp. 29-34.

Кб. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян, Я.Суликовский. Моделирование требуемых режимов работы и анализ их устойчивости для многоцелевых изохронных циклотронов. // ПИСЬМА В ЭЧАЯ, 2009. Т.6, № 6(155). Стр. 805-813.

К7.1.V.Amirkhanov, G.A.Karamysheva, I.N.Kiyan, J.Sulikowski. Mathematical Aspects of Modeling of Required Operation Modes of Multi Purpose Isochronous Cyclotrons. // The International Conference: Mathematical Modeling and Computational Physics, 2009. Dubna, Russia. Book of Abstracts D11-2009-88. Dubna 2009. Pp. 203-204.

K8.1.V.Amirkhanov, G.A.Karamysheva, I.N.Kiyan, J.Sulikowski. Mathematical Aspects of Modeling of Required Operation Modes of Multi Purpose Isochronous Cyclotrons. // BULLETIN of PFUR. Series Mathematics. Information Sciences. Physics, 2010. No 2 (2). Pp. 99-103.

К9. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян, Я.Суликовский. Методика математического моделирования режимов работы многоцелевых изохронных циклотронов. // ВЕСТНИК ТГУ. Серия: Прикладная математика, 2010. №9, Выпуск 1 (14). Стр. 55-68.

К10. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян. Я.Суликовский. Улучшенная математическая модель расчета токов в концентрических катушках коррекции основного магнитного поля многоцелевого изохронного циклотрона. // Scientific report 2008-2009, LIT JINR. 2009-196. Pp. 109-112.

Kll. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян, Я.Суликовский. Решение обратной задачи расчета токов в гармонических катушках многоцелевого изохронного циклотрона. // Информационно-телекоммуникационные технологии и математическое моделирование высокотехнологичных систем. Тезисы докладов Всероссийской конференции с международным участием. 18-22 апреля 2011 года. Российский университет дружбы народов. Москва. 2011. Стр. 253-255. К12. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян. Я.Суликовский. Математическая модель расчета токов в двух парах гармонических катушек многоцелевого изохронного циклотрона. Контроль точности расчета изохронного магнитного поля. // Scientific report 2010-2011, LIT JINR. 2011-130. Pp. 135-138. K13. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян. Я.Суликовский. Оптимизация параметров магнитной системы многоцелевого изохронного циклотрона АИЦ-144. // ПИСЬМА В ЭЧАЯ, 2012. Т.9, № 2(172). Стр. 235-246. К14. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян. Я.Суликовский. Математическое моделирование основного режима работы многоцелевого изохронного циклотрона АИЦ-144. // ПИСЬМА В ЭЧАЯ, 2012. Т.9, № 4-5(174-175). Стр. 647-653.

К15. И.Н.Киян. GaussDLL. Комплекс программ матрично-векторных операций и решения СЛАУ с вещественными коэффициентами. // Библиотека готовых программ ОИЯИ JINRLIB:

http://wwwinfo.iinr.ru/programs/iinrlib/gaussdll/index.html.

Получено 21 июня 2012 г.

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

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

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

ВВЕДЕНИЕ.

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

2. Содержание диссертационной работы.

3. Постановка общей задачи и ретроспектива ее решения.

4. Описание основных параметров и магнитной системы АИЦ-144.

ГЛАВА 1. МАТЕМАТИЧЕСКИЕ МОДЕЛИ.

1.1. Постановка задачи по моделированию основного режима работы АИЦ-144.

1.2. Порядок проведения магнитных измерений и подготовки исходных данных.

1.3. Математические модели для расчета токов в главной, в двадцати концентрических катушках АИЦ-144 и обеспечение устойчивости решения.

1.4. Математическая модель для расчета токов в двух парах гармонических катушек АИЦ-144.

1.5. Численный контроль точности расчета изохронного магнитного поля.

1.6. Выводы.

ГЛАВА 2. РАЗРАБОТАННОЕ ПРОГРАММНОЕ ОБЕСПЕЧЕНИЕ.

2.1. Описание комплекса программ Cyclotron Operator HELP Program Complex 20042012.

2.2. Описание динамически связанной библиотеки GaussDLL.

2.3. Выводы.

ГЛАВА 3. ПРИМЕНЕНИЕ КОМПЛЕКСА ПРОГРАММ CYCLOTRON OPERATOR HELP PROGRAM COMPLEX 2004-2012.

3.1. Моделирование основного режима работы АИЦ-144: протоны, кинетическая энергия вывода Ek ~ 60±0,5 МэВ, частота ВЧ-генератора Frj- = 26,15;.26,26 МГц.

3.2. Результаты запуска и эксплуатации АИЦ-144 в основном режиме его работы

3.3. Выводы.

Во введении приводятся общая характеристика и краткое содержание диссертационной работы, ставится задача на разработку новой методики математического моделирования режимов работы многоцелевого изохронного циклотрона АИЦ-144, позволяющей переходить на работу в моделируемых режимах без остановки и разборки ускорителя для проведения дополнительных магнитных измерений, дается ретроспектива решения поставленной задачи, приводится описание основных параметров и магнитной системы АИЦ-144 [7].

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

В диссертации изложены результаты разработок и исследований, выполненных соискателем в Лаборатории ядерных проблем Объединенного института ядерных исследований (ЛЯП ОИЯИ) в период с 2004 по 2012 год в рамках проблемно-тематических планов научно-исследовательских работ ОИЯИ (№ 03-2-1102-2010/2012 "Совершенствование Фазотрона ЛЯП (ОИЯИ) и разработка циклотронов для физических и прикладных исследований") и протоколов о международном сотрудничестве, осуществляющемся между ОИЯИ и Институтом ядерной физики Польской академии наук (ИЯФ ПАН) (№3970-2-10/12 "Протокол о выполнении совместной научно-исследовательской работы"). Работа выполнена при финансовой поддержке РФФИ, грант № 10-01-00467-а.

Актуальность. В настоящее время одним из главных направлений в области развития ускорительной техники является разработка, производство и эксплуатация универсальных изохронных циклотронов. Они находят применение как в области научных исследований (производство радиоизотопов для исследований в области ядерной спектроскопии, электроядерная энергетика), так и в области медицины (производство радиоизотопов для различных видов томографии, протонная терапия раковых опухолей): DC-60 [1], AIC-144 [2], VINCY [3]. Циклотроны данного вида являются многоцелевыми и могут использоваться для ускорения частиц с различным отношением числа стандартных единиц атомной массы к числу элементарных зарядов (A/Z). Энергия вывода в таких машинах может плавно варьироваться для каждого типа ускоряемых частиц.

Диссертационная работа посвящена созданию и применению на практике новой методики математического моделирования режимов работы многоцелевого изохронного циклотрона АИЦ-144, расположенного в ИЯФ ПАН и предназначенного как для проведения протонной терапии меланомы глаза у пациентов (при работе в основном режиме), так и для производства радиоизотопов, используемых в научных исследованиях в области ядерной спектроскопии. АИЦ-144 был разработан группой сотрудников ИЯФ ПАН под руководством д-ра Е.Схвабэ [7]. Система вывода АИЦ-144 была разработана группой сотрудников ОИЯИ при участии сотрудников ИЯФ ПАН [9]. Рассматриваемая методика позволяет моделировать новые режимы работы АИЦ-144 и переходить на работу с ними без остановки и разборки циклотрона для проведения дополнительных магнитных измерений. Новая методика содержит в себе набор алгоритмически связанных методов для расчета токов в главной катушке, в двадцати концентрических катушках, в двух парах гармонических катушек, а также для численного контроля точности расчета изохронного магнитного поля, рассчитываемого по аналитическим формулам М.М. Гордона (Прил., п.1), [4]. Общий алгоритм моделирования новых режимов работы АИЦ-144 включает в себя расчет токов во всех катушках циклотрона [К6]-ь[К13], численный контроль точности расчета изохронного магнитного поля [К14], аналитический [15], [КЗ] и численный [27] расчет фазового движения и частот свободных колебаний ускоряемых ионов. В общий алгоритм входят рабочие алгоритмы, позволяющие по наибольшему значению тока выведенного пучка проводить:

• Улучшение фазового движения ускоряемых частиц, которое проводится за счет регулировки амплитуды центрального бампа результирующего магнитного поля циклотрона (от токов в главной катушке и в двадцати концентрических катушках) [К13].

• Оптимизацию параметров магнитной системы циклотрона (регулировку амплитуды и фазы первой гармоники рабочего магнитного поля ускорителя (от токов в главной катушке, в двадцати концентрических катушках и в двух парах гармонических катушек), проводимую на радиусе заброса пучка ускоренных ионов в электростатический дефлектор) [К 13].

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

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

Намного более удобным представляется подход, когда новый режим работы моделируется без остановки и разборки циклотрона для проведения дополнительных магнитных измерений. Моделирование в этом случае осуществляется на основе уже имеющихся измеренных и пересчитанных карт исходных магнитных полей. Такой подход требует разработки и развития, как ряда новых математических моделей, так и ряда программных и рабочих алгоритмов. При использовании нового подхода, время перехода на моделируемый режим работы значительно сокращается. Например, переход на основной режим работы АИЦ-144, осуществленный в июне 2009 года, занял всего шесть часов. Сюда вошло время на аналитические и численные расчеты, измерение тока внутреннего и внешнего пучка ускоренных ионов, а также время на улучшение фазового движения ускоряемых частиц. Еще два часа ушло на оптимизацию параметров магнитной системы циклотрона, проведенную на ускорителе в апреле 2011 года с помощью разработанного комплекса программ Cyclotron Operator HELP Program Complex 2011.

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

Цель диссертационной работы заключается:

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

2. В создании на языке программирования С++ в среде MS Visual Studio .NET комплекса программ Cyclotron Operator HELP Program Complex 2004-2012 и динамически связанной библиотеки GaussDLL, реализующих разработанные математические модели и алгоритмы.

3. В практическом применении разработанного программного обеспечения при компьютерном моделировании новых режимов работы АИЦ-144 (В том числе основного режима работы циклотрона: протоны, кинетическая энергия вывода Ек ~ 60±0,5 МэВ, частота ВЧ-генератора Frf = 26,15;.26,26 МГц).

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

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

1. Разработана новая методика математического моделирования режимов работы многоцелевого изохронного циклотрона АИЦ-144, позволяющая переходить на работу в моделируемых режимах без остановки и разборки ускорителя для проведения дополнительных магнитных измерений [К1]-^[К14].

2. Впервые созданы математическая модель и алгоритм, которые позволяют проводить расчет тока в главной катушке АИЦ-144 на основе карт исходных магнитных полей, измеренных для трех уровней тока в катушке возбуждения главного магнита [К9].

3. Разработан новый метод расчета токов в двадцати концентрических катушках для рассчитанного уровня тока в главной катушке АИЦ-144, основанный на минимизации квадратичного функционала, в состав которого введена штрафная функция. Данная функция предназначена для направленного покомпонентного ввода вектора решения (вектора нормированных значений оптимизированных токов в концентрических катушках) в границы области допустимых значений. Ввод найденного вектора решения в рамки заданных граничных значений осуществляется в направлении, определяемом в соответствии с расположением зон влияния концентрических катушек на основное магнитное поле АИЦ-144 (от центра к краю циклотрона) [К6]н-[К10].

4. Впервые созданы математическая модель и алгоритм, которые позволяют проводить расчет токов в двух парах гармонических катушек для рассчитанного уровня тока в главной катушке АИЦ-144. Расчет проводится путем решения неоднородной системы нелинейных алгебраических уравнений, решаемой итерационным способом [К11]-н[К13].

5. Усовершенствован метод контроля точности расчета изохронного магнитного поля, рассчитываемого по аналитическим формулам М.М.Гордона. Показано, что для основного режима работы АИЦ-144 значение модуля относительной ошибки частоты обращения протонов во всем диапазоне изохронизации формируемого магнитного поля составляет \ aF0\<3E-4, что сопоставимо по своей величине со значением модуля относительной погрешности измерения карт исходных магнитных полей, измеряемых с помощью датчика Холла, | ams \ =1+2Е-4 [К14].

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

1. На языке программирования С++ в среде MS Visual Studio .NET в стандартах DBA (Dialog-based Application) и SDI (Single Document Interface) создано программное обеспечение, реализующее разработанные математические модели и алгоритмы:

• Комплекс программ Cyclotron Operator HELP Program Complex 2004-2012, установленный на сервере АИЦ-144 и используемый для моделирования новых режимов работы циклотрона.

• Динамически связанная библиотека GaussDLL, включающая в свой состав комплекс программ матрично-векторных операций и решения СЛАУ с вещественными коэффициентами, размещенная на сайте ОИЯИ в библиотеке готовых программ HNRLIB [Kl5].

2. С помощью комплекса программ Cyclotron Operator HELP Program Complex 2009 для АИЦ-144 было выполнено компьютерное моделирование основного режима его работы, предназначенного для ускорения пучка протонов с необходимыми для протонной терапии меланомы глаза параметрами: кинетической энергией вывода Ek.ext ~ МэВ, значением и стабильностью тока выведенного пучка h„ = 20+40 нА, cjlh„~±5 % соответственно [К7МК10].

3. Смоделированный с помощью созданного программного обеспечения основной режим работы АИЦ-144 был установлен на циклотроне в июне 2009 года без эмпирической подстройки токов в концентрических катушках и частоты ВЧ-генератора. Благодаря эффективной работе сотрудников Циклотронного отдела ИЯФ ПАН и высокой стабильности установленного режима работы, подразделение протонной терапии ИЯФ ПАН впервые на протяжении ряда лет получило возможность бесперебойной работы и использования пучка выведенных протонов:

• Для настройки медицинских приборов и аппаратуры, которая проводилась в течение около двух лет с момента установки смоделированного режима работы на АИЦ-144.

• Для проведения успешной протонной терапии меланомы глаза у первой партии пациентов, впервые проведенной в Польше на АИЦ-144 в период с февраля по апрель 2011 года [К12НК14].

4. С помощью комплекса программ Cyclotron Operator HELP Program Complex 2011 для основного режима работы АИЦ-144 в апреле 2011 года была проведена оптимизация параметров магнитной системы циклотрона [К 13]. Выведенный пучок ускоренных протонов был успешно использован в ИЯФ ПАН для производства радиоизотопов, нарабатываемых с целью проведения научных исследований в области ядерной спектроскопии.

5. С помощью комплекса программ Cyclotron Operator HELP Program Complex 2012 в июле 2012 года на АИЦ-144 было проведено моделирование двух вариантов основного режима его работы: протоны, кинетическая энергия вывода Ек~ 60/60,5 МэВ, частота ВЧ-генератора Frj = 26,155/26,26 МГц. Моделирование было выполнено с учетом результатов последних магнитных измерений, проведенных на АИЦ-144 в сентябре и октябре 2011 года. Измерение тока пучка выведенных протонов проводилось в ионоводе после прохождения пучком двух корректирующих магнитов и двух квадрупольных линз. Коэффициент вывода пучка ускоренных ионов рассчитывался по измеренным значениям как отношение тока внешнего пучка в ионоводе к току внутреннего пучка на конце плато его значений перед электростатическим дефлектором (на радиусе R = 61 см). После проведения оптимизации параметров магнитной системы АИЦ-144 коэффициент вывода пучка ускоренных ионов для указанных вариантов основного режима работы циклотрона составил Кех, = 15+2 % и Кех1 = 20±2 % соответственно (Прил., п.2), что полностью соответствует установленным техническим требованиям. Практическая значимость разработанных методов и комплексов программ документально подтверждена: отчетами о проведенных на АИЦ-144 в период с 2004 по 2012 год физических экспериментах, подписанными начальником ЦО ИЯФ ПАН (Прил., п.2); письмом от дирекции ИЯФ ПАН к дирекции ЛИТ ОИЯИ BD/47/2011 от 27 апреля 2011 года об успешных результатах лечения раковой опухоли глаза у первых девяти пациентов (Прил., п.З). Новая методика и разработанное программное обеспечение, используемые в настоящее время для моделирования режимов работы

АИЦ-144, могут найти широкое практическое применение в процессе эксплуатации аналогичных ускорителей, расположенных в различных российских и зарубежных научно-экспериментальных центрах.

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

• VI Международный семинар по проблемам ускорителей заряженных частиц памяти В.П. Саранцева, 8-10 сентября, 2005. Алушта, Украина.

• XXI Russian Particle Accelerator Conference 2008 (RuPAC 2008), September 28 -October 3, 2008. Zvenigorod, Russia.

• The International Conference Mathematical Modeling and Computational Physics 2009 (MMCP 2009), July 7-11, 2009. Dubna, Russia.

• Всероссийская конференция (с международным участием): Информационно-телекоммуникационные технологии и математическое моделирование высокотехнологичных систем. 18-22 апреля, 2011. Москва, Россия.

• IX Международный семинар по проблемам ускорителей заряженных частиц памяти В.П. Саранцева 17-21 сентября, 2011. Алушта, Украина.

Результаты, включенные в диссертацию, докладывались на семинарах ИЯФ ПАН, Лаборатории информационных технологий и Лаборатории ядерных проблем Объединенного института ядерных исследований (ЛИТ ОИЯИ, ЛЯП ОИЯИ). Публикации. Материалы диссертации опубликованы в 14 печатных работах, из которых 7 статей в рецензируемых научных журналах, рекомендуемых ВАК (в одном зарубежном и трех российских) [К2], [К5], [Кб], [К8], [К9], [К13], [К14], 3 статьи в научных сборниках [К10]н-[К12], 2 статьи в сообщениях ОИЯИ [КЗ], [К4] и 2 аннотации докладов в материалах конференций [Kl], [К7], а также в виде 1 распространяемого программного продукта [К 15].

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

Достоверность представленных в диссертации результатов моделирования новых режимов работы АИЦ-144 подтверждается:

• Результатами тестирования разработанного программного обеспечения: комплекса программ Cyclotron Operator HELP Program Complex 2004-2012 и динамически связанной библиотеки GaussDLL.

• Сопоставлением результатов компьютерного моделирования основного режима работы АИЦ-144, выполненного с помощью различных версий разработанного программного обеспечения.

• Методическими результатами контроля точности расчета изохронного магнитного поля, полученными путем сопоставления аналитических и численных оценок, выполненных на основе использования расчетной карты результирующего магнитного поля АИЦ-144.

• Успешными результатами научно-методических натурных экспериментов, проведенных на АИЦ-144 в период с 2004 по 2012 год. (Отчеты, удостоверенные подписью начальника ЦО ИЯФ ПАН (Прил., п.2)).

• Успешными результатами проведения протонной терапии меланомы глаза у первой партии пациентов, впервые проведенной в Польше на АИЦ-144 в период с февраля по апрель 2011 года. (Письмо от дирекции ИЯФ ПАН к дирекции ЛИТ ОИЯИ BD/47/2011 от 27 апреля 2011 года (Прил., п.З)).

Структура и объем диссертации. Диссертация состоит из введения, трех глав, заключения и приложения. Список работ, опубликованных по теме диссертации, содержит 15 наименований. Полный объем диссертации составляет 135 страниц машинописного текста, включая 8 таблиц и 47 рисунков.

Результаты работы, представленные в автореферате, своевременно опубликованы в статьях, полностью раскрывающих его содержание. Полученные результаты неоднократно докладывались И.Н. Кияном на российских и международных конференциях. В целом, проделанная автором работа отличается завершенностью, имеет научную новизну и большую практическую значимость, подтвержденную документально: отчетами о проведенных на АИЦ-144 в период с 2004 по 2012 год физических экспериментах, подписанными начальником ЦО ИЯФ ПАН; письмом от дирекции ИЯФ ПАН к дирекции ЛИТ ОИЯИ BD/47/2011 от 27 апреля 2011 года об успешных результатах проведения протонной терапии меланомы глаза у первых девяти пациентов. Разработанная методика и созданное программно§а5$^пенение используются в настоящее время для мо^—*—**"■■— ••—-----------. ami.lWt1 угемо

Заместитель директора ИЯФ ПАН по научной работе профессор

Начальник ЦО ИЯФ ПАН ul. Radzikowskii'K0152,31-342 Krakftw. POLAND www.ifj.ctlu.pl switchboard: 12 662 8000

СПИСОК РАБОТ, ОПУБЛИКОВАННЫХ ПО ТЕМЕ ДИССЕРТАЦИИ

К1. И.Н.Киян. Новый метод моделирования режимов работы многоцелевого изохронного циклотрона. // 6-ой Международный семинар памяти профессора В.П.Саранцева, 2005. Алушта, Украина. Сборник аннотаций D9-2005-123. Дубна 2005. С. 16.

К2. И.Н.Киян. Новый метод моделирования режимов работы многоцелевого изохронного циклотрона. // ПИСЬМА В ЭЧАЯ, 2006. Т.З, N7(136). Стр. 49-54.

КЗ. И.Н.Киян, Р.Тарашкевич. Анализ области применения аналитических формул для расчета частот свободных колебаний в изохронных циклотронах. // Сообщение ОИЯИ Р9-2005-114, Дубна 2005.

К4. И.Н.Киян, Р.Тарашкевич. Методы расчета режимов работы и параметров изохронного циклотрона. // Сообщение ОИЯИ Р9—2006^40, Дубна 2006.

К5. Igor N. Kiyan and others. Simulation of operation modes of isochronous cyclotron by a new iterative method. // NUKLEONIKA, 2007. International Journal of Nuclear Research, V.52 (1). Pp. 29-34.

Кб. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян, Я.Суликовский. Моделирование требуемых режимов работы и анализ их устойчивости для многоцелевых изохронных циклотронов. // ПИСЬМА В ЭЧАЯ, 2009. Т.6, № 6(155). Стр. 805813.

К7. I.V.Amirkhanov, G.A.Karamysheva, I.N.Kiyan, J.Sulikowski. Mathematical Aspects of Modeling of Required Operation Modes of Multi Purpose Isochronous Cyclotrons. // The International Conference: Mathematical Modeling and Computational Physics, 2009. Dubna, Russia. Book of Abstracts D11-2009-88. Dubna 2009. Pp. 203-204.

K8. I.V.Amirkhanov, G.A.Karamysheva, I.N.Kiyan, J.Sulikowski. Mathematical Aspects of Modeling of Required Operation Modes of Multi Purpose Isochronous Cyclotrons. // BULLETIN of PFUR. Series Mathematics. Information Sciences. Physics, 2010. No 2 (2). Pp. 99-103.

K9. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян, Я.Суликовский. Методика математического моделирования режимов работы многоцелевых изохронных циклотронов. // ВЕСТНИК ТГУ: Серия: Прикладная математика, 2010. №9, Выпуск 1 (14). Стр. 55-68.

К10. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян, Я.Суликовский. Улучшенная математическая модель расчета токов в концентрических катушках коррекции основного магнитного поля многоцелевого изохронного циклотрона. // Scientific report 2008-2009, LIT JINR. 2009-196. Pp. 109-112.

Kll. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян, Я.Суликовский. Решение обратной задачи расчета токов в гармонических катушках многоцелевого изохронного циклотрона. // Информационно-телекоммуникационные технологии и математическое моделирование высокотехнологичных систем. Тезисы докладов Всероссийской конференции с международным участием. 18-22 апреля 2011 года. Российский университет дружбы народов. Москва. 2011. Стр. 253-255.

К12. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян, Я.Суликовский. Математическая модель расчета токов в двух парах гармонических катушек многоцелевого изохронного циклотрона. Контроль точности расчета изохронного магнитного поля. // Scientific report 2010-2011, LIT JINR. 2011-130. Pp. 135-138.

K13. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян, Я.Суликовский. Оптимизация параметров магнитной системы многоцелевого изохронного циклотрона АИЦ-144. // ПИСЬМА В ЭЧАЯ, 2012. Т.9, № 2(172). Стр. 235-246.

К14. И.В.Амирханов, Г.А.Карамышева, И.Н.Киян, Я.Суликовский. Математическое моделирование основного режима работы многоцелевого изохронного циклотрона АИЦ-144. // ПИСЬМА В ЭЧАЯ, 2012. Т.9, № 4-5(174-175). Стр. 647653.

К15. И.Н.Киян. GaussDLL. Комплекс программ матрично-векторных операций и решения СЛАУ с вещественными коэффициентами. // Библиотека готовых программ ОИЯИ JINRLIB: http://wwwinfo.iinr.ru/programs/iinrlib/gaussdll/index.html.

ЗАКЛЮЧЕНИЕ

На защиту выносятся следующие основные результаты диссертационной работы:

1. Разработана новая методика математического моделирования режимов работы многоцелевого изохронного циклотрона АИЦ-144, позволяющая переходить на работу в моделируемых режимах без остановки и разборки ускорителя для проведения дополнительных магнитных измерений [К1]-ь[К14].

2. Впервые созданы математическая модель и алгоритм, которые позволяют проводить расчет тока в главной катушке АИЦ-144 на основе карт исходных магнитных полей, измеренных для трех уровней тока в катушке возбуждения главного магнита [К9].

3. Разработан новый метод расчета токов в двадцати концентрических катушках для рассчитанного уровня тока в главной катушке АИЦ-144, основанный на минимизации квадратичного функционала, в состав которого введена штрафная функция, предназначенная для направленного покомпонентного ввода вектора решения в границы области допустимых значений [К6]н-[К10].

4. Впервые созданы математическая модель и алгоритм, которые позволяют проводить расчет токов в двух парах гармонических катушек для рассчитанного уровня тока в главной катушке АИЦ-144. Расчет проводится путем решения неоднородной системы нелинейных алгебраических уравнений, решаемой итерационным способом [К11]-^[К13].

5. Усовершенствован метод контроля точности расчета изохронного магнитного поля, рассчитываемого по аналитическим формулам М.М.Гордона. Показано, что для основного режима работы АИЦ-144 значение модуля относительной ошибки частоты обращения протонов во всем диапазоне изохронизации формируемого магнитного поля составляет | <7F0 \<ЗЕ-4, что сопоставимо по своей величине со значением модуля относительной погрешности измерения карт исходных магнитных полей, измеряемых с помощью датчика Холла, | сгш \ =1+2Е-4 [К14].

6. На языке программирования С++ в среде MS Visual Studio .NET в стандартах DBA (Dialog-based Application) и SDI (Single Document Interface) был создан комплекс программ Cyclotron Operator HELP Program Complex 2004-2012, реализующий разработанные математические модели и алгоритмы.

7. С помощью комплекса программ Cyclotron Operator HELP Program Complex 2009 для АИЦ-144 было выполнено компьютерное моделирование основного режима его работы, предназначенного для ускорения пучка протонов с необходимыми для протонной терапии меланомы глаза параметрами [К7]-ь[К10].

8. Смоделированный с помощью созданного программного обеспечения основной режим работы АИЦ-144 был установлен на циклотроне в июне 2009 года без эмпирической подстройки токов в концентрических катушках и частоты ВЧ-генератора. Установленный режим работы был использован при проведении успешной протонной терапии меланомы глаза у первой партии пациентов, впервые проведенной в Польше на АИЦ-144 в период с февраля по апрель 2011 года [К12ИК14].

9. С помощью комплекса программ Cyclotron Operator HELP Program Complex 2011 для основного режима работы АИЦ-144 в апреле 2011 года была проведена оптимизация параметров магнитной системы циклотрона [К 13]. Выведенный пучок ускоренных протонов был успешно использован в ИЯФ ПАН для производства радиоизотопов, нарабатываемых с целью проведения научных исследований в области ядерной спектроскопии.

10. С помощью комплекса программ Cyclotron Operator HELP Program Complex 2012 в июле 2012 года на АИЦ-144 было проведено моделирование двух вариантов основного режима его работы: протоны, кинетическая энергия вывода Ек~ 60/60,5 МэВ, частота ВЧ-генератора Frf= 26,155/26,26 МГц. Указанные магнитной системы АИЦ-144 коэффициент вывода пучка ускоренных ионов для указанных вариантов основного режима работы циклотрона составил Кт = 15 ±2 % и Кех1 = 20±2 % соответственно (Прил., п.2), что полностью соответствует установленным техническим требованиям.

Благодарности

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

Особо хочу поблагодарить дирекцию ИЯФ ПАН, сотрудников Циклотронного Отдела ИЯФ ПАН и, лично, начальника ЦО г-на Я. Суликовского за организацию и помощь в проведении на АИЦ-144 ряда физических экспериментов.

Отдельно хочу поблагодарить дирекцию ЛЯП ОИЯИ за организацию международного научного сотрудничества с ИЯФ ПАН и дирекцию ЛИТ ОИЯИ за оценку выполненной работы.

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

1. Ввод исходных данных и расчет численных констант на основе обозначений, принятых в статье М.М. Гордона [4]: а = со, б = С о тхс

Ъ = т, гАл /л 1

С, г = Я мап > где с - это скорость света, со0 - круговая частота обращения частицы, т0 - масса частицы, тх - стандартная атомная единица массы, д - заряд частицы, е -элементарный электрический заряд, г - средний радиус замкнутой равновесной орбиты, - его стартовое значение.

2. Гармонический анализ введенной карты поля для расчета коэффициентов Фурье: Нп(г), Сп(г) и параметров азимутальной вариации магнитного поля: Кп(г) и К'п(г): г42 С

К(г) = 2

Кп(г)+г3 ■ чСу с1г с1г п = 1- N, г=1,2,.5, где N - периодичность магнитной структуры. 3. Расчет параметра ро в азимутально-симметричном магнитном поле: г

Д>= а ■ ^

4. Расчет параметров среднего магнитного поля: 1

4=4-.' М' =

5. Расчет поправочного коэффициента §

4-1.

1 8

1-2:

Ьп2+\)-Кп п2-М')2 = /'=7,2,.5

6. Оптимизация параметра /? по методу Ньютона:

А,

1 - + — ¿Л

1 +

2М

7. Если |£/?|<106, то переход к пункту 8. Если неравенство не соблюдается, то возврат к пункту 4.

8. Расчет релятивистского импульса частицы:

V, 2

7+1

9. Пересчет параметров среднего магнитного поля: М' = \ + р2

М" = 3 • (1 + /?2)- (1 + 2/?2)

10. Расчет поправочного коэффициента £: е^2

I 1 мх |«2 -М' 2 („2-М'Н п = 1-Ы, / =1,2,.5

11. Расчет изохронного магнитного поля с учетом азимутальной вариации магнитного поля (флаттера): Р

В(г) = С

V*- У

12. Расчет очередного значения среднего радиуса замкнутой равновесной орбиты: г = г + 5г

13. Если г<К/тп11, то возврат к пункту 2. Если неравенство не соблюдается, то переход к пункту 14.

14. Конец расчета изохронного магнитного поля с учетом флаттера.

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

16.12.2004 года.

16.12.2004 года в Институте Ядерной Физики Польской Академии Наук (ИЯФ ПАН) в Кракове на изохронном циклотроне AIC144 был проведен эксперимент по проводке пучка протонов от источника ионов до радиуса вывода. Эксперимент был проведён для двух режимов работы под руководством начальника Циклотронного Отдела (ЦО) МТалаха на основе расчётов, выполненных Кияном И.Н. с помощью комплекса программ, разработанного Кияном И.Н. на С++ и Transact SQL. Эксперимент прошел успешно. Токи пучков были получены после установки токов катушек изохронного циклотрона в течение 10 минут каждый при минимальной эмпирической подстройке тока главной катушки без регулировки токов катушек коррекции и без регулировки частоты RF генератора. Результаты эксперимента приведены в приложении N1 и N2

Директор ИЯФ: профессор Ежабек М. <паллий

АИЦ-144 подтверждение о проведении эксперимента

1. http ://fls2. i inr.ru/flnr/dc60 .html.2. http ://www. i fj. edu. p 1/str/dc/index. html.3. http://www.vinca-at.org/tesla.html.

2. M.M. Gordon. Calculation of Isochronous Fields for Sector-Focused Cyclotrons // Particle Accelerators, 1983. Vol. 13. Pp. 67-84.

3. В.Н. Аносов, Ю.Н. Денисов, В.П. Дмитриевский, В.В. Кольга. Исследование движения и коррекции фазы сгустка в изохронном циклотроне. // Препринт ОИЯИ, Р9-8496, Дубна, 1974.

4. Е. Схвабэ и др. Магнитная система автоматического изохронного циклотрона AIC-144. // Proceedings of the International Seminar on Isochronous Cyclotron Technique. Institute of Nuclear Physics, Krakow, 13th—18th November, 1978. Pp. 237256. Krakow, 1978.

5. C.M. Davisson, R.B. Theus. Linear Programming Vs Least Square Method in Determing the Best Trim Coil Currents for the NRL cyclotron. // Particle Accelerators, 1976. Vol. 7. Pp. 103-110.

6. Nikolaj A. Morozov and others. A new extraction system for the upgraded AIC-144 cyclotron. // NUKLEONIKA, 2001. International Journal of Nuclear Research, V.46 (2). Pp. 51-57.

7. В.Ф. Формалев, Д.Л. Ревизников. Численные методы. // Учебное пособие. ФизМатЛит, Москва, 2004.

8. Н.С. Бахвалов, Н.П.Жидков, Г.М. Кобельков. Численные методы. // Учебное пособие. БИНОМ. Лаборатория знаний. Москва. 2006.

9. Н.С. Пискунов. Дифференциальное и интегральное исчисление. // Учебники для втузов. В 2-х т. T.I, Т.Н. Интеграл-Пресс. Москва. 2004.

10. А.А. Самарский. Введение в численные методы. // Учебники для втузов. Лань. Санкт Петербург. 2009.

11. A.A. Самарский, П.Н. Вабищев, Е.А. Самарская. Задачи и упражнения по численным методам. // Учебное пособие. Москва. КомКнига. 2007.

12. Т. Stammbach. Cyclotrons. // Joint Universities Accelerator School. Archamps-France. 1997.

13. Дж. Ливингуд. Принципы работы циклических ускорителей. // Издательство иностранной литературы. Москва. 1963.

14. Е.Г. Комар. Основы ускорительной техники. // Атомиздат. Москва. 1975.

15. Н. Schildt. С++ Nuts & Bolts for Experienced Programmers. // Osborne McGraw-Hill, Berkeley. 1995.

16. Г. Шилдт. Самоучитель С++. 3-е издание, переработанное и дополненное. // bhv, Санкт-Петербург. 2001.

17. О. Gurevich, N. Gurevich. Teach Yourself Visual С++ 4 in 21 Days. // SAMS Publishing, Indianapolis. 1996.

18. K. Gregory. Special Edition Using Visual С++ 6. // Que, Indianapolis. 1998.

19. А. Черносвитов. Курс MCSD Visual С++ 6 и MFC. // Питер, Санкт-Петербург.

20. А. Черносвитов. Visual С++ 7. Учебный курс. // Питер, Санкт-Петербург. 2001.

21. Е.В. Самсонов Расчет динамики пучка протонов в магнитном поле карты 26p25Ms циклотрона АИЦ-144. // Внутренний отчет. Краков. 2007. Стр. 1-7.

22. И.М. Ситник. Модификация пакета программ минимизации FUMILI. Памяти проф. И.Н. Силина. // Библиотека готовых программ ОИЯИ JINRLIB: http://wwwinfo.iinr.ru/programs/iinrlib/fumilim/index.html.

23. С.А. Костромин, Г.А. Карамышева, Е.В. Самсонов, И. Йонген. Влияние магнитного поля от ускоряющей системы на динамику пучка в циклотроне. // ПИСЬМА В ЭЧАЯ, 2012. Т.9, № 4-5(174-175). Стр. 635-639.

24. И.Н. Киян, С.А. Костромин. Отчет по работе с АИЦ-144 (Июль-август 2012, Краков-Дубна). // Внутренний отчет. Дубна. 2012. Стр. 1-12.

25. JI.A. Сена. Единицы физических величин. Учебно-справочное руководство. // Наука. Москва. 1988.

26. X. Кухлинг. Справочник по физике. // Мир. Москва. 1982.

27. Г. Корн, Т. Корн. Справочник по математике для научных работников и2000.