Лекции теплотехника часть1
.pdfдеформация рассматриваемого объема, связанная с изменением температуры, является очень малой величиной по сравнению с самим объемом;
внутренние источники теплоты в теле, которые в об-
щем случае могут быть заданы как qυ = f (x,у,z,τ), распределены равномерно. (В ряде случаев внутри объектов исследования могут протекать процессы, в результате которых будет выделяться или поглощаться теплота. Например, выделение джоулевой теплоты при прохождении электрического тока по проводникам; выделение или поглощение теплоты при протекании ряда химических реакций и т.д. Дей-
ствие внутренних источников количественно характеризуется их мощностью qυ, Вт/м3, т.е. количеством теплоты, которое выделяется в единице объема среды в единицу времени).
В основу вывода дифференциального уравнения теплопроводности положен закон сохранения энергии, который в рассматриваемом случае может быть сформулирован следующим образом: количество теплоты dQ, введенное в элементарный объем извне за время dτ вследствие теплопроводности, а также от внутренних источников, равно изменению внутренней энергии вещества содержащегося в элементарном объеме:
dQ1 dQ2 |
dU , |
(1.12) |
где dQ1 – количество теплоты, введенное в элементарный объем путем теплопроводности за время dτ, Дж;
dQ2 – количество теплоты, которое за время dτ выделилось в элементарном объеме dυ за счет внутренних источников;
dU – изменение внутренней энергии вещества, содержащегося в элементарном объеме dυ, за время dτ.
Выделим в теле элементарный параллелепипед со сторонами dx, dy, dz (рис. 1.3), причем параллелепипед расположим так, чтобы его грани были параллельны соответствующим координатным плоскостям.
11
Количество теплоты, которое подводится к граням элементарного объема за время dτ в направлении осей Ох, Оу, Oz, обозначим соответственно dQx, dQy, dQz.
Рис. 1.3 К выводу дифференциального уравнения теплопроводности
Количество теплоты, которое будет отводиться через противоположные грани в тех же направлениях, обозначим
соответственно dQx+dx, dQy+dy, dQz+dz.
Количество теплоты, подведенное к грани dydz в направлении оси Ох за время dτ, согласно закону Фурье составляет
dQx dydzd t x , |
(1.13) |
Количество теплоты, отведенное через противоположную грань элементарного параллелепипеда в направлении оси Ох, за тот же промежуток времени запишется как
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
dQx dx |
dydzd |
|
t |
|
dx , |
(1.14) |
|||
|
|
|
x |
||||||||
|
|
t |
|
|
x |
|
|
|
|||
|
|
|
|
|
|
|
|
||||
где |
t |
|
dx |
– температура второй |
грани, |
а величина |
|||||
x |
|||||||||||
|
|
|
|
|
|
|
|
|
|
||
xt dx определяет изменение температуры в направление оси Ox. Уравнение (1.14) можно представить в следующем виде:
12
dQx dx dydzd |
t |
dxdydzd |
2t |
. |
(1.15) |
|
x |
x2 |
|||||
|
|
|
|
Разница между количеством теплоты, подведенного к элементарному параллелепипеду и отведенного от него за время dτ в направлении оси Ох
dQx1 dQx dQx dx dydzd |
t |
|
dydzd |
t |
|
||||||
x |
x |
||||||||||
|
2t |
|
|
2t |
|
(1.16) |
|||||
dxdydzd |
dxdydzd |
|
|
|
|||||||
|
. |
|
|
||||||||
x |
2 |
|
x2 |
|
|
||||||
|
|
|
|
|
|
|
|
||||
Аналогичным образом можно найти количество теплоты, подводимое к элементарному объему и в направлениях двух других координатных осей Оу и Oz:
dQy1 |
dxdydzd |
2t |
; dQz1 dxdydzd |
2t |
. |
(1.17) |
||
y |
2 |
z 2 |
||||||
|
|
|
|
|
||||
Количество теплоты, подведенное в результате теплопроводности к рассматриваемому объему, будет равно:
dQ1 dQx1 dQy1 dQz1
|
|
2 |
t |
|
|
2 |
t |
|
|
2 |
t |
|
(1.18) |
|
|
|
|
|
|
|
|
|
|||||||
dxdydzd |
x |
2 |
y |
2 |
z |
2 |
. |
|
||||||
|
|
|
|
|
|
|
|
|
||||||
Определим вторую составляющую уравнения (1.12). Учитывая, что мощность внутренних источников равна qυ, количество теплоты, выделенное в элементарном объеме dυ за время dτ будет составлять
dQ 2 q d d q dxdydzd |
. |
(1.19) |
Изменение внутренней энергии выделенного элемента, массой ρdυ, при изменении его температуры на величину dt выразится как
dU c d dt c dxdydz |
t |
d , |
(1.20) |
|
|
||||
|
|
|
где с – массовая теплоемкость вещества, Дж/(кг·К); ρ – плотность вещества, кг/м3;
13
dt |
t |
|
d . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Подставляя полученные выражения (1.18), (1.19) и (1.20) в |
||||||||||||||||||||||||||||||||||||||
уравнение (1.12), получим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
2t |
|
|
2t |
|
|
|
2t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
||||||||||||
dxdydzd |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q |
dxdydzd c dxdydz |
|
d . |
(1.21) |
||||||||||||||
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|||||||||||||||||||||||
|
|
x |
|
|
|
y |
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Сократив все члены равенства (1.23) на dxdydzd |
и разделив |
|||||||||||||||||||||||||||||||||||||
на сρ получим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2t |
|
|
|
|
|
||||||||
|
|
|
|
t |
|
|
|
|
|
|
2t |
|
|
|
2t |
|
|
|
q |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
(1.22) |
||
|
|
|
|
|
|
|
c |
|
x |
2 |
|
|
y |
2 |
|
|
z |
2 |
|
c |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|||||||
Величину |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
называют оператором Лапласа |
|||||||||||||||
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|
|
2 |
||||||||||||||||||||||||
|
|
|
|
|
x |
|
|
|
|
y |
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
и обозначают 2 (знак читается “набла”); величину c – ко-
эффициентом температуропроводности и обозначают буквой
а, м2/с. При указанных обозначениях дифференциальное уравнение теплопроводности принимает вид
t |
a 2t |
q |
. |
(1.23) |
|
|
|||
|
c |
|
||
При отсутствии внутренних источников теплоты (qυ = 0) дифференциальное уравнение теплопроводности примет вид
t |
a 2t. |
(1.24) |
|
|
|||
|
|
Дифференциальное уравнение теплопроводности выражает связь между изменением температуры во времени и ее распределением в пространстве. Действительно, левая часть уравнения (1.24) характеризует скорость изменения температуры некоторой точки тела во времени, правая – пространственное распределение температуры вблизи этой точки.
Дифференциальное уравнение теплопроводности в цилиндрической системе координат r-φ-z при наличии внутренних источников теплоты имеет вид:
14
t |
|
2t |
1 t |
|
1 2t |
2t |
|
q |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
r |
2 |
r |
r |
|
r |
2 |
|
|
2 |
z |
2 |
|
|
c |
. |
(1.25) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
где r – радиус-вектор; φ – полярный угол;
z – аппликата.
Коэффициент температуропроводности а является физическим параметром вещества. Он существен для нестационарных тепловых процессов и характеризует скорость изменения температуры. Если коэффициент теплопроводности характеризует способность тел проводить теплоту, то коэффициент температуропроводности является мерой теплоинерционных свойств тела. Из уравнения (1.23) следует, что изменение температуры во
времени t
для любой точки пространства пропорционально
величине а. Иначе говоря, скорость изменения температуры в любой точке тела будет тем больше, чем больше коэффициент температуропроводности. Поэтому при прочих равных условиях выравнивание температур во всех точках пространства будет происходить быстрее в том теле, которое обладает большим коэффициентом температуропроводности. Коэффициент температуропроводности зависит от природы вещества. Например, жидкости и газы обладают большой тепловой инерционностью и, следовательно, малым коэффициентом температуропроводности. Металлы обладают малой тепловой инерционностью, так как они имеют большой коэффициент температуропроводности.
2 Условия однозначности
Так как дифференциальное уравнение теплопроводности выведено на основе общих законов физики, то оно описывает явление теплопроводности в самом общем виде. Поэтому можно сказать, что полученное дифференциальное уравнение описывает целый класс явлений теплопроводности. Чтобы из бесчисленного количества выделить конкретно рассматриваемый процесс и дать его полное математическое описание, к дифференциальному уравнению необходимо присоединить математическое описание всех частных особенностей рассматриваемого процес-
15
са. Эти частные особенности, которые совместно с дифференциальным уравнением дают полное математическое описание конкретного процесса теплопроводности, называются условиями однозначности или краевыми условиями.
Условия однозначности включают в себя:
геометрические условия, характеризующие форму и размеры тел, в которых протекает процесс;
физические условия, характеризующие физические свойства среды и тела;
начальные условия, характеризующие распределение температур в изучаемом теле в начальный момент времени;
граничные условия, характеризующие взаимодействие рассматриваемого тела с окружающей средой.
Геометрическими условиями задаются форма и линейные размеры тела, в котором протекает процесс.
Физическими условиями задаются физические параметры тела (λ, с, ρ и др.) и может быть задан закон распределения внутренних источников теплоты.
Начальные условия необходимы при рассмотрении нестационарных процессов и состоят в задании закона распределения температуры внутри тела в начальный момент времени. В общем случае начальное условие аналитически может быть записано следующим образом (при τ = 0):
t f x, y,z . |
(1.26) |
В случае равномерного распределения температуры в теле начальное условие упрощается (при τ = 0):
t to const . |
(1.27) |
Граничные условия могут быть заданы несколькими способами.
Граничные условия первого рода. При этом задается распределение температуры на поверхности тела для каждого момента времени:
tп f x, y,z, . |
(1.28) |
где tп – температура на поверхности тела; х, у, z – координаты поверхности тела.
16
В частном случае, когда температура на поверхности является постоянной на протяжении всего времени протекания процессов теплообмена, уравнение (1.28) упрощается и принимает вид:
tп const . |
(1.29) |
Граничные условия второго рода. При этом задаются значения теплового потока для каждой точки поверхности тела и любого момента времени.
Аналитически это можно представить следующим образом:
qп f x, y,z, . |
(1.30) |
где qп – плотность теплового потока на поверхности тела; х, у, z – координаты поверхности тела.
В простейшем случае плотность теплового потока по поверхности и во времени остается постоянной:
qп qo const . |
(1.31) |
Граничные условия третьего рода. При этом задаются температура окружающей среды tо.с. и закон теплообмена между поверхностью тела и окружающей средой. Граничное условие третьего рода характеризует закон теплообмена между поверхностью и окружающей средой в процессе охлаждения и нагревания тела.
Например, если теплообмен между поверхностью тела и окружающей средой осуществляется конвективной теплоотдачей, плотность теплового потока, отводимого от поверхности тела (или подводимого к ней) согласно закону Ньютона-Рихмана
q α tп tо.с. ,. |
(1.32) |
где tп – температура поверхности тела, оС ;
tо.с. – температура окружающей среды, оС ;
α – коэффициент пропорциональности, называемый коэффициентом теплоотдачи, Вт/(м2ּК).
Коэффициент теплоотдачи характеризует интенсивность теплообмена между поверхностью тела и окружающей средой. Он численно равен количеству теплоты, отдаваемой (или воспринимаемой) единицей поверхности в единицу времени при разно-
17
сти температур между поверхностью тела и окружающей средой, равной одному градусу.
Согласно закону сохранения энергии, количество теплоты, которое отводится с единицы поверхности тела в единицу времени вследствие теплоотдачи, должно равняться количеству теплоты, подводимому к единице поверхности в единицу времени вследствие теплопроводности из внутренних объемов тела, т.е.
|
t |
|
|
|
α tп tо.с. λ |
|
|
, |
(1.33) |
|
||||
|
n |
пов |
|
|
где индекс “пов” показывает, на то, что температура и градиент относятся к поверхности тела.
Равенство (1.33) является математической формулировкой граничного условия третьего рода и является действительным для каждого момента времени.
Следует заметить, что в общем случае коэффициент теплоотдачи α и температура окружающей среды tо.с. могут быть переменными, но обязательно заданными функциями времени и координат.
Дифференциальное уравнение (1.23) совместно с условиями однозначности дают полную математическую формулировку конкретной задачи теплопроводности. Поставленная таким образом задача разрешается аналитическим, численным или экспериментальным методом.
18
Лекция № 3
ТЕПЛОПРОВОДНОСТЬ ЧЕРЕЗ ПЛОСКУЮ СТЕНКУ ПРИ СТАЦИОНАРНОМ РЕЖИМЕ
При установившемся, или стационарном, тепловом режиме температура тела во времени остается постоянной, т.е.t / 0 . В этом случае при отсутствии внутренних источников теплоты (qυ = 0) дифференциальное уравнение теплопроводности будет иметь вид:
или |
|
2t 0 |
|
(2.1) |
|||
|
|
|
|
|
|
||
|
2t |
|
2t |
|
2t |
0. |
(2.2) |
|
x2 |
y2 |
z2 |
||||
|
|
|
|
|
|||
Рассмотрим однородную изотропную стенку, высота и ширина которой бесконечно велики по сравнению с толщиной. В этом случае температура будет изменяться только в направлении, перпендикулярном плоскости стенки (т.е. по толщине), значит,
t |
|
t |
0, |
(2.3) |
|
y |
z |
||||
|
|
|
а дифференциальное уравнение теплопроводности примет вид
2t |
0. |
(2.4) |
|
x2 |
|||
|
|
Граничные условия первого рода
Пусть толщина стенки равна δ, коэффициент теплопроводности материала стенки постоянен и равен λ. На наружных поверхностях стенки поддерживаются постоянными температуры tc1 и tc2 (рис. 2.1).
Граничные условия зададим следующим образом: |
|
|||
при |
х = 0 |
t = tc1, |
(2.5) |
|
при |
х = δ |
t = tc2. |
||
|
||||
19
Уравнение (2.4) и условия (2.5) дают полную математическую формулировку рассматриваемой задачи.
Рис. 2.1 Теплопроводность однородной плоской стенки
В результате решения поставленной задачи должно быть найдено распределение температуры в плоской стенке, т. е. функция t = f(x), и получена формула для определения плотности теплового потока через стенку.
Закон распределения температур по толщине стенки найдем в результате двойного интегрирования уравнения (2.4).
Первое интегрирование дает:
dt |
C , |
(2.6) |
|
||
dx |
1 |
|
|
|
После второго интегрирования получим:
t = C1x + C2. (2.7) Из уравнения (2.7) следует, что при постоянном коэффици-
енте теплопроводности температура в стенке изменяется по линейному закону (рис. 2.2).
Постоянные интегрирования C1 и С2 в уравнении (2.7) определим из граничных условий:
при |
х = 0 |
t = tc1 |
и |
С2 |
= tc1, |
(2.8) |
|
при |
х = δ |
t = tc2 |
и |
С1 = - (tc1 - tc2) / δ. |
|||
|
|||||||
20
