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

Савандер В.И., Увакин М.А. Физическая теория ядерных реакторов. Часть 1. Однородная размножающая среда и теория гетерогенных структур

.pdf
Скачиваний:
179
Добавлен:
21.07.2024
Размер:
1.45 Mб
Скачать
V0 Σt0

%

бл

1

 

бл r

r

 

 

(E) =

 

 

(r , E)dr .

Для вычисления потока резо-

Φ

 

V

Ф

 

 

 

бл V

 

 

 

 

 

 

 

бл

 

 

 

нансных нейтронов в блоке воспользуемся уравнением баланса нейтронов с применением вероятностей первых столкновений (всюду в дальнейшем индекс «0» относится к блоку, а «1» – к замедлителю):

0 % 0

 

 

 

E α0

0

%

0

 

 

 

 

 

 

V0 Σt Φ

 

 

 

 

 

 

 

+

(E) = P00 V0 Σs (E ) Φ

 

(E )dE

 

 

 

 

 

E

 

 

 

 

 

 

(8.2.7)

E

 

 

 

 

 

 

 

 

 

 

α

 

%1

 

 

 

 

 

 

 

 

+P10 V1

1 1

′ ′

 

 

 

 

 

 

 

Φ

 

 

 

 

 

 

Σs (E )

(E )dE .

 

 

 

 

 

 

E

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

Φ% яч(E) = 1 . Интегрируя по всей ступеньке замедления,

ξΣp Vяч E

получим

V0 Σt0 Φ% 0 (E) = P00 V0 Σ0p Φ% яч(E) + +P10 V1 Σ1s Φ% яч(E),

Φ% 0(E) = (P00 V0 Σ0p + P10 V1 Σ1s ) Φ% яч(E).

Подставляя это выражение в формулу для вероятности нейтрону поглотится в резонансе, получим:

181

P =

V0 ρП

J

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

r

ξΣ

p

V

 

ef

 

 

 

 

 

 

 

 

 

яч

 

 

 

 

 

 

 

(8.2.8)

 

 

 

 

σar (E) (P00 V0 Σ0p + P10

V1

Σ1s )

Jef =

 

dE.

 

 

 

 

V

Σ0

E

 

 

 

E

 

 

 

 

 

 

 

 

 

 

0

t

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

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

P10 V1 Σ1s = P01 V0 Σt0 , откуда следует:

P

= P

 

V0

Σ0

,

P

=1P .

(8.2.9)

V

Σ

10

01

 

 

01

00

 

 

 

 

1

1

 

 

 

 

Модель тонкого блока. В модели «тонкого блока» принимается, что рассеянием нейтронов в блоке можно пренебречь. Это значит, что основное количество резонансных нейтронов падает на

блок из замедлителя. В этом случае можно положить Σ0p = 0 , тогда

% 0

 

V0 Σt0

 

V1 Σ1s

% яч

%

яч

 

Φ

(E) = (1P00 )

V Σ1

 

V Σ0

Φ

(E) = (1P00 ) Φ

 

(E) . (8.2.10)

 

 

1

s

 

0

t

 

 

 

 

Используя рациональное приближение для вероятностей первых

столкновений 1P00

=

 

1

 

 

 

 

, получим

 

 

+

 

 

0

Σt0

 

 

 

1

L

 

 

 

 

 

%

0

(E) =

 

 

1

 

%

яч

(E) .

 

Φ

 

 

 

 

0

0

Φ

 

(8.2.11)

 

 

 

 

 

 

 

1

+ L

Σt

 

 

 

В результате этих преобразований выражение для эффективного резонансного интеграла будет иметь следующий вид:

Jef =

 

 

σar (E)

dE .

(8.2.12)

 

 

 

 

 

 

 

(1

+ L0

Σt0 ) E

E

r

 

 

 

 

 

 

 

 

 

 

182

Для вычисления резонансного интеграла применим стандартную замену переменных x = ((E E0 )Г2). Представляя формулу Брейта – Вигнера для формы резонансной линии в новой перемен-

ной σar = 1+σ0x2 , получаем

Jef

=

Гn σr0

+∞

dx

 

 

 

 

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1+ x2

