Учебная работа. Расчет гидравлических характеристик реакторного контура

1 Звезда2 Звезды3 Звезды4 Звезды5 Звезд (Пока оценок нет)
Загрузка...

Расчет гидравлических характеристик реакторного контура

Введение

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

В данной курсовой работе нам нужно вычислить мощность главного циркуляционного насоса, необходимую для перекачки теплоносителя по петле первого контура ядерной энергетической установки с реактором ВВЭР-1000. Это главного циркуляционного насоса.

Исходные данные для расчета

тепловая мощность реактора Q=2800 [МВт]

Давление в нижнем коллекторе реактора Pн.к.=157 [бар]

температура на входе в активную зону tвх. = 290[°С].

подогрев на активной зоне Dt = 30[°C].

Исходные данные к курсовой работе:

Реактор:

Размер кассеты «под ключ», H238 ммДиаметр твэлов, dтв.9,1 ммКластерные каналы12,6 x 0,85 ммЦентральная трубка10,3 x 0,65 ммКластерные стержни8,2 ммШаг решетки твэл, S12,75 ммТолщина стенки шестигранного чехла1,5 мм

Длина трубопроводов первого контура:

От реактора до парогенератора11.5мОт ПГ до ГЦН13мОт ГЦН до реактора25,3мВнутренний диаметр трубопровода0,87м

Коэффициенты местных сопротивлений

Активная зона:

Вход в демпферные трубы, ξдт1,47Вход в хвостовик кассеты, ξвх.х1,0Концевые решетки, ξк.р.1,15Дистанционирующая решетка, ξд.р.0,96На сужение в верхней части кассеты, ξсуж2Выход из кассеты, ξвых.к2,18Контур:

ГЗЗ на входе в реактор, ξвх.гзз0,7ГЗЗ на выходе из реактор, ξвых.гзз0,76Выход из реактора, ξвых.р0,8Вход в реактор, ξвх.р0,6изгиб трубопровода на 90˚, ξ900,3Потери в ГЦН, ξгцн2,5

Парогенератор:

Площадь поверхности теплообмена, F6486 м2Количество трубок ПГ, n14000 шт.диаметр трубок ПГ16 x 1,5 мм

Местные сопротивления:

Вход в коллектор ПГ, ξвх.к.пг1,0Вход теплоносителя из коллектора в трубки ПГ, ξвх.тр.пг0,5Выход из коллектора, ξвых.к.пг0,9Выход из трубок в коллектор, ξвых.тр.пг1,0

Определение расхода через реактор и через все ТВС

Определим Сp через среднюю температуру и давление в нижнем коллекторе из таблицы Ривкина.

tср

Pн.к.=157 бар

Сp=5648

Расход по одной ТВС

Расчёт геометрии ТВС.

участок

Длина участка

Гидравлический диаметр участка

Площадь

участок

Этот участок является диффузором, на котором происходит переход от круглого к шестиугольному сечению.

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

,

— эквивалентный диаметр, — площадь шестиугольника.

На этом участке происходит расширение, найдем угол .

, тогда .

По таблице найдем поправочный коэффициент на плавное расширение:

Рассчитаем местный коэффициент сопротивления на плавное расширение:

участок

Длина ,

Площадь поперечного сечения ,

Гидравлический диаметр

Площадь проходного сечения =

=

где r-радиус вписанной окружности

Смоченный периметр

Где а — сторона шестиугольника

участок

Длина .

Гидравлический диаметр равен гидравлическому диаметру ячейки ТВС:

где S — шаг решетки твэла

участок

Длина ,

гидравлический диаметр- ,

участок

гидравлический диаметр— ,

площадь поперечного сечения

.

Расчет площадей на 4 и 5 участках для различных типов ТВС:

Кластерные открытые: N1 = 80[шт.]

Площадь поперечного сечения на участках 4 и 5 одинакова и равна:

Кластерные стержни опущены на 40%: N2 = 29[шт.]

Площадь поперечного сечения на участке 4: , площадь поперечного сечения на участке 5:

С выгорающим поглотителем (СВП): N3 = 54[шт]

Для СВП кассет площадь поперечных сечений будет одинаковой, также она останется неизменной на участках 4 и 5 и будет равна

Определение количества дистанционирующих решеток:

Шаг решеток равен следовательно, количество решеток равно

Вывод зависимостей для определения расхода

Для определения расходов воспользуемся соотношениями:

заменим скорость ;

;

Вынесем за скобку и оставшуюся часть обозначим через

Получим

Где

Для нашего случая запишем

где С1, С2, С3 , — постоянные зависящие от коэффициентов сопротивления трения и местных коэффициентов сопротивления, а также геометрии ТВС.

Из этой системы

, ,

.

Потери давления на различных ТВС ; заменим среднюю скорость ; получим

Перенесем G2 в противоположную сторону, ; обозначим

Тогда можно записать

где — коэффициент местного сопротивления.

Возьмем необходимые для выполнения работы значения удельных объемов и динамической вязкости, из таблиц Ривкина.

t0CV,м3/кгΡ,кг/м3m2900,001339746,825992.724*10-63050,0013966716,0246386.702*10-6

