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

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

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

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

от свойств теплоносителя и его температуры;

от температурного напора;

от скорости движения и режима течения теплоносителя;

от вида конвекции;

от геометрии тела и условий обтекания.

Всвязи с этим коэффициент теплоотдачи в отличие, например, от коэффициента теплопроводности не является теплофизической характеристикой вещества, и его значения не приводятся

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

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

фициентов теплоотдачи.

2.1. Конвективный тепловой поток

Исходя из определения процесса теплоотдачи, тепловой поток, определяемыйпо законуНьютона– Рихмана, складываетсяиз двух слагаемых: конвективная составляющая потока qк итепловой поток за счет теплопроводности теплоносителя qт:

q qк qт .

(2.3)

60

Тепловой поток, обусловленный теплопроводностью, определяется в соответствии с законом Фурье:

qт f gradt,

(2.4)

где λf – теплопроводность теплоносителя.

Для определения конвективной составляющей потока рассмотрим простой пример. Пусть в канал поперечным сечением F со скоростью w втекает теплоноситель плотностью ρ, теплоемкостью ср, температурой t (рис. 2.1). За время τ он продвинется по каналу на расстояние l и займет объем V = F · l = F · w · τ. Масса

поступившего в канал теплоносителя

. Вместе с потоком

теплоносителя переносится энтальпия i = ср·ρt.

 

Количество теплоты,

 

 

 

 

внесенное в канал теплоно-

 

 

 

 

сителем определяется так:

 

 

 

 

W m i V i

 

 

 

 

F w cpt, Дж .

 

 

 

 

Тепловой поток, про-

 

 

 

 

шедший через входное се-

 

 

 

 

чение канала,

 

 

 

 

Q W F w cpt, Вт ,

 

Рис. 2.1. Определение

 

конвективного теплового потока

 

 

 

 

 

а плотность конвективной составляющей теплового потока

 

q Q / F w c

t,[Вт / м2

].

(2.5)

к

p

 

 

 

В результате плотность теплового потока в условиях конвек-

ции с учетом (2.4) и (2.5) будет определяться выражением

 

q w cpt f gradt .

 

(2.6)

61

Закон Ньютона – Рихмана (2.1) позволяет определить только абсолютную величину плотности теплового потока. Однако, поскольку плотность теплового потока является величиной векторной, она может быть определена по компонентам

q

 

w c

t

 

t

;

 

 

f x

 

 

x

x

p

 

 

 

q

 

w

c

t

 

t

;

(2.7)

 

f y

 

y

y

p

 

 

 

qz wzcpt f zt .

Для этого в рассматриваемом объеме теплоносителя необхо-

димо знать поля скоростей

w fw x,y,z, и температур

t ft ( x, y,z, ) (в общем случае нестационарные). Эти функции

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

уравнение переноса энергии в подвижной среде;

уравнение неразрывности;

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

как поле температур и компоненты поля скоростей wx ,wy ,wz

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

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

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

62

2.2. Гипотеза пограничного слоя

Вместе с тем существует приближённый метод решения данной системы уравнений, позволяющий существенно её упростить. Он основывается на использовании теории пограничного слоя, разработанной Л. Прандтлем в 1904 г. В основу теории положена гипотеза о прилипании жидкости к стенке при движении теплоносителя вдоль нее. Вследствие действия сил вязкого трения вблизи поверхности образуется тонкий заторможенный слой жидкости, в пределах которого её скорость изменяется от нуля (на стенке) до скорости невозмущенного потока на удалении от стенки. Этот слой получил название гидродинамического (динамического) по-

граничного слоя. Толщина этого слоя, как правило, возрастает вдоль по потоку.

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

В зависимости от скорости и условий движения теплоносителя режим течения может быть ламинарным или турбулентным. Рассмотрим процесс продольного омывания стенки невозмущенным потоком теплоносителя, движущегося с постоянной скоростью w0 (рис. 2.2). При набегании потока на стенку сначала образуется ламинарный пограничный слой 1, толщиной δл. Если скорость потока достаточно большая, то в дальнейшем образуется турбулентный пограничный слой 2 толщиной δт. Следует отметить, что даже при турбулентном режиме течения вблизи стенки всегда имеется тонкий слой, где течение остается ламинарным 3.

Его называют ламинарным подслоем.

63

Рис. 2.2. Схема возникновения пограничного слоя

Еслитемпературыстенкиижидкостинеодинаковы, то вблизи стенки образуется тепловой пограничный слой, в котором происходит изменение температуры от её значения на стенке до температурывневозмущенномпотоке, гдеонапостоянна. Толщинатеплового слоя зависит как от теплофизических характеристик теплоносителя, так и характера его течения. В общем случае толщины динамического и теплового слоев не совпадают. Это зависит от физических свойств теплоносителей. Соотношение толщин слоев определяется безразмерным числом Прандтля:

Pr f ,

аf

где νf – кинематическая вязкостьтеплоносителя, м2/с; аf – температуропроводность теплоносителя, м2/с. Например, для вязких жидкостей с низкой теплопроводностью Pr 1 и толщина динамического слоя больше теплового. И наоборот, у жидких металлов температуропроводность высокая и Pr 1 . В этом случае тепловой слой превышает толщину динамического. Для газов толщины слоев примерно одинаковы и Pr 1 .

