- •«Томский политехнический университет» а.Г. Коротких, и.В. Шаманин основы гидродинамики и теплообмена в ядерных реакторах
- •Содержание
- •Глава 1. Основные дифференциальные уравнения 7
- •Глава 2. Теплопроводность в ядерном реакторе 26
- •Глава 3. Конвективный теплообмен в ядерном 62
- •Глава 4. Расчет теплофизических 91
- •Введение
- •Глава 1. Основные дифференциальные уравнения тепломассообмена
- •1.1. Вязкость
- •1.2. Расход жидкости
- •1.3. Уравнение неразрывности потока
- •1.4. Статическое давление
- •1.5. Уравнение теплопроводности
- •1.6. Условия однозначности для процессов теплопроводности
- •1.7. Уравнение движения
- •Глава 2. Теплопроводность в ядерном реакторе
- •2.1. Температурное поле
- •2.2. Температурный градиент
- •2.3. Тепловой поток. Закон Фурье
- •2.4. Коэффициент теплопроводности
- •2.6. Теплопроводность при наличии внутренних источников тепла
- •Глава 3. Конвективный теплообмен в ядерном реакторе
- •3.1. Факторы, влияющие на теплоотдачу
- •3.2. Моделирование процессов конвективного теплообмена
- •3.3. Теплоотдача в однофазной среде при свободном движении жидкости (естественная конвекция)
- •3.4. Теплоотдача в однофазной среде при вынужденном течении жидкости
- •3.5. Теплоотдача при вынужденном поперечном омывании труб и пучков труб
- •3.6. Теплоотдача жидких металлов
- •3.7. Теплоотдача при кипении
- •Глава 4. Расчет теплофизических и гидродинамических параметров элементов активной зоны
- •4.1. Расчет удельных значений тепловыделения по длине технологического канала
- •4.2. Расчёт параметров теплоносителя по длине тк и количества теплоты, выделяющейся на отдельных участках и в центральном тк
- •4.3. Расчёт коэффициентов теплоотдачи с поверхности твэл
- •4.4. Расчёт распределения температуры в твэл
- •4.5. Распределение температуры в блоке замедлителя
- •V зам тк внеш
- •4.6. Распределение температуры в органах регулирования
- •4.7. Гидравлический расчёт технологического канала
- •4.8. Оценка критических тепловых потоков
- •4.9. Оценка напряжённого состояния элементов конструкции тк
- •Заключение
- •1. Дементьев б.А. Ядерные энергетические реакторы: Учебник.
- •3. Исаченко в.П., Осипова в.А., Сукомел а.С. Теплопередача. –
- •10. Гидродинамика и теплообмен в атомных энергетических уста-
- •11. Тепловые и атомные электрические станции: Справочник / Под
- •117 634050, Г. Томск, пр. Ленина, 36, рио тгу
3.3. Теплоотдача в однофазной среде при свободном движении жидкости (естественная конвекция)
Свободное движение возникает за счет изменения в объеме рас- сматриваемой жидкости массовых (объемных) сил. Такими силами яв- ляются сила тяжести, центробежная сила или силы, наведенные в жид- кости электромагнитным полем высокой напряженности. Наиболее хо- рошо изучено свободное движение жидкости, вызванное гравитацион- ными силами.
При теплообмене температура жидкости переменна. Вследствие этого возникает разность плотностей и, как следствие, разность грави- тационных сил, представляющих собой подъемную силу.
Рассмотрим свободное гравитационное течение в большом объеме вблизи поверхности пластины и цилиндра.
Пусть вертикальная пластина с неизменной температурой поверх- ности, равной tc, находится в жидкости или газе. Жидкость вдали от пластины неподвижна (вынужденное течение отсутствует), температура жидкости вдали от пластины постоянна и равна tж < tc. При этом у пла- стины появляется подъемное движение нагретого слоя жидкости. Вдали от пластины скорость равна нулю.
Слой нагреваемой движущейся жидкости является одновремен- но гидродинамическим и тепловым пограничным слоем, в пределах это-
го слоя изменяется скорость течения от нуля до максимума и темпера-
тура от tc до tж.
Рис. 3.2. Изменение коэффициента теплоотдачи
при свободном движении жидкости
Вначале толщина пограничного слоя мала и ее течение носит ла- минарный характер. Постепенно в движение вовлекается все большое
количество жидкости, толщина пограничного слоя растет, затем он раз- рушается (переходный режим) и возникает турбулентный режим тече- ния жидкости. При ламинарном режиме течения коэффициент теплоот- дачи с увеличением толщины пограничного слоя уменьшается (см. рис. 3.2), при переходном режиме – резко возрастает, при турбу- лентном режиме – остается постоянным по высоте пластины.
На основе математического описания процесса конвективного те- плообмена при естественной конвекции выявлена структура уравнения подобия:
Nu = f(Gr, Pr). (3.4)
В результате экспериментального исследования теплоотдачи ус-
тановлено, что при
103 Gr
Pr
109
– ламинарный режим течения
жx ж
жидкости в пограничном слое;
109 Gr
Pr
6 1010
– переходный
жx ж
режим течения жидкости; Gr
Pr
6 1010
– турбулентный режим те-
чения жидкости.
жx ж
Число Грасгофа рассчитывается по формуле:
g t x3
Grжх 2 , (3.5)
где – температурный коэффициент объемного расширения;
t tc tж
– разность температур. Для капельных жидкостей значения
=
f(t)
приводятся
в
справочниках,
для
газов
рассчитывается
по
фор-
муле
1
Tж
, полученной на основе совместного решения уравнений
1 v
и pv RT .
v T p
Значения числа Прандтля, зависящие только от теплофизических свойств жидкостей, приводятся для различных теплоносителей (жидко- стей, газов) в справочной литературе.
Для расчета коэффициентов теплоотдачи
Nu жx
x
рекоменду-
ются следующие уравнения при свободном движении жидкости вдоль вертикальной поверхности:
при ламинарном режиме локальные значения , описываемые кривой = f(x) (рис. 3.2) в ламинарной области пограничного слоя, рас-
считываются по уравнению
Nu 0.60(Gr
Pr )0.25 Prж
0.25
; (3.6)
жх жх ж
Prc
средние коэффициенты теплоотдачи на участке поверхности с высо- той h с ламинарным течением в пограничным слое, рассчитываются по уравнению
Nu 0.75(Gr
Pr )0.25 Prж
0.25
; (3.7)
жh жh ж
Prc
при турбулентном режиме коэффициенты теплоотдачи рассчи- тываются по уравнению
Nu 0.15(Gr
Pr )1 / 3 Prж
0.25
; (3.8)
жх жх ж
Prc
при переходном режиме средний коэффициент теплоотдачи можно определить по формуле
пер
тур лам . (3.9)
2
Здесь определяющей является температура жидкости за предела- ми движущегося слоя (Prc выбирается по местной температуре стенки). Определяющий размер – высота пластины, отсчитываемая от места на- чала теплообмена.
Коэффициент теплоотдачи зависит от направления теплового по- тока, и обусловлено это неодинаковыми средними температурами жид- кости вблизи поверхности при нагреве (tс > tж) или охлаждении (tс < tж) жидкости, а также зависимостью теплофизических свойств жидкости от ее температуры. Как следствие этого, коэффициент теплоотдачи ка- пельных жидкостей при нагреве больше, чем при охлаждении. Влияние указанного фактора учитывается в уравнениях подобия множителем
Prж
0.25
. Так при нагреве жидкости
Prж
0.25
> 1, при охлаждении
Prc
0.25
0.25
Prc
Prж
< 1. Отношение
Prж
при течении определенной капель-
Prc
Prc
ной жидкости тем более отличается от единицы, чем больше темпера- турный напор.
Для газов с достаточной точностью можно считать, что
0.25
Prж
= 1.
Prc
Форма поверхности при естественной конвекции жидкости играет
второстепенную роль (важна ее протяженность), поэтому по вышепри- веденным формулам рассчитывается теплоотдача от плоских, цилинд- рических или вертикальных поверхностей иной формы.
Для горизонтальных плоских поверхностей (пластин) приведен- ные выше формулы применимы для расчета, но в этом случае вычис- ленный коэффициент теплоотдачи надо увеличить на 30 %, если тепло- отдающая поверхность обращена вверх, или уменьшить на 30 %, если теплоотдающая поверхность обращена вниз. В качестве определяющего
размера берется меньшая сторона пластины.
Для горизонтальных труб при ламинарном режиме течения жид-
кости
для
расчета
среднего
коэффициента
теплоотдачи
пользуется следующее
уравнение:
Nu жd
d
ис-
Nu 0.5(Gr
Pr )0.25 Prж
0.25
. (3.10)
жd жd ж
Prc
где определяющий размер – наружный диаметр трубы d.
Естественная конвекция в ограниченном объеме характеризуется наличием восходящих и нисходящих потоков, когда условия свободно- го движения жидкости значительно отличаются от ее движения в неог- раниченном пространстве.
а)
Рис. 3.3. Условия возникновения свободного движения в ограниченном объеме
Через газовые прослойки передача теплоты между поверхностями осуществляется тремя способами: теплопроводностью, конвекцией и излучением. Через прослойки капельной жидкости – двумя: теплопро- водностью и конвекцией.
Во всех случаях передачу теплоты рассчитывают по формулам теории теплопроводности, но коэффициент теплопроводности среды
заменяют эквивалентным, учитывающим перенос теплоты другими спо- собами.
Для плоских прослоек тепловой поток рассчитывают по уравне-
нию
для цилиндрических
Q экв
t1 t2 F , (3.11)
2 l
Q
экв t1 t2 . (3.12)
ln d2
d1
ж
экв к , (3.13)
где – коэффициент теплопроводности жидкости; к – коэффициент, учитывающий перенос тепла конвекцией.
Для прослоек любой формы при рассчитывается по формуле
(Grж
Pr ) 103
коэффициент к
к 0.18(Grж
ж
g (t
t ) 3
t t
где Grж
1 2 , определяющая температура t
2 ж
1 2 .
2
ж
Pr ) 103 принимают, что к = 1.
Для газовых плоских прослоек
qл
, (3.15)
экв к
t t
1 2
где qл – плотность теплового потока, передаваемого излучением через газовую прослойку.
Для газовых цилиндрических прослоек
л
Q ln d2
d1
. (3.16)
экв к
2 l (t1 t2 )
