Курсовое
.pdf2. Тепловой расчет теплообменного аппарата
Определим количество теплоты, отводимое от воздуха:
Q = W1(t1′ − t1′′) = W2(t′′2 − t′2) = 30 · 1,015 · (107,8 − 37,6) = 2138 кВт
Для упрощения расчетов теплофизические характеристики воздуха и
воды возьмем при t = 100 °С и t = 20 °С соответственно.
Уравнение теплопередачи имеет вид:
Q = kFΔtлог,
где k – коэффициент теплопередачи; F – площадь поверхности
теплообмена; tлог – среднелогарифмический температурный напор.
Дальнейший расчет будет направлен на нахождение компонентов вышеупомянутой формулы и нахождению из нее площади поверхности
теплообмена. |
|
|
|
|
|
||
|
Произведем расчет температурного напора |
tлог: |
|||||
|
|
Δtб − Δtм |
|
77,8 |
− 17,6 |
||
|
Δtлог = |
|
|
= |
|
|
= 40,5 °С, |
|
ln |
Δtб |
ln |
77,8 |
|||
|
|
|
|
||||
|
|
Δtм |
|
17,6 |
|
||
|
|
|
|
|
|
||
где |
Δtб = t1′ − t2′′ = 107,8 |
− 30 = 77,8 °С, |
Δtм = t1′′ − t2′ = 37,6 − 20 = |
||||
17,6 °С. Основанием для выбора данных формул является выбор противоточной схемы движения теплоносителей. Она является более эффективной при прочих равных условиях. Также по расчетам лог для противоточной схемы оказалась больше, чем для прямоточной (при прямотоке
лог = 32,8 °С), что также указывает на использование противоточной схемы.
11
|
120 |
|
|
100 |
|
|
80 |
|
°С |
60 |
|
t, |
|
|
|
|
|
|
40 |
|
|
20 |
|
|
0 |
F |
|
|
|
|
Воздух |
Вода |
Рисунок 6 – Зависимость температуры теплоносителей от теплоемкости массовых |
||
|
расходов для прямоточной схемы |
|
|
120 |
|
|
100 |
|
|
80 |
|
°С |
60 |
|
t, |
|
|
|
|
|
|
40 |
|
|
20 |
|
|
0 |
F |
|
|
|
|
Вода |
Воздух |
Рисунок 7 – Зависимость температуры теплоносителей от теплоемкости массовых расходов для противоточной схемы
На рисунках 7 и 8 представлены зависимости для полученных выше
данных.
12
Коэффициент теплопередачи k определяется по формуле:
k = [( |
1 |
+ Rз) |
d |
+ |
1 |
]−1, |
|
|
|
α ɳφ |
|||||
|
α |
d |
в |
|
|||
|
2 |
|
|
1 |
|
||
где α1, α2 – коэффициент теплоотдачи воздуха и воды |
|||||||
соответственно, Вт/(м2 · К); |
Rзагр – |
|
термическое сопротивление слоя |
||||
загрязнения, (м2 · К)/Вт; d, dв – наружный и внутренний диаметр трубы соответственно, мм; ɳ – коэффициент эффективности оребрения; φ –
коэффициент оребрения.
Внутренний диаметр трубы:
dв = d − 2δ = 19,6 − 2 · 1,5 = 16,6 мм.
Примем начальное приближение α2 = 6 кВт/(м2 · К).
Найдем значение термического сопротивления слоя загрязнения.
Примем время работы теплообменного аппарата τ=5000 часов. Тогда, исходя из значений, снятых с графика на рисунке 9, термическое сопротивление будет равно:
Rзагр = 0,0005 (м2 · К)/Вт.
Рисунок 8 – Изменение во времени термического сопротивления слоя загрязнения
13
Для дальнейшего расчета предварительно найдем ряд необходимых величин, таких как: площадь поверхности ребер Fр; площадь поверхности межреберных промежутков Fс; площадь полной поверхности оребренной трубы Fрс; площадь поверхности гладкой трубы, несущей оребрение F1.
Используем формулы:
F |
|
= |
π(D2 − d2) |
+ πDδ |
= |
π · (37,22 − 19,62) |
+ π · 37,2 · 0,5 = 1628,7 мм2; |
||||||
р |
|
|
|||||||||||
|
2 |
|
|
|
|
|
1 |
|
2 |
|
|||
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
F |
с |
= |
πd(S |
р |
− δ |
) = π · 19,6 · (2,8 − 0,1) = 166,3 мм2; |
||||
|
|
|
|
|
|
|
|
2 |
|
|
|
||
Fрс = Fр + Fс = 1795 мм2;
1 = π р = · 19,6 · 2,8 = 172,4 мм2
В расчетах принимаются значения δ1 = 0,5 мм и δ2 = 0,1 мм исходя из
среднего значения толщины ребра δр = 0,3 мм. То есть, δр = δ1+δ2 .
2
Коэффициент оребрения рассчитывается по формуле:
|
|
|
|
φ = |
рс |
= |
1795 |
|
= 10,41. |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
172,4 |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Найдем |
коэффициент |
теплоотдачи |
воздуха α1 исходя из формулы |
|||||||||||||||||||||
α = |
Nu·λ1 |
, где Nu – число Нуссельта, λ |
– коэффициент теплопроводности |
||||||||||||||||||||||
|
|||||||||||||||||||||||||
1 |
l0 |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
воздуха, l0 – линейный размер. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
Линейный размер рассчитывается по формуле: |
|
|||||||||||||||||||||||
|
|
|
|
|
d Fc |
|
|
|
|
|
|
|
|
|
|
|
|
Fp |
|
|
|||||
|
|
|
l0 = |
+ √ |
π |
(D2 − d2) |
|
|
|||||||||||||||||
|
|
|
|
Fpc |
4 |
Fpc |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
19,6 166,3 |
|
|
|
|
|
|
1628,7 |
|
|||||||||||||||
|
l0 = |
+ |
√ |
π |
(37,22 − 19,62) |
= 27,2 мм. |
|||||||||||||||||||
|
1795 |
|
|
1795 |
|
||||||||||||||||||||
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
14
Критерий Нуссельта:
Nu1 = 0,36 ∙ Re1n ∙ Pr10,33 ∙ Cz ∙ Csφ−0,5,
Где Cz и – поправочные коэффициенты на число рядов в пучке и коэффициент формы пучка соответственно. При числе рядов больше четырех поправочный коэффициент принимаем Cz = 1.
Коэффициент формы пучка рассчитываем по формуле:
|
σ |
− 1 |
|
0,1 |
3,31 − 1 |
|
0,1 |
= ( |
1 |
|
) |
= ( |
|
) |
= 1,092, |
|
|
|
|||||
|
σ′ |
− 1 |
|
|
1,95 − 1 |
|
|
|
2 |
|
|
|
|
|
|
где σ1 = 1 = 64,819,6 = 3,31;
σ2 = 2 = 20,319,6 = 1,036;
σ′ |
= √ |
σ12 |
+ σ2 = √ |
3,492 |
+ 1,0362 = 1,95. |
|
||||||||||||
|
|
|
|
|||||||||||||||
2 |
|
|
|
|
4 |
|
2 |
|
4 |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Число Рейнольдса: |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
c l |
0 |
|
21 27,2 10−3 |
||
|
|
|
|
|
|
|
|
Re = |
1 |
= |
|
|
= 59749, |
|||||
|
|
|
|
|
|
|
|
|
|
ν |
|
9,56 10−6 |
||||||
Где ρ = |
|
p1′ |
= |
|
|
250000 |
|
= 2,29 кг/м3; |
||||||||||
|
′ |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
RT1 |
|
287∙(107,8+273,15) |
||||||||||
ν = |
|
μ |
= |
21,9∙10−6 |
= 9,56 ∙ 10−6 м2/с . |
|||||||||||||
|
|
|
||||||||||||||||
|
|
ρ |
|
|
2,29 |
|
|
|
|
|
|
|
|
|
|
|||
Показатель степени n:
n = 0,6 ∙ φ0,07 = 0,6 10,410,07 = 0,707.
Подставим значения в формулу для нахождения числа Нуссельта:
Nu1 = 0,38 597490,707 0,6880,33 1 1,092 10,41−0,5 = 270,7.
15
Отсюда коэффициент теплоотдачи составляет
α = |
Nuλ1 |
= |
270,7 0,0321 |
= 319,5 Вт/(м2 · К). |
|
|
|||
1 |
l0 |
|
27,2 10−3 |
|
|
|
|||
Определим коэффициент эффективности оребрения ɳ. Для этого
воспользуемся следующей формулой:
η = 1 − (1 − ψ ∙ ∙ ξ) ( ),
где ψ – поправочный коэффициент; – коэффициент эффективности ребра; ξ – поправочный коэффициент для
трапециевидных ребер. |
|
|
|
|
|
|
|
|
||
Формула для вычисления поправочного коэффициента: |
|
|||||||||
|
|
|
ψ = 1 − 0,058m′h |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где m′ = √ |
2αр |
= √ |
2·319,5 |
= 74,5; h = h |
|
+ |
δр |
= 8,8 + |
0,3 |
= 8,95 мм. |
|
0,3·10−3·384 |
р |
|
|
||||||
|
δрλр |
|
|
2 |
|
2 |
|
|||
Тогда:
ψ = 1 − 0,058 74,5 8,95 10−3 = 0,96.
Коэффициент эффективности круглого ребра постоянного сечения и поправочный коэффициент для трапециевидных ребер определим исходя из найденного значения m′h = 0,67 по графику (рисунок 9).
16
Рисунок 9 – К определению коэффициента эффективности круглых ребер постоянного сечения (а) и поправочного коэффициента для трапециевидных ребер (б)
При |
D |
= |
37,2 |
= 1,9: =0,85; при √ |
δ2 |
= √ |
0,1 |
= 0,45: ξ=1,05. |
d |
19,6 |
δ |
|
|||||
|
|
|
0,5 |
|
||||
|
|
|
|
1 |
|
|
|
|
Подставим полученные значения в формулу коэффициента
эффективности оребрения:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1628,7 |
|
||||
η = 1 − (1 |
− ψ ∙ ∙ ξ) ( |
|
|
|
) = 1 − (1 − 0,96 0,85 1,05) |
( |
|
|
|
) = 0,87 |
|||||||||||
|
1795 |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
После нахождения всех необходимых величин, рассчитываем |
|||||||||||||||||||||
коэффициент теплопередачи: |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
k = [( |
|
1 |
+ 0,0005) |
∙ |
19,6 |
+ |
1 |
|
|
]−1 |
= 882,8 |
|
Вт |
. |
|||||||
|
|
|
|
319,5 0,87 10,41 |
|
|
|||||||||||||||
6000 |
|
|
|
|
16,6 |
|
|
|
|
|
м · К |
||||||||||
Поверхность теплообмена: |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
2138 103 |
= 63,7 м2. |
|
|
|
|
|
|
|||||
|
|
|
= |
|
|
|
= |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
1 |
|
|
|
лог |
|
882,8 40,5 |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
17
3. Определение характеристик трубного пучка
Расчет площади фронтального сечения производится по формуле:
f = fϰ1,
где f1 – площадь проходного сечения, а ϰ – коэффициент загромождения.
Рассчитаем площадь проходного сечения:
f1 = |
М1 |
= |
30 |
|
= 0,62 м2 . |
|
21 2,29 |
||||
|
c1ρ1 |
|
|||
Коэффициент загромождения для шахматного пучка бывает двух видов:
фронтальный и диагональный. Необходимо рассчитать оба варианта и выбрать наибольший из них.
Расчет фронтального коэффициента загромождения:
ϰф = |
S1 − d |
− |
2hрδр |
= |
68,4 − 19,6 |
− |
2 |
8,8 0,3 |
= 0,69. |
||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
S1 |
|
|
|
S1Sр |
68,4 |
|
68,4 2,8 |
|
|
|
|||||||||
Расчет диагонального коэффициента загромождения: |
|
|
|
||||||||||||||||||
|
2(S′ − d) |
|
|
4h |
δ |
2 (39,8 − 19,6) |
|
4 8,8 0,3 |
|
||||||||||||
ϰд = |
2 |
|
− |
|
|
р р |
= |
|
|
|
|
|
− |
|
|
|
|
= 0,54. |
|||
S1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
S1Sр |
68,4 |
|
|
|
|
68,4 2,8 |
|
||||||||||
Очевидно, что из двух значений коэффициента загромождения выберем фронтальный, так как по условию нам необходимо большее значение.
Подставим полученные значения в формулу площади фронтального сечения:
= 0,620,69 = 0,9 м2
18
Ширина a и высота L пучка в данном теплообменном аппарате одинаковы. Так как a = √f, следовательно, обе величины можно рассчитать по формуле:
a = L = √f = √0,9 = 0,95 м.
Число труб в ряду:
|
n1 = |
a |
= |
0,95 |
= 14. |
|
|
|
|
|
|
||||
|
S |
68,4 ∙ 10−3 |
|
||||
|
|
|
1 |
|
|
|
|
Число рядов в пучке: |
|
|
|
|
|
|
|
z = |
F1 |
= |
|
|
63,7 |
|
= 78. |
|
|
||||||
πdLn |
π 19,6 10−3 0,95 14 |
||||||
|
1 |
|
|
|
|
|
|
Определим глубину пучка:
b = z ∙ S2 = 78 20,3 10−3 = 1,6 м.
Объем теплообменного аппарата определяется формулой:
V = a2b = 0,952 1,6 = 1,44 м3
Обращаясь к изначально заданным параметрам теплообменного аппарата, проведем проверку полученного значения объема:
|V − V11| |
100% = |
|1,44 − |
1,5| |
100% = 4 % |
|
|
|
||
V11 |
1,5 |
|
|
|
Полученное значение объема меньше заданного в варианте на 4%.
Рассчитаем массу теплообменного аппарата. Для этого воспользуемся формулой:
M = ρMVM.
Объем металла рассчитывается по формуле:
|
π(d − dв)2 |
|
π(D − d)2 |
|
|
L |
|
V = ( |
|
∙ L + ( |
|
∙ δ |
) ∙ |
|
) ∙ n ∙ z |
|
|
|
|||||
M |
4 |
|
4 |
р |
|
1 |
|
|
|
|
|
Sp |
|||
19
|
|
π(19,6 ∙ 10−3 |
− 16,6 ∙ 10−3)2 |
|
|
|
|||
VM = |
( |
|
|
|
|
∙ 0,95 |
|
|
|
|
|
4 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
+ ( |
π(37,2 ∙ 10−3 |
− 19,6 ∙ 10−3)2 |
∙ 0,3 ∙ 10−3) ∙ |
0,95 |
) ∙ 14 ∙ 78 |
|||
|
|
|
4 |
|
2,8 ∙ 10−3 |
||||
|
|
|
|
|
|
|
|
||
|
|
= 0,097 м3. |
|
|
|
|
|
||
Вернемся к нахождению массы теплообменного аппарата. Так как материал аппарата – алюминий, то плотность берем соответствующую
(ρ=8900 кг/м3). Находим:
M = 8900 0,097 = 863,3 кг.
Проверка массы относительно начальных данных:
|M − M11| |
100% = |
|863,3 − 1470| |
100% = 41% |
|
|
|
|
||
M11 |
1470 |
|
|
|
Полученное значение массы меньше заданного в варианте на 41%.
Далее необходимо определить коэффициент теплоотдачи со стороны воды. Ранее он был принят как начальное приближение, однако на данном этапе необходимо уточнить α2.
20
