автореферат диссертации по энергетике, 05.14.03, диссертация на тему:Методы и программные комплексы формирования систем константного обеспечения расчетов активации материалов и выходов газообразных продуктов в электро-ядерных установках

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

00501365^

Л

Т'

КУПЦОВ ИЛЬЯ СЕРГЕЕВИЧ

МЕТОДЫ И ПРОГРАММНЫЕ КОМПЛЕКСЫ ФОРМИРОВАНИЯ СИСТЕМ КОНСТАНТНОГО

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

Специальность 05.14.03 - «Ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации»

АВТОРЕФЕРАТ

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

ОБНИНСК 2012

1 5 [;Ь\? Ши

005013652

Работа выполнена в Обнинском институте атомной энергетики - филиале федерального государственного автономного образовательного учреждения высшего профессионального образования «Национальный исследовательский ядерный университет «МИФИ»

Научный руководитель:

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

Ведущая организация

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

доктор физико-математических наук Титаренко Юрий Ефимович

кандидат физико-математических наук Блохин Анатолий Иванович

Институт Физики Высоких Энергий (ГНЦРФ-ИФВЭ)

Защита состоится «30» марта 2012 г. в 10 часов на заседании диссертационного совета Д.201.003.01 при ГНЦ РФ Физико-энергетический институт имени А.И. Лейпунского в конференц-зале по адресу: 249033, г. Обнинск, пл. Бондаренко, 1.

С диссертацией можно ознакомиться в библиотеке ГНЦ РФ-ФЭИ. Автореферат разослан февраля 2012 г.

Ученый секретарь диссертационного совета

доктор технических наук, Т.Н. Верещагина

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

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

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

Концептуальные исследования свидетельствуют о перспективности электро-ядерного способа производства нейтронов для трансмутации дол-гоживущих радиоактивных отходов (ДРАО). Этот способ заключается в использовании подкртического реактора (называемого бланкетом) с внешним источником нейтронов в виде нейтронопроизводящей мишени и ускорителя протонов. В российской технической литературе этот тип установок называют электро-ядерными установками (ЭЛЯУ), в зарубежной - Accelerator Driven System (ADS).

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

• Развитие методов расчета сечений реакций глубокого расщепления в области промежуточных энергий протонов 20-2000 МэВ;

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

• Развитие и создание компьютерных программ расчета полной и парциальной активноти, остаточного энерговыделения в материалах мишени;

• Создание тестированных библиотек оценненых ядерных данных о взаимодействии протонов и нейтронов с энергиями до 3 ГэВ с материалами мишени, ускорителя и бланкета.

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

Объектом исследования являются методы и программные комплексы формирования систем константного обеспечения расчетов активации материалов и выходов газообразных продуктов в электро-ядерных установках материалы в энергетическом диапазоне от 50 МэВ до 1 ГэВ.

Цель диссертации - разработка и совершенствование методов и программных комплексов формирования систем константного обеспечения расчетов активации материалов и образования легких газообразных продуктов в перспективных ЯЭУ на основе методов теории принятия решений при многих критериях, усовершенствование расчетного кода CASCADEX, включающего развитие методов расчета сечений реакций глубокого расщепления, эмиссии кластеров и выходов легких ядер в мишени электроядерных установок (ЭЛЯУ), обновлению и расширению файлов библиотек оцененных ядерных данных HEPAD, HEAD.

Цель диссертации определяет постановку и решение следующих задач:

• Изучение особенностей, областей применения каскадно-испарительной модели CASCADEX;

• Реализация процедуры расчета выходов легких ядер на базе модели CASCADEX;

• Разработка алгоритма оценки высокоэнергетических ядерных данных на основании многофакторного анализа;

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

• Оценка моделей высокоэнергетических ядерно-ядерных взаимодействий на основании разработанных программных средств и методик.

На защиту выносятся следующие положения:

• Методы оценки активационных ядерных данных в высокоэнергетической области с использованием методов теории принятия решений при многих критериях;

• Программный комплекс CASCADEX 1.2, включающий в себя учет выходов легких ядер, определяющих накопление газообразных продуктов в мишени ЭЛЯУ и продуктов реакций глубокого расщепления;

• Интерактивные программные комплексы подготовки, верификции и валидации баз активационных ядерных данных в области энергий от 0,1 до 10 ГэВ;

• Результаты сравнительного анализа моделей ядерно-ядерных взаимодействий в области энергий от 50 МэВ до 1 ГэВ и данных натурного эксперимента для основных конструкционных материалов ЭЛЯУ (56Fe, Со, l84W, 209Bi ипа1РЬ).

Научная новизна данного диссертационного исследования заключается в том, что в работе:

• Впервые разработан интерактивный программный комплекс для подготовки, верификации и валидации активационных ядерных данных в области энергий от 50 МэВ до 10 ГэВ на основе каскадно-испаритаельных моделей, необходимых для расчетного обоснования проектов ЭЛЯУ;

• Впервые предложен подход к оценке активационных ядерных данных в области энергий от 50 МэВ до 10 ГэВ на основе методов теории принятия решений при многих критериях, позволяющий формировать самосогласованный набор систем константного обеспечения расчетов ЭЛЯУ;

• Впервые предложена модифицированная версия ПК CASCADEX 1.2, включающая в себя учет выходов легких ядер, определяющих накопление газообразных продуктов в мишени ЭЛЯУ и продуктов реакций глубокого расщепления на основе модифицированного подхода Хаузе-ра-Фешбаха, реализованного в TALYS 1.2;

• Впервые представлены результаты сравнительного анализа моделей ядерно-ядерных взаимодействий в области энергий от 50 МэВ до 1 ГэВ и данных натурного эксперимента для основных конструкционных материалов ЭЛЯУ (56Fe, 59Со, l84W, 209Bi и nalPb).

Практическая ценность работы состоите следующем:

• Проведена верификация и валидация программного комплекса CASCADEX на известных экспериментах и бенчмарках. Даны рекомендации по свободным параметрам данной модели для изотопов необходимых для константного обеспечения ЭЛЯУ;

• Подготовлены секции в библиотеки оцененных ядерных данных для ряда изотопов необходимых для расчета характеристик при проектировании и оптимизации эксплуатационных режимов ЭЛЯУ;

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

Научные положения, выводы и рекомендации, сформулированные в диссертации, обоснованы поскольку были апробированы в международном проекте «Benchmark of Spallation Models» проводимым МАГАТЭ, в работе использовались современные экспериментальные данные и теоретические представления о механизмах ядерных реакций в рассматриваемом энергетическом диапазоне, а также современные программы расчета ядерно-ядерных взаимодействий в высокоэнергетической области.

Личный вклад соискателя в представленную работу заключается:

• В разработке и реализации справочно-информационного интерактивного комплекса, подготовке библиотек ядерно-физических данных;

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

• Проведении всех верификационных и валидационных расчетов на известных экспериментах и бенчмарках с их последующим анализом для выявления возможности применения к решению задач разработки ЭЛЯУ;

• Оценке моделей на основании разработанных программных средств и методик;

• В обновлении и расширении файлов библиотек оцененных ядерных данных HEPAD, HEAD.

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

• Безопасность АЭС и подготовка кадров. X Международная конференция, Обнинск, 1 - 4 октября 2007 г.

• XV семинар по проблемам физики реакторов («Волга-2008»), Актуальные проблемы физики ядерных реакторов - эффективность, безопасность, нераспространение. Москва, 2-6 сентября 2008 года.

• Workshop on Accelerator Radiation Induced Activation (ARIA'08), Villigen, Switzerland, October 13 - 17,2008.

• Безопасность АЭС и подготовка кадров. XI Международная конференция, Обнинск, 29 сентября - 2 октября 2009 г.

• Научная сессия НИЯУ МИФИ-2009, Москва, 25-30 января 2009

• ND-2010 - International Conference on Nuclear Data for Science and Technology. Jeju Island, Korea, 26 - 30 April, 2010.

• Научная сессия МИФИ - 2010, Москва, 25-31 января, 2010.

• Научная сессия МИФИ - 2011, Москва, 02 - 05 февраля, 2011.

• Научная сессия МИФИ- 2012, Москва, 30 января - 4 февраля, 2012.

Результаты исследования представлены в статьях: по теме диссертации опубликовано 16 работ, в том числе 4 статьи в ведущих рецензируемых научных журналах, определенных ВАК, 12 - в материалах всероссийских и международных конференций и семинаров. Получено 2 авторских свидетельства.

Структура и объем диссертации. Диссертация состоит из введения, четырех глав основного текста, заключения, библиографического списка, включающего в себя 117 наименований. Работа изложена на 121 страницах с 30 иллюстрациями и 15 таблицами.

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

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

Глава I посвящена краткому описанию потребностей в высокоэнергетических ядерных данных в энергетике, медицине и физике твердого тела, основным современным программно-расчетным комплексам для получения, обработки, представления высокоэнергетических ядерных данных и обзору существующих библиотек ЯФД в высокоэнергетической области.

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

В таблице 1 дается краткий обзор основных приложений высокоэнергетических реакций.

Таблица 1. Применения реакций глубокого расщепления

Применения Существующие проекты

Высокоинтенсивный источник продуктов глубокого расщепления SINQ (Швецария), Ер = 0,6 ГэВ, мощность = 1-1,5 МВт

SNS (США), Ер = 1 ГэВ, мощность = 1-1,5 МВт

J-PARC (Япония), Ер = 3 ГэВ, мощность = 1 МВт, и 50 ГэВ адроный ускоритель

ESS (Европа), Ер = 1,5 ГэВ, мощность = 10 МВт

ADS (УУС), ADTT, трансмутация (эксперименты под бенчмарки) MEGAPIE, MUSE, RACE-ISU, Y ALINA - бустер и т.д.

Новые установки GSI-FA1R, SPIRAL, EURISOL и т.д.

Облучение материалов высокоинтенсивным потоком частиц Установки по облучению, производство редких изотопов и т.д.

Безопасность и радиационная защита Физическая защита, радиационные измерения, дозиметрия и т.д.

Разработка устройств детектирования Для текущих и будущих проектов по ускорительным системам

Астрофизика, космос

Деятельность МАГАТЭ в сфере малых протонных ускорителей и их применений

Обсуждаются проблемы, связанные с требованиями к материалам в системах ввода пучка и мишенях. Приведены перечни международных проектов - бенчмарков' для подготовки констант существующих и инновационных проектов ЭЛЯУ2.

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

• невозможность получения всех констант экспериментальным путем в связи с длительностью и большими материальными затратами на эксперимент;

• трудности, связанные с расчетной обработкой эксперимента;

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

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

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

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

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

Бенчмарк (англ. benchmark) — контрольная задача, необходимая для определения сравнительных характеристик производительности компьютерной системы. ЭЛЯУ — ЭЛектроЯдерная Установка.

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

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

• теоретический расчет ядерных данных;

• корректировка рассчитанных ядерных данных с доступными экспериментальными данными;

• выбор критериев сравнения данных;

• сравнение рассчитанных ядерных данных с экспериментальными ядерными данными с помощью выбранных критериев;

• формирование файлов оцененных ядерных данных на основании проведенного сравнения в установленном формате.

В §1.4 приведены примеры интерактивных систем оценки и хранения ядерных данных.

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

Программа TALYS предназначена для моделирования взаимодействий легких частиц (протонов, нейтронов, дейтронов, тритонов, 3Не- и альфа-частиц) с веществом. В TALYS заложены строгие физические модели для описания ядерных реакций, в связи с чем область применения программы ограничена диапазоном от ] кэВ до 200 МэВ. Для расширения энергетического диапазона используется статистическая модель программы CASCAD/INPE, способная моделировать высокоэнергетические взаимодействия от 100-150 МэВ до 3 ГэВ. Таким образом, основа программы CASCADEX представляет собой соединительный модуль между TALYS и CASCAD/INPE.

В главе приведены усовершенствования CASCADEX, включающие в себя интеграцию с последней версией кода TALYS 1.2, реализацию модулей расчета дважды дифференциальных сечений выходов легких ядер (р, п, d, t, 3Не, 4Не) на основе коалесцентной модели и модели движущегося источника, учет процесса образования кластеров. Дано описание интерактивной среды IViS CASCADeX, обеспечивающей возможность проведения многовариантных серийных расчетов, осуществления обоснованного выбора параметров моделей по результатам количественного сопоставления с экспериментом и включающей в себя широкий спектр средств обработки и представления результатов расчетов. Приведены результаты использования разработанного программного обеспечения для задач совершенствования константного обеспечения.

В §2.1 дается краткое описание общего алгоритма программ, рассчитывающих высокоэнергетические ядерно-ядерные взаимодействия. Приводятся основные характеристики и параметры программы CASCADEX. Описаны этапы развития программного комплекса CASCADEX.

В §2.2. рассматривается каскадно-испарительная модель на примере программы CASCADE, проходящая в две стадии. Первая, быстрая стадия, называется внутриядерным каскадом и описывается в рамках модели взаимодействия двух облаков ферми-газа. Данная стадия модели является статистической, поэтому необходим метод Монте-Карло для моделирования взаимодействий частиц внутри ядра мишени. «Каскадная» стадия модели завершается после того как, налетевшая частица полностью потеряла свою энергию или вылетела из ядра мишени, ядро при этом переходит в предравновеснос состояние. В предравновесии ядро перераспределяет избыточную энергию между нуклонами ядра, при этом существует большая вероятность, что вся энергия сконцентрируется на одном нуклоне, и он вылетит из ядра. Если этого не произойдет, то ядро перейдет в равновесное возбужденное состояние и снимет возбуждение путем испарения вторичных частиц. Процесс испарения частиц представляет собой вторую (медленную) стадию каскадного процесса.

Описание испарения частиц в рамках каскадно-испарительной модели связано с проблемой создания эффективного алгоритма для моделирования этого процесса методом Монте-Карло.

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

Принцип работы программы CASCADEX заключается в следующем. Информация о возбужденных ядрах, образующихся после стадии внутриядерного каскада, накапливается после каждой истории взаимодействия первичной частицы с ядром-мишенью. Таким образом, после прослеживания всех историй первичных частиц образуется «популяция» составных ядер, причем ее размеры и содержимое зависят от «пороговой» энергии, ниже которой каскадная стадия считается завершенной. Эта энергия является параметром модели и ограничена значением 200 МэВ.

В иллюстративных целях проведено сравнение расчетных и экспериментальных данных по выходам продуктов в реакциях под действием протонов с энергиями от 0,8 до 1,2 ГэВ. Для сравнения использовались также и результаты расчетов остаточных ядер, полученные по программам СЕМ03 и INCL4/ABLA, наиболее популярным в настоящее время. Полученные результаты выходов продуктов реакций умножались на одно и то же полное сечение неупругих взаимодействий протонов с ядрами для получения сечений образования ядер. Кривые массовых распределений остаточных ядер для различных ядер-мишеней и энергий первичных протонов приведены на рис.1: ядро 197Аи, облучаемое протонами с энергией 800 МэВ (а); ядро 208РЬ, облучаемое протонами с энергией 1 ГэВ (б); ядро 65Си, облучаемое протонами с энергией 1.2 ГэВ (в); ядро 9 Тс, облучаемое протонами с энергией 800 МэВ (г). Кривые, полученные по CASCADEX, обозначены пунктиром, сплошной линией и штрихом представлены расчеты по моделям СЕМ03 и INCL4/ABLA. Экспериментальные данные взяты из работ Ю.Е. Титаренко, F. Rejmund, J. В (2001-2002).

40 бо 80 10« 120 140 1б0 160 20с

Массовое число продуктов

Массовое число продуктом

б)

