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

8185

.pdf
Скачиваний:
0
Добавлен:
24.11.2023
Размер:
1.43 Mб
Скачать

60

Производная 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.Дрейцер Г.А., Дзюбенко Б.В., Якименко Р.И. Интенсификация теплообмена и анализ методов сравнения теплогидравлической эффективности

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]