Курс лекций по ТМО_ Абузова Ф[1].Ф
..pdfЛинейная плотность теплового потока для диаметра d (d1 ≤d ≤d2 ):
ql |
= π k d (tж1 −tж2 ); |
(3.38) |
||||||
q |
|
= |
Q |
= |
π (tж1 −tж2 ) |
. |
(3.39) |
|
l |
l |
Rц |
||||||
|
|
|
|
|
||||
|
|
|
|
|
kΣ |
|
|
Критический диаметр цилиндрической стенки
Термическое сопротивление теплопередаче однослойной цилиндри-
ческой стенки
Rцk |
= |
1 |
|
+ |
|
1 |
ln |
d2 |
+ |
1 |
. |
(3.40) |
|
α1 d1 |
|
|
d1 |
|
|||||||||
|
|
2 |
λ |
|
α2 d2 |
|
|||||||
|
|
|
|
|
|
|
|||||||
|
|
Rα |
|
|
|
R λ |
|
|
Rα |
2 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
Если d2 возрастает, то Rα2 уменьшается, а Rλ – возрастает.
Примем, что α1 , d1, α2 и λ – постоянные величины, и найдем экстре-
мум:
d |
(Rцk )= |
1 |
|
1 |
− |
1 |
=0 . |
(3.41) |
|
|
|
2 |
|||||
dd2 |
2 λ d2 |
α2 d2 |
|
|
Решая (3.41), получаем критический диаметр цилиндрической стенки:
d |
|
=d |
|
= |
2 λ |
. |
(3.42) |
2 |
кр |
|
|||||
|
|
|
α2 |
|
|||
|
|
|
|
|
|
||
Полученное значение критического диаметра соответствует |
мини- |
мальной величине Rцk .
Для цилиндрической стенки существует минимальное значение Rцk и
максимальное значение ql .
33
Rцk
Rц |
Rц |
α2 |
|
|
λ |
Rцα1
Рис. 3.3. Зависимость термических сопротивлений от d2
С учётом слоя изоляции имеем
RцkΣ |
= |
1 |
+ |
1 |
ln |
d2 |
+ |
1 |
ln |
dиз |
+ |
1 |
. |
(3.43) |
|
2 λст |
d1 |
2 λиз |
|
|
|||||||||
|
|
α1 d1 |
|
|
|
d2 |
α2 dиз |
|
Примем, что d1, d2 , λст , α1 , α2 – постоянные величины, а диаметр изоляции dиз изменяется.
Рис. 3.4. Зависимость линейной плотности теплового потока ql от dиз
Из рис. 3.4 видно, что при возрастании dиз от d2 до dкр наложение изоляции увеличивает потери. Линейная плотность теплового потока ql у
голой трубы ( dиз =d2 ) и трубы с dиз = dиз.эф одинакова. При возрастании
34
диаметра изоляции от dкр и дальше значение ql уменьшается. Для рацио-
нальной эффективно работающей изоляции необходимо выполнение ус-
ловия
dкр <d2 . |
(3.44) |
Таким образом, соответствующее значение d2 расположено на рис. 3.4 правее dкр . В случае несоблюдения условия (3.44) выбирают материал
изоляции так, чтобы λ = |
α2 d2 |
. |
|
|
|
|
(3.45) |
|
|
|||||||||||
|
|
|
|
|
|
|
||||||||||||||
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В промышленности λ |
|
|
≤0,2 |
|
Вт |
|
. Если принять,α |
|
= 4 |
Вт |
то |
|||||||||
|
|
|
|
|
|
|
м2 град |
|||||||||||||
|
|
|
|
из |
|
|
|
м град |
2 |
|
|
|||||||||
|
|
|
|
|
|
|
|
2 0,2 |
|
Вт |
|
|
|
|
|
|||||
критический диаметр d |
|
= |
2 |
λ |
= |
м град |
=0,1 м. |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|||||||||||||
кр |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
α2 |
|
|
|
|
|
Вт |
|
|
|
|
|
|||||||
|
|
|
|
4 |
|
|
|
|
|
|||||||||||
|
|
|
|
м2 град |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
Как правило, для большинства диаметров в промышленности усло-
вие (3.44) выполняется и оптимальная изоляция выбирается исходя из минимума приведённых годовых расходов:
Рпр =Э+E К, |
(3.46) |
где Е – нормативный коэффициент сравнительной экономической эф-
фективности;
ΔК – капитальные вложения, [ К]= руб.;
Э – эксплуатационные затраты, [Э]= руб. .
|
|
год |
|
Коэффициент Е находится по формуле |
|
||
E = |
1 |
, |
(3.47) |
|
τсл
где τсл – срок службы, [τсл ]= лет.
Заметим, что для плоской стенки оптимальной изоляции нет.
35
Тепловая изоляция
Вкачестве тепловой изоляции используют волокнистые, порошковые
ипористые материалы, заполненные воздухом. Основа препятствует кон-
векции, воздух создаёт термическое сопротивление.
Коэффициент теплопроводности в значительной степени зависит от
увлажнения, т.к. λ20° |
=0,0296 |
Вт |
, λ20° |
=0,599 |
Вт |
, поры заполня- |
|
м град |
м град |
||||||
возд |
|
вода |
|
|
|||
|
|
|
|
|
ются водой, есть капиллярный эффект.
Тепловую изоляцию рассчитывают по СНиП 2.04.14–88 « Тепловая
изоляция оборудования и трубопроводов», Госстрой, 1989 г.
Особенно низкий коэффициент теплопроводности имеет вакуумно-
−4 |
Вт |
|
|
пористая изоляция при высоких температурах: λ =10 |
|
. |
При этом |
м град |
следует учитывать предельные температуры, при которых изоляция раз-
рушается (по СНиП).
1.3.3. Нестационарная теплопроводность
Физические представления о процессах нагревания и охлаждения тел
Рис. 3.5. Нагрев тела в печи
При нагревании и охлаждении тел температурное поле зависит ещё и от времени – это нестационарные процессы теплопроводности.
Различают 2 группы процессов:
36
∙тело стремится к тепловому равновесию со средой, и на бесконечности эти температуры выравниваются, но на начальных стадиях температу-
ра оси отстает от температуры поверхности;
∙температура периодически изменяется – это процессы с периодиче-
ским нагреванием и охлаждением (регенеративные теплообменные ап-
параты).
Лекция 4
Охлаждение (нагревание) неограниченной пластины
Рис. 4.1. Охлаждение неограниченной пластины
Задача симметрична относительно оси, qv =0 . По (1.39)
∂t = a ∂2t . ∂τ ∂x2
Введём избыточную температуру
ϑ = t −tж .
Из (4.1) и (4.2) можно получить
∂ϑ =a ∂2ϑ . ∂τ ∂x2
(4.1)
(4.2)
(4.3)
37
Условия однозначности:
1.Геометрические: ширина и длина пластины намного больше толщины.
2.Физические: λ, с, ρ , a – const.
3. |
Временные: |
τ =0 ; |
t(x,0) = f (x) ; J0 = f (x) −tж . |
||
4. |
Граничные: |
пусть тело остывает. |
|||
|
|
x =0 ; |
∂ϑ =0 ; |
||
|
|
|
∂x |
|
|
|
|
x = δ ; |
∂ϑ =− |
α |
ϑ; |
|
|
|
|||
|
|
|
∂x |
λ |
−λ |
∂ϑ =α (tc −t |
ж ). |
|
∂x |
|
Цель – получить зависимость вида: ϑ = f (x,τ) .
Введём безразмерные координаты и числа подобия:
∙ относительное превышение температуры
Θ = |
t −tж |
= |
ϑ |
. |
(4.4) |
to −tж |
|
||||
|
|
ϑo |
|
∙ относительные координаты
Χ = |
x |
; |
(4.5) |
|
|||
|
δ |
|
∙число Био – число краевого подобия, характеризует связь между по-
лем температур твёрдого тела и условиями теплоотдачи на его поверх-
ности:
Bi = |
α δ |
; |
(4.6) |
|
|||
|
λ |
|
∙число Фурье – безразмерное время, характеризует связь между скоро-
стью изменения температурного поля и физическими характеристиками и размерами тела:
Fo = |
α τ |
. |
(4.7) |
|
|||
|
δ2 |
|
Цель: получить безразмерное превышение температуры над темпе-
ратурой окружающей среды, т.е. перейти от ϑ = f (x,τ) к θ = f (Χ, Bi, Fo) .
38
Необходимо решить уравнение (4.3). Решение ищем в виде
|
|
ϑ =ϕ(τ) ψ(x) . |
|
(4.8) |
||||
Продифференцируем (4.8) и подставим в (4.3): |
|
|
||||||
|
|
∂ϕ(τ) ψ(x) =a ∂2ψ(x) ϕ(τ) . |
|
(4.9) |
||||
|
|
∂τ |
|
|
∂x2 |
|
|
|
|
|
ϕ (τ) |
=a |
ψ′′(x) |
|
|
||
|
|
ϕ(τ) |
|
. |
|
(4.10) |
||
|
|
ψ(x) |
|
|||||
Уравнение (4.10) интегрируют с учётом условий однозначности, по- |
||||||||
сле преобразований, перейдя к безразмерным переменным, имеем. |
|
|||||||
∞ |
2 sin µn |
|
|
|
|
|
|
|
Θ = ∑ |
|
2 |
|
|
||||
|
|
|
cos(µn Χ) exp(−µn |
Fo) , |
(4.11) |
|||
|
|
|
||||||
n=1 µn |
+sin µn cos µn |
|
|
где µn – корень характеристического уравнения:
|
ctg µn = |
µn |
. |
(4.12) |
|
|
|||
|
|
Bi |
|
|
|
При Fo ≥0,3 берут 1-ый член ряда (µn =µ1 ), получают |
|
||
∙ |
на оси: |
|
||
|
ΘΧ=0 =N(Bi) exp(−µ12 Fo); |
(4.13) |
||
∙ |
на поверхности: |
|
||
|
ΘΧ=1 =P(Bi) exp(−µ12 Fo). |
(4.14) |
Функции N(Bi) и P(Bi) затабулированы и для зависимостей (4.13) и (4.14) имеются графики.
Количество теплоты, отданное пластиной в процессе охлаждения
Количество внутренней энергии тела относительно среды (в началь-
ный момент времени):
Qп = 2 δ f ρ c (tо −tж ). |
(4.15) |
Рассмотрим процесс от τ =0 до τ = τ1 .
К моменту времени τ = τ1 внутренней энергии осталось:
39
|
|
|
|
|
Q1 =2 δ f ρ c ( |
|
|
|
|
|
|
|
|
|
(4.16) |
|||||
|
|
|
|
|
t1 −tж ), |
|
|
|
|
|
||||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
t1 – средняя температура тела. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Изменение внутренней энергии тела: |
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
−t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
( |
|
|
|
ж ) |
1 |
ж |
п ( |
1) |
|
||||||
|
|
Q = Q |
п |
1 |
t |
о |
−t |
|
|
|
|
|
, (4.17) |
|||||||
|
|
|
|
|
|
|||||||||||||||
|
|
|
−Q |
=2 δ f ρ c |
|
|
1− |
|
|
|
|
= Q 1−Θ |
||||||||
|
|
|
|
|
|
|
|
|
|
|
tо −t |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
ж |
|
|
|
|
где Θ1 – средняя безразмерная избыточная температура к моменту вре-
мени τ = τ1 .
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t1 −tж |
|
; |
|
|
|
(4.18) |
|||||
|
|
|
|
|
|
|
|
Θ |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
tо −tж |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
Χ |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
= |
|
∫ Θ dΧ . |
|
|
(4.19) |
||||||||
|
|
|
|
|
|
|
|
Θ |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
Χ |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Подставим (4.11) в (4.19) и проинтегрируем: |
|
|
|
||||||||||||||||||||
|
|
|
|
|
∞ |
2 sin2 µn |
|
|
2 |
|
|
|
|||||||||||
|
|
|
Θ = ∑ |
|
|
|
|
|
|
||||||||||||||
|
|
|
2 |
+µn sin µn cos µn |
|
exp(−µn |
Fo) . |
|
(4.20) |
||||||||||||||
|
|
|
|
|
n=1 µn |
|
|
|
|
|
|
||||||||||||
При Fo ≥0,3 берут 1-й член ряда (µn =µ1 ), получают |
|
|
|||||||||||||||||||||
|
|
|
|
2 sin2 µ1 |
|
|
|
|
|
|
|
2 |
|
|
|
|
2 |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Θ = |
|
exp(−µ1 Fo) |
=M(Bi) exp(−µ1 |
Fo) . |
(4.21) |
||||||||||||||||||
µ12 +µ1 sin µ1 cos µ1 |
Значения функции М(Bi) затабулированы.
(4.22)
Аналогично решены задачи для бесконечного цилиндра и шара [1].
Влияние числа Вi на процессы нестационарной теплопроводности
1) Bi →∞ (Bi >100 ).
40
Рис. 4.2. Влияние числа Bi на процессы нестационарной теплопроводности Bi →∞
Здесь Fo0 <Fo1 <Fo2 <Fo3 .
Bi = α δ = 1 δ = Rλ . λ 1 α λ Rα
(4.23)
Если Rλ Rα Bi →∞.
Значение коэффициента теплоотдачи большое, охлаждение (нагре-
вание) тела определяется физическими свойствами и размерами тела.
2) Bi →0 (Bi <0,1).
Рис. 4.3. Влияние числа Bi на процессы нестационарной теплопроводности Bi →0
41
Здесь Fo0 <Fo1 <Fo2 <Fo3 .
Если Rα Rλ Bi →0 .
Процесс охлаждения (нагрева) определяется интенсивностью тепло-
отдачи на поверхности пластины.
3) 0,1<Bi <100 .
Рис. 4.4. Влияние числа Bi на процессы нестационарной теплопроводности 0,1<Bi <100
Здесь Fo0 <Fo1 <Fo2 <Fo3 .
Если Rα ≈Rλ 0,1<Bi <100 .
Охлаждение (нагрев) определяется как внутренним, так и внешним термическим сопротивлением.
Охлаждение тел конечных размеров
Эти расчёты основываются на теореме перемножения решений.
Решением задачи представляется произведение безразмерных тем-
ператур для тел неограниченных размеров, в результате пересечения ко-
торых образовалось рассматриваемое тело.
Задачи нестационарной теплопроводности сведены в табл. 1.
42