" Ркщ цыви аіі^ртм* гцвтп рісцсла

Н

X «

и

А

•хр Тігвівпісс «1 ві, 2С01 .- ссм-оі(мсмі'кг.ои)

-їмси*

-СА5САРЕХІСА5САРЕ-»ТА1.УЗ*

їй 20 зо 4и ьо Массовое число продуктов

В)

100

»

м ©

г

о

0.1

Р(Я°0

т

<хр. Тіаг«пісо «і «і,

- С-ЕиоЗОіС^РХ 2.вс> -ІИСІ-* (УСИРХ 3.вє» СА5СА0еХ (СД.гСАРЕ^ТАї'

10 20 30 40 50 60 70 00

Массовое число продуктов

Рис. 1, Массовые распределения сечений накопления ядер-продуктов реакции

Как видно из графиков на рис.1, сечения, рассчитанные по программе CASCADEX, находятся ближе к экспериментальным данным, чем сечения, рассчитанные по другим программам. Однако имеется существенное расхождение в области продуктов глубокого расщепления, особенно отчетливо это наблюдается для ядер-мишеней 6SCu и 99Тс, а также для легких фрагментов, образующихся при взаимодействии высокоэнсргетических протонов с тяжелыми ядрами.

В §2.3-2.7 приведено подробное описание всех модернизаций CASCADEX, необходимых для расширения функциональных и технических возможностей программы. Представлены примеры расчета дважды дифференциальных сечений в программе CASCADEX 1.2.

Основные пункты можно отметить как:

Оптимизация под Windows-платформу

Для улучшения производительности и стабильности CASCADEX код был протестирован на компиляторах Compaq Fortran 6.Х, Watcom Fortran, CygWin G95, Intel Visual Fortran 11 .X. Расчет скорости показал, что с опциями по умолчанию компилятор Intel Visual Fortran увеличивает производительность программы примерно на 10-15% на компьютерах одинаковой конфигурации.

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

Модели выходов легких ядер (коалесцентная и модель движущегося источника)

В CASCADEX был включен модуль расчета дважды дифференциальных сечений вылета легких ядер (d, t, 3Не, 4Не) с использованием двух моделей (модель движущегося источника и коалесцентная модель). Обе модели, несмотря на то, что являются феноменологическими, показали неплохое согласие с экспериментальными данными, и в настоящее время широко используются в подобных комплексах программ.

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

( V*"1

(4-я/З )-Р0]

f <тГ

d Edil Z\N\

' p /

Д2-тъ-{Е-Ес)У;

A

d2°7 ] dzdQ. J

где 2, и 2Р - атомные номера ядра-мишени и налетающей частицы; N,иNp-число нейтронов в ядре-мишени и налетающей частицы соответственно;

клона; Ес - кулоновский потенциал для а-частиц; dloprp' jdtpd€l - дважды дифференциальное сечение эмиссии неравновесных нейтронов или протонов; Еа - энергия испускаемой а -частицы, связанная с энергией вылетающего нуклона соотношением е0 = 4 • Ер - 2 • Ес; Р0 — радиус «коалисценции», параметр модели.

Равновесные спектры налетающих частиц (апж ) и дважды дифференциальное сечение эмиссии неравновесных нейтронов (d2csp"jdzpdQ.) вычислялись с помощью программы CASCADeX 1.2.

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

d2N А, г— [[-£-Z'£t + Е, • yftT-cos(Q)]\

dQdE ~ 0 'ехр { т J'

где N0 - нормировочная константа; Еж = E-Z-Ес- энергия отталкивания, обусловленная кулоновским взаимодействием; Ес - энергия кулоновского отталкивания, приходящаяся на единичный заряд; Z - заряд испущенной частицы; =m-V2j2 - кинетическая энергия частицы массой т в системе центра масс при скорости V; 7(G) = 70 ■ ехр(- 9/45,5) - температура ядра; 8 - угол детектирования.

Учет процессов образования кластеров

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

ст (а, А, Е) = о(А,Е) • Ща, А) ■ Р(а, А, Е), где g(A,E) - сечение неупругого взаимодействия протонов, имеющих энергию Е, с ядром А ; N(a, А) - эффективное число кластеров а в ядре А ; Р(а,А,Е) - вероятность нарушения связи кластера а с ядром А при энергии падающих протонов Е.

Разработана интерактивная система

Для целей автоматизации расчетов и обработки результатов, а также для подбора параметров моделей, была разработана система визуального моделирования IViS (Interactive Visual System). Данного рода программная система нацелена на решение проблем обработки, оценки, согласования, вали-дации, верификации и наглядного представления полученных данных как

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

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

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

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

В §3.1.1 описываются функциональные возможности СИИК, основными из которых являются

• поиск оптимальных параметров модели в каждой из областей;

• получение расчетных данных по каждой из моделей;

• выбор оптимальной модели;

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

В §3.1.2 обосновывается выбор языка программирования, в качестве которого была выбрана среда разработки Microsoft Visual Studio .Net 2008. В §3.2 проводится описание компонентов СИИК, включающих в себя - интерактивную среду (оболочку), объединяющую библиотеки высокоэнергетических ядерных данных, модели высокоэнергетических ядерных реакций, экспериментальные данные в высокоэнергетической области (EXFOR) и средства подключения данных из файлов;

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

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

стадии испарения, который наравне с другими моделями высокоэнергетических ядерных реакций входит в инструментальный набор СИИК;

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

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

В §33 предлагается детализированная структурная схема блоков н модулей СИИК.

§3.3.1. Блоки СИИК

Функционально СИИК можно разделить на три равноправных блока: управления, хранения, обработки.

Блок управления - основной блок СИИК, включающий в себя MDI-менеджер (MDI - Multiple Document Interface), менеджер решений и менеджер внутренней базы данных.

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

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

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

Блок хранения включает в себя две структурные единицы: база данных в формате ENDF/B (HEPAD, 1EAF) и XML-хранилище.

База данных в формате ENDF/B позволяет подключать файлы с актива-ционными данными.

XML-хранилище реализует функции хранения и быстрого доступа к файлам проектов и Решений.

Блок обработки данных

Для обработки данных СИИК включает в себя пять структурных блоков:

• подключение данных из базы EXFOR;

• расчет;

• подключение произвольных данных;

• обработка данных ENDF/B-формата;

• статистика.

Блок визуализации предназначен для наглядного представления и сравнения полученных данных. Визуализатор доступен практически из любой части СИИК. Обмен данными с визуализатором происходит по принципу клиент-сервер. Те части СИИК, из которых доступен визуализатор, отправляют соответствующие запросы и информацию по протоколу (TCP/IP).

§3.3.2. Модули СИИК

Модуль подключения данных из базы EXFOR позволяет пользователю использовать экспериментальные данные для последующего их применения в задачах сравнительного анализа. При этом экспериментальные данные могут быть получены и загружены в СИИК как через удаленную базу данных, расположенную на официальной сайте МАГАТЭ3 (www-nds.iaea.org/exfor), так и из файлов, расположенных на локальном компьютере пользователя.

Модуль расчета. В модуле расчета происходит работа с программными комплексами, включенными в СИИК. Пользователю необходимо выбрать модели, с помощью которых будут произведены соответствующие вычисления, а также задать общие для всех моделей входные данные. Кроме задания общих для всех моделей параметров необходимо задать параметры, характерные для каждой модели.

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

Модуль подключения произвольных данных необходим для подключения данных, полученных из других источников (например, данные, взятые из интернета или рассчитанные на программных комплексах, к которым у СИИК нет доступа).

Модуль обработки данных ENDF/B-формата. Его функцией является формирование файлов оцененных библиотек ядерно-физических данных и считывания данных из файлов готовых оцененных библиотек ядерно-физических данных в формате ENDF/B.

Модуль статистики. Для сравнения результатов расчетов по моделям с экспериментальными данными в настоящее время используются пять рекомендованных факторов согласия (F-, Н-, R-, L-, D-факторы) и два фактора,

3 МАГАТЭ — Международное агентство по атомной энергии.

предложенных авторами: корреляция и модифицированный Л-фактор Кт. Их описание приведено в табл. 2.

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

Таблица 2. Факторы согласия

N )| N Оценка интегральной близости к эксперименту, при условии, что данные могут сильно различаться

н аГр - аТ!с' Асг;хр 2 5 Показатель степени отражает допустимую степень компенсации малых значений одних слагаемых большими значениями других. Чем больше значение показателя, тем больше степень возможной компенсации

и ч* п аТ? ~ а?к о-;хр 3

ь ь = 1 N сак 2 Го-;*"-о-™"-" г Оценка интегральной близости к эксперименту при условии обеспечения примерно одинакового вклада различных областей

и 1 І-] Г- -|7 > — СЧІС 1 <Т1 .Дс7Н

Л 1 N са/с Оценка интегральной относительной близости к эксперименту

Сог — (=0 Оценка тенденций поведения сечений (Ковариа-ция)

Лот Лот = тах< СÙР'о?*]' Оценка максимального относительного отклонения сечений

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

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

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

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

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

В главе IV представлено применение разработанных комплексов для решения задач получения высокоэнергетических ядерных данных и последующей их оценки. Разработанное программное обеспечение было использовано для проведения расчетов в рамках проекта, организованного МАГАТЭ, по сравнительному анализу моделей реакций глубокого расщепления («Benchmark of Spallation Models»), а также обновлению библиотек активационных данных.

В §4.1 проводится определение оптим&пьных параметров модели CASCADeX 1.2 с целью дать рекомендации по областям применимости CASCADeX 1.2 и оценки значений свободных параметров модели в диапазоне энергий от 50 МэВ до 1 ГэВ и массовых чисел от 13 до 240. Была проведена серия расчетов, по результатам которых были определены значения свободных параметров обеспечивающих наилучшее согласование с экспериментом.

В качестве примера представлены оптимальные значения параметра пороговой энергии для выбранных пяти нуклидов. Экспериментальные значения были взяты из базы данных EXFOR, а также экспериментальных значений, представленных в проекте МАГАТЭ «Benchmark of Spallation Models» (табл. 3). Значение пороговой энергии равное нулю, означает, что в данной реакции предпочтительно использовать подход Вайтскопа-Ивинга вместо формализма Хаузера-Фешбаха.

Таблица 3. Оценка значений пороговой энергии

Реакция Ядро выхода Пороговая энергия (МэВ) Эксперимент (авторы) Год

58Ni(p, я) "Fe 10 G.F. Stein, S.J. Mills, F.M. Nortier, B.R.S. Simpson, B.R. Meyer 1990

"Fefc,*) 2lNe 0 K. Ammon, 1. Leya, B. Lavielle, E. Gilabert, J.C. David, U. Herpers, R. Michel 2008

55Mn(p, x) 36CI 0 Th. Sehiekel, F. Sudbrock, U. Herpers, M. Gloris 1996

mPb(p,x) 2WPb 50 Ю.Е. Титарснко 2006

204B|. 50 2006

72Ga 100 - 2006

,иІг 100 /» 2006

200-j-j 100 H 2006

,6Nb 100 2006

205Bi 50 ft 2006

0 2006

На основании количественных сопоставлений с экспериментами проведена оценка свободных параметров моделей выходов легких ядер. Оптимальные значения коэффициентов 1\ - для коалесцентной модели и Т0, N0 и Р - для модели движущегося источника для указанных реакций представлены в табл. 4.

Сравнительный анализ факторов согласия по моделям, входящих в «Benchmark of Spallation Models», для реакций na'Pb(p,x) представлен в табл. 5. Также, в таблице выделены штриховкой рекомендованные экспертами МАГАТЭ модели (incMabla, incl4gemini, inc!4smm) и сплошным цветом - модели, рассчитанная совокупность факторов согласия которых удовлетворяет условию парето-оптимальности (CASCADE/ASF, СЕМ02, g4bic, phits4jam, CASCADeX 1.2). Разница указанных совокупностей, возможно, связана с оценкой экспертами эффективности моделей в условиях традиционного однокритериалыюго подхода, расчеты по которым были осуществлены в ходе реализации проекта.

§4.2 в рамках многофакторного анализа, предоставляемого средой IViS CASCADeX был проведен сравнительный анализ CASCADeX 1.2 с моделями расчета высокоэнергетических нуклон-нуклоных реакций, входящих в программный комплекс MCNPX. В работе было проведено сравнение факторов для различных пороговых энергий программы CASCADeX 1.2 (0, 50, 100, 150, 200 МэВ) с факторами, полученными в работах под руководством С.Н.М. Breeders (2006). В таблицах 6-7 представлены результаты данного сравнения.

Таблица 4. Оценка параметров моделей выходов легких ядер

Мишень Реакция Энергия протона, МэВ Модель движущегося источника Коалес цент-ная модель 7 ^-сс Ссылка

No Р То Ро

РЬ РЫр,а) 63 0,15 0,035 5,5 19,837 F.E. Bertrand and R.W. Pelle, Phys. Rev. С 8 (1973) 1045

Pb(p,3He) 63 0,011 0,068 8,5 170 19,804

РЬШ 63 0,1 0,06 9,5 160 10,024

Pb (p,d) 63 0,205 0,068 10,5 180 10,024

Au Au (р,а) 2500 7 0,015 7 19,458 A. Bubak et al., Phys. Rev. С 76 (2007)014618

Aü(p,3Hc) 2500 0,4 0,004 11 210 19,804

Au (p,t) 2500 2,2 0,02 13,5 10,024

Au (p,d) 2500 3,3 0,015 11,5 10,024

Ni Ni(p,o) 175 0,21 0,33 7 240 10,298 S.V. Fôrtsch et al., Phys. Rev. С 43 (1991)691

Ni(p,3He) 175 0,039 0,03 13 10,235

Nifoi) 175 0,06 0,028 9 105 5,307

Nifocf) 175 0,4 0,028 8 250 5,275

