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

книги из ГПНТБ / Кутателадзе, С. С. Пристенная турбулентность

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

где

| т ~

.

 

 

 

 

 

 

 

 

Вычисляя

i g

из (3.3.5),

(6.3.2)

и (6 .3.5), для точки от­

рыва

( с / = 0 )

найдем, что

 

 

 

 

 

 

 

 

StK P = ( - f6)ll

V'2

О -

I) -4-

 

 

(8.4.6)

Ограничиваясь

приближением

( 1 — | ) ? А у ~ х | ,

что пра­

вильно только в окрестности стенки

( t i - ^ l - C ô ) ,

но соответ­

ствует оценкам (6 .2 .8), получаем

 

 

 

 

 

 

StK P « - f ( - /бУкр т 1 1

^ -

 

 

(8.4.7)

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

\

= Ц± = StPrRe** - ^ .

 

 

(8.4.S)

Для рассматриваемых условий

( Р Г А ; 1 , критические

парамет­

ры по оценке (6.2.8)) уже при Re** Ä; 10'j

| і % Р 1. Принимая

во внимание это обстоятельство

и решая

систему

уравнений

(8.4.7)

и (8.4.8) при параметрах

(6.2.8) п (6.2.10),

получаем

оценку

 

 

 

 

 

 

 

 

 

 

 

 

 

StK P -

-^ßr-

 

 

 

 

(8.4.9)

Еще раз следует отметить, что эта оценка относится к усло­ виям газового потока ( Р г » 1 ) при совпадающих толщинах теплового и динамического пограничных слоев и критических параметрах (6.2.8).

Ниже приводятся результаты некоторых расчетов по фор­ муле (8 . 4 . 9) .

>,2

 

Re'

103

10"

105

 

 

2 5t,.„

 

 

 

1.0-

 

1,06

0,73

0

 

 

Ы 1,3

0,8: с ° ° °

°

0

 

 

 

 

0,025 0,050 0,075

"Re#* 0,25

Рис. 8.5. Влияние градиента давления па теплообмен по опытам А. И. Леонтьева, А. Н. Обливина, П. Н. Ромаиеико [102].

Как видно, в практически наиболее часто встречающемся диапазоне чи­ сел Рейнольдса теплоотдача слабо зависит от градиента давления (приведенные на рис. 8.5 и 8.6 экс­ периментальные данные подтвер­ ждают этот вывод). Однако в об-

 

 

 

с

Cä'CÄg

OSS*

 

 

 

 

 

 

 

 

e

S

 

 

 

CE

 

 

 

 

 

 

 

 

 

о

/

 

 

 

 

 

в

2

0,6

IIt

 

 

 

в

3

 

 

 

 

в

4

 

 

 

 

ѳ

5

0,4

 

 

 

 

 

 

!

 

 

 

 

 

0,2

0,2

 

0,4

0, ff

 

0,8

 

 

 

Рис.

8.6. Распределение

температур в пограничном слое

 

при различных

значениях формпараметра:

 

/ — / . R e * * ° ' 2 s

= -0.or40;

2 — / R e * * ° . 2

5 = -0,0550;

3 — ; Т ! е * * 0 > 2 ^ = -0,02В2;

4—/Не**0^25

-0,0049; 5

— / Н е * * 0 , 2 5 - 0 . 0 1 2 І .

ласти весьма больших чисел Re зависимость St(/) становится вполне ощутимой.

8.5 Тепловой

пограничный

слой

при обтекании непроницаемой

пластины

металлической

жидкостью

(Рг^;1)

Вследствие того, что число Прандтля намного меньше еди­ ницы, тепловой пограничный слой проникает в область пото­ ка, гидродинамически невозмущенную. Кроме того, в окрест­ ности обтекаемого тела даже вне вязкого подслоя À^>ÀT и особенности распределения турбулентного переноса при до­ статочно малых расстояниях от стенки несущественны. По­ этому для зоны динамического пограничного слоя с большой степенью точности можно написать:

О • . у <8, q ~[%-Ѵ

Cppl-^j

—-,

(8.5.1)

где / определяется по формуле

(3.5.6).

 

Полагая, что турбулентность

внешнего потока

близка к

нулю, имеем

 

 

 

 

6 < y < ô T ) ? = *

ï

-

(8-5.2)

121

Рис. 5.7. Сопоставление расчетов с опытными дан­ ными по теплообмену в жидких металлах [96].

Продольное обтекание одиночной трубы:

/ — турбулентный слоіі

по (7.5.3),

Р г - 0 , 0 1 ;

2 ламинар­

ный слон, Рг=0,01 .

Опытные

точки дл я

Рг=*0,007.

Распределение теплового потока при аппроксимации по граничным условиям не зависит от режима течения и может быть описано формулой (8.4.5). Распределение скорости опре­ деляется формулами третьей главы, поскольку перепады тем­ ператур в потоке жидкого металла обычно невелики из-за большой теплопроводности и физические свойства меняются слабо. Эти уравнения численно решены автором и Е. Д. Фе­ доровичем. Результаты расчетов удовлетворительно описыва­ ются интерполяционной формулой (0,005^Рг<0,05; 10 3 < < Р е < 2 - 1 0 5 )

— = 0,38 ^ а J

(8.S.3)

На рис. 8.7 формула (8.5.3) сопоставляется с опытными данными, а также с расчетом для ламинарного пограничного слоя. Как видно, в жидких металлах турбулентность влияет на теплообмен, хотя и значительно слабее, чем в газах. При течении в трубах влияние турбулентности несколько больше, чем при обтекании пластины.

8.6. Тепловой пограничный слой при больших числах Прандтля

При Рг>-1 тепловой пограничный слой «тонет» в вязком подслое динамического пограничного слоя. Последний вслед­ ствие его малой толщины молено считать плоским как при те­ чении на пластине, так и в круглой трубе.

122

Пусть в области 0 < g < g i (£/' = 0)

 

 

 

^И- =

Рч" .

 

(8.6.1)

Тогда

 

 

 

 

 

 

 

 

 

Ят =

X Рг ßT Ti",

 

(8.6.2)

где Рт = ß РгГ1 -

 

 

 

 

 

 

 

Тепловой поток с точностью до знака

 

 

 

Я =

(Я +

*т)

 

 

(8.6.3)

причем

вследствие

малой

теплоемкости

вязкого подслоя

qœqcT

и, следовательно,

 

 

 

 

 

- = ( 1 + Р т Р г т і " ) а = - -

(8.6.4)

Интегрируя в пределах

 

г) =

0,

•&=0,

г) = т)і, т>=1, имеем

 

 

