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

книги из ГПНТБ / Зысина-Моложен, Л. М. Теплообмен в турбомашинах

.pdf
Скачиваний:
10
Добавлен:
22.10.2023
Размер:
15.66 Mб
Скачать

Из этой формулы видно, что для случая обтекания пластины, когда скорость на внешней границе пограничного слоя равна скорости набегающего потока и, соответственно, а 0 = ат и М 0 = Мо,,

должно иметь место формальное равенство Nu** и Nu$c* при отнесении физических констант к температуре торможения.

В работе [150] решается задача о теплообмене в турбулентном пограничном слое при обтекании поверхности потоком несжимае­ мого газа для Рг = 1. Рассматривается задача о теплообмене и сопротивлении при малых и больших температурных напорах. Решение строится аналогично решению [149] для ламинарного пограничного слоя. Рассматриваются интегральные соотношения импульсов и энергии. Вводятся переменная Дородницына

СО у

^= \ f 0dy==№ dy'

о о

условные толщины динамического и теплового пограничного слоя:

со оо

* 4

t*_

t* Лг, «Г = ] г ( 1 - Я « п

 

(для случая малых температурных напоров гидродинамическая и тепловая задачи разделяются и решаются отдельно) и параметры:

 

 

и д:

x = - ^ - G 1 ( R r ) ;

f t

и ■с, д а .

 

PCp t ои

 

 

Для функции Gi (/?;*) принимается то же значение, что и в ги­ дродинамической задаче по методу Лойцянского [115], только

вместо R в формулу вводится RT, и получается

Gi (r ;*)= 153,2 (r ;*)v*.

Для пластины принимается

 

(IV. 123)

 

X = Хо =

1-

Вводятся

нормированные параметры:

 

 

/т =

Z = JL

(IV. 124)

 

(/т)отр

К-- v !

 

 

Хо

 

ГДб (/т)отр

значение /т в точке отрыва.

 

Эти форм-параметры вводятся в интегральное соотношение

энергии, которое после преобразований

и упрощений

сводится

к виду

 

 

i | = - ^ F T(/T; R;*) +

£ ? «

(IV. 125)

141

где

1-f- т

М /т; ИГ): г — 2/т.

(/т)отр

Далее принимается допущение о наличии аналогии между

нормированными функциями %и /т в ламинарном и турбулентном пограничных слоях и рекомендуется эти функции определять по табл. 4 в зависимости от динамического форм-параметра. Для определения (/т)отр принимается допущение

 

 

 

 

 

(Ы отр . турб

/отр.

турб

 

 

 

 

. .

 

 

 

 

 

 

(/т)отрлам

/отр.

лам

 

 

 

 

 

 

 

 

Отсюда при /отр. турб =

2,

/ отр. л ам

=

 

0,0681

ДЛЯ

( / т ) отр. лам

=

= —0,22

и

( / т)отр. лам = — 0,033

получается

соответственно

(/т)отр . турб = —0,650’

при

Рг =

1

и

(/т) отр. турб =

—0,970 при

Рг = 0,73.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Kux = ~ % - f ^ .

 

 

 

 

 

(IV. 127)

 

 

 

 

 

 

 

Gi(RT)

 

 

 

 

 

 

 

 

При значительных температурных напорах вводятся параметры:

W

 

 

T'w

/ гг **\

f

Я

 

бт

 

 

Gi (r D ;

 

 

 

U

/отр. турб

Т 0 u ( R

);

/т = ^

(/т)оТр. турб

^ о

(IV.128)

 

 

Try

Т га

 

 

 

 

Цш

 

 

g i(r ;*).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£оРс,Я2 T'o G(R");

 

 

 

 

 

 

 

 

 

 

 

 

PoVo^Xo Т°

 

 

 

 

 

 

 

Предполагается,

что влияние теплообмена полностью учиты­

 

 

 

 

 

 

 

вается:

 

1)

введением

в

уравне­

 

 

 

 

 

 

 

ния T J T 0;

2)

изменением

форм-

 

 

 

 

 

 

 

па!раметра /т; 3) выражением вели­

 

 

 

 

 

 

 

чин 6 ,6

, 6, и т. п. в

новой

пере­

 

 

 

 

 

 

 

менной.

На основе

 

этих

допущений

 

 

 

 

 

 

 

принимается,

что

для

пластины

 

 

 

 

 

 

 

(Н =

Н 0 =

1,4 =

const)

/отр не за­

 

 

 

 

 

 

 

висит от

T J T о, ^ вид

нормирован­

 

 

 

 

 

 

 

ных

функций

1(f),

Н (/),

Ж (7)

не

 

 

 

 

 

 

 

зависит от TwlTo и может быть взят

 

 

 

 

 

 

 