Та Ta (p,a) 1200 8,6 0,01 4,5 220 18,452 C.-M. Herbach et al., Nucl. Phys. A 765 (2006) 426

Ta(p,3He) 1200 0,17 0,021 12,5 18,418

Ta(p,i) 1200 1,7 0,02 7,5 170 9,337

T&(p,d) 1200 2,9 0,025 7,5 9,319

Bi Bi(p,3He) 62 0,0186 0,088 5,92 19,984 F.E. Bertrand and R.W. Pelle, Phys. Rev. С 8 (1973) 1045

Bife/) 62 0,11 0,07 6,5 150 10,2

Bi(p,d) 62 0,2 0,07 10,5 200 10,1

Fe F e(p,a) 62 0,39 0,02 4,5 240 9,656 F.E. Bertrand and R.W. Pelle, Phys. Rev. С 8 (1973) 1045

Fe(p,3He) 62 0,042 0,045 5,8 110 9,595

FefeO 62 0,06 0,045 5,4 120 4,989

F e[p,d) 62 0,45 0,055 6,5 260 4,958

Al Мір,a) 160 0,09 0,04 10 6,083 A. Cowley et al., Phys. Rev. С 54 (1996)778

Таблица 5. Факторы согласия для реакций "а'РЬ(р, х)