L0

ρП σr0

 

 

Г E0 −∞

1+

 

 

 

 

 

 

 

 

 

 

 

1+ x2

 

Вводя обозначение βгет =

 

0 ρП σr0

и интегрируя в пределах от -

L

до +, получим аналитическое выражение для эффективного резонансного интеграла для приближения «тонкого блока:

 

Jef

=

 

Гγ σr0

+∞

dx

=

JR

,

(8.2.13)

 

 

 

 

 

 

2 E0

1+ x2 гет

1гет

 

 

 

 

−∞

 

 

 

где

JR =

 

π Гγ σr0

– обычный резонансный интеграл.

 

 

2 E0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как в гомогенной среде замедлителя и резонансного погло-

тителя

βгом =

σr0 ρП , то в гетерогенной среде роль сечения по-

 

 

Σp

тенциального рассеяния играет величина Σгетp = L10 .

Слабый резонанс:

βгет <<1

 

Jef

 

= J R .

 

 

 

 

 

Сильный резонанс:

βгет >>1

 

Jef

 

=

 

J R

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

βгет

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

4V0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку β

гет

= L

ρ

П

σ

r0

=

 

 

ρ

П

σ

r0

, а

ρ

П

V = M

0

, то

 

 

 

 

 

 

 

 

 

 

S0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

имеем выражение:

183

Jefгет

S0

.

(8.2.14)

 

 

M0

 

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

Jef = A + B

S

.

(8.2.15)

 

 

M

 

Модель толстого блока. Случай «толстого» блока характеризуется тем, что рассеянием резонансных нейтронов в блоке пренеб-

речь нельзя. Это следует из того, что L0 Σ0p >>1. Уравнение баланса нейтронов в ячейке для этого случая имеет вид:

 

 

0

%

0

 

 

 

 

0

% яч

(E)

 

 

1

% яч

(E) . (8.2.16)

 

V0 Σt

Φ

 

(E) = P00 V0 Σs

Φ

 

+ P10 V1 Σs Φ

Используя теорему взаимности:

 

P

 

V Σ0

= P

V Σ1 , запишем

 

 

 

 

 

 

 

 

 

 

 

 

 

01

 

0 t

10

1

s

 

 

V Σ0

 

 

 

 

 

 

 

 

 

 

 

 

P

= P

 

 

0

 

t

.

Следовательно,

поток резонансных

нейтронов в

 

 

 

 

 

10

01 V Σ1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

s

 

 

 

 

 

 

 

 

 

 

 

 

блоке имеет вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

Σ0

 

 

 

 

яч

 

 

 

 

 

 

 

 

%

 

 

s

P + P

 

%

(E) .

 

 

(8.2.17)

 

 

 

 

Φ (E) =

0

Φ

 

 

 

 

 

 

 

 

 

 

 

 

 

00

01

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Σt

 

 

 

 

 

 

 

 

 

 

Поскольку P00 =1P01 , то

184

 

0

 

 

Σ0

 

 

 

 

 

 

 

 

яч

 

%

 

 

s

(1

 

 

 

 

 

 

%

 

Φ (E) =

Σ0

P01) + P01 Φ (E) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

(8.2.18)

 

Σ0

 

 

 

Σ0

 

 

 

 

 

 

 

 

 

 

 

яч

 

 

 

 

=

s

+ P

1

s

 

 

%

 

 

 

 

Σ0

Σ0

Φ (E).

 

 

 

 

01

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

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

 

 

0

 

 

 

0

 

 

dE

 

Jef =

σar (E)

Σs

+ P01

1

Σs

 

 

.

 

E

 

Σ0

 

 

 

Σ0

 

 

E

r

 

t

 

 

 

t

 

 

 

 

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

интеграла. Вводя стандартную замену переменных:

 

 

 

 

 

 

 

 

 

 

x =

E E0r

,

 

dx =

2

dE,

dE =

Г

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г

2

 

 

 

Г

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и распространяя интегрирование по х от - до +получим:

 

Jef(1)

 

 

σ0r Гγ

+∞

 

σ

 

 

 

 

 

σ0r Гγ

+∞

