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

книги / Теоретические основы теплотехники. Теория теплообмена

.pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
2.94 Mб
Скачать

 

 

 

t

 

 

 

 

t

 

 

grad t x

 

 

 

cos n,x

 

 

 

 

;

 

n

x

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

t

 

 

grad t y

 

 

 

cos n, y

 

 

 

;

(1.2)

n

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

t

 

 

grad t z

 

 

 

 

cos n,z

 

 

 

 

.

 

 

n

 

z

 

 

 

 

 

 

 

 

 

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

Тепловым потоком Q называется количество теплоты W [Дж], проходящее через некоторую контрольную поверхность площадьюF вединицувремениτ[с]. ВсистемеСИтепловойпоток измеряется в ваттах:

Q

W

Дж

 

(1.3)

 

 

с

= Вт.

 

 

 

 

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

Плотность теплового потока q – это тепловой поток, прохо-

дящий через элементарную площадку dF изотермической поверх-

ности в данной ее точке:

 

dQ

Вт

(1.4)

q

dF

 

м

2 .

 

 

 

 

Плотность теплового потока является векторной величиной, направленной в сторону, противоположную направлению градиента температуры, т.е. в сторону уменьшения температуры.

ТепловойпотокQ черезизотермическуюповерхностьиплот-

ность теплового потока q связаны очевидным соотношением

Q q dF .

(1.5)

F

 

10

1.2. Основной закон теплопроводности. Коэффициент теплопроводности

Основным законом теплопроводности является закон Фурье. Он устанавливает связь между неоднородностью температурного поля, обусловленную величиной температурных градиентов, с возникающими тепловыми потоками. Согласно гипотезе Фурье, количество теплоты δW, проходящее через элемент изотермической поверхности dF за промежуток времени , пропорционально температурному градиенту:

W dF d grad t Дж .

(1.6)

Для теплового потока и плотности теплового потока выраже-

ние закона Фурье будет иметь вид

 

 

 

 

 

dQ dF

 

grad t

 

Вт ;

 

 

 

 

 

 

 

 

 

 

2

 

(1.7)

 

 

 

q grad t

 

 

 

Вт/ м

 

.

 

Минус в правой части математической записи закона учитывает разнонаправленность векторов теплового потока и градиента температуры.

Коэффициентом пропорциональности в данном уравнении служит коэффициент теплопроводности – это важный тепло-

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

ной величине температурного градиента:

 

 

 

 

 

 

 

 

q

 

 

Вт / м К .

(1.8)

 

 

 

grad t

 

 

 

 

Коэффициент теплопроводности определяется экспериментально и зависит от множества факторов:

температуры;

давления;

11

влажности;

природы тела (однородное или неоднородное, изотропное или анизотропное и т.д.).

Большинство методов экспериментального определения теплопроводности основано на измерении теплового потока и градиентатемпературвисследуемомвеществе. Коэффициенттеплопроводности при этом определяется из соотношения (1.8).

Ниже приведен порядок значений коэффициентов теплопроводности различных материалов, Вт/(м К):

0,05–0,1 – строительные теплоизоляционные материалы;

1–10 – пластмассы;

20–50 – железо, стали;

393 – медь;

410 – серебро.

Материалы с низким значением коэффициента теплопроводности [меньше 0, 25 Вт/(м К)], применяемые для тепловой изоля-

ции, называются теплоизоляционными.

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

 

 

,

(1.9)

0 1 b t t0

 

где 0 – значение коэффициента теплопроводности материала при температуре t0 ; b – постоянная, определяемая при обработке дан-

ных экспериментального исследования теплопроводности для данного материала.

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

12

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

Так, например, для приближённого нахождения коэффициента теплопроводности чистого (однокомпонентного) газа можно использовать следующую формулу:

0,25 9k 5 cv ,

(1.10)

где k – показатель адиабаты; μ – коэффициент динамической вязкости, Па с; – изохорная теплоемкость, Дж/(кг К).

В формуле (1.10) учитывается природа, структура и количество степеней свободы молекул газа, а также его реологические характеристики.

Теплопроводность газов изменяется в пределах 0,006...

...0,1 Вт/(м К) и с ростом температуры увеличивается.

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

Несколькопрощетеплопроводностьопределяетсядляжидкостей, так как их вязкость существенно выше, и быстрое развитие конвекции затруднено. Коэффициент теплопроводности жидкостей лежит в пределах от 0,07 до 0,7 Вт/(м К). Опыты показывают, что для большинства жидкостей с повышением температуры коэффициент теплопроводности убывает, исключение составляют вода и глицерин.

Для изотропной однородной среды коэффициент теплопроводности const и одинаков по направлениям осей координат, а в неоднородной среде в зависимости от температуры его

13

значение x, y,z может быть принято в соответствии с фор-

мулой (1.9). В большинстве практических случаев значение коэффициентатеплопроводности такжеможет быть принято средним в рассматриваемом диапазоне температур.

Для кристаллов и некоторых искусственных материалов, например композитов, теплопроводность зависит от направления

(анизотропия теплопроводности). Для таких материалов коэффи-

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

xx xy xz

 

 

 

 

(1.11)

yx yy yz .

 

 

 

zx zy zz

 

 

Физически здесь λij означает теплопроводность среды в направленииi приградиентетемпературывнаправлении j. Тензор (1.11) является симметричным тензором, то есть выполняются условия i j j i .

Для ортотропных тел, имеющих различную теплопроводность в трех взаимно перпендикулярных направлениях, тензор теплопроводности при соответствующем выборе направлений осей координат имеет шаровой вид:

 

xx 0 0

 

 

 

 

yy 0

 

(1.12)

0

.

 