Интегральный расчет для реакции па'РЬ(р,х), количество экспериментальных точек 279

Факторы согласия Количество

Модели H D R F полученных расчетных точек

Cas4 6,17 0,69 0,91 5,14 276

CasASF .Л'ЛЛ 4,62 0,49 0,91 ; '2,57 278 '.■■i'i'

CASCADeX 1.2 5,82 0,71 0,46 10,98 261

СЕМ02 -j 4,84 '0,51 1,05 ' 2,44 279

СЕМ03 5,21 0,56 1,06 2,46 279

g4bert 14,80 1,02 1,40 4,00 266

g4bic 4,39 0,53 0,69 3,73 258

'л ''-'incftabla.

1 < mcMgemmi'" Млб?/ <>У:2'А% ■■'/, ,

' m'(;f4smm --,/, ■'¿9,'.57,;--- Vaw-y v

BertDres 7,37 0,72 1,15 2,59 279

IsabelaAbla 13,13 1,08 1,77 2,29 250

IsabelGemini 30,30 1,70 2,49 2,79 240

Isabelasmm 10,04 0,92 1,35 4,04 250

Phitsjqmd 42,86 2,23 2,26 6,43 270

phits4jam 5,63 0,54 0,93. 2,12 277

phits/bertini 6,75 0,61 1,16 2,08 277

Обработка всей совокупности данных для l84W на основе методов многокритериального анализа показала преимущества моделей Bertini/Dresner CASCADE, CASCADE/ASF, CASCADcX 1.2. Анализ параметра пороговой энергии показал наилучшее совпадение с экспериментом при пороговой энергии 0 и 50 МэВ для 184W. Более того, при условии, что все факторы одинаково значимы и в качестве агрегированного показателя эффективности используется их линейная свертка, CASCADE/ASF представляет собой наилучшую модель для расчета сечений реакции l8"W(р, х) при соответствующих энергиях налетающей частицы.

В случае с 59Со эффективными оказались модели, заложенные в программы Bertini/Dresner, INCL4/Dresner, CEM2k, CASCADeX 1.2 с пороговой энергией, равной 50 МэВ. При этом наилучшей моделью при условии, что все факторы равнозначны, стала CASCADeX 1.2. Необходимо отметить' что модели Bertini/Dresner и CASCADeX 1.2 для данных реакций всегда входят в совокупность^парето-эффективных решений. Результаты интегрального расчета Сои W для количества экспериментальных точек, равного 218, представлены в табл. 8.

Таблица 6. Факторы согласия для различных моделей для реакции¿н-|84\У

Энергия налетающего протона 0,8 ГэВ, количество экспериментальных точек 67, сгпоп= 1636 мб Энергия налетающего протона 1,6 ГэВ, количество экспериментальных точек 91, стпо„= 1687 мб Интегральный расчет, количество экспериментальных точек

Модель 158

Факторы согласия

Я £> Я Р Н й Я Г Н £> Л

ВегйП1/Т)ге5ПеГ 5,08 0,38 0,83 1,76 6,89 0,48 0.87 1,87 6,19 0,44 0,85 1.82

