Часть 1 Основы теплопередачи
.pdf
тенсивному переносу механической энергии, результатом которого является интенсивное перемешивание жидкости.
Если в потоке имеет место разность температур, то перемешивание приводит к выравниванию температуры за счет переноса теплоты.
Режим движения оказывает существенное влияние на характер течения жидкости, так как он определяет механизм переноса теплоты в направлении нормали к стенке (рис 2.7). Так, при ламинарном режиме перенос теплоты в жидкости осуществляется в основном путем теплопроводности. При турбулентном режиме такой способ имеет место лишь у стенки, в ядре потока тепло переносится за счет перемешивания. В этом случае интенсивность теплоотдачи определяется термическим сопротивлением пограничного слоя,
которое по сравнению с термиче- |
|
|
ским сопротивлением ядра потока |
|
|
является определяющим. |
Рис. 2.7 |
|
2.1.1. Гидродинамический пограничный слой |
|
|
Допустим, что мы имеем продольное обтекание поверхности |
||
твердого тела, причем вдали от стенки скорость частиц жидкости |
и |
|
температура постоянны. Вследствие действия сил внутреннего трения у стенки образуется тонкий слой заторможенной жидкости, в пределах которого скорость частиц у стенки изменяется от W=0 на стенке до скорости в невозмущенном потоке . Как было обнаружено, в этом слое сохраняется ламинарный характер течения. Этот слой был назван гидродинамическим пограничным слоем. В этом слое скорость
частиц жидкости асимптотически приближается к скорости |
, а сама |
толщина гидродинамического пограничного слоя обозначается |
(рис. |
2.7). Другими словами, – это расстояние от стенки, на котором скорость частиц жидкости отличается от скорости частиц в невозмущен-
ном потоке |
на 1 |
|
- 53 - |
Толщина гидродинамического пограничного слоя δ зависит от многих переменных:
( ж ж ж )
2.1.2. Тепловой пограничный слой
Советским теплофизиком Г.Н. Кружилиным было обнаружено существование у стенки теплового пограничного слоя. Тепловой пограничный слой – это слой жидкости у стенки, в пределах которого сосредотачивается основное изменение температуры жидкости. При этом температура жидкости изменяется от значения до значения
в невозмущенном потоке (рис. 2.8).
Толщина теплового пограничного слоя (K) – это расстояние от стенки, на котором температура частиц жидкости в слое отличается от температуры в невозмущенном потоке вдали от стенки на 1 .
Установлено, что толщина теплового пограничного слоя K также зависит от большого числа переменных:
( |
ж |
ж ж |
ж ж |
) |
В общем случае толщины K и |
неодинаковы, однако считает- |
|||
ся, что они одного порядка. Можно показать, что между толщиной теплового пограничного слоя K и коэффициентом теплоотдачи существует определенная связь.
Если учесть, что на стенке существует тонкий слой неподвижной жидкости, то для него справедлив закон Фурье:
Если посчитать, что изменение температуры в этом слое имеет место в направлении одной координаты, а закон изменения температуры в слое посчитать линейным, то получим
Рис. 2.8
- 54 -
С другой стороны, в процессе теплоотдачи
( ж)
Сопоставление этих соотношений позволяет сделать вывод:
( )
Таким образом, коэффициент теплоотдачи является функцией коэффициента теплопроводности жидкости и толщины теплового пограничного слоя K:
( ж ж ж ж ж )
Одной из задач конвективного теплообмена является определение коэффициента теплоотдачи
2.2. Аналитическое описание процесса конвективного теплообмена
2.2.1. Дифференциальное уравнение теплоотдачи
Учитывая, что процесс теплоотдачи определяется не только теплофизическими характеристиками текущей жидкости, но и связан с режимом движения среды у стенки, то этот процесс должен быть описан системой дифференциальных уравнений, полученных на основе известных физических законов сохранения: энергии, количества движения и массы.
Поскольку в процессе теплоотдачи на стенке имеется слой неподвижной жидкости, то для этого слоя может быть применен закон Фурье:
С другой стороны, в процессе теплоотдачи
( с ж)
Таким образом,
( |
|
) |
(2.1) |
|
ж
( )
ж
- 55 -
Сопоставление этих соотношений позволяет сделать вывод о том, что для расчета необходимо знание градиента температуры
в движущейся жидкости у стенки. Последнее может быть получено лишь путем решения дифференциального уравнения, которое описывает температурное поле в движущейся жидкости у стенки и называется дифференциальным уравнением энергии.
2.2.2. Дифференциальное уравнение энергии
Для упрощения вывода предположим, что текущая жидкость
несжимаемая, |
то есть |
Физические параметры |
λ, |
Будем считать, что жидкость однородна и изотропна, а |
|
внутренние источники теплоты отсутствуют. В основу вывода дифференциального уравнения энергии положен закон сохранения энергии, который для данного случая может быть сформулирован следующим образом: количество теплоты dQ (Дж), введенное в элементарный объем dV движущейся жидкости за промежуток времени d извне, идет на изменение энтальпии жидкости, протекающей через этот объем:
Проверим размерность: |
|
|
|
|
|
|
|
||||
|
кг Дж |
м с |
кг Дж К |
м с Дж |
|||||||
|
|
|
|
|
|
|
|
|
|
||
|
м кг с |
м кг К с |
|||||||||
Начало вывода дифференциального уравнения энергии в точности повторяется, как и при выводе дифференциального уравнения теплопроводности.
В потоке движущейся жидкости выделяется элементарный объем, и подсчитывается количество теплоты, которое подводится и отводится от элементарного объема жидкости по всем направлениям
0X, 0Y, 0Z и идет на изменение |
жидкости. В результате избыток |
|||||
теплоты, подсчитанный ранее, запишется |
|
|
||||
( |
|
+ |
|
+ |
|
) |
|
|
|
||||
Сдругой стороны, согласно закону сохранения энергии
-56 -
|
|
|
|
|
( |
|
|
+ |
|
+ |
|
) |
|
|
|
|
|
|
|
|
|||||||
|
Поскольку теплоотдача является процессом при совместном |
||||||||||||
действии конвекции и теплопроводности, то |
|
|
|||||||||||
|
|
|
|
|
|
т п + |
ко |
|
|
||||
|
|
т п |
|
(по закону Фурье), |
|||||||||
|
|
|
|||||||||||
|
|
|
|
ко |
|
|
, |
|
|
||||
где |
- массовый расход жидкости через единицу площади сечения |
||||||||||||
поперечного направлению движения жидкости. Но вместе с массой
жидкости переносится и энтальпия: |
|
|
|
|
||
|
Дж |
|
Дж |
|||
( |
|
|
К |
|
|
) |
кг К |
кг |
|||||
Тогда |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
( |
|
|
|
|
|
|
|
|
|
|
) |
|
|
|
|
|
( |
|
|
|
|
|
+ |
|
|
|
|
|
|
) |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
Аналогично |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
( |
|
|
|
|
+ |
|
|
|
|
|
|
) |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
( |
|
|
|
+ |
|
|
|
|
|
) |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
( |
|
+ |
|
|
|
+ |
|
|
) |
( |
|
|
|
|
|
+ |
|
|
|
|
+ |
|
|
) + |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
+ |
( |
|
|
|
+ |
|
|
|
|
|
|
+ |
|
|
|
|
) |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
Для несжимаемых жидкостей при
( + + )
- 57 -
Тогда
+ |
|
+ |
|
+ |
|
|
|
( |
|
+ |
|
+ |
|
) |
|
|
|
|
|
|
|
Здесь (мс2), тогда правая часть полученного уравнения энер-
гии характеризует изменение температуры в жидкости при переносе теплоты путем теплопроводности.
Производная |
|
характеризует изменение температуры в дан- |
|
ной точке жидкости во времени, тогда сумма слагаемых в левой части
++ представляет собой изменение температуры
при движении жидкости при переходе от точки к точке, то есть является конвективным изменением температуры, поэтому
( ).
Дифференциальное уравнение энергии запишется
(2.2)
где a – коэффициент температуропроводности жидкости. На основании понятия о полной производной:
+ + +
где
Подсчитаем число неизвестных: |
. Число неиз- |
вестных равно пяти, то есть система незамкнута. Для ее решения необходимо добавить к этим двум уравнениям уравнения движения, которые описывают поле скоростей в движущейся жидкости у стенки.
- 58 -
2.2.3. Дифференциальные уравнения движения
Эти уравнения выводятся в курсах гидродинамики. Приведем упрощенный вывод дифференциального уравнения движения для случая одномерного течения несжимаемой вязкой жидкости, в которой скорость частиц жидкости изменяется в направлении одной координаты. Закон изменения скоростей для трехмерного течения уравнений движения будет записан без вывода.
Выделим в потоке жидкости элементарный объем dV (рис. 2.9) и подсчитаем силы, действующие на него. Движение жидкости осуществляется вдоль оси x.
Скорость меняется в направлении оси Y (рис. 2.10).
Вывод дифференциального уравнения движения основан на применении второго закона механики, согласно которому сила равна массе, умноженной на ускорение:
F=mg.
Определим силы, действующие на выделенный объем в потоке жидкости. Разделим их на массовые силы и силы поверхностные. Массовые силы – это центробежные силы, сила тяжести, силы электрической природы. При выводе будем учитывать действие лишь силы тяжести, поверхностные силы – это силы давления и силы трения.
Подсчитаем проекции всех сил на ось движения. Проекция силы тяже-
Рис. 2.9
Рис. 2.10
- 59 -
сти на ось 0X: |
Проверим размерность: |
|||
|
кг |
м |
||
|
|
м |
|
|
|
м |
с |
||
В сечении X действует сила давления, приложенная к грани dydz: |
||||
X →pdydz ;
Для сечения X + dX : |
|
+ → – ( + |
) |
Равнодействующая этих сил равна их алгебраической сумме:
На рис. 2.9 показана сила трения S, которая противодействует движению. В сечении Y действует касательная сила, противодействующая движению:
Y→ -Sdxdz ;
Y+dY→ ( + ) .
Равнодействующая силы трения
d
Но по закону Ньютона
Тогда
( )
|
Проекция равнодействующих всех сил согласно |
второму зако- |
|
ну механики равна массе элемента |
, умноженной |
на ускорение |
|
|
|
|
|
+
- 60 -
где - – полная производная от скорости по времени, учитыва-
ющая инерционные силы. Раскрывая значение полной производной, получим
( |
|
|
+ |
|
|
|
|
|
+ |
|
|
|
|
|
+ |
|
|
|
|
) |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.3) |
|||
|
|
|
|
|
|
+ |
( |
|
|
+ |
|
|
|
|
+ |
|
|
|
|
) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
Без вывода аналогично: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
( |
|
|
+ |
|
|
|
|
|
+ |
|
|
|
|
|
|
+ |
|
|
|
|
|
) |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.4) |
|||
|
|
|
|
+ |
( |
|
|
+ |
|
|
|
+ |
|
|
|
) |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
Аналогичное соотношение для проекции сил на ось 0Z:
( |
|
+ |
|
|
|
+ |
|
|
|
+ |
|
|
|
) |
|
||
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.5) |
|||
|
|
+ |
( |
|
|
+ |
|
|
+ |
|
|
) |
|||||
|
|
|
|
|
|||||||||||||
В векторной форме
+
Появилась шестая неизвестная –p:
( )
Чтобы сделать систему замкнутой, необходимо добавить еще одно уравнение, которое называется дифференциальным уравнением неразрывности.
- 61 -
2.2.4. Дифференциальное уравнение неразрывности
Дифференциальное уравнение неразрывности выводится на основе применения закона сохранения массы. Подсчитывается масса жидкости, которая втекает и вытекает из противоположной грани элементарного объема в направлении трех осей.
Например, d - масса жидкости в кг:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
||
Сделаем проверку размерности для этого выражения: |
||||||||||||||||||||||
|
|
|
|
кг м |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
м с |
кг |
|
|
|
|
|
||||||
|
|
|
|
м |
|
с |
|
|
|
|
|
|||||||||||
Масса, которая вытекает из противоположной грани в направ- |
||||||||||||||||||||||
лении оси 0X: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Избыток массы в направлении всех трех осей, в частности 0X: |
||||||||||||||||||||||
|
|
( |
|
|
|
|
|
|
|
|
|
|
|
|
) |
|
|
|||||
Для получения этого избытка массы функцию |
|
|
разложим в ряд |
|||||||||||||||||||
Тейлора: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
|
|
|
dx+…, |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
( |
|
|
|
+ |
( |
|
|
|
) |
|
|
) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
( |
|
) |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( |
|
) |
. |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
* |
( |
) |
+ |
|
|
|
( |
|
) |
+ |
( |
|
|
) |
+ |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
На основе применения закона сохранения массы этот избыток массы идет на уменьшение плотности жидкости в объеме dV и равен изменению массы жидкости во времени:
…
+ ( ) + ( ) + ( )
- 62 -
