- •Севастопольский институт ядерной энергии и промышленности
- •Основы теории ядерных реакторов Курс для эксплуатационного персонала аэс
- •Содержание
- •Перечень сокращений
- •Тема 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. Расчёт времени подпитки первого контура концентрированным раствором борной кислоты до достижения безопасной стояночной концентрации
- •Литература
6.4. Геометрический параметр цилиндрического реактора без отражателя и поле тепловых нейтронов в нём
Большинство энергетических тепловых реакторов имеют цилиндрическую форму активной зоны или очень близкую к ней. Среди многих соображений при выборе формы активной зоны побеждает стремление сделать её симметричной, технологичной и удобной для организации теплосъёма.
Геометрический параметр Вг2 критической цилиндрической активной зоны может быть с равным успехом найден и из решения уравнения критичности, и из решения волнового уравнения, но первое является трансцендентным и не разрешается относительно В2 аналитически, следовательно, при известных величинах k, L2 и т из уравнения критичности можно найти только величину В2 (методом последовательных приближений), но нельзя получить удобной аналитической зависимости Вг2 от размеров активной зоны (радиуса и высоты её).
Такое выражение В2 = f(Rаз,Hаз) можно получить только после решения волнового уравнения и разрешения его относительно Вг2:
Bг2 = - [ 2Ф(r)] / [Ф(r)], (6.4.1)
Но волновое уравнение - дифференциальное уравнение второго порядка в частных производных, поэтому для получения конкретного (частного) его решения надо обязательно располагать парой граничных условий.
6.4.1. Граничные условия для решения волнового уравнения. Предположим вначале для простоты, что речь идёт о гомогенной цилиндрической активной зоне, окруженной пустотой (вакуумом). Почему именно пустотой?
Дело в том, что вакуум, кроме того, что он даёт возможность сравнивать различные критические активные зоны в одинаковых условиях, является в некотором смысле абсолютным поглотителем утекающих из активной зоны нейтронов, поскольку он не содержит в себе объектов, с которыми нейтрон может столкнуться, изменить направление движения и вернуться обратно в активную зону.
Е динственной точкой цилиндрической гомогенной активной зоны, о величине плотности потока тепловых нейтронов мы можем хоть что-то сказать, является центр её (середина её высоты по оси симметрии). И единственное, что мы можем сказать относительно плотности потока тепловых нейтронов в этой точке, - то, что величина Ф(r) в ней максимальна, поскольку это - наиболее удалённая от всех границ активной зоны точка, и возможности для утечки тепловых нейтронов из неё за пределы активной зоны минимальны. z
0 r
Рис.6.6. Размещение начала цилиндрической системы координат в
геометрическом центре цилиндрической активной зоны.
И если поместить начало цилиндрической системы координат в центр активной зоны (рис.6.6), то первое граничное условие
Ф(r=0,z=0) = Фо = Фmax
- выглядит неопределённо, так как неясна конкретная величина этого максимума. Это же граничное условие (как условие максимума функции Ф(r,z)) можно записать более определённо:
(dФ/dr)z=o = 0 и (dФ/dz)r=o = 0. (6.4.2)
Второе граничное условие в такой ситуации должно быть обязательно нетривиальным, то есть должно указывать на любое конкретное значение функции Ф(r,z) в какой-либо точке активной зоны. Здесь мы выдыхаемся: при всем желании указать такую точку в пределах активной зоны мы не в состоянии. На действительных границах активной зоны (при r = Rаз или z = Hаз/2) величина плотности потока тепловых нейтронов - явно не нулевая.
Поэтому в качестве второго граничного условия вводится искусственное условие, состоящее в следующем.
Предположим, что распределения Ф(r) и Ф(z) в пределах активной зоны от центра к периферии имеет характер нелинейного уменьшения. Но, если вообразить (рис.6.7), что функция Ф(r), переходя через границу активной зоны, продолжает уменьшаться линейно, причём, с тем же угловым коэффициентом, что и на границе активной зоны, то на некотором отстоянии d от границы активной зоны линейно-экстраполированная таким образом функция Ф(r) уменьшается до нуля.
Отстояние (d) от границы активной зоны в вакуум, на котором линейно-экстраполированная на границе активной зоны функция распределения плотности потока тепловых нейтронов обращается в нуль, называется длиной линейной экстраполяции.
Условная цилиндрическая поверхность, эквидистантно отстоящая от поверхности реальной цилиндрической активной зоны на длину линейной экстраполяции, называется экстраполированной границей активной зоны.
Величины полуразмеров реальной цилиндрической активной зоны, увеличенные на длину линейной экстраполяции, называют экстраполированными полуразмерами активной зоны:
R' = Rаз + d (6.4.3)
H'/2 = Hаз/2 + d (6.4.4)
Ф ( r )
Действительные границы активной зоны
Касательная к графику Ф(r) на границе активной зоны
r
Длина линейной экстраполяции d
Экстраполированные границы активной зоны
Действительный радиус
активной зоны Rаз Rаз + d = R - экстраполированный радиус активной зоны
Рис.6.7. К понятию длины линейной экстраполяции, экстраполированных границ
и экстраполированных размеров активной зоны.
Кинетическая теория даёт простую формулу для длины линейной экстраполяции:
d = 0.7104 tr = 0.7104 / tr (6.4.5)
Учитывая сказанное, второе граничное условие звучит просто:
На экстраполированных границах активной зоны величина плотности потока тепловых нейтронов равна нулю: Ф(r=R') = 0 и Ф(z=H'/2) = 0 (6.4.6)
6.4.2. Результат решения волнового уравнения для цилиндрической гомогенной активной зоны. Если записать волновое уравнение в цилиндрической системе координат, начало которой совпадает с центром активной зоны, и решить его при обозначенных выше граничных условиях, то интеграл этого уравнения будет иметь вид:
(6.4.7)
Выражение (6.4.7) означает, что:
- распределение величины плотности потока тепловых нейтронов по высоте цилиндрической гомогенной активной зоны (в точках равноудаленных от оси симметрии на расстояние r) подчиняется закону косинуса:
Ф(z) r=idem = Фоr cos(z/H'), (6.4.8)
где Фоr = Ф(z=0, r) - значение плотности потока тепловых нейтронов на цилиндрической поверхности радиуса r на середине высоты активной зоны (рис.6.8):
d
d d
r
d
Рис.6.8. Эпюры распределения плотности потока тепловых нейтронов по высоте цилиндрической
гомогенной активной зоны по оси симметрии и на разных отстояниях от оси.
z
d
d d
r
d
Рис.6.9. Эпюры распределения плотности потока тепловых нейтронов по радиусу
цилиндрической гомогенной активной зоны на разных уровнях по её высоте.
распределение плотности потока тепловых нейтронов по радиусу активной зоны (в плоских круговых поверхностях на любой фиксированной высоте z над (или под) центром активной зоны) подчиняется закону функции Бесселя первого рода нулевого порядка:
Ф(r) z=idem = Фоz Io(2.405r/R'), (6.4.9)
где Фоz = Ф(z,r=0) - значение плотности потока тепловых нейтронов на оси симметрии активной зоны на высоте z (рис.6.9).
Функция Бесселя первого рода нулевого порядка Io(x) для действительного аргумента x появляется при решении волнового уравнения в цилиндрической системе координат. Начальный участок графика этой функции (при изменении x в пределах от 0 до 2.405) напоминает график функции косинуса в пределах от 0 до /2: при x = 0 Io = 1, а при x = 2.405 Io = 0 (рис.6.10). Более того, значения этих функций при значениях аргумента x в указанных интервалах их с точностью до + 2% совпадают.
I0(x) 1.0
0.5
0
- 0.5
0 1 2 3 4 5 х
Рис.6.10. График функции Бесселя первого рода нулевого порядка Io(x) для действительного аргумента.
В связи с тем, что график Io(x) пересекает ось абсцисс при xo = 2.405, это значение аргумента называют первым корнем (или первым нулём) функции Бесселя первого рода нулевого порядка.
Функция Io(x), наряду с другими бесселевыми функциями, приводится в графическом и табличном виде в справочниках по специальным функциям.
Характер косинусоидально-бесселевского распределения плотности потока тепловых нейтронов в цилиндрической гомогенной активной зоне действителен (совпадает с реальным) для любых точек активной зоны, исключая точки, лежащие в пределах относительно тонкого приграничного слоя толщиной ~ 2tr среды активной зоны, где действительный характер распределения Ф(z,r) несколько отклоняется от аналитического в сторону увеличения.
Учитывая, что транспортные макросечения сред активных зон ВВЭР не превышают нескольких см -1, соответствующие им величины длины линейной экстраполяции d оказываются не выше 1 см. Поэтому распределение Ф(z,r) в цилиндрических гомогенных активных зонах с размерами более 1 м фактически определяется не столько величиной d, сколько действительными размерами активной зоны.
Этот вывод справедлив и для гетерогенных тепловых реакторов.
6.4.3. Выражение для геометрического параметра цилиндрической активной зоны. Это выражение получается путём подстановки решения волнового уравнения (6.4.7) в равенство (6.4.1). После преобразований получается:
Bг2 = (/H')2 + (2.405/R')2 (6.4.10)
Как видим, геометрический параметр имеет размерность см-2, а его величина обратно пропорциональна квадрату линейных размеров активной зоны реактора.
О величине геометрического параметра говорят такие цифры:
- для реактора космической спутниковой электростанции (R' 6 см, H'~ 11 см) величина Вг2 0.2422 см-2;
- для реактора морского атомохода (R' 50 см, H' 100 см) Вг2 3.3 10-3 см-2;
- для реактора ВВЭР-1000 (R' = 158 см, H' = 355 см) Вг2 = 3.1 10-4 см-2;
- для реактора РБМК-1000 (R' = 590 см, H' = 700 cм) Вг2 = 3.7 10-5 см-2.
Падающий характер изменения величины Вг2 с ростом линейных размеров активной зоны позволяет качественно разрешить вопрос о соотношении величин геометрического и материального параметров в некритических реакторах (в критических реакторах, как уже отмечалось, Вг2 = Вм2).
Величина материального параметра для любого реактора определяется только составом материалов, входящих в его активную зону. Следовательно, для гетерогенного реактора, активная зона которого состоит из одинаковых ячеек, величина материального параметра для всей активной зоны уже определена составом материалов одиночной ячейки: ведь соотношение материалов в одиночной ячейке и во всей активной зоне, составляемой из определённого числа таких ячеек, одинаково. Значит, величина материального параметра от числа размещаемых в его активной зоне ячеек не зависит и в процессе загрузки топливных ячеек в активную зону не меняется.
Теперь представим себе процесс зарядки активной зоны и доведения её до критического состояния: в загруженный замедлителем реактор вначале вставляется центральная ТВС, затем вокруг неё размещается первый слой из 6 таких же ТВС, затем последовательно ставятся на свои места 12 ТВС второго слоя, затем - 18 ТВС третьего слоя и т.д., - до тех пор, пока не будет набрано критическое количество ТВС, при котором в активной зоне начинается самоподдерживающаяся цепная реакция деления.
Ясно, что в процессе доведения активной зоны до критического состояния растёт радиус набора активной зоны, а значит, величина геометрического параметра:
Bг2 = (/H')2 + (2.405/R')2
в процессе набора критической массы будет уменьшаться. И когда активная зона достигнет критичности, величина геометрического параметра снизится до величины материального параметра.
Таким образом, в подкритическом реакторе величина геометрического параметра больше величины материального параметра, а в надкритическом (который получился бы, если бы в активную зону добавили еще одну ТВС сверх критического их количества) - наоборот - величина материального параметра стала бы больше величины геометрического параметра.
6.4.4. Оптимальное соотношение размеров цилиндрической активной зоны. Из всего сказанного в принципе должно быть ясно, как рассчитывать критические размеры активной зоны цилиндрического гомогенного реактора по заданному составу материалов его активной зоны:
а) по составу материалов активной зоны рассчитать величины их эффективных микросечений и средних макросечений для всей среды активной зоны;
б) рассчитать , , , , т и L2, то есть получить k, т и L2;
в) методом последовательных приближений решить уравнение критичности реактора
k exp(-B2т)/(1 + B2L2) = 1
относительно величины В2, являющейся в критическом реакторе и материальным, и геометрическим параметром;
г) подставляя найденную величину В2 в её выражение:
(/H')2 + (2.405/R')2 = B2, (6.4.11)
можно было бы искать экстраполированные критические размеры активной зоны реактора (Н' и R'), но одно уравнение с двумя неизвестными имеет бесчисленное множество пар решений. Иными словами, одному и тому же значению В2 удовлетворяют и блинообразные активные зоны (с малым отношением Н'/R'), и, наоборот, колоннообразные активные зоны (с большим отношением Н'/R'). Следовательно, для получения определённого решения уравнения (6.4.11) необходимо задаться величиной соотношения размеров активной зоны (Н'/R'). Из каких соображений?
- Из соображений экономии нейтронов: из стремления при данной величине объёма активной зоны сделать минимальной утечку тепловых нейтронов. При одинаковой плотности тока утечки тепловых нейтронов по всей поверхности активной зоны решение задачи на минимум утечки сводится к решению задачи на минимум поверхности цилиндрической активной зоны при заданном ее объёме. Это имеет место при соотношении (Н'/R') = 2, то есть когда высота цилиндра равна его диаметру.
Но на цилиндрической части поверхности активной зоны градиент плотности потока тепловых нейтронов получается немного выше, чем на плоских поверхностях верхнего и нижнего торцов активной зоны, а, значит, величины плотности тока утечки тепловых нейтронов на цилиндрической поверхности будут выше, чем на плоских торцах.
Поэтому для нахождения минимально-возможной общей утечки тепловых нейтронов из активной зоны необходимо решать задачу на экстремум для величины общего тока утечки тепловых нейтронов через всю поверхность активной зоны (S):
Iобщ = I(S) dS
(S)
Решение этой задачи дает оптимальное соотношение размеров цилиндрической активной зоны
(H'/R')opt = 1.948 (6.4.12)
по соображениям экономии тепловых нейтронов в активной зоне.
Цилиндрические активные зоны с (Н'/R')<1.948 принято называть уплощёнными (т.е. более плоскими по сравнению с активными зонами с оптимальным соотношением размеров), а зоны с (H'/R')>1.948 - удлинёнными.
Например, активная зона РБМК-1000 (Наз = 7м, Dаз = 11.8 м) характеризуется отношением Н'/R' 1.19, т. е. является сильно уплощённой, а активная зона ВВЭР-1000 (Наз = 3.55 м, Rаз =1.58 м, Н'/R' 2.25) - является явно сильно удлинённой. В той и другой активных зонах экономия тепловых нейтронов оказалась принесённой в жертву иным соображениям.
В ВВЭР-1000 уменьшение отношения Н'/R' привело бы к увеличению диаметра активной зоны за счёт сокращения её высоты, а вместе с этим - и к увеличению диаметра корпуса реактора, а, значит, - к увеличению толщины стенки корпуса (корпус - сосуд, работающий под большим давлением), материалоёмкости реактора и к увеличению его стоимости. Именно поэтому (главным образом) активная зона ВВЭР-1000 выполнена удлинённой.
У РБМК-1000 (канального реактора) таких проблем нет: активная зона находится под незначительным давлением азотно-гелиевой смеси, охлаждающей графитовую кладку; высокое давление имеет место только внутри труб технологических каналов; уменьшение высоты активной зоны (или высоты технологических каналов) за счёт увеличения диаметра активной зоны оказывается даже благотворным делом: с точки зрения укорочения технологических каналов и увеличения численности параллельно работающих каналов, при котором снижается гидравлическое сопротивление активной зоны, а, значит, - и энергетические затраты на циркуляцию теплоносителя в контуре МПЦ.