Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Прикладная гидрогазодинамика..pdf
Скачиваний:
90
Добавлен:
15.11.2022
Размер:
19.64 Mб
Скачать

Глава 17 ТУРБУЛЕНТНЫЕ СТРУИ

В природе и технике широко распространены течения жидко­ стей с тангенциальным разрывом скорости. Течения жидкостей оди­ накового агрегатного состояния по обе стороны поверхности тан­ генциального разрыва называются струйными. Свободные, струи не ограничены стенками. Характерным примером свободной турбу­ лентной струи является поток газов, истекающих из сопла реактив­ ного двигателя в атмосферу. На рис. 17.1 схематично показаны ус­ ловные границы такой струи и «мгновенная фотография» ее турбу­ лентной структуры (под осью) с некоторыми ее характеристиками.

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

метры торможения Т0*ф Тн*, ро*фрп*, Qo*¥:Qa*- При этом

важно

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

разрыва

Ро=Рш в дозвуковых и расчетных сверхзвуковых

струях,

которые

являются изобарными dp/dx=dpfdy=0.

 

 

Поверхность

тангенциального

разрыва скорости при Re =

= е°^ц°—

0 >

Ю3 неустойчива

и генерирует

турбулентность

W)

 

 

 

 

 

(см. п. 6.2), которая называется свободной, так как ее возникнове­

ние не связано с присутствием стенок. Интенсивность

свободной

турбулентности &= и'/ит (ит— скорость на оси струи)

достигает

20% и существенно превосходит естественную степень

турбулент­

ности в трубах (2 ... 5%) и в развитом пристеночном турбулентном пограничном слое (8... 12%). Струйная турбулентность определя­ ет смешение сред, структуру струй и их название — свободные тур­ булентные струи.

К л а с с и ф и к а ц и я т у р б у л е н т н ы х струй*:

I.По числу М0 = Ыо/йо истечения из сопла:

1)сверхзвуковые нерасчетные (см. п. 13.4) и изобарные расчет­

ные;

*При Re<103 могут существовать ламинарные струи, редко встречающиеся

втехнике и здесь не рассматриваемые.

2)изобарные дозвуковые.

II. По относительному движению активной струи и окружающей

среды:

1) спутные — струи и окружающая среда движутся водном

направлении с разными скоростями и0 и ив (см. рис. 17.1); 2) затопленные — струи вытекают в неподвижную окружающую

среду «н=0; 3) встречные — струи и окружающая среда движутся в противо­

положных направлениях; 4) струи в сносящем потоке — направления потоков пересека­

ются;

5)закрученные струи. III. По форме сопла:

1)плоско-параллельные;

2)осесимметричные и другие.

IV. По соотношению и состоянию параметров газа:

1)изотермические однородные T0 = TBi QO = Q H> Со= сн;

2)подогретые и охлажденные Т0^ Т Я\ QO^Q H*,

3)содержащие примесь

4)двухфазные — содержащие примеси в виде мельчайших

взвешенных жидких или твердых частиц; 5) струйные течения, сопровождающиеся химическими и дру­

гими реакциями и т. д.

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

Теория струй является частью более общей теории погранично­ го слоя. Однако практическая значимость струйных течений и мно­ гочисленные исследования, обобщенные в капитальных моногра­ фиях Г. Н. Абрамовича [2], А. С. Гиневского [10] и др. * выделили теорию турбулентных струй в самостоятельный раздел гидрогазо­ динамики.

З а д а ч а т е о р и и с т р у й . Д а н о : 1) ширина 260 плоского сопла или радиус Ro осесимметричного; 2) равномерные поля пара­ метров потоков на срезе сопла е0, Ро, Т0, Q0, и0 и в окружающей сре­ де ен, Рн= Ро, Т1ЪQ„, с„, ин (см. рис. 17.1).

О п р е д е л и т ь параметры струи, т. е. поля скоростей, темпера­ тур, плотностей и концентраций: и= и(х, у); Т=Т(х, у); Q= Q(X, у);

с=с(х, у).

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

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

При анализе турбулентных струй используются:

а) полуэмпирическая теория турбулентности (см. пп. 6.1 6.5); 6) основные уравнения гидрогазодинамики; в) экспериментальные данные.

17.1. ОБЩИЕ СВОЙСТВА И СТРУКТУРА ТУРБУЛЕНТНЫХ СТРУЙ

Рассмотрим схему осесимметричной неизотермической струи в спутном потоке (см. рис. 17.1). Такая струя характеризуется степе-

нью подогрева 0= 7'о/Тн, отношением плотностей n = Q B/Q0,

спутно-

стью m = u ju 0 и разностью концентраций избыточного

элемента

Со—£н»

 

Как уже упоминалось, на границе струи у кромок сопла имеет

место тангенциальный разрыв скорости и0—ив, dufdy=oo

и возни*

* Абрамович Г. Н., Крашенинников С. Ю., Секундов А. Н. Турбулентные те*

чения при воздействии объемных сил и неавтомодельности. М., «Машинострое­ ние», 1975, с. 96.

кает специфическая, свободная турбулентность *. Поперечный пульсации переносят моли активной струи за пределы цилиндри­ ческой поверхности RQ. Эти моли сообщают часть своего количест­ ва движения внешней среде, вовлекая ее в движение вдоль оси х и одновременно сами подтормаживаются. Вместо жидкости, выне­ сенной пульсациями из активной струи, поперечные пульсации вно­ сят в активную струю через поверхность, радиус которой R0y моли жидкости из окружающей среды, которые подтормаживают актив­ ный поток и ускоряются сами. Механизм обмена теплом и избы­ точным компонентом аналогичен. Подторможенные частицы ак­ тивного потока, перемешанные с вовлеченными в движение части­ цами окружающей среды, образуют турбулентный пограничный слой струи (слой смешения), толщина которого у кромки сопла равна нулю и возрастает в направлении течения. Турбулентный об­ мен ликвидирует тангенциальный разрыв, так что в любом сечении

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

распределе­

ние осредненных скорости, температуры, плотности

и концентра­

ции.

 

Рассмотрим сначала поля скорости в спутной струе, а затем

поля температуры и концентрации.

О—Л совпа­

В н е ш н я я г р а н и ц а п о г р а н и ч н о г о с л о я

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

а

#„= 0,01 {ит — ин) при y = Rrр,

(17. 1)

где ит— скорость

на оси струи в данном сечении х.

на внеш­

П о п е р е ч н а я

с о с т а в л я ю щ а я с к о р о с т и v

ней границе направлена к оси струи

(линия тока v—и на рис. 17.1),

так как в струю вовлекается масса

окружающей жидкости. Это

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

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