dx

 

=

 

 

 

 

 

 

 

s0

 

 

dx =

 

 

 

 

=

 

 

2E0r

 

 

 

 

 

 

2E0r

 

+ x2 гет

 

 

 

 

 

−∞ σs0 (1+ x2 ) + σ0r

 

−∞1

 

=

 

JR

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1гет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где βгет

 

 

σ

 

 

σs

 

 

 

Σ0

= σp + σ0 . Эта составляющая резо-

=

 

0r ;

0 = σp +

s

 

ρП

 

 

 

 

 

σs0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

185

Теперь рассмотрим

второе

 

слагаемое

резонансного интегра-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

dE

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ла

σ

ar

(E) P

1

Σs

 

 

. Воспользуемся рациональным при-

 

 

 

 

 

 

 

 

 

 

01

 

 

 

 

 

Σ0

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ee

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ближением для вероятности P01:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P01 =

 

 

1

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1+

 

 

0 Σt0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

Σar

 

 

 

 

P

 

1Σs

=

 

 

 

 

 

 

 

 

 

 

1Σs

 

 

 

Σt

 

 

− Σs

 

 

=

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

01

 

 

 

0

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

 

0

)

2

 

 

 

 

 

0

 

 

 

 

0

 

2

 

 

 

 

 

 

 

 

Σt

 

1+ L

Σt

 

 

 

 

 

Σt

 

 

 

L

 

