Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

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

.pdf
Скачиваний:
43
Добавлен:
12.11.2022
Размер:
1.91 Mб
Скачать
(4.58/)

min

E

,

E

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

+

 

 

α Σs (Е)Ф(E)dE

 

1

+Sδ(E E0 ) = 0 .

(4.55)

 

 

−α)E

 

 

E

(1

 

 

В заключение дадим несколько новых определений.

 

Плотностью рассеяния нейтронов будем называть скорость

рассеяния нейтронов в единичном фазовом объеме dVdE

 

Fs (Е) = Σs (E)Ф(E) ,

(4.56)

т.е. в единичном объеме пространства dV и в единичном интервале энергий dE около энергии Е.

Аналогично плотностью столкновений нейтронов будем на-

зывать величину

 

F(Е) = Σ(E)Ф(E) =[Σs (E) + Σa (E)]Ф(E) ,

(4.57)

являющуюся скоростью столкновений нейтронов в единичном фазовом объеме.

Плотностью замедления нейтронов будем называть функцию j(Е), определяющую количество нейтронов в единичном объеме пространства, энергия которых в единицу времени изменяется от значений энергии выше, чем Е, до значений энергии ниже, чем Е. Используя рассуждения, примененные выше для получения выражения скорости процесса (4) в уравнении (4.55), с учетом dE = (E −αE) и формулы (4.35) можно получить

 

 

min

E

,

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

j(E) =

 

 

 

 

α

Σs (Е)Ф(E)dE

 

dE

 

=

 

 

 

 

 

 

 

−α)E

 

 

 

 

E

 

 

 

 

 

 

 