,

41,

 

 

 

 

2

 

/ т ^ Г

^

Г ' .

 

 

S t T / ^ = Рг С

(8.6.5)

о

Таким образом, независимо от интенсивности затухания тур­

булентности в вязком подслое диффузия

теплоты

и массы

при больших числах Прандтля пропорциональна

квадратно­

му

корню

коэффициента

трения. При Рг -*- ooßT Pr"^> 1,

и решение

интеграла

(8.6.5) не зависит

от значения

гц. При

ß T = c o n s t имеем:

 

 

 

 

 

 

 

 

„. = 3, St = ^ =

ß T ' 3 c ; ' 2 P r - 2

' 3 ,

 

 

(8.6.6)

 

 

 

-

у 6

 

 

 

 

 

 

п. = 4,

St = -j-

ßT '4 c}'2 Pr~3 / 4 .

 

 

 

(8.6.7)

Как

было

показано

в третьей и четвертой

главах,

коэффи­

циент турбулентной вязкости затухает в вязком подслое про­

порционально расстоянию

от стенки

в степени,

не

менее

 

Т

ѵ'

 

 

третьей. При этом ѵг2

и, еслизтртд

= '» т 0

П Р И

" = 3

 

оі lay

ди/ду

 

 

• ^ - « Р ^ / з р г л 4 .

 

 

(8.6.8)

Соответственно турбулентное число Прандтля в вязком под­ слое будет возрастать обратно пропорционально расстоянию от стенки.

123

ю

i l

л

а

в

 

 

 

 

 

 

m г

 

 

 

 

I

 

!

! I

 

 

 

ж ;

 

Е

ё

 

 

 

 

 

 

L U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ME

3-/0

 

 

 

 

 

 

 

/0"

 

 

 

 

'0 е P г

Рыс. Ä.c?. Теплоil массообмсп

при больших

числах Рг:

/ — по 1301): 2 — п о

[272]; 3 — по

[343]; 4 — по [275];

5 — по [313]; б -

экспери­

ментальные данные М. Е. Кишиневского и Т. С. Корниенко.

 

Подставляя

значение

Ат из

(8.6.8)

в

формулу

(8.6.7),

имеем

 

 

 

 

 

 

 

St

V -pr

> — ^ - ^ г т -

 

 

(8.6.9)

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

