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

Теплопередача учебное пособие

..pdf
Скачиваний:
38
Добавлен:
15.11.2022
Размер:
4.77 Mб
Скачать

ривается изучаемый процесс, с математической точки зрения являются величинами бесконечно малыми, а с физической точки зрения– величинами еще достаточно большими, чтобы в их пределах можно было игнорировать дискретным строением среды ирассматривать ее как сплошную. Полученная таким образом зависимость является общим дифференциальным уравнением рассматриваемого процесса.

Для облегчения вывода этого уравнения сделаем следующие допущения:

1)тело однородно и изотропно;

2)физические параметры постоянны;

3)деформация рассматриваемого объема, связанная с изменением температуры, является очень малой величиной по сравнению с самим объемом;

4)внутренние источники теплоты в теле распределены равномерно.

Воснову вывода дифференциального уравнения теплопроводности положен закон сохранения энергии, который в рассматриваемом случае может быть сформулирован следующим обра-

зом: количество теплоты dQ1 , введенное в элементарный объем dV извне за время dτ вследствие процессов теплопроводности, а также количество теплоты dQ2 , полученное от внутренних ис-

точников теплоты, равно изменению внутренней энергии вещества dQ, содержащегося в элементарном объеме за время dτ:

dQ1 + dQ2 = dQ.

(1.6)

Для нахождения составляющих уравнения (1.6) выделим в

теле элементарный параллелепипед со сторонами dx,

dy и dz

(рис. 1.3).

Количество теплоты, которое подводится к граням элементарного объема за время dτ в направлении осей 0x, 0 y, 0z,

обозначим соответственно 0 y, dQy , dQz .

11

z

dQz+dz

 

 

dQy

 

 

dQx

 

dQx+dx

dz

 

dy

 

 

 

dx

 

dQy+dy

dQz

x

0

 

 

y

Рис. 1.3. К выводу дифференциального уравнения теплопроводности

Количество теплоты, которое будет отводиться через противоположные грани в тех же направлениях, обозначим dQx+dx ,

dQy +dy , dQz +dz . Количество теплоты, подведенное к грани dy dz в направлении оси 0x за время dτ, составляет

dQx = qx dy dz dτ,

где qx – проекция плотности теплового потока на направление

нормали к указанной грани.

Количество теплоты, отведенное через противоположную грань элементарного параллелепипеда в направлении оси 0x, запишетсяввиде

dQx+dx = qx+dx dydzdτ.

Тогда

dQx1 = dQx dQx+dx ,

12

или

dQx1 = qx dydzdτqx+dx dydzdτ .

(а)

Функция qx+dx является непрерывной в рассматриваемом интервале dx и может быть разложена в ряд Тейлора:

 

 

 

qx

 

2 qx

 

dx2

qx+dx

= qx

+

 

dx +

 

 

 

 

 

 

+ ...

x

 

x

2

2!

 

 

 

 

 

 

 

Если ограничиться двумя первыми членами ряда, то уравнение (а) запишется в виде

dQx1

= −

qx

dxdydzdτ.

(б)

 

 

 

x

 

Аналогичным образом можно найти количество теплоты, подведенное к элементарному объему и в направлениях двух других координатных осей 0 y и 0z.

Количество теплоты dQ1 , подведенное теплопроводностью

к рассматриваемому объему извне,

 

 

qx

 

qy

 

 

qz

dQ1

= −

 

+

 

+

 

 

x

y

z

 

 

 

 

dx dy dz dτ. (в)

Определим вторую составляющую уравнения (1.6). Обозначим количество теплоты, выделяемое внутренними источниками в единице объема среды в единицу времени и называемое мощностьювнутренних источников теплоты, через qV , Вт/м3, тогда

dQ2 = qV dVdτ. (г)

Третья составляющая в уравнении (1.6) найдется в зависимости от характера термодинамического процесса изменения системы.

13

При изохорном процессе вся теплота, подведенная к элементарному объему, уйдет на изменение внутренней энергии вещества, заключенного в этом объеме, т.е.

dQ = dU = cρ

t

dτdV ,

(д)

 

 

τ

 

где c – удельная теплоемкость, Дж/ (кг °С); ρ – плотность ве-

щества, кг/м3.

Подставляя полученные выражения (в), (г) и (д) в уравнение (1.6), получаем

 

t

 

qx

 

 

qy

 

qz

 

 

 

cρ

 

= −

 

 

+

 

 

+

 

 

+ qV .

(1.7)

τ

x

 

z

 

 

 

y

 

 

 

Выражение (1.7) является дифференциальным уравнением энергии для изохорного процесса переноса теплоты.