таким, как для ламинарного погра­

 

 

 

 

 

 

 

ничного слоя.

Эти

функции

приве­

Рис. 45.

Зависимость

функций

 

дены

на

рис.

45.

Далее

решение

 

I

я, я ; от /

 

 

строится

аналогичным

образом

и

деления локальных

 

 

 

получаются соотношения для опре­

коэффициентов сопротивления итеплоотдачи:

 

 

 

 

 

Cl-

21(f) Г0.

 

 

 

 

 

 

 

 

 

 

 

 

 

G(R**) Тшг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(IV. 129)

 

 

 

 

 

Nux = х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т"

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

< V №

 

 

 

 

 

 

 

 

 

142

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

22 Турбулентное движение вязкой жидкости в гладкой цилиндрической трубе

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

Логарифмический профиль скорости, описываемый уравне­ нием (IV.52), как уже говорилось выше, имеет место во всех уча­ стках трубы, кроме тонкого пристенного слоя, названного лами­ нарным подслоем.

Рассмотрим круглую трубу радиусом а. На оси трубы и —

= Umax,

следовательно, уравнение (IV.52) можно записать в виде

 

= 5,75 lg

+ 5,5.

(IV. 130)

Если произвести почленное вычитание формул

(IV. 130) и

(IV.52),

получим

 

 

 

^ ™ * ^ = 5,751g — ;

(IV. 131)

 

0*

У

 

если вместо и взять среднюю скорость

 

 

a

 

 

 

и ср= i Jи2п (а ~ у) йУ<

 

то

о

 

 

и„ ■-U,ср =

 

 

 

3,75.

(IV. 132)

В отличие от ламинарного движения в трубе, когда UmaJUср — 2, при турбулентном движении это отношение зависит от R и из меняется от 1,3 (при R я» 5- Ю3) до 1,15 (при R *=« 3-108). Объяс­ няется это отмечавшейся выше зависимостью от R степенных профилей скорости при турбулентном движении. :

Формулы для логарифмических профилей скорости содержат динамическую скорость и*, связанную с касательным напряжением трения Тц,. Для того чтобы формулы были определенными, необ­ ходимо найти связь между rw, 0 тт и Ucp. Эта связь определяется формулой сопротивления трубы турбулентному движению жидкости.

Определяемое этим сопротивлением падение давления Ар на участке длиной L трубы диаметром d может быть выражено, формулой

Ар ==%- т ^ ¥ - >

<1УЛЗЗ)

где к — коэффициент сопротивления.

 

143

При равномерном движении жидкости по трубе перепад Др уравновешивается сопротивлением трения xwn dL, т. е.

A n d 2

, ,

£±р

= тшя dL.

Отсюда следует, что

Ар — tw ^ .

Подставляя значение Др в (IV. 133), получим

 

 

 

Ьдо-- Л

8

 

 

 

 

 

 

 

 

откуда

 

 

 

и2

 

 

 

 

 

 

 

 

 

 

 

v] = k

ср

 

 

 

 

 

8

 

или

 

 

 

 

 

 

 

 

 

U c р

2 V2

 

(IV. 134)

 

 

 

 

к

 

 

 

 

 

Найдем связь между к и Rcp = UcPd/v. Для этого из (IV. 130)

с помощью тождественного преобразования получим

 

U

■ U , ср

и,ср

 

 

+ 5,5.

(IV. 135)

 

 

 

 

 

2U,ср

 

Отсюда, используя (IV. 132) и (IV. 134), найдем

 

1

5,75 lg

Ук

р

+ 5,5 — 3,75

 

V I

4V2

Кср

 

 

=

21g(Rcp^ ) - 0 , 8 .

(IV. 136)

Экспериментальные исследования ряда авторов подтверждают эту формулу. Неудобство ее в том, что она дает связь между к и Rcp в неявном виде.

Никурадзе на основе обработки опытных данных получил для

коэффициента сопротивления формулу

 

* А, = .0,0032+ ! g - .

(IV. 137)

Сопоставление формул (IV.136) и (IV.137) с формулой (III.121) для коэффициента сопротивления при ламинарном течении в тру­ бах показывает, что при турбулентном течении коэффициент со­ противления выше.

Логарифмический профиль скорости и логарифмическая фор­ мула сопротивления являются полуэмпирическими соотношениями,

144

справедливыми в той мере, в какой возможно использование чис­ ленных значений полуэмпирических констант турбулентности х = 0,4 и г|0 = 11,5 и в какой справедлива идея о существовании ламинарного подслоя и турбулентного ядра при течении в трубе. Наряду с этими формулами — так называемыми универсальными формулами — до настоящего времени широкое распространение имеют чисто эмпирические соотношения, соответствующие эмпи­ рическим степенным профилям скорости

