Реферат: Тепловой расчет реактора
МЭИ (ТУ)
Кафедра парогенераторостроения
Типовой расчёт по курсу:
Генераторы тепловой энергии
Тепловой расчёт ВВЭР
Студент: Иванов А.А.
Группа: С-2-95
Преподаватель: Двойнишников В.А.
Москва 2000 год
Аннотация.
В данной работе решались следующие задачи:
расчёт реактора при m = 1 и qv = 100 и определение его экономичности и надёжности при учёте наложенных ограничений: 1.6 < n < 2.2,
2 < Wт <10 м/с, tоб < 350 оС, tc < 2300 оС.
нахождение области допустимых значений относительной высоты активной зоны m и удельного энерговыделения qv (m = 0.8 … 1.6,
qv = 50 … 150) при учёте наложенных ограничений: 1.6 < n < 2.2,
2 < Wт <10 м/с, tоб < 350 оС, tc < 2300 оС.
для выбранного варианта расчёт температуры сердечника, оболочки и теплоносителя по высоте активной зоны.
Содержание:
Введение Исходные данные Тепловой расчёт реактора при m = 1 и qv = 100 МВт/м33.1. Определение размеров активной зоны реактора и скорости теплоносителя
3.2. Определение коэффициента запаса по критической тепловой нагрузке
3.3. Расчёт максимальных температур оболочки ТВЭЛа и материала
топливного сердечника
3.4. Определение области допустимых значений m и qv
3.5. Расчёт распределения температуры теплоносителя, оболочки и топливного
сердечника по высоте активной зоны реактора
4. Выводы
1. Введение
Назначение и виды тепловых расчётов реакторов.
Тепловой расчет ядерного реактора является одной из необходимых составных частей процесса обоснования и разработки конструкции. Без него невозможны ни предварительные поисковые проработки, ни определение оптимальных проектных решений.
Тепловые расчеты обычно выполняются одновременно с гидравлическим и нейтронно-физическим расчетами реактора. В зависимости от задач, решаемых на том или ином этапе проработки конструкции, различают поисковые и поверочные расчеты
Поисковые тепловые расчеты проводятся в период определения основных конструктивных решений. При их выполнении, как правило, известны тепловая мощность реактора, распределение плотности энерговыделения, вид теплоносителя и его параметры все эти данные получают в результате нейтронно-физического расчета, а также тип и конструкция ТВЭЛов и кассет, определяемых техническим заданием на основе накопленного опыта проектирования, изготовления и эксплуатации. В результате определяются размеры активной зоны и других элементов реактора, находятся, а при необходимости уточняются параметры теплоносителя, определяются характерные температуры, выбираются конструкционные материалы и топливные композиции.
По мере разработки конструкции тепловые расчеты выполняются снова, но более детально, с учетом выбранных конструктивных решений, как для номинального режима, так и для работы на частичных нагрузках. Также обсчитываются тепловые режимы работы оборудования при переходных процессах при пуске, останове, изменении нагрузки, характерных как для штатных ситуаций, так и в аварийных случаях. Во всех этих случаях тепловой расчет носит характер поверочного, и его основной задачей является определение термодинамических характеристик теплоносителя и тепловых параметров характеризующих условия функционирования элементов ядерного реактора. Обеспечение надежной работы реактора в целом и его отдельных элементов, достижение высокой экономичности реакторной установки требует высокой точности определения теплотехнических параметров, что ведет к существенному усложнению всех видов расчетов, в том числе и теплового. Необходимость же их автоматизации приводит к созданию сложных программных комплексов, объединяющих тепловые, Гидравлические, нейтронно-физические и прочностные расчеты.
Настоящий метод ориентирован на использование несколько упрощенного теплового расчета, базирующегося на одномерном представлении протекания процессов тепло - и массообмена в одной ячейке активной зоны реактора.
2. Исходные данные.
Для выполнения теплового расчета водо-водяного энергетического реактора (ВВЭР) в соответствии с упрощенной методикой требуются исходные данные, условно подразделяемые на режимные и конструктивные,
Данные режимного типа:
Тепловая мощность ВВЭР N = 1664.87 МВт
Конструктивные данные:
Характеристики кассеты:Число ТВЭЛов в кассете nТВЭЛ = 331
Шаг решётки а¢ ¢ = 12.75·10-3 м
Размер кассеты “под ключ” а¢ = 0.238 м
Толщина оболочки кассеты δ = 1.5·10-3 м
Характеристика ТВЭЛа:Радиус топливного сердечника r1 = 3.8·10-3 м
Внутренний радиус оболочки r2 = 3.9·10-3 м
Внешний радиус оболочки rq = 4.55·10-3 м
Размер ячейки а = 0.242 м Материал оболочки ТВЭЛов и кассет: 99% циркония и 1% ниобия Топливная композиция: двуокись урана
3.Тепловой расчёт реактора при qv= 100 МВт/м3 и m= 1
Определение размеров активной зоны реактора и скорости теплоносителя. Температура теплоносителя на выходе из реактораtвых = 314 ° C
Принимаем из расчёта парогенератора
Температура теплоносителя на входе в реакторtвх = 283 ° C
Принимаем из расчёта парогенератора
Перепад температур теплоносителя между входом и выходомΔtт = tвых - tвх = 314 – 283 = 31 ° С
Температура воды на линии насыщенияЗапас до температуры кипения δt = 30 ° C
ts = tвых + δt = 314 + 30 = 344 ° C
Давление в реактореP = 15.2 МПа
Расход воды (теплоносителя) на один реакторсредняя температура воды в реакторе tср = = 298.5 ° C
средняя теплоёмкость воды Cp = 5.433 кДж/кг
Gт = =9885.05 кг/с
Принимаем из расчёта парогенератора.
Объём активной зоны реактора.Средняя плотность тепловыделения АЗ реактора qv = 100 МВт/м3
VАЗ = = 16.648 м3
Диаметр активной зоны реактораПараметр m* = = 1
DАЗ = = 2.767 м
Число кассет в активной зонеПлощадь поперечного сечения ячейки: Sяч = 0.866·a2 = 5.072·10-2 м2
= 178.2 шт.
т.к.дробное, то округляем его до ближайшего большего целого числа
Nкас = 179 шт. с последующим уточнением величин:
DАЗ== 3.4 м
m = = 0.993
Высота активной зоны реактораHАЗ = m·DАЗ = 0.993·3.4 = 3.376 м
Тепловыделение в ТВЭЛахДоля теплоты выделяемая в ТВЭЛах κ1 = 0.95
Qт = κ1·N = 0.95·3064 = 2910.8 МВт
Суммарная поверхность ТВЭЛF = 2·π·rq·HАЗ·nТВЭЛ·Nкас = 2·π·4.55·10-3·3.376·331·179 = 5719 м2
Расход теплоносителя через одну кассетуGтк = = 90.22 кг/с
3.1.14. Скорость теплоносителя в активной зоне реактора
сечение для прохода теплоносителя около одного ТВЭЛа SвТВЭЛ = 0.866·(a¢ ¢ )2-
-π·rq2 = 0.866·(12.75·10-3)2 – π·(4.55·10-3)2 = 7.574·10-5 м2
сечение для прохода теплоносителя в кассете Sвкас = SвТВЭЛ·nТВЭЛ = 7.574·10-5·331 = 2.507·10-2 м2
плотность воды при средней температуре и давлении в реакторе ρв = 713.2 кг/м3
Wт = = 5.046 м/с
Определение коэффициента запаса по критической тепловой нагрузке.3.2.1. Коэффициенты неравномерности тепловыделения
Эффективная добавка отражателя δ0 = 0.1 м
Эффективная высота активной зоны Hэф = HАЗ + 2·δ0 = 3.376 + 2·0.1 = 3.576 м
по оси реактора: Kz = = 1.489
по радиусу активной зоны: Kr = = 2.078
3.2.2. Коэффициент неравномерности тепловыделения в объёме АЗ
Kv = Kz·Kr = 1.489·2.078 = 3.094
Максимальная величина тепловой нагрузки на единицу поверхности ТВЭЛаСредняя тепловая нагрузка на единицу поверхности ТВЭЛа qF = = =0.509 МВт/м2
qmax = qF·Kv = 0.509·3.094 = 1.575 МВт/м2
Критический тепловой поток кризиса первого рода для трубы d = 8 ммТеплота парообразования теплоносителя R = 931.2 кДж/кг
Температура воды на линии насыщения ts = 347.32 ° C
Величина паросодержания теплоносителя в центральной точке реактора xкр = = = -0.2782
qкр(8) =
=
= 1.347·3.5990.5549·е0.4173 = 4.161 МВт/м2
Критический тепловой поток кризиса первого рода для труб диаметром 2rqqкр(2rq) = = 3.901 МВт/м2
Коэффициент запаса по критической нагрузке.nзап = = 2.477
Расчёт максимальных температур оболочки ТВЭЛа и материала топливного сердечника.3.3.1. Максимальное тепловыделение в центре реактора приходящееся на единицу высоты ТВЭЛа.
ql,0= = 4.503·10-2 МВт/м
Коэффициент теплоотдачи от стенки к теплоносителю.Коэффициент теплопроводности теплоносителя λ = 548.3·10-3 Вт/(м·К) при температуре tcр
Эквивалентный диаметр сечения для прохода воды dэкв = = 6.851·10-3 м
Кинематическая вязкость воды. Для её определения необходимо найти динамическую вязкость. μ = 8.936·10-5 Па/с. ν = = 1.253·10-7 м2/с
Критерий Рейнольдса Re = = 2.759·105
Число Прандтля Pr = 0.9217
α==3.685·104 Вт/м2К
Перепад температуры между оболочкой ТВЭЛа и теплоносителем в центре реактора.Δθа0 = = 40.61 ° С
Координата в которой температура на наружной поверхности оболочки ТВЭЛа максимальна.Z*==0.4287м
Максимальная температура наружной поверхности оболочки ТВЭЛаt= 351.7 ° C
Температурный перепад в цилиндрической оболочке ТВЭЛаКоэффициент теплопроводности материала оболочки λоб = 24.1 Вт/(м·К)
Δθоб0 = = 43.55 ° С
Температурный перепад в зазоре ТВЭЛаКоэффициент теплопроводности газа в зазоре λз = 30 Вт/(м·К)
Δθз0 = = 18.52 ° С
Температурный перепад в цилиндрическом сердечникеКоэффициент теплопроводности в цилиндрическом сердечнике λс = 2.7 Вт/(м·К)
Δθс0 = = 1261 ° С
Перепад температур между теплоносителем и топливным сердечникомΔθс = Δθа0 + Δθоб0 + Δθз0 + Δθс0 = 42.46 + 43.55 + 18.52 + 1261 = 1366 ° С
Максимальная температура топливного сердечникаt = 1674 ° C
3.4 Определение области допустимых значений m и qv
Исходные данные для расчёта по программе WWERTR
Тепловая мощность реактора [МВт] Давление в реакторе [МПа] Перепад температур воды [° C] Радиус топливного сердечника ТВЭЛа [м] Внутренний радиус оболочки ТВЭЛа [м] Внешний радиус оболочки ТВЭЛа [м] Шаг решетки [м] Размер кассеты “под ключ” [м] Размер ячейки [м] Толщина оболочки кассеты [м] Эффективная добавка отражателя [м] Число ТВЭЛов в кассете [шт] Температура воды на линии насыщения [° С] Теплота парообразования [кДж/кг] Теплоемкость воды [кДж/кг·К] Теплопроводность воды [Вт/м·° С] Кинематическая вязкость воды [м2/с] Число Прандтля Плотность воды [кг/м3] Теплопроводность оболочки ТВЭЛа [Вт/м·° С] Теплопроводность газа в зазоре ТВЭЛа [Вт/м·° С] Теплопроводность двуокиси урана [Вт/м·° С] Удельное энерговыделение [кВт/л] Относительная высота активной зоны Расч. скорость воды [м/с] Расч. коэффициент запаса Расч. координата точки с мак. темп. оболочки [м] Расч. мак. температура оболочки ТВЭЛа [° С] Расч. мак. температура сердечника ТВЭЛа [° С] |
N = 1664.84 P = 15.2 Δt = 31 r1 = 3.8·10-3 r2 = 3.9·10-3 rq = 4.55·10-3 а¢ ¢ = 12.75·10-3 а¢ = 0.238 а = 0.242 δ = 1.5·10-3 δ0 = 0.1 nТВЭЛ = 331 ts = 344 R = 1020.9 Cp = 5.433 λ = 556.658·10-3 ν = 1.21·10-7 Pr = 0.905 ρв = 724.4 λоб = 23.9 λз = 30.5 λс = 2.7 qv = 100 m = 0.995 Wт = 4.345 nзап = 2.699 Z*= 0.333 t= 343.957 t = 1623.37 |
Результаты расчёта по программе WWERTR.
№ |
m* |
DАЗ |
Wт |
nзап |
Z* |
t |
t |
- |
м |
м/с |
- |
м |
° С |
° С |
|
qv = 50.0 кВт/л |
|||||||
1 2 3 4 5 |
0.800 1.004 1.203 1.409 1.608 |
4.602 4.267 4.018 3.812 3.647 |
2.754 3.204 3.614 4.015 4.386 |
3.433 3.731 3.990 4.234 4.451 |
0.546 0.699 0.850 1.007 1.160 |
345.5 342.1 339.7 337.8 336.3 |
1016.8 1013.4 1010.5 1007.7 1005.2 |
qv = 75.0 кВт/л |
|||||||
1 2 3 4 5 |
0.802 1.006 1.201 1.405 1.611 |
4.018 3.726 3.512 3.333 3.184 |
3.614 4.202 4.730 5.253 5.755 |
2.707 2.941 3.141 3.332 3.510 |
0.413 0.530 0.645 0.766 0.889 |
351.4 347.2 344.2 341.9 340.0 |
1343.5 1339.7 1336.2 1332.7 1329.3 |
qv = 100.0 кВт/л |
|||||||
1 2 3 4 5 |
0.804 1.001 1.209 1.405 1.604 |
3.647 3.390 3.184 3.028 2.897 |
4.386 5.076 5.755 6.362 6.950 |
2.290 2.482 2.662 2.817 2.962 |
0.339 0.433 0.533 0.630 0.729 |
356.2 351.5 347.9 345.3 343.2 |
1662.9 1659.0 1654.7 1650.7 1646.7 |
qv = 125.0 кВт/л |
|||||||
1 2 3 4 5 |
0.801 1.005 1.213 1.411 1.605 |
3.390 3.143 2.953 2.807 2.689 |
5.076 5.905 6.692 7.405 8.067 |
2.009 2.183 2.341 2.479 2.602 |
0.289 0.372 0.459 0.543 0.627 |
360.5 355.1 351.1 348.2 346.0 |
1976.9 1972.6 1967.8 1963.2 1958.7 |
qv = 150.0 кВт/л |
|||||||
1 2 3 4 5 |
0.806 1.010 1.206 1.412 1.609 |
3.184 2.953 2.784 2.641 2.528 |
5.755 6.692 7.529 8.365 9.126 |
1.812 1.969 2.102 2.231 2.345 |
0.256 0.330 0.402 0.479 0.555 |
364.0 358.3 354.2 350.9 348.4 |
2286.2 2281.8 2276.9 2271.4 2266.1 |
m = 0.8
m = 1.0
m = 1.2
m = 1.4
m = 1.6
Границы возможного диапазона значений qv
для каждого параметра (по графикам).
m параметры |
0.8 |
1.0 |
1.2 |
1.4 |
1.6 |
Wт |
— |
— |
— |
— |
— |
— |
— |
— |
— |
— |
|
nзап |
— |
— |
— |
— |
— |
108.1 |
123.6 |
139.9 |
— |
— |
|
t |
68.83 |
91.04 |
116.4 |
141.6 |
— |
— |
— |
— |
— |
— |
|
t |
— |
— |
— |
— |
— |
— |
— |
— |
— |
— |
|
Диапазон допустимых значений |
— |
— |
— |
— |
— |
— |
— |
— |
— |
— |
Прочерк в таблице означает, что максимальное или (и) минимальное значение величины находится за границами рассматриваемой области.
Знак "*" означает, что ни одно значение не входит в накладываемые ограничения.
Анализ таблицы показывает, что при заданных начальных условиях не существует значений m и qv, которые удовлетворяли бы наложенным ограничениям.
3.5. Расчёт распределения температуры теплоносителя, оболочки и топливного сердечника по высоте активной зоны реактора. m = 1.4, qv = 125 кВт/л.
№ |
Координата, м |
Температура теплоносителя, ° С |
Температура сердечника, ° С |
Температура оболочки, ° С |
1 2 3 4 5 6 7 8 9 10 11 |
-1.981 -1.585 -1.188 -0.792 -0.396 0.000 0.396 0.792 1.188 1.585 1.981 |
292.0 293.1 295.5 299.0 303.3 308.0 312.7 317.0 320.5 322.9 324.0 |
416.8 898.6 1328.0 1666.8 1885.0 1963.2 1894.4 1684.9 1353.1 928.5 448.8 |
294.8 306.6 318.5 329.5 338.5 344.9 348.0 347.5 343.6 336.4 326.8 |
Выводы по проведённой работе.
При m = 1 и qv = 100 получено, что данный пример не удовлетворяет условию экономичности n = 2.477 (1.6 < n < 2.2) и незначительно условию надёжности tоб = 351.7 oC (tоб < 350 oC).
При заданных начальных условиях характеристики теплоносителя и реактора, и поставленных ограничениях на скорость теплоносителя, коэффициент запаса, максимальную температуру оболочки и теплоносителя; области допустимых значений относительной высоты активной зоны m и удельного энерговыделения qv (m = 0.8 … 1.6, qv = 50 … 150) не существует. Во всех случаях кроме последнего (m = 1.6 и qv = 150, здесь n > 2.2) не проходит по надёжности.
При расчёте температур по высоте активной зоны получено для m = 1.4 и qv = 125: температура сердечника максимальна в середине высоты ТВЭЛа, температура оболочки максимальна на высоте z = 0.5, а температура теплоносителя максимальна в верхней части ТВЭЛа. Максимальный градиент температуры теплоносителя в середине высоты ТВЭЛа.