краткий курс тепломассообмена
.pdfЦилиндр конечных размеров
Аналогично находятся температуры и для других тел конечных размеров путем перемножения температур их образующих.
Такие задачи часто возникают в печах, где происходит конвективный нагрев заготовок. Из-за не равномерности их нагрева, внутри заготовок возникают термические напряжения, которые могут привести к повреждению заготовок.
Регулярный тепловой режим
В стадии регулярного режима начальные условия начинают играть второстепенную роль, и изменение избыточной температуры любой точки объема подчиняется экспоненциальному закону
J = AU |
×e−μ12Fo , |
|
1 |
1 |
|
откуда
ln ϑ = −μ2Fo + ln AU |
= −μ2 |
a |
τ + const |
||
R2 |
|||||
1 |
1 1 |
1 |
|
или, сокращенно,
ln ϑ = −mτ + const
41
Тепловой режим, при котором натуральный логарифм избыточной температуры изменяется по закону прямой линии, называется регулярным тепловым режимом. Множитель m характеризует скорость протекания явления и носит название темпа процесса
|
ln J - ln J |
|
|
ln |
T1 |
− Tж |
|||||
m = |
|
= |
|
T - T |
|||||||
1 |
|
2 |
|
|
|
2 |
ж |
. |
|||
- t |
t |
|
|
|
|||||||
|
t |
2 |
|
|
2 |
- t |
|||||
|
|
1 |
|
|
|
|
|
1 |
|
График изменения логарифма избыточной температурыот времени.
Для неограниченной пластиныприBi → ∞ темп охлаждения прямо пропорционален коэффициенту температуропроводности:
a =ψ × m.
Здесь y-размерный коэффициент формы.
æ |
R |
ö2 |
|
|
1 |
|
|
|
yн.пл. = ç |
÷ |
= |
|
|
|
. |
||
|
|
p |
ö |
2 |
||||
è m1 |
ø |
æ |
|
|
||||
|
|
|
ç |
|
÷ |
|
|
|
|
|
|
|
|
|
|||
|
|
|
è |
2R ø |
|
|
Аналогично определяетсяразмерный коэффициент и для тел иной формы:
yн.цил. |
= |
|
1 |
|
|
; |
|
|
|
p |
ö |
2 |
|||
|
æ |
|
|
|
|||
|
ç |
|
|
÷ |
|
|
|
|
|
|
|
|
|||
|
è |
1,31R.ø |
|
|
yшара = |
1 |
|
|
; |
||
æ p |
ö |
2 |
||||
|
ç |
|
|
÷ |
|
|
|
|
|
|
|
||
|
è R |
ø |
|
|
42
КОНВЕКТИВНЫЙ ТЕПЛООБМЕН
Понятие конвективного теплообмена охватывает процесс теплообмена при движении жидкости или газа. При этом процесс переноса теплоты осуществляется одновременно конвекцией и теплопроводностью.
Конвекция. Конвекция возможна только в текучей среде, где перенос теплоты неразрывно связан с переносом самой среды.
Конвективный теплообмен между потоком жидкости или газа и поверхностью соприкасающегося с ним тела называется теплоотдачей. При ее расчетах используют закон Ньютона–Рихмана:
dQ = α (tc − tж ) dF .
В этом уравнении наибольшей трудностью является определение коэффициента теплоотдачи α, который зависит от множества факторов и часто определяется опытным путем. Для расчета коэффициента теплоотдачи используют критерии теплового подобия. Тогда, коэффициент теплоотдачи можно найти используя число Нуссельта, которое характеризует отношение термического сопротивления теплопроводности слоя жидкости толщиной l к термическому сопротивлению теплоотдачи:
Nu = α ×l .
λж
Число Нуссельта зависит от геометрии тела, условий теплообмена, свойств и характера движения среды (ламинарный, турбулентный). Различают также свободную и вынужденную конвекцию.
Свободная конвекция
Явление конвекции можно объяснить законом Архимеда и явлением теплового расширения тел. При повышении температуры объем жидкости возрастает, а плотность уменьшается. Под действием архимедовых сил менее плотная нагретая жидкость поднимается вверх, а более плотная холодная жидкость опускается вниз.
С явлением конвекции связаны процессы глобальной циркуляции атмосферных масс воздуха. Все ветры вызваны конвекционными потоками, возникающими из-за того, что большая часть энергии Солнца попадает на Землю вблизи экватора. Когда воздух нагревается, он расширяется и поднимается, а взамен к экватору устремляется поток более холодного и плотного воздуха. Так образуется ветер.
Движение среды в рассматриваемом объеме возникает за счет неоднородности распределения температуры, и как следствие, неоднородности распределения плотности. Такое движение называется гравитационным
43
движением, а конвекция свободной. Данная задача распространена для отопительных установок, так как в помещении движение воздуха создается в основном за счет разности температур воздуха возле отопительной установки и вдали от нее.
В данном случае число Нуссельта будет зависеть от числа Грасгофа и числа Прандтля:
Nu = f (Gr,Pr) .
Число Грасгофа характеризует отношение гравитационныхсил к силам вязкости:
Gr = g ×β (tc -tж )l3 ,
νж2
g – ускорение свободного падения, м/с2, νж – кинематическая вязкость среды, м2/с,
β – коэффициент объемного расширения (для жидкости берется по таблице, для газов рассчитывается по формуле β = T1 , Т – температура газов, К), 1/К.
Число Прандтляхарактеризует отношение скорости распространения количества движения и скорости распространения температуры:
Pr = νж .
aж
Число Нуссельта также зависит от геометрии тел, тогда:
вдоль вертикальной пластины
при 103 < Gr ×Pr <109 (ламинарный режим течения) |
|||||
Nulж = 0.6×(Grlж ×Prж ) |
0.25 |
æ |
Pr ö0.25 |
||
|
ç |
ж |
÷ |
, |
|
|
|
||||
|
|
è |
Prc ø |
|
|
|
|
|
|
|
0.25 |
æ |
Pr |
ö0.25 |
|
|
||
|
|
Nulж = 0.75×(Grlж ×Prж ) |
, |
|
|||||||||
|
|
|
ç |
ж |
÷ |
|
|||||||
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
è |
Prc |
ø |
|
|
||
при 6×1010 < Gr ×Pr |
(турбулентный режим течения) |
|
|||||||||||
|
|
|
|
|
|
|
|
0.33 æ |
Pr |
ö0.25 |
|||
|
Nulж = Nulж = 0.15×(Grlж ×Prж ) |
||||||||||||
|
|
ç |
|
ж |
÷ |
, |
|||||||
|
|
|
Prc |
||||||||||
|
|
|
|
|
|
|
|
|
è |
|
ø |
|
44
около горизонтальных труб
|
|
0.25 |
æ |
Pr |
ö0.25 |
||
Nudж = 0.5×(Grdж ×Prж ) |
|||||||
|
ç |
ж |
÷ |
, |
|||
|
Prc |
||||||
|
|
|
è |
ø |
|
æ |
ö0.25 |
||
ç |
Prж |
÷ |
– поправка Михеева, для газов принимаемая равной единице, так как |
|
|||
è |
Prc ø |
|
число Прандтля для газов слабо зависит от температуры.
при Gr ×Pr < 200 Nudж =1.18×(Grdж ×Prж )0.125 ,
в этой формуле определяющей температурой является нетемпература жидкости, а t = 0.5(tc + tж ) .
Вынужденная конвекция
Движение среды в рассматриваемом объеме происходит под действием внешних сил, (например за счет насоса или вентилятора).
В данном случае число Нуссельта будет зависеть от числа Рейнольдса и числа Прандтля:
Nu = f (Re,Pr) .
Число Рейнольдса определяет соотношение между силой инерции и силой внутреннего трения:
Re = ω ×l .
ν ж
ω – скорость движения потока, м/с.
|
|
Для продольно омываемой плоской поверхности |
|||||||
при |
Re <105 |
(ламинарный режим течения) |
|
||||||
|
|
|
|
|
æ |
ö0.25 |
|
||
|
|
|
Nulж = 0.33×Rel0.5ж Prж0.33 |
ç |
Prж |
÷ |
, |
||
|
|
|
|||||||
|
|
|
|
|
è |
Prc ø |
|
||
|
|
|
|
|
æ |
ö0.25 |
|
||
|
|
|
|
|
|
||||
|
|
|
Nulж = 0.66×Rel0.5ж Prж0.33 |
ç |
Prж |
÷ |
, |
||
|
|
|
|
||||||
|
|
|
|
|
è |
Prc ø |
|
||
при |
105 < Re |
(турбулентный режим течения) |
|
45
|
|
|
|
|
|
|
æ |
ö0.25 |
|
||||
|
Nulж = 0.0296×Rel0.8ж Prж0.43 ç |
|
Prж |
÷ |
, |
||||||||
|
|
|
|
||||||||||
|
|
|
|
|
|
|
è |
|
Prc ø |
|
|||
|
|
|
|
|
|
|
æ |
ö0.25 |
|
||||
|
|
|
|
|
|
|
|
||||||
|
|
Nulж = 0.037 ×Rel0.8ж Prж0.43 |
ç |
Prж |
÷ . |
||||||||
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
è |
|
Prc ø |
|
|||
Для движенияжидкости внутри трубопровода |
|||||||||||||
при Re < 2300 Gr ×Pr < 8×105 |
(вязкостный режим течения) |
||||||||||||
|
|
|
|
|
æ |
d ö0.33 |
æ |
|
μc ö−0.14 |
, |
|||
|
|
|
|
|
|
||||||||
|
Nudж =1.55×ç Redж |
÷ |
ç |
|
|
÷ |
|||||||
|
|
||||||||||||
|
|
|
|
|
è |
l ø |
è |
|
μж ø |
|
при |
Re < 2300 |
8×105 < Gr ×Pr |
(гравитационно-вязкостный режим течения) |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.1 æ |
Pr |
ö0.25 |
||||
|
|
|
Nudж = 0.15×Re0.33dж Prж0.33 (Grdж ×Prж ) |
||||||||||||||||||||
|
|
|
|
|
|
ç |
ж |
÷ |
, |
||||||||||||||
|
|
|
|
|
|
Prc |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
è |
ø |
|
|
при |
104 < Re |
(турбулентный режим течения) |
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
æ |
|
Prж |
|
ö0.25 |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
Nudж |
= 0.021×Red0.8ж Prж0.43 |
ç |
|
|
÷ . |
|
|
|||||||||||
|
|
|
|
|
|
Prc |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
è |
|
ø |
|
|
|
||||||
|
|
При поперечномомывании одиночного цилиндра |
|||||||||||||||||||||
при |
Re <1000 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
æ |
|
Prж |
|
ö0.25 |
, |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
Nudж = 0.5×Re0.5dж Prж0.38 ç |
|
|
÷ |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
Prc |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
è |
|
ø |
|
|
|
|
||||||
при |
103 < Re < 2×105 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
æ |
Prж |
ö0.25 |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
Nudж |
= 0.25×Red0.6ж Prж0.38 ç |
|
÷ |
, |
|
|
|||||||||||
|
|
|
|
|
|
Prc |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
è |
ø |
|
|
|
||||||||
при |
3×105 < Re < 2×106 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
æ |
Prж |
ö0.25 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
Nudж |
= 0.023×Re0.8dж Prж0.38 |
ç |
÷ . |
|
|
|||||||||||||
|
|
|
|
|
Prc |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
è |
ø |
|
|
|
46
При поперечномомывании пучка труб
Различают два вида расположения труб в пучке: коридорное и шахматное. Расстояние между двумя соседними центрами труб поперек движения среды S1, вдоль – S2, тогда
Коридорное |
|
|
|
|
Шахматное |
||
расположение труб. |
|
|
|
|
расположение труб. |
||
для коридорного расположения труб |
|
|
|
|
|
||
|
|
|
|
æ |
ö0.25 |
||
|
|
|
|
||||
|
Nudж = 0.26×Red0.65ж |
Prж0.33 |
ç |
Prж |
÷ |
εsεi , |
|
|
|
||||||
|
|
|
|
è |
Prc ø |
|
−0.15
ε = æ S2 ö – поправочный коэффициент, учитывающий влияние s ç ÷
è d ø
относительных шагов для глубинных рядов пучка.
для шахматного расположения труб
|
|
|
|
|
|
|
|
|
|
æ |
ö0.25 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
Nudж = 0.41×Re0.6dж Prж0.33 |
ç |
Prж |
÷ |
εsεi , |
|
||||
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
è |
Prc ø |
|
|
|
|
|
S1 |
|
|
æ |
S1 |
ö |
1 |
|
|
|
|
|
S1 |
|
|
при |
|
εs |
6 |
, |
|
при |
|
εs =1.12 . |
|||||||
< 2 |
= ç |
÷ |
|
|
³ 2 |
||||||||||
|
|
|
|||||||||||||
S2 |
S2 |
|
|
S2 |
|||||||||||
|
|
|
è |
ø |
|
|
|
|
|
|
|
|
Поправочный коэффициент εi учитывает изменения теплоотдачи в начальных рядах труб (при не высокой степени турбулентности набегающего потока для первого ряда труб можно принять равным 0.6).
47
Теплообмен при конденсации
Конденсация представляе т собой процесс перехода пара (газа) в жидкое состояние.Конденсация пара часто встречается на практике в ко нденсаторах паровых турбин, опреснительн ых установках и многих других теп лообменных аппаратах.
Различают пленочную ко нденсацию, когда жидкая конденсированная фаза образуется на поверхности теплообмена в виде устойчивой пленки и смачивает всю поверхность, а также капел ьную конденсацию, когда конденсат не смачивает всю поверхность, а происходит образование капель на п оверхности теплообмена. Искусственно капельную конденсацию можно получить, если смазать поверхность маслом ил и воском.
Капельная конденсация на покрашенной поверхности
При пленочной конденсации неподвижного пара при лам инарном течении пленки
на вертикальной стенке
4 |
|
|
rρ2 gλ3 |
|
|
|||
α = 3 4 4μж (tн − tc )h , |
|
|||||||
|
|
|
ж ж |
|
||||
μж =νж ρж – динамическая вязкость, кг/м∙с. |
|
|||||||
на горизонтальной трубе |
|
|
|
|
|
|||
|
|
|
|
|
|
|||
α = 0.7284 |
rρж2 gλж3 |
|
|
, |
||||
μж (tн − tc )d |
||||||||
|
|
|
|
48
Теплообмен при кипении
Кипением называется процесс интенсивного парообразования, происходящего во всем объеме жидкости, находящейся при температуре насыщения или несколько перегретой относительно температуры насыщения, с образованием паровых пузырей.
Различают два основных режима кипения: пузырьковый –кипение, при котором пар образуется в виде отдельных периодически зарождающихся, растущих и отрывающихся паровых пузырей, и пленочный – характеризуется наличием на поверхности пленки пара, обволакивающую эту поверхность и отделяющей ее от жидкости.
Пузырьковое кипение
в условиях свободного движения
Для определенного рода жидкостей (например, для воды) коэффициент теплоотдачи при развитом кипении зависит лишь от тепловой нагрузки и давления насыщения
α= 3×q0.7 p0.15 = 38.7Dt2.33 p0.5 ,
вкоторые q и p следует подставлять в Вт/м2 и бар соответственно.
Пленочное кипение при ламинарном движении вдоль вертикальной стенки
|
|
|
|
|
|
α = 0.9434 |
rρп (ρж - ρп ) gλп3 |
, |
|||
μп (tп - tc )h |
|||||
|
|
|
|
при ламинарном движении на наружной поверхности горизонтального цилиндра
|
|
|
|
|
|
α = 0.7284 |
rρп (ρж - ρп ) gλп3 |
, |
|||
μп (tп - tc )d |
|||||
|
|
|
|
при турбулентном движении
|
|
|
|
|
α = 0.254 |
(ρж - ρп ) g ×λж2 срп |
. |
||
|
||||
|
|
νж2 |
49
МАССООТДАЧА
Вдвижущейся однокомпонентной среде теплота переноситься теплопроводностью и конвекцией. Этот процесс называется конвективным теплообменом. По аналогии перенос вещества в многокомпонентной среде совместно происходящими процессами молекулярной диффузии и конвекции называют конвективным массообменном.
Воднородной по температуре и давлению макроскопически неподвижной двухкомпонентной смеси плотность потока массы одного из компонентов за счет молекулярной диффузии определяется законом Фика
ji = -D ¶¶ρni .
j – плотность потока массы, кг/(м2К); D – коэффициент диффузии.
Аналогично теплоотдаче конвективный массообмен между жидкой или твердой поверхностью и окружающей средой называют массоотдачей.
Для расчетов массоотдачи используется уравнение:
jc = β (ρc − ρ0 ),
β – коэффициент массоотдачи, м/с; ρс, ρ0 – концентрации диффузионного вещества на поверхности раздела фаз и вдали от нее соответственно, кг/м3.
Исходя из аналогии процессов теплообмена и массообмена, можно записать:
Nuдиф = f (Re,Prдиф ).
Prдиф = νD – диффузионное число Прандтля;
Тогда, коэффициент массоотдачи можно найти, используя число Нуссельта диффузионное:
Nuдиф = βD×l .
50