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

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

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

Несмотря на различие стержней по назначению, действие их на АЗ одинаково, а методы их расчета одни и те же. Их эффективность зависит не только от дополнительного захвата нейтронов, но и от увеличения утечки нейтронов из АЗ, обусловленного деформацией нейтронного поля. На рис. 4.1 качественно показано распределение Ф(r) по оси цилиндрической АЗ с поглощающим стержнем (2) и без него при одинаковой средней мощности (1). Увеличение потока Ф вблизи границы приводит к увеличению вероятности утечки и падению Кэф.

H

Рис.4.1. Зависимость Ф(r) для цилиндрической зоны

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

91

4.3. Материалы и форма поглощающих стержней

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

К легким органам СУЗ относятся стержни, выполненные из сильно поглощающего материала, но малого поперечного размера. Эти стержни собирают в пучки (кластеры) и используют в реакторах ВВЭР.

По аналогии с оптикой, стержень, поглощающий все

падающие

на

него тепловые

нейтроны,

называют

« черным».

Но

плотность потока,

проходящего через

вещество, определяется законом Ф = Фoe

− ∑a x

(рис.4.2).

 

 

 

 

 

Ф

 

 

 

 

 

 

 

 

Фo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

Рис.4.2. Зависимость Ф(х)

92

Стержень считается « черным», если он поглощает 95 % падающих нейтронов. Отсюда определяется толщина слоя, при которой стержень можно считать « черным»:

Фпогл = 1 − e−∑a d = 0,95 ,

Фo

т.е. a d ≈ 3 .

В табл.4.1 приведены значения a и d для некоторых веществ.

Таблица 4.1. Характеристики « черных» стержней

Вещество

-1

d, см

 

a , см

 

Бор

96

0,03

Кадмий

112

~0,03

Железо

0,2

15

Ртуть

15

0,2

При выборе материалов нужно учитывать механическую прочность и коррозионную стойкость. Всем требованиям удовлетворяет бористая сталь. 3%-я бористая сталь может работать более 10 лет без заметного выгорания изотопа бора. Широко применяют карбид бора В4С. Он коррозионно стоек в натрии, гелии, азоте, диоксиде углерода. В воде до температуры 250 °С он слабо взаимодействует с ней, а при 360 °С уже необходима защитная оболочка. Также используют гафний. Он коррозионно стоек в воде, и до 400 °С его механические свойства позволяют его использовать в качестве стержней. В зависимости от конструкции АЗ стержни могут быть круглыми, пластинчатыми, крестообразными и пр. Стержень некруглой формы эквивалентен круглому, если равны их эффективные поверхности, на которых плотность потока тепловых нейтронов обращается в нуль.

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

93

часть выполняется из поглощающего материала, а нижняя представляет собой обычную ТВС. В начале кампании стержни своей поглощающей частью находятся в АЗ и компенсируют избыточную реактивность. По мере выгорания топлива стержни извлекаются, и на их место вводится топливная часть.

Недостатком такой конструкции является то, что подвижные ТВС работают в более тяжелых условиях по сравнению с неподвижными и под АЗ требуется наличие свободного пространства для размещения подвижных ТВС, что приводит к увеличению размера корпуса реактора.

Указанные недостатки можно устранить, если поглотители сделать полыми, а внутри установить неподвижные ТВС.

4.4. Теория поглощающего стрежня. Одногрупповое приближение

Рассмотрим цилиндрический гомогенный реактор без отражателя с эффективными размерами R и Н, внутрь которого по оси введен цилиндрический поглощающий стержень радиусом Rст. Нашей задачей является расчет эффективности управляющего стержня. Задачу будем решать в предположении, что R>>М (длина миграции) и Rст<<R. Кроме того, вначале найдем решение задачи в одногрупповом приближении, когда все нейтроны считаем тепловыми.

Запишем уравнение диффузии

 

 

D Ф − ∑a Ф + S = 0 ,

(4.1)

где S = ka Ф – источник нейтронов.

 

Подставив значение S, получим

 

 

 

 

Ф

 

kФ

 

 

 

Ф

 

+

 

= 0 .

(4.2)

 

 

2

2

 

 

 

L

 

L

 

Здесь L2 =

D

– квадрат длины диффузии.

 

 

 

a

 

 

 

 

 

94

 

 

 

Группируя два последних слагаемых, получим

k − 1

= B

2

или

k

= 1.

(4.3)

L2

 

 

1+ B2L2

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

Подставив (4.3), получим

 

 

Ф + B2Ф = 0 .

 

 

(4.4)

Найдем решение волнового уравнения в отсутствие

стержня. Запишем уравнение (4.4)

в развернутом виде:

1 ∂

(r

Ф

) +

2Ф

+ Bo2Ф = 0 ,