ВеЛт1/АВЬА 5.04 0,39 0,83 2,28 5,67 0,44 0,89 2,63 5,41 0,42 0,86 2,48

^АВЕЬ/ОгеБпег 5,05 0,37 0,78 2,13 5,45 0.44 0,8 2,6 5,28 0,41 0,79 2.40

¡БАВЕЬ/АВЬА 5,35 0,38 0,78 2,13 5,91 0.4 0.83 2,83 5,68 0,41 0,81 2,58

тОЛ/Огевпег 5,56 0,41 0,75 2,2 5,25 0,38 0,79 2,73 5,38 0,39 0,77 2,51

ШСЬ4/АВЬА 6,18 0,43 0,75 2,54 6,08 0,4 0,85 2,57 6,12 0,41 0,81 2,56

СЕМ2к 4,85 0,43 0,8 2,89 5,88 0,45 0,89 3,6 5,47 0,42 0,85 3,30

САБСАБЕ 4,72 0,39 0,78 1,65" 4,9 0,35 0,83 2,85" 4,82 _0,37 0,81 2,39е

САЗСАОЕ/АБР 4,34 0,33 0,86 1,57 4,51 0,33 0,87 1.69 4,44 0,33 0,87 1,64

СА5САОеХ-1.2 (200 МэВ) 5,07 0,39 0,73 2,47° 7,52 0,44 0,86 5,86'" 6,64 0,42 0,81 4,40е

СА8СЛОеХ-1.2 (150 МэВ) 4,96 0,40 0,72 1,99" 6,49 0,40 0,83 7,05* 5,92 0,40 0,79 4,81е

СА5САПеХ-1.2 (100 МэВ) 5,39 0,41 0,71 1,82" 5,15 0.37 0,80 4,37* 5,25 0,39 0,76 3,33е

СА8САХ)еХ-1.2 (50 МэВ) 5,64 0,43 0,69 1,65" 4,89 0,36 0,75 4,63' 5,20 0,39 0,73 3,40е

СА8САОеХ-1.2 (0 МэВ) 4,88 0,38 0,70 1,68" 4,93 0,35 0,79 3,27* 4,82 0,36 0,75 2,29е

"Количество точек (ЛО равно 58, N = 86, 'Л' = 144.

Таблица 7. Факторы согласия для различных моделей для реакции р+кСо

Модель Энергия налетаюшего протона 1,2 ГэВ, количество экспериментальных точек 20, ст„оп = 772 мб Энергия налетающего протона 1,6 ГэВ, количество экспериментальных точек 20, ап[ш = 773 мб Энергия налетающего протона 2,6 ГэВ, количество экспериментальных точек 20, о„0п = 770 мб Интегральный расчет, количество экспериментальных точек 60

Н D R F Н D R F Я D R F Я D R F

Bertini/Dresner 4,87 0,32 0,7 1,74 4,51 0,33 0,71 1,96 4,29 0,32 0,71 1,76 4,56 0,32 0,71 1,82

Bertini/ABLA 15,85 0,81 1,5 2,07 13,79 0,81 1,43 2,11 13,71 0,8 1,42 2,11 14,71 0,81 1,45 2,1

ISABEL/Dresner 4,58 0,28 0,89 1,58 5,66 0,33 1,01 1,65 5,78 0.34 1,15 1,55 5,37 0,32 1,02 1,59

ISABEL/ABLA 21,17 1,1 1,83 2,31 23,16 1.27 2,04 2,45 28.15 1,63 2.38 2,75 24,34 1,33 2,08 2,5

INCL4/Dresner 4,16 0,25 0,91 1,56 4,30 0,25 0,95 1,51 4,42 0,27 1,03 1,47 4,29 0,26 0,96 1,51

JNCL4/ABLA 20,35 1,02 1,78 2,21 20,05 1,13 1,87 2.30 26 1,45 2.23 2,55 22.3 1,2 1,96 2,35

CEM2k 6,52 0,41 0,93 1,78 5.82 0.37 0,84 1,78 5.23 0,36 0,8 1,86 5,88 0,38 0,86 1,81

CASCADE 12,79 0,6 1,12 2,52 11,80 0,59 U1 2,38 10,26 0,58 1,08 2,31 11,66 0,59 1,1 2,4

CASCADE/ASF 6,02 0,36 1.1 1.5 5.51 0.37 1,09 1,48 5.51 0,37 1,08 1,49 5,69 0,37 1,09 1,49

CASCADeX-1.2 (200 МэВ) 10,43 0,59 1,15 2,14 9,33 0,59 1,14 2,16 9,57 0,62 1,16 2,13 9,79 0,6 1,15 2,14

CASCADeX-1.2 (150 МэВ) 10,13 0,61 1,19 1,96 9,17 0,6 1,19 1,93 8,48 0,55 1,15 1,87 9,29 0,59 1,18 1,92

CASCADeX-1.2 (100 МэВ) 6,27 0,38 1,09 1,52 6,37 0,43 1,08 1,7 6,08 0,4 1,08 2,04 6,24 0,4 1,08 1,76

CASCADeX-J.2 (50 МэВ) 3,5 0,22 0,99 1,32 3,47 0,23 1,01 1,38 3,85 0,27 1,03 1,39 3,61 0,24 1,01 1,36

CASCADeX-1.2 (0 МэВ) 16,0 0,74 1,44 2 14,5 0,73 1,43 1,98 14,1 0,76 1,44 2 14,9 0,74 1,44 1,99

N)

Таблица 8. Интегральный расчет 5,Со и l84W

Модель Интегральный расчет 5,Co и l84W, количество экспериментальных точек 218

Н D R F

Bertini/Dresner 5.79 0,41 0,81 1,82

' ;;::•;> ^Шиш&ЩНй- г Ч ЧП 1 ' " 1

