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

Крючков Теория переноса нейтронов 2007

.pdf
Скачиваний:
2975
Добавлен:
16.08.2013
Размер:
7 Мб
Скачать
☆

E∫0

E

 

 

h(E

′

)

 

E0

 

 

 

1

 

 

 

 

 

 

 

 

′

)

 

 

 

 

 

 

 

E0

 

 

 

E0

d E

′

 

 

′

)

 

d E

′

 

 

 

d E

′

 

 

 

 

 

Σa (E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Σa (E

 

 

 

 

=

∫

 

 

 

−

 

 

 

 

 

 

 

= ln

 

 

 

 

−

 

 

∫

 

 

 

 

 

 

 

 

;

 

 

E′

 

 

 

E

 

 

 

E′

 

 

Σtot (E′)

 

 

 

 

 

 

E

 

 

 

 

E

E′

 

 

Σtot (E′)

 

 

 

 

 

E0

d E′

h(E′)

 

 

 

 

 

 

 

 

E

0

 

 

E0

 

d E′ Σ

a

(E′)

=

 

 

 

 

exp ∫

 

 

 

 

 

 

 

= exp ln

 

 

 

 

− ∫

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E′

 

 

 

 

 

 

 

 

 

 

 

E′

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

E

 

 

 

 

 

 

Σtot (E′)

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

E0

d E′

Σ

a

(E′)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

0

exp −

E∫

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E′

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

Σtot (E′)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(E)=

 

q h(E

0

)

 

 

 

 

E0

d E′ Σ

a

(E′)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp −

∫

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

(3.69)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E′

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

E

 

 

 

Σtot (E′)

 

 

 

 

 

 

 

Найдем плотность замедления j(E) в рассматриваемой среде

 

E

 

FS (E′)E

E

 

 

ΣS (E′)Φ(E′)E

 

 

j(E)= ∫0 d E′

= ∫0 d E′

 

=

E′

 

 

 

 

 

 

E

 

E

 

 

 

 

E′

 

 

E

ΣS (E′)E

 

 

 

E

 

 

 

 

E

 

 

 

= ∫0d E′

Σtot (E′)Φ(E′)= ∫0 d E′ h(E′)

 

( (E′)+

Σtot (E′)E′

E′

E

 

 

 

E

 

 

 

 

 

 

 

 

E0

 

(E′)

 

 

q h(E0 )E

 

+ qδ(E′−E 0 ))= E ∫ dE′ h(E′)

 

 

 

+

 

 

 

 

 

.

 

 

E′

 

E0

 

 

 

 

E

 

 

 

 

 

 

 

При выводе было использовано выражение для h(E) и формула (3.65). Функция (E) удовлетворяет уравнению (3.66), следова-

тельно, для плотности замедления окончательно получим:

E0

′

 

 

′

 

 

qh(E0 )

 

 

 

 

(E )

 

 

 

 

 

j(E)= E ∫ dE h(E )

 

 

 

+

 

 

 

 

= E (E),

 

E

′

 

 

E0

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или, используя (3.69):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E0

 

 

Σ

a

(E′)

 

 

j(E)= E (E)= qh(E0 ) exp −

∫

dE′

 

 

.

(3.70)

 

 

 

 

 

 

 

E

 

 

Σtot (E′)E′

 

141

В свою очередь, на основе вышеприведенных выкладок можно записать:

 

 

 

 

(E)=

j(E)

.

(3.71)

 

 

 

 

 

 

 

 

 

 

 

 

E

 

Если ∑a = 0 , то для бесконечной среды

j(E)= q и из (3.71) сле-

дует, что

(E)=

j(E)

=

q

, т.е. в среде

устанавливается спектр

E

 

 

 

 

 

E

 

Ферми.

Используя выражение для плотности замедления (3.70), найдем вероятность избежать поглощения при замедлении от энергии Е0 до энергии Е:

 

 

j(E)

 

E

0

 

Σ

 

(E′)

 

 

 

 

ϕ(E0

, E)=

 

 

dE′

a

 

1

 

 

= exp −

∫

 

 

 

 

.

(3.72)

J (E0 )

 

 

 

 

 

 

 

E

 

Σtot (E′) E′

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В прил. 4 приведен пример использования формулы (3.72) для расчета избежать поглощения на первом резонансе урана-238 при замедлении на водороде.

3.12.Замедление на ядрах с атомной массой больше единицы

сучетом поглощения

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

вданной постановке задачи без дополнительных приближений (как

вслучае замедления на водороде) решить не представляется возможным.

По отношению к энергии источника (это область мегаэлектронвольт) область резонансов (до сотен килоэлектронвольт) – асимптотическая область энергий, т.е. область, отстоящая от энергии источника более, чем на три ступеньки замедления на замедлителе.

142

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

Случай слабого поглощения. Случай слабого поглощения оз-

начает, что ∑∑as((EE))<<1, и реализуется в области неразрешенных

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

Рассмотрим энергетический интервал ∆E (рис. 3.19) в области [E1, E2], где реализуется случай слабого поглощения, и баланс нейтронов в этом энергетическом интервале.

 

 

 

 

 

j(E)

∆E

 

j(E+∆E)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E1

 

 

 

 

 

 

E2

Поглощение

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3.19. Схема баланса нейтронов в энергетическом интервале

Исходя из определения плотности замедления и основываясь на законе сохранения нейтронов, можно утверждать, что изменение плотности замедления на энергетическом интервале ∆Е равно скорости реакции поглощения нейтронов в интервале ∆Е, т.е.

143

j(E + ∆E)− j(E) = ∑a (E)Φ(E)∆E .

(3.73)

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

Φ(E)≈

j(E)

.

(3.74)

ξ ∑S (E) E

 

 

 

В случае отсутствия поглощения j(E) = q и (3.74) представляет

собой спектр Ферми в асимптотической области энергий в непоглощающей среде.

Подставим (3.74) в (3.73):

j(E + ∆E)− j(E)≡ ∆j(E)= ∑a(E)j(E)∆E

ξ∑S (E) E

и перейдем к дифференциальному уравнению

dj((E)) = ∑a ((E)) dE , j E ξ∑S E E

решением которого на рассматриваемом интервале энергий [E, E2] является функция

 

E

2

Σa (E′)

 

 

 

 

 

 

dE′

j(E)= C exp −

∫

 

 

 

.

ξΣS (E′)

 

 

E

 

E′

 

 

 

 

 

144

Учитывая слабое поглощение, т.е. равенство ∑S (E′)≈ ∑ tot (E′),

получим выражение для вероятности избежать поглощения при замедлении на интервале [E1, E2]:

 

 

 

j(E )

 

 

 

E2

ϕ(E

 

, E )=

 

 

 

2

1

= exp

 

−

∫

 

 

1

j(E2 )

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

1

Σ (E′)

 

 

 

dE′

 

a

 

 

.

(3.75)

ξΣtot (E′)

E′

 

 

 

 

 

Случай узкого изолированного резонанса. Определим веро-

ятность избежать поглощения в энергетическом интервале [E1, E2] при замедлении на ядрах с атомной массой больше единицы, если в рассматриваемом энергетическом интервале в сечении поглощения (радиационного захвата) поглотителя расположен узкий изолиро-

ванный резонанс (рис. 3.20).

На рис. 3.20 выделены две характерные области энергий:

область I, где реализуется случай слабого поглощения, т.е.

∑a << ∑S ;

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

σа(E) σ°

Г

12 σ°(E)

∆E′

II I

E1

Er

E2 D

E

Рис. 3.20. Зависимость сечения поглощения от энергии

145

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

с ядром замедлителя Г<< E2r (1−αз ).

Резонанс называется изолированным, если расстояние между ним и следующим, расположенным выше, резонансом много больше ступеньки замедления для замедлителя D >> Er (1− αз ).

Узкие изолированные резонансы реально существуют. Например, рассмотрим первый резонанс 238U при замедлении на углероде.

Первый резонанс 238U имеет следующие параметры: Er = 6,7 эВ, Г = = 26,5·10-3 эВ. Второй резонанс 238U лежит при энергии 21 эВ. Для

углерода α = 0,716. Проверим выполнимость условий узости и изолированности первого резонанса урана-238 при замедлении на углероде:

Er (1− αC ) = 6,7 (1− 0,716) 0,95 эВ; 2 2

Er (1− αC ) 1,9 эВ.

Оба из условий 26,5·10–3 << 0,95 и (21–6,7) = 14,3 >> 1,9 эВ с хорошей точностью выполняются, поэтому первый резонанс урана238 можно считать узким и изолированным при замедлении на углероде.

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

Условие узости резонанса позволяет утверждать, что вероятность того, что энергия нейтрона при замедлении попадет в область резонанса, т.е. в область энергий [E – 3Г, E + 3Г], мала, поскольку эта область существенно меньше ступеньки замедления на замедлителе.

146

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

ln 21

резонанса нейтрон в среднем испытает N = 6,7 = 1,142 ≈ 7,2

ξC 0,158

столкновений с ядрами среды в области слабого поглощения. Выберем в пределах резонанса (в интервале энергий [E1, E2]) ин-

тервал ∆E′ и рассмотрим уравнение баланса нейтронов (3.75):

j(E′+ ∆E′)− j(E′)= ∆j(E′)= ∑a (E′)Φ(E′)∆E′.

Но если в случае слабого поглощения в качестве спектра нейтронов Φ(E′) можно было использовать невозмущенный спектр

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

Найдем спектр нейтронов в области узкого изолированного резонанса Φ(E′) (рис. 3.21).

Область резонанса

Область слабого поглощения (спектр Ферми)

E1

E2

E′

E

Е′

∆E′

α

 

 

 

Рис. 3.21. Схема для нахождения спектра нейтронов в области резонанса

147

Резонанс

 

–

узкий,

поэтому

выполняется равенство

(E2 − E′)<<

E′

− E2

 

 

∆E′ подавляющее боль-

 

 

 

, т.е.

в интервал

α

 

 

 

 

 

 

 

 

 

E′

 

шинство нейтронов приходит из интервала

E2

,

 

 

. А это – об-

 

 

 

 

α

 

ласть слабого поглощения, в которой установился невозмущенный спектр Ферми Φ0 (E) вне зависимости от того, есть резонанс или

его нет, поскольку наличие резонанса ни как не влияет на спектр нейтронов в более высокой области энергий. Наличие или отсутствие резонанса влияет только на спектр нейтронов в области [E ′, E2], размеры которой по сравнению со ступенькой замедления, а следовательно, и вклад в баланс нейтронов в интервале ∆E′ малы. Таким образом, в интервал ∆E′ приходит примерно одинаковое число нейтронов за счет рассеяния при более высоких энергиях вне зависимости от наличия или отсутствия резонанса. Но в случае стационарной задачи скорость генерации нейтронов в ∆E′ всегда равна скорости исчезновения нейтронов из ∆E′ . Таким образом, получается логическая цепочка, изображенная на рис. 3.22.

В изображенной на рис. 3.22 логической диаграмме все знаки равенства, за исключением жирного, обсуждены выше. А жирный знак равенства следует из логики приведенной диаграммы.

Таким образом, скорость исчезновения нейтронов из ∆E′ не зависит от наличия или отсутствия резонанса. Но в случае наличия

резонанса

эта

величина

определяется

соотношением

∑tot (E′)Φ(E′)∆E′,

поскольку любое столкновение с ядром среды

приводит или к рассеянию, или к поглощению нейтрона. А в случае отсутствия резонанса реализуется случай слабого поглощения, т.е. Σtot ≈ ΣS , и устанавливается невозмущенный спектр нейронов –

спектр Ферми – Φ0 (E′) . Поэтому в случае отсутствия резонанса

скорость

исчезновения нейтронов из ∆E′ рассчитывается как

∑S (E′)Φ0

(E′)∆E′ .

148

Скорость

генерации

нейтронов

в

 

Скорость исчезновения нейтронов

из

 

∆E′ в случае наличия резонанса

 

 

∆E′ в случае наличия резонанса

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Скорость

генерации

нейтронов

в

 

Скорость

исчезновения нейтронов

из

 

∆E′ в случае отсутствия резонанса

 

 

∆E′ в случае отсутствия резонанса

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3.22. Логическая диаграмма

Итак, получено следующее равенство:

∑tot (E′)Φ(E′)∆E′ = ∑S (E′)Φ0 (E′)∆E′ ,

(3.76)

где

′

)

 

 

 

 

 

Φ0 (E′)=

j(E

.

 

ξ∑S (E′)E′

 

 

 

 

Отметим, что поскольку сечение рассеяния замедлителя фактически не зависит от энергии, а плотность замедления слабо меняется на малом интервале [E1, E2] в случае слабого поглощения, то правая часть равенства (3.76) фактически не зависит от энергии. А это означает, что плотность столкновений (левая часть (3.76)) постоянна в области узкого изолированного резонанса. Поскольку плотность замедления слабо меняется на малом интервале [E1, E2] в случае слабого поглощения, то в формуле для невозмущенного спектра Ферми (Ф0) заменим в числителе j(E ′) на j(E2), где j(E2) – плотность замедления, «падающая» на резонанс:

Φ0

(E′)=

j(E′)

≈

j(E2 )

.

ξ∑S (E′)E′

ξ∑S (E′)E′

 

 

 

 

Подставим последнее выражение в (3.76) и выразим спектр нейтронов Ф(Е) в пределах узкого изолированного резонанса:

149

′

)=

∑S (E′)Φ0 (E′)

 

j(E2 )

 

Φ(E

∑tot (E′)

=

ξ∑tot (E′)E′ .

Таким образом, в пределах узкого изолированного резонанса устанавливается следующий спектр:

Φ(E) =

j(E2 )

,

(3.77)

ξ∑tot (E)E

который называется спектром Вигнера. Отметим, что в отличие от спектра Ферми в знаменателе спектра Вигнера стоит не макроскопическое сечение рассеяния замедлителя, которое фактически не зависит от энергии, а полное сечение, которое сильно зависит от энергии, поскольку включает сечение поглощения поглотителя:

Σtot (E) = ΣaП(Е) + ΣSП(Е) + ΣSЗ(Е) + ΣaЗ(Е) .

Если пренебречь сечением поглощения замедлителя по сравнению с сечением рассеяния замедлителя и сечением рассеяния поглотителя по сравнению с сечением поглощения поглотителя, то для полного сечения получится следующее выражение:

Σtot (E) ΣaП(Е) + ΣSЗ = NПσaП(Е) + NЗσSЗ .

(3.78)

Если в среде очень мало поглотителя (случай бесконечного разбавления), то каким бы большим не было микроскопическое сечение поглощения поглотителя первое слагаемое в (3.78) стремится к нулю, следовательно, полное сечение стремится к сечению рассеяния замедлителя, а спектр Вигнера – к спектру Ферми (кривая 1 на рис. 3.23). При повышении концентрации ядер поглотителя растет влияние первого слагаемого в (3.78) и в спектре нейтронов (3.77) начинает наблюдаться «провал», который ведет себя противоположным образом по отношению к микроскопическому сечению поглощения (кривая 2 на рис. 3.23).

150