- •Севастопольский институт ядерной энергии и промышленности
- •Основы теории ядерных реакторов Курс для эксплуатационного персонала аэс
- •Содержание
- •Перечень сокращений
- •Тема 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. Расчёт времени подпитки первого контура концентрированным раствором борной кислоты до достижения безопасной стояночной концентрации
- •Литература
5.5. Вероятность избежания утечки замедляющихся нейтронов
Вероятность избежания утечки замедляющихся нейтронов pз - это доля нейтронов, избежавших утечки при замедлении, от общего числа нейтронов поколения, начавших процесс замедления в активной зоне. Но величину этой вероятности можно переосмыслить и по-другому:
*) Имеются в виду активные зоны одинакового состава.
Выражение для скорости генерации тепловых нейтронов в реальной активной зоне qт получено в предыдущем пункте. Подходя к величине qт с теми же мерками в рассуждениях, что и к qт, несложно получить:
qт = k aФ (5.5.1)
Подставляя (5.4.14) и (5.5.1) в приведенную логическую формулу для pз, имеем:
pз = exp (- B2т) (5.5.2)
Сравнивая (5.5.2) с начальным предположением (5.1.1), мы должны согласиться, что гипотеза (5.1.1) была не лишена оснований: pз действительно определяется, во-первых, величиной параметра реактора (позже убедимся, что параметр B2 имеет и геометрический смысл), а, во-вторых, - величиной комплексной характеристики замедляющих свойств среды активной зоны реактора, каковой и является возраст тепловых нейтронов (величина, равная шестой части среднего квадрата пространственного смещения замедляющегося нейтрона, то есть пропорционально связанная с величиной квадрата средней длины замедления).
Чем выше величина возраста тепловых нейтронов в реакторе (то есть чем хуже замедляющие свойства среды активной зоны), тем меньше величина вероятности избежания утечки замедляющихся нейтронов, поскольку величина возраста определяет толщину приграничного слоя активной зоны, из которого возможна утечка замедляющихся нейтронов. Чем меньше замедляющихся нейтронов располагают возможностью для утечки, тем выше доля замедляющихся нейтронов, которые останутся к концу замедления в активной зоне (то есть выше величина pз).
Что же касается другой величины, определяющей значение pз, - параметра реактора B2, то пока можно лишь сказать, что с геометрией активной зоны эта величина каким-то образом связана, о чём свидетельствует её размерность - см-2. С непростым смыслом величины B2 ещё предстоит познакомиться подробнее.
Спектр замедляющихся нейтронов Ферми в гомогенной непоглощающей среде
Игнорируя вывод, приведём конечный вид спектра Ферми (с выводом можно познакомиться, например в [8]).
Распределение величины плотности потока замедляющихся нейтронов в непоглощающих средах (имеются в виду свойства не поглощать замедляющиеся нейтроны, а не тепловые) оказывается подчинённым закономерности:
(5.6.1)
Это выражение справедливо как для простой однородной среды, так и для сложных гомогенных сред, состоящих их нескольких сортов замедляющих ядер. В этом случае в формулу (5.6.1) подставляется сумма значений замедляющей способности k компонентов сложной среды
s = 1s1 + 2s2 + 3s3 + ... + ksk = (5.6.2)
Учитывая классическую зависимость кинетической энергии нейтрона от его скорости (Е = mv2/2) и связь плотности потока и плотности нейтронов одинаковой скорости
выражение для спектра замедляющихся нейтронов можно записать так:
(5.6.3)
Таким образом, величина плотности замедляющихся нейтронов по энергиям в непоглощающей среде распределяется по закону "Е-3/2", то есть плавно возрастает с уменьшением энергии нейтронов в процессе их замедления, и столь же плавно (без скачков и изломов) переходит при Е = Ес в максвелловский спектр тепловых нейтронов (рис.5.9).
В реальной замедляющей среде активной зоны, которой свойственно очень слабое поглощение замедляющихся нейтронов в замедлителе, и заметное поглощение их в топливе твэлов (где содержится резонансный захватчик замедляющихся нейтронов – уран-238), реальный спектр замедляющихся нейтронов проходит ниже изображённого на рис.5.9., но качественно выглядит так же.
n(E)
Спектр тепловых нейтронов (Максвелла)
Спектр замедляющихся нейтронов (Ферми)
E
Ec
Рис.5.9. Граница тепловых и замедляющихся нейтронов - энергия сшивки энергетических
спектров Максвелла и Ферми (спектры ненормированные).
Уже отмечалось, что об энергетическом спектре нейтронов в реакторе имеет смысл говорить только для критического реактора, так как любой энергетический спектр в поглощающей среде имеет динамически равновесный характер: плотность нейтронов любой энергии Е поддерживается неизменной во времени в любом единичном объёме за счёт неизменной разницы скоростей:
- прихода замедляющихся нейтронов на уровень энергии Е из области более высоких энергий;
- ухода замедляющихся нейтронов с уровня энергии Е в область более низких энергий;
- появления в единичном объёме новых нейтронов энергии Е за счёт делений ядер (если таковые имеются в рассматриваемом единичном объёме);
- поглощения нейтронов при энергии Е (если в единичном объёме наличествуют поглотители) и
- утечки нейтронов энергии Е из единичного объёма (понимая под утечкой разницу скоростей ухода и прихода нейтронов с энергией Е в этом единичном объёме).
В рамках одногруппового возрастного приближения договорились считать, что поглощение в диапазоне энергий замедления отсутствует, а заметное поглощение замедляющихся нейтронов в реальных средах учитывать с помощью вероятности избежания резонансного захвата (). В этом приближении среда активной зоны поглощает только тепловые нейтроны. И если предполагать, что энергетический спектр тепловых нейтронов - максвелловского типа, величина и положение максимума на шкале энергий в нём явно должны зависеть от поглощающей характеристики среды (a) и замедляющей способности среды (s) в области энергий перехода от замедляющихся нейтронов к тепловым (то есть в области энергии сшивки Ес) и ниже (то есть в пределах самого спектра тепловых нейтронов).
Действительно, поскольку в тепловой области энергий микросечения поглощения изменяются по закону "1/v" (или"Е-1/2"), то основное поглощение тепловых нейтронов происходит при более низких энергиях левого крыла спектра Максвелла; чем выше величина микросечения поглощения активной зоны a (а значит - и a), тем больше тепловых нейтронов поглощается при Е < Eнв, тем больше "выедание" левого крыла спектра, а это значит, что положение максимума спектра (характеризуемое Енв - наиболее вероятной энергией тепловых нейтронов) с увеличением поглощающей способности среды должно смещаться вправо, в область более высоких энергий. Иначе говоря, с увеличением поглощающих свойств среды энергетический спектр тепловых нейтронов ужестчается. Но, т.к. положение максимума в спектре тепловых нейтронов определяет температуру нейтронов Тн (ведь Енв= kTн), то можно выразиться иначе: с увеличением поглощающих свойств среды повышается температура нейтронов в ней. И чем больше величина макросечения поглощения среды (a) - тем больше величина температуры нейтронов в ней (Тн) отклоняется от её термодинамической температуры (Т) в сторону увеличения.
С другой стороны, чем выше величина замедляющей способности среды (s), тем с большей скоростью пополняется за счёт замедления нейтронов весь спектр тепловых нейтронов (в том числе и его левое крыло). Поэтому положение его максимума, наоборот, смещается влево, в область более низких энергий, т.е. с увеличением замедляющей способности среды спектр тепловых нейтронов «умягчается», и температура тепловых нейтронов в среде с лучшими замедляющими свойствами меньше отличается от термодинамической температуры этой среды по сравнению со средой с более слабыми замедляющими свойствами.
Таким образом получается, что температура нейтронов Тн находится в прямой зависимости от величины a среды активной зоны и в обратной зависимости - от s.
Это дало повод к предположению, что обе зависимости являются пропорциональными, что дает лёгкую возможность построить полуэмпирические зависимости для расчёта температуры нейтронов в тепловом реакторе:
Тн = Тз[1 + 1.8 (a/s)] - для уран-водных гомогенных сред, (5.6.4) и
Тн = Тз[1 + 0.91(Aa/s)] - при использовании других замедлителей (с массовым числом A) (5.6.5)
Таким образом, температура нейтронов в тепловом реакторе - величина, прямо пропорциональная термодинамической температуре активной зоны и величине, обратной коэффициенту замедления среды в ней.
Приведённые формулы получены, строго говоря, для гомогенных смесей топлива и замедлителя, однако, с достаточной точностью могут служить и для оценки температуры нейтронов в гетерогенных активных зонах реакторов соответствующих типов, для чего в них должны подставляться средняя термодинамическая температура замедлителя Тз и величины гомогенизированных макросечений активной зоны a и s.
*) Гипотеза насчёт пропорциональности Тн,Тз и обратного коэффициента замедления являлась бы строго-доказательной только в том случае, если бы было строго доказано, что энергетический спектр тепловых нейтронов, как бы ни деформировала его среда своими поглощающими и замедляющими свойствами, всегда сохраняет максвелловскую форму.
Разговор был посвящен фермиевскому спектру замедляющихся нейтронов, но был смещен к спектру тепловых нейтронов, чтобы подчеркнуть два важных момента:
- Во-первых, между двумя этими энергетическими спектрами в критическом тепловом реакторе существует какая-то неясная пока взаимосвязь: любое, самое малейшее, изменение в спектре тепловых нейтронов обязательно должно повлечь за собой изменение и в спектре замедляющихся нейтронов, и, наоборот, изменение в фермиевском спектре должно отразиться и на максвелловском распределении тепловых нейтронов. Эту взаимосвязь можно было бы предсказать и без углубленного экскурса в спектры, если просто (по-материалистически) уверовать в то, что всем поведением замедляющихся и тепловых нейтронов управляет одна Среда, только Среда и исключительно Среда. Это Она, благодаря присущим ей природным свойствам, единым образом, с единой закономерностью рождает с определённой скоростью быстрые нейтроны деления, замедляет их с определённой (связанной со скоростью генерации быстрых нейтронов) интенсивностью, поглощает их с определенной (связанной со скоростями генерации и замедления) скоростью, и определяет скорость утечки нейтронов из любого микрообъёма, - и все это природное предопределение свойственно нейтронам любой энергии. Иначе говоря, Среда, в силу своих природных свойств, и расставляет все нейтроны по энергиям в единый энергетический спектр, не являющийся ни спектром Уатта, ни спектром Ферми и ни спектром Максвелла.
- Во-вторых, хотя спектры Уатта, Ферми и Максвелла и отражают закономерности рождения, замедления и поглощения нейтронов, связать их в единое аналитическое выражение для единого энергетического спектра реактора не удаётся, так как спектры Уатта и Максвелла нормированы соответственно на один нейтрон деления и тепловой нейтрон, а спектр Ферми - вообще не является нормированным.
Не стоит воспринимать сказанное как критику существующей физической теории реакторов: скорее, это её очередная проблема, с которой любому человеку, связанному с ядерной энергетикой, следует быть знакомым. Потому, что её разрешение радикальным образом изменит не только систему теоретических представлений и облегчит их восприятие: должны открыться новые возможности для простых и быстрых инженерных расчётов и анализа эксплуатационной безопасности любых реакторов.
Пока же мы вынуждены довольствоваться упомянутой триадой энергетических спектров и для расчётов критических свойств реакторов пользоваться громоздкими (хотя и приближёнными) вычислениями с помощью ЭВМ.