-—-=Ç>13 —0,115- Значение этого коэффициента, соответствую­ щее вычислениям (см. гл. 4) для закона пристенной турбу­ лентной вязкости (8.6.1), равно 1,29. Совпадение теоретиче­ ских и экспериментальных данных для столь сложной про­ блемы следует считать вполне удовлетворительным.

Связь между числом Стентона и коэффициентом трения при различных числах Прандтля можно представить следую­ щим образом: при С/' = 0_для Рг-*-0, Ргя*1, Р г - ѵ с о , St^= =7^/(с/), St —^с/, St — ]/c/ соответственно.

124

8.7. Вырождение теплового пограничного слоя на адиабатической поверхности

Рассмотрим развитие теплового пограничного слоя, набе­ гающего на "адиабатическую поверхность. Схема такого те­ чения и возникающая деформация профиля температур пока­ заны на рис. 8.9.

В области формирования теплового пограничного слоя ( x O ' i ) имеем условия:

 

А

 

д

Т

9ст.

 

 

 

q = qCr,

 

 

 

у = 0;

 

 

(8.7.1)

у = 8т; q = 0; ~ = 0

 

 

 

В области

адиабатической

 

стенки

(x>xt)

соответствую­

щие условия имеют вид:

 

 

 

 

 

у

= 0;

q~0;

f

- =

0;

]

 

 

 

 

 

 

 

 

(8.7.2)

у = ÔT; q = 0; ~ = 0. |

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

среды

(р,

 

ц.).

Поэтому

 

дальнейшее

развитие теп­

 

лового

пограничного

слоя

 

локализуется

в

границах

 

динамического слоя

(при

 

Л'3>Л'І

о т ~ б ) .

Наиболь­

 

шие тепловые

потоки

и

 

наиболее

высокие

коэф­

 

фициенты

 

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

 

теплопроводности

лока- Рис. 8.9.

Схема тепловой завесы за участ-

лизованы

в

окрестности

ком теплообмена,

вязкого

подслоя, т. е. в

зо­

 

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

дТ

адиабатическую стенку. Вследствие этого область с -щ^ = 0д

125

интенсивно размывается от стенки в глубь потока (расши­ ряется зона с TœTc7). Одновременно вследствие присоеди­ нения массы из внешнего потока к развивающемуся динами­ ческому пограничному слою температура этого слоя посте­ пенно выравнивается так, что при х\-+<х> Г с т - ѵ 7 0 . Меха­ низм выравнивания температуры в окрестности адиабатиче­ ской поверхности можно записать символически следующим образом:

Т-+Т„^Та. (3.7.3)

Толщина потери энергии при постоянных физических свой­ ствах пограничного слоя определяется аналогично толщине потери импульса формулой

где от *— отношение толщины потери энергии к толщине теп­ лового пограничного слоя; | т — отношение расстояния от стен­ ки у к той же толщине бт .

Механизму выравнивания температуры (8.7.3) соответ­

ствует

предельное значение

толщины потерн энергии

 

 

 

1

 

 

бГ ,..„,.„vm

> Г cûdgT.

(8.7.5)

 

 

о

 

Таким

образом, при движении теплового пограничного

слоя

по адиабатической поверхности происходит размывание тол­

щины потери энергии

по толщине динамического

погранично­

го слоя. О

масштабе

такого дефекта можно судить по сле­

дующему примеру при условиях Р г = 1 , ô T = ô , р = 1 :

 

 

a)

f=0,

 

 

(0 =

1",

 

Ъ=1Я,

 

х<х.,

1**

 

П

 

 

 

 

 

 

,

 

'S**

1

от =

„ ^

;і) ( 1

 

 

 

, ;

х

 

< х - >

со, о.

1 + n '

 

 

(1 + п

9

+ 2 л ) '

 

х

 

 

~ >

 

 

 

И | х

 

я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.7.6)

 

 

б ) / = / к р ,

 

С О = ё 1 ; 2 ,

 

 

 

х<хх,

V*

= q ^ V ;

 

 

хг <

х-*

 

со,

7

$*"->-1-

(8.7.7)

В табл. 8.1 приведены

значения ôT * для различных значе­

ний п.

 

 

 

 

 

 

 

 

 

 

 

 

 

Как видно, относительная толщина потери энергии на адиабатической поверхности, расположенной за зоной тепло-

126

10

20

50

40

SO

У, км

Рис. 8.10. Деформация профиля температур в по­ граничном слое на адиабатической стенке за пори­ стым участком [324]:

.ѵ,=52 мм, U = 6,5 м/с, V с т =0 . 61 м/с.

