Часть 1 Основы теплопередачи
.pdf
grad t |
|
|
t |
cos(n,x ) |
t |
|
x |
n |
x |
||||
|
|
|
;
cos(n,x ) |
n |
|
n |
. |
|
x |
x |
||||
|
|
|
Коэффициент теплопроводности λ характеризует способность тел проводить теплоту, его величина характеризует количество теплоты, которое проходит в единицу времени через единицу площади поверхности при температурном grad = 1. Коэффициент теплопроводности, как правило, определяется опытным путем и зависит от температуры и давления:
= |
|q | |
|
|
grad |
t |
||
|
,Вт/(м∙К).
(1.6)
Таблица. 1.1
Ориентировочные значения λ различных веществ при t=20 °C
Вещество |
|
λ Вт/(м К) |
Воздух |
|
0,025 |
Вода |
|
0,6 |
|
Твердые материалы |
|
Асбест |
|
0,11 – 0,3 |
Кирпич |
|
0,8 |
Пенопласт |
|
0,05 |
Стекловата |
|
0,047 |
|
Металлы |
|
Углеродистая сталь |
|
45 – 60 |
Медь |
|
390 |
Количество теплоты, проходящее через единицу площади поверхности в единицу времени, называют плотностью теплового потока:
Вт
м
(1.7)
Плотность теплового потока есть вектор, положительное направление которого указывает в сторону уменьшения температуры (рис. 1.1). Тогда основной закон теплопроводности может быть записан в виде
(1.8)
- 13 -
Направления grad t и q противоположны, этим объясняется знак «–» в законе Фурье.
Для подсчета q в любой точке тела необходимо знать значение grad t = ∂t/∂n. Температурное поле внутри тела может быть определено в общем случае для конкретной задачи лишь путем решения дифференциального уравнения теплопроводности.
1.3. Дифференциальное уравнение теплопроводности
Для упрощения вывода сделаем следующие допущения: будем считать, что тело однородно и изотропно. Его физические параметры, например λ= const, P = const ρ const, не зависят от температуры. Пренебрежем также деформацией объема вследствие переноса теплоты, так как она мала по сравнению с самим объемом. Вывод дифференциального уравнения теплопроводности основан на применении закона сохранения энергии, который применительно к данному случаю может быть сформулирован следующим образом: количество теплоты dQ, введенное в элементарный объем dV за промежуток времени dτ, вследствие теплопроводности идет на изменение внутренней энергии вещества dU, содержащегося в данном объеме:
dQ =dU=ρCV |
t |
dVdτ ,кДж, (1.9) |
|
|
|||
|
|
где ρ – плотность, кг/м³; CV, – теплоемкость при постоянном объеме,
кДж/(кг∙К); |
t |
, К/с; dV, м³; τ, с. |
|
||
|
|
|
Выделим в теле элементарный объем (рис.1.2). Обозначим через dQx, dQy и dQz количества теплоты, которые подводятся к граням элементарного объема в направлении осей 0X, 0Y, 0Z. Обозначим также через dQx+dx, dQy+dy и dQz+dz количества теплоты, которые отводятся от противоположных граней элементар-
Рис. 1.2 ного объема в направлении тех же осей.
- 14 -
Подсчитаем dQx = qxdydzdτ (dydz = dF), где qx – проекция вектора плотности теплового потока на ось 0X:
dQx+dx = qx+dxdydzdτ.
Подсчитаем избыток теплоты в направлении оси 0X и всех остальных осей:
dQx– dQx+dx = (qx– qx+dx) dydzdτ = …
Для получения этого избытка теплоты функцию (qx+dx) разложим в ряд Тейлора с точностью до второго члена ряда:
qx+dx= qx+
q |
x |
||
|
|
||
x |
|||
|
|||
|
|
||
dx+…
и подставим значение функции qx+dx в избыток теплоты. Получим
dQ |
x |
|
dQ |
x dx |
|
q |
|
q |
|
|
|
q |
x |
|
dx |
|
dydzd |
|
q |
x |
|
dVd |
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
x |
|
x |
|
x |
|
|
|
x |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
,
где dxdydz = dV.
По аналогии
dQ |
|
dQ |
|
|
q |
y |
|
dVd |
;dQ |
|
dQ |
|
|
q |
z |
|
dVd . |
||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
y |
|
y dy |
|
|
y |
|
|
|
z |
|
z dz |
|
z |
|
||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Согласно закону сохранения энергии этот избыток Q идет на изменение внутренней энергии вещества, содержащегося в данном объеме:
q |
x |
|
q |
y |
|
q |
z |
|
||
dQ |
|
|
|
|
|
dVd |
||||
x |
y |
z |
||||||||
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|||
C |
|
t |
dVd |
||
v |
|
|
|||
|
|
||||
|
|
|
|
||
.
Тогда
t |
|
1 |
q |
x |
|
q |
y |
|
q |
z |
|
||
|
|
|
|
|
|
. |
|||||||
|
C |
x |
y |
z |
|||||||||
|
v |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||
По закону Фурье
q |
x |
q |
n |
cos(n,X |
|
|
|
)
|
t |
cosn,X ; |
t |
cos(n,X ) |
t |
. |
|
n |
n |
x |
|||||
|
|
|
|
q |
x |
|
|
|
|
t |
|
|
|
2t |
, |
|
|
|
|
|
|
|
|
|
|||||
x |
|
|
|
x |
2 |
|||||||
|
x |
|
x |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
- 15 - |
|
|
|
q |
|
|
|
|
|
t |
|
|
t |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
y |
|
|
|
|
|
|
|
|
|
y |
|
|
y |
2 |
||||||
|
|
|
|
|
||||||
|
y |
|
y |
|
|
|||||
|
q |
|
|
|
|
t |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
2 |
|
|
. |
|
|
|
|
|
|
|
|||
|
|
|
z |
z |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
z |
|
|
|
2 |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
t |
|
|
|
|
2t |
|
|
|
2t |
|
2t |
|
, |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
x |
2 |
|
y |
2 |
z |
2 |
||||||||
|
|
|
v |
|
|
|
|
|
|
|
|
|
||||||
,
C |
v |
|
a
,
где a – физический параметр вещества, характеризующий скорость изменения температуры в нестационарных процессах, который был назван коэффициентом температуропроводности. Это мера теплоинерционных свойств тела. Чем больше значение а, тем больше скорость изменения температуры в теле и тем меньше тепловая инерционность тела. Наивысшее значение для а характерно для металлов. Размерность коэффициента температуропроводности:
|
|
3 |
кг К |
|
2 |
|
а |
|
Вт м |
= |
м |
, |
|
|
|
|
||||
|
|
м К кг Дж |
|
с |
|
|
|
|
|
|
|||
тогда
t |
2 |
|
а t , |
|
2 оператор Лапласа.
(1.10)
Дифференциальное уравнение теплопроводности устанавливает связь между пространственным и временным изменением температуры. Рассмотрим применение дифференциального уравнения теплопроводности для решения некоторых практических задач теплопроводности.
- 16 -
1.4. Применение дифференциального уравнения теплопроводности для решения практических задач 1.4.1. Теплопроводность плоской однослойной стенки при стационарном режиме
Сформулируем задачу. При этом будем считать, что два размера стенки, например, в направлении осей Y и Z являются существенно бóльшими по сравнению с ее толщиной (рис. 1.3), и поэтому
2t 0;
у 2
Тогда
t |
|
2 |
|
z |
2 |
|
|
0;
t |
0. |
|
|
||
|
t=f(x). (1.11)
При стационарном режиме на поверхностях стенки установились во времени тем-
пературы tc1 и tc2.
Для получение формулы для плотности теплового потока необходимо решить дифференциальное уравнение теплопроводности для сформулированной нами
задачи, чтобы получить закон изменения температуры в данной стен-
ке. Поскольку t=f(X), то |
2t |
0 |
– и мы имеем дифференциальное |
|
x 2 |
||||
|
|
|
уравнение теплопроводности для стационарной одномерной задачи. Первое интегрирование дает
t |
C . |
|
x |
||
1 |
||
|
Второе интегрирование дает закон изменения температуры в такой стенке: t = C1X+C2. Это уравнение прямой. Значит, в плоской стенке температура меняется по линейному закону и не зависит от свойств материала.
Запишем граничные условия:
1. X=0, t=tc1=C2;
- 17 -
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
||
2. X= ,t=tc2=C1∙ +tс1.Отсюда |
|
|
|
|
|
|
t |
c 2 |
c 1 |
|
|
. |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Запишем закон Фурье: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
t |
|
t |
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Вт |
|
|||
|
|
q |
|
|
c 2 |
|
c 1 |
|
|
|
|
|
|
t |
|
|
|
, |
, |
|||||||
|
|
x |
|
|
|
|
|
|
t |
c 1 |
c 2 |
|
м |
2 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
где |
|
– тепловая проводимость плоской однослойной стенки. Отсюда |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
плотность теплового потока: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
q |
t |
c 1 |
t |
c 2 |
, |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где – термическое сопротивление теплопроводности плоской одно-
слойной стенки.
Если известна площадь поверхности, через которую передается
теплота, то |
|
|
|
|
|
|
|
t |
|
||
Q = qF= |
t |
c 1 |
c 2 |
|
|
|
|
, Вт. |
|||
|
|
|
|
||
|
|
|
|
|
Площадь теплового потока пропорциональна температурному напору на поверхности плоской стенки tc1 – tc2 и обратно пропорциональна термическому сопротивлению.
1.4.2. Теплопроводность плоской многослойной стенки
Допустим, что плоская стенка состоит из n разнородных слоев толщиной 1, 2… n с коэффициентами теплопроводности λ1, λ2…. λn (рис. 1.4). Слои такой стенки плотно прилегают друг к другу, поэтому температура на соприкасающихся поверхностях одна и та же.
При стационарном режиме в стенке установились во времени температуры tc1,tc2,tc3,tcn,tc(n+1). Тогда для первого слоя
- 18 -
|
|
|
t |
|
|
|
|
q |
t |
c 1 |
c 2 |
|
. |
||
|
|
|
|
||||
|
|
|
|
|
|
||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
При |
|
стацио- |
|||||
нарном режиме |
чис- |
||||||
ленно то же значение плотности теплового
потока для |
|
второго |
|||||
слоя: |
|
|
|
|
|
|
|
|
|
t |
|
|
|||
q |
t |
c 2 |
c 3 |
|
. |
||
|
|
|
|
||||
|
|
|
|
|
|
||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
Рис. 1.4
Докажем, что это так.
Для этого используем промежуточное выражение дифференциального уравнения теплопроводности
|
|
|
|
( |
|
+ |
|
+ |
|
|
) |
|
||
|
|
|
|
|
|
|
||||||||
Это объясняется тем, что при стационарном режиме |
t |
0 |
, |
|||||||||||
|
||||||||||||||
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
(1.12)
1 |
0 |
. |
|
|
|||
|
|||
v |
|
||
|
|
Тогда |
q |
x |
0 |
, qx |
= const. Для n-го слоя значение плотности теплово- |
||||||||
|
|||||||||||||
x |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||
го потока составит |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
t |
|
|
||||
|
|
|
|
|
q |
t |
cn |
c (n 1) |
|
. |
|||
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
n |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
n |
|
|
||
Записанные исходные соотношения переписываем каждое относительно температурного напора и складываем:
- 19 -
|
|
|
|
|
|
1 |
|
|
2 |
|
|
n |
|
t |
|
t |
|
q |
|
|
|
..... |
|
. |
|||
c 1 |
c (n 1) |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
1 |
|
|
2 |
|
|
n |
||
Тогда
|
t |
c 1 |
t |
c (n 1) |
|
|
Вт |
|
||
q |
|
|
, |
. |
||||||
|
|
n |
|
|
|
м |
2 |
|||
|
|
|
|
|
|
|
||||
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i 1 |
|
i |
|
|
|
|
|
Последнее соотношение может быть переписано в виде
( ( ))
|
|
|
Вт |
(1.13) |
∑ |
|
|
||
|
|
|
||
|
|
|
|
|
|
|
|
|
Заметим, что знаменатель представляет собой сумму термических сопротивлений многослойной плоской стенки. Его величина не меняется от перестановки слоев многослойной стенки, и поэтому не изменяется плотность теплового потока и сам тепловой поток.
1.4.3. Теплопроводность однослойной цилиндрической стенки
Задача решается в цилиндрических координатах (рис. 1.6). При этом уравнение температурного поля в общем виде запишется:
t=f (τ,r,φ,z) |
(1.14) |
Переход от полярных (рис. 1.5) к цилиндрическим координатам приводит к дифференциальному уравнению теплопроводности в общем виде
t |
t |
|
1 t |
|
1 |
t |
|
t |
|
|||||
|
|
2 |
|
|
|
|
|
|
2 |
|
|
2 |
|
|
|
a |
r |
|
|
r r |
|
r |
|
|
|
|
z |
|
. |
|
2 |
|
|
2 |
2 |
|
2 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Сформулируем |
|
|
задачу. |
||||||||||
|
Пусть l >> r2 – r1, тогда |
2t |
0. |
Рис. 1.5 |
z 2 |
- 20 -
Будем считать, что температура в цилиндрической стенке не меняется по окружности цилиндра, поэтому
|
2 |
t |
|
|
|
|
|
|
|
|
|
|
2 |
0 |
, |
а также |
|
|
|||||
|
|
|
|||
|
|
|
|
|
|
t |
0 |
(так как ста- |
|
|
|||
|
|||
|
|
ционарный режим). Тогда дифференциальное уравнение для сформулированной нами задачи запишется в виде
2 |
|
1 dt |
|
||
d |
t |
|
0. |
||
r |
2 |
r dr |
|||
|
|
||||
|
|
|
|||
Введем новую переменную u:
Рис. 1.6
Разделим переменные:
dt |
u |
||
dr |
|||
|
|
||
du |
|
dr |
|
u |
r |
||
|
|||
,
du |
|
u |
0. |
|
dr |
r |
|||
|
|
|||
0 |
|
|
|
и проинтегрируем: lnu+lnr = const |
1. Перепишем последнее выра- |
|||
жение в виде u∙r = 1'. Подставим значение u: |
||||
dt |
dt |
C1 ' |
|
|
dr r |
1 ', dr |
|
. |
|
r |
|
|||
Разделим переменные: dt 1 'drr .
Проинтегрируем:
t= 1'∙lnr + 2.
- 21 -
В цилиндрической стенке температура изменяется по логарифмическому закону и не зависит от свойств материала.
|
Применим граничные условия: |
|
1. |
r=r1 то t=tc1= 1'∙lnr1 + |
2 ; |
2. |
r=r2,то tc2 = 1'∙lnr2 + |
2. |
Из выражения 2 вычтем выражение 1 и получим
t |
|
t |
|
C 'ln |
r |
, |
|
|
2 |
||||
|
|
|
|
|
|
|
|
c 2 |
|
c 1 |
1 |
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
1 ' tc 2 tc 1 .
Применим закон Фурье:
Q
dt |
F |
|
dr |
||
|
,
где F – площадь боковой поверхности цилиндра.
Тогда тепловой поток через цилиндрическую стенку запишется
Q |
1 |
' |
2 r |
|
|
|
'2 . |
|||
|
|
|
||||||||
r |
|
|
1 |
|||||||
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
Подставим значение 1' : |
|
|
|
|
|
|
|
|
|
|
Q |
(t |
c 1 |
t |
c 2 |
) |
,Вт |
||||
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
||
|
|
|
1 |
ln |
d |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
d |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где
1 |
ln |
d |
|
2 |
d |
||
|
2 1
представляет собой линейное термическое сопротивле-
ние теплопроводности однослойной цилиндрической стенки.
Часто тепловой поток через цилиндрическую стенку относят к
единице трубы , тогда |
|
|
|
q |
|
Q |
, Вт/м, |
|
где qℓ - линейная плотность теплового потока.
q |
Q |
|
(t |
c 1 |
t |
c 2 |
) |
, |
Вт |
. |
|
||||
|
|
|
|
|
|
|
(1.15) |
||||||||
|
|
1 |
|
|
d2 |
|
|
м |
|||||||
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
ln |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
2 |
|
|
d1 |
|
|
|
|
|
||||
|
|
|
- 22 - |
|
|
|
|
|
|
|
|
|
|||
