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

книги из ГПНТБ / Алтухов В.А. Основы аэродинамики летательных аппаратов

.pdf
Скачиваний:
38
Добавлен:
30.10.2023
Размер:
16.73 Mб
Скачать

ним охлаждением* Обозначим суммарный поток тепла с .внутрен­ ней стороны через <70борЗаметим, что q06oP может быть как отрицательным, так и положительным.

 

 

 

 

 

 

 

 

Т а б л и ц а 5.3

М а т

е

р

и

а

л

р-Земли,

Р-Солнца

Т = 200—500°К

Т =■ 20U—600°К

 

 

 

 

 

 

Металлы:

 

 

 

 

0,04

-f-0.il

0,10 4-0,49

А1 . . . . ......................

F e ......................................

 

 

 

 

0,06

ч-0,71

0,45

N

i ......................................

 

 

 

 

0,04

4- 0,39

0.4

Сплавы:

 

 

 

 

0,04

ч -0,55

0,53

типа

дюралюминий

легированная

 

сталь

0,12 -+-0,62

0,6

Изоляционные

материалы:

 

0,85

0

стекло .

-

. „ . *

 

плексиглас

 

. . . .

 

0,89

0

Крашеные поверхности:

 

 

 

темные . . . . . . .

0,8 4-0,99

0,97

с в е т л ы е ........................

 

 

0,8 4 - 0,9

0,14 4-0,18