ISABEL/Dresner 5,30 0,39 0,85 2,19

1SABEL/ABLA 13,65 0,66 1,16 2,56

!NCL4/T)resner 5,10 0,35 0,82 2,25

1NCL4/ABLA 12,81 0,63 1,13 2,50

CEM2k 5,59 0,41 0,85 п 2,90

CASCADE (original) 7,51 0,43 0,90 2,39

,CASCADE-?ASF'> . . "r- '

' CkSCADEK-1 (ОГМэВЪ i 0>47

- - С A SCADEX-1 * "¿ЩйМ1 i" - W& '•£> h ^Го/За^:';- if

CASCADEX-1.2 (100 МэВ) 5,56 0,39 0,86 2,88

CASCADEX-1.2 (150 МэВ) 7,08 0,46 0,90 3,92

1 CASCADEX-1.2 (200 МэВ) 7,70 0,47 0,91 3,71

ЗАКЛЮЧЕНИЕ В СООТВЕТСТВИЕ С ЗАДАЧАМИ И ЦЕЛЯМИ

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

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

• Разработана методика оценки моделей высокоэнергетических ядерных реакций на основе теории принятия решения при многих критериях;

• Произведена оценка адекватности моделей высокоэнергетических ядерно-ядерных взаимодействий на основании разработанных программных средств и методик на наиболее важных для конструирования ЭЛЯУ изотопах;

• Обновленные секции оцененных библиотек высокоэнергетических данных HEPAD и HEAD 2011 для расчёта электро-ядерных установок;

• Произведена оценка внутренних параметров каскадно-испарительной модели CASCADEX, коалесцентной модели и модели движущегося источника для следующих реакций 58Ni(p, х), S6Fe(p, х), 55Мп(р, х)

РЬ(р,х) в диапазоне от 150 МэВ до 1 ГэВ, па,РЬ(р, х), 209Bi(p, х) в диапазоне от 50 МэВ до 3 ГэВ и Pb(p, a), Pb(p, 3Не), РЬ(р, d), Pb(p, t), Au(p, a), Au(p, He), Au(p, d), Au(p, t), Ni(p, a), Ni(p, 3He), Ni(p, t), Ni(p, d), Ta(p, a),

Та(р 3Не), Та(р, d), Та(р, t), Bi(p, a), Bi(p, ЗНе), Bi(p, t), Bi(p, d), Fe(p, a), Fe(p, He), Fe(p, t), Fc(p, d), Al(p, а) в диапазоне от 63 МэВ до 1.2 ГэВ;

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

СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ

Статьи в научных рецензируемых журналах

1. Станковский А.Ю., Конобеев А.Ю., Купцов И.С. Программа CASCADEX для расчета взаимодействий частиц и ядер промежуточных и высоких энергий с веществом.// Известия вузов. Ядерная энергетика, 2008. №4.-С. 65-75.

2. Коровин Ю.А., Станковский А.Ю., Андрианов A.A., Конобеев А.Ю., Купцов И.С. Справочно-информационный интерактивный комплекс подготовки и верификации ядерно-физических данных в высокоэнергетической области.// Известия вузов. Ядерная энергетика, 2009. № 3. - С. 66-76.

3. Андрианов A.A., Конобеев А.Ю., Коровин Ю.А., Купцов И.С., Станковский А.Ю. Усовершенствованный программный комплекс CASCADEX 1.2 для расчета реакций глубокого расщепления.// Известия вузов. Ядерная энергетика, 2011. № 2. - С. 5-16.

4. Andrianov A.A., Korovin Yu. A., Kuptsov ].S., Stankovskiy A.Yu. Interactive Information System for Preparation and Veri_cation of Nuclear Data in the High-energy Range.// Journal of the Korean Physical Society, Vol. 59, No. 2, August 2011,-P. 1096-1099.

Материалы конференций и тезисы докладов

1. Korovin Yu. Development of a new code to simulate radiation damage and gas accumulation in the structural materials of ADS / Yu. Korovin, I.S. Kuptsov, A. Konobeyev, A. Stankovskiy, A. Natalenko // Workshop on Accelerator Radiation Induced Activation (ARIA'08), Villigen, Switzerland, 13-17 October, 2008. - Proc. PSI ARIA 08. - P. 48-54.

2. Коровин, Ю.А. Сравнительный анализ каскадно-испарительных моделей, применяемых для расчета протонных активационных данных в энергетическом диапазоне до 1 ГэВ / Ю.А. Коровин, И.С. Купцов, A.A. Ната-ленко, С.А. Осыкин // Актуальные проблемы физики ядерных реакторов -эффективность, безопасность, нераспространение - Волга-2008: материалы XV Семинара - М.: МИФИ, 2008. - С.ЗЗ.

3. Конобеев А.Ю., Станковский АЛО. Коровин Ю.А., Купцов И.С. Модифицированная каскадно-испарительная модель для расчета взаимодействий частиц и ядер промежуточных и высоких энергий с веществом. / Материалы XV Семинара по проблемам физики реакторов «Волга-2008», Москва, 2-6 сентября 2008 г. - Москва: НИЯУ МИФИ, 2008. - С. 96-97.

4. Купцов И.С., Наталенко A.A., Коровин Ю.А., Станковский А.Ю. Интерактивная среда комплексной подготовки ядерных активационных данных в энергетическом диапазоне от 0.2 до 1 ГэВ. / Научная сессия НИЯУ МИФИ-2009, Москва, 25-31 января 2009 г. - Москва: НИЯУ МИФИ, 2009. -Т.1.-С. 157.

5. Коровин Ю.А., Станковский А.Ю., Купцов И.С. Усовершенствованная каскадно-испарительная модель на основе формализма Хаузера-Фешбаха для расчета активационных данных в высокоэнергетической области. / Научная сессия НИЯУ МИФИ-2009, Москва, 25-31 января 2009 г. -Москва: НИЯУ МИФИ, 2009. - Т. 1. - С. 157.

6. Андрианов A.A., Коровин Ю.А., Купцов И.С., Догов A.A. Сопоставление расчетов дважды дифференциальных сечений выхода легких ядер по коалесцентной модели и модели движущегося источника. / Безопасность АЭС и подготовка кадров: тез. докл. XI Международной конференции, Обнинск, 29 сентября - 2 октября 2009 г. - Обнинск: ИАТЭ, 2009. - Т. 1. -С. 117.

7. Андрианов A.A., Ачаковский О.И., Евдокимов А.Н., Коровин Ю.А., Купцов И.С. Интерактивная среда работы с комплексом программ расчета радиационных повреждений. Безопасность АЭС и подготовка кадров // XI Международная конференция, Обнинск, 29 сентября - 2 октября 2009 г. -Обнинск: ИАТЭ, 2009. - Т.1. - С. 96-97.

8. Купцов И.С., Андрианов A.A., Коровин Ю.А., Догов A.A. Сравнительный анализ расчетов дважды дифференциальных сечений выходов легких ядер в высокоэнергетическом диапазоне по различным моделям. / Научная сессия НИЯУ МИФИ-2010, Москва, 25-31 января 2010 г. - Москва: НИЯУ МИФИ, 2010. - Т.1.-С. 236.

9. Купцов И.С., Андрианов A.A. Опыт применения информационной интерактивной среды подготовки высокоэнергетических ядерных активационных данных. / Научная сессия НИЯУ МИФИ-2010, Москва, 25-31 января 2010 г. - Москва: НИЯУ МИФИ, 2010. - Т.1. - С. 237.

10. Купцов И.С., Андрианов A.A., Коровин Ю.А. Модификация кода CASCADEX для расчета дважды дифференциальных сечений выходов легких ядер в высокоэнергетической области / Научная сессия НИЯУ МИФИ-2010, Москва, 25-31 января 2010 г. - Москва: НИЯУ МИФИ, 2010. - Т.1. - С. 234.

