I. Расчет физических характеристик активной зоны
Введение Вариант 6
В качестве исходных данных заданы: - Замедлитель; - Теплоноситель; - Ядерное горючие. В процессе выполнения работы необходимо: 1) Оценить размеры активной зоны; 2) Рассчитать физические характеристики активной зоны реактора (коэффициент размножения нейтронов в бесконечной среде, возраст теплового нейтрона, длина диффузии) в свежем топливе; 3) Расчет изменения изотопного состава горючего для определения кампании топлива, кампании реактора (количество перегрузок); 4) Определение коэффициента размножения нейтронов в процессе выгорания и запаса реактивности реактора; 5) Определение температурного эффекта и коэффициента реактивности; 6) Компенсируемость реактора (число регулирующих стержней для компенсации рассчитанного запаса реактивности). Исходные данные. Тепловая мощность Qр, МВт………………………………………………………………….…5676 Электрическая мощность Nэл, МВт………………………………………………………….….2100 Шаг ячейки а, см……………………………………………………………………………………25 Высота активной зоны НАЗ, м………………………………………………………………………7 Вид топлива………………………………………………………………………………………UO2 Обогащение по U235 х235, %............................................................................................................2,6 Средняя плотность топлива γUO2, г/см3……………………………………...…………………...9,5 Замедлитель…………………………………………………………….……………………...графит Средняя плотность графита γс, г/см3……………………………………….……………………1,65 Средняя температура графита Тс, 0К…………………………………………………………….750 Давление теплоносителя на входе в активную зону p1р, МПа…….………………….…………7 Удельное энерговыделение в расчете на единицу объема активной зоны qv, МВт/м3…………………………………………………………………………….……..4,2 ТВЭЛы……………………………………………………стержневые с наружным охлаждением Наружный диаметр оболочки твэла d2, мм………………………..…………………………..13,5 Внутренний диаметр оболочки твэла d1, мм……………………..…………………………...11,7 Диаметр топливного сердечника dUO2, мм…………………………..………………...……....11,5 Наружный диаметр центральной каркасной трубки dц.тр,2, мм……………………….……….15 Внутренний диаметр центральной каркасной трубки dц.тр,1, мм………………….……………7 Количество твэлов в одной ТВС nтвэл …………………………………………….…………….18 Наружный радиус трубы технологического канала R2, мм………………….………………..44 Внутренний радиус трубы технологического канала R1, мм…………………………………40 Число технологических каналов nтк.........................................................................................3089
Физический расчет I. Расчет физических характеристик активной зоны Поскольку реактор гетерогенный, начинаем расчет с определения объемов Vi в ячейке с шагом а = 25 см, приходящихся на 1см высоты:
Остальную часть объема ячейки занимают зазоры. Содержащийся в них газ можно в расчете не учитывать. Определим ядерные концентрации веществ по формуле: Вычислим теперь макроскопические параметры гомогенизированной активной зоны, нужные для оценки температуры тепловых нейтронов, а именно Σa и ξΣs. Для этого составляем таблицу 1, в которой величины σa и σs взяты при энергии E=0,025 эВ. Таблица 1.
- среднелогарифмическая потеря энергии Замедляющая способность урана пренебрежимо мала и поэтому не учитывается. Для гомогенизированной активной зоны вычислим макроскопические сечения, нужные для расчета температуры нейтронного газа. Пользуясь табл. 1, находим: Полученное значение Σа=0,0163см -1 рассчитано для энергии нейтронов 0,025 эВ, что соответствует максимуму максвелловского спектра при температуре нейтронного газа 200С. Температуру нейтронного газа оценим по формуле: . Температура замедлителя Т0 (средняя по объему активной зоны) определяется расчетом процессов тепловыделения и теплопередачи в кладке реактора: Т0=750 0К. Пересчитаем Σа на эту температуру, предполагая, что справедлив закон : Следовательно, Найдем средние сечения для тепловых нейтронов. Границей тепловой группы Егр условно считается точка пересечения спектров Ферми и Максвелла. Эта точка определяется подбором или графически из следующего трансцендентного уравнения: , где ; Задавшись хгр = 6, по таблицам из справочника по ядерно-физическим константам при Тn = 900 0К находим . Средние сечения всех других элементов (в виду справедливости закона ) можно посчитать по формуле: , где F(xгр)=1,01 (определено по графику при хгр=6). Т. е. , i = [C, Н2О, U238, О, Zr, Er] Пересчитывая данные табл.1, получаем: Σа(Егр,ЕТ) – сечение поглощения среды, усредненное по спектру Максвелла при температуре Тn в интервале энергий от нуля до Егр. Находим: и по графику получаем, что этой величине соответствует хгр=5,9, то есть можно считать, что совпадение заданного и полученного хгр удовлетворительное. Примем в расчет , найденные для хгр=6: ; ; ; ; ; ; Найдем транспортные сечения, усредненные по спектру Максвелла, по формуле: Для этого найдем средний косинус угла рассеяния по формуле: ; ; ; ; ; ; ; ; ; ; . Для воды при расчете транспортного сечения используем особую формулу: В гетерогенной активной зоне неоднородность потока тепловых нейтронов в ячейке может оказаться весьма существенной, что скажется, в конечном счёте, на коэффициенте размножения. Учтём эту неоднородность. Ячейка разбивается на две зоны, центральную - блок, и периферийную - замедлитель. Центральная зона гомогенизируется. Удобно принять в качестве радиуса блока наружный радиус циркониевого канала: Условный радиус ячейки: Таблица 2.
Пользуясь данными табл. 2, находим: Теперь определяем коэффициент и длину диффузии: Вторая зона ячейки – чистый графит. Два штриха означают, что это зона “два”. Будем считать, что источники тепловых нейтронов распределены в каждой зоне ячейки равномерно и мощность их пропорциональна замедляющей способности зон. Мощность источников во второй зоне ячейки можно принять за единицу, тогда в первой зоне она будет равна: В диффузионном приближении формулы для средних нейтронных потоков в первой и второй зоне ячейки имеют вид: ; Не следует придавать какое-либо значение абсолютным величинам и размерности и , так как потоки определяются здесь с точностью до произвольного общего множителя. Для дальнейшего расчета важно только отношение: Теперь приступим к расчету коэффициента размножения: Коэффициент теплового использования θ вычисляем с учетом гетерогенности по формуле: Для вычисления νа находим в справочнике при хгр=6 и Тn=9000К: . Получаем: Коэффициент ε будем рассчитывать с учетом гетерогенности. Пусть блоком является та же самая область ячейки, что и при расчете θ, с радиусом = 4,4см. Для расчета (это сечения, усредненные в области энергий Е>1,4Мэв) составляем табл.3, пренебрегая U-235.
Таблица 3.
Вычисляем по данным табл.3: На рис.4 (см. стр.33 Румянцев Г.Я.) находим вероятность первого столкновения в цилиндрическом блоке: Учитывая, что шаг ячеек реактора довольно велик, принимаем:
Определяем ε: Теперь рассчитаем коэффициент j - вероятность избежать резонансный захват. Для учета эффекта Доплера нужна средняя температура урана, которая должна быть найдена при расчете теплопередачи в твэлах (в нашем случае средняя температура урана является исходным данным ТU = 10000K). Блоком теперь будет называться каждый отдельный твэл без оболочки. Тогда: Подставим все эти величины в выражение: Отсюда: По формуле четырех сомножителей находим: Определим теперь усредненные диффузионные параметры для тепловой и надтепловой группы: D2, L2, D1 и τ.
Для тепловых нейтронов: Σtr вычисляется методом простой гомогенизации: Следовательно: ; Для надтепловой группы, константы должны определяться усреднением по спектру в интервале от Eгр до 14 МэВ. В двухгрупповой теории коэффициент диффузии обычно определяют при E=1 эВ, поскольку для большинства веществ сечение в области E>1 эВ относительно слабо меняется. Для определения Σtr(1эв) составим табл.4. При этом пренебрежем поглощением всех элементов, и для простоты будем считать, что U-235 не отличается от U-238, учитывая, что вклад U-235 в величину Σtr весьма невелик. Воду можно представить как смесь ядер водорода и кислорода. Таблица 4.
В результате получим для гомогенизированной активной зоны: Квадрат длины замедления вычисляем по формуле: Для этого сначала найдем . Коэффициенты Аij берем из книги Галанина А.Д. «Теория ядерных реакторов на тепловых нейтронах». Расчет удобно вести с помощью табл.5. Таблица 5.
По таблице находим:
Это величина соответствует нижней границе замедления Емин = 0,2 эв. В нашем случае: следовательно, нужно убавить на величину: Итак, для активной зоны реактора:
Популярное: Почему стероиды повышают давление?: Основных причин три... Почему двоичная система счисления так распространена?: Каждая цифра должна быть как-то представлена на физическом носителе... Модели организации как закрытой, открытой, частично открытой системы: Закрытая система имеет жесткие фиксированные границы, ее действия относительно независимы... Почему человек чувствует себя несчастным?: Для начала определим, что такое несчастье. Несчастьем мы будем считать психологическое состояние... ©2015-2024 megaobuchalka.ru Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. (2602)
|
Почему 1285321 студент выбрали МегаОбучалку... Система поиска информации Мобильная версия сайта Удобная навигация Нет шокирующей рекламы |