Потоками тепла, возникающими вследствие теплопроводно­ сти вдоль поверхности стенки, будем (пренебрегать, так как о б ы ч ­

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

Результирующий тепловой поток к стенке, идущий на нагре­ вание (охлаждение) стенки, очевидно, будет равен:

Я к онв + *7солн + Я з г м

Ч п у ч £ Я о б о р •

(5.33)

Если представить себе внешнюю обшивку летательного аппарата в виде тонкой стенки толщиной ост и считать, что температура стенки устанавливается мгновенно по всей толщине, то количе­ ство тепла, поглощаемое единицей поверхности обшивки в еди­ ницу времени, будет равно:

Д<? = сстТсЛ т^ ,

(5.34)

at

где t — время;

сст — удельная теплоемкость материала стенки; Тст — удельный вес материала стенки; ост — толщина стенки.

180

Подставив

выражения для

отдельных

тепловых

потоков

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

баланса

(5.33), получим

 

* ( Тг- г„) +

?сРс +

д3Рз ± д0бор - asт* =

(С8Т)С Т L •

(5.35)

 

 

 

 

a t

 

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

Расчет должен проводиться для различных точек поверхно­ сти летательного аппарата. В результате решения уравнения (5.35) получается для определенной точки поверхности зависи­ мость

г 10-

Уравнение теплового баланса (5.35) показывает принципи­ альные пути снижения температуры стенки летательного аппа­

рата.

 

следует

стремиться

Для уменьшения конвективного потока

к наименьшему значению

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

теплоотдачи а. Он

в сильной степени зависит

от плотности воздуха р-

Коэффици­

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

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

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

Лучистый поток тепла от стенки можно существенно увели­ чить за счет увеличения степени черноты стенки.

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

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

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

lit

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

В

этом

случае для установившегося процесса

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

когда

d Т

= 0, уравнение (5.35)

примет более простой вид:

 

at

 

 

 

 

 

а (7\ — Тст) =

ое Тс$ . .

(5.36)

Так как все величины в этом уравнении от времени не зависят, то в каждой точке поверхности тела установится своя равновес­ ная температура Т„ = Тге, удовлетворяющая уравнению (5.36). Эту температуру часто называют р а в н о в е с н о й р а д и а ц и ­ о н н о й т е м п е р а т у р о й . На фиг. 5.16 и 5.17 приведены

Фиг. 5.16

результаты расчетов средней равновесной радиационной темпе­ ратуры для плоских пластин с хордой 1 и 20 м и коэффициен­ том черноты е = 0,8 при турбулентном пограничном слое, а на фиг. 5.18 для передней точки тупоносого тела при ламинарном слое с радиусом закругления 1 м. Из приведенных графиков сле­ дует, что равновесная радиационная тёмпература 7 ге при боль­ ших числах М значительно меньше температуры восстановления.

182

особенно

на

больших высотах

 

 

 

 

(фиг. 5.17).

Влияние длины

тела

Я®

 

 

 

(числа Re) на температуру его

 

 

у

поверхности сравнительно невелико.

т

 

 

§ 7. СВЯЗЬ МЕЖДУ КОЭФФИЦИЕНТАМИ

 

 

 

яде

 

 

 

ТРЕНИЯ

И ТЕПЛООТДАЧИ

 

 

 

 

(ЧИСЛОМ СТЭНТОНА)

 

 

 

 

 

Явления трения и теплопередачи

3X10

 

 

 

тесно .связаны между собой.

Сила

Ш

 

 

 

трения возникает в результате пере-

 

 

 

носа количества движения от одного

 

 

 

 

слоя воздуха к другому. Этот пере-

д

*

в

/г е м

нос происходит или за счет молеку­

 

Фиг.

5.18

лярного движения, или за счет тур­

 

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

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

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

Z

V

т*

2

 

 

 

 

 

1 J

SHt

 

 

Yxi

 

- —1

 

 

 

 

Фиг.

5.19

 

Проведем внутри турбулентного слоя контрольные плоскости

1— 1

и 2—2,

параллельные поверхности

обтекаемой пластины

(фиг.

5.19).

Вследствие турбулентного

течения между этими

плоскостями' через выделенный ими слой воздуха будет проис­ ходить обмен частицами. Обозначим через 0 ,_ 2 весовое количе-

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

183

Если Кд., — Vxl есть разность продольных скоростей на пло­ скостях 2—2 и 1 — 1, то произведение

- ^ ( K v a - K r ,)

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

'(5-37)

6

Вместе с переносом масс жидкости в процессе турбулентного обмена происходит перенос и тепла. Результирующий тепловой поток между плоскостями 2—2 и 1 — 1, отнесенный к единичной площадке, в единицу времени при ср = const запишется в виде

Я\ Яъ ~ ^ 1 - 2 ср (7*2

 

 

 

 

(5.38)

Исключая G,_, из последних двух уравнений, находим

 

 

Ч\ ~

<7г = (Ti - тг) S cp ту1 -----

 

 

 

(5-39)

 

Кд2 — Кд,

 

 

 

 

Обратим, далее,

b ih имаиие на следующее обстоятельство.

Вну­

три турбулентного

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

у самой

стенки

сущест­

вует ламинарный подслой толщиной h

(фиг. 5.20).

Расположим

 

тсштрольную 1ПЛоскость

1 — 1

 

на высоте

у => /г,

а

кон­

 

трольную плоскость

2—2 на

 

внешней границе пограничного

 

слоя.

 

 

 

 

 

 

На внешней границе погра­

 

ничного слоя при у '■>

 

 

Ъ = ,*я = 0 и Яг = Яй— О-

 

Для

контрольной

плоскости

 

1— >1 :

=

Vxl =

Vh.

 

 

Для рассматриваемого случая расположения контрольных плос­

костей уравнение (5.39)

перепишется в виде

 

4h

= 4Cp g

Tb ~ T h

(5.40) .

v%- v h

 

184

При линейном законе изменение скорости и температуры по ламинарному подслою получим согласно закону трения Ньютона

Ч

= *сТ=

Ц-^-'

(5.41)

и по закону Фурье

 

А

 

 

 

 

7Л=

<7сТ^ ^

^ .

(5.42)

С другой стороны, для теплорого потока по закону теплоотдачи Ньютона можно написать

?лт=^ст = а (7 'с 1 -Тъ).

(6.43)

Решая совместно' систему четырех уравнений (5.40)

(5.43) и

исключая из нее величины qh>Th. и А, получим

 

(5.44)

V,

V,

1 4- - 1 А ( Р г - 1)

 

К /

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

Для определения отношения К воспользуемся приведен­ v *

ными в предыдущей главе соотношениями

КаВ

Re,

VJi

Re* 160;

8 = 0,037 J L l .

• l/R e

С помощью этих соотношений без особых затруднений най­ дем

К

2,135 Re70,1

(5.45)

 

Уь

Полученное соотношение позволяет выражение (5.44) привести к виду

= ________ Тст Ср g

(5.46)

* * .V i[l+2,135R e;W (Pr — 01

1 8 5

Переходя здесь к безразмерному критерию Стэнтона

g?cp Vs

и учитывая, что

 

Чт

=

cf*

 

 

получим

?V*

=

2

 

 

 

 

 

 

 

s,- =

v

i

-

<5-47>

где

_

 

 

 

 

 

5 = 1 + 2,135 ReJ0'1(Рг - 1).

(5.48)

Величина 5 носит название фактора обобщенной аналогии Рейнольдса.

Уравнение (5.47) и определяет искомую зависимость между приведенным коэффициентом теплоотдачи — числом Стэнтона и коэффициентом трения.

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

обобщенной аналогии Рейнольдса 5 получается несколько отлич­ ным от соотношения (5.48). Однако и формула (5.48) дает доста­ точно хорошие для практики результаты.

Если перейти в формулах (5.47) и (5.48) к средним значе­ ниям числа Стэнтона и коэффициента трения Cj и принять чис­ ло Re*, входящее в фактор аналогии Рейнольдса, постоянным

по всей

длине пластины

и разным

среднему

значению

п

Re

Vco ftp»

,

где

.

пластины,

то получим.

Recp =

= -------—

о — длина

 

•2

р.ет2

 

 

 

 

 

следующее выражение для критерия Стэнтона:

 

 

 

St =

JC________ 1

 

(5.49)

 

 

 

2

1 +

2,2Re-°-1(Pr — 1)

 

Для

ламинарного

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

формула

(5.47) для

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

п

 

ЗГ= РгТ .

(5.50)

Это же значение S дает неплохие результаты и для ламинар,- ного слоя. Учитывая вышесказанное, в дальнейшем во всех слу­ чаях будем принимать как наиболее простое выражение для фактора аналогии Рейнольдса по формуле (5.50).

(8 6

Формула (5.47) для числа Стэнтона оказывается применима

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

ифактор аналогии Рейнольдса определять по некоторой средней температуре пограничного слоя.

Коэффициент теплоотдачи в этом случае может быть опреде­ лен по формуле

 

 

a = St*p*Cp* Кос.

(5.51)

Значения

р*, ср* и St*

определяются по некоторой средней ‘по

толщине

температуре

слоя. Эта средняя температура

носит

название

о п р е д е л я ю щ е й т е м п е р а т у р ы .

 

§ 8. ОПРЕДЕЛЯЮЩАЯ ЭНТАЛЬПИЯ. ОПРЕДЕЛЯЮЩАЯ ТЕМПЕРАТУРА

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

В случае наличия теплопередачи меняются температуры по пограничному слою и, следовательно, изменяется и средняя — определяющая температура.

■При температурах < 700°К, когда удельные теплоемкости воз­ духа можно считать постоянными, значение определяющей тем­ пературы дается следующей формулой:

Т * = у (Г5 +

Тст) + 0,22 (Tr - Tt).

(5.52)

Эта формула получена на

основании точных решений

задачи

о теплопередаче для ламинарного пограничного слоя. Но, как приближенную, ее можно применять и для турбулентного слоя.

Структура формулы (5.52) достаточно ясна. Первый член равеи среднеарифметическому значению температур на стенке и на внешней границе пограничного слоя. Второй член учитывает в среднем по толщине пограничного слоя повышение темпера­ туры вследствие кинетического нагрева. •

Учитывая формулу (5.3) и формулу для температуры тормо­

жения

(1.32),

определяющую температуру при ср = const

мож­

но выразить через число М на внешней

границе пограничного'

слоя

-

 

 

 

 

 

j*

1_

 

 

+ 0.22

г Мг\

(5.53>

 

2 1 +

( 1 +

^ 52 ) ^ci

 

 

 

 

где —, а г — коэффициент восстановления температуры. Го

1 S T

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

t* = — (г8 + г’ст) + 0,22 (ir it).

(5.54)

При отсутствии диссоциации энталыпия

i* i(T*)

с*„, Т*,

гд? с*т — среднее'значение удельной теплоемкости

ср в интер­

вале температур от 0 до Т* °К.

 

Для диссоциированного воздуха г* =

i(p, Т*). Разность

гг — h = г * ~

>

(5-55)

где г* — коэффициент восстановления температуры находится

по формуле

 

,- * = т |/р + .

(5.56)

Величина Рг* находится по определяющей температуре Т*.

При определении коэффициента восстановления V*

следует

принимать п = 2 для ламинарного слоя и п —- 3 для турбулент­ ного слоя.

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

(5.5/)

и учесть, что энтальпия торможения

'о = *6 +

2g

то с помощью (5.55), (5.56) и (5.57) из уравнения (5.54) можно получить следующее выражение:

 

 

г *

 

W

+ 0,22 г*

(5.58)

 

 

h

 

 

 

 

 

 

 

В

случае

теплоизолированной

стенки, когда

/•ст == ir и

Тсх=

Тг ,

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

 

 

 

 

;*

 

v 2

(5.59)

 

 

 

- ^ = 1

+ 0,36г* Д 5- ,

 

 

 

h

 

g'h

 

а при

=

+4

== const и Гст — Тг

из формулы

(5.53) будем

иметь

 

 

 

 

 

 

 

 

 

— =

1 + 0,144 гМ г2.

(5.60)

При температуре стенки, равной температуре торможения, ■формула (5.53) переходит в формулу (4.57), полученную нами для адиабатического пограничного слоя в § 8 гл. IV.

188

§ в. ТРЕНИЕ ИТЕПЛООТДАЧА НА ПЛАСТИНЕ ПРИ НАЛИЧИИ КИНЕТИЧЕСКОГО НАГРЕВА

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

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

Действительно, подставив в формулу (4.44) значение коэф­ фициента В и показателя степени а для ламинарного и турбу­

лентного слоев, будем иметь .

 

 

= 0,054 Re»-*

(у-)" '.

(5.61)

С/ л

Р \ Ps /

 

где число Re определяется по параметрам на внешней границе пограничного слоя.

' При степенной зависимости вязкости от температуры с пока­ зателем степени п = 0,76 получим

 

-j j- =

0 ,0 5 4 R e°’3 ^ | r j

.

 

(5 .6 2 )

При малых

числах

М отношение

Т~

а

с

—-- = 1,0,

> 1,0

 

 

 

Т*

 

cfa

(см. гл. IV, §

8). С ростом чисел М величина отношения ~Т-

уменьшается и

приближается к единице.

 

 

cf л

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

На фиг. 5.21 показано влияние точки перехода на коэффици­

ент трения -пластины иля чисел Re=>107 и R e =

106 в зависимо­

сти от числа М.

Как и следовало ожидать, влияние положения

точки перехода

на коэффициент трения резко

снижается при-

больших числах М.

Обратимся теперь к теплопередаче при больших числах М. Согласно методу определяющей температуры, критерий Стэнтона

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

(5.47) и

(5.50) с коэффициентом трения зависимостью

 

q

2

 

St* = - j -

Рг*~т .

(5.63)

189

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