П е р е х о д н ы й у ч а с т о к с т р у и имеет относительно не­ большую протяженность хп/*н~ 1,2 1,5. В области, примыкающей

коси струи переходного участка, происходит интенсивная турбули-

*В реальных течениях поля скоростей неравномерны из-за присутствия на стенках сопла пограничного слоя. Струйная турбулентность возникает из-за на­ личия конечного, но большого градиента скорости dujdy, а расчеты необычайно

усложняются.

зация и перемешивание жидкости, только что вышедшей из ядра постоянной скорости, т. е. степень турбулентности еи^и'/ит увели­ чивается (см. рис. 17.1). Это приводит к уменьшению осевой ско­ рости ит и деформации полей скорости (а также полей температу­ ры и концентрации) в поперечных сечениях переходного участка. Деформация полей заканчивается к сечению хи, где степень турбу­ лентности на оси струи достигает максимума.

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

Задача 17.1. Нарисуйте схему совмещенных полей скорости u= f(y) для различных сечений основного участка струи (см. рис. 17.1).

Поля абсолютной скорости u=f(y) в различных сечениях основ­ ного участка струи различны: чем дальше отстоит сечение от среза сопла, тем шире поле и меньше ит (см. рис. 17.1).

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

точках любых сечений

/,

2 yi/Rr[A =

безразмерные

скорости

равны

Ul~~и* =

U2~~UH ,

т. е. безразмерные поля

скорости

и

Ыт1 — ин

ит2—ин

 

 

и

 

и.

Аи

* / /

\

для

 

 

---------- = ------= J (УЩп)

всех сечении совпадают и описыва-

Ит

Нн

^и т

 

 

 

 

 

 

ются формулой Шлихтинга

 

 

 

 

 

 

 

_и - и н

= ( 1 _

Т1з/2)21

(17. 2)

 

 

 

 

“ т — и н

 

 

 

где ын, «т, и — скорости: спутного потока, на оси струи, на рассто­ янии у от оси струи в данном сечении, соответственно, r\ = ylRrp — расстояние от оси струи до точки со скоростью и, выраженное в до­ лях RTV, или полуширины Ьгр данного сечения струи.

Асимптотичность пограничного слоя (см. рис. 17.1) затрудняет экспериментальное определение его толщины Rrp или бгр. Поэтому, обычно, нормирование ординаты поля ведется не по Rгр или Ьт а

по ординате ус при которой

— = 0,5. Ордината ус достаточ-

но точно определяется по экспериментальным

данным,

a yJRjv

рассчитывается по формуле (17.2)

 

 

 

yJRT? = 0AM

Rrp = 2,27ye

 

(17.3)

и расчетная формула (17.2) приобретает вид

 

 

Поле скоростей Au/Aum= f(yfyc) приведено на рис. 17.2.

и (17.4)

Опыты показывают, что формулы Шлихтинга

(17.2)

применимы для самых различных

струй: осесимметричных, пло­

ских, затопленных, в спутном потоке и в противотоке, подогретых и охлажденных (0,25<QH/eo<4), для течений в турбулентном слое

ззо

за плохообтекаемым телом

и т. д.

 

 

 

 

 

 

Для

всех

этих

случаев

безразмер­

Ч

ч

 

 

!

 

ные

поля

скорости

универсальны,

 

 

 

I

 

V s \

 

 

 

т. е. не зависят

ни от

числа Рей­

 

^

\

 

 

 

нольдса Re0 ни от числа Маха М0,

0,75

\

\

 

 

 

 

v

 

 

 

ни от других индивидуальных

осо­

 

 

V v V / / ^/

 

 

бенностей струи.

 

 

 

 

 

0,5

 

 

 

 

 

 

 

 

 

\

\

 

 

На основном

участке

течение

 

 

 

 

 

йи/АитЪк

\

\

 

приобретает такую

же

структуру,

 

 

 

 

 

 

 

\ \

 

как при истечении

из точечного

ис­

0,25

 

 

\

\

 

точника Ro-+-0 или из плоской щели

 

 

 

 

S

 

 

 

 

 

ч \

2Ь0-+0 — из полюса струи

(см. рис.

 

 

 

 

 

 

17.1). Полюс струи может распола­

 

0 ,5

1,0

1,5

2 ,0

у /у с

 

 

 

 

 

 

гаться как елева, так и справа от

Рис.

17.2.

Безразмерные

поля

среза сопла, на расстоянии х0

или

температуры и скорости на основ­

совпадать с ним, в зависимости

от

ном

участке струи

 

 

 

условий истечения

(рис.

17.3).

 

 

 

 

 

 

 

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

Б е з р а з м е р н о е п о л е с к о р о с т е й в п о г р а н и ч н о м

слое н а ч а л ь н о г о у ч а с т к а с т р у н е

небольшой погрешно­

стью можно рассчитать по формулам (17.2),

(17.3) и (17.4), в кото­

рых вместо ит следует подставлять и0, а вместо RTр или ЬТр — тол­ щину пограничного слоя в данном сечении начального участка, т. е.

расстояние от

границы

ядра постоянной

скорости до грани-

цы струи.

 

 

 

П о д о б и е

поле й

т е м п е р а т у р ы и

к о н ц е н т р а ц и и

и з б ы т о ч н о й

примеси . Подобие полей скорости в поперечных

сечениях пограничного слоя струи предопределяет подобие полей температуры и концентрации.

Опыты показывают, что при одинаковых граничных условиях

Ц

цн_

7’

7*н __

£

_0 01

и т

Ин

Т т

Т н

с m

сн

тепловой и диффузионный пограничные слои толще динамического #гр.т=^гр.с>#гр.и, что соответствует более интенсивному турбулент­ ному переносу тепла и вещества по сравнению с переносом количе­ ства движения (см. рис. 6.5). Это значит, что турбулентное число Прандтля меньше 1.

Для осесимметричной струи 7?Гр.т = # г р .с ~ 1,18/?гр.„ и Ргт»0,8'. Для плоскопараллельной струи Rrp.T = Rrp.cS:: l,45#rp.u и Ргт«

«0,5.

При этих условиях безразмерные поля скорости, температуры и концентрации сливаются

ди

= /

(у/%гР.и)= т £ - = / (У/Ягр.т)

Ас

д«т

f{y!RTVJ-

 

Л* тп

Аст

Однако, для упрощения исследования обычно принимают, что границы струи по температуре и концентрации совпадают с грани­ цами по скорости Ягрл= ^.ро—Ягр.и=Ягр.и(х)- При этом граничные условия несколько изменяются

В этом случае безразмерные поля температуры и концентрации хорошо аппроксимируются формулой

7*

Т н

с Сн

|

/ у \3 /2

Тт

7'н

ст — с„

 

\2,27ус!

(Л ^ин _ у'г

п 7 .5)

ит Ин '

 

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

постоянных температур и концентраций

короче

ядра постоянной

скорости

и безразмерные

температура

Т _т

и

концентрация

— ------

Q __ Q

 

 

 

7"о — 7’н

 

 

вдоль оси х уменьшается быстрее, чем

безразмерная

-----51

с 0 — Сн

 

 

 

 

 

 

скорость

Um

и" (см. рис.

17. 1).

 

 

 

 

U.Q

и н

 

 

 

 

Уравнения (17.2) ... (17.5) замечательны тем, что не зависят от индивидуальных свойств струи и позволяют определить искомые поля скорости, температуры, плотности и концентрации, если оире* делены границы струи Rrp.u=Rт р . и ( х ) или Ьтр.и= Ьг р,и ( х ) и значения параметров на оси струи ит, Тт, ст (в изобарных струях поле плотности определяется полем температуры Q=PlRT).

Итак, струйные течения обладают свойствами, характерными для пограничного слоя: поперечные размеры течения малы по сравнению с продольными, поперечная скорость мала по сравИе*

нию с продольной v^tiy поперечный градиент продольной скорости больше продольного ди/ду^ди/дх, безразмерные поля скорости, температуры и концентрации универсальны.

О с о б е н н о с т и с в о б о д н о й т у р б у л е н т н о с т и . В при­ стеночном турбулентном пограничном слое стенка гасит пульсации. Поэтому величина пути смещения или масштаба турбулентности в поперечных сечениях пропорциональна расстоянию от стенки /= = ку = 0,4у (см. п. 8.1). Струйный турбулентный пограничный слой не ограничен стенками. Поэтому величина пути смешения в дан­ ном сечении имеет постоянное значение. Из подобия полей скоро­ сти следует, что отношение характерных линейных размеров сече­ ний вдоль оси х величина постоянная, т. е. величина пути смешения пропорциональна толщине пограничного слоя

//# Гр = const; l/br? = const.

(17.6)

Опыты показывают, что свободная турбулентность имеет двоя­ кую структуру. Основная часть пульсаций имеет сравнительно ма­ лый масштаб и высокие частоты от нескольких килогерц до 200 Гц и содержат основную часть турбулентной энергии. На эту структу­ ру налагается система больших вихрей с частотой пульсаций по­ рядка 20.... 30 Гц. Расширение свободных турбулентных струй оп­ ределяется движением этих вихрей, для которых справедлива за­ висимость (17.6). Большие вихри искривляют границы погранич­ ного слоя с ядром постоянной скорости и с окружающей средой и осуществляют захват нетурбулентной жидкости. Эта модель пред­ полагает наличие сравнительно резкой границы между турбулент­ ной и нетурбулентной жидкостью, что подтверждается опытом. В тонком слое, в месте соприкосновения турбулентной и нетурбу­ лентной жидкостей, должна проявляться вязкость, так как пере­ дача завихренности может происходить только за счет сил сдвига. Этот тонкий слой называется ламинарным надслоем, по аналогии с ламинарным подслоем в турбулентном пограничном слое на твер­ дой поверхности. Очевидно, что в области границ струйного погра­ ничного слоя течение имеет перемежающийся характер, так как через данную точку пространства хаотически во времени проходят моли жидкости различной степени турбулентности. На рис. 17.1 со­ поставляются поле скорости и коэффициент перемежаемости у (см. п. 6.1) в сечении основного участка струи. Вблизи оси струи коэф­ фициент перемежаемости равен единице, а в области границы он резко падает до нуля. Характерно, что ширина струи, определенная по пульсациям скорости, т. е. по у, всегда превышает ширину, опре­ деленную по осредненной скорости. График распределения степени турбулентности е и = и'/ит по сечению основного участка струи по­ казывает неравномерность этого распределения. Максимум интен­ сивности примерно соответствует максимуму du[dy.

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

ззз

ского струйного пограничного слоя пропорциональна пульсацион- ■ной составляющей поперечной скорости,

db

db

d x

db

___

(17. 7)

dt

d x

d t

d x

 

 

 

По (6.16) v'= u' Ttilduldy. Ввиду подобия полей скорости гради­ ент du/dy во всех сечениях пограничного слоя пропорционален мак­ симальной разности скоростей, отнесенной к его толщине dufdy~ ~ т —ии)/Ь, где ит и ын скорости на внутренней и внешней гра­

ницах пограничного слоя, тогда

 

 

v '---- —1 (ип —ии), а с учетом (17. 6) — v' ~ ( « m— ы„).

 

П

 

 

 

Подставляя значение v'

в (17.7), получим, что нарастание

толщи­

ны пограничного слоя по длине струи

 

db

___ v ' _

\и т и н[

^ | у

d x

и

\и\

 

пропорционально интенсивности турбулентности v'/u, характерной для всего данного сечения пограничного слоя. Так как v'ju всегда положительно, то струя расширяется — db/dx>0. При этом величи­ на v' является фактором, увеличивающим Ъ во времени, а и — уменьшающим, сносящим приращение Ь вдоль течения.

Х а р а к т е р н у ю с к о р о с т ь и в уравнении (17.8),

с учетом

■сжимаемости M= QH/QO^ 1, рационально определять как среднемас­

совую для данного сечения пограничного слоя. Однако,

скорость

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

ьь

«С1,= \Q4dy/^Qdy.

(17.9)

о

6

 

Вследствие того, что и и Q изменяются как вдоль, так и поперек

■струи, уравнение (17.8) с использованием иср получается

чрезвы­

чайно сложным не только для решения, но

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

анализа.

 

 

Избыточные скорости Дит = ит—мн, температура АТт = Тт—Гп,

а следовательно и плотность рт QH по длине

основного

участка

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

нута. Поэтому сложную зависимость (17.9)

для характерной ско­

рости в сечении можно заменить упрощенной:

и = Qm^-m ~b QU^ H

(17. 10)

Qm-b QI

 

Подставляя удвоенное значение и из (17.10) в (17.8), заменяя знак пропорциональности знаком равенства и константой С, получим за­ кон нарастания толщины струи

d b __ £

и т

и»________ Qm

с -

1—

(17.11)

dx

2

,

QH

2

 

 

 

ит + « н

Qm

Qm Mm

 

 

 

 

Здесь константа С может быть определена только эксперимен­ тальным путем.

А н а л и з з а к о н а н а р а с т а н и я т о л щ и н ы п о г р а н и ч ­ ного с ло я струй.

1. Затопленная изотермическая струя ын = 0, Т0 = ТВ и Q0 = (>H-

Для этого простейшего случая из (17.11) получим линейный за­

кон расширения струи

 

V 1 = V 1 = C ; Яз=Сх\ Ь3=Сх.

(17.12)

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

Для

начального

участка

0,27;

ан^15°30'.|

^ ^

Для

основного

участка

С ^0,22;

12°30' J

 

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

Эксперименты показывают, что полюс, в котором пересекаются прямолинейные границы основного участка изотермических затоп­ ленных струй, практически совпадает со срезом сопла (х0 = 0). По­ этому отсчет абсцисс сечений, как для начального и переходного, так и для основного участков, в формулах (17.12) следует произво­ дить от среза сопла (см. рис. 17.3).

Расчеты струй, с использованием интегрального уравнения ко­

личества движения и значений константы С (см. п. 17.2), показы­

вают, что вследствие пространственности течения, начальный и пе­

реходный участки осесимметричной струи короче соответствующих

участков плоской при одинаковых

законах расширения границ.

В реальных течениях прямолинейные границы начального и основ­

ного участков плавно сопрягаются

криволинейной границей пере­

ходного участка.

 

2.

Два полубесконечных потока при Q0 = QH характерны тем, что

на границах пограничного слоя um = u0 = const и ип = const, поэтому

границы прямолинейны (см. 17.11)

 

-**- = С

= const и b= const х.

 

dx

\ит\ + |и„|

Сопоставляя эти выражения с (17.12), получим

Ь _

!ц ч ~ “ иI

(17. 14)

h

|u.„| + 1«„|

 

При спутном движении двух полубесконечных потоков ит и «н

имеют одинаковые знаки. Поэтому, с увеличением

спутности т=

= и„/и0, утолщение пограничного слоя уменьшается. В пределах

0 < т < 0 ,5 эксперименты подтверждают закономерность

(17.14),

что указывает на стабилизирующее действие спутности,

приводя­

щее к уменьшению тангенциального разрыва скорости и к сниже­ нию интенсивности турбулентности (см. рис. 17.1). При 0 ,5 < т ^ 1 величина 6/63» 0,3 и формула (17.14) не применима.

Задача 17.2. Объясните, почему возникает турбулентный пограничный слон при т = 1 „ когда тангенциальный разрыв скорости формально отсутствует.

При распространении струи во встречном потоке, скорости на границах пограничного слоя имеют разные знаки, тогда из (17.14) имеем

(17. 15)

Ь3 ит + ин

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

3. Струя конечной толщины в спутном потоке ( QO = QH) -

В на­

чальном участке um=u0=const и un=const. Поэтому все

выводы

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

Для основного участка формула (17.11) принимает вид

 

= c\Um~ ц"1-.

 

(17.16)

dt

ип + и„

 

 

При цн<ыо величина осевой скорости

um вдоль оси основного

участка уменьшается и стремится к ив.

Вследствие этого

db[dx

вдоль оси также уменьшается и граница струи в спутном

потоке

криволинейна (рис. 17.4).

4. Затопленная струя, плотность газа которой не равна плотно­ сти газа в окружающей среде гс = ен/бо^1- Разница в плотностях может быть достигнута за счет разницы в температурах, так как для изобарной струи Qn/Qo — To/TH, либо за счет использования раз­ личных газов. При и,,= 0 формула (17.11) принимает вид

Если струя подогрета 0= 7’о/7’н>1 и ен/(>о>1, то она расширяет­ ся в большей степени, чем изотермическая. Длины начального и пе­ реходного участков сокращаются, границы их прямолинейны. В ос-

Рис. 17.4. Границы осесимметричных за­ топленных струй при m=UiJuo<\

новном участке gm увеличивается вдоль оси, приближаясь к дн, поэто­ му границы криволинейны. Однако, на расстоянии более 50 Ь0 от сопла, они приобретают направление, па­

раллельное

границам

изотермичес­

кой струи

С= 0,22. Холодная струя

0 = 7о/7,п< 1 > QH/QO<1

напротив, расширяется медленнее изотер­

мической, но к х>50Ь0 ее границы также принимают направле­

ние, параллельное границам изотермической затопленной

 

струи.

Дальнобойностью струи называется расстояние от среза

сопла,,

при котором скорость на оси достигает половины исходной

ит =

= 0 ,5 w 0-

Задача 17.3. При одинаковых Ь0 и и0 качественно сравните границы и даль­ нобойность следующих струй: затопленных — изотермической, подогретой, ох­ лажденной; изотермических — затопленной и в спутном потоке.

17.2. РАСЧЕТ СТРУЙ

Для решения поставленной задачи, т. е. для определения пара­ метров газа в струе и ( х , у ), Т(х, у), Q ( X , у) и с(х, у), кроме уже определенных границ струи (17.12) и (17.13) и универсальных про­ филей параметров в поперечных сечениях слоя смешения (17.4) и (17.5), необходимо определить изменение этих параметров вдоль оси основного участка струи Um(x), Tm(x), дт (я) и ст(х). Для это­ го используем следующие основные уравнения газовой динамики в интегральной форме.

1. Уравнение количества движения. Выделим произвольный ко­ нечный участок изобарной осесимметричной затопленной струи. Напряжения трения на контрольной поверхности равны нулю — от­ сутствуют поперечные градиенты скорости. Силы давления уравно­ вешиваются вследствие изобарности струи. Таким образом, проек­ ция на ось х суммы сил, действующих на элемент, равна нулю и уравнение выражает постоянство количества движения в любом се­ чении струи:

^гр

Qo«oJi/fl)= f §u22nydy,

(17. 18)

b

 

где у — текущий радиус точки сечения, в которой

скорость равна

и, а плотность — д.

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

Заменим y/Rrp = i\; d(y/RTp)= dr\;

=

JL = Д_Л* ..

Учитывая, что Rrp/Ro и um/uo

кгр Ко

U0

ani u0

не зависят от у и г\ и могут быть

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

 

 

 

 

1.

(17. 23)

Подставляя значение и/ит из (17.2), получим, что ^ ^—~)~ г1^г1=

о

= 0,067 и тогда из (17.23) найдем связь между безразмерным ра­ диусом поперечного сечения и безразмерной скоростью на оси это­ го сечения

Rr.JR0=2,72u0/um. (17.24)

Радиус переходного сечения определим из условия ит=ио

R„= RJRo= 2J2.

(17.25)

Длина переходного участка хп определяется из условий

x J R 0 = RJCRo и R JR 0 = 2 J 2 -x ,= R J C ^ 2 J 2 /0 ,2 2 ^ V 2 A .

(17. 26}

Абсцисса начального участка в соответствии с эксперименталь­ ными и расчетными данными [2]

x H= x J R 0zs8.

(17.27)

Отсюда угол между образующей ядра постоянных

скоростей и

осью струи ая~7°. Для плоской струи хи= xBfb0« 9 и ая~ 6°30' (см. рис. 17.3).

Граница основного участка осесимметричной_ и плоской струй

определяется по (17.12) и

(17.13)

одинаково Л = 0,225, 5=0,225.

Осевая скорость газа um=Um {x)

в основном участке осесиммет­

ричной струи определим, используя

(17.24), (17.12)

и (17.13)

 

Лв-=

12;1_.

(17.28)

 

и0

X/RQ

 

Подобный расчет для плоской струи дает

 

=

-М = г

и х„== 14,4.

(17. 29)

и 0

у

х / ь 0

 

 

 

Сопоставление (17.28)

и

(17.29)

подтверждает

уже отмечав­

шееся более быстрое затухание скорости в осесимметричной струе. Объемный расход газа через поперечные сечения основного уча­

стка осесимметричной струи

 

 

^гр

п

2 1

Q= \ u2nydy=2nRoU0-^IIL 1-^ -)

\

J

“ 0 \ # 0 J j и т

В этом уравнении KRQ2U0= Qo — расход газа через сопло, 1

(RJRof = 7A(u0fum и \ - ^ т|£/л = 0,128;

JО«ш

Q/Q0= l,9 u 0/um. (17.30)

Сопоставление (17.30) и (17.28) показывает, что расход газа на основном участке осесимметричной струи возрастает пропорцио­ нально х. В конце начального участка ит=ио и QH=l,9Qo-

Изменение температуры газа вдоль оси основного участка осе­ симметричной струи определим, выполнив преобразования уравне­ ния сохранения избыточной энтальпии (17.20), аналогично преоб­ разованию уравнения сохранения импульса.

При n = QH/Q0— 1, 0 = 7’о/Тн —>1 и Ср= СРо> получим

“от-АГт(— f 2 Г— — T|rfTl=l,

и0 ЛГо \ #0 /о ЬТтJ Цщ

где, в соответствии с (17.5), при РгТ = 0,8

АТ

дт„

Принимая во внимание подобие полей температур и концентра­ ций, найдем

Д^m

ДСm

Q gg

Um

10,9

(17.31)

ДТ’о

Дсо

Щ

x/R0

 

Полученные формулы дают возможность рассчитать осесиммет­ ричную и плоскую затопленные изобарные струи только для при­ нятых идеализированных условий, при которых определяющим активным параметром является только начальный импульс ео«о2яЯ02.

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

Спутность т = ив/и0 и отношение плотностей n = QulQo существен­ но влияют на границы струи и изменение параметров в слое сме­ шения.

Увеличение плотности активной струи, например за счет ее ох­ лаждения 0= 7’О/7Н< 1 , приводит к увеличений длины начального участка струи и к соответствующему смещений кривой, характери­ зующей затухание осевой скорости, вдоль осН основного участка (рис. 17.5).

Аналогичное влияние увеличения спутностИ показано на рис. 17.4.

Особенно существенное влияние на струйное течение оказывает

неравномерность распределения параметров сТРУи на срезе сопла, обусловленное пограничным слоем на его стен№х- Эта неравномер-

ность

изменяет

начальные

интег­

 

ральные

параметры

струи: рас­

 

ход,

импульс,

энтальпию,

что

 

обычно

в расчетах

учитывается

 

введением коэффициентов, оцени­

 

вающих

эту

неравномерность.

 

Кроме того, неравномерность про­

 

филя скорости приводит к повы­

 

шению уровня турбулентности и

 

к сокращению длины

начального

 

участка.

 

 

 

 

 

 

При истечении в спутный по­

Рис. 17.5. Изменение безразмер­

ток действие внутреннего и внеш­

ной скорости вдоль оси струи в

него

пограничных слоев

сумми­

зависимости от начального подо­

руется и усиливается. Этим объ-

грева

ясняется развитие конечного тур­

 

булентного слоя

смешения

п р и т = н н/ио=1*

Наконец, существенным фактором, с помощью которого можно управлять процессами турбулентного смешения, является искусст­ венная турбулизация струй, которая становится эффективной при е0>5% .

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

циссы переходных

сечений соответствующих

параметров

хии, *пт,

хпс, то изменение этих

параметров

вдоль

оси осесимметричной

струи может быть рассчитано по следующим формулам

 

 

& иГп

__ х \\и

. А7*т ___

-Уц.т .

т ___ -^п.с

 

(17.32)

 

Дно

х

Д7\)

 

х

 

ДСо

х

 

 

 

 

 

 

Для плоско-параллельной струи

 

 

 

 

 

 

Дн0

1

/ ~ х \ш

. Д7*т

у

Г

-^нт

. ^ с т

Г

х ис

(17.33)

у

х

ДГо

 

х

ДСо

| /

х

 

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

Расчет струй при более сложных начальных условиях можно найти в специальной литературе.

характеристики прямоугольного крыла бесконечного размаха. На всех участках единичной длины такого крыла действуют одинако­ вые полная аэродинамическая сила R, подъемная сила Ry, нор­ мальная к вектору скорости на бесконечности и сила лобового соп­ ротивления RXi совпадающая по направлению с направлением век­

тора скорости на бесконечности, т. е. исключаются

концевые эф­

фекты, присущие крыльям или лопаткам конечных размеров.

В соответствии с положениями теории

подобия

используются

не только силы и моменты сил, но и их безразмерные

коэффици­

енты

R.,

 

 

м

 

с„

 

С„

,__• , го х

б

w*

Qoo^o

 

 

*+ОО 01

Ьг

где R, Rv, Rx — силы, действующие на один погонный метр крыла бесконечного размаха Н/м; М — момент силы R относительно выб­ ранной оси; z — характерная длина.

Отношение величины подъемной силы профиля к силе лобового сопротивления называется качеством профиля

K= Ry/Rx = Cy/Cx.

Задача 18.1. Покажите, что качество есть отношение проекции па горизон­ таль длины планирования профиля к высоте, с которой начинается его планиро­

вание. Подсчитайте, на

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

1000 м планер

(К =38), пассажирский

самолет (К=10) и сверхзвуковой

лайнер (К = 6).

В общем случае Су, Сх, Ст зависят от геометрии

профиля или

других тел, угла атаки и от критериев подобия:

чисел Re, М и др.,

атакже от взаимодействия исследуемых тел с соседними телами.

Вдинамически подобных системах Су, Сх, Ст и К одинаковы.

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

Теорема Жуковского о подъемной силе (п. 4.9), постулат Жу­ ковского—Чаплыгина (см. ниже) с использованием метода кон­ формного отображения (см. п. 3.10) позволяют определить вели­ чину Ry и Су теоретически.

В течение ряда лет, после получения Н. Е. Жуковским формулы подъемной силы Rv= д001Т'с0Г, ученые и инженеры не могли ее ис­ пользовать при проектировании крыльев, винтов самолетов и лопа­ ток турбомашин, так как не были известны причины самопроиз­ вольного появления циркуляции скорости Г (см. п. 3.5) вокруг про­ филя и методы расчета ее величины. В 1908 г. Н. Е. Жуковский и С. А. Чаплыгин сформулировали свой знаменитый постулат.

П о с т у л а т Ж у к о в с к о г о —Ч а п л ы г и н а . При безотрыв­ ном обтекании профиля, вокруг него возникает циркуляция скоро­ сти Г такой величины, при которой струи плавно стекают с задней острой кромки с конечной скоростью. Постулат — необходимое до­ полнение к теореме Жуковского. Он определяет причину возникно­ вения и величину циркуляции скорости вокруг профиля.

В о з н и к н о в е н и е ц и р к у л я ц и и с к о р о с т и и р а с п р е ­ д е л е н и е д а в л е н и я по пр о фил ю. Циркуляция скорости и подъемная сила равны нулю при обтекании симметричного профи-

Рис. 18.2. Бесциркуляционное об­

Рис. 18.3. Обтекание профиля с

текание профиля:

циркуляцией скорости

а—при нулевом угле атаки; б—при }г-

 

ле атаки а

 

ля при нулевом угле атаки а= 0

(рис. 18.2). В этом случае задняя

острая кромка совпадает с задней критической точкой или точкой схода струи. Суммарные силы давления на нижнюю и верхнюю по­ верхности профиля одинаковы.

При увеличении угла атаки а> 0 передняя критическая точка К перемещается вниз по профилю к точке Б, а точка схода струй В — вверх по профилю к точке А (рис. 18.2, б). Если было бы возмож­ но плавное обтекание жидкостью задней острой кромки, то устано­

вилось бы новое бесциркуляционное обтекание

профиля — Г=0,

Ау= 0. Однако, при радиусе закругления острой

кромки г—0, ско-

рость безотрывно обтекающей жидкости должна беспредельно воз­ растать (см. п. 3.8), так что давление, вычисленное по уравнению Ьернулли, должно было бы неограниченно уменьшаться /?—>•(—°°), что невозможно. В действительности на верхней поверхности про­ филя самопроизвольно возникает течение жидкости к задней кри­ тической точке, где давление понижено. Это течение возвращает точку схода струй в заднюю острую кромку профиля. При этом поток жидкости срывается с острой кромки в виде начального или разгонного вихря, вращающегося против часовой стрелки с цирку­ ляцией (—Г) и сносится потоком (рис. 18.3).

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

должна остаться равной нулю r s = 0, как бы далеко не

уносился

потоком разгонный вихрь. Но это возможно только в том

случае,

если при срыве разгонного вихря вокруг профиля установится цир­ куляция, равная циркуляции разгонного вихря и направленная в об»

ратную сторону. При наложении такой циркуляций на гипотетиче­ ское бесциркуляционное течение (см. рис. 18.2, б) получаем реаль­ ное безотрывное обтекание профиля, при котором скорости над профилем становятся больше, а под профилем меньше скорости не­ возмущенного потока. В соответствии с уравнением Бернулли дав­ ление под профилем повышается, а над профилем — понижается, что приводит к возникновению подъемной силы.

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

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

О п р е д е л е н и е

ч и с л о в о й

в е л и ч и н ы

ц и р к у л я ­

ции с к о р о с т и

д л я з а д а н н ы х

п р о ф и л я

и у с л о в и й

о б т е к а н и я (роо, W7*,, а ). На основании постулата Жуковского— Чаплыгина определяется точное положение передней и задней кри­ тических точек на профиле в физической плоскости г. Отобразив профиль и течение около него на вспомогательную плоскость £ круга (см. п. 3.10), находим соответствующие критические точки и поток, что и определяет искомую циркуляцию, одинаковую для круга и профиля в обеих плоскостях и определяемую выражением

Г = mJbW* sin (а — а0)= mub W sin ал,

(18. 1)

где ти — коэффициент, приближенное значение которого для про­ извольного профиля имеет вид

 

ша = л (1+ 0.77С )]/ 1 + (7/2)2-

(18.2)

 

Подъемная сила для профиля единичной длины, Н/м

 

 

Ry=QW„T = mubQWl sin аА

(18. 3)

и

С>/ = 1 ^ — = 2/rau sin аА.

(18.4)

 

- F - *

 

 

Для обычно применяемых малых углов атаки sin а л ^ а л

и

 

Ry= mubQW ia A;

(18.5)

 

Су = 2тиаА.

(18.6)

12

950

34.5

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

Задача 18.2. Объясните, почему для обычных самолетов необходимы доста­ точно длинные взлетно-посадочные полосы.

П о д ъ е м н а я с ил а п л о с к о й п л а с т и н ы п р и б е з ­ о т р ы в н о м д о з в у к о в о м о б т е к а н и и и д е а л ь н о й ж и д ­ костью. Полагая в (18.2) С= 0 и / = 0, получим

Су = 2 п а .

(18.7)

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

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

Задача 18.3. Обоснуйте форму профиля,

приведенного_на_рис. 18.1: закруг­

ление передней и заострение задней кромок,

изогнутость f и СФ 0.

при Q=

Задача 18.4. Человек,

массой т = 80 кг поддерживается в воздухе

= 1,2 кг/м3 воздушным змеем, который буксируется со скоростью №„ =

100 км/ч.

Определить площадь змея

при угле атаки а= 3 °, если

он выполнен:

а)

в виде

профиля С=10% ; /=2%;

ао=3°; б) виде плоской

пластины. Ответ:

Sa =

= 2,34 м2; S <5= 5 м2.

 

 

 

 

 

18.2. ВЗАИМОДЕЙСТВИЕ ТЕЛ С ПОТОКАМИ РЕАЛЬНОЙ ЖИДКОСТИ

Возникновение подъемной силы при обтекании тел реальными жидкостями определяется теоремой Жуковского и постулатом Жу­ ковского—Чаплыгина так же, как в уже изученном обтекании тел идеальной жидкостью.

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

подъемной силы и силы лобового сопротивления R = R y+Rx- Силы лобового сопротивления крыла или лопатки турбомашины

состоит из сопротивлений профильного /?*р, волнового Rxb и индук­

тивного Rxi

 

 

т. е.

 

(18. 8)

 

 

П р о ф и л ь н о е

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

это сопротивление бес­

конечно длинного

крыла при обтекании

его с дозвуковой скоро-

стью. Оно слагается из сопротивлений трения и давления

Сгр= С*тр+ С*д.

(18- 9)

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

Для того, чтобы найти распределение давления по поверхности

профиля, обтекаемого реальной жидкостью,

следует рассмотреть

обтекание идеальной жидкостью профиля,

все сечения

которого

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

слоя б* =

= 6* (я) [см. п. 15.1]. В этом случае распределение давления на пе­ реднюю часть профиля практически не изменится из-за малого зна­

чения б и б* на этом начальном

участке

конфузорного течения.

В кормовой диффузорной области

течения

б и б* увеличиваются

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

Задача 18.5. Нарисуйте схемы распределения статического давления по про­ филю при обтекании его идеальной и реальной жидкостью и докажите, что в первом случае R x д = 0 , а во втором — Rx д > 0 .

Для определения силы Rx лобового сопротивления тела приме­ ним теорему количества движения к контрольному объему жидко­ сти (см. п. 4.1), включающему обтекаемое тело (рис. 18.4) * При*-

Рис. 18.4. Лобовое сопротивление профиля Сх

* Формула подъемной силы Ry= роо^«>Г также может быть получена из теоремы количества движения.

менение теоремы количества движения основано на том, что тело действует на жидкость с силой, равной силе лобового сопротивле­ ния, но противоположно ей по направлению, что вызывает соответ­ ствующее этой силе уменьшение количества движения жидкости. Это проявляется в существовании за обтекаемым телом аэродина­ мического следа, т. е. спутной струи. Частицы жидкости, подтормо­ женные в пограничном слое, попадают в этот след и постепенно, обычно турбулентно, смешиваются со сносящим потоком. Границы следа расширяются, а скорость на оси увеличивается, приближаясь в бесконечности к скорости невозмущенного потока (см. гл. 17) *. Непосредственно за плохообтекаемым телом давление в следе мо­ жет быть существенно ниже давления в невозмущенном потоке. Однако, оно очень быстро выравнивается, в то время как отличие в скорости остается на больших расстояниях за телом.

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

Применив для жидкости в контрольном объеме ABCD интег­ ральное уравнение количества движения (4.11) и учитывая, что расход жидкости через сечения АВ и CD одинаков, получим

+j p ~ d y - +\ p d y - R x= )' {Q„W „W -Q „W l) d y

или

+ 00

+ oo

RX= Q„W„ j

{W „ -W )cly+ J (p „ - p ) d y . (18.10)

—oo

— oo

Формула (18.10) показывает, что возникновение силы лобового сопротивления объясняется рассеянием механической энергии в потоке вследствие вязкости жидкости, т. е. возрастанием энтропии. Для определения Сх необходимо знать распределение параметров потока на контрольной поверхности.

Задача 18.6

Определить Сх профиля (см.

рис. 18.4), если хорда Ь = 2 м;

6 = 0,2 м; №оо=100 м/с; 1^т = 90 мс. Принять,

что распределение скоростей ли­

нейное, а давление на всех границах постоянно р = р 00. Ответ: С.х= 0,02 .

Хорошо-

и п л о х о о б т е к а е м ы е тела. Тела, обтекаемые

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

Сопротивление тонких профилей при дозвуковом безотрывном обтекании обусловлено почти исключительно трением. При малых Re^<Re.x кр и нулевом угле атаки сопротивление тонкого профиля практически равно сопротивлению тонкой пластины при ламинар­

* Формула поля скоростей (17.2) получена Шлихтингом для дальнего следа.

ном пограничном слое. При уве­ личении Re*>ReXKp сопротивле­ ние увеличивается за счет появ­ ления турбулентного погранично­ го слоя. Округление передней кромки позволяет избежать от­ рыв потока при изменении угла

Рис. 18.5. Зависимость •сопротивлений трения и давления от относитель­ ной толщины профиля v (n O T O K дозвуковой)

Рис. 18.6. Экспериментальная кривая Су(а) для единичного профиля

атаки от нулевого. Заострение задней кромки обеспечивает плав­ ный сход струи с нее в определенном диапазоне углов атаки и по­ лучение подъемной силы, а также малые dp/dx в кормовой части профиля и безотрывное его обтекание. _

С увеличением относительной величины С профиля увеличива­ ется диффузорность течения в его кормовой части и сопротивление давления (рис. 18.5).

Зависимость Су и Сх несимметричного профиля от угла атаки а приведена на рис. 18.6. Для несимметричного профиля угол нуле­ вой подъемной силы отрицателен (схо = —7°). При а< ао подъем­

ная сила направлена вниз

Ry<CО,

0. Вначале

Су возрастает

пропорционально а, что качественно соответствует

теоретической

зависимости (18.6), однако

измеренные

значения

Су получаются

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

рой кромки *. Любой из рассмотренных ранее методов ликвидации отрыва пограничного слоя (см. п. 15.6) может обеспечить увеличе-»

ние акр и Су.

К о э ф ф и ц и е н т л о б о в о г о с о п р о т и в л е н и я п л о х о ­

о б т е к а е м ы х

тел.

Рассмотрим зависимость

Cx= f(Re) на при­

мере обтекания

шара

при Мсзо= 0 и постоянной

интенсивности тур­

булентности Еоо

невозмущенного потока (см. рис. 5.2).

I. При малых числах ReooClOO отрыв пограничного слоя отсут­

ствует, картина обтекания шара близка к картине обтекания шара идеальной жидкостью и сила лобового сопротивления является почти исключительно силой сопротивления трения. Резкое сниже­ ние Сх с увеличением Reoo показывает, что в этой области сила сопротивления пропорциональна скорости №«>, что характерно для

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

II. При 100<Reoo<2,5-103 в кормовой части возникает неустой­ чивое вихревое движение. Это понижает давление в кормовой час­

ти и приводит к увеличению сопротивления давления

и замедле­

нию падения Сх с увеличением Re^.

 

 

пограничный

III. В области 2,5*

I03< R e(XJ<2* 105 ламинарный

слой отрывается от поверхности шара. Положение точки

(линии)

отрыва ламинарного

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

слоя

не зависит от

Re^ и

происходит в одной и той же точке

S при

ф»80°

при

dp/dxX

X (б*)2/цоо1^оо= 1,92 (см. п. 15.6), т. е. при dp/dx>0. Это указывает на то, что при обтекании шара вязкой жидкостью, точка минималь­ ного давления М располагается, вследствие воздействия погранич­ ного слоя на поток, при ф»70°, а не при 90°, как в случае идеаль­ ной жидкости. Ламинарный пограничный слой, оторвавшийся от поверхности, продолжает на некотором расстоянии £—Т течь в ви­ де ламинарной струйки. Образовавшаяся поверхность тангенци­ ального разрыва скорости неустойчива, поэтому в сечении Т возни­ кает интенсивная струйная турбулентность. В результате интенсив­ ного турбулентного обмена жикость подсасывается оторвавшимся (активным) потоком из кормовой области, где давление понижает­ ся до величины р = ( р —роо)/(е оо Woo/2) ^ —0,4 аналогично тому, как это происходит при поперечном обтекании цилиндра, где рж —1,2 (см. рис. 4.13). На расстоянии одного-двух диаметров шара давле­ ние в среде сравнивается с давлением в окружающей среде /) = 0 и из этой области возникает встречный завихренный поток к кормо­ вой поверхности шара — зона обратных токов.

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

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

тоянное значение для шара Cs = 0,45 ... 0,47 (для цилиндра С, «1,3) и сила лобового сопротивления Rx~ W^2. Область III называется областью локальной автомодельности по числу Рейнольдса.

В этой области происходит перестройка структуры обтекания. При увеличении числа Рейнольдса Re» увеличивается интенсивность турбулентности в области обратных токов и точка Т перехода ла­

минарного течения в турбулентное приближается к точке 5 отрыва ламинарного пограничного слоя.

Рис. 18.7. Степень тур­ булентности в зависимо­ сти от критического чис­ ла Рейнольдса

IV.

При Re,x>«3,2-105 наблюдается кризис сопротивления, при

котором Сх уменьшается в 3 ... 5 раз при незначительном увеличе­

нии числа

Рейнольдса. Это интересное явление объясняется

тем,

что при Re<x> кр точка перехода Т совпадает с точкой отрыва

5 ла­

минарного пограничного слоя, т. е. ламинарный пограничный слой

турбулизуется перед его отрывом. Турбулентный пограничный слой обладает большей сопротивляемостью отрыву. Поэтому точка от­ рыва, теперь уже турбулентного пограничного слоя, резко, кризис­ ным образом, перемещается по потоку и устанавливается при угле <р=120 140° в области большего dp/dx>0. Обтекание шара улуч­ шается, кормовая поверхность шара, на которую действует пони­ женное давление, уменьшается и Сх уменьшается за счет уменьше­ ния сопротивления давления (сопротивление трения увеличивает­

ся за счет появления турбулентного пограничного

слоя). При

Re<x,>ReooKp снова наступает зона автомодельности

по числу Re*,

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

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

1

В области 2,5- 103<Reoo<2- ю» можно снизить Сх с 0,47 до 0,1,

если искусственно турбулизовать ламинарный

пограничный слой,

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

И з м е р е н и е и н т е н с и в н о с т и т у р б у л е н т н о с т и в п о т о к а х м е т о д о м ша р а . Положение точки Т перехода лами­ нарного пограничного слоя в турбулентный зависит от степени тур-

булентности e«> набегающего потока и шероховатости поверхности шара. Зависимость ReooKp=/(eoo) при гидравлически гладком шаре представлена на рис. 18.7. Эта строго выполняющаяся зависимость используется для определения степени турбулентности потока по измерению Сх шара. При этом, за ReTOKP принимают значение числа Re», при котором C*Kp=0,3.

Плохообтекаемые тела широко используются в качестве пара­ шютов, тормозных щитков самолетов, стабилизаторов пламени в форсажных камерах сгорания ТРД и в различного рода турбули-

заторах потоков.

В л и я н и е с ж и м а е м о с т и г а з а ил и ч и с л а М» н а в е л и ч и н у Су и Сх. В о л н о в о е с о п р о т и в л е н и е . Рас­ смотрим течение идеального газа при р«> = const между верхней по­ верхностью профиля и отвердевшей поверхностью тока (рис. 18.8). Увеличение числа М» в этом случае соответствует увеличению пол-

ного давления невозмущенного потока

-

(\_\

к~ 1м 2

у

 

— /'«■I 1~г

2

При увеличении Моо от 0 до 0,3 сжимаемость газа практически не проявляется, скорость в любых сечениях струйки, например /_и /,

изменяется в одинаковой степени, коэффициент давления

ро

Р-

р™— ^°°

—( 1

^°°

I т е остается постоянным и Сх и

QX

9 X

v

si

)'

 

2

2

 

 

 

 

Су2 сохраняют неизменное значение (рис. 18.9).

При 0,3<Моо<МКр сжимаемость жидкости проявляется в том, что максимальная скорость течения, например в минимальном се-

чении 2 увеличивается в большее число раз, чем увеличивается ^оо. Это объясняется уменьшением плотности при ускорении газа. Такое увеличение скорости вызывает соответственное уменьшение давления Рг/Р»* и увеличение Су.

Теоретические расчеты и опыты показывают,

что Су профиля

2 области 0,3<Моо<Моонр следует рассчитывать по формуле*

Cy= C yj V 1 -M l,

(18.11)

где Сун — коэффициент подъемной силы данного профиля при за­ данном угле атаки в потоке несжимаемой жидкости (Мв><0,3).

При кризисе Моо=Мкр<1 скорость в минимальном сечении 2 достигает скорости звука, а давление — критического минимально-

го значения ——= ^ + ^к~1 = 31 (1)- При дальнейшем увеличении

Мсо>М кр отношение ргкр/Роо* не изменяется, поэтому

Р«л (1)

увеличивается и Су уменьшается.

 

Коэффициент лобового сопротивления Сх профиля

в области

0,3< Моо< МЬ1, сохраняет неизменное значение вследствие сохра­ нения неизменной картины течения и картины распределения дав­ ления по профилю (при обтекании профиля идеальным газом Сх=0).

Область обтекания профиля MKP<;MDO<;l качественно отлич­ на от области М»^Мкр. При М*,>Мкр, за минимальным сечением 2 , поток, становясь сверхзвуковым, продолжает ускоряться в рас­ ширяющемся канале, а давление — падать. Сверхзвуковая область замыкается скачком уплотнения. Связанное с этим понижение ста­ тического давления на кормовую часть профиля (см. п. 15.6) вызы­ вает резкое увеличение лобового сопротивления давления (см. рис. 18.8), которое в связи с причиной возникновения называется волн о­ вы м сопрот ивлением.

Увеличение Сх продолжается и при Мж> 1 за счет образования перед профилем головной ударной волны.

При дальнейшем увеличении числа М®, когда угол «о полуютана сверхзвукового профиля становится меньше ©щах (см. рис. 12.7), отошедшая ударная волна превращается в присоединенные косые скачки уплотнения н Сх начинает снижаться вместе с уменьшени­ ем угла косых скачков и волновых потерь.

При увеличении М »> 3 ... 5 углы головных косых скачков уплот­ нения изменяются незначительно (см. рис. 1217) и Сх перестает за­ висеть от числа М®.

Область М,„<Мот<1 называется област ью волн ового кри зи са сопрот ивления или звуковы м барьером , так как в ней происходит пяти-, десятикратное увеличение Сж при незначительном увеличе­

нии Мот-

* Число М„.=М«р><1 называется «ел® М жтвщщештт шитика яро кю-

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

С т р е л о в и д н о е к р ы л о — это крыло, передняя кромка ко­ торого образует с поперечной осью самолета угол стреловидности X (см. рис. 18.8, б).

Поток, обтекающий стреловидное крыло, можно представить как сумму двух потоков: 1) тангенциального №* = №00 sin %текуще­ го вдоль крыла и не влияющего на распределение давления по профилю (если пренебречь влиянием вязкости); 2) нормального к крылу {Wn = Woo cos х, определяющего распределение давления по профилю. Следовательно, эффективное число Мп, определяющее

 

характер обтекания профиля

 

Ми = Wn/doo= Moo cos х

тем

 

меньше числа Моо, чем больше

 

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

 

 

 

Результаты опытов показы­

 

вают,

что

увеличение угла

 

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

приводит

к

 

увеличению

Мкр, существенно­

 

му

снижению максимального

Рис. 18.10. К обтеканию шара

значения Сх и смещению его

в

область

более

высоких

чи­

 

сел М о о .

Поэтому все современ­

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

 

В действительности, тангенциальная составляющая потока вза­

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

Это

уменьшает положительный эффект стреловидности.

Для устране­

ния этого вредного влияния вязкости на поверхности

крыла уста­

навливаются ребра, препятствующие

перетеканию

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

слоя.

крыла приводит

к увеличению

Уменьшение удлинения Т = 1 / Ь

Мкр и к снижению Сх (см. рис. 18.8, б) за счет уменьшения разре­ жения около верхней поверхности крыла. Последнее объясняется усилением эффекта перетекания воздуха (см. ниже) и использова­ нием относительно более тонких профилей.

Звуковой барьер в авиации был преодолен за счет двух техни­ ческих достижений:

1) замены поршневых авиационных двигателей реактивными, что обеспечило получение потребных тяг при малом весе силовой установки;

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

В л и я н и е с ж и м а е м о с т и на о б т е к а н и е п л о х о о б ­

т е к а е м ы х

тел. В качестве

примера рассмотрим обтекание

шара. Для

шара М1ф«0,6 (рис.

18.10). При М«><Мкр сжима­

емость газа проявляется в увеличении абсолютных значений dp/dx^O. Это приводит к стабилизации ламинарного погранично­ го слоя и затягиванию кризиса сопротивления на большие числа ReKp, например, ReHp^4 ,5 7 -10^ при Моо = 0,5 вместо ReKP = 3,2* 105 при Моо= 0 (см. рис. 5.2).

При Моо>М1ф на поверхности шара образуется зона сверхзву­ ковой скорости, замыкающаяся скачком уплотнения, который, не­ зависимо от режима течения в пограничном слое, вызывает его от­ рыв. В связи с этим, при Мкр<Моо<1 кризис сопротивления шара со снижением Сх (см. рис. 5.2) совсем не возникает (Сх не зависит от Reoo).

При увеличении числа М*, в области М1ф< М 00< 1 -наблюдается резкое увеличение сопротивления шара, совпадающее с появлени­ ем ударных волн. Точка отрыва пограничного слоя смещается по направлению к передней критической точке (рис. 18.10) и обра­ зуется турбулентный след большего диаметра с пониженным дав­ лением на кормовую поверхность шара. При Моо>1 перед шаром устанавливается отошедшая криволинейная ударная волна и Сх продолжает повышаться и достигает величины Сд;т ах~1,05 при

Моо^ 1,7.

Как показывают эксперименты, точка отрыва пограничного слоя смещается на корму на угол 0 —110°, зона пониженного дав­ ления сокращается и в связи с этим наблюдается незначительное снижение Сх.

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

18.3. КРЫЛО КОНЕЧНОГО РАЗМАХА. ИНДУКТИВНОЕ СОПРОТИВЛЕНИЕ

Пусть крыло конечного размаха / установлено в потоке идеаль­ ной жидкости под положительным углом атаки а и имеет циркуля­ цию скорости Г, направленную по часовой стрелке и положитель­ ную подъемную силу (рис. 18.11, а). В этом случае давление жид­ кости на нижней поверхности крыла больше, чем на верхней Рпаж> > рв• Концевые эффекты крыла конечного размаха состоят в том, что возникает самопроизвольный поперечный ток жидкости из об­ ласти большего давления в область меньшего: на нижней поверх­ ности крыла — от оси симметрии к торцам; затем вокруг торцов; на верхней поверхности — от торцов к оси симметрии. Взаимодей­ ствие этого тока с невозмущенным потоком приводит к образова­ нию около торцов вихревых шнуров и вихревой пелены за задней кромкой. Вихревые шнуры вызывают отклонение невозмущенного потока вниз, уменьшая действительный угол атаки по сравнению с геометрическим, что и является причиной появления индуктивного сопротивления.

При отрицательном угле атаки изменяются на обратные нап­ равления циркуляции скорости, подъемной силы (Рниж<Рв), попе­ речного тока и направление отклонения потока. Индуктивное соп­ ротивление сохраняется. Только при установке крыла на угол ну­ левой подъемной силы (Г = 0; Ry*=0, рнпж=Рв) поперечные токи прекращаются, пропадают вихри и индуктивное сопротивление.

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

крыла.

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

Рис. 18.11. Крыло конечного размаха:

а—простейшая вихревая схема с П-образным вихрем; б—скос потока за крылом конечного размаха

не может окончиться в жидкости (см. п. 3.5), и за торцами крыла подхватывается потоком и вытягивается вдоль линии тока в беско­ нечность. Эти вихри называются свободными или вихревыми усами и имеют циркуляцию скорости Г такую же, как у присоединенного вихря. В первом приближении принимают, что крыло конечного размаха действует на поток так же, как П-образный вихрь, состо­ ящий из присоединенного вихря и двух свободных вихрей *.

Свободные вихри индуцируют за крылом вертикальное поле ско­ ростей, направленное вниз, вызывающее отклонение всего потока на угол Да (рис. 18.11, б).

Средняя скорость скоса потока определяется величинами цир­

куляции скорости и относительного размаха

крыла

(см. Юрьев

Б. Н. Экспериментальная аэродинамика,

ч. II.

Оборонгиз,

М. 1938.):

 

 

 

 

 

(18.

12)

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

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

где l = l/b — относительный размах или удлинение прямоугольного крыла.

Угол скоса потока

tg Да = v J W H= C y/n'l

весьма мал, так как оср<с№п. Поэтому tg A a^sin Д а«Д а и cos Д а » 1.

Таким образом приходим к выводу, что крыло конечного раз­ маха обтекается потоком со скоростью W=WH+ vcp, равной при­ мерно Wn и направленной под истинным аэродинамическим углом атаки a, = a—Да, где а — геометрический угол атаки.

В соответствии с теоремой Жуковского на крыло действует сила

о„Г2

/? = в й ’^ д а = с , у ! 1ь,

направленная перпендикулярно вектору скорости W, т. е. отклонен­ ная от направления, перпендикулярного вектору скорости набегаю­

щего потока WH.

Проекция силы R на ось у дает подъемную силу крыла конечно­ го размаха

Ry = Rcos Дa ^ R .

Проекция силы R на направление невозмущенного потока дает силу индуктивного сопротивления

R xi — R sin Aa — R b.a=Q nW Hrik a .

 

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

конечного

размаха определим с учетом, что Г= CvbWHf2

 

Cxi— Cy^ a = C yjnL2

(18.13)

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

Коэффициент индуктивного сопротивления пропорционален квадрату коэффициента подъемной силы и обратно пропорциона­ лен удлинению крыла. Следовательно индуктивное сопротивление отсутствует либо в случае нулевой подъемной силы (Су = 0), либо при бесконечном удлинении крыла (Тж оо).

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

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

ния увеличивают удлинение Z= 20 ... 30 вместо 1= 6

7 для пасса­

жирских самолетов.

 

Задача 18.7. Придумайте устройство, уменьшающее Cxi

заданного крыла

при заданном угле атаки.