|
число нейтронов эпитепловых энергий, |
|
|
ϕ = |
избежавших резонансного захвата в 238 U |
. |
(6.5) |
|
|||
|
все нейтроны поколения эпитепловых энергий |
|
|
В ВВЭР-1000 φ ≈ 0,71. Величина φ изменяющийся параметр, влияющий на ядерную безопасность: рост мощности ЯР обусловлен ростом числа делений в активной зоне ЯР и ростом температуры топлива. Вследствие эффекта Доплера растет сечение радиационного захвата эпитепловых нейтронов 238U (и вероятность их захвата 1 – φ). В результате поток ТН в активной зоне ЯР уменьшается, и ЯР саморегулируется – снижает мощность.
Нахождение 238U в ЯР неизбежно в силу того, что это низкообогащенный, природный уран, но в каждой бочке меда есть своя ложка дегтя, и наоборот. Резонансное поглощение, в общем случае, не полностью вредно, так как при этом из 238U может образоваться в меньших количествах вторичное ядерное топливо – делящийся изотоп 239Pu, далее – с уменьшающейся вероятностью пороговый нуклид 240Pu и из него делящийся нуклид 241Pu:
β− |
β− |
α |
23892 U + 01n → 23992 U* → |
23993 Np → 23994 Pu |
→ , |
23,5 мин |
2,3 сут |
2,4 104 лет |
23994 Pu + 01n → 24094 Pu + 01n → 24194 Pu .
Кроме того, резонансное поглощение на 238U за счёт доплерэффекта даёт замечательную способность активной зоны ЯР типа ВВЭР саморегулироваться – снижать мощность на первоначальном этапе, т.е. при росте мощности повышает безопасность (см. выше). С другой стороны, при снижении мощности необходимо будет компенсировать высвобождающуюся реактивность, чтобы обеспечить ядерную безопасность.
(4) – участие в делении топлива принимает часть нейтронов, избежавших утечки при диффузии ТН до момента их поглощения топливом: рд – вероятность избежать утечки ТН в процессе диффузии
рд = число ТН поколения, избежавших утечки при диффузии. (6.6)
В ВВЭР-1000 pд ≈ 0,99. Величина pд изменяющийся параметр, например, за счет изменения температуры и давления замедлителя.
141
Из решения волнового уравнения Гельмгольца следует развер-
нутое условие критичности реактора |
|
|
|
|
||||
|
k |
∞ |
exp(−χ2τ |
т |
) |
=1, |
(6.7) |
|
|
|
|
1+ B |
2 2 |
|
|
||
|
|
|
L |
|
|
|
|
|
С учетом (6.4) pз = exp(-χ2τт) и |
г |
|
|
|
|
|||
|
|
|
|
|
|
|||
рд = (1 + Вг2 L2)-1, |
(6.8) |
|||||||
где L – длина диффузии, τт – возраст тепловых нейтронов, χ2 – ма- |
||||||||
териальный параметр ЯР, Вг2 – геометрический параметр ЯР, |
k∞ – |
|||||||
эффективный коэффициент размножения нейтронов в бесконечной среде (решение уравнения Гельмгольца и перечисленные параметры ЯР рассмотрим далее в этой главе).
(5) – участие в делении топлива продолжает часть нейтронов, которая поглощается именно топливной композицией: Θ – вероятность поглощения ТН топливом (238U и после начала кампании 239Pu) или коэффициент использования ТН
Θ = все ТН поколения, поглощенные топливом и замедлителем. (6.9)
В ВВЭР-1000 Θ ≈ 0,79. Характерные значения в ЯР 0,8–0,9. Величина Θ изменяющийся параметр, влияющий на ядерную безопасность через саморегулируемость ЯР (изменение свойств сред активной зоны ЯР с изменением мощности ЯР и температуры топлива, т.е. температуры, давления и расхода теплоносителя, содержания в нем поглощающих добавок и т.д.).
(6) – появляется (i + 1)-поколение БН: f = |
Σf |
– вероятность |
|
||
|
Σa |
|
того, что поглощение ядром ТН закончится делением ядра (функция компоновки активной зоны ЯР).
νf – среднее число вторичных нейтронов, испускаемых каждым делящимся ядром под действием ТН. Для каждого нуклида это параметр, зависящий от энергии нейтрона. Экспериментально установленные зависимости имеют вид:
для ТН:
νf5(Е) = 2,416 + 0,1337Е [6]; νf5(Е) = 2,862 + 0,1357E [6];
142
для БН (Еп=1.1 МэВ):
νf5(Е) = 2,409 + 0,1389E [6].
Параметр
η= fv |
|
= ν |
|
Σf |
, |
(6.10) |
|
|
|
||||
|
f |
|
f Σa |
|
||
где η – среднее число БН нового поколения, образовавшихся при захвате одного ТН предыдущего поколения топливной композицией. Параметр η – не константа даже при заданной энергии нейтронов, а функция, зависящая от ядерной концентрации 238U и 239Pu:
|
|
η = |
Σ |
а топл |
= |
ν5 R5f |
+ ν9 R9f |
|
= |
|
|
|
|
|
|
|
|
||||
|
|
59 |
|
Σ |
а |
|
|
R5 |
+ R9 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
a |
a |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
ν5σ5f + ν9σ9f |
|
N |
9 |
|
|
||||||
|
5 |
|
9 |
N9Ф |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||
= |
ν5σf N5 |
Ф+ ν9σf |
= |
|
|
|
|
|
N5 |
. |
(6.11) |
||||||||||
σ5 N |
Ф+ σ9 N |
Ф |
|
|
|
|
|
||||||||||||||
|
|
|
|
N |
9 |
|
|
|
|
|
|
||||||||||
|
a 5 |
a |
|
9 |
|
|
|
|
|
σ5а + σ9a |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
N5 |
|
|
|
|
|||||
Для чистого изотопа 235U η5 ≈ 2,5, для топлива ВВЭР-1000 η ≈ 1,8. Величиной η практически сложно управлять, но она и не влияет на ядерную безопасность.
Указанная цепочка рассуждений (6.1) приводит к формуле эф-
фективного коэффициента размножения нейтронов |
|
|
|
||||||
nбi ε рз ϕ рд Θ η5 = nб(i+1) , |
|
|
|
||||||
k |
эф |
= |
nб(i+1) |
= ε р |
ϕ р |
Θ η = k |
∞ |
р р |
. (6.12) |
|
|||||||||
|
|
|
з |
д |
5 |
з д |
|
||
|
|
|
nбi |
|
|
|
|
|
|
Формула четырех сомножителей |
|
|
|
|
|
||||
|
|
|
k∞ = ε ϕ Θ η |
|
|
|
(6.13) |
||
– эффективный коэффициент размножения нейтронов в бесконечной активной зоне реактора, рз∞ = рд∞ = 1 и для ВВЭР-1000 k∞ = = 1,03·0,71·0,79·1,8 = 1,04 (это весьма условная цифра) > 1.
Для критического ЯР |
|
kэф = k∞ рз рд =1 |
(6.14) |
– эффективный коэффициент размножения нейтронов, для ВВЭР
рз = 0,97; рд = 0,99, для kэф = 1 необходимо k∞ = 1,04 > 1:
143
<1− ЯР подкритичен;
kэф =1− ЯР критичен на заданном уровне мощности; (6.15)>1− ЯР надкритичен.
Следует отметить, что формулы для каждого из сомножителей (6.1) имеют простой вид только в гомогенном ЯР. Здесь мы рассмотрели коэффициенты формулы 4-х сомножителей ознакомительно, более подробно вернемся к ним позже в рамках курса «Физика ядерного реактора».
В гомогенном ЯР в формуле 4-х сомножителей ε ≈ 1 и k∞ = ϕ Θ η5 . Если в активной зоне ЯР используется чистый изотоп, например, 235U, то φ = 1, и формула 4-х сомножителей сокращается до k∞ = Θ η5 .
6.2. Нейтронно-физические основы теории критических размеров голых реакторов
В гл. 3 было получено стационарное уравнение диффузии для неразмножающей среды (3.37) и рассмотрены его решения для
элементарных геометрий |
1 |
|
S(rG) |
|
|
|
G |
G |
|
|
|
||
Ф(r ) − |
|
Ф(r ) + |
|
= 0 |
, |
(6.16) |
2 |
D |
|||||
|
L |
|
|
|
|
где L = D / Σa – длина диффузии моноэнергетических нейтронов. Также было показано, что стационарное уравнение диффузии
для размножающей среды имеет вид (3.37/) |
|
S(rG) |
|||||||||||||||
G |
1 |
|
|
νf |
|
Σf |
|
|
|
|
G |
|
|||||
Ф(r ) + |
|
|
|
|
|
|
|
|
|
|
− |
1 Ф(r ) = − |
|
|
|||
|
2 |
|
|
Σa |
|
D |
|||||||||||
или |
L |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
S(r ) |
|
|
|
||
G |
|
|
|
|
2 |
|
|
G |
|
|
|
|
|
||||
Ф(r ) |
+χ |
Ф(r ) |
= − |
|
, |
(6.17) |
|||||||||||
D |
|||||||||||||||||
где |
|
|
|
|
|
νf Σf |
|
−1 |
|
|
|
||||||
χ2 |
= |
|
1 |
|
|
(6.18) |
|||||||||||
|
2 |
|
|
|
|
|
|
||||||||||
|
|
|
Σa |
|
|
|
|
|
|
||||||||
|
|
|
|
|
L |
|
|
|
|
|
|
|
|||||
– материальный параметр реактора (по англ. баклинг).
144
Уравнение вида |
|
|
Ф+ B2Ф = 0 |
|
|
|
||||||||
|
|
|
|
|
|
|
(6.19) |
|||||||
|
|
|
|
|
|
|
г |
|
|
|
|
|
|
|
называется волновым уравнением Гельмгольца. |
|
|
||||||||||||
Решение неоднородного дифференциального уравнения (6.17) |
||||||||||||||
ищут как сумму общего решения однородной части уравнения |
||||||||||||||
(6.19), а также любого частного решения неоднородного уравне- |
||||||||||||||
ния. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рассмотрим теперь общую методику решения уравнения (6.17) с |
||||||||||||||
определением материального и геометрического параметров для |
||||||||||||||
голого (без отражателя) критического ЯР с различной формой ак- |
||||||||||||||
тивной зоны: параллелепипед, цилиндр, сфера. |
|
|
||||||||||||
Пусть активная зона голого ЯР имеет форму |
|
|
||||||||||||
параллелепипеда c ребрами а0, b0, c0 (рис. 30). |
|
|
||||||||||||
Исходное уравнение имеет вид волнового |
|
|
||||||||||||
уравнения Гельмгольца |
Ф+χ2Ф = 0 : |
|
|
|
|
|||||||||
∂2Ф + |
∂2Ф |
+ |
∂2Ф |
+ B2Ф = 0, |
|
(6.20) |
|
|
||||||
∂x2 |
∂y2 |
|
∂z2 |
г |
|
|
|
|
|
|
|
|
|
|
где B2 имеет смысл геометрического парамет- |
Рис. 30. Активная зона |
|||||||||||||
г |
|
|
|
|
|
|
|
|
|
|
|
|
реактора в виде |
|
ра ЯР. |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
параллелепипеда |
||
Точку отсчета (начало координат) распо- |
|
|
||||||||||||
ложим в центре параллелепипеда. Граничные условия с учетом |
||||||||||||||
экстраполяции размеров для (6.20) имеют вид: |
|
|
||||||||||||
|
|
|
|
|
|
a |
, y, z) = 0, |
|
|
|||||
|
|
|
|
Ф(± |
2 |
|
|
|||||||
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, z) |
= 0, |
|
(6.21) |
||||
|
|
|
|
Ф(x,± |
2 |
|
||||||||
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
) |
= 0, |
|
|
||
|
|
|
|
Ф(x, y,± |
2 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где a, b, c – экстраполированные границы параллелепипеда (актив- |
||||||||||||||
ной зоны) критического ЯР, a = a0 +2 δ, так как поправка добав- |
||||||||||||||
ляется дважды – для параллельных граней параллелепипеда с ре- |
||||||||||||||
альными границами активной зоны а0, b0, c0. |
|
переменных |
||||||||||||
Решение |
(6.20) |
будем |
искать |
|
с |
разделением |
||||||||
Ф(x, y, z) = X (x)Y ( y)Z(z). Подставляя это выражение в уравнение |
||||||||||||||
(6.20), получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
145