Компьютерное моделирование физических явлений. Малютин В.М - 37 стр.

UptoLike

точке с абсциссой x
k+1
. Тогда разыгрываем новое направление скорости
нейтрона
1
21
k
μ
γ
+
=
и затем повторяем весь цикл снова (но, конечно, уже с другими
значениями
γ).
Все
γ написаны без индексов, так как имеется в виду, что каждое
значение
γ используется всего один раз. Для расчета одного звена
траектории нужны три значения
γ. Начальные значения для каждой
траектории:
00
0, 1x
μ
=
=
.
После того как будут сосчитаны N траекторий, окажется, что N
+
нейтронов прошли сквозь пластинку, N
нейтронов отразились от нее, а
N
0
нейтронов были поглощены в ней. Очевидно, искомые вероятности
приближенно равны отношениям:
0
0
,,
NN
ppp
NN
+−
+−
≈≈.
N
N
На рис. 2.7 приведена блок-схема программы для расчета этой
задачи. Индекс j это номер траектории, индекс k — номер
столкновения (вдоль траектории).
Отметим только, что метод Монте-Карло позволяет решать
гораздо более сложные задачи об элементарных частицах: исследуемая
среда может состоять из различных веществ и иметь любую
геометрическую структуру; энергия частицы при каждом столкновении
может меняться. Можно учитывать много других ядерных процессов
(например, возможность деления атома при столкновении с нейтроном
и образование при этом новых нейтронов). Можно рассчитать условие
возникновения и поддержания цепной реакции и т. д.
37