CC BY
37УДК 621.039.578:629.7
обзор программных комплексов и результатов расчетно-экспериментальных исследований и оптимизации характеристик систем с термоэмиссионным преобразованием энергии
© 2020 г. Алексеев П.А., Кротов А.д., Кухарчук О.ф., Пышко А.П., ярыгин в.и.
АО «Государственный научный центр Российской Федерации -Физико-энергетический институт имени А.И. Лейпунского» (АО «ГНЦ РФ-ФЭИ») Пл. Бондаренко, 1, г. Обнинск, Калужская обл., Российская Федерация, 249033,
e-mail: postbox@ippe.ru
Основной целью настоящей статьи является краткий обзор мультифизичных задач и методов их решения в области космических ядерных энергетических установок, решаемых и разрабатываемых специалистами АО «ГНЦ РФ-ФЭИ». Физико-техническое обоснование разработки и эксплуатации космических ядерных энергетических установок с прямым преобразованием тепловой энергии в электричество включает в себя электротеплофизический расчет характеристик термоэмиссионного электрогенерирующего элемента и электрогенерирующего канала в целом, нейтронно-физические расчеты и оптимизацию структуры активной зоны и радиационной защиты.
Для определения электротепло физических характеристик электрогенерирующего канала применяется его трехмерный детальный расчет (программный комплекс COMSOL-ЭГК). Для нейтронно-физического расчета характеристик активной зоны и радиационной защиты применяются единые трехмерные расчетные программные комплексы, реализующие метод Монте-Карло (в частности, MMKFK-2, MCNP), позволяющие учитывать гетерогенность и все особенности геометрии. Перечисленные в обзоре расчетные методики и комплексы системно и комплексно применяются в ГНЦ РФ-ФЭИ при обосновании характеристик проектируемых в настоящее время космических ядерных энергетических установок.
Ключевые слова: термоэмиссионный реактор-преобразователь, радиационная защита, электрогенерирующий канал, MCNP, COMSOL-ЭГК.
DOI 10.33950/spacetech-2308-7625-2020-3-114-128
review of software solutions and results of computational and experimental studies and performance optimization of systems with thermionic energy conversion
Alekseev p.A., Krotov A.D., Kukharchuk O.F., pyshko A.p., Yarygin V.I.
Joint Stock Company State Scientific Centre of the Russian Federation -Institute for Physics and Power Engineering named after A.I. Leypunsky (JSC SSC RF-IPPE)
1 Bondarenko sq, Obninsk, Kaluga region, 249033, Russian Federation, e-mail:postbox@ippe.ru
The main objective of this paper is to summarize multiphysical tasks and methods for their solution in the area of space nuclear power systems being solved and developed by specialists of JSC SSC RF-IPPE. The physical and technical justification of the development and operation of space nuclear power systems with direct conversion of thermal energy into electricity includes an electrothermophysical calculation
of characteristics of a thermionic electrical power generating cell and electrical power generating channel as a whole, neutronic physical calculations and structural optimization of the core and radiation shielding.
To determine the electrothermophysical performance of the electrical power generating channel, its 3D detailed computation (software COMSOL-EGK) is used. For the neutronic calculation of the core and radiation shielding characteristics, unified 3D computational software systems are used that implement the Monte Carlo method (in particular, MMKFK-2, MCNP), making it possible to take into account the heterogeneity and all geometry features. The computational procedures and solutions listed in the review are applied in a systematic and integrated manner at JSC SSC RF-IPPE in justification of performance of space nuclear power systems being currently designed.
Key words: thermionic re actor-converter, radiation shielding, electrical power generating channel, MCNP, COMSOL-EGK.
АЛЕКСЕЕВ П.А. КРОТОВ А.Д. КУХАРЧУК О.Ф. ПЫШКО А.П. ЯРЫГИН В.И.
АЛЕКСЕЕВ Павел Александрович — кандидат технических наук, старший научный сотрудник АО «ГНЦ РФ-ФЭИ», e-mail: palekseev@ippe.ru
ALEKSEEV Pavel Aleksandrovich — Candidate of Science (Engineering), Senior research scientist at JSC SSC RF-IPPE, e-mail: palekseev@ippe.ru
КРОТОВ Алексей Дмитриевич — старший научный сотрудник АО «ГНЦ РФ-ФЭИ», e-mail: akrotov@ippe.ru
KROTOV Aleksey Dmitrievich — Senior research scientist at JSC SSC RF-IPPE, e-mail: akrotov@ippe.ru
КУХАРЧУК Олег Филаретович — доктор физико-математических наук, заместитель генерального директора АО «ГНЦ РФ-ФЭИ», e-mail: kuh@ippe.ru
KUKHARCHUK Oleg Filaretovich — Doctor of Science (Physics and Mathematics), Deputy General Director at JSC SSC RF-IPPE, e-mail: kuh@ippe.ru
ПЫШКО Александр Павлович — кандидат физико-математических наук, начальник отдела АО «ГНЦ РФ-ФЭИ», e-mail: pyshko@ippe.ru PYSHKO Aleksandr Pavlovich — Candidate of Science (Physics and Mathematics), Head of Department at JSC SSC RF-IPPE, e-mail: pyshko@ippe.ru
ЯРЫГИН Валерий Иванович — доктор технических наук, профессор, главный научный сотрудник АО «ГНЦ РФ-ФЭИ», e-mail: ecs-yar@ippe.ru
YARYGIN Valery Ivanovich — Doctor of Science (Engineering), Professor, Chief research scientist at JSC SSC RF-IPPE, e-mail: ecs-yar@ippe.ru
введение
АО «ГНЦ РФ-ФЭИ» обладает уникальными компетенциями в области физико-технического обоснования разработки и эксплуатации космических ядерных энергетических установок (КЯЭУ) с прямым преобразованием тепловой энергии в электричество, включая методы электротеплофизического расчета характеристик термоэмиссионного электрогенерирующего элемента (ЭГЭ) и электрогенерирующего канала (ЭГК) в целом, нейтронно-физические расчеты и оптимизацию структуры активной зоны и радиационной защиты. Физические основы, заложенные при обосновании характеристик КЯЭУ первого поколения В.Я. Пупко, Ю.С. Юрьевым, В.А. Малыхом, А.А. Дубининым и другими, в настоящее время подтверждаются и углубляются разработкой и применением современных расчетных комплексов и новейших математических методов оптимизации и поиска решений. Так, вместо одномерного упрощенного расчета характеристик ЭГЭ/ЭГК (комплекс TFEDM) применяется трехмерный детальный расчет всего ЭГК (COMSOL-ЭГК), вместо гомогенного расчета активной зоны и отдельного, не учитывающего проникающие элементы расчета радиационной защиты, применяются единые трехмерные расчетные комплексы, позволяющие учитывать гетерогенность и все особенности геометрии активной зоны и радиационной защиты, реализующие метод Монте-Карло (ММК) (в частности, программные комплексы MMKFK-2, MCNP). Совместное применение перечисленных комплексов при расчетном обосновании характеристик проектируемых КЯЭУ позволяет наиболее полно исследовать пространство возможных решений и находить лучшие из них.
Для определения выходной электрической мощности КЯЭУ необходимо провести электротеплофизический расчет характеристик ЭГЭ/ЭГК. Одними из основных исходных данных для этого расчета являются значения тепловыделения в топливе. Такие исходные данные предоставляет нейтронно-физи-ческий расчет характеристик активной зоны, который, в первую очередь, определяет критические параметры активной
зоны, запас реактивности реактора на кампанию и т. д. Зная габариты активной зоны, ее тепловую мощность и расположение органов регулирования, можно проводить расчет радиационной обстановки вокруг реактора и оптимизационные исследования в обоснование массогабаритных характеристик теневой радиационной защиты. В конечном итоге масса и габариты радиационной защиты, пространственное расположение отсеков космического аппарата с оборудованием определяют облик ядерно-энергетической установки, ее конкурентоспособность при сравнении характеристик КЯЭУ с другими типами ядерных энергоустановок с различными методами преобразования тепловой энергии в электрическую, в т. ч. и солнечными энергетическим установками.
методы расчета тепло- и электрофизических характеристик термоэмиссионных электрогенерирующих систем
С 1960-х гг. начали развиваться методы расчета тепло- и электрофизических характеристик ЭГК и термоэмиссионного реактора-преобразователя (ТРП) в целом. Спецификой этих методов является то, что приходится совместно решать тепло- и электрофизические задачи, так как термоэмиссионный ЭГК является одновременно ядерным твэлом и электрогенератором.
Конечной задачей тепло- и электрофизического расчета ТРП является получение его выходных вольтамперных характеристик (ВАХ), т. е. зависимости электрического тока от напряжения при заданной тепловой мощности реактора. Знание ВАХ позволяет определить электрическое напряжение и ток, выходную электрическую мощность и эффективность (КПД) преобразования тепловой энергии ТРП в рабочем режиме. Одновременно необходимо определить значение максимального уровня температуры эмиттеров ЭГК в активной зоне реактора-преобразователя, который не должен превышать предельно допустимого значения для данного материала эмиттерной оболочки. Однако, для того чтобы получить эти важнейшие выходные характеристики ТРП, состоящего из сотен ЭГК, соединенных в сложную электрическую цепь и находящихся в активной зоне реактора
при различных уровнях плотности тепловыделения, необходимо провести целую серию теплоэлектрофизических расчетов единичного ЭГК при различных значениях тепловыделения, давления пара цезия в межэлектродном зазоре (МЭЗ) и ряда других внешних параметров. В настоящее время это стало возможным благодаря развитию информационных технологий и численных методов анализа. Широкое распространение получили интерактивные программы графического представления результатов расчета, основанные на решении краевых задач математической физики с помощью метода конечных элементов, такие как АЫЗУЗ, СОЫБОЬ, ЗоЫЦ^откв и др. [1].
Основные уравнения математической модели тепло- и электрофизических процессов, протекающих в единичном ЭГЭ термоэмиссионного ЭГК, были получены Ю.С. Юрьевым в ходе работ по созданию термоэмиссионной КЯЭУ «ТОПАЗ» [2]. Дальнейшее развитие методов численного решения системы нелинейных уравнений ЭГЭ проводилось В.А. Ружниковым [3-5], А.А. Шиманским [6], В.В. Синявским, Ю.В. Бабушкиным, В.П. Зиминым [7-10], В.А. Линником [11], Е.Г. Виноградовым [12] и др.
Первоначально из-за отсутствия достаточно мощных вычислительных средств разрабатываемые методы расчета тепло- и электрофизических характеристик ЭГК основывались на аналитическом решении системы уравнений ЭГЭ, что приводило к необходимости упрощения математической модели. В частности, приходилось вводить допущения о линейности локальной ВАХ термоэмиссионного преобразователя (ТЭП), постоянстве температуры коллектора и плотности эмиссионного тока по длине ЭГЭ [2]. При этих допущениях уравнение теплопроводности для коллектора вообще исключалось из математической модели, а уравнение теплопроводности для эмиттера при условии линеаризации члена, описывающего теплопередачу излучением, и уравнение для межэлектродного напряжения становились линейными и допускали аналитическое решение. Позднее для решения нелинейного уравнения теплопроводности для эмиттера были применены вариационные методы и метод Галеркина [11], позволяющие
более точно рассчитать распределение температуры по длине эмиттера.
Последующий прогресс в развитии расчетных методов исследования ВАХ ЭГК был связан с разработкой численных методов решения системы нелинейных дифференциальных уравнений ЭГЭ [3-6]. Тем не менее, используемая при расчетах математическая модель тепловых и электрических процессов ЭГЭ, даже при некоторых упрощающих допущениях (азимутальная симметрия, пренебрежение радиальным распределением параметров ввиду «тонкости» электродов и др.), сводится к системе одномерных нелинейных дифференциальных уравнений второго порядка, описывающих только распределение температуры эмиттера, коллектора и разности потенциалов между электродами и их распределение по длине ЭГК. Нелинейность этих уравнений обусловлена наличием теплопередачи между электродами, излучением и эмиссионным током электронов. Во все уравнения входит зависимость эмиссионного тока от температуры эмиттера и коллектора, межэлектродного напряжения и ряда других параметров, которая и замыкает эти уравнения в единую систему. Расчет этой зависимости (локальной ВАХ ТЭП) является сложной задачей, требующей рассмотрения элементарных процессов в низкотемпературной плазме МЭЗ и процессов, протекающих на электродах ТЭП. До настоящего времени эта задача в полном объеме не решена, поэтому экспериментальное определение локальных ВАХ в ходе лабораторных стендовых испытаний ТЭП/ЭГЭ с электронагревом является актуальной задачей.
Расчет ВАХ многоэлементного ЭГК обычно сводится к последовательному расчету ВАХ входящих в него ЭГЭ и их суммированию при заданной величине протекающего электрического тока, одинакового для всех ЭГЭ вследствие их последовательного соединения. Основным блоком расчетной методики является решение системы тепловых и электрических уравнений для ЭГЭ.
Известные в литературе расчетные методики в основном позволяют выполнять оценочные и инженерные расчеты [13]. Этим методикам присущи следующие характерные особенности, существенно снижающие точность расчета:
• расчет ЭГК сводится к последовательному расчету входящих в его состав ЭГЭ;
• при построении математических моделей электротеплофизических процессов в ЭГЭ используется ряд упрощающих допущений, сильно влияющих на соответствие используемой модели реальным процессам;
• используемые математические модели электротеплофизических процессов в ЭГЭ в большинстве методик являются одномерными;
• для используемых математических моделей электротеплофизических процессов в ЭГЭ характерен некоторый произвол в выборе граничных условий.
Модернизация существующих или создание новых расчетных методик должны идти по пути частичного или полного устранения указанных выше упрощений. Сформулируем минимальный перечень требований к современным расчетным методикам:
• переход от одномерного к трехмерному моделированию;
• детальный учет зависимости свойств конструкционных материалов и сред от температуры;
• переход от математического моделирования процессов в ЭГЭ к ЭГК в целом;
• возможность проведения расчетов ЭГЭ/ЭГК сложной геометрии;
• возможность «сквозного» расчета для ЭГК, а в перспективе — ТРП и ЭУ в целом.
Для решения этой задачи в ГНЦ РФ-ФЭИ была выполнена модификация стандартного пакета конечноэлементного анализа СОМБОЕ [14].
Разработанный на базе СОМБОЬ программный комплекс СОМБОЬ-ЭГК, оптимизированный под задачи расчета термоэмиссионных ЭГК, предоставляет пользователям полный спектр инструментов для решения поставленной задачи:
• построение трехмерной геометрической модели;
• описание физических процессов на языке математической физики;
• построение конечноэлементной расчетной сетки геометрической модели;
• расчет и постобработка результатов расчета.
Программный комплекс позволяет моделировать физические процессы электромагнетизма и сопряженного теплообмена
в трехмерной геометрии как в пределах одного ЭГЭ, так и внутри коммутационного пространства ЭГК. Важным достоинством описываемого комплекса является возможность прямого использования чертежей различных ЭГК, выполненных, например, в среде автоматизированного проектирования AutoCAD.
COMSOL-ЭГК реализует методику расчета тепло- и электрофизических характеристик ЭГК и термоэмиссионной электрогенерирующей системы (ТЭС) в целом на основе трехмерной математической модели с возможностью использования в качестве исходных данных атласа экспериментальных изотермических и эквипотенциальных ВАХ ТЭП в широком диапазоне изменения таких параметров, как температура электродов, давление пара цезия и ряда других. Методика предназначена для обоснования проектных решений перспективных термоэмиссионных КЯЭУ и наземных ТЭС с ядерным или неядерным нагревом эмиттеров.
В численный алгоритм комплекса внедрены решения уравнений для электрического потенциала, генерации тока, переноса энергии тепловым излучением и электронами эмиссии, джоулева тепловыделения в электропроводящих материалах. С помощью программного комплекса COMSOL-ЭГК можно моделировать практически все физические процессы, которые описываются уравнениями в частных производных. Далее задачи решаются методом конечных элементов. Имеются также дополнительные средства для моделирования: средства для геометрических построений, генераторы сетки, инструменты постобработки и т. п.
В качестве примера для трехмерного численного моделирования с помощью программного комплекса COMSOL-ЭГК была использована базовая конструкция девятиэлементного ЭГК для термоэмиссионных КЯЭУ второго поколения (так называемый унифицированный ЭГК [15]).
Расчетная модель ЭГК представлена на рис. 1. Эта модель многоэлементного ЭГК имеет 2D-осесимметричное приближение. На рисунке представлен общий вид ЭГЭ/ЭГК, используемый в расчетной модели, а также более детальный вид ЭГЭ в области его межэлектродной коммутации.
Разбиение модели ЭГК на конечные элементы производилось программой СОЫБОЬ-ЭГК в полуавтоматическом режиме с выбором типа разбиения (квадратная или треугольная форма расчетных ячеек), а также минимального и максимального размера ячеек, который варьировался для различных сред расчетной сетки.
Рис. 1. Расчетная модель девятиэлементного унифицированного ЭГК, разработанная в программной среде СОМ$ОЬ-ЭГК [14]: 1 — цезиевая среда; 2 — ядерное топливо; 3 — эмиттер; 4 — коллектор; 5 — коллекторный пакет; 6 — газоотводящее устройство; 7 — межэлектродная коммутационная перемычка
О а)
Рис. 2. Распределение температурного поля: а — многоэлементного ЭГК — коллекторная концевая серия при величине протекающего электрического тока 100 А (электродная пара Р1-ВХ2У [12]); б — многоэлементного ЭГК при величине протекающего электрического тока 100 А (электродная пара Р1-ВХ2У [12])
Основным результатом расчетов элек-тротеплофизических характеристик ЭГК с помощью программного комплекса СОЫБОЬ-ЭГК являются стационарные двумерные распределения температуры, потенциала и плотности генерируемого электрического тока. Программа также позволяет рассчитывать распределение тепловых потоков в любой рассматриваемой области задачи, например, тепловой поток через МЭЗ. По найденным распределениям определяются такие выходные характеристики ЭГК, как изомощностная ВАХ, электрическая мощность, КПД преобразователя и ряд других.
На рис. 2 представлены распределения температурного поля по конструкционным элементам отдельного термоэмиссионного ЭГЭ и девятиэлементного ЭГК в целом. На рис. 3, а показаны расчетные изомощностные ВАХ ЭГК, а на рис. 3, б — зависимости выходной электрической мощности ЭГК от тока для различных материалов электродных пар.
Использование программного комплекса СОЫБОЬ-ЭГК для расчетов электротеплофи-зических характеристик многоэлементного термоэмиссионного ЭГК в трехмерной постановке показало достаточную гибкость и эффективность этой методики, позволяющей во всей полноте учесть реальную конструкцию ЭГК и разнообразие физических свойств материалов его конструктивных элементов.
Моделирование тепловых и электрических характеристик много-
элементного ЭГК с помощью программного комплекса СОЫБОЬ-ЭГК открывает путь к полноценному исследованию
б)
термомеханических напряженно-деформированных состояний конструкционных элементов и сред ЭГК в трехмерной постановке, что является одной из приоритетных задач при обосновании ресурсных характеристик ЭГК и термоэмиссионной ЯЭУ в целом. Использование программного комплекса СОМБОЬ-ЭГК может сыграть важную роль в разработке методов прогнозирования ресурса термоэмиссионных ЭГК по результатам его петлевых реакторных испытаний на укороченной временной базе.
щ
с
о и1
а
н
о
Л
240 220 200 ISO-160 140 120 100 8060 4020
— W(::«) -W(iii)
— Pt- ВХ2У
* ЧЧ
VS
ч\
S
Ч *
\
\ х
V \
0 123 4 56789 Выходное напряжение ЭГК, В
а)
10
Электрическийг б)
Рис. 3. Изомощностные эквипотанциальные ВАХ ЭГК для различных материалов электродных пар (а) и зависимость выходной электрической мощности ЭГК от тока для различных материалов электродных пар (б) [14]
методы нейтронно-физического расчета и оптимизации активной зоны термоэмиссионного реактора-преобразователя
В связи со спецификой применения термоэмиссионных ЯЭУ на них накладываются ограничения, прежде всего — по массогабаритным параметрам.
Следовательно, задача нейтронно-физического расчета — выбрать такую конфигурацию активной зоны, которая удовлетворяла бы поставленным условиям по тепловой мощности, запасу реактивности, обеспечению безопасной эксплуатации, а с другой стороны, имела бы минимально возможные габариты и массу. Такая задача является, по сути, оптимизационной и может быть решена с применением соответствующих математических методов.
Расчет нейтронно-физических характеристик ТРП изначально проводился по одномерным программам в рамках теории возмущений, вариационного исчисления и др. [16]. В гомогенной модели невозможно было учесть многие особенности малогабаритных реакторов. В частности, не было возможности точно моделировать располагающиеся в отражателе поворотные цилиндры регулирования реактивности. Расчеты проверялись экспериментально на физических стендах [17]. Основываясь на данных экспериментов, для оптимизации эффективного коэффициента размножения нейтронов (Кэфф) определялась функция эффективности топлива, которая показывала, насколько изме-
няется K
эфф
реактора при перераспре-
делении топлива по объему активной зоны. Функция эффективности топлива определялась дифференцированием кривой зависимости Кэфф от положения топлива в активной зоне [18].
С развитием вычислительных машин и методов расчета ядерных реакторов более широко стали применяться алгоритмы, основанные на использовании ММК. Применяемые в ГНЦ РФ-ФЭИ современные расчетные программные комплексы, основанные на ММК, такие как ММКГК-2 [19] и МСЫР [20], позволяют детально воспроизводить компоновку активной зоны термоэмиссионного реактора, в состав которого также входят торцевой и боковой отражатели, поворотные органы регулирования, коллекторы теплоносителя. При моделировании учитывается изотопный состав материалов конструкционных элементов, положение органов регулирования, наполненность каналов теплоносителем и др. Пример расчетной модели реактора «ТОПАЗ», построенной в программном комплексе МСЫР, показан на рис. 4.
Рис. 4. Расчетная модель ТРП «ТОПАЗ», построенная в программном комплексе МСМР
Также достаточно детально учитываются особенности конструкции ЭГК. В состав расчетной модели ЭГК входят ядерное топливо, эмиттер, слои коллекторного пакета, система отвода газообразных продуктов деления, коммутационные перемычки между отдельными ЭГЭ, торцевые отражатели нейтронов и т. д. В некоторых случаях для ускорения счета коммутационные перемычки представляют в виде гомогенной смеси. На рис. 5 представлен фрагмент характерной расчетной модели такого ЭГК.
Рис. 5. Фрагмент расчетной модели ЭГК, построенной в программном комплексе МСМР: 1 — коллекторный пакет; 2 — топливо; 3 — эмиттер; 4 — система отвода газоподобных продуктов деления; 5 - коммутация; 6 — тракт теплоносителя
Такое детальное представление ЭГК необходимо, в т. ч., и для получения распределения тепловыделения по высоте и диаметру твэл, которое необходимо для электротеплофизического расчета ЭГК (см. выше).
Расчет реактора подразумевает изменение множества параметров, поэтому подбор компоновки активной зоны, толщины отражателя и др. параметров для удовлетворения поставленным условиям является сложной задачей.
Наряду с развитием вычислительных методов появились новые методики оптимизации. Одним из таких методов является генетический алгоритм (ГА) —
адаптивный метод поиска, основанный на механизмах и методах эволюции и генетических процессах. Подбор вариантов решений проходит с помощью математических формул, описывающих такие природные механизмы, как скрещивание и мутации, на которые, в свою очередь, влияет окружающая среда [21].
Важной особенностью ГА является то, что он не требует явного задания оптимизируемой функции и способен работать напрямую с численными решениями уравнения переноса, получаемыми из расчетной программы (например, из расчетов ММК).
Удачным подходом стало соединение комплекса для расчета реактора и оптимизационной методики, оперирующей непосредственно размерами конструктивных элементов активной зоны и реакторными функционалами, получаемыми при расчете реактора.
На основе данных, заданных пользователем, разработанная программа подготавливает набор начальных точек, далее, после применения операторов ГА, количество вариантов решений увеличивается в несколько раз. Затем программа автоматически запускает расчет всех этих вариантов с применением расчетного комплекса. По результатам расчета проводится оценка удовлетворения поставленным условиям и делается вывод о выполнении поставленной задачи. Генетический алгоритм применяется для анализа полученной в результате расчетов реактора информации и выработки на основе этого анализа рекомендаций по изменению параметров конструкции.
При разработке термоэмиссионного реактора-преобразователя важной задачей является снижение коэффициента неравномерности энерговыделения по радиусу активной зоны (К), так как высокий Кг влияет на значение характеристики ЭГК, что негативно отражается на способности реактора-преобразователя обеспечивать необходимую электрическую мощность [22].
ЭГК в ТРП располагаются, как правило, по концентрическим окружностям.
От их расположения зависит запас реактивности реактора и энерговыделение в каждом ЭГК. При этом, как показывают результаты расчетов, уменьшение радиального коэффициента неравномерности энерговыделения приводит к снижению запаса реактивности реактора, т. е. задача состоит в том, чтобы найти такое расположение ЭГК, при котором коэффициент неравномерности будет низким, а запас реактивности достаточным.
При моделировании активной зоны положение ЭГК задается двумя координатами. Эти координаты рассчитываются из значения радиуса окружности расположения ЭГК и их количества на отдельной окружности. Таким образом, оптимизационными параметрами становятся значения радиусов окружностей расположения ЭГК и количество ЭГК на каждой из окружностей. Также параметры могут быть дополнены, например, радиусом активной зоны, величиной обогащения топлива, углом поворота окружностей относительно друг друга.
На все параметры могут быть наложены ограничения, например, минимально возможное расстояние между ЭГК соседних окружностей, минимальное соотношение радиуса активной зоны и радиуса периферийной окружности расположения ЭГК и др.
Таким образом, решаемую оптимизационную задачу можно представить в следующим виде:
K ^ min
K .. ^max
эфф
K < K
r rm.n
R+ - R. > d
i+1 l
или
K фф > K фф .
эфф эффт.п
R+. - R. > d,
i+i i '
где К и К — минимальное
^ гтт эффшт
значение коэффициента неравномерности энерговыделения по радиусу реактора и минимально допустимое значение коэффициента размножения нейтронов, соответственно; Ri — значение радиуса окружности расположения ЭГК; d -минимально допустимое расстояние между ЭГК.
Результатами работы программы являются численные значения параметров, отвечающих за такую компоновку активной зоны реактора, расчетные значения функционалов которой отвечают поставленным требованиям.
Важным практическим результатом стало то, что для различных ТРП космического назначения удалось получить
снижение радиального коэффициента неравномерности энерговыделения с 1,30 до 1,09 [23], что приводит к увеличению генерируемой электрической мощности, при этом запаса реактивности достаточно для полной кампании. Распределение относительной тепловой мощности ЭГК в исходном и оптимизированном вариантах компоновки активной зоны ТРП показано на рис. 6.
Рис. 6. Распределение относительной мощности ЭГК:
♦ — исходный вариант; ■ — оптимизированный вариант
Применяемый подход позволяет наиболее полно исследовать оптимизационное пространство, а значит, получать наилучшие из возможных решений, что особенно важно для малогабаритных термоэмиссионных ЯЭУ.
методы нейтронно-физического расчета и оптимизации теневой радиационной защиты термоэмиссионного ядерного реактора-преобразователя
В общем контексте проектирования КЯЭУ радиационная защита занимает особое место, так как в ряде случаев именно ее габариты и масса определяют принципиальную возможность реализации проекта в рамках наложенных ограничений.
Определение габаритов радиационной защиты — это оптимизационная задача, так как с одной стороны необходимо обеспечить определенный уровень излучений за защитой, а с другой — ее минимально возможные габариты и массу.
Первоначально при оптимизации применялись методы, требующие наличия явного вида уравнения целевой функции, такие как, например, градиентные методы или методы множителей Лагранжа. При этом описание процессов проводилось с помощью сложных систем многогрупповых уравнений. Все это требовало внесения упрощений — ослабление
излучения описывалось через экспоненциальное приближение; вторичное у-излучение, вызванное радиационным захватом нейтронов в материалах радиационной защиты, не учитывалось; слои защиты рассматривались целыми, не имеющими каких-либо неоднород-ностей [24].
С появлением расчетных программных комплексов, основанных на ММК, стало возможным подробное описание радиационной защиты ТРП. В качестве примера на рис. 7 представлена такая трехмерная модель.
Рис. 7. Трехмерная модель реактора и радиационной защиты, построенная в программном комплексе MCNP:
1 — активная зона; 2 — орган регулирования; 3 — боковой отражатель; 4 — рессора привода органов регулирования; 5 — теневая радиационная защита; 6 — труба с теплоносителем; 7 — привод органов регулирования
В модели учитываются проникающие элементы, такие как трубы с теплоносителем, приводы органов регулирования. Кроме того, радиационная защита разбита на тяжелый, ослабляющий поток у-излучения, и легкий — ослабляющий поток нейтронов, компоненты.
Особенностью применяемого подхода является то, что для расчета реактора и радиационной защиты применяется один и тот же комплекс. Расчет радиационной защиты происходит при полноценной, без каких-либо упрощений или изменений, конечной модели реактора, что позволяет учесть множество тонких факторов и избежать гомогенизации элементов ядерной энергетической установки, вносящей неопределенности в получаемые значения радиационных функционалов [25].
Кроме того, разработан подход, позволяющий значительно снизить время расчета одного варианта компоновки радиационной защиты. Метод последовательных расчетов с весовыми окнами позволяет ускорить сходимость результатов в десятки и сотни раз.
Суть метода состоит в разделении всего объема расчетной модели регулярной независимой сеткой, для всех ячеек которой программным образом рассчитываются функции ценности, являющиеся весовыми окнами для каждой последующей итерации расчета переноса частиц.
Неполнота статистики компенсируется использованием среднего значения потоков в соседних ячейках, пренебрегая ослаблением в областях с плохой статистикой и тем самым практически гарантируя заниженную оценку ценности и устойчивость счета [26].
Детальное моделирование радиационной защиты позволяет создать методику ее оптимизации, основанную на физическом смысле тока контрибутонов [27]. Помещая поверхность Б в любом месте между ядерным реактором и защищаемым объектом, можно определить значение функционала радиационной нагрузки на защищаемом объекте, а также вклад в него элемента поверхности Б, что и является током контрибутонов [28].
Если в качестве поверхности Б выбрать поверхность слоя радиационной защиты и покрыть ее «черным телом», имеющим кольцевой вырез, как это показано на рис. 8, в защищаемые объекты будут попадать только частицы из кольцевой зоны поверхности радиационной защиты, над которой этот вырез расположен.
Рис. 8. Схема расчета тока контрибутонов с поверхности непрофилированного слоя: 1 — «черное тело»; 2 — кольцевая зона профилируемого слоя радиационной защиты
А значения радиационных нагрузок на них равны вкладу данной зоны в результат, получаемый со всей поверхности. Далее, смещая положение кольцевой зоны так, чтобы внутренний радиус выреза на следующем шаге равнялся внешнему, всю поверхность профилируемого слоя радиационной защиты разбивают на центральную окружность и кольцевые зоны, от которых вычисляется их индивидуальный вклад в дозу или флюенс на каждом из защищаемых объектов.
Полученные величины дозовых нагрузок для каждого положения кольцевого выреза интерполируются, полученная функция является распределением тока контрибутонов по поверхности радиационной защиты. Распределение тока контрибутонов С(г ), рассчитанное для одного из защищаемых объектов, приведено на рис. 9, а. Далее с использованием методов теории оптимального управления определяется профиль слоя радиационной защиты. Полученный таким образом профиль радиационной защиты обеспечивает допустимые дозовые нагрузки на защищаемых объектах (рис. 9, б).
Высокие значения вкладов в центре радиационной защиты связаны, в первую очередь, с наличием в защите проникающих конструкций, через которые происходит повышенная утечка излучений, а во вторую — с близостью этой зоны к оси симметрии ЯЭУ, на которой расположен защищаемый объект. Другое увеличение удельного вклада связано с утечками через трубы с теплоносителем.
Другим подходом к оптимизации характеристик радиационной защиты, разработанным в АО «ГНЦ РФ-ФЭИ», является, как и в случае оптимизации активной зоны (см. выше), соединение расчетного комплекса с оптимизационной методикой [29].
Этот подход позволяет максимально точно учитывать влияние на перенос излучения многочисленных неоднород-ностей в радиационной защите, а также влияние элементов космического аппарата, расположенных за защитой. Кроме того, данная методика не требует какого-либо аналитического описания зависимости целевой функции от варьируемых параметров. Как следствие, она исключает модельную погрешность, а также погрешность результатов косвенных расчетов, вызванную математическими операциями над величинами, определенными с погрешностью.
В качестве ограничений в данной задаче при поиске минимума массы защиты могут быть заданы допустимые потоки нейтронов и у-частиц на защищаемом объекте. В случае поиска минимальных уровней потоков нейтронов и у-квантов максимально допустимая масса защиты может быть задана в качестве ограничения.
О 2000 4 ООО 6 000 8 ООО 10000 t, см? а)
б)
Рис. 9. Результат расчета профиля радиационной защиты: а — пример рассчитанного распределения тока контрибутонов по поверхности слоя радиационной защиты как функция площади поверхности слоя; б — полученный профиль радиационной защиты (1 — реактор; 2 — труба с теплоносителем; 3 — приводы органов регулирования; 4 — теневая радиационная защита)
На рис. 10 представлены результаты работы программы для задачи оптимизации профиля и толщины легкой радиационной защиты ТРП при наличии крупногабаритного рассеивателя. При оптимизации конус радиационной защиты разбивается на блоки. Высоты каждого блока являются параметрами оптимизации и варьируются независимо друг от друга. Также в качестве параметра оптимизации выбран угол тени радиационной защиты. В результате
расчетов удалось снизить массу радиационной защиты на величину не менее 40 кг.
Полученные результаты свидетельствуют о хорошей применимости разработанных алгоритмов и компьютерных комплексов к задачам оптимизации радиационной защиты и возможности осуществления с их помощью оптимального поиска в автоматическом режиме.
б)
Рис. 10. Результаты работы методики оптимизации [23]:
а — 0 итерация, М = 160 кг; б — найденное решение, М = 118 кг; 1 — реактор; 2 — теневая радиационная защита; 3 — угол тени
Применение разработанной методики позволяет добиться снижения массы радиационной защиты уже после нахождения ее оптимальной компоновки и учесть требования к допустимым радиационным нагрузкам одновременно нескольких защищаемых объектов.
заключение
На протяжении более 60 лет в АО «ГНЦ РФ-ФЭИ» накапливаются компетенции в области процессов термоэмиссионного преобразования тепловой энергии в электричество, в физике активной зоны малогабаритных реакторов, прохождении излучения через вещество. Все эти знания в соединении с современными расчетными комплексами и новыми математическими методами позволяют проводить физико-техническое обоснование характеристик при разработке и эксплуатации инновационных КЯЭУ с прямым преобразованием тепловой энергии в электричество. Описанные методы расчетного обоснования характеристик КЯЭУ, несомненно, будут востребованы для решения новых задач, появившихся со вступлением мировой космической индустрии в эру, названную в литературе Космос 2.0 [30].
Список литературы
1. Официальный сайт-обозреватель CAD, CAE и CAM тематик. Обзор программных продуктов, реализующих CAD, CAE-технологии. Режим доступа: http:// ww w.pro cae.ru/article s/15/13 .html?showall= 1 (дата обращения: 17.01.2020 г.).
2. Пупко В.Я., Юрьев Ю.С. и др. Некоторые проблемы разработки термоэмиссионного реактора-преобразователя. Обнинск: Препринт ФЭИ-27, 1965. 20 с.
3. Дмитриев В.М., Ружников В.А. Оптимизация геометрического профилирования в термоэмиссионных электро-генерирующих каналах. Обнинск: Препринт ФЭИ-704, 1976. 19 с.
4. Ружников В.А. Численный метод совместного решения тепловой и электрической задач для термоэмиссионного электрогенерирующего канала. Обнинск: Препринт ФЭИ-774, 1977. 20 с.
5. Ружников В.А. Методы расчета тепловых и электрических характеристик систем прямого преобразования энергии: Уч. пос. Ч. 1. Термоэмиссионный
электрогенерирующий канал ЭГК. Обнинск: ФЭИ, 2001. 25 с.
6. Шиманский А.А. Эффективный алгоритм расчета ВАХ и температурных полей термоэмиссионного ЭГК на основе одномерной математической модели // Сб. тезисов докладов конференции «Ядерная энергетика в космосе». Обнинск, 1990. С. 316.
7. Синявский В.В. Методы определения характеристик термоэмиссионных твэлов. М.: Энергоатомиздат, 1990. 184 с.
8. Бабушкин Ю.В., Зимин В.П., Синявский В.В. Моделирующая система KOPTES для исследования тепловых и электрических процессов в термоэмиссионных системах преобразования энергии // Ракетно-космическая техника. Труды. Сер. XII. Королёв: РКК «Энергия». 1998. Вып. 1-2. С. 60-78.
9. Синявский В.В., Савинов А.П., Алимов В.И. и др. Имитационная модель взаимосвязанных нейтронно-физических, тепловых и электрических процессов для исследования статических, динамических и ресурсных характеристик термоэмиссионного реактора-преобразователя на быстрых нейтронах // Ракетно-космическая техника. Труды. Сер. XII. Королёв: РКК «Энергия». 1996. Вып. 2-3. С. 49-63.
10. Бабушкин Ю.В., Мендельбаум М.А., Савинов А.П., Синявский В.В. Алгоритм расчета характеристик термоэмиссионных электрогенерирующих сборок // Известия АН СССР. Энергетика и транспорт. 1981. № 2. С. 115-122.
11. Линник В.А. Расчетно-теоретические методы исследования выходных характеристик термоэмиссионных электрогенери-рующих элементов, электрогенерирующих сборок (каналов) и реакторов-преобразователей космических ЯЭУ. Обнинск: Препринт ФЭИ-3058, 2005. 70 с.
12. Виноградов Е.Г., Ярыгин В.И. Методика расчета электротеплофизи-ческих характеристик термоэмиссионного электрогенерирующего канала: Уч. пос. Обнинск: ИАТЭ, 2008. 40 с.
13. Ярыгин В.И., Ружников В.А., Синявский В.В. Космические и наземные ядерные энергетические установки прямого преобразования энергии / Монография. М.: НИЯУ МИФИ, 2016. 364 с.
14. Полоус М.А., Ярыгин В.И., Виноградов Е.Г. Программный комплекс для трехмерного численного расчета тепловых и электрических характеристик
многоэлементного электрогенерирующего канала термоэмиссионной ЯЭУ // Известия вузов. Ядерная энергетика. 2012. № 2. C. 151-160.
15. Выбыванец В.И., Гонтарь А.С., Еремин С.А. и др. Базовый электрогенерирующий канал двухрежимных термоэмиссионных ЯЭУ. Научно-технические проблемы разработки и создания // Сб. докладов Международной конференции «Ядерная энергетика в космосе 2005». Москва-Подольск, 2005. Т. 1. С. 79-82.
16. Артюхов Г.Я., Ионкин В.И, Кудрявцев В.П. и др. Оптимизация поля энерговыделения и Кэф в реакторах типа «ТОПАЗ» // Труды отраслевой юбилейной конференции «Ядерная энергетика в космосе». Обнинск, 1990. С. 24-25.
17. Пупко В.Я., Кузмин В.И. Использование функционалов теории возмущений для минимизации загрузки реакторов с произвольным спектром нейтронов // Атомная энергия. 1968. Т. 24. Вып. 3. С. 231-234.
18. Артюхов Г.Я., Истомина И.В., Макаренков Ю.Д. и др. Максимизация Кэф в гетерогенном реакторе перераспределением горючего в твэлах // Атомная энергия. 1974. Т. 37. Вып. 2. С. 135-138.
19. Полевой В.Б., Леонтьев В.В., Овчинников А.В. и др. Базовый пакет программ комплекса MMKFK-2 для решения методом Монте-Карло задач переноса нейтронов в физике реакторов (MMKFK-2-BASE). ОФАП ЯР, № 00371. М., 1996.
20. MCNP - General Monte Carlo N-Particle Transport code. LA-12625-M, Vers. 4B, 1997.
21. Гладков Л.А., Курейчик В.В., Ку-рейчик В.М. Генетические алгоритмы. М.: Физматлит, 2010. 368 с.
22. Алексеев П.А. Поиск оптимальной схемы расположения ЭГК в активной зоне термоэмиссионного реактора преобразователя космического назначения // Ядерная энергетика. 2011. № 2. С. 51-60.
23. Алексеев П.А. Создание цифрового помощника выполнения проектных расчетов // Сб. тезисов VIII научно-практической конференции молодых ученых и специалистов атомной отрасли. СПб.: Медиапапир, 2019. С. 96-97.
24. Орлов В.В., Абагян А.А., Федо-ренко Р.П., Дубинин А.А., Суворов А.П.
Оптимизация физических характеристик защиты от излучений / В сб. «Вопросы физики защиты реакторов». М.: Атом-издат, 1966. Вып. 2. С. 5-21.
25. Пышко А.П., Плотников А.Ю. Расчет и оптимизация радиационной защиты перспективных космических ЯЭУ // Атомная энергия. 2004. Т. 97. Вып. 1. С. 46-54.
26. Чернов C.B., Сонько A.B., Хо-ромский B.A. Расчет полей излучений методом итераций «весовых окон» в проекте АСММ 10/100 кВт // Радиационная защита и радиационная безопасность в ядерных технологиях. Сб. тезисов докладов 10-й юбилейной Российской научной конференции. Обнинск: НОУ ДПО «ЦИПК Росатома», 2015. С. 9-10.
27. Williams M.L. Generalized contributon response theory // Nuclear Science and Engineering. 1991. V. 10S. P. 355-3S2.
2S. Ехлаков И.А., Пышко А.П. Метод расчета тока контрибутонов с использованием
«черного тела» в задачах радиационной защиты КЯЭУ // Радиационная защита и радиационная безопасность в ядерных технологиях. Сб. тезисов докладов 10-й юбилейной Российской научной конференции. Обнинск: НОУ ДПО «ЦИПК Росатома», 2015. С. 15-16.
29. Ехлаков И.А., Пышко А.П. Автоматизированная методика поиска оптимальной компоновки радиационной защиты // Радиационная защита и радиационная безопасность в ядерных технологиях. Сб. тезисов докладов 10-й юбилейной Российской научной конференции. Обнинск: НОУ ДПО «ЦИПК Росатома», 2015. С. 14-15.
30. Дегтярева Г. Космические старт-апы: миф или реальность? Режим доступа: https://www.forbes.ru/tehnologii/342405-kosmicheskie -startapy-mi-realnost (дата обращения 17.01.2020 г.).
Статья поступила в редакцию 19.02.2020 г. Окончательный вариант — 29.06.2020 г.
Reference
1. Ofitsial'nyi sait-obozrevatel' CAD, CAE i CAM tematik. Obzor programmnykh produktov realizuyushchikh CAD, CAE-tekhnologii [Official web site reviewing the subject of CAD, CAE and CAM. A review of software products implementing CAD/CAE technologies]. Available at: http://www.procae.ru/articles/15/13.html?showall=1 (accessed 17.01.2020).
2. Pupko V.Ya, Yur'ev Yu.S. etc. Nekotorye problemy razrabotki termoemissionnogo reaktora-preobrazovatelya [Certain problems in development of a thermionic converter-reactor]. Preprint FEI-27. Obninsk, 1965. 20p.
3. Dmitriev V.M., Ruzhnikov V.A. Optimizatsiya geometricheskogo profilirovaniya v termoemissionnykh elektrogeneriruyushchikh kanalakh [Optimizing geometry profiling in thermionic electric power generating channels]. PreprintFEI-704. Obninsk, 1976. 19p.
4. Ruzhnikov V.A. Chislennyi metod sovmestnogo resheniya teplovoi i elektricheskoi zadach dlya termoemissionnogo elektrogeneriruyushchego kanala [Numerical method for joint solution of thermal and electrical problems for a thermionic power-generating channel]. PreprintFEI-774. Obninsk, 1977. 20p.
5. Ruzhnikov V.A. Metody rascheta teplovykh i elektricheskikh kharakteristik sistem pryamogo preobrazovaniya energii: Uch. pos. Ch. 1. Termoemissionnyi elektrogeneriruyushchii kanal EGK [Methods for calculating thermal and electrical properties for direct energy conversion systems: Textbook. Part 1. Thermionic Power-Generating Channel PGC]. Obninsk, FEIpubl, 2001.25p.
6. Shimanskii A.A. Effektivnyi algoritm rascheta VAKh i temperaturnykh polei termoemissionnogo EGK na osnove odnomernoi matematicheskoi modeli [An effective algorithm for calculating volt-ampere characteristic and thermal fields of a thermionic PGC based on a one-dimensional math]. Coll. of abstracts of paper presented at the conference Nuclear Power Engineering in Space. Obninsk, 1990. P. 316.
7. Sinyavskiy V.V. Metody opredeleniya kharakteristik termoemissionnykh tvelov [Methods for determining characteristics of thermionic fuel elements]. Moscow, Energoatomizdatpubl., 1990. 184p.
8. Babushkin Yu.V., Zimin V.P., Sinyavskiy V.V. Modeliruyushchaya sistema KOPTES dlya issledovaniya teplovykh i elektricheskikh protsessov v termoemissionnykh sistemakh preobrazovaniya energii [Simulation system KOPTES for studying thermal and electrical processes in thermionic energy conversion systems]. Raketno-kosmicheskaya tekhnika. Trudy. Ser. XII. Korolev, RKK «Energiya»publ., 1998, issue 1-2,pp. 60-78.
9. Sinyavskiy V.V., Savinov A.P., Alimov V.I. etc. Imitatsionnaya model' vzaimosvyazannykh neitronno-fizicheskikh, teplovykh i elektricheskikh protsessov dlya issledovaniya staticheskikh, dinamicheskikh i resursnykh kharakteristik termoemissionnogo reaktora-preobrazovatelya na bystrykh neitronakh [Simulation model of interrelated neutronic, thermal and electrical process for studying static, dynamic and endurance properties of a fast-neutron thermionic converter-reactor]. Raketno-kosmicheskaya tekhnika. Trudy. Ser. XII. Korolev, RKK «Energiya» publ., 1996, issue 2-3, pp. 49-63.
10. Babushkin Yu.V., Mendel'baum M.A., Savinov A.P., Sinyavskii V.V. Algoritm rascheta kharakteristik termoemissionnykh elektrogeneriruyushchikh sborok [An algorithm for calculating characteristics of thermionic power-generating assemblies]. Izvestiya ANSSSR. Energetika i transport, 1981, no. 2, pp. 115-122.
11. Linnik V.A. Raschetno-teoreticheskie metody issledovaniya vykhodnykh kharakteristik termoemissionnykh elektrogeneriruyushchikh elementov, elektrogeneriruyushchikh sborok (kanalov) i reaktorov-preobrazovatelei kosmicheskikh YaEU [Analytical methods for studying output characteristics of thermionic power-generating elements, powergenerating assemblies (channels) and converter-reactors of space nuclear power systems]. Preprint FEI-3058. Obninsk, 2005. 70p.
12. Vinogradov E.G., Yarygin V.I. Metodika rascheta elektroteplofizicheskikh kharakteristik termoemissionnogo elektrogeneriruyushchego kanala: Uch. pos. [A method for calculating electrical, physical and thermal properties of a thermionic power generating channel: Textbook]. Obninsk, IATEpubl., 2008. 40p.
13. Yarygin V.I., Ruzhnikov V.A., Sinyavskii V.V. Kosmicheskie i nazemnye yadernye energeticheskie ustanovki pryamogo preobrazovaniya energii [Space- and ground-based direct power conversion nuclear power systems]. Moscow, NIYaUMIFI publ, 2016.364p.
14. Polous M.A., Yarygin V.I., Vinogradov E.G. Programmnyi kompleks dlya trekhmernogo chislennogo rascheta teplovykh i elektricheskikh kharakteristik mnogoelementnogo elektrogeneriruyushchego kanala termoemissionnoi YaEU [A software package for three-dimensional numerical calculation of thermal and electrical properties of a multielement power-generating channel of a thermionic nuclear power system]. Izvestiya vuzov. Yadernaya energetika, 2012, no. 2, pp. 151-160.
15. Vybyvanets V.I., Gontar' A.S., Eremin S.A. etc. Bazovyi elektrogeneriruyushchii kanal dvukhrezhimnykh termoemissionnykh YaEU. Nauchno-tekhnicheskie problemy razrabotki i sozdaniya [Basic Power-Generating Channel of dual-mode thermionic nuclear power systems. Scientific and technical problems in development and construction]. Collected papers of the international conference Nuclear Power in Space 2005. Moscow-Podol'sk, 2005, vol. 1, pp. 79-82.
16. Artyukhov G.Ya., Ionkin V.I, Kudryavtsev V.P. etc. Optimizatsiya polya energovydeleniya i Kef v reaktorakh tipa «TOPAZ» [Optimizing power-generating field and Kef in TOPAZ-type reactors]. Proceedings of the sectorial anniversary conference Nuclear Power in Space. Obninsk, 1990. Pp. 24-25.
17. Pupko V.Ya, Kuzmin V.I. Ispol'zovanie funktsionalov teorii vozmushchenii dlya minimizatsii zagruzki reaktorov s proizvol'nym spektrom neitronov [The use of perturbation theory functionals for minimizing reactor load with an arbitrary neutron spectrum]. Atomnaya energiya, 1968, vol. 24, issue 3, pp. 231-234.
18. Artyukhov G.Ya., Istomina I.V., Makarenkov Yu.D. Maksimizatsiya Kef v geterogennom reaktore pereraspredeleniem goryuchego v tvelakh [Maximizing Kef in a heterogenic reactor through redistribution of fuel in fuel elements]. Atomnaya energiya, 1974, vol. 37, issue 2, pp. 135-138.
19. Polevoi V.B., Leont'ev V.V., Ovchinnikov A.V. etc. Bazovyi paket programm kompleksa MMKFK-2 dlya resheniya metodom Monte-Karlo zadach perenosa neitronov v fizike reaktorov (MMKFK-2-BASE) [A basic software package of the system MMKFK-2 for solving problems in neutron transfer in reactor physics using Monte Carlo method (MMKFK-2-BASE)]. OFAP YaR, № 00371. Moscow, 1996.
20. MCNP - General Monte Carlo N-Particle Transport code. LA-12625-M, Vers. 4B, 1997.
21. Gladkov L.A., Kureichik V.V., Kureichik V.M. Geneticheskie algoritmy [Genetic algorithms]. Moscow, Fizmatlit publ., 2010.368p.
22. Alekseev P.A. Poisk optimal'noi skhemy raspolozheniya EGK v aktivnoi zone termoemissionnogo reaktora — preobrazovatelya kosmicheskogo naznacheniya [Investigation into the optimal location of PGC in the active zone of a thermionic converter-reactor for space applications]. Yadernaya energetika, 2011, no. 2, pp. 51-60.
23. Alekseev P.A. Sozdanie tsifrovogo pomoshchnika vypolneniya proektnykh raschetov [Development of a digital assistant for design analysis]. Book of abstracts of the 8th workshop for young scientists and specialists working in atomic industry. Saint-Petersburg, Mediapapirpubl., 2019. Pp. 96-97.
24. Orlov V.V., Abagyan A.A., Fedorenko R.P., Dubinin A.A., Suvorov A.P. Optimizatsiya fizicheskikh kharakteristik zashchity ot izluchenii [Optimization of physical properties of radiation shielding]. In: Voprosy fiziki zashchity reaktorov. Moscow, Atomizdat publ., 1966. Issue 2, pp. 5-21.
25. Pyshko A.P., Plotnikov AYu. Raschet i optimizatsiya radiatsionnoi zashchity perspektivnykh kosmicheskikh YaEU [Analysis and optimization of radiation shielding for advanced space nuclear power systems]. Atomnaya energiya, 2004, vol. 97, issue 1, pp. 46-54.
26. Chernov S.V., Son'ko A.V., Khoromskii V.A. Raschet polei izluchenii metodom iteratsii «vesovykh okon» v proekte ASMM 10/100 kVt [Analysis of radiation fields using the method of iterations of 'weight windows' in the project ASMM 10/100 kW]. Radiation shielding and radiation safety in nuclear technologies. Coll. abstracts of papers presented at the 10th anniversary Russian scientific conference. Obninsk, NOU DPO TsIPKRosatomapubl, 2015. Pp. 9-10.
27. Williams M.L. Generalized contributon response theory. Nuclear Science and Engineering, 1991, vol. 108, pp. 355-382.
28. Ekhlakov I.A., Pyshko A.P. Metod rascheta toka kontributonov s ispol'zovaniem «chernogo tela» v zadachakh radiatsionnoi zashchity KYaEU [Method for calculating the current of contributons using «black body» in the problems of radiation shielding for space-based nuclear power systems]. Radiation shielding and radiation safety in nuclear technologies. Coll. abstracts of papers presented at the 10th anniversary Russian scientific conference. Obninsk, NOU DPO TsIPK Rosatoma publ., 2015. Pp. 15-16.
29. Ekhlakov I.A., Pyshko A.P. Avtomatizirovannaya metodika poiska optimal'noi komponovki radiatsionnoi zashchity [Automated procedure to find the optimal radiation shielding configuration]. Radiation shielding and radiation safety in nuclear technologies. Coll. abstracts of papers presented at the 10th anniversary Russian scientific conference. Obninsk, NOU DPO TsIPK Rosatoma publ, 2015. Pp. 14-15.
30. Degtyareva G. Kosmicheskie startapy: mif ili real'nost'? [Space startups: myth or reality?/ Available at: https://www.forbes.ru/tehnologii/342405-kosmicheskie-startapy-mif-ili-realnost (accessed 17.01.2020).
