Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Монте-Карло.docx
Скачиваний:
7
Добавлен:
23.08.2019
Размер:
50.81 Кб
Скачать

Терминология

Для подробного объяснения принципа работы метода Монте-Карло необходимо ввести некоторые термины:

  • Зона: геометрическая ячейка произвольной формы и сложности, состоящая из единственного материала. Из конечного числа зон и плоскостей, их ограничивающих, состоит вся система

  • Материал: система состоит из зон, зоны состоят из материалов. Каждый материал характеризуется материальным номером, температурой и нуклидным составом

  • Вес зоны ω: каждой зоне соответствует свой вес. Он описывает, насколько значима конкретная зона. Веса можно использовать для понижения погрешности, чтобы точнее рассчитать важные зоны

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

  • История: «путь» частицы, от источника, до ее поглощения (или для частицы с малым весом разыгрывается русская рулетка ). Для понижения погрешности, одна история может записываться для нескольких частиц (используя веса). В программе MCU используется несколько другое понятие истории – каждая промоделированная частица входит в поколение, поколения частиц объединяются в серию, а серии частиц объединяются в историю.

2) Выбор взаимодействия

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

Вероятность взаимодействия:

Вероятность частицы провзаимодействовать в среде с полным сечением (E) в пути от x0 к x1 равна

(2.1)

Случайная величина определяет, произойдет ли взаимодействие ( , или не произойдет ( Если p заменить на r, а |x0-x1| заменить на , то

(2.2)

Так как распределение величины (1-r) аналогично распределению величины r, то можно заменить (1-r) на r, чтобы получить выражение

(2.3)

Взаимодействие происходит при . Если это условие выполняется, то есть расстояние до точки взаимодействия.

Выбор изотопа

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

(2.4)

Макросечения всех реакций складываются в полное макросечение. Примеры составляющих полного макросечения можно найти в таблице 2.1

Символ

=

Вид сечения

=

полное сечение взаимодействия

=

упругое рассеяние

=

=

потенциальное рассеяние

=

резонансное рассеяние

=

интерференционное рассеяние

=

неупругое рассеяние

=

полное сечение рассеяния

=

=

поглощение

=

=

испускание гамма луча: (n,γ) реакция

=

испускание протона: (n,p) реакция

=

испускание альфа-частицы: (n,α) реакция

=

испускание двух нейтронов : (n,2n) реакция

=

деление

=

неупругое столкновение (сечение реакции)

=

таблица 2.1: выборка сечений для нейтронов

Выбор реакции, происходит снова через случайную величину. Если отдельные взаимодействия последовательно пронумеровать с индексом i, то j-ое взаимодействие выбирается по принципу:

(2.5)