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

Семенов - учебное пособие по кинетике и регулированию ЯЭУ

.pdf
Скачиваний:
337
Добавлен:
26.05.2014
Размер:
955.2 Кб
Скачать

изменить коэффициент размножения зоны на величину k . Реактор в критическом состоянии описывается

уравнением

 

 

DФ + B2Ф = 0 ,

(4.29)

где B2 = B2

+ DB2

внутри V1 и B2 = B2

+ DB2 вне объема

1

 

o

 

V1; B2 - изменение B2 всей зоны,

необходимое для

компенсации dB2 в объеме V1. Умножим уравнение (4.26) на Ф, а уравнение (4.29) на Фo , вычтем одно из другого и проинтегрируем полученный результат по размеру зоны:

(ФDФo - ФoDФ)dV + (Bo2 - B2 oФdV = 0 .

(4.30)

V

V

 

Для вычисления первого интеграла воспользуемся

тождеством

 

 

 

Ñ(ФÑФo ) = ÑФÑФo + ФÑ2Фo .

 

Отсюда

ФÑ2Фo = Ñ(ФÑФo ) - ÑФÑФo .

(4.31)

Используя (4.31), упростим подынтегральное выражение

ФÑ2Фo - ФoÑ2Ф = Ñ(ФÑФo - ФoÑФ) .

Так как

Ñ(ФÑФo - ФoÑФ)dV =

¶Фo

- Фo

¶Ф

) dS

 

¶n

¶n

 

V

 

V

 

 

 

 

и на поверхности зоны в

 

силу

граничных

условий

Ф1S = Фo1S = 0 , то

интеграл

будет равен

нулю.

Второй

интеграл разобьем на два:

 

 

 

 

 

 

 

