Моделирование распространения фотонов с помощью метода Монте-Карло

Пусть

Тогда

Уравнение условно-критического реактора

k – эффективный коэффициент размножения

Замедление нейтронов при рассеянии

E

Ej

Ei

Ei+1

Ei-1

источник мгновенных нейтронов деления;

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

Сеточная модель реактора в R-Z геометрии

Рассмотрим теперь член описывающий утечку -

.
Согласно теореме Остроградского можно написать,

, где

Конечно-разностная аппроксимация указанных выше интегралов будет выглядеть:

Поток нейтронов на границе

определяется исходя из условия равенства токов

Граничные условия

ФNR,z=0, z=1,….,NZ

Фr,0=0, r=1,….,NR

Фr,NZ=0, r=1,….,NR

Нормировка потока нейтронов