0

 

 

 

 

0 zz

 

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

14

влаги внутри пористого материала. Другой пример – рост теплопроводности пористого материала при повышении температуры по причине роста теплопроводности газа, заполняющего поры. Большоевлияниенатеплопроводностьоказываетиплотностьданных тел. Чем она выше, тем меньше поры, а значит, больше основного материала, из которого состоит пористое тело. А теплопроводность твердых тел, как правило, в разы больше теплопроводности газов.

1.3. Общая постановка задачи теплопроводности. Дифференциальное уравнение Фурье

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

определения зависимости q q x, y,z, , необходимо знатьфунк-

цию температурного поля в области:

 

t f ( x, y,z, ).

(1.13)

Дляполученияфункциональнойзависимости(1.13) используется общий метод математической физики, когда из всего пространства рассматривается лишь элементарный объем, все процессы в котором происходят за малый отрезок времени. Полученная таким образом зависимость является общим дифференциальным уравнением исследуемого процесса. Интегрируя дифференциальное уравнение, можно получитьаналитическую зависимость между величинами для всей расчетной области в течение всего времени протекания процесса. Функция температурного поля (1.13) является интегралом решения дифференциального уравне-

ния теплопроводности.

15

Впервые вывод дифференциального уравнения теплопроводности был предложен М.В. Остроградским в 1830 г. В основу вывода положен закон сохранения энергии в следующей формулировке: количество теплоты dW, полученное элементарным объемом извне за время в процессе теплопроводности (а также от внутренних источников теплоты), равно изменению внутренней энергии или энтальпии объема (соответственно для изохорного или изобарного процессов теплообмена).

При выводе приняты следующие допущения:

тело сплошное, однородное, изотропное;

физические параметры (свойства) тела постоянны;

в теле нет внутренних источников и стоков теплоты;

тело неподвижное, перенос теплоты осуществляется только посредством теплопроводности;

деформация объема, связанная с изменением температуры, пренебрежимо мала.

t

t + t

Рис. 1.3. Расчётная схема

16

Закон сохранения энергии в процессе передачи тепла теплопроводностью сформулируем для элементарного объема тела (параллелепипеда) с размерами dx,dy,dz (рис.1.3).

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

Рассмотрим процесс передачи тепла в направлении оси . В соответствии с законом Фурье теплота, подведенная к объему за время , запишется как

Wx λdy dz dτ

t

,

(1.14)

x

 

 

 

где – площадь грани.

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

Wx Wx λdy dz dτ

 

t t .

(1.15)

x

 

 

 

Примем, что температура в пределах элементарного объема изменяется линейно от температуры t на левой грани до темпера-

турыt + t направойграни. Тогда t

 

t

d x , авыражение(1.15)

x

 

 

 

 

 

 

будет иметь вид

 

 

 

 

 

 

 

 

t

 

2

t

 

 

Wx Wx λdy dz dτ

 

 

d x .

(1.16)

 

 

2

 

x

x

 

 

Количество теплоты Wx, затраченное на изменение энтальпии элементарного объема, равно разности между (1.14) и (1.16):

17

W

 

2t

dx dy dz d

2t

d dV ,

(1.17)

x2

x2

x

 

 

 

 

где dV = dxdydz – объем параллелепипеда.

Аналогичные соотношения можно получить для передачи теплоты по направлениям и .

Тогда суммарное количество теплоты W, затраченное на изменение энтальпии объема, будет равно сумме изменений по каждому направлению:

 

2

t

 

2

t

 

2

t

 

 

W Wx Wy Wz

 

 

 

 

 

dVd .

(1.18)

 

2

 

2

 

2

 

x

y

z

 

 

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

элементарного объема изменилась на величину

 

W dm cp t ,

(1.19)

где dm = ρ dV – масса элементарного объема; ср – изобарная теплоемкость.

Предполагая, что температура объема за промежуток времени

dτ изменялась линейно от температуры t

до температуры t + t,

можно записать t

t

d , а выражение (1.19) будет иметь вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W dV c

 

 

t

d .

(1.20)

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Приравнивая (1.18) и (1.20) и сокращая полученное равенство

на dV и dτ, получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

2

t

 

2

t

 

 

2

t

 

 

cp

 

 

 

 

 

 

 

 

 

2t.

(1.21)

 

 

2

 

2

 

 

 

2

 

 

 

x

y

 

 

 

z

 

 

Полученное уравнение называется дифференциальным уравнением теплопроводности. Размерность уравнения (1.21) Вт/м3.

18

Часто уравнение записывают в размерности К/с:

 

 

 

t

a 2t ,

(1.22)

 

 

 

 

 

 

 

 

 

где a

 

– коэффициент температуропроводности, м2/с. Ко-

 

 

cp

 

 

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

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

2t 0.

(1.23)

В случае, когда в исследуемой области действуют внутренние источники теплоты (например, джоулево тепло при пропускании электрического тока по проводнику), энтальпия элементарного объема изменяется как за счет проходящих через объем тепловых потоков (1.18), так и за счет теплоты источников q :

 

t

 

2

t

 

2

t

 

2

t

 

 

 

cp

 

 

 

 

 

 

 

q ,

(1.24)

 

 

2

 

2

 

2

 

 

x

 

y

 

z

 

 

 

где q – мощность внутренних источников теплоты, Вт/м3.

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

 

t

 

 

2

t

 

2

t

 

2

t

 

 

 

cp

 

x

 

y

 

z

 

 

q .

(1.25)

 

 

2

 

2

 

2

 

 

 

x

 

y

 

z

 

 

 

Проекции вектора плотности теплового потока на координатные оси Оx, Oy, Oz определяются выражениями

19