Найдем местные потери для разных типов кассет.

  • Кластерные кассеты, стержни подняты
  • Кластерные кассеты, стержни опущены
  • 3) Кассеты СВП

    Запишем основные соотношения.

    Числа Рейнольдса находятся из формулы

    где — скорость, -гидравлический диаметр.

    Заменим ;

    Получим ; обозначим

    Тогда можно записать

    Перепад давления из уравнения ; заменим среднюю скорость ; получим ;

    Перенесем в противоположную сторону, ; обозначим

    Можем записать; обозначим

    Тогда можем записать

    Коэффициент сопротивления трению , так, как в нашем случае числа Рейнольдса будут

    В активной зоне коэффициент сопротивления трению будет больше и определяется по формуле

    где , при

    — относительный шаг решетки твэл

    ;

    — шаг решетки твэл.

    Вычислим для каждого участка коэффициенты и .

    Для кластерный кассет

    1 участок

    3 участок

    4 участок

    5 участок

    Опущенные на 40%:

    6 участок

    Для ТВС СВП

    Определение расходов по различным типам кассет.

    Первоначально числа Рейнольдса будут определятся из условия одинаковости расходов по всем типам кассет.

    Из них значения и .

    Первое приближение

    Кластерные кассеты, стержни подняты

    1 участок

    4 участок

    3 участок

    5 участок

    6 участок

    Кластерные кассеты, стержни опущены

    1 участок

    4 участок

    3 участок

    5 участок

    6 участок

    Кассеты СВП

    1 участок

    4 участок

    3 участок

    5 участок

    6 участок

    Найдем общие потери по длине.

  • Кластерные кассеты, стержни подняты
  • Кластерные кассеты, стержни опущены
  • Кассеты СВП
  • Определим окончательное приближения, и определим расходы.

  • Кластерные кассеты, стержни подняты
  • Кластерные кассеты, стержни опущены
  • Кассеты СВП
  • Для определения расходов воспользуемся соотношениями:

    заменим скорость ;

    ;

    Вынесем за скобку и оставшуюся часть обозначим через

    Получим

    Где

    Для нашего случая запишем

    , ,

    .

    ; ; ;

    Вычислим значения расходов по полученным формулам.

    . Проведем второе приближение, используя полученное

    значение расходов для разных типов кассет.

    1)Кластерные кассеты, стержни подняты

    1 участок

    4 участок

    3 участок

    5 участок

    6 участок

    2)Кластерные кассеты, стержни опущены

    1 участок

    4 участок

    3 участок

    5 участок

    6 участок

    3). Кассеты СВП

    1 участок

    4 участок

    3 участок

    5 участок

    6 участок

    Найдем общие потери по длине.

    ). Кластерные кассеты, стержни подняты

    ). Кластерные кассеты, стержни опущены

    ). Кассеты СВП

    Определим окончательное ) Кластерные кассеты, стержни подняты

    ) Кластерные кассеты, стержни опущены

    ) Кассеты СВП

    Определим расходы.

    Проверим, выполняется ли баланс расходов

    -выполняется.

    Вычислим расхождения между расходами в предыдущей и последующей итерациях выраженное в процентах:

    ; ;

    ; .

    Так как значения расходов во втором и первом приближении отличаются меньше чем на 0,5%, то можно продолжать вычисления с полученными во втором приближении значениями.

    Потери давления в активной зоне:

    Определение потерь давления на петле

    горячая ветка:

    Давление и температура теплоносителя на горячей ветке:

    t=313[C°]

    Р= 15,291 МПа

    Для этих параметров определим значения плотности и динамической вязкости, используя таблицу Ривкина

    ;

    ; .

    Запишем уравнение Бернулли для сечений 4 и 3:

    , учитывая, что V3 и V4 равны, получим: .

    Расход на горячей ветке

    .

    Число Рейнольдса на горячей ветке:

    где .

    Коэффициент сопротивления трения рассчитаем по формуле Никурадзе:

    ,

    Парогенератор:

    давление и средняя температура в парогенераторе:

    t=2990 C

    P=15,150 МПа

    Для этих параметров определим значения плотности и динамической вязкости, используя таблицу Ривкина:

    ; ;

    .

    Расход через парогенератор

    Расход через одну трубку парогенератора

    Число Рейнольдса в трубке парогенератора

    где ,

    Коэффициент сопротивления трения рассчитаем по формуле Никурадзе:

    Длина трубки парогенератора:

    .

    .

    холодная ветка:

    Давление и средняя температура на холодной ветке:

    t=2850 C

    P=15,109 МПа

    Для этих параметров определим значения плотности и динамической вязкости, используя таблицу Ривкина:

    ; ;

    .

    Запишем уравнение Бернулли для сечений 1 и 2:

    , учитывая что V1 и V2 равны, получим: .

    Расход на холодной ветке

    .

    Число Рейнольдса на холодной ветке

    Коэффициент сопротивления трения рассчитаем по формуле Никурадзе:

    .

    Суммарные потери давления:

    Расчет мощности ГЦН:

    .

    циркуляционный насос реактор теплоноситель

    список литературы

    1.Авдеев Е.Ф., «Расчет гидравлических характеристик реакторного контура».

    2.Ривкин С.Л., Александров А.А. «Таблицы теплофизических свойств воды и водяного пара».

    Учебная работа. Расчет гидравлических характеристик реакторного контура