вии с физическим смыслом в конкретной задаче теории замедления нейтронов:
1)условие на внешней выпуклой экстраполированной границе среды с вакуумом для конечных или полубесконечных сред;
2)условия сшивки решений на границе раздела сред с различными физическими свойствами;
3)начальное условие или условие внешнего источника и т.д.
I. Условие на внешней выпуклой экстраполированной границе среды с вакуумом для конечных или полубесконечных сред. По аналогии с условием (3.30) Ф(rэ,t) = 0 , налагаемым на
решение односкоростного уравнения диффузии и рассмотренным ранее, это можно сформулировать как равенство нулю плотности замедления на экстраполированной границе:
j(rэ,τ) = 0 , |
(5.25) |
где rэ – экстраполированная граница среды с вакуумом.
II. Условия сшивки решений на границе раздела сред с различными нейтронно-физическими свойствами. По аналогии с условиями, налагаемыми на решение односкоростного уравнения диффузии и рассмотренными ранее, и на их базе это можно сформулировать как равенство плотностей замедления и токов замедления.
Из условия (3.31) Ф1 (r0 ,t) = Ф2 (rG0 ,t) с учетом (5.8) перехода от
плотности потока замедляющихся нейтронов Ф к плотности замедления следует j имеем
|
|
|
|
|
|
|
|
|
|
j (r ,τ) |
= |
|
j |
2 |
(rG,τ) |
|
. |
|
|
|
|
|
|
(5.26) |
|||||||||||||
|
|
|
|
|
|
|
1 |
0 |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
( |
ξΣs ) |
|
( |
|
|
|
|
) |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
ξΣs |
|
|
|
|
|
|
||||||||||||||||||||||||
Из условия (3.32) |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
G |
с учетом (5.8) сле- |
||||||||||
D Ф (r ,t) = D Ф |
2 |
(r ,t) |
|||||||||||||||||||||||||||||||||||
дует, что |
|
|
|
|
|
1 |
|
|
1 |
0 |
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
0 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
D |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
n ( j1 (rG0 ,τ)) = |
|
|
|
|
|
n ( j2 (rG0 ,τ)) |
(5.27) |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
ξΣs |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ξΣs 2 |
|
|
|
|
|
||||||||||||
или с учетом D = |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
3Σtr |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
n j(r,τ) |
|
|
|
|
|
= |
|
n |
j(r,τ) |
|
|
|
/ |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
(5.27 ) |
||||||
|
|
|
|
|
|
|
|
|
ξΣsΣtr |
|
1 |
|
r =r0 |
|
|
|
|
|
|
ξΣsΣtr |
2 |
|
r =r0 |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
136
III. Начальное условие или условие внешнего источника.
Данное условие задает распределение источников быстрых нейтронов, возраст которых равен нулю
j(r ,0) = Sδ(rG)δ(τ) |
(5.28) |
или |
|
j(r ,0) = S(rG) , |
(5.28/) |
где j(rG,0) – число нейтронов, испускаемых источником в единице
объема в единицу времени и начинающих замедляться; S – мощность источников нейтронов; δ(r ) , δ(τ) – дельта-функции; S(r ) –
распределение внешних источников быстрых нейтронов.
Часто уравнение возраста записывают в элементарной форме с использованием вероятности избежать поглощения нейтрона в процессе замедления в среде.
Если рассматривать поглощающую и непоглощающую среды, то уравнение возраста (5.12) соответственно можно записать в виде
j |
(rG,τ) − |
1 |
j |
(rG,τ) − |
∂jпогл (r ,τ) |
= 0 |
(5.29) |
|
L2 |
||||||||
погл |
|
погл |
|
∂τ |
|
|
||
|
|
τ |
|
|
|
|
|
и элементарного уравнения возраста (уравнения возраста Ферми)
jнепогл (rG,τ) − ∂jнепогл (rG,τ) = 0
∂τ
c начальными условиями (условиями внешнего источника)
jпогл (r ,0) = jнепогл (rG,0) = S(rG) .
Связь между плотностями замедления в поглощающей и непоглощающей средах можно записать через вероятность избежать
поглощения нейтрона в процессе |
замедления в среде ϕ(τ) : |
|
||
j |
(r ,0) = j |
|
(rG,0)ϕ(τ) . |
(5.31) |
погл |
непогл |
|
|
|
Подставляя функцию (5.31) в ее уравнение (5.29), получим уравнение
ϕ(τ) j |
|
(rG,τ) − |
1 |
ϕ(τ) j |
(rG,τ) − |
|
|
|
L2 |
|
|||||
непогл |
|
|
непогл |
|
|
||
|
|
|
τ |
|
|
|
|
−ϕ(τ) ∂jнепогл (rG,τ) −∂j |
|
(rG,τ) ∂ϕ(τ) = |
0, |
||||
|
∂τ |
непогл |
|
∂τ |
|
||
|
|
|
|
|
|||
которое перепишем в виде
137
G |
∂jнепогл (rG,τ) |
G |
d [ln ϕ(τ)] |
|
1 |
|
= 0 . (5.32) |
||
jнепогл (r ,τ) − |
|
− jнепогл (r ,τ) |
|
|
+ |
|
|
||
∂τ |
dτ |
2 |
|||||||
|
|
|
|
L |
|
|
|||
|
|
|
|
|
τ |
|
|||
=0 |
!!! |
|
|
|
|
|
|
|
|
Так как первые два слагаемых есть элементарное уравнение возраста (5.29/), равное нулю, то и скобка равна нулю. Отсюда получаем уравнение относительно ϕ(τ)
d [ln ϕ(τ)] |
= − |
1 |
, |
(5.33) |
dτ |
2 |
|||
|
L |
|
|
|
|
|
τ |
|
|
для которого начальное условие с учетом (5.30) определяется как
ϕ0 = ϕ(0) =1. |
(5.34) |
Физический смысл условия (5.31): возраст родившихся, быстрых нейтронов деления τ = 0, в момент рождения, т.е. до начала замедления в среде вероятность избежать поглощения таких нейтронов в среде ϕ0 = ϕ(τ = 0) =1 – достоверное событие.
Решая уравнение (5.30), можно получить простейшее дифференциальное уравнение с разделяющимися переменными
1
ϕ(τ)
ln ϕ
ϕ0
dϕ(τ) |
|
|
1 |
, |
ϕ |
|
d |
ϕ(τ) |
|
|
|
τ |
d |
τ |
|
|
||||
= − |
∫ |
|
= −∫ |
|
, |
|||||||||||||||
|
2 |
|
ϕ(τ) |
|
2 |
|||||||||||||||
dτ |
|
|
|
L |
ϕ0 |
|
|
|
|
0 |
L |
|
|
|
||||||
|
|
|
|
τ |
|
|
|
|
|
|
|
|
τ |
|
||||||
= −τ |
dτ |
, ϕ(τ) = ϕ |
|
exp |
|
− |
τ |
|
dτ |
|
, |
|||||||||
2 |
|
|
∫ |
|
|
2 |
|
|
||||||||||||
∫ |
|
|
|
|
N0 |
|
|
|
|
|
|
|||||||||
0 |
Lτ |
|
|
|
|
=1 |
|
|
|
|
0 |
|
Lτ |
|
|
|
||||
|
τ |
dτ |
|
|
ϕ(τ) = exp |
−∫ |
|
. |
(5.35) |
L2 |
||||
|
0 |
τ |
|
|
Подставляя для интеграла в экспоненте (5.35) выражение (5.11) с учетом перехода (*) равенства (5.14) (из (4.48)
|
|
|
|
|
|
|
|
|
|
|
|
|
E0 |
|
|
|
|
|
dE |
|
|
|
|
|||||
|
|
|
|
du = d ln |
|
|
|
|
= − |
|
|
), |
|
|
|
|||||||||||||
|
|
|
|
E |
E |
|
|
|
||||||||||||||||||||
имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
d |
τ |
|
|
|
|
|
|
du |
|
|
|
|
|
|
|
|
|
1 |
|
|
dE |
|
|||||
|
|
|
= |
|
|
|
|
|
|
= |
|
|
|
|
|
|
, |
|||||||||||
|
|
|
|
2 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
E |
|||||||
|
|
|
|
|
3ξΣ |
Σ |
|
|
3ξΣ |
Σ |
|
|||||||||||||||||
|
|
L |
|
|
tr |
L |
|
|
L |
|
|
|||||||||||||||||
|
|
|
|
τ |
|
|
|
|
|
s |
|
|
τ |
|
|
|
|
|
s |
|
tr τ |
|
|
|
||||
откуда с учетом |
D = |
|
|
1 |
|
|
и |
|
1 |
|
= |
Σa |
|
|
получим |
|
|
|
||||||||||
|
3Σ |
|
|
2 |
|
|
|
|
|
|||||||||||||||||||
|
|
|
tr |
|
|
|
L |
|
|
|
D |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
τ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
138
|
dτ |
= |
1 |
|
|
|
dE |
= |
|
|
Σa |
|
dE |
. |
|
||||||||||
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
3ξΣ |
Σ |
|
|
|
E |
|
|
|
|
E |
|
|||||||||||||
|
L |
tr |
L |
|
|
|
|
|
|
|
ξΣ |
s |
|
|
|||||||||||
|
τ |
|
|
s |
|
τ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Окончательно получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
E0 |
|
Σa |
|
|
dE |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
∫ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
ϕ(E0 , E) = exp |
− |
|
ξΣs E |
|
. |
(5.36) |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
E |
|
|
|
||||||||||
Это выражение для слабопоглощающих сред ( Σа (E) << Σs (E) ,
Σ(E) = Σs (E) + Σа (E) ≈ Σs (E) ) совпадает с выражениями (4.76), (4.79).
6. НЕЙТРОННО-ФИЗИЧЕСКИЕ ОСНОВЫ РАСЧЕТА ЯДЕРНЫХ РЕАКТОРОВ
6.1. Нейтронный цикл в ядерном реакторе
Цепная ядерная реакция (ЦР) – это реакция деления ядер урана и плутония нейтронами, рождающимися при делении. Необходимое условие протекания самоподдерживающейся ЦР – рождение в каждом очередном акте деления не менее одного нейтрона. Достаточность этого условия зависит от исхода конкуренции четырех процессов, имеющих место в размножающей среде (среде, содержащей нуклиды ядерного топлива):
-деление ядер топлива;
-радиационный захват в топливе;
-радиационный захват нейтронов конструкционными материалами активной зоны и неделящимися компонентами топлива;
-утечка (вылет) нейтронов за пределы размножающей среды.
ВЯР должен поддерживаться баланс скоростей генерации и поглощения, утечки нейтронов. В общем случае скорость генерации тепловых нейтронов может отличаться от скорости поглощения и утечки.
Нейтронным циклом (рис. 29) называется совокупность процессов, приводящих к появлению нового поколения быстрых нейтронов (БН) в реакторе (и, следовательно, к смене поколений при контролируемой цепной реакции деления, сопровождающейся вы-
139
делением ≈200 МэВ при каждом акте деления ядра U-235 ТН, получаемых из БН деления путем замедления через реакции рассеяния). Наиболее эффективно U-235 и Pu-239 делятся ТН.
Рассмотрим следующую цепочку
(1) |
(2) |
(3) |
(4) |
(5) |
|
nбi →nбi ε→nбi ε рз →nбi ε рз ϕ→nбi ε рз ϕ рд → |
(6.1) |
||||
|
|
|
|
η |
|
(5) |
|
(6) |
|
P |
|
→nбi ε рз ϕ рд Θ→nбi ε рз ϕ рд |
Θ f υf = nб(i+1) . |
|
|||
В данных обозначениях:
(1) – рост числа первичных БН за счет деления 238U под действием БН: ε – коэффициент размножения на БН:
ε = полное число БН , полученных при делении нейтронами всех энергий.
(6.2)
В ВВЭР-1000 величина ε ≈ 1,03. Характерные значения в ЯР ε = = (1,02–1,05) < 1,06, и его практически нельзя изменить, так как природный уран – это в основном 238U. Увеличение плотности потока нейтронов в ЯР на твердотельном топливе за счет деления на БН составляет около 2–3%.
(2) – участие в делении топлива принимает часть БН, избежавших утечки в процессе замедления до резонансных энергий, имеющихся в зависимости сечения радиационного захвата урана238: рз – вероятность избежать утечки нейтронов в процессе замедления до эпитепловых нейтронов
р = число БН поколения, избежавших утечки при замедлении. (6.3)
з |
полное число БН поколения, начавших замедление |
|
В ВВЭР-1000 pз ≈ 0,97. Величина pз – изменяющийся параметр, например, за счет изменения температуры и давления замедлителя.
Из решения уравнения возраста Ферми можно получить
pз = exp(–χ2τт), (6.4)
где χ2 – материальный параметр ЯР, τт – возраст тепловых нейтронов.
(3) – участие в делении топлива принимает часть нейтронов, избежавших резонансного поглощения ядрами 238U в конце процесса замедления до ТН: φ – вероятность избежать резонансного захвата нейтронов в конце замедления
140