U

= ( J L \ n

(IV. 138)

Umax

\ О. )

 

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

а. = с R7Pm.

(IV. 139)

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

, _

0,3164

(IV. 140)

Д

_ л ос

соответствует профилю скорости с п = 1/7 и хорошо согласуется с экспериментом до значений Rcp < 2•105.

Используя формулу (IV. 138), можно найти для степенных про­ филей скорости соотношения:

Ucp _

2

(IV. 141)

U m a x

( п + 1 ) ( п + 2 ) ’

 

т

2п

(IV.142)

л + 1’

 

 

5п+1 2 (л—1)

2

 

с = 2n+1r\on+1

[(л + l)(n + 2)]n+1.

(IV. 143)

Кроме степенного профиля скорости типа (IV. 138) часто ис­

пользуется также степенной профиль в координатах

ы/и*; yvjv,

имеющий вид

 

(IV. 144)

V*■= А

Если применить это выражение для границы ламинарного

подслоя, где

. yv* _

 

 

 

можно получить выражение для А в виде

 

 

1

(IV. 145)

А = т]о- \

или т]о = Л"-1 .

10 Л. М. Зысина-Моложен и др.

145

Экспериментально установлено, что в диапазоне 40 <

^ <

< 700 с опытом хорошо согласуется формула

 

 

 

 

i

 

_

 

8.74 ( а ) 1'-.

 

 

(IV. 146)

в диапазоне 80 <

 

<

1100 — формула

 

 

 

 

 

— =-9,60

V

V

.

 

 

(IV. 147)

 

 

V*

\

 

/

 

 

 

Применим формулу (IV. 144) к оси трубы, где у — а,

^

— ^Anaxi

тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

»(=?)*->• ( ^ ) ‘(тйгГ-

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

( ^

) "

+ ,=

л ( " = й )" =

л н ; „ .

 

(IV. 148)

Так как

 

 

 

 

 

 

 

 

 

 

 

 

 

то

 

 

 

 

 

 

/

 

Р^тах ’

 

 

 

 

 

 

 

"+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тш \

 

2

1 Р_п

 

 

 

 

Р^п

 

 

 

 

*\тах>

 

 

 

 

 

 

 

 

 

 

 

 

 

И Л И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_

 

/

 

 

 

 

—2п

1Г2

 

 

Т

О

1

\ " + 1 п

Л + 1

Umax

 

 

 

 

^

\

А )

 

 

^тах

Н

2

 

 

Если обозначить

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

*

 

 

(IV. 149)

 

 

 

 

 

 

 

 

Л «+1

 

 

 

то получим формулу сопротивления

 

 

 

 

 

 

 

 

 

 

 

 

— 2 п

о

 

 

 

 

 

^

 

 

D

л-И Р ^ тах

 

 

(IV.150)

 

 

 

 

Ъ^шах

-- 2— *

 

В частности,

при

/г =

V7

и

Л =

8,7

из (IV. 150)

получится

формула

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

=

0,0225 Rmax4.

 

(IV.151)

 

 

р(7;

 

 

 

 

 

 

 

 

 

146

аналогичная хорошо известной формуле для пластины, соответ­ ствующей закону «одной седьмой».

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

площадь поперечного сечения

(IV.152)

периметр

 

для круглой трубы аг = dl4.

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

работы Ито отношение Хиз (изогнутой трубы)

к Хпр (прямой трубы)

выражается формулой

 

 

 

Яиз

 

2-10,05

(IV. 153)

^пр

 

 

 

 

 

для ламинарного течения

 

 

2,22 -1

 

1—

11,6 0,45

 

R Vajfo

 

 

 

 

при

 

 

 

11,6 <

R (-^-)°'5 <2000,

 

Здесь г0— радиус кривизны оси трубы.

При расчете теплоотдачи в трубах до недавнего времени широко использовалась аналогия Рейнольдса, сущность которой заклю­ чается в предположении равенства коэффициентов турбулентного обмена импульса и тепла (Е = К) в уравнениях (IV.40) и (IV.85). Позже оказалось, что такое равенство имеет место только при Рг = = 1. Аналогия Рейнольдса была впоследствии в соответствии с опытными данными для очень больших и очень малых значений числа Рг уточнена Карманом, Белтером, Мартинелли, и в настоя­

щее время считается, что для трубы,

так же как и для пластины,

T = (E + V ) f - ’

+

и что коэффициенты К и Е связаны друг с другом, но не равны друг другу. Существующие опытные данные по этому вопросу освещены в п. 20.

10*

147

В настоящее время имеются аналитические решения для рас­ чета теплоотдачи в круглой трубе, основанные на той или иной модели турбулентного обмена (краткий обзор приведен в работе [108]), однако каждое из этих решений справедливо в узком диа­ пазоне изменения чисел Рг. В связи с этим для практических рас­ четов правильнее пользоваться частными критериальными урав­ нениями, апробированными экспериментально. Основные из этих уравнений:

1) для газов при 0,5 < Рг < 1,0:

Nu?— 0,022R0,8Pr0'6; NuT= 0,021 R ^Pr0.*5, (IV. 154)

где Nu9— соответствует теплообмену при постоянной плотности теплового потока q на поверхности, NuT— теплообмену при по­

стоянной температуре стенки;

2) для жидких металлов при Рг < 0,1:

