Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Предоплата всего
Подписываем
МЭИ (ТУ)
Кафедра парогенераторостроения
Типовой расчёт по курсу:
Генераторы тепловой энергии
Тепловой расчёт ВВЭР
Студент: Иванов А.А.
Группа: С-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, а температура теплоносителя максимальна в верхней части ТВЭЛа. Максимальный градиент температуры теплоносителя в середине высоты ТВЭЛа.