В твердых телах перенос теплоты осуществляется по закону Фурье (1.3):

q= −λ t n .

n

Проекции вектора плотности теплового потока на координатные оси 0x, 0 y, 0z определяются выражениями:

qx

= −λ

t

;

qy

= −λ

t

; qz

= −λ

t

.

 

 

 

 

 

x

 

 

y

 

z

Подставляя полученные выражения проекций вектора плотноститепловогопотока в уравнение (1.7), получим

 

t

 

 

t

 

 

 

 

t

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

cρ

 

 

=

 

 

λ

 

 

+

 

 

λ

 

 

+

 

 

λ

 

 

+ qV .

(1.8)

τ

 

 

 

 

 

 

 

x

x

y

 

y

 

z

 

z

 

 

Выражение (1.8) называется дифференциальным уравнением теплопроводности. Оно устанавливает связь между времен-

14

ным и пространственными изменениями температуры в любой точке тела, в которых происходит процесс теплопроводности.

Если принять, что теплофизические характеристики постоянны, то

 

 

 

 

 

t

 

 

 

 

λ

 

 

 

2t

 

 

2t

 

2t

 

 

 

qV

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

+

 

 

 

 

 

 

+

 

 

 

.

(1.9)

 

 

 

 

 

∂ τ

 

 

 

 

 

 

x

2

 

 

y

2

 

z

2

cρ

 

 

 

 

 

 

 

 

cρ ∂

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В уравнении (1.9) можно обозначить

 

 

 

 

 

 

 

 

 

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

2t

 

 

 

2t

 

2t

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= а

,

 

 

 

 

 

 

 

+

 

 

 

 

 

+

 

 

 

 

 

=

 

t

 

,

 

 

 

 

 

 

 

сρ

 

 

x2

y2

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

=

x2+

y2+

 

 

 

 

 

– выражение оператора Лапласа в декар-

 

 

z2

 

товой системе координат;

 

а – коэффициент температуропровод-

ности, являющийся мерой тепловой инерции вещества, м2 с.

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

= a

2t+

 

qV

.

 

 

 

 

 

 

 

 

(1.10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cρ

 

 

 

 

 

 

 

 

 

 

 

Выражение

 

 

 

2t

 

в

 

 

цилиндрической

системе

координат

имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2t

 

 

 

1∂ t

 

 

 

1 ∂

 

2t

 

2t

 

 

 

 

 

 

 

 

 

 

t=

 

 

 

 

 

+

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r2

 

 

 

 

r r 2

φ2

 

 

 

где r

– радиус-вектор; φ – полярный угол;

z

 

– аппликата.

 

Если система тел не содержит внутренних источников тепла

( qV = 0 ), тогдауравнение (1.10) принимает вид уравнения Фурье:

t

= a 2t.

(1.11)

∂ τ

 

 

 

 

15

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

2t

 

2t

2t

 

qV

 

 

 

+

 

 

 

+

 

+

 

= 0.

(1.12)

x2

 

z2

λ

 

y2

 

 

 

При стационарной теплопроводности и отсутствии внутренних источников теплоты уравнение (1.9) примет вид уравнения Лапласа:

2t

 

2t

2t

 

 

 

+

 

 

 

+

 

= 0.

(1.13)

x2

y2

z2

 

 

 

1.3. УСЛОВИЯ ОДНОЗНАЧНОСТИ ДЛЯ ПРОЦЕССОВ ТЕПЛОПРОВОДНОСТИ

Полученное дифференциальное уравнение теплопроводности описывает целый класс соответствующих явлений. Чтобы выделить конкретно рассматриваемый процесс и дать его полное математическое описание, к дифференциальному уравнению необходимо присоединить математическое описание всех частных особенностей рассматриваемого процесса.

Эти частные особенности, которые совместно с дифференциальным уравнением дают полное математическое описание конкретного процесса теплопроводности, называются условиями однозначности или краевыми условиями [1, 2].

Условия однозначности включают в себя:

1)геометрические условия, характеризующие форму и размеры тела, в котором протекает процесс;

2)физические условия, характеризующие физические свойства среды;

16

3)временные (начальные) условия, характеризующие распределение температуры в изучаемом теле в начальный момент времени;

4)граничные условия, характеризующие взаимодействие рассматриваемого тела с окружающей средой.

Существуют граничные условия первого, второго, третьего и четвертого рода.

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

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

tc = f (xс, yс, zс, τ) ,

(1.14)

где tс – температура на поверхности тела;

xс, yс, zс – координа-

ты поверхности тела.

 

