Полезная разность температур

3.1.3 Полезная разность температур

Общая полезная разность температур равна:

ΣΔt_п=Δt_п1+ Δt_п2

Полезные разности температур по корпусам (в °С) равны:

Δt_п1=t_г1-t_кl=142,9-122,82=20,08°С

Δt_п1=t_г1-t_кl=120, 3-63,29=57,01°С

Тогда общая полезная разность температуры равна:

ΣΔt_п=20,08+57,01=77,09 °С.

Проверим общую полезную разность температуры:

ΣΔt_п=t_г1-t_бк-(ΣΔ'+ΣΔ"+ΣΔ'")=142,9-47,42-(0,66+15,73+2)=77,09°С

3.1.4 Определение тепловых нагрузок

Расход греющего пара, производительность каждого корпуса по выпаренной воде и тепловые нагрузки по корпусам определим путем совместного решения уравнений тепловых балансов и уравнения баланса по воде для всей установки:

Q₁=D∙(I_г1-i₁)=1,03∙[G_н∙С_н∙(t_к1-t_н)+w ₁∙(Н_вп1-Св∙t_к1)+Q_конц1]; (3.6)

Q₂= w₁∙(I_г2-i₂)=1,03∙[(G_н- w₁)∙С₁∙(t_к2-t_к1)+w ₂∙(I_вп2-Св∙t_к2)+Q_конц2]; (3.7)

где 1,03 - коэффициент, учитывающий 3% потерь тепла в окружающую среду;

С_н, С₁ - теплоемкости растворов соответственно исходного и в первом корпусах. кДж/(кг*К); С_н =4.14кДж/(кг*К), С₁=3.994кДж/(кг*К).Св- теплоемкость воды, кДж/(кг*К).;t_н - температура кипения исходного раствора при давлении в корпусе,⁰С; D- расход греющего пара, кг/с;

При решении уравнения можно принять:

I_вп1 ≈ I_Г2; I_вп2 ≈ I_бк;

Q_конц1, Q_конц2-теплоты концентрирования по корпусам, кВт; так как эти величины имеют небольшое значение, то ими пренебрегаем.

Q₁=D∙(2744-104.6)=1,03∙[5∙4.14∙(122.82-122.6)+w₁∙(2711-4,19∙122.82)]

Q₂=w₁∙(2711-516.1)=1,03[(5-w₁)∙3,994∙(63.29-122.82)+w₂∙(2585-4.19∙63.29)]

W=w₁+w₂=3.95

Решение этой системы уравнений дает следующие результаты:

D=2.384 кг/с; w₁=1.859 кг/с; w₂=2.091 кг/с;

Q₁=6292 кВт; Q₂=4080 кВт

Результаты расчета сведены в табл. 1.

Таблица 1

3.1.5 Выбор конструкционного материала

Выбираем конструкционный материал, стойкий в среде кипящего раствора CuSO₄ интервале изменения концентраций от 4 %, до 19 % [6]. Легированные стали с содержанием никеля являются нестойкими в среде растворов CuSO₄. В этих условиях химически стойкой является сталь марки Х17 (5 балл стойкости). Скорость коррозии ее не менее 0,1 мм/год, коэффициент теплопроводности λ_ст=25,1 Вт/(м∙К).

3.1.6 Расчет коэффициента теплопередачи

Коэффициент теплопередачи для первого корпуса определяем по уравнению аддитивности термических сопротивлений:

К₁=(1/α₁+Σδ/λ+1/α₂)^-1. (3.8)

Примем, что суммарное термическое сопротивление равно термическому сопротивлению стенки δ_ст/λ_ст и накипи δ_н/λ_н

Термическое сопротивление загрязнений со стороны пара не учитываем. Получим:

Σδ/λ=δ_ст/λ_ст+δ_н/λ_н (3.9)

где δ_ст, δ_н - толщина стенки, толщина слоя накипи, м.

при δ_ст=0,002 м.

при δ_н=0,0005 м.

где λ_ст, λ_н - коэффициент теплопроводности стенки и накипи, Вт/(м∙К).

при λ_ст=25,1 Вт/(м*К).

при λ_н=2 Вт/(м*К).

Σδ/λ=0,002 /25,1+0,0005/2=2,87 ∙10^-4 м²*К/Вт.

Коэффициент теплоотдачи от конденсирующего пара к стенке α₁ равен:

α₁=2,04∙((r₁∙ρ_ж1²∙λ_ж1³)/(μ_ж1∙Н∙Δt₁))^1/4. (3.10)

где r₁ - теплота конденсации греющего пара, Дж/кг;

р_ж1, λ_ж1, μ_ж1 - соответственно плотность (кг/м³), теплопроводность Вт/(м*К), вязкость (Па*с) конденсата при средней температуре пленки:

t_пл=t_г1-Δt₁/2

где Δt₁ - разность температур конденсации пара и стенки, °С.

Расчет α₁ - ведем методом последовательных приближений.

В первом приближении примем Δt₁=2 °С. Тогда

t_пл=142.9-2/2=141,9°С.

α₁=2,04∙(2144∙10³ ∙1032²∙0,408³/0,19∙10^-3∙4∙2)^1/4=6484 Вт/(м²∙К)

Для установившегося процесса передачи тепла справедливо уравнение:

q=α₁∙Δt₁=Δt_ст/(ΣΔδ/λ)=α₂∙Δt₂, (3.11)

где q - удельная тепловая нагрузка, Вт/кв.м;Δt_cт - перепад температур на стенке, °С;Δt₂ - разность между температурой стенки со стороны раствора и температурой кипения раствора, °С.

Отсюда:

Δt_ст=α₁∙Δt₁∙(Σδ/λ)=6484∙2∙2,87∙10^-4 =3,72°С.

Тогда

Δt₂=Δt_п1-Δt_ст-Δt₁=20,08-3,72-2=16,36°С.

Коэффициент теплопередачи от стенки к кипящему раствору для пузырькового кипения в вертикальных кипятильных трубках при условии естественной циркуляции раствора равен:

α₂=А∙(q^0,6)=780∙(q^0,6)∙(λ₁^1,3)∙(ρ₁^0,5)∙(ρ_п1^0,06)/((с₁^0,3)∙(σ₁^0,5)∙(г_в1^0,6)∙(ρ₀^0,66)∙(μ₁^0,3))(3.12)

По справочной литературе определяем:

λ₁=0,4159 Вт/(м∙К); ρ₁=1068 кг/м³; ρ_п1=1,22 кг/м³; σ₁=0,067 Н/м; г_в1=2200∙10³ Дж/кг; ρ₀=0,529 кг/м³; с₁=4095 Дж/кг∙К; μ₁=0,265∙10^-3 Па∙с

Подставив эти значения, получим:

α₂=780∙(q^0,6)∙0,4159^1,3∙1068^0,5∙1,22^0,06/0,067^0,5∙(2200∙10³)^0,6∙0,529^0,66∙

4095^0,3∙(0,265∙10^-3)^0,3=7,408∙(6484)^0,6=1435 Вт/(м²∙К)

Проверим правильность первого приближения по равенству удельных тепловых нагрузок:

q₁=α₁∙Δt₁=6484∙2=12968 Вт/кв.м

q₂=α₂∙Δt₂=1435∙16,36=2348 Вт/кв.м

q_1≠q₂

Для второго приближения примем Δt₁=5,0 град

Пренебрегая изменением физических свойств конденсата при изменении температуры на 3,0 град, рассчитаем α₁ по соотношению:

α₁=6484∙(2/5)^1/4=5156 Вт/(м²∙К)

Получим:

Δt_ст=5156∙5∙2,87∙10^-4=7,4 град;

Δt₂=20,08-5-7,4=7,68 град;

α₂=7,408*(5156∙5)^0,6=3285 Вт/(м²∙К)

q₁=5156∙5=25780 Вт/м²

q₂=3285∙7,68=25229 Вт/м²

q_1≈q₂

Расхождение между тепловыми нагрузками не превышает 3%, следовательно, расчет коэффициентов α₁ и α₂ на этом можно закончить.

Находим К₁:

К₁=(1/5156+2,87∙10^-4+1/3285)^-1=1271 Вт/(м²∙К).

Далее рассчитываем коэффициент передачи для второго корпуса К₂.

В первом приближении примем Δt₁=4 °С. Тогда:

Δt_пл=120. 3-2/2=118.3°С

α₁=2,04∙(2210∙10³ ∙1133∙0.4265/4∙4∙0.335∙10^-3)^1/4=5164 Вт/м²К

Δt_ст=5164∙4∙2,87∙10^-4=5,93°С

Δt₂=57,01-4-5,93=47,08°С

α₂=780∙(q^0,6)∙0,4366^1/3∙1187^0,5∙0,15^0,06/0,096^0,5∙(2350*10³)^0,6 0,529^0,66 ∙3509^0,3∙(0,851∙10^-3)^0,3 = 4,34(388∙4)^0,6=1683 Вт/м²∙К

Проверим правильность первого приближения по равенству удельных тепловых нагрузок:

q₁=α₁∙Δt₁=5164∙4=20656 Вт/м²

q₂=α₂∙Δt₂=1683∙47,08=79236 Вт/м²

q_1≠q₂

Используя вышеописанный метод приближения, найдем:

Δt₁=18.65°С

α₁=5164∙(4/18,65)^1/4=3514 Вт/м²К

Δt_ст=3514∙18,65∙2,87∙10^-4=18,81°С

Δt₂=57,01-18,81-18,65=19,55°С

α₂=4,34∙(3514∙18,65)^0,6=3368 Вт/м²∙К

q₁=65536 Вт/м²

q₂=65845 Вт/м²

q₁≈q₂

Определим К₂:

К₂=(1/3514+2,87*10^-4+1/3368)^-1=1151 Вт/м²∙К