4.4 Корисна різниця температур
Корисні різниці температур по корпусам розраховуємо за формулою:
∆t кор = t г – t к ; (4.21)
∆t кор. 1 = 110,11 – 119,632 = -9,522 ºС;
∆t кор 2 = 101,96 – 90,359 = 11,601 ºС,
∆t кор 3 = 87,59 – 74,369 = 13,221 ºС.
Загальну корисну різницю температур розраховуємо за формулою:
Σ ∆ tкор = ∆t кор 1 + ∆t кор. 2 + ∆t кор 3 . (4.22)
Σ ∆ tкор = - 9,522 + 11,601 + 13,221 = 15,3 С.
Загальну корисну різницю температур перевіряємо за формулою:
∆ tкор = t г 1 – tбк – (Σ∆' + Σ∆" + Σ∆"'), (4.23)
Σ∆ tкор = 110,11 – 68,34 – (8,81 + 14,66 + 3) = 15,3 0С.
4.5 Визначення теплових навантажень
Витрати гріючої пари в 1- й корпус , продуктивність кожного корпусу по випареній воді і теплові навантаження по корпусам визначимо шляхом загального рішення рівнянь теплових балансів по корпусам і рівняння балансу по воді для всієї установки:
Q1 = D (Iг 1 – i1) = 1.03[Gп cп (t к 1 – t п) + w1 (Iвп 1 – св t к 1)]; (4.24)
Q2 = w1 (Iг 2 – i2) = 1.03[(Gп - w1) c1(t к 2 – t к 1) + w2(Iвп 2 – св t к 2)]; (4.25)
Q3 = w2 (Iг 3 – i3) = 1.03[(Gп - w1 – w2) c2 (t к 3 – t к 2) +w3(Iвп 3 – св t к 3)]; (4.26)
W = w1 + w 2+ w3, (4.27)
де 1,03 – коефіцієнт, який враховує 3% витрат в оточуюче середовище;
сп,с1,с2 – теплоємність розчину відповідно вихідного в першому і в другому корпусах, кДж/(кг·К) [9, с.157];
t п - температура кипіння вихідного розчину при тиску в першому корпусі;
tв – теплоємність води, К);
D – витрати гріючої пари, кг/с;
і – ентальпія води, кДж/кг, [2, табл.LVII, с.532].
Температуру кипіння вихідного розчину при тиску в першому корпусі розраховуємо за формулою:
tп = tвп 1 + ∆'п, (4.28)
де ∆'п – температурна депресія для вихідного розчину, ºС.
tп = 102,96+ 1,0 = 103,96 ºС
При рішенні рівняння (4.24) – (4.27) приймаємо, [1,с.170]:
Iвп 1≈Іг 2; Iвп 2≈Іг 3; Iвп 3≈Ібк. (4.29)
Отримаємо систему рівнянь:
Q1 = D·(2700 - 491,12) = 1,03·[3,5·3,207·(119,632 – 102,96) + w1·(2681 - 4,19·119,632)];
Q2=w1 (2681 - 371,77)=1,03·[(3,5 – w1)·3,544·(90,359–119,632)+w2·(2653 – 4,19·90,359)];
Q3= w2·(2653–290,7) =1.03 [(3,5- w1–w2)·2,858 (74,369–90,359) + w3(2608–4,19·74,369)].
W = w1 + w2 + w3 = 3.
Рішення цієї системи рівнянь дає наступні результати:
D = 0,595 кг/с; w1 = 0,482 кг/с; w2 = 0,61 кг/с; w3 = 1,236 кг/с;
Q1 = 1314,28 кВт; Q2 = 1113,05 кВт; Q3 = 1441,003 кВт.
Невелике відхилення навантаження по випарюваній воді в кожному корпусі від попередньо прийнятих ( , , ) не перевищє 3,5 %, том не будемо перераховувати концентрації та температури кипіння розчинів у корпусах.
Результати розрахунку заносимо в таблицю 4.3:
Таблиця 4.3 – Результати розрахунку параметрів у корпусах
Параметри |
1 корпус |
2 корпус |
3 корпус |
Продуктивність по випарюваній воді,w, кг/с |
0,94 |
0,998 |
1,061 |
Концентрація розчинів x, % |
12,22 |
16,18 |
25 |
Тиск гріючих парів Рг, МПа |
0,2 |
0,1443 |
0,0866 |
Температура гріючих парів tг,ºС |
120,119 |
109,920 |
95,656 |
Температура кипіння розчину tк,ºС |
115,862 |
103,316 |
83,722 |
Корисна різниця температур ∆tкор ,ºС |
4,257 |
6,604,58 |
11,934 |
4.6 Вибір конструкційного матеріалу
Вибираємо конструкційний матеріал, стійкий в середовищі киплячого розчину NаCl з інтервалом зміни концентрації від 10 до 25 % В цих умовах хімічно стійкою є сталь Х18Н10Т. Швидкість корозії її не більше 0,1 мм/рік, коефіцієнт теплопровідності λст = 16,4 Вт/(м·К).
4.7 Розрахунок коефіцієнтів теплопередачі
Коефіцієнт теплопередачі розраховують за рівнянням адитивності термічних опорів:
, (4.30)
де α1 і α2 – коефіцієнти тепловіддачі від пари, що конденсується, до стінки і від стінки до киплячого розчину відповідно, Вт/м2∙К;
- загальний термічний опір, м2∙К/Вт.
(4.31)
Приймаємо, що сумарний термічний опір, рівний термічному опору стінки і накипу . Приймаю для всіх корпусів товщину шару накипу δн=0,5мм, λн=3,05 Вт/(м∙К) [1, с. 172]. Термічний опір зі сторони пари не враховуємо. Отримаємо:
м2 К/Вт.
Коефіцієнт тепловіддачі від пари, що конденсується, до стінки розраховуємо за формулою:
(4.32)
де r – теплота конденсації гріючої пари, Дж/кг;
ρр – густина конденсату при середній температурі плівки, кг/м3 [2,табл. XXXIX с. 520];
μр – в’язкість конденсату при середній температурі плівки, Па∙с, [2,табл. XXXIX с. 520];
λр – теплопровідність конденсату при середній температурі плівки, Вт/(м∙К), [2,табл. XXXIX с. 520].
Середню температуру плівки конденсату розраховуємо за формулою:
, (4.33)
де Δt1 – різниця температур конденсації пари і стінки, ºС.
В першому наближенні приймаємо Δt1 = 1 ºС.
ºС
Приймаємо r = 2206,55 кДж/кг; ρр = 943,305 кг/м3; μр = 231,953∙10-6 Па∙с; λр = 68,596∙10-2 Вт/(м∙К) [2].
Розрахунок α1 ведуть методом послідовних наближень за формулою (4.32):
Вт/(м2∙К).
Для встановленого процесу передачі тепла справедливе рівняння:
, (4.34)
де q – питоме теплове навантаження, Вт/м ;
Δtст. – перепад температур на стінці, ºС;
Δt2 – різниця між температурою стінки із сторони розчину і температурою кипіння розчину, ºС
Перепад температур на стінці розраховуємо за формулою:
, (4.35)
ºС
Різницю між температурою стінки із сторони розчину і температурою кипіння розчину розраховуємо за формулою:
(4.36)
ºС
Коефіцієнт тепловіддачі від стінки до киплячого розчину для пузиркового кипіння в вертикальних киплячих трубках при умовах природної циркуляції розчину розраховуємо за формулою:
, (4.37)
де λ – теплопровідність розчину при температурі кипіння, Вт/(м∙К);
ρ – густина розчину при температурі кипіння, кг/м3;
ρп – густина пари при температурі кипіння, кг/м3;
σ – поверхневий натяг, Н/м;
rв – теплота пароутворення, Дж/кг;
с – теплоємність розчину при температурі кипіння, Дж/кг∙К;
μ – в’язкість розчину при температурі кипіння, Па∙с;
ρ0 – густина водяних парів при нормальних умовах, кг/м3, [1, с.172].
Фізичні властивості киплячих розчинів NаCl та їх парів зводимо у таблицю.
Таблиця 4.4 – Фізичні властивості киплячих розчинів NаCl та їх парів
Параметр |
Корпус |
Література |
||
І |
ІІ |
ІІІ |
||
Теплопровідність розчину λ, Вт/(м∙К) |
0,675 |
0,6719 |
0,6370 |
[10, c.330], [9, c.157] |
Густина розчину ρ, кг/м3 |
1020,3 |
1049,24 |
1178,25 |
[3, c.54] |
Теплоємність розчину с, Дж/кг∙К |
3610.16 |
3544,2 |
2858 |
[10, c.53] |
Продовження таблиці 4.4 |
||||
В’язкість розчину μ, Па∙с |
0,28556 |
0,3651 |
0,542524 |
[10, c.171] |
Поверхневий натяг σ, Н/м |
0,0743 |
0,0756 |
0,0789 |
[9, c.157] |
Теплота пароутворення rв, Дж/кг |
2207,9∙103 |
2283,7∙103 |
2322,46∙103 |
[2, табл.LXII, c.532] |
Густина пару ρп, кг/м3 |
1,1082 |
0,4305 |
0,2367 |
[2, табл.LXII, c.532], |
Підставивши числові значення, отримаємо:
;
Вт/(м2.К).
Перевіримо достовірність першого наближення по рівності питомих теплових навантажень:
Вт/м2;
Вт/м 2,
q'≠q''.
Для другого наближення приймаємо Δt1=0,36 ºС.
Нехтуючи зміною фізичних властивостей конденсату при зміні температури на 0,64 ºС, розраховуємо α1 по співвідношенню:
(Вт/(м2.К);
ºС;
ºС;
Вт/(м2.К);
Вт/м2;
Вт/м2;
Як бачимо q'≈q''. Розходження між тепловими навантаженнями не перевищує 3%, тому розрахунок коефіцієнтів та на цьому завершуємо.
Коефіцієнт теплопередачі для першого корпусу розраховуємо за формулою (4.30):
Вт/(м2.К)
Далі розраховуємо коефіцієнт теплопередачі для другого корпусу.
В першому наближенні приймаємо Δt1 = 0,8 ºС.
Середню температуру плівки конденсату розраховуємо за формулою (4.33):
ºС
Приймаємо r = 2234,95 кДж/кг; ρр = 951,336 кг/м3; μр = 257,248∙10-6 Па∙с; λр = 68,49∙10-2 Вт/(м∙К) [2].
Коефіцієнт тепловіддачі від пари, що конденсується, до стінки розраховуємо за формулою (4.32):
Вт/(м2∙К).
Перепад температур на стінці розраховуємо за формулою (4.35):
ºС
Різницю між температурою стінки із сторони розчину і температурою кипіння розчину розраховуємо за формулою (4.36):
ºС
Коефіцієнт тепловіддачі від стінки до киплячого розчину розраховуємо за формулою (4.37):
Вт/(м2.К).
Перевіримо достовірність першого наближення по рівності питомих теплових навантажень:
Вт/м2,
Вт/м2,
q'≠q''.
Для другого наближення приймаємо Δt1=0,73 ºС.
Нехтуючи зміною фізичних властивостей конденсату при зміні температури на 0,07 ºС, розраховуємо α1 по співвідношенню:
(Вт/(м2.К);
ºС;
ºС;
Вт/(м2.К);
Вт/м2;
Вт/м2;
Як бачимо q'≈q''. Розходження між тепловими навантаженнями не перевищує 3%, тому розрахунок коефіцієнтів та на цьому завершуємо.
Коефіцієнт теплопередачі для другого корпусу розраховуємо за формулою (4.30):
Вт/м2.К
Розраховуємо коефіцієнт теплопередачі для третього корпусу.
В першому наближенні приймаємо Δt1 = 2 ºС.
Середню температуру плівки конденсату розраховуємо за формулою (4.33):
ºС
Приймаємо r = 2273,86 кДж/кг; ρр = 961,741 кг/м3; μр = 299,635∙10-6 Па∙с; λр = 68,140∙10-2 Вт/(м∙К) [2].
Коефіцієнт тепловіддачі від пари, що конденсується, до стінки розраховуємо за формулою (4.32):
Вт/(м2∙К).
Перепад температур на стінці розраховуємо за формулою (4.35):
ºС
Різницю між температурою стінки із сторони розчину і температурою кипіння розчину розраховуємо за формулою (4.36):
ºС
Коефіцієнт тепловіддачі від стінки до киплячого розчину розраховуємо за формулою (4.37):
Вт/(м2.К).
Перевіримо достовірність першого наближення по рівності питомих теплових навантажень:
Вт/м2;
Вт/м2;
q'≠q''.
Для другого наближення приймаємо Δt1=1,9 ºС.
Нехтуючи зміною фізичних властивостей конденсату при зміні температури на 0,1 ºС, розраховуємо α1 по співвідношенню:
(Вт/(м2.К);
ºС;
ºС;
Вт/(м2.К);
Вт/м2;
Вт/м2;
Як бачимо q'≈q''. Розходження між тепловими навантаженнями не перевищує 3%, тому розрахунок коефіцієнтів та на цьому завершуємо.
Коефіцієнт теплопередачі для третього корпусу розраховуємо за формулою (4.30):
Вт/(м2.К)