(1

 

min

E

,

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

α

Σs

 

(Е)Ф(E)

 

(E −αE)dE

=

 

 

 

 

 

(1−α)E

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

min E

 

,

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

α Fs (Е)

(E −αE)dE

.

 

(4.58)

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

(1−α)E

 

 

 

 

В случае стационарного монохроматического источника нейтронов (с энергией нейтронов Е0 и постоянной мощностью S), расположенного в непоглощающей нейтроны среде, плотность замедления нейтронов j(Е) не зависит от энергии нейтронов и

j(Е) = S.

111

4.4. Энергетическое распределение замедляющихся нейтронов в бесконечной гомогенной непоглощающей среде

Случай 1. Легкие ядра. При замедлении на водороде (А = 1) без поглощения в среде ( Σa = 0 ) α = EE= AA +11 2 = 0 , т.е. нейтрон

отдает практически всю энергию после однократного акта рассеяния. В силу этого и определения (4.56) уравнение замедления (4.55) примет вид

Е

 

dE

 

 

Fs (Е) + 0

Fs (Е)

+Sδ(E E0 ) = 0 .

(4.59)

E

E

 

 

 

Решение данного уравнения с особенностью при Е = Е0 имеет вид

Fs (E) = ES .

Согласно (4.56) плотность потока замедляющихся на легких ядрах нейтронов (спектр Ферми) имеет вид

Ф(E) =

S

, E < E0.

(4.60)

Σs E

 

 

 

Переходя в выражении (4.60) от переменной энергии к переменной летаргии согласно (4.53) Ф(u) = EФ(E) , имеем

F(E) = S = const и Ф(u) =

S

= const .

(4.61)

 

 

Σs

 

Случай 2. Тяжелые ядра. При замедлении на тяжелых ядрах (А > 1) без поглощения в среде ( Σa = 0 ) нейтрон отдает энергию

порциями и уравнение такого ступенчатого замедления (4.55) в пределах первой ступеньки замедления имеет вид

 

1

 

Е0

 

dE

 

 

Fs1 (Е) +

 

E

Fs1 (Е)

+Sδ(E E0 ) = 0 ,

 

1

α

E

 

 

 

 

αE0 E E0 .

(4.62)

Для второй и последующих (i = 2,…) ступенек замедления слагаемое источника отсутствует

112

 

1

Е0

 

dE

 

 

Fsi (Е) +

E

Fsi (Е)

= 0 .

(4.63)

1−α

E

Решение уравнения (4.62) в пределах первой ступеньки замедления имеет вид:

 

 

 

 

S

 

E

α

 

 

 

 

 

 

 

 

 

 

 

 

F (E) =

 

 

 

1−α

 

+Sδ(E E ) ,

αE E E .

(4.64)

 

 

 

 

0

 

 

 

 

 

 

 

1

 

 

s

 

1

−α

 

 

 

 

0

 

0

 

 

0

 

 

 

 

 

E1−α

 

 

 

 

 

 

 

E0

 

 

E0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С учетом определения летаргии (4.45)

u = ln

,

т.е.

= eu .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

E

Согласно (4.53)

Fs (u) = EFs (E) и первое слагаемое функции (4.64)

с переменной летаргии имеет вид

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

(u) =

e

 

u .

 

 

 

 

(4.65)

 

 

 

 

 

 

 

 

1−α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

1−α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В пределах второй ступеньки замедления из уравнения (4.62)

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (Е) =

 

 

1

αE0

F (Е/ )

dE

+

 

1 E

/α F (Е)

dE

= 0 . (4.66/)

 

 

 

E

 

 

 

 

 

1

−α

E

 

 

 

E

s2

s2

 

1−α αE

s1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

Функция Fs2(E) терпит разрыв в точке Е = αЕ0 по причине того, что нейтроны при однократном рассеянии от энергии Е0 не могут перейти сразу во вторую ступеньку замедления, но формируют фундамент для нее при первом однократном рассеянии (на первой ступеньке замедления). В этом случае

F (αЕ ) =

 

 

1

E0

F (Е)

dE

.

(4.67)

 

 

 

1

−α αE

 

s2 0

s1

E

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

Решение для третьей и последующих решений замедления сложны.

Функция плотности замедления и ее производные в силу ступенчатого процесса замедления терпят разрывы на границах ступе-

нек замедления при Е = αЕ0, α2Е0, α3Е0 и т.д.

В диапазоне энергий, достаточно удаленном от энергии источника Е0, т.е. для нейтронов, испытавших значительное число соударений с ядрами (на k-й ступени), уравнение замедления имеет вид (4.66) с заменой 2 k. Решением такого уравнения является функция

113

F

(E) =

S

, Ф(E) =

S

.

 

 

 

(4.68)

 

 

 

 

 

 

 

sk

 

ξE

 

EΣs (E)ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

В соответствии с определением (4.56) Ф(E) =

S

. Пере-

EΣs (E)

 

 

ξ

ходя от переменной энергии к переменной летаргии согласно (4.53)

Fsk (u) = EFsk (E) , Ф(u) = EФ(E) , имеем

F

(u) =

 

S

 

= const и

Ф(u) =

S

 

.

(4.68/)

 

 

 

 

 

sk

 

 

ξ

 

Σs (u)ξ

 

 

 

 

 

 

Следует отметить для ф.(4.73/), что спектр Ферми (4.62) устанавливается при замедлении нейтронов на тяжелых ядрах (А > 1) только при отсутствии поглощения нейтронов в среде и только в области энергий, удаленной на достаточно большое число ступенек замедления от энергии внешних источников Е0 (так называемая асимптотическая область энергий).

Случай 3. Смесь ядер. При замедлении на смеси ядер без поглощения в среде ( Σa = 0 ), проводя аналогичное рассмотрение для

плотности рассеяния в области энергий удаленной от энергии внешних источников на достаточно большое число ступенек замедления, можно получить

 

 

 

 

 

 

F

(E) =

 

S

,

(4.69)

 

 

 

 

 

 

 

 

 

 

 

sk

 

ξсмE

 

 

 

 

 

Σs _ i

 

 

 

 

 

 

где

 

см =

 

 

i – среднелогарифмическая потеря энергии в сме-

ξ

ξ

 

 

 

i

Σs

 

 

 

 

 

си, ξi – среднелогарифмическая потеря энергии при соударении на ядрах i-го типа, Σs (E) = Σs _ i (E) – сечение рассеяния смеси при

i

условии одинаковой (похожей, изменяющейся по одному закону) зависимости сечений компонентов смеси от энергии.

4.5.Энергетическое распределение замедляющихся нейтронов

вбесконечной гомогенной поглощающей среде

Ранее указывалось, что возможно нерезонансное (часто слабое) поглощение нейтрона средой и резонансное (сильное) поглощение нейтрона в среде на узком изолированном резонансе или на группе

114

узких изолированных резонансов, например, резонансное поглощение в уране-238.

Рассмотрим общее уравнение замедления (4.55) для указанных случаев.

Случай 1. Замедление на легких ядрах (водороде, А = 1) с по-

глощением. Водород относительно слабо поглощает нейтроны. Этот случай будет справедлив и при сильном разбавлении тяжелых ядер водородом NН >> Nя, так как при этом для макросечений рассеяния водорода и тяжелых ядер будет выполняться

ΣНs (Е) = NHσHs (E) >> Σяs (Е) = Nяσяs (E) .

Уравнение замедления (4.55) для данного случая примет вид

 

 

E

 

dE

 

 

−Σs (Е)Ф(E) −Σa (Е)Ф(E) +0 Σs (Е)Ф(E)

+Sδ(E E0 ) = 0

 

или

 

E

 

E

 

 

 

 

 

 

 

Σ(Е)Ф(E) = [Σs (Е) a (Е)]Ф(E) =

 

E

 

(Е)Ф(E) dE

 

 

 

 

= 0

Σs

+Sδ(E E0 ) .

(4.70)

E

 

E

 

 

 

 

После ввода множителя

 

 

 

 

y(E) = Σs (E) / Σ(E) ,

0 < y(E) 1 ,

(4.71)

уравнение (4.70) примет вид

 

 

 

 

 

E

y(Е)Σ(Е)Ф(E) dE+Sδ(E E0 ) .

 

Σ(Е)Ф(E) =

0

(4.72)

 

E

 

E

 

Решение данного уравнения приводит в итоге к выражению для

плотности замедления (см. (4.58))

 

 

 

 

 

 

 

 

Е0

Σ

 

(E) dE

 

j(E) = F(E)E = y(E0 )S exp

 

а

 

 

.

(4.73)

 

 

 

 

E

Σ(E)

E

 

 

 

 

 

 

 

 

Всоответствии с определением плотности столкновений (4.57)

иформулы (4.73) можно записать

Ф(E) =

j(E)

, 0 < E < E0.

(4.74)

Σ(E)E

 

 

 

Отсюда следует, что при слабом поглощении нейтронов в среде Σа (E) << Σs (E) , Σ(E) = Σs (E) + Σа (E) ≈ Σs (E) и с учетом определения плотности замедления (4.58/) j(Е) = S.

115

Плотность потока нейтронов в первом приближении описывается спектром Ферми

Ф(E) =

j(E)

=

S

, 0 < E < E0.

(4.75)

 

Σs (E)E

 

Σs (E)E

 

 

Если ввести в рассмотрение вероятность избежать поглощения при замедлении нейтронов на легких ядрах от энергии Е0 до энер-

гии Е в виде

j(E)

 

 

 

 

 

 

 

 

ϕ(E , E) =

,

 

j(E ) = y(E )S ,

 

 

 

 

0

j(E0 )

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

то из формулы (4.73) следует, что

 

 

 

 

 

 

 

 

 

 

Е0

Σ

 

(E) dE

 

ϕ(E0 , E) = exp

 

а

 

 

.

(4.76)

 

 

 

 

 

 

E

Σ(E)

E

 

 

 

 

 

 

 

 

 

 

Из выражения (4.76) следует, что при отсутствии поглощения в среде Σа (E) = 0 и ϕ(E0 , E) =1.

Случай 2. Замедление на тяжелых ядрах (А > 1) с поглоще-

нием. Уравнение замедления в данном случае имеет вид (4.55)

−Σs (Е)Ф(E) −Σa (Е)Ф(E) +

min

E

,

E

 

 

 

 

 

 

 

 

0

 

 

 

 

 

+

 

 

α Σs (Е)Ф(E)dE

 

1

+ Sδ(E E0 ) = 0 ,

 

 

−α)E

 

 

