8185
.pdf60
Производная U(x) по направлению x:
2tgγ
U (x) uh0 (d0 2xtgγ)2 .
Находим толщину потери импульса δ** по формуле [34]:
Zδ 2 0,47 xU5dx,
νU6 0
откуда
δ** 2 |
h |
2xtgγ |
6 1 |
|
|
1 |
|
||
0,06 |
0 |
|
|
|
|
|
|
|
. |
|
uh tgγ |
|
h4 |
(h |
2xtgγ)4 |
||||
ν |
|
|
|
|
|
||||
|
|
0 |
|
0 |
0 |
|
|||
Находим модифицированный формпараметр χ: |
|
|
|||||||
|
|
δ** 2 |
|
|
|
|
|
|
|
|
χ |
U . |
|
|
|
|
|
ν
(3.92)
(3.93)
(3.94)
(3.95)
По модифицированному формпараметру χ из таблиц определяем значения f1(χ) и f2(χ) и первоначальный формпараметр λ:
|
f (χ) |
δ* |
; |
|
(3.96) |
|||
|
|
|||||||
|
1 |
|
|
δ** |
|
|
|
|
|
|
|
|
|
|
|
||
f |
|
τ |
ст |
(x)δ** |
|
|||
2(χ) |
|
|
|
, |
(3.97) |
|||
|
|
|
|
|||||
|
|
|
|
μU |
|
где τст (2 λ)μU / δ – касательное напряжение на стенке.
По уравнениям (3.96) и (3.97) находят толщину вытеснения и напряжение
трения на поверхности.
По формпараметру λ и уравнению (3.98) находим толщину пограничного
слоя δ из уравнения |
|
|
|
δ* |
|
3 λ . |
(3.98) |
δ10 120
Расчет пограничного слоя возможен до точки отрыва (λ = –12).
Пограничный слой, в свою очередь, также влияет на потенциальное тече-
ние. Это обратное влияние выражается в том, что образование пограничного слоя сопровождается оттеснением линий массового потока, вследствие чего по-
61
тенциальное течение искажается (изменяется распределение давления) и отрыв наступает раньше. Полный расчет для каждого сечения диффузора приводится в табл. 3.1. Если учесть влияние оттеснения линий тока, то толщина ламинар-
ного подслоя в точке отрыва достигает 5 – 6 мм. Согласно экспериментальным данным – около 4 мм.
Рассчитанные по изложенному методу характеристики пограничного слоя до точки его отрыва получены, кроме того, без учета влияния генерируемой в результате отрыва потока турбулентности. Необходимая для этого информация и формулы для расчета имеются в [38]. В частности, влияние турбулентности на коэффициент трения для пластины в безградиентном потоке воздуха реко-
мендуется учитывать следующей эмпирической зависимостью
cf / cf 0 1 0,25th(0,2Tu), |
(3.99) |
где cf0 – коэффициент трения при степени турбулентности набегающего потока
Tu ≤ 0,1 %. Формула справедлива для Tu до 10 %. Здесь степень турбу-
лентности рассчитана по трем составляющим пульсаций скорости, т. е.
|
u |
2 |
ν |
2 |
2 |
|
|
Tu 100 |
|
|
w |
. |
(3.100) |
||
|
|
U |
|
||||
|
|
|
|
|
|
||
При продольном обтекании цилиндра |
|
|
|
|
|||
cf / cf 0 1 0,22th(0,2Tu). |
(3.101) |
||||||
Аналогичная зависимость для учета влияния степени турбулентности по- |
|||||||
тока на интенсивность теплообмена имеет вид |
|
||||||
St /St0 1 0,41th(0,2Tu). |
(3.102) |
||||||
При продольном обтекании цилиндра |
|
|
|
|
|||
St /St0 1 0,37th(0,2Tu). |
(3.103) |
||||||
Более заметное влияние степени турбулентности набегающего потока на |
|||||||
теплообмен объясняется изменением турбулентного числа Прандтля |
|
||||||
Prт / Prт0 1 0,225th(0,2Tu). |
(3.104) |
Таблица 3.1
Результаты расчета ламинарного пограничного слоя в диффузоре
Скорость в узком сечении (на входе в диффузор) u = 5 м/с; d0 = 0,047 м; γ = 12°; b = 0,12 м
x |
h(x) |
U(x) |
U′(x) |
Z |
δ** |
χ |
λ |
f1(χ) |
f2(χ) |
δ* |
τст |
δ |
|
|
0 |
0,047 |
5 |
-45,2012 |
0 |
0 |
0 |
0 |
2,554054 |
0,234921 |
0 |
0 |
0 |
|
|
0,0024 |
0,048019739 |
4,893820869 |
-43,3018 |
0,00025 |
0,00019 |
-0,010758 |
-0,770 |
2,592143 |
0,221249 |
0,000518 |
0,019912 |
0,001692 |
|
|
0,0048 |
0,049039478 |
4,792057566 |
-41,5197 |
0,00054 |
0,00028 |
-0,022224 |
-1,570 |
2,634273 |
0,206601 |
0,00079 |
0,012396 |
0,002524 |
|
|
0,0072 |
0,050059216 |
4,694440234 |
-39,8453 |
0,00086 |
0,00036 |
-0,034428 |
-2,416 |
2,681779 |
0,190679 |
0,001073 |
0,008581 |
0,003351 |
|
|
0,0096 |
0,051078955 |
4,600720560 |
-38,2703 |
0,00124 |
0,00043 |
-0,047401 |
-3,310 |
2,73544 |
0,173446 |
0,001094 |
0,007805 |
0,00334 |
|
|
0,0120 |
0,052098694 |
4,510669677 |
-36,7898 |
0,00166 |
0,0005 |
-0,061175 |
-4,250 |
2,795916 |
0,154957 |
0,001398 |
0,005578 |
0,004168 |
|
|
0,0144 |
0,053118433 |
4,42407629 |
-35,3879 |
0,00214 |
0,00057 |
-0,075781 |
-5,265 |
2,866194 |
0,134673 |
0,00172 |
0,00404 |
0,005001 |
|
|
0,0168 |
0,054138172 |
4,340745024 |
-34,0674 |
0,00268 |
0,00063 |
-0,091254 |
-6,365 |
2,948648 |
0,112446 |
0,002064 |
0,002891 |
0,005846 |
|
|
0,0192 |
0,055157911 |
4,260494953 |
-32,8193 |
0,00328 |
0,00070 |
-0,107626 |
-7,580 |
3,048027 |
0,087772 |
0,002134 |
0,002257 |
0,005875 |
|
|
0,0216 |
0,056177649 |
4,183158286 |
-31,6387 |
0,00395 |
0,00077 |
-0,124931 |
-8,960 |
3,172656 |
0,059833 |
0,002538 |
0,001346 |
0,006774 |
62 |
|
0,0240 |
0,057197388 |
4,108579200 |
-30,5206 |
0,00469 |
0,00084 |
-0,143205 |
-10,59 |
3,338215 |
0,027332 |
0,002671 |
0,000615 |
0,006878 |
||
|
||||||||||||||
0,0257 |
0,057925773 |
4,056916067 |
-29,7579 |
0,00527 |
0,00089 |
-0,156871 |
-12,00 |
3,500000 |
0 |
0,00315 |
0 |
0,007875 |
|
63
3.8Обобщение опытных данных по сопротивлению и теплообмену
впрофилированных каналах
Наличие в модели участков безотрывного обтекания стенки в подслое и зон рециркуляции, существование которых подтверждено большим количе-
ством экспериментов при изучении течения в профилированных трубах и ка-
налах (см. в частности [15, 22]), позволило A.Л. Ефимову и O.K. Бережной предположить, что зависимости (3.31) и (3.37) могут быть использованы при обобщении опытных данных по теплообмену и сопротивлению профилиро-
ванных труб и каналов. Более того, именно зоны рециркуляции рассматри-
вались в качестве основной причины нарушения аналогии процессов переноса в таких каналах, поскольку многочисленные исследования обтекания уступов, за которыми образуются подобные зоны, показывают, что касательное напряже-
ние на стенке в этих зонах практически равно нулю, в то время как плотность
теплового потока имеет конечное значение.
В результате данные по теплообмену в диффузорно-конфузорных ка-
налах, в трубах с поперечной и винтовой накаткой, со спиральными вставками удалось обобщить зависимостями вида (3.37), в которых Reкp и показатели сте-
пени p при Re были представлены степенными функциями обобщенных гео-
метрических переменных – относительной длины продольного профиля, степе-
ни дросселирования сечения канала, закрутки профиля и др. Среднеквадратич-
ная погрешность при этом составила менее 20 % для Reкp и не более 3,7 % для
показателя степени p при Re [39].
В частности, для диффузорно-конфузорных асимметричных каналов:
Reкр 7667X10,617X20,159X30,475X40,046; |
|
|
|
|
(3.105) |
|
p 0,7024X1 0,039X20,047X30,242X40,008, |
|
|
|
|
(3.106) |
|
|
|
|
|
|
|
|
где X1 = l′ / d′; Х2 = Xd / Xk; Х3 = dmax / dmin; Х4 = L / l0; l |
l2 |
d2 |
|
l2 |
d2 ; |
|
|
d |
|
|
k |
|
|
d 4V / F;V – свободный объем, F – поверхность теплообмена; dmax, dmin – мак-
симальный и минимальный гидравлические диаметры канала; L – общая длина;
64
l0 – длина одного периода диффузорно-конфузорного канала; ld, lk – длины диффузорной и конфузорной частей канала; d = dmax – dmin; Xd = ld′ / dd′,
Xk = lk′ / dk′ – относительные длины диффузорного и конфузорного участков ка-
нала, рассчитываемые аналогично относительной длине периода Х1. Зависи-
мость (3.105) в графической форме представлена на рис. 3.3.
Однако при обобщении данных по сопротивлению в такой же форме по-
грешность оказалась существенно большей. И более приемлемый результат удалось получить при обработке данных на основе представления о течении в профилированных трубах и каналах, как в каналах с дискретно шероховатыми стенками. Для этого за основу была взята формула А.Д. Альтшуля для труб с песочной шероховатостью
68 |
|
0,25 |
|||
λ 0,11 |
|
|
|
|
, |
Re |
|
||||
|
|
d |
|
которая после соответствующей модификации приняла следующий вид:
λ A |
|
n1 |
BX |
n2 |
BX |
n3 |
BX |
n4 |
BX |
n5 |
/ Re CX |
|
0,25 |
68 / Re exp BX1 |
2 |
3 |
4 |
5 |
6 |
, (3.107) |
Рис. 3.3. Зависимость Reкp (X1, X2, X3, X4) для диффузорно-конфузорного канала.
где X1, X2,…, X5 – безразмерные геометрические переменные (симплексы), харак-
теризующие относительную длину периода канала, общее число периодов, сте-
65
пень его сжатия-расширения, кривизну и закрутку профиля и рассчитываемые
по общим для любых типов каналов правилам; Х6 = / (l + d) – относительная
высота профиля канала; l – длина периода и d = 4V / F – гидравлический диа-
метр канала; V, F – объем и площадь поверхности теплообмена периода канала;
А, В, С, п1, п2,..., п5 – эмпирические константы.
Результаты обобщения данных с использованием зависимостей (3.31) и (3.37) показали, что основной причиной изменения показателя степени при Re в
исходных эмпирических формулах вида Nu = A∙Reр и λ = В∙Req (от 0,7 до 0,87
для теплообмена и от –0,32 до 0 для сопротивления), с помощью которых были обработаны опытные данные для каждой исследованной поверхности, является проявление шероховатости, а не изменение показателя в степенных профилях
скоростей и температур турбулентного ядра. |
|
|
|
|
|
|
|
Если воспользоваться зависимостями (3.2) и |
(3.30), то |
при n = 7 и |
|||||
Reкp = 22150 универсальные координаты и y δ |
/ ν, где |
|
|
|
|
– ди- |
|
τ |
ст |
/ ρ |
|||||
0 |
* 0 |
* |
|
|
|
|
намическая скорость, имеют соответственно значения 15 и 50 (рис. 3.4). Таким образом, толщина подслоя в соответствии с предложенной моделью включает не только ламинарный подслой по Л. Прандтлю (y+ < 5), но и большую часть буферного слоя по Т. Карману (y+ ≈ 5…70) [34].
Рис. 3.4. Универсальный логарифмический закон распределения скоростей в гладкой трубе: 1 – действительная кривая перехода от ламинарного течения к турбулентному; 2 – профиль скоростей в ламинарном пограничном слое по решению Г. Блазиуса.
66
4. КОНТРОЛЬНЫЕ ВОПРОСЫ И ЗАДАНИЯ ДЛЯ
САМОСТОЯТЕЛЬНОЙ РАБОТЫ
1. Назовите цель и исходные данные проектного (конструктивного) теп-
лового расчета теплообменного аппарата.
2. Назовите цель и исходные данные поверочного теплового расчета теп-
лообменного аппарата.
3.Для каких схем относительного движения теплоносителей средний температурный напор рассчитывают как среднелогарифмический?
4.В каких случаях необходимо вводить поправку на схему движения теп-
лоносителей?
5. Когда можно использовать среднеарифметический температурный на-
пор?
6. Чем метод среднего температурного напора отличается от метода эф-
фективности?
7.Дайте определение эффективности теплообменника?
8.Что такое число единиц переноса и как его рассчитывают?
9.От чего зависит величина поправки ε t?
10.Запишите формулы для расчета конечных значений температур теп-
лоносителей, используя величину эффективности?
11. Как рассчитать тепловую мощность теплообменника по методу сред-
него температурного напора и по методу эффективности?
12.Чему равен температурный напор в паропреобразователе?
13.Перечислите наиболее распространенные методы интенсификации теплообмена в теплообменных аппаратах.
14.В какой из областей потока относительно стенки (в ламинарном под-
слое, буферной области или турбулентном ядре потока) сосредоточено основ-
ное термическое сопротивление при турбулентном течении в трубе газа и в ка-
ких при течении жидкости?
15. Во сколько примерно раз можно интенсифицировать теплообмен в
67
теплообменнике при вынужденном течении, чтобы рост гидравлического со-
противления не опережал рост интенсивности теплообмена?
16.Чем можно объяснить интенсификацию теплообмена при течении в профилированных трубах и каналах?
17.Какова длина начального участка профилированного канала?
18.Назовите основные признаки модели турбулентного течения с лами-
нарным прерывистым подслоем?
19.Напишите выражение для толщины ламинарного пограничного слоя на пластине.
20.В какой из областей наиболее эффективно применять профилирован-
ные каналы, с точки зрения интенсификации теплообмена? В ламинарной, пе-
реходной, при развитом турбулентном течении?
21.Напишите формулы для расчета локального и среднего теплообмена при ламинарном продольном обтекании гладкой пластины.
22.Напишите формулы для расчета локального и среднего трения при турбулентном продольном обтекании гладкой пластины.
23.Напишите выражение для степенного профиля скоростей.
24.Объясните физический смысл критического числа Рейнольдса для ла-
минарного прерывистого подслоя.
25. Какие из известных вам способов интенсификации теплообмена заре-
комендовали себя как наиболее эффективные и почему?
26. Какие из геометрических размеров поверхностей теплообмена прини-
мают в качестве характерных при расчете конвективного теплообмена в трубах и каналах для ламинарного, переходного и турбулентного режимов, при попе-
речном обтекании пучков труб, при конденсации пара на вертикальной стенке и на пучке горизонтальных труб?
27. Укажите правила выбора определяющих значений скоростей теплоно-
сителей для труб и каналов, пучков труб.
28. Какова особенность выбора определяющего значения температуры теплоносителя при расчете (выборе табличных значений) теплофизических
68
свойств теплоносителей при вынужденном течении в трубах и каналах, попе-
речном обтекании пучков труб, при естественной конвекции, при конденсации пара и кипении жидкостей?
29.Как учитывают влияние начального участка на средний теплообмен в трубах или каналах?
30.С какого ряда пучка труб по потоку теплоносителя, коэффициенты теплоотдачи практически остаются постоянными?
31.При каком значении числа рядов труб средний коэффициент теплоот-
дачи пучка в поперечном потоке перестает изменяться?
32. Какие факторы влияют на интенсивность теплообмена и гидравличе-
ское сопротивление при течении теплоносителей на начальных участках труб и каналов, поперечно обтекаемых пучков труб?
69
СПИСОК ЛИТЕРАТУРЫ
1.Антуфьев, В.М. Эффективность различных форм конвективных поверхностей нагрева. М.-Л.: Энергия, 1966.
2.Бажан П.И., Мунябин К.Л. Обобщенное уравнение для расчета теплоотдачи и трения в каналах // Тр. Третьей Рос. нац. конф. по теплообмену. В 8-ми т. Т. 2. Вынужденная конвекция однофазной жидкости. М.: Издательство МЭИ, 2002. С. 49-51.
3.Бережная O.K., Ефимов А.Л. Расчет развитого турбулентного течения и теплообмена на основе модели прерывистого подслоя // Тр. Третьей Рос. нац. конф. по теплообмену. В 8-ми т. Т. 2. Вынужденная конвекция однофазной жидкости. М.: Издательство МЭИ, 2002. С. 56-57.
4.Барановский Н.В., Коваленко Л.М., Ястребенецкий А.Р. Пластинчатые и спиральные теплообменники. М.: Машиностроение, 1973.
5.Бурков В.В. Алюминиевые теплообменники сельскохозяйственных и транспортных машин. Л.: Машиностроение (Ленингр. отд.), 1985.
6.Величко В.И., Пронин В.А. Интенсификация теплоотдачи и повышение энергетической эффективности конвективных поверхностей теплообмена. М.: Издательство МЭИ, 1999.
7.Величко В. И., Пронин В. А. Расчет теплоотдачи в плоском канале с отрывом и присоединением воздушного потока // Науч. тр. Межвуз. тематический сборник № 54 Интенсификация тепломассообмена в энергетических установках. М.: МЭИ, 1985. С. 84-91.
8.Гортышов Ю.Ф., Олимпиев В.В. Теплообменные аппараты с интенсифицированным теплообменом. Казань.: Изд-во Казанского ГТУ, 1999.
9.Гухман А.А. Интенсификация конвективного теплообмена и проблема сравнительной оценки теплообменных поверхностей // Теплоэнергетика. 1977. № 4. С. 5-8.
10.Дрейцер Г.А., Дзюбенко Б.В., Якименко Р.И. Интенсификация теплообмена и анализ методов сравнения теплогидравлической эффективности