Простейший случай tс = const .

 

Граничное условие второго рода. При этом задается зна-

чение теплового потока для каждой точки поверхности тела и любого момента времени:

qс = f (xс, yс, zс, τ),

(1.15)

где qс – плотность теплового потока на поверхности тела. Простейший случай qс = const .

Граничное условие третьего рода. При этом задается температура окружающей среды tж и закон теплообмена меж-

ду поверхностью тела и окружающей средой. Для описания процесса теплообмена между поверхностью тела и средой используется закон Ньютона–Рихмана.

Согласно этому закону количество теплоты, отдаваемое единицей поверхности тела в единицу времени, пропорционально разности температур поверхности тела tс и окружаю-

щей среды tж ( tс > tж ):

17

q = α(tс tж ),

(1.16)

где α – коэффициент пропорциональности, называемый коэффициентом теплоотдачи, Вт/(м2·°C).

Коэффициент теплоотдачи характеризует интенсивность теплообмена между поверхностью тела и окружающей средой. Численно он равен количеству теплоты, отдаваемому (или воспринимаемому) единицей поверхности в единицу времени при разности температур между поверхностью тела и окружающей средой, равной одному градусу.

Согласно закону сохранения энергии количество теплоты, которое отводится с единицы поверхности в единицу времени вследствие теплоотдачи, должно равняться теплоте, подводимой к единице поверхности в единицу времени вследствие теплопроводности из внутренних объемов тела, т.е.

α(tс tж ) = −λ

t

,

(1.17)

 

 

n с

 

где n – нормаль к поверхности тела; индекс «с» указывает на то, что температура и градиент относятся к поверхности тела

(при n = 0 ).

Окончательно граничное условие третьего рода можно записать в виде

 

t

= −

α

(tс tж ) .

 

 

 

 

 

(1.18)

 

λ

 

n с

 

 

 

Граничное условие четвертого рода. Это условие харак-

теризует условие теплообмена системы тел или тела с окружающей средой по закону теплопроводности.

В рассматриваемых условиях имеет место идеальный контакт (температуры соприкасающихся поверхностей одинаковы) и равенство тепловых потоков, проходящих через поверхность соприкосновения (рис. 1.4):

18

 

t

 

 

t

2

 

 

1

= λ

 

 

 

 

 

λ1

 

 

2

 

 

 

 

 

,

 

 

 

 

n с

 

 

n с

(1.19)

(t1 )c = (t2 )c .

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

t1

 

t

 

 

 

 

ϕ

1

ϕ 2

 

С

 

 

)

 

 

t2

2

 

 

t

 

 

=(

 

 

 

С

 

 

 

)

 

 

 

1

 

 

 

(t

 

 

 

 

1

2

n

Рис. 1.4. Граничное условиечетвертого рода

Дифференциальное уравнение теплопроводности совместно с условиями однозначности дает полную математическую формулировку конкретной тепловой задачи. Поставленная таким образом задача может быть решена аналитическим, численным или экспериментальным методами.

1.4. ТЕПЛОПРОВОДНОСТЬПРИСТАЦИОНАРНОМРЕЖИМЕ

При установившемся или стационарном тепловом режиме

температура тела во времени остается постоянной, т.е. ∂ t

τ= 0.

При этом дифференциальное уравнение теплопроводности

(1.10) будет иметь вид

 

λ 2t+ qV= 0 .

(1.20)

 

19

Если внутренниеисточники теплоты отсутствуют, qV = 0 то

2t= 0 .

(1.21)

1.4.1. Передача теплоты через плоскую стенку (qV

= 0)

Рассмотрим однородную стенку толщиной δ (рис. 1.5), коэффициент теплопроводности λ постоянен, на наружных поверхностях стенки поддерживаются постоянные темпера-

туры t1 и t2 [1, 2].

t

t1

t2

x

0

δ

 

Рис. 1.5. Однородная плоская стенка

При заданных условиях температура будет изменяться только в направлении, перпендикулярном плоскости стенки, т.е. в направлении оси x. А в направлении 0 y и 0z температу-

рабудет постоянной.

t

=

t

= 0 .

(а)

 

 

 

y

z

 

Таким образом, уравнение теплопроводности (1.21) примет вид

2

d t = 0 . (1.22) dx2

Граничные условия будут следующими:

при

х = 0

t = t1;

 

при

х = δ

t = t

2

.

(б)

 

 

 

 

 

Уравнение (а) и условия (б) дают полную математическую формулировку рассматриваемой задачи.

Закон распределения температуры по толщине найдется в результате двойного интегрирования уравнения (1.22).

20