обмена, при х - э - о о становится близкой к единице. На рис. 8.10 приведены измерения деформации профиля температур при обтекании адиабатической пластины набегающим турбулент­ ным тепловым пограничным слоем. Качественно картина впол­ не соответствует условиям (8.7.3). Максимальная деформа­ ция профиля температур в этих опытах соответствовала уве­ личению от в шесть раз при теоретическом пределе, равном, примерно, девяти.

Т а б л и ц а 8.1

f

 

п

 

 

1.0

1 12

0

б т

, X <

хх

0,097

0,076

0,066

0

0

 

 

 

 

 

 

 

0_

,

X - > со

0,875

0,910

0,925

1

~ * *

, X <

xt

0,058

0,041

0,035

0

о т

 

 

 

 

/кр

 

 

 

 

 

 

 

^ * »

X ->

со

0,666

0,666

0,666

0,666

0,г

,

 

 

 

 

127

8.8. Адиабатическая температура стенки при тепловой завесе

Рассмотрим дозвуковое обтекание адиабатической непро­ ницаемой пластины набегающим на нее тепловым погранич­ ным слоем среды с постоянными физическими свойствами. В таком простейшем случае интегральное соотношение энер­ гии (2.8.2) примет вид

 

,

V

 

àT„

О

(8.8.1)

dx

' ГСТ

Т0

dx

и. соответственно,

при

х>х\

 

 

 

 

ат*(Т„—Та)

=

const.

 

(8.8.2)

Константа определяется

по параметрам пограничного

слоя в

сечении Х\:

 

 

 

 

 

 

 

 

Т„

 

° т . 1

 

(8.8.3)

 

^ с т ,

1

 

ô*"

 

 

 

 

 

(Индексом 1 обозначены параметры в

сечении х = х и

т. е. в

начале адиабатического

участка.)

 

 

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

пограничного слоя газа

[102].

 

 

Используя результаты предыдущего параграфа, можем на­

писать, что в рассматриваемых

условиях

Г с т , 1 - Т0

 

 

(S.8.4)

*

1 -і- 2 я

а * * '

где б* и ох, —толщина

динамического

пограничного слоя в се­

чениях Л'>.Гі И Х[.

Для тепловой завесы, создаваемой путем охлаждения (пли

нагрева) участка

от х—0 до х=хх

непроницаемой пластины,

при условии, что

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

динамический пограничный

слой формируется с кромки обогреваемого участка и далее

развивается монотонно на всех участках

пластины, при /г = 1 /7

справедливы соотношения:

 

 

$ =

(*)•••

«.8.5)

128

где Ѳ — отношение температур в (8.8.3). Эту величину обычно называют эффективностью тепловой завесы.

В турбулентном пограничном слое интенсивность теплооб­ мена очень велика, и предельное соотношение (8.8.4) реали­ зуется на сравнительно малом расстоянии от сечения хь По-

т*~т

ст 'о

4

2

0,1

8

6

4

2

 

 

 

 

 

 

 

 

 

10 '

1

2

4

6

8

10°

2

4 £fj

 

 

 

 

 

 

 

 

 

*t

Рис. 8.11.

Э ф ф е к т и в н о с т ь

т е п л о в о і і

з а в е с ы

з а

 

 

у ч а с т к о м т е п л о о б м е н а :

 

 

/ — расчет

по

формуле

(8.8.6);

 

опытные

данные

[301:

2 —1,2< Т "Т° <Г\А;

1 4 < У < 4 0

м/с;

 

 

 

 

Т

 

 

 

 

 

 

 

 

3 4 < - ^ - < 6 , 2 ;

U = 75

м/с.

 

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

простой интерполяционной

формулой

 

Ѳ . ^ ( і +

1 5 , 5 ^ ) - 0 , 8 ,

(8.8.7)

которая при х~>Х\ переходит в (8.8.6). Такого типа интер­ поляция оказывается эффективной и в более сложных тече­ ниях.

На рис. 8.11 показано сопоставление расчета по формуле (8.8.6) с опытными данными.

S.9. Тепловая завеса на неадиабатической поверхности

Если предварительно сформированный тепловой погранич­ ный слой набегает на поверхность с заданным тепловым по­

током так, что

при х<.Х\ не равен дс т

при х>Х\,

то го­

ворят о тепловой завесе на неадиабатической

поверхности.

Вычтем из

уравнения теплопереноса

для

случая

сстф0

это

же уравнение, записанное для адиабатической поверхно­

сти.

Имеем

 

 

 

 

9 З а к а з № 42n

129

Соседние файлы в папке книги из ГПНТБ