Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Часть 1 Основы теплопередачи

.pdf
Скачиваний:
80
Добавлен:
12.03.2015
Размер:
2.96 Mб
Скачать

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 за промежуток времени , вследствие теплопроводности идет на изменение внутренней энергии вещества 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.3

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, 2n с коэффициентами теплопроводности λ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 -

lnd2 d1

В цилиндрической стенке температура изменяется по логарифмическому закону и не зависит от свойств материала.

 

Применим граничные условия:

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 -