(4.5)

 

 

 

 

 

 

 

 

r r

r

z

2

 

 

 

 

 

где Bo2 – геометрический параметр в отсутствие стержня. Решение уравнения ищем методом разделения

переменных Ф(x,z) = ϕ(r)f(z) :

 

1 1 d

(r

dϕ

) +

 

1 d2f

+ (α2 + β2 ) = 0

,

(4.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ r dr

dr

 

f dz

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где α2 + β2 = Bo2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2f

 

+ β2f = 0 .

 

(4.7)

 

 

 

 

 

 

 

dz2

 

 

 

 

 

 

 

 

 

 

 

dϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

d

(r

) + α2ϕ = 0 .

 

(4.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r dr

 

 

 

 

dr

 

 

 

 

 

 

Решение уравнения (4.7) имеет вид

 

 

 

 

 

 

 

 

 

f = A1 cos βz + B1 sinβz .

 

(4.9)

Так как решение симметрично относительно

середины f(z) = f(z)

 

 

АЗ (z= 0), то В1 = 0.

 

 

На краях при z =

H

 

Ф = 0, т.е.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

π

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

β

H

=

, отсюда β =

.

 

(4.10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

H

 

 

Уравнение (4.8) после подстановки α2 принимает

вид

95

1 d

 

dϕ

 

π2

 

 

 

 

 

(r

 

) + (Bo2

 

)ϕ = 0 .

(4.11)

 

 

 

 

 

 

 

r dr

 

dr

 

H2

 

Это уравнение Бесселя. Его решение выражается

через цилиндрические гармоники – функции Бесселя:

ϕ = AIo (αr) + BYo(αr) .

(4.12)

Здесь Io , Yo соответственно функции Бесселя нулевого

порядка первого и второго рода. Их качественный вид показан на рис.4.3.

Рис.4.3. Зависимость Io , Yo от r

Io (0) = 1 , Yo (0) = −∞ . Так как Yo (0) = −∞ , то второе решение нужно отбросить, т.е. В = 0. Окончательно решение записывается в виде

Ф(r,z) = CIo

(α r) cos

πz

.

 

 

 

 

(4.13)

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Первый корень

Yo (α r) = 0 отвечает α = 2,405, т.е.

экстраполированный

радиус

 

зоны

R =

2,405

.

Таким

 

 

образом, зависимость Ф(r,z) имеет вид

 

 

 

α

 

 

πz

 

 

 

Ф(r,z) = CI

 

(

2,405 r

) cos

,

(4.14)

 

 

 

 

 

 

 

o

 

R

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Bo2 =

π2

+

2,405 2

.

 

 

 

 

 

 

 

(4.15)

H2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2

 

 

 

 

 

 

 

 

 

96

 

 

 

 

 

 

 

 

 

 

 

Поместим теперь по оси зоны центральный поглощающий стержень и определим α и β в этом случае. При наличии стержня, расположенного по всей высоте зоны, граничные условия по z сохраняются, поэтому распределение по оси z остается прежним и β = βo . Уравнение для ϕ(r) тоже остается прежним, только граничные условия будут другими. Внутри стержня поток

будет равен нулю на экстраполированной границе Rэфст . Поэтому ϕ(r) = 0 при r = R и r = Rэфст ,

 

AIo (αR) + BYo (αR) = 0,

(4.16)

 

AI (αRэф ) + BY (αRэф ) = 0.

 

 

 

 

 

o

ст

 

 

 

o

 

ст

 

Из приведенных условий находится параметр α:

 

I (αRэф )

 

=

I (αR)

 

 

 

o

ст

 

o

 

.

(4.17)

 

Y (αRэф )

 

 

 

 

 

 

 

Yo (αR)

 

 

 

 

o

ст

 

 

 

 

 

 

 

Компенсирующая

 

 

 

 

способность

стержня

определяется соотношением

 

 

 

 

 

 

 

 

 

k = k k o

= (α 2 − αo2 )L2 ,

(4.18)

(α 2 + β2 ) = B2 .

 

 

 

 

 

Поскольку Rст << RАЗ ,

 

 

то

вытеснением

среды

поглотителем можно пренебречь, а L2 считать величиной

постоянной. Кроме того,

в больших реакторах Rстэф << R ,

так что изменение α при введении стержня мало:

 

α = αo + α и

α << 1.

 

 

 

 

 

 

 

 

 

 

αo

 

 

 

Функцию Io (αR) разложим в ряд Тейлора:

 

Io (αR) Io (αoR) +

 

dIo

 

αR = Io (αoR) I1(αoR) αR ,

 

d(αoR)

 

 

 

 

 

 

 

 

 

 

Io (αoR) = 0 при αoR = 2,405 , а

 

 

dIo

= −I1 o R) = −0,519 .

d(αR)

Таким образом,

97

Io (aR) = -0,519DaR .

Аналогично

 

Y (aR) @ Y (a R) +

 

dYo

 

DaR = Y (a

o

R) = 0,51,

 

 

 

 

 

o

 

o o

 

 

 

 

d(aoR)

 

 

 

 

 

o

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

так как

dYo

 

= -Y (a

o

R) @ 0,1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d(aR)

 

 

o

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее при aRэф

<< 1

 

 

I (aRэф ) @ 1 , а

 

 

 

 

 

 

 

 

 

ст

 

 

 

 

 

 

 

 

o

 

ст

 

 

 

 

 

 

 

Y (aRэф ) @ -

2

(0,116 + Ln

 

 

R

 

 

) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

o

 

 

ст

 

 

 

 

p

 

 

 

2,405 Rстэф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя полученные оценки в (4.17), получим

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

=

0,519DaR

 

 

 

2

(0,116 + Ln

 

 

 

R

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

0,51

 

 

 

 

 

 

 

 

p

2,405 Rстэф

 

 

 

 

 

 

 

 

 

 

 

и

 

Da =

 

 

 

 

 

 

 

0,51p

 

 

 

 

 

 

 

 

 

.

 

 

(4.19)

 

2 × 0,519(0,116 + Ln

 

 

 

R

 

 

)R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,405 R стэф

 

 

 

 

 

Так как Da = a - ao мало, то формулу (4.18) можно

переписать в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dk = (a2 - ao2 )L2 = (a + ao )(a - ao )L2 = 2a0 DaL2 ,

и,

подставляя

 

α , получим для эффективности стержня

следующую формулу:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dk =

 

 

 

 

 

 

7,5L2

 

 

 

 

 

 

.

 

 

 

 

 

(4.20)

 

 

 

 

R2 (0,116 + Ln

 

 

 

R

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,405 Rстэф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.5.Двухгрупповое приближение

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

98

индексом 2. Для реактора в стационарном состоянии имеем:

D11 - 1 Ф1 + k 2a Ф2 = 0 ,

(4.21)

D22 - 2a Ф2 + 1 Ф1 = 0 ,

(4.22)

D11c - 1c Ф1c = 0 ,

(4.23)

D22c - 2c Ф2c + 1c Ф1c = 0 ,

(4.24)

Первые два уравнения описывают диффузию в АЗ, а вторые – в поглощающем стержне. 1 - сечение

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

Все нейтронные потоки должны быть конечными и положительными. Граничные условия на поверхности стержня следующие. Быстрые нейтроны материалом стержня не поглощаются, т.е. число быстрых нейтронов, входящих в стержень, равно числу нейтронов, выходящих из стержня. Отсюда следует, что на поверхности стержня ÑФ1 = 0 , а на внешней границе Ф1(Rэ ) = 0 . Для тепловых

нейтронов сохраняется прежнее граничное условие. Поток тепловых нейтронов равен нулю на экстраполированной

границе: Rэфст = R - d , где d – длина линейной

экстраполяции, и Ф2(Rэ ) = 0 .

 

 

 

 

 

 

 

 

Решение

 

системы уравнений (4.21) –(4.24) с учетом

граничных условий записывается в виде

 

 

 

 

 

 

Dk =

7,5M2

 

 

 

 

1

 

 

 

 

 

 

.(4.25)

R2

 

 

 

 

 

 

 

 

 

 

 

 

 

t

R

 

t

L t

 

 

 

 

 

 

 

0,116 (1 +

 

) + Ln

 

+

 

Ln

 

 

 

 

 

 

2

эф

2

эф

 

 

 

 

 

L

2,4 Rст

 

L

MRст

 

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

99

4.6. Теория возмущений

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

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

Рассмотрим теорию возмущений в одногрупповом приближении. Пусть до введения стержня реактор находился в стационарном состоянии:

 

k− 1

 

Фo + Bo2Фo = 0 ,

(4.26)

где B2 =

 

. Иногда для уточнения

одногруппового

L2

 

 

 

приближения заменяют L2 M2 . Этот результат можно получить из двухгруппового приближения.

Введем в малый объем реактора V1 поглощающий стержень, что приведет к изменению поглощения нейтронов в этом объеме:

(a )1 = (a )o + δ ∑a ,

(4.27)

Здесь (a )o сечение поглощения вне объема V1; δ ∑a

возмущение в объеме V1.

Материальный параметр В2 в объеме V1 тоже

изменится:

 

 

 

 

B

2

= B2

+ δ B2 .

(4.28)

1

o

 

 

В остальном объеме параметр В2 будет иметь первоначальное значение. Понятно, что после внесения возмущения реактор выйдет из стационарного состояния. Чтобы он остался в критическом состоянии, нужно

100