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

Система автоматического управления Митенков Ф.М., Чирков В.А

..pdf
Скачиваний:
50
Добавлен:
20.01.2021
Размер:
5.49 Mб
Скачать

рого на долл ниже критического значения, скорость прироста реактивности

вызовет увеличение плотности нейтронов с периодом

сек.

Практическое значение скачкообразной реакции плотности нейтронов

при внезапном изменении

заключается в возможности организации защиты

на случай неконтролируемой большой скорости перемещения регулирующих органов.

Из уравнения (3.74) следует, что контроль по периоду может обеспечить безопасность не только в отношении установившегося периода, но и в отношении такой скорости возрастания реактивности, которая может привести реактор к мгновенной критичности в очень короткое время.

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

,

где индекс "0" - относится к исходному равновесному значению. Подстановка этого выражения в (3.29) приводит к уравнению (S* не учитывается):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.77)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

. Вводя

интегрирующий множитель получим общее решение уравнения (3.77)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(3.78)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для постоянной скорости прироста реактивности

 

это уравнение при-

водится к виду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

(3.79)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- значение при . Уравнение (3.79) справедливо вплоть до значения ρ, соответствующего мгновенной критичности.

Для реактора, первоначально находившегося в подкритическом состоянии, можно использовать уравнения (3.77) и (3.78), отсчитывая время от момента достижения критичности, при котором

60

(3.80)

Величина

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

критичностью .

Малые значения

. Для низких скоростей изменения реактивности

скорость генерации запаздывающих нейтронов

любой

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

(3.81)

где - начальный уровень плотности нейтронов; - установившийся обратный период, соответствующий реактивности,

имеющийся в момент времени t и

. (3.82)

Из уравнения (3.81) следует, что величина обратная периоду реактора, определяется выражением

.

(3.83)

Соотношение (3.82) можно представить в одногрупповом приближении, пренебрегая временем жизни мгновенных нейтронов, следующим образом

 

.

(3.84)

 

Тогда из уравнений (3.75) и (3.84) следует:

,

что находится в согласии с выражением (3.74).

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

61

ф (реактивность ρ),

ГЛАВА IV

ОБРАТНЫЕ СВЯЗИ В РЕАКТОРЕ

§4-1. Общие сведения.

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

Из-за этого изменяются ядерные сечения поглощения, рассеяния, деления

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

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

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

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

Примером отрицательной обратной связи является ухудшение свойств воды как замедлителя в водо-водяных реакторах при увеличении температуры.

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

Реактивность в уравнениях кинетики можно представить в виде:

,

(4.1)

где введенная извне реактивность в может быть задана, а ос

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

внутренними процессами в данном реакторе и характером возмущения.

62

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

Для исследования динамического поведения реактора наибольший интерес представляют температурные обратные связи (из-за их "быстродействия").

§ 4-2. Температурные обратные связи.

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

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

:

 

 

 

 

,

 

(4.2)

 

 

 

 

где - температура i-го компонента активной зоны.

 

При это изменения реактивности

 

, обусловленное изменением тем-

пературы i-го элемента от до , будет:

 

 

 

 

 

 

 

 

 

.

(4.3)

 

 

 

Получение аналитического выражения для температурного коэффици-

ента реактивности

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

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

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

Гомогенные реакторы. Рассмотрим точечную модель гомогенного реак-

тора без отражателя. Используя взаимосвязь и ф , выражение для

можно

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

 

 

 

 

ф

.

(4.4)

 

 

 

Поскольку для всех практически интересных случаев ф близко к 1, то хорошим приближением для инженерных расчетов может служить выражение

ф

.

(4.5)

 

ф

63

В рассматриваемом приближении

ф

,

(4.6)

где η - среднее число нейтронов, испускаемых при захвате одного теплового нейтрона ядрами топлива;

ε- коэффициент размножения на быстрых нейтронах;

-вероятность избежать резонансного поглощения;

-коэффициент теплового использования;

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

ветственно.

В соответствии с (3.5) и (3.6) получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(4.7)

 

 

 

 

 

 

 

 

 

 

 

 

Обозначая

 

 

 

, можем записать

 

 

 

 

 

 

 

 

 

,