11. Kuptsov LS., Andrianov A.A., Korovin Yu.A. Stankovskiy A.Yu. Interactive Information System for Preparation and Verification of Nuclear Data in the Hugh-Energy Range. / International Conference on Nuclear Data for Science and Technology Proceedings of Nuclear Data conference 2010 April 26-30, 2010 Jeju Island, Korea, ND 2010-1443.

12. Купцов И.С., Андрианов A.A., Догов A.A. Оценка параметров ускорительно управляемых систем для мишеней из Fe, Pb, W в диапазоне энергий 0.1-1.5 ГэВ. / Научная сессия НИЯУ МИФИ-2012, Москва, 25-31 января 2012 г. - Москва: НИЯУ МИФИ, 2012. - Т.1. - С. 207.

Авторские свидетельства

1. Система визуального моделирования для программы СА8САОеХ. Свидетельство о государственной регистрации программы для ЭВМ №20116112258, дата регистрации 18 января 2011 г.

2. Справочно-информационный интерактивный комплекс подготовки и верификации ядерно-физических данных в высокоэнергетической области. Свидетельство о государственной регистрации программы для ЭВМ №20116112259, дата регистрации 17 марта 2011 г.

Компьютерная верстка И.С. Купцов

ЛР№ 020713 от 27.04.1998

Подписано к печати 2. ? ■ О 2.2.012 г- Формат бумаги 60x84/16

Печать ризограф. Бумага МВ Заказ № 2 /3 Тираж 100 экз. Печ. л. 1,5 Цена договорная

Отдел множительной техники ИАТЭ 249035, г. Обнинск, Студгородок, 1

61 12-1/657

Федеральное государственное автономное образовательное учреждение высшего профессионального образования

ОБНИНСКИЙ ИНСТИТУТ АТОМНОЙ ЭНЕРГЕТИКИ -филиал Национального исследовательского ядерного университета

«МИФИ»

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

КУПЦОВ Илья Сергеевич

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

Специальность 05.14.03 -«Ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации».

ДИССЕРТАЦИЯ

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

Научный руководитель

доктор физико-математических наук,

профессор Ю.А. Коровин

Обнинск 2012

СОДЕРЖАНИЕ

ВВЕДЕНИЕ..................................................................................................................6

ГЛАВА I. Современное состояние в области создания интерактивных

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

1.1. Потребности в высокоэнергетических ядерно-физических данных: состояние и перспективы.......................................................................................17

1.2. Современные программные средства расчета высокоэнергетических

ядерно-физических данных...................................................................................20

1.3. Библиотеки экспериментальных и оцененных ядерных данных в области высоких энергий......................................................................................29

1.3.1. Оцененные ядерно-физические данные.............................................................29

1.3.2. Экспериментальные ядерно-физические данные............................................33

1.4. Интерактивные системы оценки и хранения ядерных данных...............36

1.4.1. Ядерно-информационный комплекс ЗАШБ.......................................................36

1.4.2. Справочно-информационный интерактивная система ядерно-физических данных............................................................................................................................38

1.5. Выводы к главе 1................................................................................................39

Глава II. Модифицированная каскадно-испарительная модель для

расчета взаимодействий частиц и ядер промежуточных и высоких энергий с веществом............................................................................................40

2.1. УСОВЕРШЕНСТВОВАННАЯ ВЕРСИЯ ПРОГРАММНОГО КОМПЛЕКСА

САБСАБеХ 1.2..........................................................................................................40

2.2. Основные физические предположения модели..........................................42

2.2.1. Стадия внутриядерного каскада......................................................................42

2.2.2. Испарение и деление...........................................................................................47

2.2.3. Примеры расчетов с использованием предлагаемой модели.........................51

2.3. Интегрирование СА8САЕ>еХ и ТАЬУБ 1.2....................................................55

2.4. Дважды дифференциальные сечения вылета легких ядер........................57

2.4.1. Выход легких ядер в программе СА8САОЕХ....................................................63

2.5. Образования кластеров...................................................................................64

2.6. Некоторые особенности программной реализации...................................65

2.7. Интерактивный визуальный интерфейс 1Уй СА8САБеХ.........................66

2.8. Выводы к II главе..............................................................................................70

Глава III. Справочно-информационный интерактивный комплекс подготовки и верификации ядерно-физических данных в высокоэнергетической области.......................................................................71

3.1. Обоснование выбора языка программирования........................................72

3.2. Функциональные возможности СИИК.........................................................74

3.2.1. Блок схема СИИК................................................................................................75

3.2.2. Модуль расчета данных.....................................................................................77

3.2.3. Модуль обработки данных.................................................................................80

3.2.4. Модуль графического представления результатов и подготовки файлов библиотек ядерных данных..........................................................................................80

3.2.5. Модуль статистического анализа...................................................................81

3.3. Описание модулей СИИК................................................................................82

3.3.1. Детализированная структурная схема СИИК................................................83

3.3.2. Блок управления...................................................................................................84

3.3.3. Блок хранения......................................................................................................85

3.3.4. Блок обработки данных.....................................................................................87

3.3.5. Модуль подключения данных из базы ЕХЕОЯ..................................................87

3.3.6. Модуль расчета...................................................................................................88

3.3.7. Модуль подключения произвольных данных (Апу-с1а1а менеджер)................91

3.3.8. Модуль обработки данных ЕИОР/В формата................................................92

3.3.9. Модуль статистики...........................................................................................94

3.3.10. Блок визуализации...........................................................................................100

3.4. Выводы кIII главе...........................................................................................101

Глава VI. Модернизация константной базы расчетов перспективных

ядерно-энергетических установок................................................................102

4.1. Бенчмаркинг каскадных испарительных моделей..................................102

4.1.1. Международный бенчмаркинг реакций глубокого расщепления.................102

4.1.2. Расчет сечений каналов реакции па'РЬ(р,х) в диапазоне от 0,1 до 1 ГэВ.... 102

4.2. Определение оптимальных параметров модели САЗСАБеХ 1.2..........105

4.3. Сравнительный анализ с моделями, входящими в МСКРХ...................108

4.4. Выводы к IV главе...........................................................................................111

ЗАКЛЮЧЕНИЕ.....................................................................................................112

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ.......................................114

ОБОЗНАЧЕНИЯ И СОКРАЩЕНИЯ

ADS Ускорительно управляемые системы

ЭЛЯУ Электроядерная установка

Бенчмарк (англ. benchmark) — контрольная задача, необходимая для определения сравнительных характеристик производительности компьютерной системы.

СИИС справочно-информационная интерактивная система

СИИК справочно-информационный интерактивный комплекс

ЯФД Ядерно-физические данные

ГНЦИТЭФ Государственный научный центр Институт экспериментальной и теоретической физики, г. Москва

QMD (КМД) Quantum-Molecular Dynamics (Квантово-Молекулярная Динамика)

ООП Объектно-ориентированный язык Программирования

CINDA Computer Index of Neutron Data - Компьютерный указатель нейтронных данных, содержит библиографические ссылки на измерения, расчеты, обзоры и оценки нейтронных сечений и других микроскопических нейтронных данных

BROND-3 Russian neutron evaluated data files - Библиотека оцененных нейтронных данных (Россия)

МАГАТЭ Международное агентство по атомной энергии

CENDL-3 Chinese Evaluated Nuclear Data Library - Китайская библиотека оцененных ядерных данных

ENDF/B-VI Evaluated Nuclear Reaction Cross Section Libraries - Библиотеки оцененных ядерных данных (США)

ENDF-6 Evaluated Nuclear Data Format - формат представления оцененных ядерных данных

EXFOR Experimental Nuclear Reaction Data File; EXchange FORmat - Международная библиотека экспериментальных ядерных данных;

Формат представления экспериментальных ядерных данных

GEM Generalized evaporation model — Обобщенная испарительная мо-

дель

INC, В Ж Intranuclear Cascade - Внутриядерный каскад

JENDL-3.3 Japanese Evaluated Nuclear Data Library - Японская библиотека оцененных ядерных данных

MCNPX Monte-Carlo N-Particle - Программа для моделирования транспор-

та частиц методом Монте-Карло

