tot_book
.pdf141
Рис. 77.
(во всех одномерных стационарных задачах следовало бы использовать обозначения dT, dх и т. д.; сохраним знак частной производной лишь для однообразия и в силу традиции). Задача (2.20) допускает прямое интегрирование:
|
∂T |
= C |
(2.21) |
|
|
||
|
∂x |
1 |
|
|
|
|
|
T = C1 x + C2 |
(2.22) |
||
|
где С1, С2 — постоянные интегрирования. |
|
|||||||
При х = 0 из равенства (2.22) следует, что Тw1 = С2; при х = δ |
|
|||||||
Тw2 |
= С1 |
+ Тw1, откуда C |
= |
Tw2 −Tw1 |
. Получаем |
|
||
|
|
|
||||||
|
|
1 |
|
|
δ |
|
||
|
|
|
|
|
|
|||
|
|
T = |
Tw2 −Tw1 |
x +T |
(2.23) |
|||
|
|
|
||||||
|
|
|
|
δ |
|
|
w1 |
|
|
|
|
|
|
|
|
|
Ранее были известны лишь температуры Тw1 и Тw2 уравнение (2.23) позволяет определить температуру в любом сечении пластины и показывает, что по координате х температура меняется линейно.
Возможен и другой случай: если Т(х) задана, то постоянную во всех сечениях пластины плотность теплового потока найдем,
подставляя в уравнение (2.4) значение ∂T ∂x и уравнения (2.21):
|
|
|
|
142 |
|
|
|
|
|
|
q = λ |
Tw1 −Tw2 |
|
(2.24) |
|
|
|
|
δ |
|
|||
|
|
|
|
|
|
|
|
Величину |
δ |
= R |
|
называют |
внутренним |
термическим |
|
λ |
|
||||||
|
λ |
|
|
|
|
сопротивлением4 ; ее размерность — м2 К/ Вт; очевидно, что
q = Tw1 −Tw2 Rλ
Решения (2.23) и (2.24) можно обобщить для многослойной пластины (рис. 78).
|
|
|
Рис. 78. |
|
4 Заметим, что |
R |
= |
Tw1 −Tw2 |
— отношение разности |
|
||||
|
λ |
|
q |
|
|
|
|
температур на двух изотермических поверхностях к плотности теплового потока. Для пластины q = idem, в более общем случае q может задаваться для одной из изотермических поверхностей.
143
Пусть в пластине n слоев с толщинами δ1, δ2,...,δn и теплопроводностями λ1, λ2,...,λn, соответственно; температуры границ слоев равны Тw1, Тw2,...,Тwn+1. Физически ясно, что в стационарном режиме плотность теплового потока q во всех сечениях пластины остается одинаковой:
q = Tw1 −Tw2 ;
δ1 λ1
q= Tw2 −Tw3 ;
δ2 λ2
.......................
q= Twn −Twn+1
δn λn
Разрешим эти уравнения относительно разностей температур, а затем сложим почленно:
Tw1 −Tw2 = qδ1 λ1 ;
Tw2 −Tw3 = qδ2 λ2 ;
....................... |
|
|
|
|
|
|
|
|
|
|
|||||
T |
|
−T |
|
= qδn |
; |
|
|
|
|
|
|
|
|
||
wn |
|
wn+1 |
|
|
|
λn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
T |
−T |
|
|
|
δ |
+ |
δ |
2 |
+... + |
δ |
n |
|
= |
||
|
= q |
1 |
|
|
|
|
|||||||||
w1 |
|
wn+1 |
|
|
|
λ1 |
|
|
|
λ2 |
|
|
λn |
|
|
= q(Rλ1 + Rλ2 +... + Rλn ), поэтому |
|
|
|
|
|||||||||||
q = |
|
|
Tw1 −Twn+1 |
|
|
|
|
|
|
|
|
||||
|
R |
+ R |
|
+... |
+ R |
|
|
|
|
|
|
|
(2.25) |
||
|
|
λ1 |
λ2 |
|
|
λn |
|
|
|
|
|
||||
Величину |
Rλ1 + Rλ2 +... + Rλn = RλΣ |
|
|
назовем суммарным |
внутренним термическим сопротивлением многослойной пластины (здесь видна аналогия с сопротивлением электрических цепей постоянного тока: суммарное сопротивление
144
последовательно соединенных проводников равно сумме составляющих). В осях Т–х распределение температур в сечении пластины представляет ломаную линию; наклон любого ее звена тем больше, чем меньше теплопроводность соответствующего слоя. Усложним задачу (рис. 79): пусть обе поверхности пластины омывают жидкости с различными температурами (Тf1 > Тf2), а коэффициенты теплоотдачи на этих поверхностях равны соответственно α1 и α2. Рассмотрим, таким образом, задачу теплопроводности при граничных условиях III рода.
Поскольку в любой стационарной задаче теплопроводности плотность теплового потока во всех сечениях пластины остается неизменной, определим значение q: на левой границе
q |
|
|
= α |
(T |
|
−T |
), внутри пластины |
q |
|
|
= |
Tw1 −Tw2 |
и на |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
||||||||
|
|
x=0 |
1 |
|
f 1 |
w1 |
|
|
|
0<x<δ |
|
δλ |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
правой границе q x=δ = α2 (Tw1 −Tf 2 ).
Рис. 79.
Выполнив те же преобразования, что и в предыдущем случае, получим
q = |
|
Tf 1 −Tf 2 |
|
= |
Tf 1 |
−Tf 2 |
, |
(2.26) |
||
1 |
|
+ δ + 1 |
|
R |
+ R + R |
|||||
|
|
|
|
|
|
|||||
|
|
α1 |
λ |
α2 |
|
α1 |
|
λ α2 |
|
|
|
|
|
|
|
|
|
|
145
где |
R |
= 1 |
α1 |
, R |
= 1 |
α2 |
— внешние |
термические |
||||
|
α1 |
|
|
α2 |
|
|
|
|
|
|
||
сопротивления, (м2 К)/Вт — величины, |
обратные коэффициентам |
|||||||||||
теплоотдачи α1 и α2. |
|
|
|
|
|
|
|
|
|
|||
|
Величину |
k = |
|
1 |
|
|
, |
которая |
характеризует |
|||
|
1 α +δ λ +1 α |
2 |
||||||||||
|
|
|
|
1 |
|
|
|
|
|
|
интенсивность передачи теплоты, называют. Из формулы (2.26) следует, что значение k численно равно отношению плотности теплового потока, передаваемого через пластину от одной жидкости к другой, к разности между температурами жидкостей
(теплоносителей)5.
Если величина α характеризует только конвективный теплообмен на поверхности пластины, то коэффициент
теплопередачи k является характеристикой всей теплопередающей системы, поскольку кроме коэффициентов теплоотдачи α1 и α2 она включает параметры δ и λ!
Величина 1k = Rk — ясно, что Rk = Rα1 + Rλ + Rα2 . Распределение температур в многослойной пластине при
граничных условиях III рода в пределах каждого слоя остается линейным, а в жидкости или газе вблизи границ пластины приобретает нелинейный характер (рис. 80); плотность теплового потока
q = |
|
|
Tf 1 −Tf 2 |
|
|
= |
Tf 1 −Tf 2 |
= k(T |
f 1 |
−T |
f 2 |
). |
(2.27) |
||
|
|
n |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
Rk |
|
|
|||||
1 |
|
+ ∑δ |
|
+ |
1 |
|
|
|
|
|
|
||||
α1 |
λ |
α2 |
|
|
|
|
|
||||||||
|
|
i=1 |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 Подробнее о теплоносителях см. разд. 2.12.
146
Рис. 80.
Нередко требуется определить температуру в произвольном сечении пластины 0 < х < δ. Запишем значения перепадов температур, используя прежнее условие q = idem:
Tf 1 −Tw1 = q α1 ;
Tw1 −Tw2 = qδ1 λ1 ;
.......................
Twn −Twn+1 = qδn λn ;
Twn+1 −Tf 2 = q αn+1 .
Температуры на границах слоев
Tw1 =Tf 1 −q 1α1
Tw2 =Tf 1 −q 1α1 +δ1 λ1
.....................................................
Twn+1 =Tf 1 −q 1α1 +δ1 λ1 +...+δn λn ;
Tfn+1 =Tf 1 −q 1α1 +δ1 λ1 +... +δn λn + 1αn+1 .
Пусть сечение х находится в j-м слое стенки. Тогда расстояние от этого сечения до границы с температурой Тwj будет
равно xj = х – (δ1 + ... + δj–1), а температура
147
T (x)=T |
|
|
1 |
|
+ |
δ |
+...+ |
δj−1 |
+ |
x j |
|
(2.28) |
f 1 |
−q |
α1 |
1 |
|
|
. |
||||||
|
|
|
|
|
λ1 |
|
λj−1 |
|
λj |
|
Рассуждения остаются в силе, если расчет начинать не с Тf1, а с Тf2 перед q поставить знак “+”, а все члены суммировать в обратном порядке.
Теплопроводность цилиндрической стенки при граничных условиях I рода рассчитать чуть сложнее, поскольку в этом случае плотность потока q зависит от текущего радиуса r (рис. 81).21.10
Рис. 81.
Обычно цилиндрические стенки достаточно “длинные” (трубы, каналы и т. п.); их поле температуры вполне точно описывают одномерным уравнением Фурье, т. е. считают, что температура вдоль оси цилиндра постоянна.
Краевая задача Фурье для осесимметричной цилиндрической стенки без внутренних источников имеет вид
∂2T |
+ |
1 ∂T |
= 0; |
||
∂r2 |
|
|
|||
r ∂r |
|||||
|
(2.29) |
||||
Tr=r1 =Tw1; |
|||||
|
Tr=r2 =Tw2.
Введем новую переменную U = ∂∂Tr , чтобы понизить порядок производных, тогда
|
|
|
|
|
148 |
|
|
|
|
∂2T |
= |
∂U |
; |
1 ∂T |
= U |
; |
∂U |
+U |
= 0; |
∂r2 |
|
∂r |
|
r ∂r |
r |
|
∂r |
r |
|
разделим переменные и проинтегрируем; получим последовательно
lnU +ln r = C1, |
|
||
∂ = |
∂r |
, |
(2.30) |
T |
C1 ∂r |
|
T = C1 ln r +C2.
Постоянные С1 и С2 найдем из граничных условий задачи
(2.29):
при |
r = r1 |
|
T =Tw1; |
Tw1 = C1 ln r1 +C2; |
|
|
|
||||||||
при |
r = r2 |
|
|
T =Tw2; |
Tw2 = C1 ln r2 +C2; |
|
|
|
|||||||
C = |
Tw1 −Tw2 |
|
; C |
|
=T |
−(T |
−T |
) |
ln r1 |
|
|
. |
|||
|
|
|
ln(r r |
) |
|||||||||||
1 |
ln(r r |
) |
|
|
2 |
|
w1 |
w1 |
w2 |
|
|
||||
|
1 |
2 |
|
|
|
|
|
|
|
|
|
1 |
2 |
|
|
После подстановки получим значение температуры Т на произвольном радиусе r:
T (r)=T |
−(T |
−T |
) |
ln r r1 |
|
. |
(2.31) |
||
ln(r r |
) |
||||||||
w1 |
w1 |
w2 |
|
|
|
||||
|
|
|
|
1 |
2 |
|
|
|
Уравнение (2.31) описывает убывающую логарифмическую функцию, что физически вполне объяснимо: с увеличением значения r плотность теплового потока снижается, хотя тепловой поток, подводимый на поверхности r =r1, и отводимый с
поверхности r = r2, остается неизменно равным Q = −λ ∂∂Tr 2πrl
(здесь l — длина расчетного участка).
Поскольку в соответствии с равенством (2.30)
∂T |
= |
C1 |
= |
(Tw1 |
−Tw2 ) |
, |
Q = |
2πλl(Tw1 −Tw2 ) |
|
, Вт. |
|||
∂r |
r |
|
|
r ln(r r ) |
|||||||||
|
|
r ln(r r |
) |
|
|
|
|||||||
|
|
|
|
|
1 |
2 |
|
|
|
1 |
2 |
|
|
Для цилиндрической стенки плотность теплового потока удобнее рассматривать либо на поверхности 2πr1l, либо на поверхности 2πr2l:
|
|
|
|
|
149 |
|
|
|
|
q1 = |
|
Q |
= |
(Tw1 −Tw2 ) |
, |
Вт м2 , |
(2.32) |
||
|
2πr1l |
||||||||
|
|
|
|
(r1 λ)ln(r2 r1 ) |
|
|
|
|
|
q2 = |
|
Q |
|
= |
(Tw1 −Tw2 ) |
|
, |
Вт м2 , |
(2.33) |
2πr2l |
|
|
|||||||
|
|
|
(r2 λ)ln(r2 r1 ) |
|
|
|
Можно также рассчитать тепловой поток, протекающий через единицу длины цилиндра,
ql = |
Q |
= |
π(Tw1 −Tw2 ) |
, Вт м, |
(2.34) |
|
l |
||||||
|
|
(1 2λ)ln(r2 r1 ) |
|
|
эту величину назовем линейной плотностью теплового потока.
Легко заметить, что
ql = 2πr1q1 = 2πr2q2; |
(2.35) |
иначе говоря, линейная плотность теплового потока ql , не зависит от текущего радиуса r. Для n-слойной цилиндрической стенки все рассуждения проведем так же, как для многослойной пластины; при этом линейная плотность теплового потока окажется равной
ql = |
π(Tw1 −Tw2 ) |
, Вт м. |
|
||||
n |
|
1 |
ln(ri+1 ri ) |
(2.36) |
|||
|
∑ |
|
|
|
|||
|
i=1 |
2λi |
|
|
|
|
|
Величину |
R |
|
= |
1 ln(r |
r ) |
назовем |
линейным |
|
lλi |
|
i+1 |
i |
|
|
|
|
|
|
|
2λi |
|
|
|
внутренним термическим сопротивлением i-го слоя, а сумму
n |
1 ln(ri+1 ri ) — суммарным линейным внутренним |
RlλΣ = ∑ |
|
i=1 |
2λi |
термическим сопротивлением многослойной цилиндрической стенки, (м·К)/Вт.
Для граничных условий III рода (рис. 82)
ql = |
|
π(Tf 1 −Tf (n+1)) |
|
|
, Вт м, |
(2.37) |
|
n |
|
|
|
|
|||
2α1 r |
21λ |
ln(ri+1 ri )+ |
|
1 |
|
||
+ ∑ |
|
|
|
||||
2α |
r |
|
|||||
1 1 |
i=1 |
|
i |
|
n+1 n+1 |
|
150
Рис. 82/
где |
1 |
= R |
|
, |
1 |
= R |
|
— линейные |
внешние |
|
|
|
|
||||||
|
|
lα |
|
|
lα |
|
|
|
|
|
2α1r1 |
1 |
|
2α2rn+1 |
2 |
|
|
||
|
|
|
|
|
|
|
|
термические сопротивления, (м·К)/Вт. Равенство (2.37) можно представить в виде
|
|
|
|
ql |
= kl π(Tf 1 −Tf (n+1)), |
|
Вт м, |
|
|
|
(2.38) |
||||
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
kl = |
|
|
|
1 |
|
|
|
|
= |
1 |
|
, |
Вт |
— |
|
|
|
n |
ln(ri+1 ri )+ 2α |
|
1 r |
|
Rlα +RlλΣ+Rlα |
|
м К |
||||||
|
|
2α1 r + ∑21λ |
|
1 |
n+1 |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
1 1 i=1 i |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
n+1 n+1 |
|
|
|
|
|
|
|
||||
линейный |
коэффициент |
|
теплопередачи, |
а |
величина |
||||||||||
1 |
= R , |
м К |
— полное линейное термическое сопротивление. |
||||||||||||
|
Вт |
||||||||||||||
|
|
lk |
|
|
|
|
|
|
|
|
|
|
|
|
|
kl |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Если |
расчет ведут для |
плотностей теплового |
потока |
||||||||||
q1 = k1(Tf 1 −Tf 2 ) |
и q2 = k2 (Tf 1 −Tf 2 ), |
|
где k1, k2 — коэффициенты |
теплопередачи, Вт/(м2·К), рассчитанные для соответствующих поверхностей однородного цилиндра, то справедливо соотношение kl = 2r1k1 = 2r2k2 .
Температура на радиусе r, отдаленном на расстояние ∆rj от границы (j – 1)-го и j -го слоев (см. рис. 82), равна