Nu9 = 6,3 + 0,003RPr; NuT = 4,8 + 0.003R Рг; (IV. 155) 3) для воды и легких жидкостей при 1 < Рг < 20

Nu9 = Nut=

0,0155R0,83Pr0'5;

(IV. 156)

4) для тяжелых жидкостей и масел при Рг >

20

Nu„= NuT=

0,0118R°'9Pr°'3.

(IV. 157)

Теплоотдача в трубах, так же как и сопротивление, мало за­ висит от формы поперечного сечения, и практически с удовлетво­ рительной точностью формулы (IV. 154)—(IV. 157) могут приме­ няться и для труб некруглого сечения при условии замены диа­ метра трубы в формулах гидравлическим диаметром

dr = 4аг.

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

Эксперименты показывают, что при входе потока в трубу в ней не сразу устанавливается стабильное течение. Вследствие влияния условий входа происходят возмущения потока, которые затухают на некотором участке стабилизации. Для ламинарного течения длина участка стабилизации LBX= 0,06aR, для турбулентного течения LBX 2а (25-^40). Опыты показывают, что интенсивность теплообмена на начальном участке Nu^ выше интенсивности при стабилизированном'течении. На рис. 46 приведены эксперимен­ тальные значения интенсивности теплообмена при различных усло­ виях входа, полученные в работе [183] для течения воздуха. Обработка этих данных применительно к средним по длине числам Нуссельта позволила получить формулу

NtJcp = 1 л__ £_

(IV. 158)

Num '

2а

148

При этом коэффициент с определяется условиями на входе:

Развитый

турбулентный

профиль ск о р о ст и

................................... 1,4

Внезапное с у ж е н и е ..............................................................................

 

6

Колено с изгибом на 9 0 ° .................................................................

7

Колено с

поворотом на

180° ........................................................

6

Рис. 46. Локальная теплоотдача на начальном участке круг­ лой трубы при различных условиях входа

В работе [126 ] в экспериментах с воздухом получено следующее соотношение средних коэффициентов теплоотдачи в изогнутых

(Nucp.из) и прямых (Nucp.np) трубах:

Nucp-из =

] + з

го

(IV. 159)

Nucp. пр

 

 

 

Здесь а — радиус трубы; г0— радиус поворота.

23. Турбулентное течение вдоль шероховатых труб

ипластин

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

случаев шероховатую трубу можно считать аэродинамически глад­ кой и влияние шероховатости не учитывать.

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

Подробное исследование турбулентного течения в шероховатых трубах было проведено Никурадзе. На рис. 47 приведены полу­ ченные в этих опытах значения коэффициентов сопротивления А.

Из рассмотрения этих кривых можно сделать следующие вы­ воды.

149

1. Относительная шероховатость не влияет на значение RH

ина развитие процесса в переходном режиме.

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

шероховатости к/а и не зависеть от R, причем чем больше к1а, тем раньше наступает этот режим, чем меньше к/a, тем до больших значений числа R сохраняется обычная для гладких труб зависи­ мость к = к (R). В автомодельном относительно R режиме ве­ личина к тем больше, чем больше к/a. Физически такой характер зависимости был объяснен из сопоставления величин абсолютной

шероховатости к и толщины ламинарного подслоя 8Л. Очевидно, могут быть рассмотрены три предельных случая:

1) к 8Л— бугорки шероховатости тонут в ламинарном под­ слое, и шероховатая труба практически является гладкой;

2) к > 6Л — бугорки шероховатости выходят за пределы ла­ минарного подслоя, обтекание их становится отрывным, и сопро­ тивление поверхности трубы движению жидкости внутри нее прак­ тически определяется сопротивлением плохообтекаемых тел, не зависящим от числам и пропорциональным скоростному напору; такой режим обтекания называется режимом развитой шерохова­ тости;

3) к «=! 8Л — промежуточный режим, когда наблюдается со­ вместное влияние и шероховатости, и трения.

При построении полуэмпирической теории турбулентного течения в шероховатых трубах вводится допущение

ТШ

(IV. 160)

Ри2к

 

150