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

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

вии с физическим смыслом в конкретной задаче теории замедления нейтронов:

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

(5.29/)
(5.30)

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

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