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

Методическое пособие 769

.pdf
Скачиваний:
23
Добавлен:
30.04.2022
Размер:
8.47 Mб
Скачать

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

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

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

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

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

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

wx wx x

wx

1

 

rwy 0;

 

 

 

 

(7.21)

 

 

 

 

 

x

r r

 

 

 

 

 

 

 

 

 

 

 

wy

 

 

1 dp

 

 

 

wy

 

(7

wy

 

 

 

 

 

 

 

 

 

 

 

r

 

 

;

r

dx

 

 

 

 

r

 

 

 

 

 

r r

 

 

 

w

wz w

 

wz a

 

r

wz

(7.23)

y

 

x

x

r

 

 

 

r

 

 

 

 

r r

 

 

при соответствующих граничных условиях:

x 0,

0 r r0,

wx wx0,

wy wy0,

wy wy0;

(7.24)

80

x 0,

r 0,

wx 0, w

y

0,

wz 0;

 

(7.25)

 

 

r

 

r

 

 

 

 

 

 

 

 

x 0, r r0, wx

wy 0, wz 0 t tc

или qc wz

r r r0 .

(7.26)

Указанная система уравнений справедлива при следующих допущениях: - жидкость является несжимаемой средой; -теплофизические свойства жидкости постоянны;

-рассматривается стационарный процесс теплоотдачи;

-не учитываются сила давления, массовые сил и вязкая диссипация;

-длина трубы существенно превышает ее диаметр;

-осевой теплопроводностью потока пренебрегают.

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

а б Рис. 7.6. Изменение температуры жидкости по длине трубы при граничных условиях

первого (а) и второго рода (б) [10]: а. 1 – rr0 0,8; 2 – rr0 0;

б: 1 – rr0 1; 2 – rr0 0,9; 3 – rr0 0

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

81

где ки;

t

 

 

tж2 tж1

 

tвх

tвых

tб

tм ,

(7.27)

 

 

 

 

лог

 

ln

tc tж1

 

ln

tвх

 

ln

tб

 

 

 

 

tc tж2

tвых

tм

 

 

 

 

 

 

 

 

 

tж1,tж2 – температуры жидкости на входе и выходе; tc – температура стен-tвх, tвых – разности температур жидкости на входе и выходе; tб, tм

большая и меньшая разности температур соответственно.

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

быть определена исходя из среднеарифметических значений tа t1 t2 2.

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

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

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

Рис. 7.7. Теплопередача через стенку трубы

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

82

противотоком. Если течение осуществляется в перекрестных направлениях, то схема течения носит название перекрестного тока (рис. 7.8).

 

 

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

б

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г

д

Рис. 7.8. Схемы течения теплоносителей: а – прямоток; б – противоток; в – перекрестный ток; г – смешанная схема; д – многократный перекрестный ток

С точки зрения энергетической эффективности, наиболее предпочтительной является схема течения с противотоком теплоносителю. Однако в условиях большого перепада температур, прямоток является наиболее предпочтительным с точки зрения сохранения поверхности теплообмена [7].

Из уравнения теплового баланса можно определить водяные эквиваленты:

G1cp1 t1' t1" G2cp2 t2" t2'

W1 t1 W2 t2,

(7.28)

где W1 G1cp1,W2 G2cp2,

– водяные

эквиваленты

теплоносителей;

t1 t1' t1"; t2 t2" t1' – среднемассовые температуры теплоносителя в теплооб-

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

На рис. 7.9 представлены изменение средней температуры теплоносителя при различных схемах течения.

Например, при W2 W1 при обеих схемах для теплоносителя с индексом 1 условия теплообмена близки к граничным условиям первого рода. При W1 W2 для противотока строго выполняется условие постоянства теплового потока (граничные условия второго рода).

83

а

б

Рис. 7.9. Изменение средней температуры теплоносителей при прямотоке (а) и противотоке (б)

Уравнение (7.28) может быть использовано в равной степени как для прямоточной, так и противоточной схемы течения. Оно справедливо для условий постоянства массового расхода теплоносителей и коэффициентов теплообмена вдоль поверхностей. Для перекрестной схемы тока выражение для разности температур может быть определено как

t tлог t ,

(7.29)

где t – поправка, определяемая из графика (рис. 7.10), а ее составляющие определяются согласно зависимостям

P t2" t2'

t1' t2' t2

tmax ; R t1' t1"

t2" t2' t1

t2 .

(7.30)

84

Рис. 7.10. Н омограмма для определения поправки t

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

7.3. Рас етные зависимости для определения коэффициент ов теплоотдачи при течении теплоносителя внутри труб ы

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

в работах [7, 8, 10, 11], так

и вновь полученные для

различного класса материа-

лов (пористые, с тчатые и

др.) [4, 5]. Теплофизика

е стоит на месте, поэтому

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

течении в микро- и наноканалах [1 , 13].

 

Теплоотдача при ламинарном режиме течения. Ламинарный

режим

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

85

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

Для ламинарного движения потока теплоносителя Re 2300 в условиях

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

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

При вязкостно-гравитационном режиме на вынужденный конвектив-

ный теплообмен оказывает влияние свободная конвекция, а силы вязкости и подъемные силы при этом соизмеримы.

Критерием того или иного режима является число Рэлея. При Ra 8 105 наблюдается вязкостный режим течения, а при Ra 8 105 – вязкостногравитационный.

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

При стабилизированном ламинарном течении в круглой трубе безразмер-

ный коэффициент теплоотдачи принимает следующие значения [10]:

 

Nu 3,66

при Tc const;

(7.31)

Nu 4,36

при qc const.

 

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

Для начального термического участка при qc const для определения ло-

кального безразмерного коэффициента теплоотдачи возможно использовать следующую зависимость [7]:

Nu

x

0,33Re0,5

Pr0,43

Pr

Pr

0,25 x d 0,1.

(7.33)

 

x

ж

ж

с

 

 

В качестве определяющей температуры теплоносителя выбрана средняя температура, а для стенки – локальная температура.

Осредненный коэффициент теплоотдачи определяется как

 

1,4 x l .

(7.34)

86

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

 

 

d 0,33

 

 

 

0,14

 

 

Nu 1,55

Pe

 

 

 

c

 

l ,

(7.35)

 

 

 

 

l

 

ж

 

 

где

 

1 l 1,7

 

 

1 l

– поправка на гидродиамический началь-

l 0,1

 

 

 

 

 

1 2,5

 

 

 

 

 

 

 

 

 

 

Re d

 

 

Re d

 

ный участок. Величину поправки рассчитывают при

1

 

l

0,1,

в остальных

 

 

 

 

Re d

1

 

 

l

 

 

случаях принимают

l 1. Формула (7.35) справедлива

 

при

 

 

0,01

и

 

 

 

 

 

 

 

 

 

 

Pe d

ж

0,07 с ж 1500.

Если расчет ведут для газов, то величину поправки c

принимают равной 1.

При вязкостно-гравитационном режиме течения средний безразмерный коэффициент теплоотдачи может быть определен как

Nu 0,15Re0,15

Pr0,33

Ra0,1 Pr

Pr

0,25

,

(7.36)

ж

ж

ж

c

l

 

 

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

Таблица 7.6 Значение поправки на гидродинамический начальный участок при

вязкостно-гравитационном режиме течения

l d

1

 

2

5

10

30 15

40 20

l

1,90

1,70

1,441,28

1,18

1,131,05

1,02

1,0

Теплоотдача при турбулентном режиме течения в трубах. Наиболее удобное практически и физически ясное уравнение для безразмерного коэффициента теплоотдачи было получено в 1970 г. Б.С. Петуховым:

 

 

RePr

 

 

 

 

 

Nu

 

 

8

 

 

,

(7.37)

 

 

 

 

 

1 900

Re 12,7

Pr2 3 1

 

8

 

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

1,82lgRe 1,64 2.

(7.38)

87

Диапазон применимости формулы 5 103 Re 5 106; 0,7 Pr 200. Частным случаем данной формулы является

 

Nu 0,0119Re0,9 Pr1 3

 

при Re 2 104;

Pr 20.

 

 

 

 

Для низких чисел Рейнольдса Re 2300 5000 может быть использо-

вана зависимость

 

 

 

 

 

 

 

 

Re 1000 Pr

 

 

Nu

 

 

8

.

(7.40)

 

 

 

 

 

 

 

 

 

 

1 12,7

8 Pr2 3 1

 

При этом наряду с формулой (7.40) возможно использование следующих

зависимостей, которые показывают хорошую сходимость:

 

 

Nu 0,023Re0,8

при Pr 1,

(7.41)

 

Nu 0,023Re0,8 Pr0,4

при Pr 0,7 7,

(7.42)

 

Nu 0,0118Re0,9

Pr0,33при Pr 20.

(7.43)

Теплоотдача при течении теплоносителя в шероховатых трубах. Теп-

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

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

Зависимость гидравлического сопротивления от режима течения, числа Рейнольдса для различной относительной высоты шероховатости представлена на рис. 7.11 [10].

На практике применяют следующие значения kэ , мкм: для стеклянных и

цельнотянутых труб из цветных металлов – 1,5 ;для бесшовных труб из углеродистой стали – 14; для сварных труб из углеродистой стали – 50.

88

Рис. 7.11 . Опытны данные коэффициента гидравлического сопротивления в зависимости от режима течения и шероховатости [10]:

 

1 – ламинарное стабил

изированное течение 64 Re;

2 – формула Блазиуса 0,316Re 0,25; 3 – форм ула Прандтля 2lgRe ,8 2 ;

kэ – эквивалентная шероховатость; ks фактическая шероховатость

Влияние

переменности те плофизических свойств на коэффициент

теплоотдачи. В

ольшинстве практических теплофизических задач принимает-

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

ее температуры от 20 до 0°С, вязкость сн

зится почти втрое, что повлечет за

собой пропорциональное увеличение числа Рейноль дса в 3 раза и коэ

фициен-

та теплоотдачи примерно в 2,3 раза. Если

качестве теплоносителя

спользо-

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

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

89