Научная статья на тему 'Подход к математическому моделированию трансмутационных процессов в ядерных энергетических установках'

Подход к математическому моделированию трансмутационных процессов в ядерных энергетических установках Текст научной статьи по специальности «Математика»

CC BY
388
83
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЗАДАЧА ЯДЕРНОЙ ТРАНСМУТАЦИИ / МАТЕМАТИЧЕСКАЯ МОДЕЛЬ / ОРИЕНТИРОВАННЫЙ МУЛЬТИГРАФ / СИСТЕМА ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ / NUCLEAR TRANSMUTATION PROBLEM / MATHEMATICAL MODEL / DIRECTED MULTIGRAPH / ORDINARY DIFFERENTIAL EQUATIONS SYSTEM

Аннотация научной статьи по математике, автор научной работы — Белозерова Алла Равильевна, Мельников Борис Феликсович

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

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

Похожие темы научных работ по математике , автор научной работы — Белозерова Алла Равильевна, Мельников Борис Феликсович

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

APPROACH TO MATHEMATICAL SIMULATION OF TRANSMUTATIONAL PROCESSES IN NUCLEAR POWER PLANTS

Background. At the present time the study and simulation of nuclear kinetics or nuclide transformations in the process neutron irradiation in power reactor facilities are especially topical for solution of the problem of nuclear fuel cycle abridgement to provide profitability of nuclear power. The purpose of this article consists in describing and investigating essentially new mathematical models in simulation of nuclear-physical processes in reactor materials irradiation. Materials and methods. The initial model of the problem is similar to a linear system of the ordinary differential equations. Relations of nuclide reciprocity and accumulation on a nuclide set of nuclide strong-coherent subsets were entered. With the help of relations a transition from the initial continuous model to a discrete model was made. The discrete model describes nuclear kinetics as an oriented multigraph. The authors averaged the estimation of neutron-physical characteristics of reactor irradiation by means of averaging a stream at certain micro-campaigns in the process of power reactor facility functioning. Results. The problem of discrete optimization of the nuclide transformations was formulated in terms of the graph theory. The discrete problem is original in the discrete optimization theory. The problem concerns a class of NP-difficult problems. The authors author a method of averaging the estimation of neutron-physical characteristics of reactor irradiations for the period of several micro-campaigns in the problem of nuclear transmutations. Conclusions. The priority direction of modernization and technological development of the Russian economy consists in development of nuclear technologies, namely in realization of the closed fuel cycle with mixed fuel in fast power reactors. The considered model of calculation of nuclide transformations into materials at reactor irradiation can be successfully used for mathematical modelling of nuclear-physical processes with accumulation and burning out actinides in a fast reactor core.

Текст научной работы на тему «Подход к математическому моделированию трансмутационных процессов в ядерных энергетических установках»

УДК 519.176, 519.622.2, 519.688

А. Р. Белозерова, Б. Ф. Мельников

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

Аннотация.

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

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

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

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

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

A. R Belozerova, B. F. Mel'nikov

APPROACH TO MATHEMATICAL SIMULATION OF TRANSMUTATIONAL PROCESSES IN NUCLEAR POWER PLANTS

Abstract.

Background. At the present time the study and simulation of nuclear kinetics or nuclide transformations in the process neutron irradiation in power reactor facilities are especially topical for solution of the problem of nuclear fuel cycle abridgement to provide profitability of nuclear power. The purpose of this article consists in describing and investigating essentially new mathematical models in simulation of nuclear-physical processes in reactor materials irradiation.

Materials and methods. The initial model of the problem is similar to a linear system of the ordinary differential equations. Relations of nuclide reciprocity and accumulation on a nuclide set of nuclide strong-coherent subsets were entered. With the help of relations a transition from the initial continuous model to a discrete model was made. The discrete model describes nuclear kinetics as an oriented multigraph. The authors averaged the estimation of neutron-physical characteristics of reactor irradiation by means of averaging a stream at certain micro-campaigns in the process of power reactor facility functioning.