NEA Nuclear Energy Agency - Агентство по ядерной энергии (АЯЭ)

IEAF-2001, Intermediate Energy Activation Files - Библиотеки оцененных ак-

IEAF-2005 тивационных ядерных данных

ИАТЭ Обнинский государственный технический университет атомной

энергетики

РОСФОНД РОСсийская библиотека Файлов Оцененных Нейтронных Данных

no программное обеспечение

MBK Модель внутриядерного каскада

MSDN Microsoft Developer Network

BNL Brookhaven National Laboratory - Брукхейвенская Национальная

Лаборатория

КПЖСР Каскадный Подкритический Жидкосолевой Реактор

MA Минорные актиниды

РАО Радиоактивные отходы

ЛУП Линейный ускоритель протонов

БР Быстрый реактор

УУС Ускорительноуправляемая система

ЛВР Легководный реактор

ВТГР Высокотемпературный газовый реактор

ВВЕДЕНИЕ

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

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

Концептуальные исследования свидетельствуют о перспективности электроядерного способа производства нейтронов для трансмутации долгоживущих радиоактивных отходов (ДРАО). Этот способ заключается в использовании подкртического реактора (называемого бланкетом) с внешним источником нейтронов в виде нейтронопроизводящей мишени и ускорителя протонов. В российской технической литературе этот тип установок называют электроядерными установками (ЭЛЯУ), в зарубежной - Accelerator Driven System (ADS).

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

W "ГЛ U

задач по созданию неитронопроизводящеи мишени. В этой связи специалистами данной предметной области был сформулирован перечень следующих НИОКР:

• Развитие методов расчета сечений реакций глубокого расщепления в области промежуточных энергий протонов 20-2000 МэВ;

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

• Развитие и создание компьютерных программ расчета полной и парциальной активноти, остаточного энерговыделения в материалах мишени;

• Создание тестированных библиотек оценненых ядерных данных о взаимодействии протонов и нейтронов с энергиями до 3 ГэВ с материалами мишени, ускорителя и бланкета.

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

Объектом исследования являются методы и программные комплексы формирования систем константного обеспечения расчетов активации материалов и выходов газообразных продуктов в электро-ядерных установках материалы в энергетическом диапазоне от 50 МэВ до 1 ГэВ.

Цель диссертационного исследования заключается в разработке методов и программных комплексов формирования систем константного обеспечения расчетов активации материалов и образования легких газообразных продуктов в перспективных ЯЭУ на основе методов теории принятия решений при многих критериях, усовершенствованию расчетного кода CASCADEX, включающего развитие методов расчета сечений реакций глубокого расщепления, эмиссии кластеров и выходов легких ядер в окне и мишени электро-ядерных установок (ЭЛЯУ), обновлению и расширению файлов библиотек оцененных ядерных данных HEPAD, HEAD.

Цель диссертации определяет постановку и решение следующих задач:

1. Изучение особенностей, областей применения каскадно-испарительной модели СА8САЭЕХ;

2. Реализация процедуры расчета выходов легких ядер на базе модели СА8САБЕХ;

3. Разработка алгоритма оценки высокоэнергетических ядерных данных на основании многофакторного анализа;

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

5. Оценка моделей высокоэнергетических ядерно-ядерных взаимодействий на основании разработанных программных средств и методик.

На защиту выносятся следующие положения:

1. Методы оценки активационных ядерных данных в высокоэнергетической области с использованием методов теории принятия решений при многих критериях;

2. Программный комплекс СА8САБЕХ 1.2, включающий в себя учет выходов легких ядер, определяющих накопление газообразных продуктов в мишени ЭЛЯУ и продуктов реакций глубокого расщепления;

3. Интерактивные программные комплексы подготовки, верификции и вали-дации баз активационных ядерных данных в области энергий от 0,1 до 10 ГэВ;

4. Результаты сравнительного анализа моделей ядерно-ядерных взаимодействий в области энергий от 50 МэВ до 1 ГэВ и данных натурного эксперимента для основных конструкционных материалов ЭЛЯУ (56Ре,59Со, 184\У, 209В1 и паИ).

Научная новизна данного диссертационного исследования заключается в том,

что в работе:

1. Впервые разработан интерактивный программный комплекс для подготовки, верификации и валидации активационных ядерных данных в области энергий от 50 МэВ до 10 ГэВ на основе каскадно-испаритаельных моделей, необходимых для расчетного обоснования проектов ЭЛЯУ;

2. Впервые предложен подход к оценке активационных ядерных данных в области энергий от 50 МэВ до 10 ГэВ на основе методов теории принятия решений при многих критериях, позволяющий формировать самосогласованный набор систем константного обеспечения расчетов ЭЛЯУ;

3. Впервые предложена модифицированная версия ПК CASCADEX 1.2, включающая в себя учет выходов легких ядер, определяющих накопление газообразных продуктов в мишени ЭЛЯУ и продуктов реакций глубокого расщепления на основе модифицированного подхода Хаузера-Фешбаха, реализованного в TALYS 1.2;

4. Впервые представлены результаты сравнительного анализа моделей ядерно-ядерных взаимодействий в области энергий от 50 МэВ до 1 ГэВ и данных натурного эксперимента для основных конструкционных материалов ЭЛЯУ (56Fe, 59Со, 184W, 209Bi и natPb).

Практическая ценность работы состоит в следующем:

1. Проведена верификация и валидация программного комплекса CASCADEX на известных экспериментах и бенчмарках. Даны рекомендации по свободным параметрам данной модели для изотопов необходимых для константного обеспечения ЭЛЯУ;

2. Подготовлены секции в библиотеки оцененных ядерных данных для ряда изотопов необходимых для расчета характеристик при проектировании и оптимизации эксплуатационных режимов ЭЛЯУ;

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

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

Научные положения, выводы и рекомендации, сформулированные в диссертации, обоснованы поскольку были апробированы в международном проекте «Benchmark of Spallation Models» проводимым МАГАТЭ, в работе использовались современные экспериментальные данные и теоретические представления о механизмах ядерных реакций в рассматриваемом энергетическом диапазоне, а также современные программы расчета ядерно-ядерных взаимодействий в высокоэнергетической области.

Личный вклад соискателя в представленную работу заключается:

1. В разработке и реализации справочно-информационного интерактивного комплекса, подготовке библиотек ядерно-физических данных;

2. Реализации модификаций модели CASCADEX, а именно, в создании интерактивной среды работы с кодом, в осуществлении процедуры расчета выходов легких ядер по коалесцентной модели и модели движущегося источника, в соединении CASCADEX и Talys 1.2;

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

4. Оценке моделей на основании разработанных программных средств и методик;

5. В обновлении и расширении файлов библиотек оцененных ядерных данных HEPAD, HEAD.

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

1. Безопасность АЭС и подготовка кадров. X Международная конференция, Обнинск, 1-4 октября 2007 г.

2. XV семинар по проблемам физики реакторов («Волга-2008»). Актуальные проблемы физики ядерных реакторов - эффективность, безопасность, нераспространение. Москва, 2-6 сентября 2008 года.

3. Workshop on Accelerator Radiation Induced Activation (ARIA'08), Villigen, Switzerland, October 13 - 17, 2008.

4. Безопасность АЭС и подготовка кадров. XI Международная конференция, Обнинск, 29 сентября - 2 октября 2009 г.

5. Научная сессия НИЯУ МИФИ-2009, Москва, 25-30 января 2009

6. ND-2010 - International Conference on Nuclear Data for Science and Technology. Jeju Island, Korea, 26 - 30 April, 2010.

7. Научная сессия МИФИ - 2010, Москва, 25 - 31 января, 2010.

8. Научная сессия МИФИ - 2011, Москва, 02 - 05 февраля, 2011.

9. Научная сессия МИФИ - 2012, Москва, 30 января - 4 февраля, 2012.

Результаты исследования представлены в статьях: по теме диссертации опубликовано 16 работ, в том числе 4 статьи в ведущих рецензируемых научных журнал