Материал: Баранник Лекции по курсу Теория переноса нейтронов 2012

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам
полное число ТН поколения

 

число нейтронов эпитепловых энергий,

 

ϕ =

избежавших резонансного захвата в 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

Источник: https://studfile.net/preview/16708734/