Расчёт коэффициентов теплопередачи
Коэффициент теплопередачи для первого корпуса определяют по уравнению аддитивности термических сопротивлений:
(4)
Примем, что суммарное термическое сопротивление равно термическому сопротивлению стенки δст/λст и накипи δн/λн. Термическое сопротивление загрязнений со стороны пара не учитываем. Получим:
Коэффициент теплоотдачи от конденсирующегося пара к стенке α1 равен:
где r1 - теплота конденсации греющего пара, Дж/кг; ρж, λж, μж - соответственно плотность (кг/м 3), теплопроводность Вт/(м·К), вязкость (Па·с) конденсата при средней температуре плёнки, tпл= tг.1. – ∆t1/2, где ∆t1 - разность температур конденсации пара и стенки, град.
Первый корпус
Расчёт ∆t1 ведут методом последовательных приближений. В первом приближении примем ∆t1=1,45 град. Тогда:
r1 =2017,06 кДж/кг ρж = 887,87 кг/м3
μж = 15,087∙10-5 Па∙с λж = 0,6737 Вт/(м∙К)
Для установившегося процесса передачи тепла справедливо уравнение:
где q - удельная тепловая нагрузка, Вт/м 2; Δtст - перепад температур на стенке, град; Δt2 - разность между температурой стенки со стороны раствора и температурой кипения раствора.
Отсюда:
Тогда:
Коэффициент теплоотдачи от стенки к кипящему раствору для пузырькового кипения в вертикальных кипятильных трубках при условии принудительной циркуляции раствора равен:
(5)
(6),
(7), (8)
Отсюда выразим α2:
(9)
Физические свойства кипящего раствора NaCl сведём в таблицу 5:
Таблица 5
Параметр |
Корпус |
||
1 |
2 |
3 |
|
Теплопроводность раствора λ, Вт/(м·К) |
0,7005 |
0,6935 |
0,644 |
Плотность растворов ρ, кг/м3 |
955,25 |
968,1 |
1072,2 |
Теплоемкость раствора с, Дж/ (кг·К) |
4301 |
4200 |
3659 |
Вязкость раствора µ·10-3, Па·с |
0,3054 |
0,4637 |
0,5820 |
Скорость движения раствора принимаем равной v=2 (м/с). Подставив численные значения, получим:
Проверим правильность первого приближения по равенству удельных тепловых нагрузок:
q' = α1∙∆t1
q' = 8 949,6064∙1,45= 12 976,9293 Вт/м2
q'' = α2∙∆t2
q'' = 1 839,2087∙9,08=16 700,015 Вт/м2
Как видим, q' ≠ q''
Расхождение между тепловыми нагрузками превышает 3%. Для второго приближения примем ∆t1=3 град.
Пренебрегая изменением физических свойств конденсата при изменении температуры на несколько градусов, рассчитаем α1 по соотношению:
Получим:
q' = 7 462,1783∙3= 22 386,5349 Вт/м2
q'' = 1 839,2087∙3,69=6 786,6801 Вт/м2
Как видим, q' ≠ q''
Для расчёта в третьем приближении строим графическую зависимость удельной тепловой нагрузки от разности температур между паром и стенкой в первом корпусе (рис. 1) и определяем ∆t1= 1,75 град. Получим:
Рис. 1. Зависимость удельной тепловой нагрузки q от разности температур t в корпусе
q' = 8 538,5952∙1,75=14 942,5416 Вт/м2
q'' = 1 839,2087∙6,36=14 697,3673/м2
Как видим, q' ≈ q''
Находим коэффициент теплопередачи:
Второй корпус
В первом приближении примем t = 7 ˚С. Тогда
tпл = tг.2 – ∆t1 /2
tпл = 163,2 – 7/2 = 159,7 °С
r2 =2103,72 кДж/кг ρж = 913,99 кг/м3
μж = 17,84∙10-5 Па∙с λж = 0,6815 Вт/(м∙К)
Тогда:
Проверим правильность первого приближения по равенству удельных тепловых нагрузок:
q' =2 357,4∙7 = 16 501,8 Вт/м2
q'' = 1 942,16∙8,41 = 16 333,57 Вт/м2
Как видим, q' ≈ q''
Расхождение между тепловыми нагрузками не превышает 3%, расчёт коэффициентов теплоотдачи на этом заканчиваем.
Находим коэффициент теплопередачи:
Третий корпус
В первом приближении примем t = 16 ˚С. Тогда
tпл = tг.3 – ∆t1 /2
tпл = 137,9 – 16/2 = 129,9 °С
r3 =2173,99 кДж/кг ρж = 934,92 кг/м3
μж = 21,31∙10-5 Па∙с λж = 0,6837 Вт/(м∙К)
Тогда:
Проверим правильность первого приближения по равенству удельных тепловых нагрузок:
q' =4 761,78∙16 = 76 188,48 Вт/м2
q'' = 3 094,73∙24,5 = 75 820,89 Вт/м2
Как видим, q' ≈ q''
Расхождение между тепловыми нагрузками не превышает 3%, расчёт коэффициентов теплоотдачи на этом заканчиваем.
Находим коэффициент теплопередачи:
