CC BY
20
УДК 533.72
DOI: 10.18384/2310-7251-2018-4-127-139
АНАЛИЗ ДИНАМИКИ ПОВЕДЕНИЯ РАДИОНУКЛИДОВ В ТЕРМОЭМИССИОННОМ ПРЕОБРАЗОВАТЕЛЕ ЯЭУ НА ОСНОВЕ РЕШЕНИЯ УРАВНЕНИЯ БОЛЬЦМАНА НА КЛАСТЕРНОЙ АРХИТЕКТУРЕ
Басалаев АА.4, Клосс Ю.Ю.12, Любимов Д.Ю.13, Квасов И.Е.4, Шувалов ПВ.4, Щербаков ДВ.1, Захаров АА.4
1 Московский физико-технический институт (государственный университет) 141701, Московская область, г. Долгопрудный, Институтский переулок, д. 9, Российская Федерация
2 Национальный исследовательский центр «Курчатовский институт» 123182, г. Москва, пл. Академика Курчатова, д. 1, Российская Федерация
3 Научно-производственное объединение «Луч»
142100, Московская область, г. Подольск, ул. Железнодорожная, д. 24, Российская Федерация
4 ООО «Лаборатория открытого кода и приложений»
123060, г. Москва, ул. Берзарина, д. 19, Российская Федерация Аннотация. Данная работа посвящена моделированию поведения газообразных продуктов деления: Kr и Xe в межэлектродном зазоре ядерной энергетической установки. Для проведения анализа используется метод численного решения уравнения Больцмана с использованием двумерных сеток. Для моделирования процессов был реализован программный комплекс, позволяющий проводить вычисления на кластерной архитектуре за счёт разделения пространственных ячеек между параллельными узлами. Были получены данные по распределениям ГПД в полости межэлектродного зазора в состоянии термодинамического равновесия с заданными граничными условиями. Ключевые слова: разреженный газ, смесь газов, уравнение Больцмана, консервативный проекционный метод, численное моделирование
ANALYSIS OF THE BEHAVIOR OF RADIONUCLIDES IN THE THERMO-EMISSION CONVERTER OF A NUCLEAR POWER PLANT BASED ON THE SOLUTION OF THE BOLTZMAN EQUATION ON CLUSTER ARCHITECTURE
A. Basa^ev4, U. Kloss12, D. Lubimov13,I. Kvasov4, P. Shuvalov4, D. Sherbakov1, A. Zaharov4
1 Moscow Institute of Physics and Technology (State University) Institutskii per. 9,141701 Dolgoprudniy, Moscow region, Russian Federation
2 National Research Center "Kurchatov Institute"
pl. Akad. Kurchatova 1,123182 Moscow, Russian Federation
© CC BY Басалаев А.А., Клосс Ю.Ю., Любимов Д.Ю., Квасов И.Е., Шувалов П.В., Щербаков Д.В., Захаров А.А., 2018.
3 "Luch" Scientific Production Association
ul. Zheleznodorozhnaya 24,142100 Podolsk, Moscow region, Russian Federation
4 Open Source Lab and Applications Ltd
ul. Berzarina 19,123060 Moscow, Russian Federation Abstract. The paper is devoted to modeling the behavior of gaseous fission products: Kr and Xe in the interelectrode gap of a nuclear power plant. For analysis, the method of numerical solution of the Boltzmann equation using two-dimensional grids is used. To simulate the processes, a software package was implemented that allows calculations to be performed on cluster architecture due to separation of spatial cells between parallel nodes. Data were obtained on the distribution of the GPA in the cavity of the interelectrode gap in the state of thermodynamic equilibrium with the given boundary conditions.
Key words: rarefied gas, gas mixture, Boltzmann equation, conservative projection method, numerical simulation.
Введение
Одним из перспективных направлений в современной физике является создание космических аппаратов, способных развивать высокую энергетическую мощность, например, для круглосуточного наблюдения за техногенными и природными процессами в различных регионах земли. Уровни необходимой мощности на данный момент не могут быть достигнуты с использованием традиционных солнечных энергетических установок, поэтому создание ядерных установок является востребованной задачей. В проектных разработках космических ядерных энергетических установок получили развитие термоэмиссионные реакторы-преобразователи на тепловых и промежуточных нейтронах [2, с.20]. Одним из элементов такого реактора является электрогенерирующий канал (ЭГК), содержащий вентилируемый твэл с ядерным топливом на основе диоксида урана. На работоспособность многоэлементного ЭГК с оксидным ядерным топливом с сообщающимися полостями твэла и межэлектродного зазора (МЭЗ) и, в частности, на изменение выходных электрических характеристик влияют многочисленные факторы, например, наличие в МЭЗ продуктов деления, образующихся в топливе при работе реактора [7, с.34]. Металлические продукты деления и их оксиды попадают в МЭЗ диффундируя через ловушку и газоотводящий тракт, конденсируясь на стенках, однако газообразные продукты деления (ГПД), например, Kr и Xe не конденсируются на ловушке и газоотводящем тракте попадая в МЭЗ ЭГК и далее вакуумно-цезиевую систему [4, с.216]. В процессе диффузии радионуклидов через вакуумно-цезиевую систему происходит их бета-распад, образующиеся в результате Ba и Sr оказывают существенное влияние на выходную электрическую мощность ЭГК из-за адсорбции на стенках коллектора [3, с. 323]. Для обоснования и отработки схемы реактора проводятся испытания в петлевом канале [9, с. 50], с помощью у - спектрометрической системы, размещенной на вакуумном трубопроводе линии откачки, регистрируются радионуклиды и измеряются объемные активности откачиваемых газов. Проблемной задачей в такой схеме ЭКГ является определение давлений и потоков ГПД в ва-
куумно-цезиевой системе, поскольку оценки выхода радиоактивных газов из топлива трудоемки и не позволяют с приемлемой точностью определить потоки ГПД через газоотводное устройство в МЭЗ.
Целью работы является: определение потока ГПД в МЭЗ из газоотводящего тракта по заданному потоку ГПД из МЭЗ в вакуумно-цезиевую систему методом решения прямой задачи по выбранным начальным потокам и последующей аппроксимации результатов; и оценка парциальных давлений 8г и Ва, образующихся в полости МЭЗ в результате бета-распада.
Для анализа газообразной смеси в полости ЭГК используется численный метод решения уравнения Больцмана [5, с. 34], которое описывает поведение разряженных смесей газов. Численно уравнение Больцмана локально по пространству, что позволяет эффективно проводить расчеты на многоядерных системах, поскольку потери производительности происходят только из-за обмена информацией между соседними ячейками, распределенными в разные процессы, специальные алгоритмы разбиения пространства позволяют уменьшить число соседних ячеек и добиться почти линейного прироста производительности от числа используемых параллельных узлов.
Многоэлементный ЭГК представляет собой сложную инженерную конструкцию, которая содержит последовательно соединенные ЭГК с внутренним расположением твэлов, заключенные в герметичную оболочку, коммутационные переходники, соединяющие эмиттеры с коллекторами соседних электрогенери-рующих элементов. Внутри коммутационных переходников выполнены каналы, соединяющие полости МЭЗ и твэла (см. рис. 1).
Пространство МЭЗ заполнено газообразным цезием - рабочим телом, через каналы в эмиттере в зазор поступают криптон и ксенон, в процессе диффузии они распадаются по цепочкам бета-распада:
38Кг ^ 38ЯЬ ^ з^г;
15348 Хе ^ ^ ^Ва.
Оценочные давления Cs в полости МЭЗ составляется 150 Па, а Kr и Xe находятся в полости МЭЗ под давлением не выше 1,5 Па. Характерные размеры МЭЗ варьируются от 0,4 мм до 520 мм, это дает разброс чисел Кнудсена в диапазоне от 0,004 до 100. Такой диапазон позволяет рассматривать движение ГПД в зазоре как движение разряженных газов и позволяет моделировать их поведение методом численного решения уравнения Больцмана.
Метод решения уравнения Больцмана
Для анализа поведения ГПД используется двумерное уравнение Больцмана, а при вычислении интеграла столкновений используется приближение потенциала твердых сфер, использование двумерной формы уравнения Больцмана возможно, поскольку ширина МЭЗ мала относительно радиуса эмиттерного пакета, а сама установка осесимметрична относительно центральной оси:
f+m f+^ f = \(ffi-f )bdbd edpi.
dt m dx m dy
Для учета бета-распада ГПД уравнение Больцмана дополняется уравнением радиоактивного распада, в котором количество молекул заменено функцией распределения, поскольку при бета-распаде импульс частицы меняется слабо:
f = —ffi + 'Lfjfj
где г - индекс газа, испытывающего распад, а индексы ) - индексы газов, которые распадаются на г газ.
Для решения системы уравнений Больцмана, используется монотонная разностная схема с использованием метода расщепления по физическим процессам [1, с.49]: на первом этапе решается уравнение переноса, поочередно вдоль направлений распространения газа, на втором этапе решается уравнение релаксации, рассматриваются процессы межмолекулярного взаимодействия в газах, на третьем этапе решается уравнение бета-распада.
д- + т= 0,/* (х,$,0) = /(х,Ь) дt т дх
^г + тдУ- = 0,/* (х,$,0) = /' (х,т/2)
дt т ду
г тЛ
3t = I(r,r),r м,0)=Г I x,
,t = tj+1 — tj
df *
dt
. = — ff - +JXfr, f - (x, 5,0) = f m (x, 5, t)
f (x, 5,0) = f - (x, 5, t)
Для решения уравнения переноса используется разностная схема второго порядка точности TVD-схемой (Total Variation Diminishing Scheme). Такая схема обладает свойствами монотонности и неотрицательности решения, что отвечает поведению функции распределения, которая как физическая величина представляет из себя суть концентрацию газа.
fi — fi , p fi+1/2 — fi-1/2
-+-
m
h
= 0
fi+1/2 = '
fj +^ Afj, p > 0,
f+1 - i-IÄfLi, p < 0.
В уравнениях среднее отклонение функции распределения А/, соответствует
одному из нескольких ограничителей второго порядка точности, наиболее оптимальным выбором является зирегЬее ограничитель, при котором среднее отклонение определяется следующим образом:
Af\ = max(min (2 f+1 - fj , f/ - f -1
,min
((+1 - f/\ ,2 f' - f - 1I)
Для того, чтобы ТУБ-схема не разваливалась в процессе счета необходимо выполнение условия, исходя из которого определяется временной шаг в задаче:
Y =
p т
m h
, 0 <у< 1.
Для решения уравнения релаксации используется схема:
/ (+ т) = / () + т1 () + 0(т2).
Интеграл столкновений вычисляется с помощью проекционного метода, более подробно описанного в [10, с.53], при котором интеграл представляется в виде:
I у =
nbmaxVsph N0 N
N
X { ,Y + 5ßv ,Y - (1 - rv ) (xv ,Y + ,Y ) - rv (5xv+Sv ,Y + 5^-sv ,Y )}
v v=0
^v = {(( fHv ) (./Xv+Sv fHv-Sv ) fav fßv }
Pav Pßv m1 m
rv =
E0 v E1v -»2 ->2 -, En v = p
E2 v - Ei-
E0v = pav + Pßv , E1v = P2v + P^v , E2v = P2v+ sv + P^v-
т
В формулах выше из импульсной сетки выбираются узлы а„ в, - отвечающие за импульсы молекул до столкновения и узлы XV, ц, - отвечающие за импульсы молекул после столкновения соответственно. Причем узлы X, X, + 8^, выбираются так, что бы выполнялось неравенство:
Е < Ео < Е2
а узел ц выбирают, чтобы выполнялся закон сохранения импульса:
= «V +Pv .
Данный численный интеграл консервативен по импульсу, энергии и веществу, а также обращается в ноль на максвелловском распределении. Окончательно, для уравнения релаксации, после аппроксимации функций их дискретными аналогами, получается разностная схема:
/у' — /у + .
Эта разностная схема явная, поскольку вклады в интеграл столкновений вычисляются через функции распределения из нижнего временного слоя. Для повышения точности вычисления интеграла столкновений узлы скоростной сетки для расчета выбираются с помощью метода сеток Коробова [9, с.163].
Разностная схема дополняется начальными и граничными условиями, начальное условие - распределение Максвелла с заданным давлением Ро и температурой То:
( )_ Po ехР(-p2 / 2mkT0)
f \p' _ IT
кТо Xехр(-р2/2ткТ0)
Для определения граничного условия на границе с нормалью п нужно доопределить функцию распределения для молекул с импульсами (р, п) > 0, молекулам задается максвелловское распределение с константой, определяемой из типа граничного условия. Для условия диффузного отражения задается температура стенки Т1, а константа определяется из равенства потоков налетающих и отраженных от стенки молекул:
,„)<о1(р) Г .2 ^
f ((p,n) > 0) _ ^-W(p'w)<0--------r-exp
f 4P j ) S(p,„)>oP-exP(-p2 /2тктя ) Р
p2
2mkTi
При определении граничного условия заданного потока ql константа определяется из суммарного потока газа в ячейке и температуры газа на границе Т1:
*-Е(р,„)<0(р) ( ^
f ((n) > 0) ^ ^-r-exp
[[Р ' ' EM>op • ехР(-p2 /2mkT® ) P
p2
2mkT1
ISSN 2072-8387
2018/№ 4
При определении условия постоянного давления задается давление Р1 и температура газа Т1, функция распределения дополняется по формуле:
f ((p, n)> 0) =
kTi-1(P ,n)<f (p )
S(p,n)>0eXP (-p2/2mkT*>)
exp
,2 Л
2mkTi
Программная реализация
Ввод параметров
Ввод начальных параметров задачи 1
Обезразмериванне введенных параметров 1
Конфигурация сетки GMSH 1
Распределение ячеек между узлами
Солвер
Синхронизация
данных ячеек
1
Решение уравнения
переноса п релаксации
1
Решение уравнения
оета-распада
!
Проверка
неотрицательности
9
Вывод результатов
Вычисление
макропараметров
1
Передача данных в
основной узел
1
Восстановление
размерности
А
Вывод результатов в формате VKT
Рис. 2. Схема модулей программы.
Программный комплекс, использующий численное решение уравнения Больцмана, реализован на языке C++, программа решает задачу эволюции газовых смесей в области с заданными начальными и граничными условиями. Задание параметров пространства, степени дискретизации и разбиение сетки на равные части производится с помощью пакета GMSH. В процессе работы программы на каждом шаге итерации создаются файлы с маркопараметрами (давление, концентрация, температура, поток) ячеек в формате vtk - стандартном формате для приложений gnuplot и paraview. Схема модулей программы изображена на рисунке 2.
В программе используется одновременное исполнение на кластерных узлах с помощью средств MPI, данные пространственных ячеек разделяются по доступным параллельным узлам поровну. При проведении расчетов на нескольких узлах в пространственную сетку добавляются фиктивные ячейки, которые содержат значения и информацию о пространственной ячейке соседнего узла, перед каждой операцией по вычислению уравнения переноса идет синхронизация значений функции распределения соседних ячеек между узлами, для решения уравнения релаксации значения соседних ячеек не используются. Для
ввода и вывода данных используется узел с номером 0, который запускается на основной машине, этот узел создает пространственную и скоростную сетки и распределяет их по остальным узлам. Для получения результатов расчетов данные всех узлов собираются на основной машине и сохраняются в файл. Прирост производительности при разбиении пространственной сетки по параллельным узлам близок к линейному от числа используемых узлов, потери времени происходят из-за передачи больших массивов функции распределения между узлами в процессе синхронизации.
Рис. 3. Двумерный вид многоэлементного ЭГК: 1 - коллекторный пакет; 2 - эмиттеры; 3 - поступление ГПД через систему газоотвода.
Методом, описанным выше, решалась задача моделирования поведения ГПД в пространстве МЭЗ. Поскольку МЭЗ обладает осевой симметрией, то для упрощения счета он рассматривался в двумерном виде, представленном на рисунке 3. Граничными условиями на коллекторе и эмиттере задается диффузное отражение с температурой 900К и 1500К соответственно. Граничным условием на левом и правом концах МЭЗ является диффузное отражение с градиентной температурой и условие постоянного давления соответственно, для С8 давление устанавливается в значение 150Па, для Кг и Хе задается граничное условие вакуума - давление ОПа. Начальным условием является распределение Максвелла с давлением 150Па для С8 и градиентной температурой, Кг и Хе в начальный момент отсутствую в области счета. Поскольку оценочные концентрации Кг и Хе малы по сравнению с концентрацией С8, то при вычислении интеграла столкновений рассматривались только столкновения: С8 - С8, С8 - Кг и С8 - Хе, столкновениями ГПД между собой можно пренебречь. Расчет производится до момента установления состояния термодинамического равновесия.
Были проведены серийные расчеты для разных потоков Кг и Хе в МЭЗ, и были получены данные по распределениям давлений, температур, потоков цезия и ГПД в полости зазора. На рисунке 4 приведен график зависимости потока Кг и Хе на правом конце МЭЗ от начального потока из газоотводящего тракта. На выходе из МЭЗ задавался ожидаемый конечный поток дк.
0.5 1.0 1.5 2.0 2.5 3.0 2 4 6
ЯоДОЧмх)1 ЧоД015 (м2с)'
Рис. 4. Конечный поток из МЭЗ криптона (а) и ксенона (б) в зависимости от их начального потока в МЭЗ.
Таблица 1.
Конечные и начальные потоки ГПД в МЭЗ
Газы qk, 1016(м2 • с)-1 qo, 1015(м2 • с)-1
Кг 1,01 1,29
Хе 1,73 3,74
Используя метод аппроксимации, т.е. минимизации квадратичной ошибки:
Е = £ | РX)-Уi |2
;=0
где р(х) - это аппроксимирующий полином, а у - это значения дискретной функции в точках ), определяется зависимость соответствующая полученным данным, из этой зависимости определяется значение начального потока ГПД в МЭЗ для заданного ожидаемого потока на выходе. Аппроксимация производилась с помощью пакета питру на локальной машине. Расчетные начальные и ожидаемые конечные потоки указаны в таблице 1.
Для оценки парциальных давлений Ва и Бг были проведены расчеты с начальным потоком, определенным методом аппроксимации, при расчетах рассматривался бета-распад Кг и Хе до Бг и Ва. При вычислении интеграла столкновений рассматривались только столкновения: Сэ - сб, Сэ - Кг и Сэ - Хе. Всего рассматривалось 7 газов: Сэ, Кг, Хе, ИЬ, Бг, С8(138), Ва. На рисунке 5 приведены графики распределения давления Кг, Хе и Бг, Ва относительно длинны МЭЗ.
Рис. 5. Давление Kr, Xe (а) и их продуктов бета-распада Sr, Ba (б) в МЭЗ.
Заключение
В ходе работы было рассмотрено применение метода численного решения уравнения Больцмана к расчету поведения газообразной смеси в полости МЭЗ, разработана программа и получены распределения потоков и давлений ГПД для разных начальных потоков. Рассмотрена зависимость потока ГПД на конце МЭЗ от потока из газоотводящего тракта, получено значение начального потока по заданному конечному. Оценены парциальные давления продуктов бета-распада Kr и Xe в полости зазора. Полученные результаты, согласуются с результатами, полученными методом решения линейного уравнения диффузии [8, с. 324].
Статья поступила в редакцию 20.08.2018 г.
БЛАГОДАРНОСТИ
Работа выполнена при финансовой поддержке РФФИ, гранты 17-08-00606, 18-08-00438,18-07-00323,17-08-00339,17-08-00252.
ACKNOWLEDGMENTS
This work was supported by the Russian Foundation for Basic Research (grant Nos 17-08-00606, 18-08-00438, 18-07-00323, 17-08-00339, and 17-08-00252).
ЛИТЕРАТУРА
1. Аристов В.В., Черемисин Ф.Г. Расщепление неоднородного кинетического оператора уравнения Больцмана // Доклады АН СССР. 1976. Т. 231. № 1. С. 49-52.
2. Васильковский В.С., Андреев П.В., Зарицкий Г.А. Проблемы космической энергетики и роль ядерных энергетических установок в их решении // Международная конференция «Ядерная энергетика в космосе - 2005». Москва - Подольск, 1-3 марта 2005 г.: сборник докладов в 3 т. Т. 1. М.: Научно-исследовательский и конструкторский институт энерготехники им. Н.А. Доллежаля, 2005. С. 20-25.
3. Выбыванец В.И., Любимов Д.Ю. Корюкин В.А. Моделирование работы долгоресурс-ного термоэмиссионного преобразователя с ядерным нагревом // Атомная энергия. 2015. Т. 118. Вып. 4. С. 233-236.
4. Гонтарь А.С., Гриднев А.А., Любимов Д.Ю. Анализ физико-химических процессов в многоэлементном ЭГК с сообщающимися полостями твэла и межэлектродного зазора // Атомная энергия. 2008. Т. 104. Вып. 4. С. 216-224.
5. Коган М.Н. Динамика разреженного газа. М.: Наука, 1967. 440 с.
6. Коробов Н.М. Тригонометрические суммы и их приложения. М.: Наука, 1989. 240 с.
7. Любимов Д.Ю., Николаев Ю.В., Шумилов А.А. Влияние продуктов деления на фазовый состав достехиометрического диоксида урана в тепловыделяющих элементах термоэмиссионных электрогенерирующих каналов // Материаловедение. 2008. № 3(132). С. 34-42.
8. Любимов Д.Ю., Федик И.И., Шумилов А.А. Влияние продуктов деления на выходную мощность термоэмиссионных ЭГК с сообщающимися и разделенными полостями твэла и межэлектродного зазора // Атомная энергия. 2011. Т. 110. Вып. 6. С. 321-327.
9. Синявский В.В. Методы и средства экспериментальных исследований и реакторных испытаний термоэмиссионных электрогенерирующих сборок. М.: Энергоатомиздат.
10. Черемисин Ф.Г. Консервативный метод вычисления интеграла столкновений Больцмана // Доклады РАН. 1997. Т. 357. № 1. С. 53-56.
1. Aristov V.V., Cheremisin F.G. [Separation of the inhomogeneous kinetic operator in the Boltzmann equation]. In: Doklady AN SSSR [Doklady Mathematics], 1976, vol. 231, no. 1, pp. 49-52.
2. Vasil'kovskii V.S., Andreev P.V., Zaritskii G.A. [Problems of space energy and the role of nuclear power plants in their decision]. In; Mezhdunarodnaya konferentsiya «Yadernaya energetika v kosmose - 2005». Moskva-Podolsk, 1-3 marta 2015g.; sbornik dokladov v 3 t. T. 1 [International conference "Nuclear power in space - 2005". Moscow-Podolsk, March 1-3, 2005: a collection of papers in 3 volumes. Vol. 1]. Moscow, Dollezhal Research and Development Institute of Power Engineering Publ., 2005. pp. 20-25.
3. Vybyvanets V.I., Lyubimov D.Yu., Koryukin V.A. [Simulation of the operation of a long-life thermionic converter with nuclear heating]. In: Atomnaya energiya [Atomic Energy], 2015, vol. 118, no. 4, pp. 233-236.
4. Gontar' A.S., Gridnev A. A., Lyubimov D.Yu. [Analysis of physico chemical processes in multicomponent electricity generating channels with communicating fuel-element and interelectrode gaps]. In: Atomnaya energiya [Atomic Energy], 2008, vol. 104, no. 4, pp. 216224.
5. Kogan M.N. Dinamika razrezhennogo gaza [Rarefied gas dynamics]. Moscow, Nauka Publ., 1967. 440 p.
6. Korobov N.M. Trigonometricheskie summy i ikhprilozheniya [Trigonometric sums and their applications]. Moscow, Nauka Publ., 1989. 240 p.
7. Lyubimov D.Yu., Nikolaev Yu.V., Shumilov A.A. [Influence of fission products on phase composition of prestoichiometric uranium dioxide in fuel elements of thermal-emission energy-generating channels]. In: Materialovedenie [Inorganic Materials: Applied Research], 2008, no. 3(132), pp. 34-42.
8. Lyubimov D.Yu., Fedik I.I., Shumilov A.A. [Effect of fission products on the output power of thermionic EGC with communicating and separated fuel-element cavities and
2000. 355 с
REFERENCES
an interelectrode gap]. In: Atomnaya energiya [Atomic Energy], 2011, vol. 110, no. 6, pp. 321-327.
9. Sinyavskii V.V. Metody i sredstva eksperimental'nykh issledovanii i reaktornykh ispytanii termoemissionnykh elektrogeneriruyushchikh sborok [Methods and means of experimental studies and reactor tests of thermionic power generating assemblies]. Moscow, Energoatomizdat Publ., 2000. 355 p.
10. Cheremisin F.G. [A conservative method for computing the Boltzmann collision integral]. In: Doklady RAN [Doklady Mathematics], 1997, vol. 357, no. 1, pp. 53-56.
ИНФОРМАЦИЯ ОБ АВТОРАХ
Басалаев Артем Викторович - начальник группы ООО «Лаборатория открытого кода и
приложений»;
e-mail: bz@lokip.ru;
КлоссЮрий Юрьевич - доктор физико-математических наук, начальник отдела моделирования физических процессов Национального исследовательского центра «Курчатовский институт»; заведующий кафедрой моделирования ядерных процессов и технологий Московского физико-технического институт (государственного университета); e-mail: kl@lokip.ru;
Любимов Дмитрий Юльевич - кандидат технических наук, доцент кафедры моделирования ядерных процессов и технологий Московского физико-технического институт (государственного университета); ведущий научный сотрудник научно-производственного объединения «Луч»; e-mail: ld@lokip.ru;
Квасов Игорь Евгеньевич - кандидат физико-математических наук, начальник группы ООО «Лаборатория открытого кода и приложений»; e-mail: kvs@lokip.ru;
Шувалов Павел Вадимович - кандидат физико-математических наук, начальник лаборатории ООО «Лаборатория открытого кода и приложений»; e-mail: shv@lokip.ru;
Щербаков Дмитрий Валерьевич - аспирант кафедры моделирования ядерных процессов и технологий Московского физико-технического института; e-mail: dmitrii.shcherbakov@phystech.edu;
Захаров Александр Алексеевич - начальник группы, ООО «Лаборатория открытого кода и приложений»; e-mail: shv@lokip.ru.
INFORMATION ABOUT THE AUTHORS
Artem V Basalaev - Head of the Group, Open Source Lab and Applications Ltd; e-mail: bz@lokip.ru;
Yury Yu. Kloss - Doctor in Physical and Mathematical Sciences, Head of the Department of Physical Processes Modeling, National Research Center "Kurchatov Institute"; Head of the Department of Modeling Nuclear Processes and Technologies, Moscow Institute of Physics and Technology (State University); e-mail: kl@lokip.ru;
Dmitry Yu. Lubimov - PhD in Engineering Sciences, Associate Professor at the Department of Modeling Nuclear Processes and Technologies, Moscow Institute of Physics and Technology (State University); Leading Researcher of "Luch" Scientific Production Association; e-mail: ld@lokip.ru;
Igor Ye. Kvasov - PhD in Physical and Mathematical Sciences, Head of Group, Open Source Lab and Applications Ltd; e-mail: kvs@lokip.ru;
Pavel V. Shuvalov - PhD in Physical and Mathematical Sciences, Head of the Laboratory Open Source Lab and Applications Ltd; e-mail: shv@lokip.ru;
Dmitry V. Sherbakov - postgraduate student at the Department of Modeling Nuclear Processes and Technologies, Moscow Institute of Physics and Technology (State University); e-mail: dmitrii.shcherbakov@phystech.edu;
Alexander A. Zaharov - Head of the Group, Open Source Lab and Applications Ltd; e-mail: shv@lokip.ru.
ПРАВИЛЬНАЯ ССЫЛКА НА СТАТЬЮ
Басалаев А.А. Клосс Ю.Ю., Любимов Д.Ю., Квасов И.Е., Шувалов П.В., Щербаков Д.В., Захаров А.А. Анализ динамики поведения радионуклидов в термоэмиссионном преобразователе ЯЭУ на основе решения уравнения Больцмана на кластерной архитектуре // Вестник Московского государственного областного университета. Серия: Физика-Математика. 2018. № 4. С. 128-139. DOI: 10.18384/2310-7251-2018-4-127-139
FOR CITATION
Basabev A.A., Kloss Yu.Yu., Lubimov D.Yu., Kvasov I.Ye., Shuvalov P.V., Sherbakov D.V., Zaharov A.A. Analysis of the behavior of radionuclides in the thermo-emission converter of a nuclear power plant based on the solution of the Boltzmann equation on cluster architecture. In: Bulletin of Moscow Region State University. Series: Physics and Mathematics, 2018, no. 4, pp. 128-139.
DOI: 10.18384/2310-7251-2018-4-127-139
