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

UptoLike

где
γ
случайная величина на интервале (0,1).
Остается еще выяснить, как выбирать случайное направление
нейтрона после рассеяния. Так как задача симметрична относительно
оси х, то направление вполне определяется одним углом
ϕ
между
направлением скорости нейтрона и осью Ох. Можно доказать, что
требование равной вероятности любого направления в этом случае
равносильно требованию, чтобы косинус этого угла
μ
= cos
ϕ
был
равномерно распределен в интервале (-1,1) и
μ
=2
γ
1 .
2.3.3. Схема моделирования траекторий нейтронов
Предположим, что нейтрон испытал k-e рассеяние внутри
пластинки в точке с абсциссой x
k
и поcле этого начал двигаться в
направлении
μ
k
. Разыграем длину свободного пробега
λ
= –(1/
Σ
)ln
γ
и вычислим абсциссу следующего столкновения (рис. 2.6)
x
k+1
= x
k
+
λ
k
μ
k
.
Проверим условие прохождения сквозь пластинку:
x
k+1
>h.
Если это условие выполнено, то счет
траектории нейтрона заканчивается и
добавляется единица к счетчику прошедших
частиц. В противном случае проверяем
условие отражения:
x
k+1
<0.
Если это условие выполнено, то счет
траектории заканчивается и добавляется
единица к счетчику отраженных частиц.
Если же и это условие не выполнено, т. е. 0
x
k+1
h, то, значит,
нейтрон испытал k+1-е столкновение внутри пластинки, и надо разыг-
рать «судьбу» нейтрона при столкновении.
Рис.2.6
Выбираем очередное значение
γ и проверяем условие
поглощения:
c
γ
<
ΣΣ
Если последнее неравенство выполнено, то счет траектории
заканчивается и добавляется единица к счетчику поглощенных частиц.
В противном случае мы считаем, что нейтрон испытал рассеяние в
36