2.2.2. Коэффициент теплоотдачи конденсирующегося бензола определяют по зависимости [1].
альфа
б=А/
4√d*∆t
б
где d-гидравлический диаметр канала;
∆t
б – температурный напор от конденсирующихся паров бензола к стенке канала со стороны бензола;
А – коэффициент, зависящий от физических свойств конденсирующегося бензола и от скрытой теплоты его парообразования. Зависимость для определения коэффициента А приведена далее.
Дальнейший
расчёт проводим последовательными приближениями, задаваясь в каждом новом приближении температурой стенки канала со стороны бензола. Каждое приближение заканчиваем сравнением заданного и полученного значений температуры стенки канала со стороны бензола.
Приближение №1.
Принимаем температуру стенки канала со стороны бензола равной tст1 = 57,90C
Тогда, ∆t
б = t
к – t
ст1 = 80,1 – 57,9 = 22,2
0C – температурный напор от конденсирующихся паров бензола к стенке канала, а средняя температура плёнки конденсирующегося бензола равна
t
п = (t
к + t
ст1)/2 = (80,1+57,9)/2 = 69
0C
В соответствии с [1]
А=С
0,75*r
0,25, где С - коэффициент, зависящий от физических параметров конденсата бензола (от температуры его насыщенных паров).
В нашем случае С=3423 [1]. В курсовой работе допускается принимать это значение для всех вариантов заданий.
Тогда,
а
льфа
б=А/
4√d∆t
,б = 3423
0,75*94,5
0,25/
4√0,0196*22,2 = 1395/0,812 = = 1718ккал/м
2часК = 1998Вт/м
2К
Для определения коэффициента теплоотдачи к воде определим её режим течения.
Число Рейнольдса
Re=wdρ/μ=wdγ/μg=(0,525*0,0196*1000)/(0,854*1,02*10
-4*9,81) = 12042
В этой зависимости
μ = 0,854сП=0,854*1,02*10
-4 = 87,108*10
-6кг*с/м
2-динамический коэффициент вязкости воды при её средней температуре t=27,2
0C (таблица №1 или [2]; перевод из единиц “сантипуаз” в “кг*с/м
2 технической системы единиц” см. в [3]).
Режим течения турбулентный.
Для турбулентного режима течения воды коэффициент теплоотдачи определим по зависимости [1].
альфа
в=А
5w
0,8в/d
0,2,
где А
5=1860 (таблица№3, Приложение).
альфа
в=1860*0,525
0,8/0,0196
0,2=2439ккал/м
2часК=2838 Вт/м
2КК
Термические сопротивления загрязнений канала приведены в задании. Термическое сопротивление загрязнения со стороны бензола равно R
загр.б =
=0,0001м
2часК/ккал, а со стороны воды -
R
загр.в = 0,0007м
2часК/ккал.
Определим
термическое сопротивление стальной спирали, принимая её толщину равной δ = 2,5мм, а коэффициент теплопередачи стали равным λ = 40 ккал/м*час*К
R
ст= б/λ = 0,0025/40 = 0,0000625м
2часК/ккал
Тогда
, коэффициент теплопередачи спирального теплообменника равен
к= 1/(1/альфа
б+R
затр.б+R
ст+R
затр.в+1/альфа
в) =
= 1/(1/1718+0,0001+0,0000625+0,0007+1/2439) =
=1/(0,000582+0,0001+0,0000625+0,0007+0,00041)= =539,2ккал/м
2часК=627Вт/м
2К
Необходимая площадь поверхности нагрева спирального теплообменника
F=Q/к∆t= 94500/(539,2*52,9)=3,25м
2
Определяем температуру стенки канала в первом приближении
t
ст1 = t
к – Q/(F*альфа
б) = 80,1 – 94500/(3,25*1718) = 63,18
0C
Так как полученное значение температуры стенки канала заметно отличается от заданного, проводим расчёт во втором приближении.
Приближение №2.
Принимаем температуру стенки канала со стороны бензола равной tст1 = 63,20C
Тогда,
температурный напор на плёнке бензола равен
∆t
б = t
к – t
ст1 = 80,1 – 63,2 = 16,9
0C,
а
средняя температура конденсирующегося бензола равна
t
п = (t
к + t
ст1)/2 = (80,1+63,2)/2 = 71,65
0C
А=С
0,75*r
0,25, где С - коэффициент, зависящий от физических параметров конденсата бензола (от температуры его насыщенных паров). Принимаем С=3423
Тогда,
альфа
б=А/
4√d∆t
,б=3423
0,75*94,5
0,25/
4√0,0196*16,9 =
= 1839ккал/м
2часК = 2140Вт/м
2К
Коэффициент теплоотдачи к воде был определён в приближении №1
альфа
в = 2838 Вт/м
2КК
Термические сопротивления загрязнений канала и стальной спирали принимаем такими же, как и в приближении №1.
Тогда
, коэффициент теплопередачи спирального теплообменника равен
к = 1/(1/альфа
б+R
затр.б+R
ст+R
затр.в+1/ альфа
в) =
= 1/(1/1839+0,0001+0,0000625+0,0007+1/2439) =
=1/(0,000544+0,0001+0,0000625+0,0007+0,00041)= =550,6ккал/м
2часК=641Вт/м
2К
Необходимая площадь нагрева спирального теплообменника
F = Q/(к*∆t) = 94500/(550,6*52,9)=3,24м
2
Определяем температуру стенки канала во втором приближении
t
ст1 = t
к – Q/(F*альфа
б) = 80,1 – 94500/(3,24*1839) = 64,2
0C
Так как полученное значение температуры стенки канала со стороны бензола близко к заданному во втором приближении, то дальнейшие приближения не проводим.
3. По данным таблицы №9 (Приложение) минимальная поверхность серийных спиральных теплообменников составляет 15м
2. Поэтому,
определяем размеры проектируемого теплообменника расчетом.
3.1. Геометрический и компоновочный расчет матрицы спирального теплообменника.
3.1.1. Из рис.5 (Приложение) видно, что
количество витков наружной спирали на 0,5 витка больше, чем количество витков внутренней спирали, т.е.
n
2=n
1+0,5
где:
n
1-количество витков внутренней спирали;
n
2-количество витков наружной спирали.
Длина внутренней спирали
L
1=π(D
1ср/2)*n
1=π((D+d-2t)/2)*n
1,
где:
L
1 – длина внутренней спирали;
D-2t – наружный диаметр внутренней спирали;
D –
наружный диаметр наружной спирали;
D
1ср – средний диаметр внутренней спирали;
t = в + δ = 10 +2,5 = 12,5мм – шаг спирали;
δ – толщина спирали;
d – внутренний диаметр наружной и внутренней спиралей (внутренний диаметр матрицы теплообменника).
Длина наружной спирали
L
2 = π((D+d)/2)*n
2 = π((D+d)/2)*(n
1+0,5))
3.1.2. Определим количество витков внутренней спирали - n1
Общее количество витков обоих спиралей
n
1+n
1+0,5 = (D-d)/(2t);
Откуда,
n
1 = [((D-d)/2t)-0,5]/2 = (D-d-t)/4t
Уравнение поверхности нагрева матрицы
В
к(L
1+L
2) = F;
L
1+L
2 = F/В
к
Здесь:
F – площадь поверхности нагрева (охлаждения) матрицы теплообменника;
В
к, L
1, L
2 – ширина каналов и длина спиралей матрицы.
Подставим в последнее уравнение полученные выше уравнения для длин наружной и внутренней спиралей.
π ((D+d-2t)/2)*n
1+ π((D+d)/2)*(n
1+0,5) = F/В
к;
(D+d-2t)*n
1+(D+d)(n
1+0,5) = 2F/(πВ
к);
(D+d)*n
1-2tn
1+(D+d)*n
1+ (D+d)/2 = 2F/(πВ
к);
2(D+d)*n
1 + (D+d)/2-2tn
1 = 2F/(πВ
к);
2n
1(D+d-t) + (D+d)/2 = 2F/πВ
к;
Подставим в это соотношение полученную ранее зависимость для количества витков внутренней спирали - n
1 = (D-d-t)/4t.
(D+d-t)*(D-d-t)/(2t)+(D+d)/2 = 2F/(πВ
к);
(D-t+d)(D-t-d) + (D+d)*t = 4Ft/(πВ
к);
D
2 - 2Dt+t
2-d
2+Dt+dt = 4Ft/(πВ
к);
D
2-Dt+t
2-d
2+dt-4Ft/(πВ
к) = 0
3.1.3. Выполним расчет, принимая внутренний диаметр матрицы равным d=150мм.
Наружный диаметр наружной спирали определим из последнего полученного уравнения
D
2-12,5D+12,5
2-150
2+150*12,5-(4*3,25*10
6*12,5)/(π500) = 0;
D
2-12,5D+156,25-22500+1875-103450,8 = 0;
D
2-12,5D-123919,55 = 0;
D=358,5мм;
Количество витков внутренней и наружной спиралей
n
1=(D-d-t)/4t=(358,5-150-12,5)/(4*12,5) = 3,92;
n
2=n
1+0,5=4,42
Длина внутренней спирали
L
1=(π(D+d-2t)/2)*n
1=π(358,5+150-25)/2*3,92=2977,2м.
Длина наружной спирали
L
2=(π(D+d)/2)*n
2=(π(358,5+150)/2)*4,42=3530,5м.
3.1.4 Проверка.
Поверхность нагрева матрицы
F=(L
1+L
2)В
к=(2977,2+3530,5)*10
-3*0,5=3,254м
2.
Т.к. необходимое значение поверхности равно 3,25м
2, то расчет проведен верно.
4. Эскизный проект рассчитанного спирального теплообменника привести на чертеже. Для разработки эскиза использовать полученные геометрические параметры спирального теплообменника, а также рис.5 и рис.6.