E

(1

 

технология решения данного уравнения «по ступенькам замедления» такая же, как и в случае замедления на тяжелых ядрах в отсутствие поглощения. Так как функция Σa (Е) может быть весьма

сложна, то решение данного уравнения аналитически для общего случая не представляется возможным. В связи с этим рассмотрим частные случаи для области энергии, достаточно удаленной на большое число ступенек замедления от энергии внешних источников (в асимптотической области энергий):

-случай 2а – случай слабого поглощения нейтронов при замедлении на тяжелых ядрах;

-случай 2б – случай поглощения нейтрона на узком изолированном резонансе (на группе узких изолированных резонансов) тяжелых ядер.

116

Указанная область и перечисленные случаи практически интересны тем, что именно в данной области энергий уже сильно замедлившиеся, эпитепловые нейтроны часто попадают в область резонансов тяжелых ядер, например, урана-238.

Случай 2а. Слабое поглощение нейтронов при замедлении на тяжелых ядрах.

Так как рассматривается слабое поглощение Σа (E) << Σs (E) , Σ(E) = Σs (E) + Σа (E) ≈ Σs (E) , то в первом приближении подобно

рассмотрению на водороде с поглощением (4.56) и рассмотрению замедления на тяжелых ядрах без поглощения (4.43) можно предположить, что плотность потока нейтронов будет совпадать со спектром Ферми

Ф(E) =

j(E)

.

(4.77)

 

 

EΣs (E)ξ

 

 

 

Действительно, в предыдущих лекциях было показано, при замедлении на тяжелых ядрах нейтрон отдает свою энергию растянуто, малыми порциями, уменьшающимися с ростом массового числа А, в то время как на водороде даже при единичном столкновении нейтрон может потерять практически всю кинетическую энергию. Отличие при замедлении нейтронов на тяжелых ядрах от их замедления на водороде для слабого поглощения заключается в том, что в знаменателе появляется множитель – среднелогарифмическая

потеря энергии ξ . Сопоставьте формулы (4.58/), (4.60), (4.61),

(4.68), (4.75) и (4.77).

Изменение плотности замедления в некотором интервале энергии будет равно скорости поглощения нейтронов в данном интервале энергии

 

 

 

 

 

 

j = j(E +

E) j(E) = Σа (E)Ф(Е) Е .

(4.78)

 

Подставляя

в (4.78)

выражение

(4.77),

получим при

j→0

dj

=

Σa (E)

 

dE

с решением в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

Σs (E)ξ E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Е0

Σ

 

(E)

dE

 

 

 

 

 

 

 

j(E) = j(E0 ) exp

 

а

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

Σ(E)

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По аналогии с формулой (4.76) выражение для вероятности избежать поглощения при замедлении нейтронов на тяжелых ядрах

117

от энергии Е0 до энергии

Е

 

при

слабом

поглощении

( Σа (E) << Σs (E) , Σ(E) = Σs (E) а (E) ≈ Σs (E) )

 

 

 

Е0

 

Σ

а

(E) dE

 

ϕ(E0

, E) = exp

 

 

 

 

 

 

 

.

(4.79)

 

 

 

 

 

 

 

 

 

Σ(E)ξ

 

 

E

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

Чем больше массовое число А, тем меньше среднелогарифмическая потеря энергии при одном соударении ξ и тем больше веро-

ятность поглощения нейтрона тяжелыми ядрами при замедлении.

Случай 2б. Поглощение нейтрона на узком изолированном резонансе (на группе узких изолированных резонансов) тяжелых ядер (приближение узких резонансов – УР-приближение).

Введем в рассмотрение ряд характеристик узкого изолированного резонанса (рис. 27): энергию максимума резонанса Еi, высоту резонанса σа0i , ширину Гi, измеряемую на полувысоте резонанса, рас-

стояние до соседнего резонанса Di.

В областях резонанса поглощение сильное

Σа (E) / Σs (E) >>1,

а в межрезонансных областях поглощение относительно слабое

Σа (E) / Σs (E) <<1 ( Σа (E) = σa (E)Nядер ).

Резонанс называется узким, если его ширина Г много меньше ступеньки замедления

Рис. 27. Характеристики резонансов

 

E2

 

E

 

A 1

2

 

α =

 

=

 

=

 

 

 

 

E1

E

A +1

 

 

 

 

 

(Е1 и Е2 – энергии до и после однократного соударения нейтрона с ядром, см. формулы (4.20), (4.21)):

Гi << (1−α)E . (4.80)

Резонанс называется изолированным, если расстояние до соседнего резонанса D много больше ступеньки замедления:

Di >> (1−α)E .

(4.81)

118

В силу слабого поглощения в областях, примыкающих к узкому изолированному резонансу (УИР) падающий на него поток замедляющихся нейтронов можно в первом приближении считать распределенным по спектру Ферми подобно вышерассмотренным случаям (4.68), (4.75), (4.77). Изменение плотности замедления в малой области энергий узкого резонанса Е = Г по аналогии с (4.78) будет равно скорости поглощения нейтронов в данном интервале энергии

j = Σа (E)Ф(Е)Г.

(4.82)

Основная трудность – определение плотности потока нейтронов Ф(Е), так как происходит отклонение в меньшую сторону от распределения спектра Ферми в связи с интенсивным поглощением («выеданием») нейтронов в области энергий УИР. Для оценки Ф(Е) определим скорость прибыли нейтронов в область энергий УИР Е = Г в результате упругого рассеяния из вышележащих интервалов энергии. Для простоты рассмотрим стационарную задачу, в которой данная скорость равна скорости убыли нейтронов из интервала Е = Г в результате упругого рассеяния нейтронов, но при отсутствии их поглощения в рассматриваемом интервале энергии

 

 

 

 

j = Σs (Eф (Е)Г,

(4.83)

где Фф (E) =

j(E)

 

– спектр Ферми (4.77). Учтем наличие по-

 

 

 

EΣs (E)ξ

 

 

 

 

глощения в интервале энергий УИР путем домножения выражения

(4.83) на вероятность поглощения

нейтронов на

резонансе

Σа (E) / Σ(E) . Таким образом, выражение для скорости поглощения

нейтронов в области энергий резонанса

Σа (E)

 

 

j = Σs (Eф (Е

.

(4.84)

 

 

Σ(E)

 

Запишем тождество для равенства скоростей поглощения (4.84) и истинной, где неизвестна истинная функция плотности потока нейтронов Ф(Е) в области энергий УИР

Σs (Eф (Е

Σа (E)

= Σа (E)Ф(Е)Г ,

(4.85)

Σ(E)

 

 

 

откуда истинная функция плотности потока нейтронов Ф(Е) в области энергий УИР

119

Ф(Е) = Фф (Е)

Σs (E)

=

j(E)

.

(4.86)

 

 

Σ(E)

EΣ(E)ξ

 

 

 

 

Выражение (4.86) называется спектром Вигнера.

Подставляя спектр Вигнера (4.86) в формулу (4.82), получим истинное выражение для скорости поглощения нейтронов в области энергий УИР

j = Σа (E)Ф(Е= j1

Σа (E) Г

 

 

 

 

 

 

.

(4.87)

Σ(E)

 

 

Е

ξ

С учетом (4.87) вероятность избежать резонансного поглощения

на УИР

j2

 

j1 j

 

j

 

Σа (E

)

 

Г

 

 

ϕ =

=

=1

=1

 

.

(4.88)

j

j

j

 

 

 

 

 

 

 

Σ(E)ξ Е

 

 

1

 

1

 

1

 

 

 

 

 

 

 

Как бы ни было велико сечение поглощения в резонансе, вероятность φ избежать резонансного поглощения нейтрона на УИР отлична от нуля, если средняя потеря энергии нейтрона при однократном рассеянии больше ширины резонанса. Однако в лекциях ранее было показано (см. формулу (4.19)), что при рассеянии на легких ядрах нейтрон отдает энергию большими порциями, в предельном случае – на водороде при числе соударений, стремящихся к единице. В силу этого вероятность φ избежать резонансного поглощения нейтрона на УИР легких ядер стремится к единице φ 1. При рассеянии на тяжелых ядрах нейтрон отдает энергию малыми порциями, тем меньше, чем больше массовое число ядра А (приближение бесконечной массы – БМ-приближение) со ступенькой

 

E2

 

E

 

A 1

2

замедления α =

 

=

 

=

 

 

 

(Е1 и Е2 – энергии до и после од-

E1

E

A +1

 

 

 

 

 

нократного соударения нейтрона с ядром (см. (4.20), (4.21)). В результате этого вероятность φ избежать резонансного поглощения нейтрона на УИР тяжелых ядер стремится к нулю φ 0.

Поглощение нейтронов на группе из N узких изолирован-

ных резонансов УИР. Учитывая вышеприведенное рассмотрение для одного УИР, можно считать следующее:

- перед УИР № 1 плотность потока замедляющихся нейтронов описывается распределением Ферми в силу относительно слабого поглощения;

120