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

Osnovy_teploperedachi_i_massoobmena_2015

.pdf
Скачиваний:
73
Добавлен:
23.03.2016
Размер:
12.13 Mб
Скачать

тем диффузии свободных электронов и имеет место тогда, когда в различных точках тела температура не одинакова.

1.2. Основной закон теплопроводности. Закон Фурье

Перенос теплоты теплопроводностью (который в чистом виде имеет место только в твердых телах с малым коэффициентом термического расширения) выражается эмпирическим законом Био–Фурье:

(1.5)

Количество теплоты dQ, проходящее через элемент изотермической поверхности dF (м²) за промежуток времени , пропорционально градиенту температуры ∂t/∂n. Коэффициент пропорциональности λ (коэффициент теплопроводности) – это справочная характеристика, его значение различно для различных материалов, характеризует способность тел проводить теплоту.

Градиент температуры ∆n – это наикратчайшее расстояние

между соседними изотермами, тогда lim

t

=

t

, К/м, который назы-

n

n

n 0

 

 

вается градиентом температуры grad t. Градиент температуры характеризует интенсивность возрастания температуры по нормали к изо-

термической поверхности. nt – векторная производная, положитель-

ное направление которой указывает в сторону возрастания температуры. Часто требуется подсчитать изменение температуры в произвольном направлении, например в направлении оси x:

Рис. 1.1

- 13 -

grad t

 

t

cos(n,x )

t

; cos(n,x )

n

 

n

.

 

n

x

x

 

 

x

 

 

 

x

Коэффициент теплопроводности λ характеризует способность тел проводить теплоту, его величина характеризует количество теплоты, которое проходит в единицу времени через единицу площади поверхности при температурном grad = 1. Коэффициент теплопроводности, как правило, определяется опытным путем и зависит от температуры и давления:

=

|q |

, Вт/(м·К).

(1.6)

grad t

 

 

 

 

 

 

 

 

Таблица. 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)

- 14 -

Направления 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 ного объема в направлении тех же осей.

- 15 -

Подсчитаем dQx = qxdydzdτ (dydz = dF), где qx – проекция вектора плотности теплового потока на ось 0X:

dQx+dx = qx+dxdydzdτ.

Подсчитаем избыток теплоты в направлении оси 0X и всех остальных осей:

dQx– dQx+dx = (qx– qx+dx) dydzdτ = …

Для получения этого избытка теплоты функцию (qx+dx) разложим в ряд Тейлора с точностью до второго члена ряда:

qx+dx= qx+

qx

dx+…

x

 

 

и подставим значение функции

dQx dQx dx = qx

qx

qx

x

 

 

где dxdydz = dV.

По аналогии

qx+dx в избыток теплоты. Получим

dx dydzd

qx

dVd ,

x

 

 

dQ

 

dQ

 

qy

dVd ;dQ

 

dQ

 

qz

dVd .

y

y dy

 

z

z dz

z

 

 

y

 

 

 

 

 

 

 

 

 

 

Согласно закону сохранения энергии этот избыток Q идет на изменение внутренней энергии вещества, содержащегося в данном объеме:

dQ

 

 

q

x

 

 

qy

 

 

q

z

 

dVd

 

C

 

 

t

dVd .

 

 

 

 

x

 

 

 

y

 

z

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

1

 

 

 

 

q

x

 

 

 

qy

 

 

 

 

q

z

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Cv

 

 

 

 

x

 

y

 

 

 

 

z

 

 

 

 

 

 

 

 

По закону Фурье

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qx qn

 

cos(n,X )

 

 

 

t

 

cosn,X ;

t

cos(n,X )

t

.

 

 

 

n

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

qx

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

2t

,

 

 

 

 

 

 

 

 

x

 

 

x

 

 

 

 

 

 

x

 

 

 

 

 

 

 

x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- 16 -

 

 

 

 

 

 

 

 

 

 

qy

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

2t

,

 

 

 

 

y

 

 

y

 

 

 

 

 

y

 

 

 

 

 

y 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qz

 

 

 

 

 

2t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

z 2

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

2t

 

 

2t

 

2t

 

,

 

 

 

a ,

 

 

 

 

 

 

 

C

 

 

 

 

 

Сv

x

2

 

 

 

y

2

 

z

2

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где a – физический параметр вещества, характеризующий скорость изменения температуры в нестационарных процессах, который был назван коэффициентом температуропроводности. Это мера теплоинерционных свойств тела. Чем больше значение а, тем больше скорость изменения температуры в теле и тем меньше тепловая инерционность тела. Наивысшее значение для а характерно для металлов. Размерность коэффициента температуропроводности:

аВт м3 кг К=м2 ,

мК кг Дж с

тогда

t

а

2t ,

2 оператор Лапласа.

 

 

(1.10)

 

 

 

 

Дифференциальное уравнение теплопроводности устанавливает связь между пространственным и временным изменением температуры. Рассмотрим применение дифференциального уравнения теплопроводности для решения некоторых практических задач теплопроводности.

- 17 -

Рис. 1.3

1.4. Применение дифференциального уравнения теплопроводности для решения практических задач 1.4.1. Теплопроводность плоской однослойной стенки при стационарном режиме

