- •Севастопольский институт ядерной энергии и промышленности
- •Основы теории ядерных реакторов Курс для эксплуатационного персонала аэс
- •Содержание
- •Перечень сокращений
- •Тема 1.
- •1.1. Строение вещества
- •1.2. Строение и характеристики атомов
- •Атомная теория раскрывает физический смысл этих характеристик в следующих основных положениях:
- •1.3. Строение ядер и свойства ядерных сил
- •1.4. Энергия связи и устойчивость ядер атомов
- •1.4.5. Энергия связи, приходящаяся на один нуклон ядра
- •1.5. Закономерность и характеристики радиоактивного распада
- •Тема 2 нейтронные ядерные реакции
- •2.2. Особенности реакции деления и их практическое значение
- •2.3. Основные характеристики нейтронных полей
- •2.4. Скорости нейтронных реакций и их характеристики
- •Тема 3 критичность реактора и условия её реализации
- •3.1. Условия осуществления критичности реактора
- •3.1.2. Эффективный коэффициент размножения и реактивность реактора
- •3.2. Нейтронный цикл в тепловом ядерном реакторе.
- •3.2.2. Нейтронный цикл и характеристики его физических процессов
- •4.1. Ядерное топливо.
- •4.2. Замедлитель.
- •4.3. Теплоноситель
- •4.4. Параметры структуры активных зон гетерогенных эяр.
- •Тема 5 замедление нейтронов в реакторе и его размножающие свойства
- •5.1. Общие начальные рассуждения
- •Вероятность избежания утечки замедляющихся нейтронов - это доля нейтронов, избежавших утечки из активной зоны при замедлении, от всех нейтронов поколения, начавших процесс замедления в активной зоне.
- •5.2. Характеристики замедляющих свойств веществ
- •5.3. Возраст нейтронов в среде
- •Величину, обратную величине транспортного смещения
- •Возраст нейтронов с энергией е - это шестая часть среднего квадрата пространственного смещения нейтрона в среде при замедлении от начальной энергии Ео до данной энергии е.
- •5.4. Уравнение возраста Ферми и его решение
- •5.5. Вероятность избежания утечки замедляющихся нейтронов
- •Спектр замедляющихся нейтронов Ферми в гомогенной непоглощающей среде
- •5.7. Время замедления нейтронов в среде активной зоны
- •Краткие выводы
- •Тема 6. Диффузия и размножающие свойства теплового реактора
- •6.1. Закон диффузии тепловых нейтронов и длина диффузии
- •6.2. Скорость утечки тепловых нейтронов из единичного объёма активной зоны
- •6.3. Волновое уравнение, уравнение критичности реактора и величина вероятности избежания утечки тепловых нейтронов
- •6.4. Геометрический параметр цилиндрического реактора без отражателя и поле тепловых нейтронов в нём
- •Краткие выводы
- •7.1. Константа
- •7.2. Коэффициент использования тепловых нейтронов
- •7.2.6. Зависимости величины от определяющих её факторов.
- •Краткие выводы
- •Тема 8 уран-238 и размножающие свойства реактора
- •8.1. Коэффициент размножения на быстрых нейтронах
- •8.1.2. Величина в цилиндрическом блоке из природного металлического урана.
- •8.2. Вероятность избежания резонансного захвата
- •Тема 9 критические размеры и нейтронное поле в реакторе с отражателем
- •9.1. Отражатель теплового реактора
- •9.2. Эффективная добавка (э)
- •9.3. Геометрический параметр и поле тепловых нейтронов в гомогенной цилиндрической активной зоне с отражателем
- •9.4. Особенности нейтронного поля в гетерогенном реакторе
- •9.5. Показатели неравномерности нейтронного поля в реакторах
- •Тема 10 температурные эффекты реактивности реактора
- •Температурный эффект и температурный коэффициент реактивности
- •Температурный эффект реактивности реактора
- •Три характерных для ввр типа кривых тэр
- •Температурный коэффициент реактивности реактора (ткр)
- •Условие устойчивости работы энергетического реактора на мощности
- •10.3. Чем определяется форма кривой тэр реактора?
- •Условные составляющие тэр и ткр
- •Мощностной тэр (ткр) реактора
- •Тэр и ткр теплоносителя
- •Раздел 3 кинетика реактора
- •Тема 11 элементарная кинетика теплового реактора
- •10.1. Элементарное уравнение кинетики реактора
- •Среднее время жизни поколения нейтронов в тепловом реакторе
- •Следовательно, время жизни запаздывающих нейтронов любой группы
- •11.3. Период реактора, период удвоения мощности и их взаимосвязь
- •Тема 12 кинетика реактора с учётом запаздывающих нейтронов
- •Система дифференциальных уравнений кинетики реактора с учётом
- •Уравнение обратных часов.
- •Переходные процессы при сообщении реактору отрицательной
- •Переходные процессы при сообщении реактору положительных реактивностей
- •Особенности переходных процессов при сообщении реактору малых и больших реактивностей
- •Как управляют реактором на малых уровнях мощности?
- •Тема 13 основы кинетики подкритического реактора при его пуске
- •Источники нейтронов в подкритическом реакторе
- •Что это за источники?
- •Устанавливающаяся в подкритическом реакторе плотность нейтронов
- •Переходные процессы при изменениях степени подкритичности реактора
- •Учитывая, что отношение начальной и конечной плотностей нейтронов
- •Время практического установления подкритической плотности
- •Процедура ступенчатого пуска и ядерная безопасность реактора
- •Краткие выводы
- •Раздел 4. Изменения запаса реактивности при работе реактора
- •Тема 14.
- •Понятия общего и оперативного запаса
- •Тема 15 уменьшение запаса реактивности с выгоранием ядерного топлива
- •15.2. Энерговыработка реактора
- •15.4. Основные характеристики выгорания
- •Тема 16 уменьшение запаса реактивности за счёт шлакования ядерного топлива
- •Кинетика роста потерь запаса реактивности за счёт шлакования
- •Тема 17 рост запаса реактивности с воспроизводством ядерного топлива
- •17.2. Система дифференциальных уравнений воспроизводства плутония-239
- •Рост запаса реактивности с воспроизводством плутония-239.
- •17.4. Коэффициент воспроизводства ядерного топлива
- •Тема 18 использование выгорающих поглотителей
- •18.1. Характеристики наиболее распространённых выгорающих поглотителей
- •18.2. Факторы, определяющие скорость выгорания вп
- •18.4. Кривая энерговыработки активной зоны реактора
- •Тема 19 отравление реактора ксеноном
- •Отравления реактора ксеноном
- •Стационарное отравление реактора ксеноном.
- •19.3. Переотравление после останова реактора («йодная яма»)
- •Переотравления реактора ксеноном после изменения уровня мощности
- •19.5. Расчёт изменений потерь реактивности за счёт переотравлений реактора.
- •Тема 20 отравления реактора самарием-149
- •20.1. Схема образования-убыли 149Sm и дифференциальные уравнения отравления реактора самарием
- •20.1. Схема образования и убыли самария-149 и сопутствующих продуктов деления и их распада
- •20.2. Потери реактивности при стационарном отравлении реактора самарием
- •20.3. Закономерность роста потерь реактивности от отравления самарием до выхода реактора на стационарный уровень отравления.
- •20.4. Нестационарное переотравление реактора самарием после останова («прометиевый провал»)
- •20.5. Переотравление самарием после пуска длительно стоявшего реактора
- •20.6. Нестационарное переотравление реактора самарием после перевода реактора на более высокий или более низкий уровень мощности
- •Раздел 5.
- •Действие вводимого в активную зону стержня-поглотителя
- •Характеристика положения стержня-поглотителя в активной зоне
- •Понятия об интегральной и дифференциальной эффективности
- •Эффективный радиус стержня-поглотителя
- •Физический вес центрального стержня-поглотителя полной длины
- •21.6. Физический вес нецентрального подвижного поглотителя
- •Характеристики поглотителей – кривые интегральной и дифференциальной эффективности
- •Изменение реактивности реактора при перемещении стержня
- •Особенности характеристик укороченных поглотителей
- •Интерференция подвижных стержней-поглотителей
- •21.11. Простейшие методы градуировки подвижных поглотителей
- •Тема 22 борное регулирование ввэр
- •22.1. Сущность борного регулирования
- •22.2. Характер изменения концентрации борной кислоты в первом контуре
- •Эффективность борной кислоты
- •Факторы, определяющие величину дифференциальной эффективности борной кислоты
- •Тема 23 расчётное обеспечение ядерной безопасности ввэр при его эксплуатации
- •Расчёт пусковой критической концентрации борной кислоты
- •Расчёт предельно допустимого расхода подпитки первого контура чистым дистиллатом при пуске ввэр
- •Время снижения концентрации борной кислоты до заданной величины
- •Расчёт безопасного значения стояночной концентрации борной кислоты
- •23.5. Расчёт времени подпитки первого контура концентрированным раствором борной кислоты до достижения безопасной стояночной концентрации
- •Литература
Уравнение обратных часов.
а) Уравнение обратных часов как характеристическое уравнение системы дифференциальных уравнений кинетики реактора. Развёрнутый вид уравнения:
свидетельствует о том, что это - алгебраическое уравнение седьмой степени относительно Т. Для того, чтобы понять это, достаточно мысленно представить, что получится в числителе выражения правой части уравнения после приведения его к общему знаменателю: седьмая степень уравнения становится очевидной. А это значит, что уравнение обратных часов в самом общем случае должно иметь семь корней.
В разделе математики “Решение дифференциальных уравнений” говорится, что вид, величины и знаки корней характеристического уравнения определяют вид решения дифференциальных уравнений. В частности, если характеристическое уравнение имеет действительные корни, то решение дифференциального уравнения (или системы дифференциальных уравнений) имеет экспоненциальный вид. Это побуждает нас заняться анализом корней уравнения обратных часов.
Но вначале следует уточнить одну немаловажную деталь: отыскивая с самого начала именно экспоненциальное решение для дифференциального уравнения изменения плотности нейтронов в виде одной экспоненты - n(t) = no exp (t/T), - мы невольно совершили небольшой просчёт, который следует исправить по ходу дела. Если уравнение обратных часов имеет семь корней, то общее решение системы дифференциальных уравнений кинетики реактора будет не одной экспонентой, а будет представлять собой сумму семи экспонент, показатели которых определяются величинами этих семи корней уравнения обратных часов:
или в более краткой форме:
(12.19)
где То, T1, T2, ... , T6 - значения семи корней уравнения обратных часов, а Ao, A1, A2, ... , A6 - величины постоянных интегрирования, находимые путём подстановки в общее решение (12.19) конкретных начальных условий.
Знаки корней уравнения обратных часов наиболее наглядно видны, если показать его решение в графическом виде (см. рис.12.3).
r
1/T6 1/T5 1/T4 1/T3 1/T2 1/T1 1/T0
-l6 -l5 -l4 -l3 -l2 -l1 0 1/T
Рис,12.3. График зависимости корней уравнения обратных часов при положительных и
отрицательных реактивностях разной величины.
На этом графике показано решение уравнения обратных часов в зависимости не от величины самого периода Т, а от обратной ему величины 1/Т: так удобнее выполнять решение уравнения аналитически.
Как видим, функция r = f(1/T) имеет шесть точек разрыва (второго рода), и именно благодаря этой разрывности отдельные корни уравнения обратных часов отображаются на графике достаточно наглядно: области изменения каждого из семи корней по оси 1/T лежат между соответствующими точками разрыва; например, нулевой обратный корень 1/То лежит правее первой точки разрыва (- l1), первый обратный корень 1/Т1 - между первой и второй точками разрыва (-l1 и -l2), второй обратный корень 1/T2 - между второй и третьей точками разрыва (-l2 и -l3 ) и так далее; значения последнего, седьмого, обратного корня 1/Т7 - располагаются левее последней, шестой, точки разрыва (-l6).
При этом знаки всех семи корней уравнения обратных часов определяются самым наглядным образом: в точках пересечения соответствующих участков графика с горизонтальной прямой, отсекающей на оси ординат рассматриваемое значение реактивности r.
Из графика видно, что, если величина сообщаемой реактору реактивности положительна, то нулевой обратный корень 1/То , а, значит, и сам корень То, - положителен (так как располагается в правой полуплоскости, правее оси О - r). Остальные шесть корней (Т1 ¸ Т6) уравнения обратных часов - отрицательны (лежат в левой полуплоскости). Если же величина сообщаемой реактору реактивности отрицательна, то все семь корней уравнения обратных часов отрицательны (лежат в левой, отрицательной, полуплоскости).Что касается величин самих корней, то они, как следует из графика, определяются только величиной сообщаемой реактору реактивности r.
Теперь о знаках постоянных интегрирования (Ао ¸ А6). Здесь не приводится полный (и очень громоздкий) аналитический вывод общего выражения для любой из постоянных интегрирования, которое имеет вид:
(12.20)
Квадраты в выражении в скобках под знаком суммы, независимо от знака корня Тi , всегда имеет положительный знак, и, поскольку все остальные величины знаменателя положительны, то весь знаменатель (12.20) - всегда положителен. А раз так, то знак постоянной интегрирования Аi всегда определяется знаком произведения rТi в числителе правой части этого выражения. То есть, если нулевой корень То при r > 0, как говорилось выше, положителен (То > 0), то произведение rТо > 0, а, следовательно, Ао > 0. При отрицательной же величине реактивности (r < 0) произведение отрицательного нулевого корня То на отрицательную величину реактивности r даёт положительную величину произведения rТо, и, следовательно, величина нулевой постоянной интегрирования будет иметь положительный знак. Проделав такой микроанализ со всеми величинами постоянных интегрирования, можно прийти к общему выводу:
при r>0: То > 0 и Ао > 0, а остальные корни (Т1¸Т6)<0 и (А1¸А6) < 0 (12.21)
при r<0: все 7 корней (То¸Т6)<0, а постоянные интегрирования (Ао¸А6)>0 (12.22)
б) Самостоятельное практическое значение решения уравнения обратных часов.
Анализ знаков корней уравнения обратных часов и постоянных интегрирования безусловно важен, так как он даёт возможность выявить качественную структуру и дать физическое толкование характеру переходных процессов при сообщении реактору реактивности того или иного знака.
Но значение уравнения обратных часов не исчерпывается только этим. Забегая немного вперёд, заметим, что из семи корней уравнения обратных часов старший корень То (старший - в математическом понимании этого слова, то есть - наибольший по абсолютной величине) имеет простое физическое толкование: это - тот самый установившийся период, определяющий “чисто экспоненциальное” изменение плотности нейтронов в реакторе в развитой стадии переходного процесса n(t). (Почему это так, станет ясно немного позже).
,% r, %
а) График
0.3
0.3
0.2
0.2
0.1
0 10 20 30 40 50 Т2 с
б) Таблица
r, % |
0.2925 |
0.2152 |
0.1745 |
0.1482 |
0.1294 |
0.1152 |
0.1039 |
0.0948 |
0.0872 |
Т2, с |
5 |
10 |
15 |
20 |
25 |
30 |
35 |
40 |
45 |
Рис. 12.2. Две наглядные формы отображения взаимосвязи реактивности реактора r и установившегося периода удвоения мощности реактора Т2 , вытекающие из решения уравнения обратных часов.
Эту величину установившегося периода, как уже говорилось, легко измерить практически с помощью самого обычного секундомера. Вместе с тем, величина установившегося периода (То) для конкретного реактора в рассматриваемый момент кампании определяется только величиной сообщённой реактору реактивности r, следовательно, уравнение обратных часов для конкретного реактора (с конкретной величиной bэ) устанавливает жёсткую однозначную взаимосвязь величин реактивности r и установившегося периода То (или r - с величиной установившегося периода удвоения мощности реактора Т2, которая, как мы знаем, пропорциональна величине установившегося периода То).
А это значит, что по величине измеренного установившегося периода удвоения мощности можно находить величину сообщённой реактору реактивности, и, наоборот, - по величине сообщённой реактору реактивности можно предсказывать, с каким установившимся периодом удвоения будет происходить установившийся разгон (или спад) мощности ядерного реактора. Это, согласитесь, практически очень важно для оператора реакторной установки. Эту взаимосвязь можно занести в программу компьютера, можно выразить в форме таблицы или в форме графика.
Пользуясь приведенными таблицей или графиком, оператор имеет возможность быстро оценить величину реактивности реактора по измеренному периоду удвоения мощности или предсказать величину установившегося периода разгона реактора по величине реактивности, которую он собирается сообщить реактору.