Results. The problem of discrete optimization of the nuclide transformations was formulated in terms of the graph theory. The discrete problem is original in the discrete optimization theory. The problem concerns a class of NP-difficult problems.

The authors author a method of averaging the estimation of neutron-physical characteristics of reactor irradiations for the period of several micro-campaigns in the problem of nuclear transmutations.

Conclusions. The priority direction of modernization and technological development of the Russian economy consists in development of nuclear technologies, namely in realization of the closed fuel cycle with mixed fuel in fast power reactors. The considered model of calculation of nuclide transformations into materials at reactor irradiation can be successfully used for mathematical modelling of nuclear-physical processes with accumulation and burning out actinides in a fast reactor core.

Key words: nuclear transmutation problem, mathematical model, directed multigraph, ordinary differential equations system.

1. Постановка задачи

Влияние процессов ядерной трансмутации ([1] и др.) на основные ядерно-физические характеристики материалов ядерных энергетических установок (ЯЭУ) требует детального изучения и количественной оценки изменения нуклидного состава этих материалов при реакторном облучении. При этом значимость трансмутационных процессов наиболее ярко выражена для тех материалов, в которых изменение химического состава может непосредственно влиять на изменение эксплуатационных свойств материалов и приводить к ограничению ресурса изделий, ответственных за безопасную эксплуатацию ЯЭУ, что в первую очередь связано с материалами отражателя, замедлителя и поглотителя нейтронов. К таким материалам относится бериллий, гафний и их соединения.

Накопление газообразных продуктов ядерной трансмутации при длительном реакторном облучении этих материалов может приводить к изменению их структурно-фазового состояния, что также необходимо учитывать при обосновании работоспособности конструктивных элементов ЯЭУ. Для конструкционных материалов (типа коррозионно-стойких сталей и никельсодержащих материалов [1]) важно определение динамики накопления радиогенных изотопов гелия и водорода в зависимости от условий реакторного облучения.

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

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

В нашем случае общая проблематика такова. Согласно [4] трансмутации в ядерных реакторах происходят в результате ядерных реакций, вызываемых нейтронами, а также в результате спонтанных распадов радионуклидов. То есть под воздействием нейтронного облучения происходит трансмутация нуклидного состава материала, протекают процессы активации и газонакоп-ления - и все это приводит к изменениям нуклидного состава. Задача расчета этих изменений (трансмутаций) возникает в широком круге теоретических исследований и прикладных работ, связанных с ядерными реакторами.

2. Основные понятия. Исходные уравнения образования нуклидов

Под скоростью реакции Ап превращения нуклида / в нуклид п понимается доля ядер нуклида /, превратившихся в нуклид п в единицу времени. Скорость реакции Ап определяется интегрированием по энергии нейтронов:

Ап = |а1п (Е)Ф(Е)ёЕ, где а1п (Е) - дифференциальное сечение реакции,

о

Ф(Е) - дифференциальная плотность потока нейтронов. В данной работе

радиоактивные распады рассматриваются наряду с ядерными реакциями как каналы нуклидных превращений. Скорость распада радионуклида / в нуклид

л* Л 1п(2) 1 т

п записывается так: Ап = Л =----, где Л - постоянная распада; Т1/2 - период

ТИ2

полураспада.

Пусть х* - концентрация ядер нуклида /. Исходные уравнения нуклид-

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

\ddj_=- £А;х,, * = 1,2...^. (1)

Физический смысл отдельных членов в уравнении (1) таков:

- член слева - изменение концентрации ядер х* нуклида / в единицу времени;

- первый член справа - концентрация ядер нуклида /, образовавшегося путем всех возможных превращений к ^ / за единицу времени (суммирование по к распространяется на все случаи, для которых Ак ф 0);