(4.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*.

(4.10)

Полученные выражения (4.7) (4.10) позволяет определить температур-

ный коэффициент реактора

по известным температурным коэффициентам

физических параметров этого реактора.

 

Температурный коэффициент для величины η обусловлен различной за-

висимостью сечений деления (

) и радиационного захвата ( ) для делящихся

материалов от скорости нейтронов. Температурные коэффициенты

при

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

 

 

 

 

 

 

 

 

 

Таблица 4.1

 

 

 

 

 

 

 

Изотоп

U235

Pu239

 

U233

 

Естественный

 

 

уран

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффициент размножения на быстрых нейтронах ε в гомогенных тепловых реакторах не зависит от температуры, поэтому температурный коэффи-

циент

равен нулю.

 

Температурная зависимость вероятности избежать резонансного захвата

 

64

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

ский эффект Допплера, т. е.

.

 

 

 

Коэффициент теплового использования для гомогенного реактора опре-

деляется выражением

 

 

 

 

 

 

 

,

(4.11)

 

 

где индексы "г" и "з" относятся соответственно к топливу и замедлителю; N - ядерная плотность.

Из (3.11) следует, что

.

(4.12)

Учитывая зависимости ядерной плотности N и сечения поглощения

от темпе-

ратуры, выражение (4.12) при некоторых упрощающих посылках можно привести к виду

 

 

 

 

 

,

(4.12а)

где

- фактор, учитывающий отношение от закона

.

 

 

Значения коэффициента f обычно заключены в интервале

. В

табл. 4.2 приведены характерные значения

при температуре

.

 

 

 

 

 

 

Табл. 4.2.

 

 

 

 

 

 

 

 

Изотопы

 

 

 

 

 

 

 

 

 

 

 

 

 

U235

 

 

 

 

 

 

Pu239

 

 

 

 

 

 

U238

 

 

 

 

 

 

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

Pu239, коэффициент

мал.

 

 

 

 

Соответствующая оценка вклада параметров L2, τ, B2 в температурный

коэффициент реактивности приводит к следующим выражениям вместо (4.9) и

(4.10):

,

(4.13)

65

(4.14)

Здесь - температурный коэффициент объемного расширения; - константа 0.1, характеризующая температурную зависимость коэффициента

диффузии

.

(4.15)

Если коэффициент

мал, а

, то

, (4.16)

где использовано общепринятое обозначение

Если , , и малы, то , и температурный коэффициент реактивности гомогенного реактора без отражателя описывается выражением

 

 

 

 

.

(4.17)

 

 

 

 

Анализ выражения (4.17) показывает, что

отрицателен и уменьшается

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

эффициента

будет

возрастать. Численные значения температурного коэффи-

циента реактивности

для гомогенных тепловых

реакторов изменяются от

- для больших графитовых реакторов до

- для небольших

реакторов с водяным замедлителем.

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

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

66

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

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

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

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

Вгетерогенном реакторе коэффициент размножения на быстрых нейтронах е может изменяться с температурой по двум причинам: тепловое расширение ТВЭЛ изменяет вepoятность вылета быстрых нейтронов из топлива (этот эффект практически мгновенный; с изменением температуры топлива несколько изменяется пространственное распределение потока тепловых нейтронов в ячейке, а значит, и пространственное распределение источников быстрых нейтронов, что снова приводит к изменению вероятности вылета быстрых нейтронов из топлива. Выравнивание потока тепловых нейтронов в объеме элементарной ячейки с ростом температуры замедлителя происходит в основном за счет увеличения L2, следовательно, этот эффект изменения ε можно считать запаздывающим. Оба эффекта обычно очень малы.

Вгетерогенном реакторе с низкообогащенным топливом влияние температуры на вероятность избежать резонансного поглощения р (эффект Допплера) весьма существенно. Приближенная формула для р в гетерогенном реакторе имеет вид:

67

,

(4.18)

где - плотность ядер резонансного поглотителя в объеме топлива;

,- объемы топлива и замедлителя соответственно;

-замедляющая способность замедлителя;

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

 

 

 

 

 

,

(4.19)

где - температура топлива в К;

;

 

 

 

 

 

- линейная функция от отношения поверхности топлива к его массе.

(4.20)

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

успевает заметно измениться, принимая, что

 

 

 

 

 

 

 

 

 

 

, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

(4.21)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

(4.22)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку

, то заключаем, что мгновенный температурный

 

 

 

 

 

 

коэффициент отрицателен и обратно пропорционален

 

.

 

 

 

 

Запаздывающие слагаемые температурного коэффициента

обусловле-

ны зависимостью от температуры произведения

(4.18).

 

 

 

 

Учитывая, что

, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(4.23)

 

 

 

 

 

 

 

 

 

 

 

Так как

почти не зависит от ,

. В случае твердого замедлителя в

ячейке

. Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

, т. е. совпадает с мгновенной

 

 

 

 

 

 

 

 

 

 

 

 

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

, то

,

(4.24)

где - объемный коэффициент расширения замедлителя. Следовательно

68

. (4.25)

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

Для коэффициента теплового использования в гетерогенном реакторе справедлива формула:

,

(4.26)

где - коэффициент проигрыша, определяемый выражением

.

(4.27)

Величины с чертой - средние по объему замедлителя и топлива потеки нейтронов.

. (4.28)

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

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

изменению f. В случае твердого замедлителя

з з

. Поскольку замедли-

тель является поглотителем типа

, получим

 

 

 

 

 

.

 

(4.29)

Уравнение (4.28) принимает вид

 

 

 

 

 

 

 

.

(4.30)

Для жидкого замедлителя с фиксированным объемом з, учитывая (4.24), получим

 

.

(4.31)

В случае урана коэффициент

, как правило, мал.

Расчет

очень сложен, и знак его всегда отрицателен (увеличение Т приводит к выравни-

69