(Σt

 

 

 

 

 

L

 

(Σt )

 

 

 

В результате получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

σ0r Гγ

 

 

 

 

 

 

 

 

1

 

 

ρП σ0r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dE

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Jef

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ x2 1+ x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

Г E0r Er

1

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

σ0r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

(ρП)

 

σs0

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

+ x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Σ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

x =

 

 

0

 

;

 

 

 

 

σ

s0

= σ

p

 

+

 

s

.

 

 

 

Путем несложных преобразова-

 

 

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

 

 

 

 

 

Г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ний представим этот интеграл в следующем виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

(σ0r )2

 

 

 

 

+∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Jef

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 E0r L0

 

ρП

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ0r

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

(1+ x

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

σs0

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

(σ0r )2

 

 

 

 

 

 

 

 

 

+∞

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

ρП

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

E0r L

 

 

(σs0 )

 

 

 

−∞

 

+ x2

+

0r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

σ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π(σ0r )2

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ0r

 

3

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4E

 

 

0

 

 

 

 

(σ

 

 

 

)

2

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L ρ

П

s0

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σs0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

186

Главная часть этого интеграла, отвечающая за пространственную неравномерность потока резонансных нейтронов, есть

ρП

 

0 = ρП

4V

M

. Следовательно, этот интеграл

пропорцио-

L

S

 

 

 

 

 

M

S

 

 

 

 

нален величине

 

, а поэтому полный интеграл по всем резонан-

 

S

 

 

 

 

 

 

 

 

 

 

сам можно представить в виде

 

 

 

 

 

 

 

 

 

Jef = A + B

S

.

(8.2.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

В отличие от модели тонкого блока, для толстых блоков эффективный резонансный интеграл пропорционален отношению поверхности блока к его массе.

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

интеграла величину P00 на величину Q00 = P(БЛ БЛ) , которая

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0Σt0

 

 

 

 

 

 

 

 

0Σt0

 

 

L

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

,

 

Q = P(БЛ БЛ) =

 

 

 

 

=

 

 

1C

 

(8.2.20)

 

 

 

 

 

 

 

 

 

 

 

 

 

00

1

+ L0Σt0 С

1

+

 

L0Σt0

 

 

 

 

 

 

 

 

 

 

 

 

1C

 

 

 

где C = Pзам(S S) . Для удобства дальнейших преобразований переобозначим величину C (1C) , так что выражение для вероятности Q00 будет иметь вид:

 

 

 

 

 

0Σt0

 

 

L

 

 

Q00 =

 

 

 

 

C

 

 

.

(8.2.21)

 

 

 

 

 

 

L0

Σt0

1+

 

 

 

 

 

 

C

 

 

187

Таким образом, в тесных решетках как бы возрастает величина эффективной средней хорды:

 

 

 

 

 

 

0

 

4V0

 

4V0

 

 

 

 

 

 

 

 

0

 

L

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

L

=

 

 

 

=

 

=

 

; S

 

= С S

 

.

(8.2.22)

 

 

 

S0 C

Sэф0

 

 

 

эф

 

 

C

 

 

 

эф

 

0

 

 

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

Как в случае тонкого блока, так и в случае толстого блока, эффективный резонансный интеграл зависит от отношения площади поверхности топливного блока S к его массе M . Для тесной решетки эффективная поверхность блока уменьшается и поэтому при всех прочих равных условиях

Jef (тесной решетки) < Jef (широкой решетки) . (8.2.23)

Принцип эквивалентности для резонансного поглощения нейтронов. Как видно из полученных результатов, выражение для эффективного резонансного интеграла как для гомогенной среды, так и для гетерогенной имеют одинаковый вид:

Jefгет =

 

JR

 

,

βгет =

ρП σr0

,

 

 

1гет

 

 

 

 

 

 

 

 

Σp

(8.2.24)

Jefгом =

 

JR

 

 

βгом =

ρП σr0 .

 

 

,

 

 

1гом

 

 

 

 

 

 

 

 

Σp

 

188

Таким образом, можно утверждать, что эффективные резонансные интегралы в гомогенной и гетерогенной среде совпадают, если имеет место равенство:

Σ

p

 

Σ

 

 

 

=

p

.

(8.2.25)

ρП

ρблП

 

 

 

Из соотношений (8.2.24) и (8.2.25) следуют два принципа эквивалентности (иногда называют теоремой эквивалентности) резонансного поглощения:

1. Независимо от природы замедлителя и геометрии блока, ячейки с одинаковым значением величины Σp ρблП имеют одинаковые

резонансные интегралы.

2. Для любой гетерогенной среды можно подобрать такую эффективную гомогенную среду, в которой эффективный резонансный интеграл будет таким же, как и в гетерогенной среде.

Второй принцип эквивалентности часто используется для расчета резонансного поглощения в сложных решетках с двойной гетерогенностью. Так, например, графитовые реакторы типа РБМК и тяжеловодные реакторы CANDU, которые относятся к реакторам канального типа, имеют сложную структуру топливного канала. Топливный канал представляет собой ТВС, как в реакторах водоводяного типа, омываемые теплоносителем. Поэтому такие топливные сборки можно представить в виде решетки твэлов в канале (первый уровень гетерогенности). В свою очередь, топливные каналы образуют решетку каналов (второй уровень гетерогенности). Для расчета резонансного поглощения в такой системе вначале рассматривают решетку на основе топливной сборки в теплоносителе. Заменяя ее гомогенной средой и добавив к величине Σp не-

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

189

8.3.Коэффициент использования тепловых нейтронов

вдвухзонной ячейке

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

 

 

 

 

drr

Eгр

 

 

 

 

 

 

 

 

Σ0a (E) Φ0 (rr, E)dE

 

 

 

θ =

 

 

 

V0

 

 

0

 

 

 

. (8.3.1)

 

rEгр

 

 

 

 

rEгр 1

 

 

 

0

0

 

r

1

r

 

dr

Σa (E) Φ

 

(r , E)dE +

dr Σa (E) Φ

(r , E)dE

 

V0

 

0

 

 

 

V1

0

 

 

 

Введем в рассмотрение вспомогательную величину

q , опреде-

ляемую соотношением:

 

 

 

 

 

 

 

rEгр 1

1

r

 

 

dr

Σa (E) Φ

(r , E)dE

 

q =

V1

 

0

 

 

,

(8.3.2)

drr

Eгр

 

 

 

 

 

 

 

 

Σ0a (E) Φ0 (rr, E)dE

 

 

V0

 

0

 

 

 

 

которая есть относительное значение скорости поглощения нейтронов в замедлителе, взятое по отношению к скорости поглощения нейтронов в топливной зоне. С помощью величины q коэффи-

циент теплового использования примет довольно простой вид

θ =

 

 

1

.

(8.3.3)

1

+ q

 

 

 

190