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

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

Подставляя сюда (1.43), найдем величину, обратную периоду реактора:

Сопоставляя первое слагаемое в (1.45) со вторым, видим, что при b ¹ r имеет место оценка

Таким образом, в первый момент после возмущения реактора скачком реактивности реакция реактора обусловлена мгновенными нейтронами:

С ростом ρ увеличивается скорость изменения концентрации нейтронов и мощности реактора:

Она равна нулю, когда реактор находится в критическом состоянии, и по мере возрастания реактивности она тоже растет.

1.6. Мгновенно-критический реактор

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

21

мгновенно-критического реактора, т.е. ρ = β. Подставляя это условие в (1.52), получим

22

Таким образом, при мгновенной критичности мощность в течение 1÷2 с увеличится до значения

и далее будет нарастать с периодом Т1. Для приведенной

выше оценки это составит n = 13,5 . no

1.7. Введение отрицательной реактивности

При введении отрицательной реактивности характер решения сохраняется, только в формуле (1.43) нужно у реактивности изменить знак:

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

Например, при мгновенном введении реактивности ρ = -3β, характерной для аварийных стержней, мощность скачком упадет на 75 %, однако 25 % тепловыделения, обусловленного запаздывающими нейтронами, остается. Оно становится малым в течение нескольких минут

(рис.1.5).

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

23

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

Эффективность стержня в зависимости от глубины погружения в зону изменяется по закону ∫ sin2 xdx , т.е.

вначале погружения эффективность стержня незначительна. Точный расчет эффективности стержня в зависимости от времени падения показывает, что в течение первых 200 мс практически никакая реактивность не вводится. Такое запаздывание введения реактивности эквивалентно появлению дополнительного времени задержки в цепях АЗ. В итоге после срабатывания АЗ на уровне 1,2 Nо мощность будет продолжать увеличиваться в течение времени (τср+200) мс с тем периодом, который был у реактора при выходе в критическое состояние по мгновенным нейтронам. За это время мощность реактора может превысить номинальную мощность в несколько сотен раз. Оставшаяся после полного введения стержней медленно изменяющаяся мощность (25 %) может быть на два порядка выше номинальной. Опасность представляет не большая достигнутая мощность, а освобожденная в результате этого процесса энергия. Действительно, в результате рассмотренного быстрого процесса теплоотвода от реактора практически не будет, т.е. температура зоны будет увеличиваться по закону

24

Cp o

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

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

Для ВВЭР γ = 1 + 20β2 1,05 ÷1,1.

Рис.1.5. Зависимость мощности реактора от времени при отрицательном скачке реактивности:

1 – вклад всех нейтронов; 2 – вклад мгновенных нейтронов

25

1.8. Кинетика реактора при линейном изменении реактивности

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

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

ρ = αt .

Анализ кинетики реактора в этих условиях проводим качественно, представляя реактор с сосредоточенными параметрами (“точечный” реактор) с учетом одной группы запаздывающих нейтронов:

Записанную систему уравнений легко свести к уравнению второго порядка. Для этого сложим исходные уравнения:

а первое уравнение продифференцируем по времени:

Начальные условия имеют вид: при t = 0 N = Nо,

dN = 0 . Последнее условие следует из того, что до dt

26

переходного процесса реактор был критичен и dC = 0 . dt

Отсюда βτ No = λCo , k(0) = 0 и из (1.59) следует, что

( dN)o = 0 . dt

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

Приведем качественный анализ уравнения (1.63).

Если бы отсутствовало слагаемое с dN , то уравнение dt

имело бы вид

d2N − b(t)N = 0 . dt2

Решением такого уравнения, как известно, является экспонента. По аналогии с уравнением затухающих колебаний

d2x + γ dx + ω2 x = 0 , dt2 dt

наличие в уравнении (1.63) слагаемого, содержащего dN , dt

играет роль тормозящей силы, если это слагаемое входит в уравнение с положительным знаком. Поскольку с течением времени это слагаемое может переменить знак, то вместо торможения получится разгон и N(t) начнет расти очень круто. Нас интересует решение этого уравнения до того момента, когда реактор станет критичным по мгновенным нейтронам, т.е. до t = β/α. Из проведенного качественного рассмотрения ясно, что решение уравнения (1.63) будет близко к экспоненциальному:

27

Результаты численного интегрирования уравнения (1.62), представленные на рис.1.6, подтверждают приведенный здесь качественный анализ.

Рис.1.6. Зависимость мощности реактора от времени:

β=0,0064, λ=0,077 1/с, τ+10-3 с, α=5 .10-4 1/с

1.9. Кинетика холодного реактора с источником

Для пуска холодного реактора необходимы начальные нейтроны, чтобы осуществить цепной процесс. Эти нейтроны могут быть обусловлены космическим излучением. Другим источником является спонтанное деление урана. Известно, что 1 т природного урана дает около 1,5×104 1/с. Если реактор запускается после останова, то в нем есть уже свободные нейтроны. Их появление обусловлено фотонейтронными источниками и излучением осколков деления.

Однако удобнее и безопаснее пользоваться внешними источниками, помещаемыми в АЗ. Дело в том, что при пуске реактора с малой мощностью внутреннего источника чувствительность пусковой аппаратуры может оказаться недостаточной для контроля выхода на критическое состояние и перехода в надкритическое состояние. Безопасность такого пуска “вслепую” обеспечивается только выбранной программой высвобождения реактивности: шагом освобождения

28

реактивности и временем выдержки между шагами. Нужно, чтобы при переходе через критическое состояние мощность достигла контролируемого уровня раньше, чем будет высвобождена опасная реактивность. Чтобы избежать пуска реактора “вслепую”, используют внешний источник нейтронов. Использование внешнего источника дает такой же эффект, как и повышение чувствительности пусковой аппаратуры. Применяют (Ra, Po, Pu-Be) источники с мощностью Q = 106÷107 1/с.

Рассмотрим пуск холодного реактора на модели “точечного” реактора в одногрупповом приближении.

называемая подкритичностью реактора; q – объемная плотность источника нейтронов. При >> β

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

Решение этого уравнения имеет вид

т.е. реактор ведет себя как размножающая среда с

1

коэффициентом умножения . Установившееся

ρпод

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

29

При достижении критического состояния (kэф = 1) плотность нейтронов будет изменяться по линейному закону

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

Если необходимо учитывать

запаздывающие нейтроны, решая систему уравнений (1.65) и (1.66). Однако при решении практических задач используют решение (1.68), заменив время жизни мгновенных нейтронов τ на среднее время жизни с учетом запаздывающих нейтронов

На практике для оценки переходного процесса используют понятие времени установления подкритической плотности потока. Оно равно времени достижения потоком значения 0,95Фуст. Используя это уравнение, найдем (1.68) при no = 0:

tуст ρ3τ . (1.72)

под

Из уравнения (1.69) следует, что скорость увеличения мощности в подкритическом состоянии определяется подкритичностью и скоростью ее изменения

т.е. при постоянной скорости высвобождения реактивности скорость роста мощности обратно пропорциональна ρпод2 .

Основным условием обеспечения ядерной безопасности при пуске реактора является выведение ЯР в контролируемое надкритическое состояние с допустимым периодом разгона и надежным исключением

30