Волков / Конспекты (в LaTeX) / 17_VolkovYN_TPN
.pdfНациональный исследовательский ядерный университет
«МИФИ»
Институт ядерной физики и технологий
Конспект лекции №17
Основы замедления нейтронов
Преподаватель: Волков Ю.Н. Группа: ТПН Дата: 19.11.2025
Москва, 2025
Содержание
1 |
Введение |
2 |
2 |
Средний косинус угла рассеяния |
2 |
3 |
Уравнение замедления |
3 |
4 |
Плотность замедления |
5 |
5 |
Замедление в непоглощающей среде на водороде |
6 |
1
1Введение
[00:01] На предыдущей лекции мы рассмотрели основные характеристики замедления нейтронов. Были введены такие понятия как коэффициент замедления и замедляющая способность, которые показывают эффективность того или иного замедлителя. Хороший замедлитель должен не только эффективно замедлять нейтроны, но и обладать малым сечением поглощения.
[01:35] Из сравнительного анализа следует, что тяжелая вода является наилучшим замедлителем по соотношению замедляющей способности и поглощения. Кроме того, мы рассмотрели закон упругого рассеяния, который описывает распределение энергии нейтрона после столкновения с ядром.
[01:57] В рамках ступеньки замедления мы установили, что плотность вероятности нахождения энергии нейтрона после столкновения в интервале от αE1 до E1 является константой. Средняя энергия после столкновения находится ровно посередине этого интервала, а средняя потеря энергии составляет половину этого интервала.
[02:11] Мы также ввели понятия целогарифмической потери энергии и литаргии. Однако один важный аспект остался нерассмотренным — средний косинус угла рассеяния.
2Средний косинус угла рассеяния
[02:57] Продолжим рассмотрение четвертого параграфа, который относится к третьему разделу. Этот параграф называется "Средний косинус угла рассеяния".
[03:28] Рассмотрим сложение векторов скоростей va + vцм = v2, где θ — угол между вектором скорости центра масс и вектором скорости нейтрона после столкновения, а ψ — угол между вектором скорости центра масс и вектором относительной скорости.
[04:24] На диаграмме скоростей:
•θ — угол рассеяния в лабораторной системе отсчета (ЛСО)
•ψ — угол рассеяния в системе центра масс (СЦМ)
[05:27] Найдем выражение для энергии нейтрона после столкновения:
E2 |
= E1 |
1 + A2 + 2A cos ψ |
|
(1 + A)2 |
|||
|
|
где A — массовое число ядра-мишени.
[06:46] Плотность вероятности рассеяния по углу ψ в СЦМ:
P (ψ) = 12 sin ψ
[07:07] Рассчитаем средний косинус угла рассеяния в СЦМ:
(2.1)
(2.2)
cos¯ ψ = ˆ π cos ψ |
· |
P (ψ) dψ = |
1 |
ˆ π cos ψ sin ψ dψ = 0 |
(2.3) |
|
2 |
||||||
0 |
|
0 |
|
[09:25] Равенство нулю среднего косинуса угла рассеяния в СЦМ означает, что рассеяние изотропно в этой системе отсчета. Это подтверждает наше предположение о том, что в среднем нейтрон рассеивается под углом π/2.
2
[10:23] Более интересен средний косинус угла рассеяния в лабораторной системе отсчета, который обозначается буквой µ с крышечкой. Эта величина важна при рассмотрении диффузии нейтронов и используется при введении транспортного сечения.
[11:39] Для нахождения среднего косинуса угла рассеяния в ЛСО рассмотрим геометрические соотношения на диаграмме скоростей. Длина отрезка OB равна:
OB = v2 cos θ = vцм + va cos ψ = |
1 + A cos ψ |
v1 |
(2.4) |
|||
|
|
|||||
|
|
|
1 + A |
|
||
[13:33] Выразим косинус угла θ: |
|
|
|
|
||
cos θ = |
|
1 + A cos ψ |
|
|
|
(2.5) |
|
|
|
||||
p |
|
|||||
1 + A2 + 2A cos ψ |
|
|
||||
[15:13] Таким образом, косинус угла рассеяния в ЛСО является функцией угла рассеяния в СЦМ.
[16:09] Для расчета среднего косинуса угла рассеяния в ЛСО необходимо вычислить интеграл:
µ¯ = ˆ π cos θ |
· |
P (ψ) dψ = |
1 |
ˆ π |
|
1 + A cos ψ |
sin ψ dψ |
(2.6) |
2 |
|
|
||||||
2 |
||||||||
0 |
|
0 |
p1 + A + 2A cos ψ |
|
|
|||
[17:23] Этот интеграл является табличным. Его вычисление дает следующий результат:
µ¯ = |
2 |
|
(2.7) |
|
3A |
||||
|
|
|||
[19:37] Этот результат используется при рассмотрении диффузии нейтронов, где средний косинус угла рассеяния играет важную роль.
3Уравнение замедления
[20:30] Пятый параграф называется "Уравнение замедления". Уравнение замедления представляет собой уравнение баланса нейтронов.
[21:34] Рассмотрим случай, когда поток нейтронов зависит только от энер-
Φ(r, E, Ω, t)
гии, то есть Φ(E). Для этого необходимо сделать следующие допущения:
•Среда бесконечная, гомогенная и неразмножающая
•Источник равномерно распределен по среде, изотропный
•Мощность источника q [н/(с·см³)]
•Источник испускает нейтроны с детерминированной энергией E0
•Задача стационарная
[25:40] Запишем уравнение баланса нейтронов на интервале энергий (E, E + ∆E): увод из ∆E равен приходу в ∆E.
3
Рис.:
[27:40] Причины увода нейтронов из интервала ∆E: |
|
|
|
|
|
• Поглощение A(E) |
|
|
|
|
|
• Замедление Rs(E) (рассеяние с потерей энергии) |
|
|
|
|
|
[29:00] Причины прихода нейтронов в интервал ∆E: |
|
|
|
|
|
• Рассеяние сверху Ps(E) |
|
|
|
|
|
• Источник Q(E) |
|
|
|
|
|
[31:10] Выражения для различных процессов: |
|
|
|
|
|
A(E) = Σa(E)Φ(E)∆E |
|
|
|
(3.1) |
|
Rs(E) = Σs(E)Φ(E)∆E |
|
|
|
(3.2) |
|
Q(E) = qδ(E − E0)∆E |
|
|
|
(3.3) |
|
[33:00] Рассеяние сверху Ps(E) вычисляется следующим образом: |
|
||||
Ps(E) = ˆ min(E0;E/α) Σs(E′)Φ(E′)dE′ |
· |
∆E |
|
(3.4) |
|
E′(1 − α) |
|||||
E |
|
||||
где первый множитель — количество нейтронов, рассеявшихся в интервале dE′, а второй множитель — вероятность того, что после рассеяния нейтрон попал в интервал ∆E.
[39:39] Уравнение замедления в общем виде:
Σt(E)Φ(E) = ˆ min(E0;E/α) Σs(E′)Φ(E′) |
dE′ |
|
+ qδ(E |
− |
E0) |
(3.5) |
|
E′(1 − α) |
|||||||
E |
|
|
|
||||
[45:55] Введем следующие обозначения:
• Σt(E)Φ(E) = F (E) — плотность столкновений
• Σs(E)Φ(E) = Fs(E) — плотность рассеяния
4
4Плотность замедления
[47:15] 6 параграф – "Плотность замедления". Плотность замедления j(E) определяется как среднее число нейтронов за единицу времени, пересекающих данное значение энергии в процессе замедления.
[48:15] Физически это означает, что на оси энергий мы рассматриваем энергию E и считаем, сколько нейтронов в единицу времени пересекает это значение при движении слева направо (уменьшение энергии).
Рис.:
[49:40] Математическое выражение для плотности замедления:
j(E) = ˆ min(E0;E/α) Σs(E′)Φ(E′)dE′ |
· |
E − αE′ |
(4.1) |
|
E′(1 − α) |
||||
E |
|
[52:30] Важное свойство: если среда непоглощающая (Σa = 0), то плотность замедления равна мощности источника:
j(E) = q |
(4.2) |
[53:20] Если среда поглощающая (Σa ̸= 0), то отношение плотностей замедления на разных энергиях:
ϕ(E1 → E2) = |
j(E2) |
(4.3) |
j(E1) |
называется вероятностью избежать поглощения при замедлении от E1 до E2.
5
Рис.:
[53:55] Эта вероятность также называется вероятностью избежать резонансного поглощения, так как основное поглощение при замедлении происходит на резонансах урана-238. В тепловом реакторе эта вероятность составляет примерно 0.75, то есть около 25
5Замедление в непоглощающей среде на водороде
[57:45] Седьмой параграф рассматривает замедление в непоглощающей среде, состоящей из ядер водорода (A = 1).
[58:50] Уравнение замедления для этого случая:
Fs(E) = ˆ E0 Fs(E′) |
dE′ |
+ qδ(E |
− |
E0) |
|
E′ |
|||||
E |
|
|
[1:00:03] Решение этого уравнения ищем в виде суммы двух функций:
˜ −
Fs(E) = Fs(E) + qδ(E E0)
(5.1)
(5.2)
где ˜ — плотность рассеяния нейтронов, испытавших хотя бы одно рассеяние, а −
Fs(E) qδ(E E0) — плотность рассеяния нейтронов источника, которые не испытали столкновений.
[1:02:10] Подставляя это решение в уравнение замедления, получаем:
F˜s(E) = ˆ E0 F˜s(E′) |
dE′ |
+ |
q |
|||
E′ |
E0 |
|||||
E |
|
|
|
|||
[1:04:00] Решение этого интегрального уравнения: |
|
|
||||
˜ |
|
|
|
|
|
|
ln Fs(E) = −ln E + C |
|
|||||
или |
C |
|
|
|
||
˜ |
|
|
||||
Fs(E) = |
E |
|
|
|||
(5.3)
(5.4)
(5.5)
6
[1:06:48] Постоянную C находим из условия при E = E0:
˜ |
C |
|
q |
C = q |
(5.6) |
Fs(E0) = |
E0 |
= |
E0 |
[1:07:35] Таким образом, плотность рассеяния:
Fs(E) = |
q |
+ qδ(E − E0) |
(5.7) |
E |
[1:07:40] В области энергий, отличных от E0, спектр нейтронов при замедлении на водороде в непоглощающей среде описывается зависимостью 1/E. Этот спектр называется спектром Ферми.
[1:09:14] Рассмотрим этот спектр в переменных литаргии. Литаргия определяется как:
u = ln |
E0 |
(5.8) |
|
E |
|||
|
|
Рис.:
[1:10:40] Переход от переменной энергии к переменной литаргии осуществляется по формуле:
|
|
|
|
Φ(u) = Φ(E) |
|
dE |
|
(5.9) |
|
|
|
|
du |
||||
|
|
|
|
|
|
|
|
|
где |
dE |
|
= E. |
|
|
|
|
|
du |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
[1:12:53] |
Таким образом, в переменных литаргии спектр Ферми является константой: |
||||||
|
|
|
|
Fs(u) = Fs(E)E = const |
(5.10) |
|||
[1:14:13] Проверим, что в непоглощающей среде плотность замедления равна мощности источника: Σa = 0, пусть также среда бесконечная (Будем доказывать, что действительно j(E) = q), источники равномерно распределены и т.д. (j ̸= j(E, r)). Тогда:
E0 |
Fs(E′)dE′ E′ |
E0 |
(E′)2 dE′ + |
E0 |
= qE |
E |
− E0 |
|
+ |
|
E0 |
= q (5.11) |
|
j(E) = ˆE |
= E ˆE |
|
|||||||||||
|
|
E |
|
q |
qE |
|
1 |
1 |
|
|
qE |
|
|
что и требовалось доказать.
7
Рис.:
8