На величину теплового потока между стенкой и теплоносителем основное влияние оказывает тепловое сопротивление погра-

64

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

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

2.3.Дифференциальное уравнение теплоотдачи

впограничном слое

Какужебылоотмечено, вламинарномпограничномслоечастицы жидкости перемещаются параллельно стенке, а перенос теплоты по нормали к стенке в пределах теплового пограничного слоя толщиной δосуществляетсятолькопутёмтеплопроводности(рис. 2.3).

0

Рис. 2.3. Тепловой пограничный слой

65

Тогда, имея в виду, что плотность теплового потока является векторной величиной, уравнение (2.6) принимает вид

q

 

t

.

(2.8)

 

 

 

f n

 

С другой стороны, в соответствии с законом Ньютона – Рихмана плотность теплового потока при теплоотдаче

q tw t f ,

(2.9)

где tw и tf – температуры стенки и потока жидкости за пределами пограничного слоя.

Приравниваявыражения(2.8) и(2.9), получаемуравнениедля расчета коэффициента теплоотдачи

 

f

 

t

.

(2.10)

tw t f

 

n

 

 

 

 

Уравнение (2.10) называют уравнением теплоотдачи в пограничном слое.

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

 

 

f

 

 

t

w

t

f

 

,

 

 

 

 

 

 

 

tw t f

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

f

.

 

 

 

 

(2.11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вместе с тем при кажущейся простоте уравнения (2.10) расчет коэффициента теплоотдачи является задачей достаточно сложной. Длянахождениячисленногозначенияαнеобходимознатьвеличину

66

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

2.4. Система дифференциальных уравнений конвективного теплообмена

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

 

 

(2.12)

q f grad t wi ,

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

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

di cp dt и

i cp d t .

(2.13)

 

t

 

В этом случае поле энтальпий определяется, если известно

поле температур в исследуемой области.

Поля температур

 

находят из решения системы

t(x, y,z, ) и скоростей w(x,y,z, )

дифференциальных уравнений конвективного теплообмена.

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

уравнение переноса энергии в подвижной среде;

уравнение неразрывности;

уравнение движения вязкого теплоносителя.

67

2.4.1. Уравнение переноса энергии

Уравнение переноса энергии представляет собой уравнение теплопроводности, записанное для подвижной среды. Оно описывает температурное поле в движущейся жидкости при известном поле скоростей.

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

Уравнение теплопроводности для твердого тела с учетом

(2.13)

 

t

 

2

t

 

2

t

 

2

t

 

cp

 

 

 

 

 

 

 

 

x

2

y

2

z

2

 

 

 

 

 

 

 

 

запишемчерезкомпонентыплотноститепловогопотокавсоответствии с законом Фурье:

cp

t

 

 

 

t

 

 

 

 

t

 

 

 

t

 

q

qy

 

q

. (2.14)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

z

 

 

 

 

 

 

 

y

 

 

x

 

x

 

y

 

y

 

z

 

z

 

x

 

z

 

Компонентыплотноститепловогопотокавусловияхконвекции:

qx wxcpt f xt ; qy wy cpt f yt ;

qz wz cpt f zt .

Подставляя (2.15) в уравнение (2.14), получим

 

 

 

t

 

 

2

t

 

 

2

t

 

2

t

 

 

 

 

 

cp

f

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

z

2

 

 

 

 

 

x

 

y

 

 

 

 

 

t

 

 

t

 

t

 

 

 

 

 

 

w

 

 

wy

cp wx

 

wy

 

 

wz

 

 

 

cp t

 

x

 

x

y

 

 

 

y

 

 

 

 

z

 

 

 

 

x

 

 

(2.15)

w (2.16)

zz .

68

Последнее слагаемое в правой части согласно закону сохранения массы для несжимаемых жидкостей ( const ) равно нулю

[см. уравнение (2.30)]:

 

 

 

 

 

 

w

wy

 

w

0.

 

 

 

 

 

 

 

 

 

 

div w

 

x

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

z

 

 

 

 

 

 

 

 

 

С учетом (2.17) представим уравнение (2.16) в виде

 

 

 

t

 

 

t

 

t

 

 

t

 

 

 

2

t

 

2

t

 

2

t

 

cp

wx

wy

wz

 

 

f

 

 

 

 

 

.

 

 

 

 

 

2

 

2

 

2

 

 

x

y

 

z

 

 

x

y

z

 

(2.17)

(2.18)

Выражение в скобках в левой части (2.18) представляет собой полную производную от температуры по времени:

dt

 

t

 

t x

 

t y

 

t z

 

t

w

t

w

 

t

w

t

.

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

x

y

z

 

 

y y

 

 

 

 

 

 

x x

 

z z

Здесьслагаемое t характеризуетизменениетемпературы

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

В итоге уравнение переноса энергии (2.18) (уравнение переноса тепловой энергии Фурье – Кирхгофа) можно записать в виде

 

d t

2

Вт

(2.19)

cp

 

f

t ,

м

3 .

d

 

 

 

 

 

С учетом выражения для температуропроводности жидкости af f cp уравнение (2.19) можно записать в другой размерно-

сти:

d t

2

K

(2.20)

 

af

t,

.

d

 

c

 

69