Расчет процессов теплопроводности и конвективного теплообмена. методические указания к выполнению курсовой работы. Мелькумов В.Н., Петрикеева Н.А
.pdfпоглощательной способностью обладают многоатомные газы в продуктах сгорания топлива, а именно СО2 и Н2О.
Определяем парциальные давления трехатомных газов в продуктах сгорания:
P |
= Pr r |
co |
2 |
≈ 1 r |
co 2 |
= r |
co 2 |
, |
(1.24) |
co 2 |
|
|
|
|
|
||||
Pн2о |
= Pr rн2 |
о |
≈1 rн2о |
= rн2о |
, |
(1.25) |
где rСО2 , rН2О - соответственно объемное содержание углекислоты и водяных паров в дымовых газах (см. состав дымовых газов в задании), доли;
Pr - давление дымовых газов, бар.
Определяем эффективную толщину излучающего слоя:
а) для эффективных экономайзеров и теплообменников с кипящим теплоносителем для гладкотрубных пучков:
4 |
|
|
, |
(1.26) |
||
lэф = 0,9Dн |
|
σ1σ2 |
−1 |
|||
π |
||||||
|
|
|
|
|
б) для трубчатых воздухоподогревателей, при течении газов внутри труб:
lэф = 0.9Dвн , |
(1.27) |
где σ1 , σ2 - соответственно относительный поперечный и продольный шаг трубного пучка.
Вычисляем параметр Plэф :
Pco2 lэф |
= rco2 lэф , |
(1.28) |
Pн2о lэф |
= rн2о lэф . |
(1.29) |
По параметру Plэф и температуре t1с с помощью номограмм рис. П.1.1 и рис. П.1.2 определяем степень черноты объемов излучающих газов εСО2 и εН2О.
Степень черноты излучающих газов определяется по формуле
εr = εco2 + β εн2о , |
(1.30) |
где β - поправочный коэффициент, определяемый по рис. П.1.3. Выбираем степень черноты поверхности труб εС в зависимости от t1с по
табл. П.1.6.
Приведенная степень черноты системы « дымовые газы - труба»:
εпр = |
εr εc |
|
εc +εr (1−εc ). |
(1.31) |
11
Удельный тепловой поток излучением от дымовых газов к стенке трубы:
|
|
|
|
|
4 |
|
|
Т |
|
|
4 |
|
|
|
q1л |
T |
|
|
1с |
|
|
, |
(1.32) |
||||||
= 5,67εпр |
|
1ж |
|
|
− |
|
|
|
|
|
||||
100 |
|
100 |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где q1л - удельный тепловой поток излучением, Вт/м2; Т1ж - средняя температура греющей среды, К;
Т1с - средняя температура стенки трубы со стороны дымовых газов, К. С учетом вышеизложенного определяем величину лучистой
составляющей коэффициента теплоотдачи со стороны греющей среды:
α1л = |
q1л |
|
. |
(1.33) |
Т1ж − |
|
|||
|
Т1с |
|
1.3. ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТА ТЕПЛООТДАЧИ СО СТОРОНЫ НАГРЕВАЕМОЙ СРЕДЫ
1.3.1. Водяной экономайзер и воздухоподогреватель
Коэффициент теплоотдачи конвекцией α2 при вынужденном движении потока в трубах для водяного экономайзера и воздухоподогревателя определяется в зависимости от режима течения нагреваемой среды (воды или воздуха) из критериальных уравнений (1.12) - (1.13). В зависимости от схемы движения может быть определен по уравнениям (1.14) –(1.15).
При этом
Nu2ж = |
α2 D |
, |
(1.34) |
|||
|
|
|
λ |
|
|
|
|
|
|
2ж |
|
|
|
Re2 |
|
= |
W2ж D |
, |
(1.35) |
|
|
|
|||||
|
ж |
|
ν2ж |
|
|
|
где D = Dвн - при течении среды внутри труб;
D = Dн - при поперечном обтекании наружной поверхности труб. Физические свойства нагреваемой среды (λ2ж , Pr2ж , ν2ж) определяются
по табл. П.1.2 - для водяного экономайзера; по табл. П.1.1 - для
воздухоподогревателя при средней температуре |
|
t2ж = 0,5(t2′ +t2′ ). |
(1.36) |
Число Прандтля Pr2с определяется при температуре стенки |
t2с = t1с из |
вышеуказанных табл. П.1.1 и П.1.2. |
|
12
1.3.2. Теплообменник с кипящим теплоносителем
Коэффициент теплоотдачи α2 от внутренней стенки трубы к кипящему двухфазному пароводяному потоку определяется в следующей последовательности. Вначале определяем α2 w для случая однофазного потока некипящей воды из критериальных уравнений (1.12) или (1.13), выбираемых в зависимости от режима течения среды (числа Re2ж).
При этом:
Re2ж |
= |
W2ж Dвн |
, |
(1.37) |
||
|
|
|
||||
|
|
|
ν2ж |
|
||
Nu2ж |
= |
α2w Dвн |
, |
(1.38) |
||
|
||||||
|
|
|
λ2ж |
|
где λ2ж , ν2ж - определяются для воды по табл. П.1.2.
Число Прандтля Pr2ж определяется при температуре t2ж , Pr2с – при температуре стенки t2с = t1с по табл. П.1.2.
После этого вычисляем α2 q для случая развитого пузырькового кипения воды из формулы
α2q = |
3.4 P0.18 |
23 |
, |
(1.39) |
1−0.0045 P q |
|
|||
|
s |
|
|
|
s
где α2 q - коэффициент теплоотдачи к кипящей воде, Вт/(м2 ·град); Ps - давление насыщения, бар;
q- удельный тепловой поток через стенку трубы, Вт/м2.
Врезультате возможны следующие случаи:
а) при |
α2q |
|
≤ 0.5 |
|
α2 |
=α2w |
, |
|
|
|
|
|
(1.40) |
|||
α |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
2w |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
б) при |
α2q |
|
> 2 |
|
α2 |
=α2q , |
|
|
|
|
|
|
(1.41) |
|||
α |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
2w |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.5 < α2w |
|
|
|
|
4α2w +α2q |
|
|
|
|||||||
в) при |
< 2 |
α2 |
=α2w |
|
|
, |
(1.42) |
|||||||||
|
||||||||||||||||
|
|
|
α |
2q |
|
|
|
|
α |
2w |
+α |
2q |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
13
При определении α2 q по формуле (1.39) неизвестной является величина удельного теплового потока q. Поэтому вначале в первом приближении определяем α2 по формулам (1.37) – (1.42), приняв q = (1÷2) ·105 Вт/м2.
Затем, после вычисления коэффициента теплопередачи К и расчетной теплообменной поверхности Нр , производим уточнение значения q по формуле
q = |
Q |
, |
(1.43) |
|
Hp |
||||
|
|
|
После этого вновь определяем α2 с новым значением q по формулам
(1.37) – (1.42) и уточняем величины К и Нр.
Цикл последовательных приближений завершается при очередной поправке величины Нр , не превышающей 3 %.
1.3.3. Теплообменник с оребренной поверхностью теплообмена (отопительный прибор)
Внешняя теплообменная поверхность данного типа теплообменника представляет собой наружную поверхность трубы, снабженную квадратными ребрами. Коэффициент теплоотдачи со стороны оребренной поверхности теплообменника определяется по формуле
|
α2пр = |
α2 |
(E Fp + Fc ), |
(1.44) |
|
|
|||
|
|
Fрс |
|
|
где α2 - коэффициент теплоотдачи внешней поверхности трубы без учета |
||||
оребрения, Вт/(м2 ·град); |
|
|
|
|
Fpc = F1 χp - |
суммарная площадь |
оребренной поверхности |
теплообменника, м2; χр - коэффициент оребрения;
F1 - площадь внутренней (неоребренной) поверхности трубы, м2; Fр - суммарная площадь поверхности ребер, м2;
Fс – площадь гладкой части наружной поверхности трубы в промежутках между ребрами, м2;
Е - коэффициент эффективности ребра.
Коэффициент теплоотдачи α2 может быть определен из критериального
уравнения для свободной конвекции с внешней (горизонтальной поверхности) трубы диаметром Dн:
0,25 |
|
Pr |
0.25 |
|
|
|
Nu2ж = 0,5 (Gr2ж Pr2ж ) |
|
2ж |
|
, |
(1.45) |
|
Pr |
||||||
|
|
|||||
|
|
2c |
|
|
|
14
Nu2ж = |
α2 Dн , |
|
|
|
(1.46) |
|
|
λ |
|
|
|
|
|
|
2ж |
|
|
|
|
|
Gr2ж = g β ∆t |
|
Dн3 |
, |
(1.47) |
||
ν 2 |
||||||
|
|
|
|
|||
|
|
|
2ж |
|
|
где λ2ж - коэффициент теплопроводности воздуха, Вт/(м·град); ν2ж - кинематический коэффициент вязкости воздуха, м2/с; Pr2ж , Pr2с - число Прандтля для воздуха;
β = 1/Тв - температурный коэффициент объемного расширения воздуха, 1/К;
∆t = t2с – tв - температурный напор, 0С.
Физические свойства воздуха (λ2ж , ν2ж , Pr2ж) определяются при темпера-
туре t2с = t1с = t1ж по табл. П.1.1. Коэффициент эффективности ребра:
|
|
h |
|
|
|
|
|
|||||
|
th |
|
|
|
|
|
2B |
|
|
|
||
|
|
|
|
|
|
|
|
|||||
E = |
|
δ |
|
|
|
|
i |
, |
(1.48) |
|||
|
|
h |
|
|
2B |
|
||||||
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
δ |
|
|
|
i |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|||
Bi |
= |
α2 δ |
, |
|
|
(1.49) |
||||||
|
|
|
λ |
p |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
где h = 0.5(Dp − Dн ) - высота ребра, м; δ - толщина ребра, м;
λр - коэффициент теплопроводности материала ребра, Вт/(м·град), [1, табл. 6,7,8; 3, табл. П.1, П.2].
Зависимость Е от числа Вi и размеров ребра приведена в табл. 1.2.
Таблица 1.2 Значение коэффициента эффективности Е для квадратных ребер
(h/ δ) |
2ВΙ = 0 |
|
(h/ δ) |
2ВΙ |
= 0,5 |
(h/ δ) |
2ВΙ |
= 1,0 |
||||||
|
|
|
|
|
|
|
|
|
|
|||||
|
Dр / Dн |
|
|
Dр / Dн |
|
Dр / Dн |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.6 |
|
2.0 |
|
3.0 |
1.6 |
|
2.0 |
|
3.0 |
1.6 |
|
2.0 |
|
3.0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.0 |
|
1.0 |
|
1.0 |
0.84 |
|
0.85 |
|
0.8 |
0.57 |
|
0.57 |
|
0.54 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
15
Рекомендуемая последовательность решения задачи:
а) выбираем участок оребренной трубы с длиной l, кратной целому числу
шагов b; |
|
б) вычисляем значения Fр , Fс , Fрс и χр |
для участка трубы длиной l ; |
в) определяем площадь внутренней |
поверхности F1 участка трубы |
длиной l; |
|
г) находим значение α2 из уравнения (1.46);
д) находим значение числа Вi из уравнения (1.49), а затем коэффициент эффективности ребра Е из табл.1.2;
е) определяем по формуле (1.44).
1.4. ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТА ТЕПЛОПЕРЕДАЧИ ТЕПЛООБМЕННОГО АППАРАТА
1.4.1.Водяной экономайзер, воздухонагреватель и теплообменник
скипящим теплоносителем
Коэффициент теплопередачи характеризует процесс теплопередачи от потока к потоку через разделяющую их стенку.
Коэффициент теплопередачи Т.А. при переносе тепла через разделяющую стенку трубы (тонкостенный цилиндр) может быть приближенно определен по формуле для плоских теплообменных поверхностей:
K = |
1 |
|
|
, |
(1.50) |
||
|
1 |
+ε + |
1 |
|
|||
где К - коэффициент |
α1 |
α2 |
теплообменного аппарата, |
||||
|
теплопередачи |
Вт/(м2 град); ε - термическое сопротивление слоя отложений со стороны дымовых
газов, (м2 ·град)/Вт (см.задание); α1 , α2 - соответственно коэффициент теплоотдачи со стороны
греющей и нагреваемой сред, Вт/(м2 ·град).
В формуле (1.50) не учтено термическое сопротивление металлической стенки трубы ввиду его малости.
16
1.4.2. Теплообменник типа « отопительный прибор»
Коэффициент теплопередачи Т.А. без учета термического сопротивления металлической стенки трубы:
K = |
|
|
|
1 |
. |
(1.51) |
|
|
1 |
+ |
|
1 |
|||
|
α1 |
α2пр χр |
|
||||
|
|
|
1.5. ОПРЕДЕЛЕНИЕ СРЕДНЕГО ТЕМПЕРАТУРНОГО НАПОРА
Средний температурный напор для аппаратов с перекрестным или смешанным током рабочих сред (водяные экономайзеры, воздухоподогреватели и аппараты с кипящим теплоносителем) определяется выражением
∆t = ε∆t ∆tσ − ∆tм |
, |
(1.52) |
||
ln |
∆tσ |
|
|
|
∆tм |
|
|
||
|
|
|
где ∆tб , ∆tм - соответственно наибольшее и наименьшее значение температурного напора между средами на входе и на выходе из Т.А., 0С;
ε∆t - поправочный коэффициент, определяемый в зависимости от параметров Р и R с помощью номограмм рис. П.2.1 – П.2.5 в зависимости от схемы движения теплоносителей.
.
P= t2′′′ −t′2′ , t1 −t2
R = |
t1′ −t1′′ |
|
|
. |
|
t2′′ −t2′ |
(1.53)
(1.54)
Значения ∆tб и ∆tм могут быть выбраны из расчетных величин ∆tвх и ∆tвых, как большее и меньшее из них, определяемых:
а) для прямоточной схемы (водяные экономайзеры и воздухоподогрева-
тели):
∆tвх = t1′ −t2′ , |
∆tвых = t1′′−t2′′; |
(1.55) |
б) для противоточной схемы:
∆tвх = t1′ −t2′′ , |
∆tвых = t1′′−t2′ ; |
(1.56) |
17
в) для теплообменников с кипящим носителем: |
|
|
∆tвх = t1′′−ts , |
∆tвых = t1′ −ts ; |
(1.57) |
г) для теплообменников типа « отопительный прибор»: |
|
|
∆t = t1ж −t2ж = t1ж −tв . |
(1.58) |
1.6. ОПРЕДЕЛЕНИЕ РАСЧЕТНОЙ ТЕПЛООБМЕННОЙ ПОВЕРХНОСТИ
Расчетная теплообменная поверхность теплообменного аппарата Нр , м2, определяется из уравнения теплопередачи (1.2) при известных значениях Q , К и ∆t :
Нр = Q/(К ∆t) . |
(1.59) |
Для теплообменников типа « отопительный прибор» расчетная поверхность отопительного прибора определяется по наружной поверхности трубы без учета оребрения, так как оребрение учтено формулой (1.51) в величине коэффициента теплопередачи данного типа теплообменника.
18
2. ЗАДАНИЯ
ЗАДАНИЕ № 1
Определить расчетную теплообменную поверхность водяного экономайзера для исходных данных, приведенных в табл. 2.1.
|
|
|
|
|
|
|
Таблица 2.1 |
|
Параметр, |
|
|
Номер варианта |
|
||||
размерность |
1 |
2 |
3 |
4 |
5 |
|||
% |
|
|
CO2, % |
8,62 |
7,98 |
8,14 |
9,07 |
8,84 |
Состав газовдымовых, |
|
|
H2O, % |
18,3 |
18,45 |
18,7 |
18,15 |
18,84 |
|
|
|
||||||
|
|
|
O2 , % |
1,72 |
1,88 |
1,91 |
1,83 |
2,02 |
|
|
|
N2 , % |
71,36 |
71,69 |
71,25 |
70,95 |
71,21 |
|
|
|
|
2,7 |
4,4 |
7,1 |
10,0 |
6,2 |
q5 |
, |
|
% |
|||||
t’1 |
, |
°C |
350 |
450 |
410 |
380 |
300 |
|
t’’1 |
, |
°C |
200 |
180 |
160 |
190 |
120 |
|
t’2 |
, |
°C |
30 |
45 |
65 |
80 |
20 |
|
V1 , нм3/ч |
10000 |
8500 |
20000 |
6800 |
15000 |
|||
G2 |
, |
т/ч |
4,0 |
6,0 |
25,0 |
4,7 |
7,1 |
|
, (м2К)/Вт |
0,015 |
0,021 |
0,007 |
0,01 |
0,005 |
|||
DH × S, мм |
32,9 × 2,8 |
41,8 × 2,8 |
41,8 × 2,8 |
26,8 × 2,3 |
32,9 × 2,8 |
|||
W1ж , м/с |
10,0 |
9,3 |
8,4 |
8,0 |
8,8 |
|||
W2ж , м/с |
0,5 |
0,68 |
0,39 |
1,1 |
1,27 |
|||
σ1 |
|
|
|
1,5 |
1,5 |
1,5 |
1,5 |
1,5 |
σ2 |
|
|
|
2,0 |
2,1 |
2,2 |
2,4 |
2,5 |
Схема течения |
противоток |
прямоток |
противоток |
противоток |
прямоток |
|||
Расположение |
шахмат. |
коридор. |
шахмат. |
коридор. |
шахмат. |
|||
труб |
|
|
|
|
|
|
|
|
19
20
Продолжение табл.2.1
Параметр, |
|
|
|
Номер варианта |
|
|
|
||||
размерность |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
|||
, % |
CO2, % |
8,57 |
8,33 |
9,12 |
8,44 |
8,27 |
8,98 |
8,75 |
8,68 |
||
H2O, % |
|
|
|
|
|
|
|
|
|||
Состав дымовыхгазов |
19,01 |
18,69 |
17,98 |
18,85 |
18,6 |
18,44 |
18,91 |
18,55 |
|||
O2 , % |
1,68 |
1,77 |
1,74 |
1,81 |
1,89 |
1,93 |
1,87 |
1,78 |
|||
|
|||||||||||
|
N2 , % |
70,74 |
71,21 |
71,16 |
70,9 |
71,4 |
70,65 |
70,47 |
70,99 |
||
|
|
|
3,5 |
5,8 |
9,3 |
8,1 |
8,4 |
7,0 |
6,0 |
5,0 |
|
q5 |
, |
% |
|||||||||
t’1 |
, |
°C |
500 |
425 |
390 |
420 |
440 |
475 |
490 |
400 |
|
t’’1 |
, |
°C |
200 |
170 |
130 |
170 |
185 |
190 |
180 |
200 |
|
t’2 |
, |
°C |
25 |
35 |
60 |
40 |
50 |
55 |
30 |
25 |
|
V1 , нм3/ч |
4500 |
7900 |
11700 |
12900 |
5700 |
4000 |
20400 |
19200 |
|||
G2 |
, |
т/ч |
3,5 |
6,2 |
12,8 |
9,0 |
8,0 |
3,4 |
15,0 |
13,0 |
|
, (м2 К)/Вт |
0,035 |
0,03 |
0,009 |
0,012 |
0,016 |
0,02 |
0,025 |
0,03 |
|||
DH × S, мм |
41,8×2,8 |
32,9×2,8 |
41,8×2,8 |
26,8×2,3 |
32,9×2,8 |
32,9×2,8 |
41,8×2,8 |
32,9×2,8 |
|||
W1ж , м/с |
9,6 |
9,2 |
9,5 |
8,2 |
8,4 |
8,6 |
9,1 |
9,0 |
|||
W2ж , м/с |
1,5 |
0,84 |
1,3 |
0,3 |
0,6 |
0,7 |
0,75 |
0,8 |
|||
σ1 |
|
|
1,5 |
1,5 |
1,5 |
1,5 |
1,5 |
1,5 |
1,5 |
1,5 |
|
σ2 |
|
|
2,3 |
2,0 |
2,3 |
2,5 |
2,4 |
2,0 |
2,2 |
2,1 |
|
Схема течения |
противот. |
противот. |
пртивот. |
прямот. |
противот. |
противот. |
прямот. |
противот. |
|||
Расположение |
шахмат. |
коридор. |
шахмат. |
коридор. |
шахмат. |
коридор. |
шахмат. |
шахмат. |
|||
труб |
|
|
|
|
|
|
|
|
|
|
20