(Bo2 - B2 oФdV =

(Bo2 - B2 oФdV + (Bo2 - B2 oФdV .

V

V V1

 

 

 

V1

 

 

 

Подставляя сюда B2 , получим

-DB2ФoФdV - (DB2 + dB2 oФdV = 0

V V1

V1

или

 

- DB2ФoФdV - dB2ФoФdV = 0 ,

V

V1

 

101

 

 

 

ФoФdV

 

отсюда

B2 = −δB2

V1

 

 

.

(4.32)

 

 

 

 

 

 

ФoФdV

 

 

 

 

V

 

 

 

 

Учитывая, что возмущение является малым, можно

положить

Ф Фo . При этом

условии

полученное

выражение примет вид

 

 

 

 

 

 

 

Фo2dV

 

 

B2 = −δB2

V1

 

.

 

(4.33)

 

 

 

 

 

 

Фo2dV

 

V

Отсюда видно, что влияние возмущения пропорционально невозмущенной плотности потока в области возмущения. Так как V1 << V , то в пределах V1

можно Фo считать постоянным:

 

 

 

 

 

 

 

 

B2 = −δB2

 

Ф2V

 

 

 

 

 

 

 

 

 

 

 

 

 

o

1

 

,

(4.34)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

< Фo2 > V

 

 

где < Фo2 >=

1

 

Фo2dV среднеквадратичный поток.

 

 

 

 

 

 

 

V V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для определения изменения реактивности нужно

связать B2 с

k и δ B2 с δ ∑a .

 

 

 

 

 

 

Поскольку

 

 

B2 =

k− 1

,

 

 

 

 

 

 

 

 

L2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

то

 

 

 

 

 

 

 

 

B2 = B2 Bo2 =

k

,

(4.35)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L2

 

 

 

 

 

 

 

 

δ k

 

(k

− 1) δ L2

δ k

 

δ L2

 

 

 

иδ B2 =

 

 

 

 

=

2 B2

2 .

(4.36)

 

2

 

 

4

 

 

 

 

 

L

 

 

 

L

 

 

 

L

 

 

L

aU

 

 

Используя

 

выражения

для

 

k= ηεϕ

и

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L2 =

1

 

,

 

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 ∑tr a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

102

 

 

 

 

 

 

 

δ k = −k δ ∑a ,

δ L2 = −L2 δ ∑a ,

(4.37)

 

a

 

 

 

 

a

 

δ B2 = −

k

 

δ ∑a +

B2L2

δ ∑a .

(4.38)

L2

L2

 

a

a

 

Так как B2L2 << 1 , то

δ ∑a .

 

δ B2 = −

k

 

(4.39)

L2

 

 

 

 

a

 

 

Подставляя (4.35) и (4.39) в формулу (4.34), получим

ρ =

k

=

δ ∑

a

 

Ф2V

 

 

 

 

o

1

.

(4.40)

k

a

 

< Фo2

> V

4.7.Примеры по применению теории возмущений

4.7.1.Эффективность эксцентрично расположенного стержня

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

постороннего тела. Используя закон k ~ Фo2V1 , можно

определить эффективность полностью погруженного стержня на расстоянии r от оси зоны (рис.4.4).

Рис.4.4. Эксцентрично расположенный стержень

103

Его

эффективность

выражается

через

эффективность центрального стержня:

 

 

 

k

=

k(r)

.

 

(4.41)

 

 

 

 

 

 

 

Фo2(0) Фo2(r)

 

 

Отсюда

 

 

 

 

 

 

 

 

 

k(r) =

koФo2(r) .

(4.42)

 

 

 

 

Фo2(0)

 

 

Если в зоне имеется много стержней, находящихся на различном расстоянии от оси, то

k

ko

Фi2(r) .

(4.43)

 

 

Фo2(0) i

 

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

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

Рис.4.5 иллюстрирует интерференцию двух стержней. Кривая 1 показывает зависимость эффективности одного стержня от удаленности его от оси реактора. Кривая 2 показывает половину эффективности

двух стержней 1 k2 (r) в зависимости от расстояния

2 k1(r)

между ними при симметричном расположении относительно оси цилиндрической активной зоны.

104

Рис.4.5. Интерференция двух стержней в зависимости от их расположения

4.7.2. Зависимость эффективности поглощающего стержня от глубины погружения

Оценку эффективности частично погруженного стержня проведем на основе теории возмущений, согласно которой

z

Ф2(z)dz

 

k(z) = ko

o

 

,

(4.44)

 

H

 

 

 

Ф2(z)dz

 

 

 

o

 

 

 

где

ko – эффективность

полностью

погруженного

стержня; Н – высота реактора. Вблизи конца стержня поток искажен, так что, строго говоря, длину стержня нужно увеличить на длину диффузии L, но при L << H концевым эффектом можно пренебречь. Если

Ф = Фo sin( πz ) , то H

105

 

 

 

 

 

z

2

 

 

 

πz

 

 

 

 

 

 

 

 

 

 

 

sin

 

(

 

 

 

 

)dz

 

k(z) =

 

 

 

 

 

 

H

 

 

ko

o

 

 

 

 

 

 

 

.

 

(4.45)

 

 

 

 

 

 

 

 

 

H

2

 

 

 

πz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

(

 

 

 

)dz

 

 

 

 

 

 

 

 

H

 

 

 

 

 

 

 

o

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как sin2 α =

1 − cos

, то

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

z

 

 

k = k

 

 

1

(z

H

 

 

sin

) .

(4.46)

 

 

 

 

 

 

o H

 

 

 

 

 

H

 

Выражение (4.46) перепишем в более удобном виде:

 

 

 

 

 

 

 

 

sin

2πz

 

 

 

k =

k

 

z

(1 −

H

) .

(4.47)

 

 

 

 

 

 

 

 

 

 

 

o H

 

 

2πz

 

 

 

 

 

 

 

 

 

 

 

H

 

При z = H / 2

 

k =

1

 

ko , а при малых z / H

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

sin x x2

1

x6

иk меняется в зависимости от z по закону кубической параболы:

k =

ko

2

(

z

)3 .

(4.48)

6

H

 

 

 

 

 

На рис.4.6 и 4.7 приведены графики интегральной и дифференциальной эффективности управляющего стержня в зависимости от его глубины погружения в активную зону для реактора ВВЭР-1000.

106

Рис.4.6. Интегральная эффективность управляющего стержня

Рис.4.7. Дифференциальная эффективность управляющего стержня

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

107

4.7.3. Стержень типа нейтронной ловушки

Эффективность стержня может быть увеличена, если наряду с тепловыми нейтронами он поглощает и быстрые. Если стержень выполнить полым с оболочкой из сильного поглотителя, а внутри полости расположить хороший замедлитель, быстрые нейтроны, проходя сквозь оболочку, будут замедляться внутри стержня и поглощаться оболочкой (нейтронная ловушка). При этом внутренний диаметр стержня должен удовлетворять условию Rвн > τ . Здесь τ − возраст нейтронов. Такие

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

4.7.4. Калибровка управляющих стержней

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

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

108

скорости нарастания реактивности реактор становится критическим во время пуска.

Существуют различные методы калибровки стержней. Один из основных методов основывается на уравнении Нордхейма:

ρ =

τ

+

βiτi

 

 

 

.

(4.49)

T kэф

T + τi

Здесь Т – установившийся период реактора.

Этот метод используется при незначительной мощности реактора, когда можно пренебречь температурными эффектами. Реактор приводится в критическое состояние при определенном положении всех стержней. Далее калибруемый стержень (группа) смещается в направлении положительной реактивности, реактор выдерживается в течение нескольких поколений нейтронов и измеряется установившийся период. Реактивность определяется по формуле (4.49). Если Т велик, т.е. τ << T , то первым слагаемым можно пренебречь. Далее реактор приводится в критическое состояние, и процедура повторяется. При этом геометрия системы должна оставаться неизменной. Компенсация избыточной реактивности наилучшим образом может быть осуществлена за счет равномерного распределения поглотителя по всему объему вследствие его растворения в замедлителе или теплоносителе. Если количество поглотителя, равномерно введенного в объем реактора, известно, то этим методом можно осуществить калибровку. Калибровочную кривую, определяющую реактивность стержня в зависимости от его положения, можно получить последовательным добавлением поглотителя. При этом реактивность стержня определяется соотношением

ρ

k

~

f

a ,

(4.50)

 

 

k

f

a

 

где f – коэффициент использования тепловых нейтронов; a – полное сечение поглощения.

109

Калибровка по методу изменения периода или по методу введения поглотителя может проводиться одновременно.

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

Как следует из (1.41),

 

N

=

 

β

 

 

N0

β +

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

и

ρ

=

N0

− 1 .

(4.51)

β

 

 

 

 

N

 

При этом необходимо, чтобы сигнал от датчиков нейтронного потока был пропорционален средней плотности потока нейтронов.

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

110