Сформулируем задачу. При этом будем считать, что два размера стенки, например, в направлении осей Y и Z являются существенно бóльшими по сравнению с ее толщиной (рис. 1.3), и поэтому

 

2t

0;

2t

0;

t

0.

 

 

z 2

 

 

 

у 2

 

Тогда

 

 

 

 

 

 

 

 

t=f(x).

 

 

(1.11)

При стационарном режиме на поверхностях стенки установились во времени тем-

пературы tc1 и tc2.

Для получение формулы для плотности теплового потока необходимо решить дифференциальное уравнение теплопроводности для сформулированной нами

задачи, чтобы получить закон изменения температуры в данной стен-

ке. Поскольку t=f(X), то

2t

0

– и мы имеем дифференциальное

x

2

 

 

 

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

t

C

 

.

 

1

x

 

 

 

Второе интегрирование дает закон изменения температуры в такой стенке: t = C1X+C2. Это уравнение прямой. Значит, в плоской стенке температура меняется по линейному закону и не зависит от свойств материала.

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

1. X=0, t=tc1=C2;

- 18 -

2. X= ,t=tc2=C1∙ +tс1.Отсюда С

 

 

 

tc 2

tc 1

.

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Запишем закон Фурье:

 

 

 

 

 

 

 

 

 

 

 

q

t

 

tc 2 tc 1

 

 

 

t

 

t

 

 

,

Вт

,

x

 

 

 

 

 

c 1

c 2

м2

 

 

 

 

 

 

 

 

 

 

где – тепловая проводимость плоской однослойной стенки. Отсюда

плотность теплового потока:

q tc 1 tc 2 ,

где – термическое сопротивление теплопроводности плоской одно-

слойной стенки.

Если известна площадь поверхности, через которую передается теплота, то

Q = qF= tc 1 tc 2 , Вт.

Тепловой поток пропорционален температурному напору на поверхности плоской стенки tc1 – tc2 и обратно пропорционален термическому сопротивлению.

1.4.2. Теплопроводность плоской многослойной стенки

Допустим, что плоская стенка состоит из n разнородных слоев толщиной 1, 2n с коэффициентами теплопроводности λ1, λ2…. λn (рис. 1.4). Слои такой стенки плотно прилегают друг к другу, поэтому температура на соприкасающихся поверхностях одна и та же.

При стационарном режиме в стенке установились во времени температуры tc1,tc2,tc3,tcn,tc(n+1). Тогда для первого слоя

- 19 -

q tc 1 tc 2 .

1

1

При стационарном режиме численно то же значение плотности теплового потока для второго слоя:

q tc 2 tc 3 .

2

 

2

Рис. 1.4

 

Докажем, что это так.

Для этого используем промежуточное выражение дифференциального уравнения теплопроводности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

Это объясняется тем, что при стационарном режиме

0,

1

0.

 

 

 

 

Сv

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

qx

0 , qx = const.Для n-го слоя значение плотности теплового

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

потока составит

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

tcn tc (n 1)

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

n

Записанные исходные соотношения переписываем каждое относительно температурного напора и складываем:

- 20 -

tc 1 tc (n 1)

q

 

1

 

 

 

2

 

.....

n .

 

 

 

 

1

 

 

 

2

 

 

 

n

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

tc 1 tc (n 1)

 

,

Вт

.

 

 

 

 

n

 

м2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

i 1 i

Последнее соотношение может быть переписано в виде

(1.13)

Заметим, что знаменатель представляет собой сумму термических сопротивлений многослойной плоской стенки. Его величина не меняется от перестановки слоев многослойной стенки, и поэтому не изменяется плотность теплового потока и сам тепловой поток.

1.4.3. Теплопроводность однослойной цилиндрической стенки

Задача решается в цилиндрических координатах (рис. 1.6). При этом уравнение температурного поля в общем виде запишется:

t=f (τ,r,φ,z)

(1.14)

Переход от полярных (рис. 1.5) к цилиндрическим координатам приводит к дифференциальному уравнению теплопроводности в общем виде

 

 

 

=a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сформулируем

 

задачу.

Рис. 1.5

 

 

 

 

 

 

 

 

 

 

 

 

 

2t

 

Пусть l >> r2

– r1, тогда

z

2 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- 21 -

Будем считать, что температура в цилиндрической стенке не меняется по окружности цилиндра, поэтому

2t

2 0 , а также

t

0 (так как ста-

ционарный режим). Тогда дифференциальное уравнение для сформулированной нами задачи запишется в виде

d 2t

1 dt

0.

 

 

 

Рис. 1.6

r 2

r dr

 

 

 

 

 

 

Введем новую переменную u:

 

 

 

 

 

dt

u ,

du

u

0.

 

 

 

dr

 

dr

r

 

Разделим переменные:

du

dr

0

 

 

 

 

 

u

r

 

 

 

и проинтегрируем: lnu+lnr = const = С1. Перепишем последнее выражение в виде u∙r = С1'. Подставим значение u:

dt

r

С1

', dt

C1 '

.

 

dr

 

 

 

 

dr

r

 

 

 

 

 

 

 

dr

Разделим переменные: dt С1 '

r .

Проинтегрируем:

t=С1'∙lnr + С2.

- 22 -

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]