- второй член справа - концентрация ядер нуклида /, уничтоженных за единицу времени (суммирование по I распространяется на все случаи, для которых А[ ф 0);

- N - общее число нуклидов в цепочке трансмутации [4].

При этом х* - концентрация ядер, в качестве ее размерности принято рассматривать млн -1. Таким образом, она обозначает количество определенного вида ядер на 106 ядер состава. Начальные условия для задачи трансмутации определяются стартовым нуклидным составом материала и записываются в виде

Гх* о, если / содержится в исходном веществе,

х(0) = ь

[0, если нет,

где хг-,0 - концентрация ядер нуклида / в начальный момент времени.

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

3. Задача ядерной трансмутации в терминах теории графов

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

Итак, сначала приведем несколько определений относительно структуры схемы нуклидных превращений.

Определение 1. Набором нуклидов называется множество всех нуклидов, включенных в схему превращений.

Определение 2. Два нуклида из набора называются взаимными, если в графе существуют направленные пути из первого нуклида во второй и из второго нуклида в первый (при этом рассматриваются и пути нулевой длины, т.е. всякий нуклид взаимен самому себе).

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

Определение 3. Класс взаимных нуклидов называется компонентой, или блоком.

Определение 4. Компонента В накапливается из компоненты С тогда и только тогда, когда для некоторого нуклида с из С и нуклида Ь из В в графе существует направленный путь из с в Ь.

Бинарное отношение накопления на множестве компонент транзитивно. Если компонента В накапливается из компоненты С и С накапливается из В, то В = С. Следовательно, отношение накопления антисимметрично и является на множестве блоков отношением частичного порядка.

Очевидно, что каждая компонента является сильно связанным подмножеством вершин ориентированного мультиграфа О.

4. Примеры представления графов в задаче ядерной трансмутации

В рассматриваемой нами задаче ядерной трансмутации изучается система связей между различными объектами (нуклидами). Объекты являются вершинами и отмечаются точками, а связи между вершинами являются дугами и отмечаются стрелками с соответствующими точками (рис. 1,а).

(n,t)

а)

9Be

11Be

(3-decay

Все применяемые далее понятия теории графов согласованы с [5]. Графом называется алгебраическая система О = (М, Я, где Я - двухместный предикатный символ. Элементы носителя М называются вершинами графа О, а элементы бинарного отношения Я с М - дугами. Таким образом, дугами являются пары вершин (а,Ь)е Я . При этом дуга (а,Ь) называется исходящей из вершины а и заходящей в вершину Ь.

Пример 1. Изображение графа О с множеством вершин М ={ 7 Ь/,9 Ве}

и множеством дуг Я = {(п,пё),(п,t)} представлено на рис. 1,а, где (п, пё) -

захват нейтрона п с испусканием нейтрона п и дейтерия ё (ё = 2Н); (п, 0 - захват нейтрона п с испусканием трития t (t = 3Н). В задаче представления трансмутаций в некотором материале при облучении в виде графа имеется несколько дуг: (п, 0, (п, пё), исходящих из вершины 9Ве и заходящих в вершину 7Ы, - такие дуги называются кратными (рис. 1,а). При этом мы используем понятие мультиграфа.

Мультиграфом О называется тройка {М,и,Р), в которой М - множество вершин, и - множество дуг, а Р с М Xи XМ - 3-местный предикат, называемый инцидентором и представляемый следующим образом: (а,и,Ь)е Р тогда и только тогда, когда дуга и исходит из вершины а и заходит в вершину Ь. Граф О = (М, Я^ называется ориентированным (орграфом), если найдется дуга (а,Ь)е Я , такая что (Ь,а)й Я .

Пусть Бм , Бр - множества меток. Пометкой или распределением меток мультиграфа О = (М,и,Р) называется пара функций /:М ^ Бм (распределение меток вершин) и g : Р ^ Бр (распределение меток дуг). Пятерка (М ,и, Р, /, ^ называется взвешенным или помеченным графом. Для вершины а е М элемент /(а) называется весом вершины а, а для дуги (а/,и,аj) е Р элемент g (,и,аj) - весом дуги (,и,аj). Каждая вершина

графа из множества М соответствует определенному нуклиду; при этом в дальнейшем изложении термины «нуклид» и «вершина» являются синонимами. Также синонимами являются термины «дуга графа» из множества и и «канал превращения нуклида» (распадам радионуклидов и реакциям нейтронов с ядрами), для которых (а,и,Ь)е Р тогда и только тогда, когда дуга и исходит из вершины а и заходит в вершину Ь. Для вершины а е М элемент Да) - вес вершины а, соответствующий концентрации определенного для вершины а нуклида, т.е. количество атомов нуклида на 106 атомов состава (млн -1). Каждая дуга графа соответствует каналу превращения, обусловленному ядерной реакцией или распадом, т.е. для дуги ие и элемент g(u) - вес дуги и, который определяется значением скорости превращения. Размерность скорости - с-1, обычный диапазон значений - [10 -30...10 -15 ].

Пример 2. Пусть

М ={ 9Ве, 10Ве,11Ве,10 в} ,

U = {(n,у),(n,2n),(Р-decay)} ,

P = {(9Be (n, y)10Be), (10Be (n,2n) 9Be),

(10Be (n, y) 11Be), (10Be (p-decay )10b)} , f: M ^ SM,

где Sm - множество меток со значениями концентраций вершин-нуклидов для природной композиции бериллия; g : P ^ Sp , где Sp - множество меток со значениями скоростей ядерных реакций в бериллии в условиях облучения ячейки 52 активной зоны реактора СМ с газовым заполнением;

f (9Be ) = 106, f (10Be ) = 0, f (Be ) = 0, f (10B ) = 0, g ((9Be (n, y) 10Be )) = 1,84-10-12, g ((10Be (n,2n )9Be )) = 1,03-10-12,

g ((10Be (n, y) 11Be )) = 2,82-10-13, g ((10Be (p-decay )10b)) = 1,37-10-14.

Помеченный граф (M,U,P, f, g} изображен на рис. 1,б и представляет

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

n

^ f (mi) = 106 - const, mt е M.

i=1

Пример 3. Граф, представленный на рис. 1,б, имеет три сильно связанные компоненты, задаваемые множествами вершин {9Be,10Be}, {11Be} и

{10B} - поскольку матрицы достижимости C и контрдостижимости Q графа имеют следующий вид:

(bj ) = E + Ag + 4 +... + AG =

f 1 0 0 0 > f 0 1 0 0 " 1 0 1 1 ^ f 0 1 0 1 ^

0 1 0 0 1 0 1 1 0 1 0 0 1 0 1 1

+ + +

0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0

V 0 0 0 1V V 0 0 0 0 V 0 V 0 0 0 V V 0 0 0 0 V

f 1 0 1 1Л f 3 2 2 3

0 1 0 0 2 3 2 2

+ =

0 0 0 0 0 0 1 0

V 0 0 0 0 V V 0 0 0 1 /

C=

Г11111 1111 0 0 10 0 0 0 1

Q = ct =

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

Гі l 0 01 110 0 1110 110 1

, C x Q =

Г і і 0 01 110 0 0 0 10 0 0 0 1

5. Физические основы стационарной задачи ядерной трансмутации

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

В расчетах стационарной задачи ядерной трансмутации изменения нейтронного спектра не учитываются, поэтому согласно [4] и приведенным выше обозначениям выпишем цепочку накопления произвольного радионуклида за весь период облучения произвольного конструкционного материала, предполагая, что период облучения состоит из временных промежутков работы реактора на мощности и его остановок:

Xi = X0 • e X2 = Xi • e-

„—Xt2

X3 = X2 • e

X4 = X3 • e-

xn = x0 • e'

Rlti —Xt^ +R9t3 —Xt4 +...—Rkt~> к—1 —Xt'

,i2^2l3

к 2k—1■

■4 к

(2)

„—Xt4

При этом для решения стационарной задачи ядерной трансмутации нужно усреднить значение скорости накопления радионуклида. Для этого вводится среднее значение скорости в виде Я; таким образом, формула определения концентрации искомого радионуклида переписывается в виде Яt

хп = Х0 • е , где Я - усредненное значение скорости образования радионуклида с учетом ядерного накопления и убыли искомого нуклида. Выразим скорость образования произвольного радионуклида через среднее сечение и плотность потока нейтронов реакторного облучения, получим

- Г х 1

а-ф-1 = 1п — . С другой стороны, из (2) аналогично проведенным рассуж-

I х0 )

дениям получим

0 •ф • t = R1 • ti — X • t2 + R2 • t3 ■

t4 + ...— Rk • t2k—1 —t2k .

Заметим, что общее время расчета t складывается из периодов облучения t\, ^, t5, t7, ..., ^к_ь где к = (п + 1)/2 - количество так называемых микро-

кампаний [4], в течение которых искомый материал находился под реакторным облучением, а также периодов остановки реактора ^, ^, t6, ^, ..., ^к между соответствующими микрокампаниями, т.е.

t = ^ +1:2 + ^ + ... + ^2 к-1 + Ьк = ^ Ч .

1

Поэтому

О- ф-1 = Я1 • ^ + Я2 -1з + ... -Як - t2k-1 -Х ^2 -Х ^4...-Х ^2к , или О- ф-1 = °- ф1 ^1 +а- ф2 -1з + ... + О- фк - Ьк-1 -Х- Ь -Х-14...-Х- ^2к ,

О- ф1 ^1 + °- ф2 - 1з + ... + О- фк - t2k-1 - Х- ^ - Х- ^^4... -Х- Ьк

Ф = -

а-г1

Фном '( ( + г2 -t3 + ••• + jk-t2k-1 ) Мг2 + г4 + ••• + {2к ) (3)

t а-1 ’

E

л max

где а = — f а(E) ■ f (E)dE • m J

Ф

Е0

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

фном -(1 (1 + *2 ^3 +... + 1к^2к-1) t ,

где /ь /2, /3, /4, ..., 4 - доля отклонения плотности нейтронного потока от номинала по микрокампаниям.

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

6. Специфика методической и вычислительной погрешностей

Диапазон скоростей реакций весьма широк, это приводит к системам дифференциальных уравнений с высоким коэффициентом жесткости - до ~10 . Для численного решения применен метод полного вычисления собственных значений и собственного базиса матрицы системы. Матрица системы дифференциальных уравнений приводится посредством подходящей нумерации нуклидов к блочно-треугольному виду с размерностью диагональных блоков не более заданной величины Вттах. Для большинства схем с числом нуклидов, превышающим Б1ттах, требуется упростить схему, чтобы удовлетворить этому требованию. Размерность системы в расчетах может достигать 50 и более. Отсюда возникает задача дискретной оптимизации: выполнить необходимые упрощения схемы так, чтобы дополнительная методическая погрешность, обусловленная сделанными упрощениями, была минимальной. Чтобы обеспечить приемлемые значения вычислительных погрешностей, на сложность схемы накладываются ограничения. При выполнении

данных ограничений матрица системы имеет блочно-треугольный вид, размерность диагональных блоков не превышает установленного максимального значения (в последней версии компьютерной программы максимальная допустимая размерность блока Dimmax = 10) Задача вычисления собственных значений и собственного базиса факторизуется по диагональным блокам • Такой подход преследует цель достичь баланса методической и вычислительной погрешностей в расчете трансмутаций Блочно-треугольный вид матрицы позволяет снизить вычислительную погрешность, но для этого необходимо упростить схему нуклидных превращений

7. Математическая постановка задачи ядерной трансмутации

Схема ядерной трансмутации в облучаемом материале представляет взвешенный (помеченный) ориентированный мультиграф G = (М,U,P, f,g •

Путем операции добавления к графу G = (М,U,P, f,g вершины Other образуется граф G' = (M ujOther} ,U, P, f, • Таким образом, в схему ядерной

трансмутации вводится дополнительная коллекторная вершина, в которую перенаправляются некоторые каналы превращений Значение переменной в системе уравнений (1), соответствующей коллекторной вершине, дает динамическую гарантированную оценку методической погрешности [6, 7] Задача дискретной оптимизации схемы ядерной трансмутации состоит в том, чтобы минимизировать методическую погрешность, обусловленную упрощениями схемы для обеспечения блочно-треугольного представления матрицы •

Иными словами, задачу ядерной трансмутации можно описать следующим образом • Удаляя из заданного ориентированного мультиграфа G = (М,U,P,f,g дуги, нужно прийти к такому подграфуG', в котором количество вершин в любом блоке не превышало бы заданного числа Dimmax, и при этом P/ с P • С учетом того, что для каждой отбрасываемой дуги (a,u,b)

потеря концентрации равна интегралу превращения по соответствующему каналу за рассматриваемое время облучения материала, необходимо отброшенные дуги перенаправить в коллекторную вершину Other • Из введенных нами определений следует, что операции удаления дуги (a,u,b) и добавления дуги (a,u,Other) к графу G' состоят в образовании графа

G* = ^М ujOther} ,U, P u{(a, u, Other)} \ {(a, u, b)}, f, gj,

причем элемент g(u) - вес дуги u сохраняется^ Цель задачи ядерной трансмутации состоит в минимизации значения концентрации коллекторной вершины Other для графа G" = (м ujOther} ,U, P uj(a, u, Other)} \ {(a, u, b)}, f, gj •

Заключение

Возникает вопрос: что необходимо сделать для возможности практического решения описанных в данной статье NP-трудных проблем с помощью детерминированных алгоритмов • Для этого мы в следующей публикации предполагаем рассмотреть эвристический подход ([8, 9] и др-), при котором

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

Например, обычно условно можно считать, что:

- для химии, металлургии сплавов и смежных с ними наук в качестве малой методической погрешности принимается значение 1000 млн -1;

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

Итак, в следующей публикации предполагается:

- привести описание точного алгоритма для ЗДО малых размерностей;

- описать подходы к построению эвристических алгоритмов решения сформулированной нами ЗДО; прежде всего применение незавершенного метода ветвей и границ, а также алгоритмов локального поиска;

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

Список литературы

1 Ма, Б. М. Материалы ядерных энергетических установок / М^ Ма^ - М^ :

Энергоатомиздат, 1987^ - 408 с 2^ Громкович, Ю. Теоретическая информатика^ Введение в теорию автоматов, теорию вычислимости, теорию сложности, теорию алгоритмов, рандомизацию, теорию связи и криптографию / Ю^ Громкович^ - СПб^ : БХВ-Петербург, 2010^ -336^

3^ Hromkovic, J. Algorithmics for Hard Problems: Introduction to Combinatorial Optimization, Randomization, Approximation, and Heuristics / I Hromkovic - 2nd Ed^ -Springer, 200L

4^ Герасимов, А. С. Справочник по образованию нуклидов в ядерных реакторах / А^ С Герасимов, Т С Зарицкая, А П Рудик - М^ : Энергоатомиздат, 1989^ -575 с

5^ Харари, Ф. Теория графов / Ф^ Хараре - М^ : Мир, 1973^ - 301 с 6^ Shimansky, G. A. TRANS_MU computer code for computation of transmutant formation kinetics in advanced structural materials for fusion reactors /

A^ Shimansky, - J NucL Mater - 1999^ - № 271-272• - P^ 30-34^

7^ Белозерова, А. Р. Оптимизация схемы расчета трансмутации методом ветвей и границ / А^ Р^ Белозерова, Г А^ Шиманский // Физическое моделирование изменения свойств реакторных материалов в номинальных и аварийных условиях : сб^ рефератов семинара^ - Димитровград : ГНЦ РФ НИИАР, 2005^ - С 75-77^

8^ Мельников, Б. Мультиэвристический подход к задачам дискретной оптимизации / Б^ Мельников // Кибернетика и системный анализ (НАН Украины) - 2006^ -№ 3^ - С 32-42•

9^ Мельников, Б. Параллельная реализация мультиэвристического подхода в задаче сравнения генетических последовательностей / Б^ Мельников, А^ Панин //

Вектор науки Тольяттинского государственного университета^ - 20121 - № 4^ -С 83-86!

References

1 Ma M^ Materialy yadernykh energeticheskikh ustanovok [Materials of nuclear power plants) Moscow: Energoatomizdat, 1987, 408 p^ h Gromkovich Yu Teoreticheskaya informatika. Vvedenie v teoriyu avtomatov, teoriyu vychislimosti, teoriyu slozhnosti, teoriyu algoritmov, randomizatsiyu, teoriyu svyazi i kriptografiyu [Theoretical informatics• Introduction into the theory of automata, the theory of computability, the theory of complexity, the theory of algorithms, randomization, the theory of communication cryptography] • Saint Petersburg: BKhV-Peterburg, 2010, 336^

3^ Hromkovic I Algorithmics for Hard Problems: Introduction to Combinatorial Optimization, Randomization, Approximation, and Heuristics• 2nd Ed^ Springer, 200L 4^ Gerasimov A^ S^ Zaritskaya T S^, Rudik A^ P^ Spravochnikpo obrazovaniyu nuklidov v yadernykh reaktorakh [Reference book on nuclide formation in nuclear reactors) Moscow: Energoatomizdat, 1989, 575 p^

5^ Kharari Teoriya grafov [Graph theory) Moscow: Mir, 1973, 301 p^

6^ Shimansky G^ A^ J. Nucl. Mater. 1999, no^ 271-272, pp^ 30-34^

7^ Belozerova A^ R^, Shimanskiy G^ A^ Fizicheskoe modelirovanie izmeneniya svoystv reaktornykh materialov v nominal’nykh i avariynykh usloviyakh: sb. referatov seminara [Physical simulation of reactor material properties changes in nominal emergency conditions: collected reports of the seminar) Dimitrovgrad: GNTs RF NIIAR, 2005, pp^ 75-77^

8^ Mel'nikov B^ Kibernetika i sistemnyy analiz (NAN Ukrainy) [Cybernetics and system analysis (National Academy of Sciences of Ukraine)) 2006, no^ 3, pp^ 32-42.

9^ Mel'nikov B^ Panin A^ Vektor nauki Tol’yattinskogo gosudarstvennogo universiteta [Science vector of Togliatti State University) 2012, no^ 4, pp^ 83-86^

Белозерова Алла Равильевна

кандидат физико-математических наук, старший научный сотрудник, Государственный научный центр Научно-исследовательского института атомных реакторов (Россия, Ульяновская область, г. Димитровград-10)

E-mail: lmni@niiar.ru

Мельников Борис Феликсович

доктор физико-математических наук, профессор, кафедра прикладной математики и информатики, Тольяттинский филиал Самарского государственного университета (Россия, г. Тольятти, ул. Юбилейная, 31Г)

E-mail: B.Melnikov@tltsu.ru

Belozerova Alla Ravil'evna Candidate of physical and mathematical sciences, senior staff scientist, State Scientific Center of the Research Institute of Nuclear Reactors (10-Dimitrovgrad, Ulyanovsk region, Russia)

Mel'nikov Boris Feliksovich Doctor of physical and mathematical sciences, professor, sub-department of applied mathematics and informatics, Togliatti branch of Samara State University (31g Yubileynaya street, Togliatti, Russia)

УДК 519.176, 519.622.2, 519.688 Белозерова, А. Р.

Подход к математическому моделированию трансмутационных процессов в ядерных энергетических установках / А. Р. Белозерова, Б. Ф. Мельников // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2014. - № 2 (30). - С. 88-100.

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