- •«Томский политехнический университет» а.Г. Коротких, и.В. Шаманин основы гидродинамики и теплообмена в ядерных реакторах
- •Содержание
- •Глава 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, рио тгу
1.4. Статическое давление
Любое статическое давление в жидкости вызывает усилие, пер- пендикулярное стенкам сосуда. Такие силы могут быть очень большими и должны учитываться при проектировании. Если обозначить через z высоту над поверхностью постоянного давления, то разность давлений между высотами z2 и z1 в покоящейся жидкости (при отсутствии трения и ускорения) равна
p2 p1
z1 gdz
g z1
z2
. (1.18)
z 2 g c gc
Выражение (1.18) полезно для вычисления основных поправок на неизотермичность манометров и т.п.
Интегрирование по замкнутому контуру дает
g z1
g z1
pB
a dz
b dz
a b dz
0 0 tb ta dz . (1.19)
gc z 2
z1
g c z 2
g c z 2
1.5. Уравнение теплопроводности
Использование при расчетах закона сохранения энергии чрезвы- чайно полезно в тех случаях, когда скорость потока жидкости, его рас- ширение и теплой поток весьма высоки, что имеет место в активной зо- не ядерных реакторов. Рассмотрим элемент массы, мгновенно зани- мающий объем dxdydz с центром в точке x, y, z за элементарный проме- жуток времени d (см. рис. 1.3).
Для облегчения вывода дифференциального уравнения сделаем следующие допущения: тело однородно и изотропно; физические пара- метры постоянны; деформация рассматриваемого объема, связанная с
изменением температуры, является очень малой величиной по сравне- нию с самим объемом; макроскопические частицы тела неподвижны друг относительно друга; внутренние источники тепла в теле, которые в
общем случае могут быть заданы как равномерно.
qv
f x, y, z,
, распределены
b) Элемент объема dxdydz
Рис. 1.3. К выводу уравнения теплопроводности
В основу вывода дифференциального уравнения теплопроводно- сти положен закон сохранения энергии, который в рассматриваемом случае может быть сформулирован следующим образом: количество те- пла dQ, введенное в элементарный объем извне за время dτ вследствие теплопроводности, а также от внутренних источников, равно измене- нию внутренней энергии вещества, содержащегося в элементарном объ- еме:
dQ1 dQ2 dQ , (1.20)
где dQ1 – количество тепла, введенное в элементарный объем путем те- плопроводности за время dτ; dQ2 – количество тепла, которое за время dτ выделилось в элементарном объеме dV за счет внутренних источни- ков; dQ – изменение внутренней энергии вещества, содержащегося в элементарном объеме dV, за время dτ.
Для нахождения составляющих уравнения (1.20) выделим в теле элементарный параллелепипед со сторонами dх, dу, dz (см. рис. 1.3b). Параллелепипед расположим так, чтобы его грани были параллельны соответствующим координатным плоскостям.
Количество тепла, которое подводится к граням элементарного объема за время dτ в направлении осей ох, оу, оz, обозначим dQx, dQy, dQz, соответственно.
Количество тепла, которое будет отводиться через противопо- ложные грани в тех же направлениях, обозначим соответственно dQx+dx, dQy+dy, dQz+dz. Количество тепла, подведенное к грани dydz в направле- нии оси ох за время dτ, составляет dQx qx dydzd , где qх – проекция плотности теплового потока на направление нормали к указанной грани. Количество тепла, отведенного через противоположную грань элемен-
тарного параллелепипеда в направлении оси ох, запишется как
dQxdx qxdx dydzd .
Разница количеств тепла, подведенных к элементарному паралле- лепипеду и отведенных от него за время dτ в направлении оси ох пред- ставляет собой количество тепла dQx1:
dQx1 dQx dQxdx . (1.21)
Функция qх+dx является непрерывной в рассматриваемом интерва- ле dх и может быть разложена в ряд Тейлора:
q q
q
x dx
x
x2
... .
xdx x x
x2 2!
Если ограничиться двумя первыми членами ряда, то уравнение
(1.21) запишется в следующем виде:
dQx1
qx
dxdydzd
. (1.22)
x
Аналогичным образом можно найти количество тепла, подводи- мое к элементарному объему и в направлениях двух других координат- ных осей оу и оz.
Количество тепла dQ1, подведенное теплопроводностью к рас- сматриваемому объему, будет равно:
dQ qx qy
qz
dxdydzd
. (1.23)
1
x y
z
Определим вторую составляющую уравнения (1.20) dQ2. Обозна- чим удельную производительность внутренних источников тепла через qv, [Вт/м3]. Удельную производительность внутренних источников на- зывают также «объемной плотностью тепловыделения». Объемная плотность – это количество тепла, которое выделяется в единице объема вещества в единицу времени. Тогда
dQ2 qv dVd . (1.24)
Третья составляющая уравнения (1.20), характеризующая измене- ние внутренней энергии, может быть найдена по известному уравнению
dQ
c
t
dVd . (1.25)
Подставляя полученные выражения (1.23), (1.24) и (1.25) в урав- нение (1.20), получаем:
t
1 qx qy
qz
qv
. (1.26)
c
x y
z c
Проекции вектора плотности теплового потока на координатные оси ох, оу, оz определяются выражениями (закон Фурье)
q
t
;
x x
q
t
;
y y
q
t
.
z z
Подставляя
полученные
выражения
проекций
вектора
плотности
теплового потока в уравнение
(1.26),
получаем:
2 2 2
t
t
t
t qv
. (1.27)
c
x2
y2
z 2 c
Если в уравнении (1.27) обозначить
a c
и
2
2
2
t
x2
t
y 2
t
z 2
2t ,
где а – коэффициент температуропроводности, м2/с;
2t
– оператор Ла-
пласа в декартовой системе координат, то получим уравнение тепло- проводности в общем виде
t a2t
qv . (1.28)
c
Уравнение теплопроводности является частным случаем уравне- ния энергии:
Dt
a2t
qv
d c
, (1.29)
в котором полная (субстанциональная) производная по времени D
d
учитывает процесс переноса тепла, вызванный движением имеющих
различную температуру элементов потока жидкой (газообразной) сре- ды:
D
w
w
w
,
d
x x
y y
z z
где wx, wy, wz – проекции вектора скорости потока на соответствующие оси декартовой системы координат.
Выражение 2t
в цилиндрической системе координат имеет вид:
2 2 2
2t
t
1 t 1
t t ,
r 2
r r r 2 2
z 2
где r – радиальная, – угловая и z – аксиальная (осевая) координаты, соответственно.
Выражение 2t
в сферических координатах имеет вид:
2t
2 t
1 t
t
1 2t
2t
1 2
,
r 2
r r r 2
r 2 1 2
2
где
cos ; и – угловые координаты.
Уравнение (1.28) называется дифференциальным уравнением теп-
лопроводности. Оно устанавливает связь между временным и простран- ственным изменением температуры в любой точке тела, в котором про- исходит процесс теплопроводности.
Коэффициент температуропроводности а является физическим параметром вещества, используется при описании нестационарных теп-
ловых процессов и характеризует скорость изменения температуры. Ес- ли коэффициент теплопроводности характеризует способность тел про- водить тепло, то коэффициент температуропроводности является мерой теплоинерционных свойств тела. Из уравнения (1.27) следует, что изме-
нение температуры во времени
t /
для любой точки пространства
пропорционально величине а. Иначе говоря, скорость изменения темпе- ратуры в любой точке тела будет тем больше, чем больше коэффициент температуропроводности а. Поэтому, при прочих равных условиях, вы- равнивание температур во всех точках пространства будет происходить быстрее в том теле, которое характеризуется большим коэффициентом температуропроводности. Величина коэффициента температуропровод- ности зависит от природы вещества. Например, жидкости и газы обла-
дают большой тепловой инерционностью и, следовательно, малым ко- эффициентом температуропроводности. Металлы обладают малой теп- ловой инерционностью, так как они имеют большой коэффициент тем- пературопроводности.
