Скачки уплотнения

ГЛАВА ЧЕТВЕРТАЯ СКАЧКИ УПЛОТНЕНИЯ

4-1. ОБРАЗОВАНИЕ СКАЧКОВ УПЛОТНЕНИЯ

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

Рис. 4-1. Сверхзвуковое течение газа в область повышенного давления.

С этой целью рассмотрим движение сверхзвукового потока вдоль плоской стенки АВ, вытекающего в среду с повышенным давлением (рис. 4-1). Слева от точки В скорость будет Cj, давление р_х и температура 7\. Правее точки В (за линией ВС) поддерживается давление р₂, более высокое, чем р_t. Если разность давлений /?₂—р_х мала, то в точке В возникает слабая волна сжатия В Ki

128

6П2.2 Дейч Михаил Ефимович.

Д 27 Техническая газодинамика. Изд. 2-е, переработ. М.—Л.

Госэнергоиздат, 1961

с черт. и илл.    6П2.2

design pashaok

Редактор Б. #. Шумяцкий    Техн.    редактор А. М. Фридкин

ГЛАВА ЧЕТВЕРТАЯ СКАЧКИ УПЛОТНЕНИЯ

4-1. ОБРАЗОВАНИЕ СКАЧКОВ УПЛОТНЕНИЯ

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

Рис. 4-1. Сверхзвуковое течение газа в область повышенного давления.

С этой целью рассмотрим движение сверхзвукового потока вдоль плоской стенки АВ, вытекающего в среду с повышенным давлением (рис. 4-1). Слева от точки В скорость будет с_г, давление и температура 7\. Правее точки В (за линией ВС) поддерживается давление р_г, более высокое, чем р_г. Если разность давлений р_г—р_г мала, то в точке В возникает слабая волна сжатия ВК_г-

Ёсли изменение давления в точке В станет конечным, то, как показывает эксперимент, волна переместится в положение ВК и будет обладать не бесконечно малой, а конечной интенсивностью. По мере увеличения давления р_г линия ВК будет поворачиваться относительно точки В влево (ВК', ВК" и т. д.). При перех'оде через волну ВК газ сжимается и поток отклоняется на некоторый угол 8 вверх от направления невозмущенного потока АВ. С ростом

Рис. 4-2. Обтекание сверхзвуковым потоком внутреннего угла.

р₂ сжатие газа в волне ВК и угол отклонения 8 увеличиваются.

Волна ВК называется плоским косым скачком уплотнения или плоской ударной волной. При переходе через такую ударную волну поток испытывает скачкообразные изменения давления, скорости и других параметров. Положение скачка определяется углом (3 между плоскостью скачка ВК и первоначальным направлением потока АВ (рис. 4-1).

Образование косых скачков уплотнения можно проследить также на простейшем примере обтекания стенки ABC, повернутой в точке В на некоторый конечный угол 8 навстречу потоку (рис.4-2).

Благодаря такому повороту стенки сечение струйки уменьшается и она суживается. В сверхзвуковом потоке это приведет к повышению давления (р₂ >¦ p_t). Причем повышение давления происходит скачкообразно при переходе через поверхность' ВК, являющуюся поверхностью скачка. Можно показать, что при обтекании рассматриваемой стенки непрерывный переход от параметров в области АВК к параметрам в области КВС физически невозможен.

Действительно, границей возмущения для области АВК_Х должна быть звуковая волна ВК_и угол наклона которой

к вектору скорости будет а == arc sin — . Вторая гра-

ница возмущения ВК₂ имеет угол наклона а ₂ = arc sin —.

Так как с_г<С_с₁ и    то    a_m2>a_ml.    Характеристика

ВК* оказывается в невозмущенной области АВКj. и линии тока должны были бы иметь форму, показанную пунктиром, что физически совершенно нереально.

Можно предположить, что косой скачок занимает среднее положение между волнами ВК_г и BK_Z\ тогда угол косого скачка р связан простым приближенным соотношением с углами a_mV а_т2 и 6:

Р ~ Т~ К>1 + ^ат2 + ⁵)'

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

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

Проанализируем условия образования таких движущихся ударных волн. Пусть в трубе постоянного сечения находится поршень (рис. 4-3). Толчок поршня вызывает слева слабую волну разрежения т' — п', а справа слабую волну уплотнения т — я. Продолжая увеличивать ско-рость'поршня толчками, мы создадим ряд слабых волн возмущений (/и,— n_v m_l — «j и т. д), перемещающихся в потоке газа в противоположных направлениях от поршня, каждая со своей скоростью, соответствующей скорости звука в данной области.

Нетрудно видеть, что справа каждый толчок повышает давление газа на малую величину, а слева — снижает его. Следовательно, в области III давление и температура будут выше, чем в областях II и /, 8 следовательно, и скорость звука «ni > ^a_[I^^>ai- Наоборот, в областях II', III' скорости звука будут меньше, чем в области Г (ящ < <[.вл<а₁). Следовательно, справа от поршня слабые волны уплотне.

Ння нагоняют друг друга, слева волны разрежения отстают друг от друга.

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

Рис. 4-3. Распространение слабых возмущений ч в трубе.

4-2. УРАВНЕНИЯ КОСОГО СКАЧКА

Как и ранее, будем рассматривать установившееся течение газа без теплообмена с окружающей средой и без трения. Предположим, что в некоторой точке сверхзвукового штока возник косой скачок уплотнения (рис. 4-4). Параметры газа до скачка обозначены индексом 1, а за скачком — соответственно индексом 2.

Рассмотрим движение газа по линии тока- ABC, пересекающей плоскость косого скачка в точке В. Как указывалось, при переходе через косой скачок линия тока деформируется, отклоняясь на некоторый угол 5. Скорость до и после косого скачка можно представить составляющими, нормальными к плоскости скачка {с_пХ и с_п2) и касательными к ней (с_п и c_t2) и, таким образом, построить треугольники скоростей до и после скачка.

Очевидно, что

nl

= с²

^Сп2 \

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

Рис. 4-4. К выводу основных уравнений косого скачка.

Закон сохранения массы — уравнение неразрывное и—для двух сечений трубки тока до и после скачка может быть записано в следующем виде:

(4-1)

Pl^C«l = P^C„2-

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

Pi — Pt = ?iCnl(^Cn2 — ^Cnl)>

или

Pi + 9i^c\x = Рг + Ра °²_пг В проекции на плоскость скачка получим:

(4-2)

Pl^Cnl (^Ct2 ~~ ^Ct]) ⁼ Q’

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

^ct\~^cn~^cr    (4^_3)

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

Закон сохранения энергии — уравнение Бернулли— может быть использовано в любой известной нам форме. Рассматриваемое течение происходит без теплообмена с окружающей средой и, следовательно, полная энергия потока сохраняется неизменной:

^ci I k р, _ ^с\ I    к р_г ^а* k + I

(4-4)

Найдем связь между скоростями. до и после скачка. Преобразуем уравнение (4-2); с учетом формулы (4-1) получим:

Р± + _с² ^l(Pz+_c²\.    (4-2а)

Pi ^{Т п1 с}я2    ^п2)

Из уравнения энергии (4-4) выразим отношение давления к плотности:

2

-4)    <⁴’^4а)

2/? ^- \ к — 1

И .

h ь

-*ir («г iir-*'--'!)-    ⁽⁴-^4б>

Подставив (4-4а) и (4-46) в уравнение (4-2а), получим: «²4-г²    )    ^Сд1    — ^к~~^{{ c}«' J

Формула (4-5) устанавливает связь между нормальными составляющими скоростей при переходе через косой скачок

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

Заменим в формуле (4-5) al по уравнению

тогда

к — 1 / ^cn! 1    2 \    2    гл а\

°n\ ^Cn2~k ₊ 1    ^2    ^а\ ’

k — 1 / ^cn2 I 2    \    2    /Д й \

^С,г\^Сп2— Г+\[    ?    l)^a2>    (^4_ба)

или, выражая скорость звука через давление и плотность а^% = k , получим:

(4-7)

Воспользуемся уравнением неразрывности (4-1); возведя в квадрат обе части этого уравнения:

г 2    2 2

Pi ^сп\    Рг ^сп2

и заменив

получим:

,2 ,2

/г!    „    ^ьпЪ

С

АРх-^=ЛР.-^--    (4-8)

°1 ^а2

Исключая последовательно из уравнения (4-7) с помощью уравнения (4-8) и р₂ или jt?j и р_а, приходим к выражениям

Из уравнений (4-6) и (4-ба) можно получить:

Л

^Ln\

к ¦

2k

п\

Из треугольников скоростей на скачке следуют очевидные соотношения:

¦¦ с_х sm р;

* п\

(4-12)

^сп2=^с. ^sinT; с == С_х COS Р =

где

Т —Р~⁵-

Тогда из уравнения (4-11) совместно с формулами (4-12)

можно получить:

k — 1 о

(4-1 la)

¦ 2~Mi sm²ft+l

k — \ sin² p — 1

или

^l~k+l^ ^cos^ sin² p

(4-116)

Уравнения (4-9) и (4-10) совместно с уравнениями (4-11) и (4-12) после несложных преобразований дают искомые зависимости:

Се k — ¦" 1 / 2Й | >2    • п л    4    \

-'- м sin² 8 — 1 —

k+ 1

1    I²    cin!    ft

Pi 2    1    k—    1,    ,

-j- • —^-4-1    1    —¦ 7ПГТ ^-l ^cos2 P

^k — ¹ Л4² sin² p ^T    fe    +    1    i    ^f

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

7*2 Р2Р1 JP1P2 '

После подстановки ^ ш получаем:

uotvn — И

Pi

Г₂ k — \^(2k

T-=\rzif    -2---ИУ

Г, -|-1у уб    / у^—*.M[Sin²p J

(4-15)

«Формулы (4-5), (4-11) — (4-15) выражают изменения параметров газа при переходе через косой скачок уплотнения в зависимости от показателя k, скорости потока до скачка М_г и угла косого скачка р. Эти формулы вместе с тем вскрывают основные физические свойства косого скачка. Так, из формулы (4-13) можно заключить, что нормальная составляющая скорости до скачка больше критической скорости. Действительно, так как —>1 и

Р1

_с2

.    «    о ^ьп 1

М, sm² р = -у-, заключаем, что

^а\

2k ( 4i k-

>1

*+1 V а? 2й и, следовательно,

п\

>1.

, Нормальная составляющая скорости за скачком меньше критической скорости:

^С_П2 ^ ^а*'

Из формулы (4-13) следует, кроме того, что угол косого скачка больше угла характеристики а_т1.

Рг Рг __ 1

Р~ Pi Т₁

В этом случае косой скачок вырождается в слабую (звуковую) волну уплотнения (слабый скачок) и угол отклонения потока стремится к нулю.

Связь между углами р и 5 устанавливается по уравнению (4-6), которое можно преобразовать к виду:

с„₂ _А — 1 f ,    2    а?    ^

k+l^-Tk-l 4, /

Но на основании формул (4-12)

^сп2 _ ^Сп1    ^ct  _ tgf _ tg (Р — с

(4-16)

°п\    ^ct    ^сп\

Следовательно,

tg(p -S) __k-_\

‘gP Й + Ч k—\M\ sin²1

Вспоминая, что

получим:

sin² p —:

(4-17)

-i- - sin² f +

tg

На рис. 4-5 представлены графики 8({3) при различных значениях для 6 = 1,3. Отметим, что с увеличением скорости невозмущенного течения увеличивается максимальный угол отклонения потока 8 . Следует подчеркнуть, что в соответствии с двойственным решением уравнения (4-17) одному значению угла отклонения потока S соответствуют два различных значения р. Опыт показывает, что плоскому косому скачку отвечает только меньшее значение р.

Щ

ВьШе в предположении, что скаЧоК занимает среднее положение между характеристиками невозмущенного и возмущенного потоков, была дана формула

Р ~ 4"(^а_т1+^ат₂ + ⁵)-    (4-17а)

Сравнение этой формулы с точным выражением (4-17) для нескольких значений приведено также на рис. 4-5. Кривые показывают, что с увеличением скорости перед скачком расхождение между результатами расчета по

Рис. 4-5. Зависимость между углом отклонения потока и углом скачка при различных скоростях невозмущен-ного потока (для k = 1,3).

--точная формула (4-17);---приближенная формула

(4 17а).

точной и приближенной формулам увеличивается. Величина погрешности зависит также от 8.

Из уравнения (4-17) следует, что 8 = 0 при Р = а_т1 и

при р = -^-. Таким образом, кривая 8 = 8(fi) имеет "Макси-

кум, положение которого определяется обычным способом. Продифференцировав уравнение (4-17) и приравняв производную нулю, после преобразования получим:

ь±1 _м\

sin²p„

4 "^Jl

kMi

(4-18)

+

16 ^ivJl

у^r (*+!)(

где — угол косого скачка, соответствующий максимальному углу отклонения потока 6

Отсюда следует, что при Mi = 1 угол Р_т = ~, а при

М₁ = оо угол P_m = arcsin {^f j • Для промежуточных значений угол Р_т с ростом М_г вначале уменьшается, а затем несколько возрастает.

Уравнение (4-11а) позволяет проследить изменение скорости потока за косым скачком М₂ в зависимости от М_г и р. С увеличением р (при постоянном М_г) М₃ уменьшается: перепад скоростей в скачке увеличивается.

При некотором значении = скорость за скачком становится звуковой (М₂=1). При дальнейшем увеличении (3 течение за скачком будет дозвуковым.

Величину р* можно определить по уравнению (4-11 а), подставляя /И₂=1. Тогда после преобразований получим:

sin² Р,

kMf

k +1 16

< + k

')

М

(4-19)

¦М,

16

+ J/^ (^+^J)

Заметим, что при М_± = 1 угол    = [З_т =    при М_х-

= оо угол p* = p_m = arcsin(}/^r-).

максимальный

Последнему значению р* отвечает угол 8 , определяемый из (4-17):

‘<7Р=г)'

:агс

Для значений А^<оо угол |i_m>[5* и, следовательно, М_2т<1. Это означает, что при максимальном угле отклонения потока скорость за скачком будет дозвуковой. Так как, однако, для всех М_г углы р_т и fL весьма близки, то в первом приближении можно считать, что максимальный угол отклонения для каждого значения скорости невозмущенного потока достигается при звуковой скорости за скачком (М_г^ 1).

Мы установили, что параметры течения за косым скачком зависят от угла косого скачка р. При увеличении [В давление, температура и плотность газа за скачком увеличиваются (параметры потока до скачка предполагаются неизменными), а безразмерная скорость уменьшается. Угол отклонения потока, как было показано, вначале увеличивается (при (3<^_т), а затем уменьшается (при pi>-j5_m).

В частном случае [3 = -^- изменения параметров в скачке

оказываются максимальными, а угол отклонения 8 = 0. Такой скачок расположен нормально к направлению скорости невозмущенного потока и называется прямым скачком.

Прямой скачок является частным случаем косого скачка; основные уравнения прямого скачка получаются из формул

(4-11) — (4-15) после подстановки Р =

Изменения давлений и плотностей в прямом скачке найдем из формул (4-13) и (4-14):

h

Pi

(4-20)

отношение температур — из формулы (4-15):

Безразмерная скорость за прямым скачкок кожет быть получена по формуле (4-11а):

^l + jTTT ^2=-2k - ^

к— 1/^и1 или по формуле (4-116)

КК = 1.    (4-24)

т. е.

с₁с₂=а\    (4-24а)

Мы видим, что произведение скоростей до и после прямого скачка равно квадрату критической скорости. Отсюда прежде всего следует, что скорость газа за прямым скачком всегда меньше критической скорости (с₂<<а*).

Формулы (4-20) — (4-23) показывают, что интенсивность прямого скачка увеличивается с ростом скорости невозмущенного потока M_t (или >Ц). Отношение плотностей при максимальной скорости стремится к конечному пределу

Ит—=-

_Pl к —г

1 к—1

а отношения давлений и температур возрастают безгранично.

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

Формулы косого скачка могут быть преобразованы к виду, удобному для анализа влияния физических свойств газа (показателя k).

с

С этой целью введем безразмерную скорость ? = - и выразим

?макс

показатель изоэнтропического процесса через максимальную скорость ₂ k + 1

макс k — 1 '

^макс + ¹

^{й =} х² “ГГ'

макс

Заменив k в уравнениях (4-13) и (4-14) его значением, получим:

Рг _^макс ~~ ¹ Sin² g _    1

п    т» 2    1    tl2

^    ^Лмакс    —^_1    макс

(4-13а

Ь

X²

р2 __ макс

(4-14а

~р7^_    1 - f,

⁺ Si sin²

Скорость за скачком выражается уравнением [формула (4-116)]

}² {-2 ^к2 «2

-Г2---I2-—COS⁵⁰?+/(?„ Р).    (4-116)

h ?1

где

—2 cos-’ji

/>. Г \    1

f («1> Ю    sin² ji X⁴

макс

Как видно, каждая из приведенных формул содержит два сомножителя, один из которых зависит только от ^ и р и не зависит от k, а второй является функцией только k.

Такая структура формул скачка позволяет приближенно оценить влияние изменения физических свойств газа и производить расчет параметров косого скачка при различных постоянных значениях k*.

Для определения других параметров скачка можно пользоваться очевидными соотношениями- Отношение температур

Г» 1-5*

Угол скачка

1-5?

tgP=| 1+-^"ЖТГ 1*8»

sm²B ^ьчг

^г /    макс

Рис. 4-6. Влияние показателя k на интенсивность косого скачка.

Для расчета скачков при различных k могут быть построены графики p₂/Pi в зависимости от g,, (I (или 6) я k (рис. 4-6). Влияние k можно оценить при одинаковых р или 3. Сравнение при одинаковых р показывает, что с уменьшением k интенсивность скачка возрастает.

4-3. УДАРНАЯ ПОЛЯРА

Зависимость между параметрами на скачке можно в весьма удобной форме представить графически. С этой целью рассмотрим треугольники скоростей на скачке (рис. 4-7).

Расположим вектор скорости до скачка с_г по оси х (отрезок OD). Отрезки OF и FD представляют собой соответственно касательную c_t и нормальную с_п1 составляющие скорости до скачка. Зная угол отклонения потока 8,

проведем линию вектора скорости за скачком с_г до пересечения с отрезком FD. Точка пересечения (точка Е)определяет величину вектора с_а, а отрезок EF выражает нормальную составляющую скорости за скачком.

Вектор скорости с₂ можно представить двумя другими составляющими: «_а и v_t. Компоненты ы₂ и v₂ являются проекциями с₂ на направление скорости потока перед скачком и на нормаль к этому направлению.

Найдем уравнение кривой, описываемой концом вектора скорости за скачком с₂, при постоянном значении вектора скорости перед скачком с_г и переменных значениях угла поворота потока за скачком 8.

Выражая это уравнение в форме связи между ы_а и v_s, мы получим кривую скорости за скачком в плоскости годографа скорости.

Для получения искомой зависимости используем основное уравнение косого скачка (4-5). Подставив в это уравнение значения с_п1 и c_t из формул (4-12), получим:

с, sin р ^ sin р -    =    а-щ±с] cos² р, (4-25)

так как (рис. 4-5)

= ^a_t l^V2 + (^Cl — “*) Ч - k+1 ^Cl^Vr

Окончательно

a²

*

o! = (Ci —--- •    (⁴‘²⁶)

К ²

77⁺F+l^Ci~~ ^u*

Кривая, отвечающая уравнению (4-26), представленная на рис. 4-8, называется ударной полярой. Кривая при надлежит к классу гипоциссоид. Ударная поляра может быть широко использована для расчета косых скачков графоаналитическим методом и для выяснения некоторых особенностей таких скачков.

Обратимся прежде всего к предельным значениям и₂, даваемым уравнением (4-26).

Легко видеть, что у₂ = 0 при

И₂-С \

и

г

Первый случай (и₂ = с_t) отвечает бесскачковому процессу; косой скачок уплотнения переходит в волну слабого возмущения (характеристику). Касательные к гипоциссоиде в точке D расположены под углом a_ml = arcsin

Рис 4-8. Ударная поляра в плоскости годографа.

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

Второй случай (u_ic₁=a²) характеризует переход косого

скачка в прямой скачок, угол которого = Этому случаю на гипоциссоиде отвечает точка А.

Из уравнения (4-26) следует, что и₂ может обратиться в бесконечность при

_    ²    I

—6 + l^Cl + с, •

Очевидно, что ветви гипоциссоиды асимптотически приближаются к прямой, проведенной параллельно оси v на

— °² 2

расстоянии ОН = ~-г]^ci ^от начала координат. Эти

ветви не имеют физического смысла, так как они дают значения скорости за скачком (точка Е₃ на рис. 4-8), большие, чем до скачка. Таким условиям отвечал бы скачок разрежения, но скачков разрежения существовать Me Mokef. Отбрасывая внешние ветви гипоциссоиды как физически нереальные (см. ниже), мы замечаем, что ударная поляра в пределах между крайними точками А и D дает два значения для вектора скорости за скачком.

Обычно плоские скачки реализуются при значениях вектора скорости потока за скачком, отвечающим точкам Е₂ (рис. 4-9,а). Второе значение скорости с₂, соответствующее

Рис 4-9 Отрыв и искривление скачка уплотнения

точкам Еj, в плоском скачке может быть реализовано лишь при специальных условиях.

Рассмотрим сверхзвуковое течение газа вдоль стенки LBC (рис. 4-9,а), постепенно увеличивая угол отклонения потока 8 (угол поворота стенки в точке В). При малых значениях 5, близких к нулю, возмущение потока невелико и скорость за скачком с₂ близка к скорости до скачка с_г. По мере увеличения § точка Е_г (рис. 4-9,6) перемещается вдоль ударной поляры от D к г, где точка г дает скорость за скачком М₂ = 1. Дальнейшее весьма небольшое увеличение 8 приводит поток за скачком к состоянию, определяемому точкой К. Здесь течение за скачком уже дозвуковое (М₂<^ 1) и § достигает максимального значения 5_т>

На рис. 4-10 представлено обтекание клина сверхзвуковым потоком. Если половинный угол раствора клина 8 меньше Ь_т для данной скорости М_г, то на носике клина возникают два прямолинейных косых скачка: АВ и АВ,, образующих так называемую головную ударную волну клина.

При дальнейшем увеличении угла 8>8_т скачок отходит от передней точки и искривляется (рис. 4-9,б). Это объясняется тем, что скорости распространена возмуЩё-ний становятся больше скорости потока. Действительно, увеличивая угол поворота стенки 8, мы тем самым увеличиваем поджатие потока, т. е. его давление, плотность и температуру. Вместе с тем растет и скорость распространения возмущений, равная скорости звука возмущенного потока a^^^kgRT^. При о]>о_п эта скорость становится больше скорости потока и поэтому возмущения проникают вперед по потоку. Однако при удалении от стенки ВС (рис. 4-9,б) давление, плотность и температура

a)    ^N    б)

Рис. 4-10 Обтекание клина сверхзвуковым потоком.

будут уменьшаться; вместе с тем будет падать и скорость распространения возмущений. На некотором расстоянии от стенки возникнет геометрическое место точек PQ (рис. 4-9,в), в которых скорость распространения возмущений снизилась до скорости потока. Очевидно, за пределы этой поверхности возмущения, вызванные стенкой, не смогут проникнуть, так как они будут сноситься потоком назад. Поверхность PQ отделяет зону невозмущенного потока от зоны возмущенного потока и представляет собой отошедший скачок уплотнения.

Следовательно, если 8 >¦ Ъ_т, то плоский косой скачок сменяется криволинейным скачком (рис. 4-10), который располагается не на носике клина, а на некотором расстоянии перед ним. Это расстояние зависит от скорости невозмущенного потока УМ, и 8. С ростом М₁ скачок приближается к носику тела. С увеличением угла отклонения при 8>8_т екачок удаляется от тела. Обтекание скругленного но-'€ика тела сверхзвуковым потоком всегда будет происходить с образованием криволинейной головной волны, отор-

вавшейся от носика, а расстояние между волной и носиком для центральной линии тока будет зависеть от скорости Мj и от формы носика.

Заметим, что при 8 = §_т поток за скачком дозвуковой и М_г несколько меньше 1 (точка К на рис. 4-9,6). Так как для нейтральной линии тока, разветвляющейся в точке А

(рис. 4-10), р = -^- и 2 = 0, то элемент скачка, пересекающий эту линию, должен быть прямым. Скорость потока за элементом прямого скачка определится точкой А на ударной поляре (рис. 4-9б,). Поток за скачком на этой линии тока всегда дозвуковой.

Все участки скачка, кроме центрального, расположены под различными углами к вектору скорости невозмущенного потока    •

Рассматривая такую искривленную головную волну, ' состоящую из большого числа малых прямолинейных элементов, можно убедиться, что по мере удаления от центральной линии тока уменьшается и угол наклона элементов скачка р . При этом можно воспользоваться ударной полярой для расчета потока за скачком для каждой линии тока в отдельности Участку головной волны KL отвечают точки ударной поляры от Л до г, в которой скорость М₃ = 1. На этом участке поток за криволинейным скачком будет дозвуковым. Следовательно, если головная волна отрывается от тела, то в некоторой области, прилегающей к носику тела, течение будет дозвуковым (эта область на рис. 4-10,6 заштрихована), а линии тока здесь будут иметь разную кривизну. В различных точках за скачком давления будут различными.

По мере удаления от точки К уменьшается наклон элементов скачка и вместе с тем уменьшается интенсивность скачка. В некоторой точке L скорости за скачком становятся звуковыми. Выше этой точки состояние за скачком определяется отрезком ударной поляры от г до D. На бесконечном удалении от тела криволинейный скачок вырождается в слабую волну уплотнения, которой отвечает точка D на ударной поляре.

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

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

В табл. 4-1 приведены значения Ь_т для двух значений &=1,4 и 1,3. Значения 8_т зависят от числа М₁(Х₁) и физических свойств газа (k) и могут быть определены по уравнению (4-17) при условии подстановки в него Р = Р_т из уравнения (4-18).

Таблица 4-1

Значения максимальных углов отклонения потока в плоском косом скачке уплотнения

4-4. ИЗМЕНЕНИЕ ЭНТРОПИИ В СКАЧКЕ

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

Интегрируя это уравнение, получим

Для обратимого (изоэнтропического) процесса Д« = 0 и

Посмотрим, как изменяется энтропия при переходе через скачок уплотнения. Исключив из уравнений (4-13) и (4-14) iWjSin²p, получим:

Рг — k— Г Pi

k -f- 1 р₂ 1 — U    1    ’

(4-28)

Р1    h__k + 1

Pi k — 1

Произведя расчет, легко убедиться, что для скачка

уплотнения, для которого —>>1, всегда

Pi

и, следовательно, согласно (4-27) при переходе через скачок энтропия газа возрастает.

Увеличение энтропии при переходе через скачок объясняется необратимым характером изменения состояния газа в скачке, „ударным* характером процесса. В результате такого процесса часть кинетической энергии газа необратимо переходит в тепло; при отсутствии энергетического обмена с внешней средой внутренняя энергия потока необратимо возрастает. Кривую, характеризующую процесс, протекающий по уравнению (4-28), называют ударной адиабат о^й.

Для скачка разрежения —<<1 и    <    1    получаем

Pi    Р1 \ Рг

As < 0, что невозможно при отсутствии энергетического обмена с внешней средой, так как противоречит второму закону термодинамики. Отсюда следует, что скачки разрежения в энергетически изолированном течении не могут возникнуть, так как это противоречило бы второму началу термодинамики. Однако, как показано в гл. 3, в сверхзвуковом течении волна разрежения с непрерывным падением давления представляет собой вполне устойчивое явление; энтропия при переходе через такую волну сохраняется неизменной.

Изменение энтропии в скачке уплотнения может быть выражено через параметры М_г и В путем замены в уравнении (4-27) — и -г по формулам (4-13) и (4-14); тогда

(4-29)

Уравнение (4-29) позволяет заключить, что изменение энтропии в косом скачке зависит от скорости невозмущенного потока Мj и угла косого скачка р.’ При неизменной скорости М, энтропия при переходе через скачок меняется в соответствии с изменением угла р. Если скачок прямолинейный и, следовательно, вдоль скачка {3 сохраняет постоянное значение, то для всех линий тока, пересекающих скачок, изменение энтропии будет одинаковым. Если же скачок криволинейный, то увеличение энтропии для каждой линии тока будет различным, так как вдоль скачка угол р меняется. Это означает, что за криволинейным скачком поток будет вихревым; за прямолинейным скачком течение остается потенциальным.

Пользуясь уравнением (4-28), можно рассмотреть изменение состояния г^за при переходе через слабый скачок. Полагая р_г = р, Pj = р и считая, что давление и плотность в скачке изменяются на бесконечно малую величину, т. е.

P»=P + dp; p_a = p-fdp;

из (4-28) получим:

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

Р    Р    '

Таким образом, изменение состояния 6 сКаЧке бесконечно малой интенсивности (слабый скачок) является изоэнтропи-ческим.

4-5. ПОТЕРИ В СКАЧКАХ УПЛОТНЕНИЯ. ПОСТРОЕНИЕ ПРОЦЕССА В ДИАГРАММЕ i-s. СКАЧКИ УПЛОТНЕНИЯ В РЕАЛЬНОМ ГАЗЕ

Рассмотрим более подробно энергетические преобразования в скачках. Предполагая энергетически изолированное течение,    заключаем,    что полная    энергия потока при переходе через    скачок не меняется.    Это значит,    что

2    2

^CI I    k Р] _ ^с2    I k р_г

2 ' k — 1 2 ¹ k — 1 р₂ ’ или, пользуясь параметрами полного торможения,

Еп ₌ Р™ .    (4-30)

Pol Ро2

Условие (4-30) можно заменить эквивалентным условием постоянства энтальпии торможения при переходе через скачок:

^г‘о; = ^г'о₂ =    (4-30а)

или при

с.

, = const; 1

-т -т    (⁴-³⁰⁶)

01-¹ 02-¹ 0‘ »

Имея в виду эти условия, рассмотрим процесс скачка в диаграмме is (рис. 4-11). Зная давление торможения до скачка р₀₁ и энтальпию торможения г₀, найдем в диаграмме is точку О j, характеризующую состояние изоэнтропи-чески заторможенного газа до скачка. По известной скорости потока до скачка c_t или давлению р_х находим точку D, которая определяет состояние движущегося газа перед скачком. В скачке статическое давление потока увеличивается до р₂. Если известен угол отклонения потока 8 и, следовательно, р, то Состояние газа за скачком определено (точка Е_г на рис. 4-11), так как по формуле (4-29) можно найти приращение энтропии As⁵¹. Заметим, что линия, соединяющая точки D и Е₂ на рис. 4-11, не характеризует изменения состояния газа в скачке, так как в диаг-

в тепловой диаграмме.

Рис. 4-12. Коэффициенты потерь и отношение давлений торможения на скачке в зависимости от угла р и скорости невозмущенного потока.

рамме is неквазистатические процессы могут быть представлены только начальной и конечной точками процесса.

Если поток за скачком изоэнтропически затормозить, то состояние полного торможения характеризуется точкой 0₂, в которой легко находится значение р₀₂. Если предоставить теперь потоку возможность изоэнтропически расшириться до давления перед скачком, то его состояние будет определено точкой Е'₂ . Скорость газа при этом может быть вычислена по уравнению энергии:

Aci²    .,

~2g~    ^2    ^02

где Я₀з — изоэнтропический перепад энтальпий за скачком. Величину Я₀₂ можно рассматривать как сумму

^02=-^Ок +^0П’

где Я_0к — кинетическая энергия потока за скачком:

тт __^2_ .

^П0к — ₂g ’

я_0п — изменение потенциальной энергии потока в скачке: Очевидно, что

Я₀ ^{2 Я}о 1 >

где Hoi ⁼ ~2^~—изоэнтропический перепад энтальпий до скачка.

Потери кинетической энергии в скачке будут:

Ah = H₀₁ — H_0z = ^ (cf — с\ ),

где Д/г легко определяется по диаграмме is как разность энтальпий i₂ — г\.

Потерю энергии нетрудно связать с основными параметрами скачка. Выразим

к-i г..

так как

(4-32)

Я₀

ft-i р±\ к

Ро2 Рп_Роз

Yol Yo2

Введем понятие коэффициента потерь энергии в скачке-A h , Я₀

Я₀₁    ¹

ft—1

ft

тогда после подстановки Я₀₁ и Я₀₂ получим:

,1

(4-33)

7^ ^_lj

'к— ¹ м\

или с помощью (2-21)

Г k+l J__1 1 Г / Pol',⁵²-¹ _ 11 .

(4-ЗЗа)*

^ L*-1 _Х2 \pj

Отношение ~—е₀ характеризует изменение давления

полного торможения в скачке. Эту величину можно представить в зависимости от параметров скачка М₁ и р.

Так как

/⁾02__ Рг Ро2

Pi ~ Pi Pi’ то- из (4-13) и табл. 2-1 получаем:

^^sin’ ') (' + ^^м* Г - ⁽⁴-^34>

Отсюда с помощью (4-11а) получаем:

6+1    _А_

^?0 — k

Представим в зависимости от тех же параметров скачка коэффициент потерь в скачке ?_е, используя формулы (4-33) и (4-35):

Заменив число M_t через Я₁ по формуле (2-21), получим зависимость в виде:

(4-36)

С* =?(*», р).

На рис. 4-12 представлен график е₀ и в косом скачке в зависимости от угла р и безразмерной скорости невозмущенного потока 4 для k = l,3. Из графика видно, что коэффициент потерь интенсивно возрастает с увеличением угла косого скачка и с ростом скорости Я,, достигая максимальных значений при р = -^- (прямой скачок).

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

Таким образом, в сверхзвуковом потоке идеальной жидкости появляется особый вид сопротивления — волновое сопротивление, зависящее от потерь в скачках, а следовательно, от формы и интенсивности скачков. Как мы видели, форма скачка и его интенсивность зависят от формы тела и скорости обтекания. Учитывая, что при уменьшении угла отклонения 8 (а, следовательно, и (3) потери в скачке уменьшаются, можно заключить, что остроконечные тела в сверхзвуковом потоке должны обладать меньшим сопротивлением, чем тела, имеющие скругленную форму.

Изменение потерь в скачках в зависимости от их интенсивности можно проследить в тепловой диаграмме. Построе-

Мне „ударной поляры* в тёпловой диаграмме удобно выполнить следующим образом¹.

По параметрам до скачка р_х и 7\ находим точку D (рис. 4-13) и при известной скорости с_г точку О_г Задаемся

- о    _п    .1

рядом значении р в пределах от р = —. Для каждого значения ,3 определяем

arcsra-

до

Я,

Ок

Рис. 4-13. Построение в тепловой диаграмме линия изменения состояния газа, соответствующей ударной поляре.

М,

~. Для каждого значения ,3 определяем — « ^Т*

Z    р    j

носим в диаграмме is точки Е'₂ , Е” и т. д. до точки А, которая соответствует прямому скачку. Геометрическое место этих точек дает состояние газа в координатах is, соответствующее ударной поляре. Заметим, что полученная

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

Для каждой точки кривой (например, Е") легко определяются:    кинетиче

ская энергия за скачком

, изменение’ по-

¹ Соответствующая кривая в диаграмме is может быть названа „ударной полярой" лишь условно, т. к. она не является векторной кривой. Приводимое построение справедливо для идеального газа.

тенциальной энергии в скачке #_оп и потери кинетической энергии Д/г. Вместе с тем здесь можно определить и все параметры скачка: р₂, Т⁵³, р₂, р₀₂, As и тем самым существенно дополнить ударную поляру, построенную в плоскости годографа. В плоскости годографа ударная поляра дает кинематическую картину течения (поле скоростей); в тепловой диаграмме мы получаем характеристики энергетических преобразований на скачке и изменения параметров потока. На линии DA легко найти также точку г, соответствующую М₂ = Я₂ = 1. Для этой цели подсчитываем критический перепад энтальпий:

гг    ^Aal    k—\.

Н —=•, , , i_a.    2 g    k + 1 ⁰

и откладываем его от линии i₀= const.

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

Для расчета скачков оказывается весьма удобным пользоваться специальными диаграммами. Такие диаграммы позволяют легко определить характеристики скачка по двум заданным параметрам. В приложении даны диаграммы косых скачков для k = l,3 и & = 1,4*.

Способ пользования диаграммой косых скачков поясняетси иа рис. 4-14. В правом верхнем квадранте диаграммы нанесены графики Рз —

5 = 8 (Р) н -у- = р (Р\ отвечающие различным, но постоянным значениям скорости перед скачком Xj. На каждой кривой надписано значение X,; в скобках указывается та величина, которую изображает данная кривая. В левом верхнем квадранте представлены графики В = Pi —

= д(Ха) и =/>(Х₂) для различных, но постоянных значений Xj. Pi

В левом- иижнем квадранте дана зависимость отношения температур

Т —

на скачке от скорости за скачком -=г- = Т (Х₂). В правом нижнем квад-

¹1

раите нанесены кривые коэффициентов потерь энергии в скачке К_с — = ?_СФ и коэффициентов восстановления давлений s₀ = e₀(ji). Таким образом, в качестве параметра для всех кривых диаграммы выбрана скорость потока до скачка X_v

Поясним на примере способ пользования диаграммой. Допустим, что ная известны угол отклонения линии тока S и скорость потока до скачка X_t. На кривой 8 = 8 ((1), соответствующей заданному значению Лц находим точку А. Проектируя эту точку на горизонтальную

Рг —

ось, найдем в точке А_г угол косого скачка Ца кривой — = р(?),

отвечающей тому же значению А,, получим точку А₂, которая определяет отношение плотностей p₂/pi- Зная А, и д, на кривой д = 5, (Х₂)

Рис. 4-14. Способ пользования диаграммой косых скачков.

в левом квадранте находим точку B_v которая определяет безразмерную скорость за скачком Х₂. Перейдя при том же значении Х₂ на кривую

~ — ~р (А₂), получим в точке В₂ отношение давлений на скачке — . Pi    Pi

Т₂ —

В точке С на кривой у- = Т (Х₂) определяем отношение темпера-Т

тур -JT- . Проектируя точку А_г на линии = С (Р) и е₀ = е₀(Р) в точках D, и П_г, найдем значения коэффициента потерь энергии и коэффициента восстановления давления. Рассматриваемая диаграмма позволяет производить расчет скачков по любым двум параметрам. Так,

Рг

например, заданными могут быть: X, и    — -и д; 8 и р и т. д.

Тепловая диаграмма удобна для расчета скачкой упЛотнёнйя в реальном газе и, в частности, во влажном паре и диссоциирующем воздухе Параметры потока перед скачком и угол отклонения на скачке обычно известны. Задаваясь рядом значений угла скачка Р₍, нетрудно найти соответствующие значения нормальных составляющих скорости.

с j cos

^с«1е ^=С1 ^SIn Рг' ^cn4i = ~7t_g (J_t _ 8) •

По основным уравнениям (4-1)—(4-4) определяются параметры за скачком 1₂₍, p_2l и удельный объем v_2l, отвечающие текущему значению Искомое решение можно наити в диаграмме i-s в точке пересечения двух кривых, одна из которых построена по параметрам i_2(i p_2t, а вторая—по i₂₍, v_9l (рис 4-13,а) С диаграммы снимаются значения (₂, р_г, и_г, действительный угол наклона скачка определяется по формуле

t_g р = Ini + ,/ (jlzlLY-7,

^    2 tg 5 - у \2Xgb)

-

где v = — .

°2

Эта же задача решается с помощью вспомогательного графика, на котором наносятся кривые v_2t (Р₍) и v_2[ (f₍) (ррс 4-13,б). Для построения второй кривой необходимо по диаграмме i-s перейти от значений р% и г₂ к значениям v₂ Решение получается в точке пересечения этих кривых, где определяются и ( и соответственно все остальные параметры за скачком В частности, скорость за скачком

сj sin р * ^Сг=1г sin(P-fi)-

4-6. ПЕРЕСЕЧЕНИЕ СКАЧКОВ

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

Два последовательных поворота стенки ABCD (рис. 4-15,а) на угол 8 приводят к образованию двух косых скачков: ВК и СК. Угол второго скачка p_a>pi, так как после первого скачка поток имеет скорость Х₂<^Х_1ШВ результате скачки пересекаются в точке К• За точкой пересечения оба скачка сливаются в один скачок KF. Линия тока, пересекающая систему двух скачков, деформируется, поворачиваясь в точках b и d на угол 8; при пересечении скачков скорости потока падают, а давления растут скачкообразно.

Рассматривая линию тока КН и предполагая, что проникновения частиц газа из области 4 в области 3' и, наоборот, из области 3' в область 4 не происходит, можно заключить, что в указанных областях давления и направления скоростей должны быть одинаковыми (поперечный градиент давления отсутствует). Но если предположить направление линии тока за скачком KF таким же, как и за вторым скачком СК, т. е. что суммарное отклонение линии тока составляет 28, то давления в областях 3 и 3' будут различными, так как линия тока bd прошла че-

расположенйых косых скачков

рез два скачка, а линия тока КН — только через один скачок; следовательно, потери в первом случае будут меньше, чем во втором (при сжатии с бесконечным числом скачков бесконечно малой интенсивности процесс будет изоэнтропическим—без потерь), причем р₃ может быть меньше или больше р₃. Отсюда можно заключить, что

области течения 3 и 3' разделены слабой волной разрежения или слабым скачком уплотнения KL, при пересечении которого поток приобретает давление p_i = p₃.    Ин

тенсивность волны KL в этом предположении легко определяется. Действительно, зная угол 8 и находим давление р_А за скачком KF. Давление р₃ нам известно в результате расчета скачков ВК и СК. Отношение — дает

Рз

интенсивность отраженной волны KL.

В    общем случае    углы отклонения потока в    точках

В    я    С могут быть    неодинаковыми. При этом в    зависи

мости от соотношения углов 8j и §₂ (8_t — отклонение в первом скачке ВК и 8, — отклонение во втором скачке СК) и суммарного угла отклонения ^ -j- Ь₂ меняются интенсивность отраженной волны KL, а также суммарные потери в рассматриваемой системе скачков. Расчеты показывают, что интенсивность отраженной волны KL, как правило, невелика и поэтому поворот потока на этой волне пренебрежимо мал, что и делает обоснованным предположение о повороте потока на волне KF на угол S_t —{— S_a.

В    зависимости от скорости невозмущенного    потока

и    суммарного угла    отклонения Sj —{— 8_а меняется    и знак

волны KL-

Характерно, что скорость за скачком KF всегда меньше скорости за скачком С К (Я₄<[Я₃); отсюда следует, что линия КН является линией тангенциального разрыва скорости. В вязкой жидкости вдоль КН развивается вихревое движение.

Изменение состояния газа вдоль линии тока при переходе через рассматриваемую систему скачков может быть представлено в тепловой диаграмме (рис. 4-15,6). В точке 2 определяется состояние газа после первого скачка, а в точке 3— после второго скачка, так как углы и известны. Легко определяются также все параметры за скачками: р₃. Т₃, Я₃ и параметры торможения р₀₃, р₀₃.

Точка 3' на изобаре р._л дает состояние газа за скач-ком KF⁵⁴. В точке 0₃ находим давление торможения за скачком KF р_(а <Zp₀₃. Потери энергии в скачке KF выше суммарных потерь в скачках ВК и СК, т. е. As' ]>

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

Предельным случаем является торможение потока вдоль плавной вогнутой стенки, в каждой точке которой поток испытывает отклонение на малый угол d§ (рис. 4-16).

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

Рис 4-16 Сжатие потока при обтекании плавной вогнутой стенки

Поток за скачком вихревой, так как скорости в разных точках за линией ВК различны.

Взаимодействие волны сжатия со скачком уплотнения качественно происходит так же, как это было указано при рассмотрении взаимодействия двух косых скачков (рис. 4-15). В точках пересечения слабых волн и скачка возникают отраженные слабые волны (рис. 4-16), так как изменения давления и направления потока при переходе через волну сжатия и скачок будут различными.

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

Отличный от рассмотренного случай пересечения двух косых скачков показан на рис. 4-17. Косые скачки возникают в результате поворота двух противоположных стенок канала на разные углы ^ и 8_а. Направления потока в зонах II и III будут неодинаковыми-, в зоне II! отклонение будет больше на угол &₂ —8_Х. Параметры течения за косыми скачками АВ и А_ХВ легко могут быть определены по известным параметрам до скачков р_х, 7\ и углам 8j и 8_а, если эти углы меньше соответствующего максимального значения 8_т для данного вектора скорости 1₁.

Рис. 4-17. Схема пересечения двух косых скачков.

а__нор\тлььое пересечение; б — процесс в скачках в тепловой диаграмме.

Параметры потока в области IV можно найти, исходя из граничных условий для линии тока, проходящей через точку В. Примем, что направления скоростей и давления во всех точках области IV будут одинаковыми. Отсюда определяется угол между вектором скорости в зоне IV и вектором 2j. Действительно, если результативное отклонение потока в зоне IV обозначить 8_4> то из рассмотрения рис. 4-17 легко установить, что отклонение потока при пересечении скачка BC_t равно S_t —{— 8₄, а при пересечении ВС оно равно 8₂ — 8₄. Задаваясь различными значениями давления в зоне W (p_lv) по формулам (4-13) и (4-17) или по диаграммам косых скачков, находим углы наклона скачков ВС и BC_t и углы отклонения потока Sj —|—    и

8₂ — 8₄. Значение р_]у, при котором величины 8₄, определенные по параметрам зон II и III, будут одинаковыми,

можно найти, построив зависимости ^_ABCt и 8_4вс от р_п. Точка пересечения этих кривых даст искомое значение 8₄. Зная величину 6₄, находим углы косых скачков р_йс и (З_вс.

Изменение состояния газа вдоль двух линий тока, пересекающих зоны // и III, в тепловой диаграмме показано на рис. 4-17,6. Суммарное приращение энтропии для этих линий тока будет одинаковым только в тех случаях, когда 8₁ = 8₂. Если интенсивность скачков АВ и А_гВ различна, то приращение энтропии вдоль рассматриваемых линий тока будет различным (точки 4 и 4' на рис. 4-17,6).

Рис. 4-18 Мостообразный скачок.

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

Устойчивое существование системы двух пересекающих косых скачков возможно не при всех условиях. Если углы вторых скачков р_вс и |3_ВС будут больше соответствующих значений Р_от, характер течения меняется. Вблизи

нейтральной линии тока, проходящей через точку В, образуется криволинейный скачок. Система пересекающихся прямолинейных косых скачков переходит в мостообразный скачок (рис. 4-18). Потери энергии в потоке при этом увеличиваются.

4-7. СТУПЕНЧАТОЕ ТОРМОЖЕНИЕ ПОТОКА

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

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

Торможение потока в различных системах скачков было подробно исследовано Г. И. Петровым и Е. П. Уховым. Следуя основным выводам этой работы, рассмотрим частный случай торможения потока в двух скачках — косом и прямом.

Рассматриваемая задача формулируется так: определить угол наклона первого — косого — скачка, при котором переход от заданной сверхзвуковой скорости к дозвуковой происходит с минимальными потерями (рис. 4-19). Расчет такой системы можно осуществить последовательно, применяя диаграммы косых скачков (см. приложение). При заданной скорости невозмущенного потока Х_г и выбранном значении угла 8 (или р,) легко определяются скорость Я_а и давление р₂ за косым скачком. Соответствующая потеря энергии или изменение давления полного торможения е₀₁ также определяется по диаграмме косых скачков (или по формулам § 4-5).

Аналогично можно найти скорость и статическое давление за прямым скачком (Я₃ и р₃) и коэффициент потери энергии (или s₀₂).

'В качестве примера на рис. 4-19 показано изменение скорости потока и коэффициентов потерь в системе косого и прямого скачков в зависимости от угла ^ для Я₁ = = 2,0 (? = 1,3). Кривые показывают, что для данной скорости имеется такое наивыгоднейшее сочетание косого и прямого скачков, при котором суммарные потери будут наименьшими.

Действительно-, с увеличением угла косого скачка р, растет коэффициент потерь С_с| в косом скачке и уменьшается скорость за косым скачком Я₂. Очевидно, что при 8 = 0 и Pi=a_ml=22°45'(дляЯ^г.О) косой скачок переходит в характеристику. В этом случае С_с1 = 0. Предельное значение угла р_х, при котором еще возможно существование плоского косого скачка, составляет Pj = P_ml = 65°40\ При этом значении р, поток за косым скачком имеет дозвуковую скорость. В указанных пределах изменений угла [5, (22°45' — 65°40') за косым скачком может существовать прямой скачок. При j5₁ = 22°45' существует только прямой скачок, а при p_t = 65°40^г— только косой.

Меняя в этих пределах угол косого скачка и подсчитывая Я₂ и /7₂ (параметры перед прямым скачком), можно найти параметры газа за прямым скачком. Скорость за

Рис. 4-19. Изменение скоростей и потерь в системе двух скачков (косой + прямой) в зависимости от угла косого скачка р _t при = 2,0; k = 1,3.

прямым скачком Я₃ = -^- увеличивается с ростом в пределах от Л₃ = 0,5 при Pi = 22°45' до Л₃ = Я, = 1 при j3_x= = 65°40\

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

где С_с2 определяется по диаграмме косых скачков для Х_й. На рис.' 4-19 нанесены значения Можно видеть, что С'₂ уменьшается с ростом В_г При Р₁ = 22°45' С^,₂ = С_с2 = = 0,227, а при р^ббЧО' С₂ = °.

Коэффициент суммарных потерь в системе двух скачков, очевидно, будет равен:

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

Кривая суммарных потерь в системе двух скачков имеет минимум при Pi = 45°. Очевидно, что значение ^ является оптимальным в отношении потерь энергии в системе скачков.

Аналогичные расчеты можно выполнить для различных скоростей Я,, определяя наивыгоднейшее значение Pi- Результаты таких расчетов даны на рис. 4-20, где суммарный коэффициент потерь С_с представлен в зависимости от угла p_t для различных значений Х_г. Жирными линиями проведены кривые в диапазоне углов fi„ при которых возможно существование системы косого и прямого скачков. Пунктирная линия ABCD соединяет точки, отвечающие p_i = a_ml. Для этих точек косой скачок имеет бесконечно малую интенсивность и, следовательно, торможение потока осуществляется только в одном прямом скачке. Точки FGHJ отвечают углу = р , при котором поток за плоским скачком имеет звуковую скорость. Для р,> >Р,. кривые С_с проведены тонкими линиями. В этом случае расчет может быть приведен в предположении существования скачка, скорость за которым дозвуковая. При Pj = 90° он становится прямым. Легко видеть, что при Pj = а_т1 и р₁ = 90° коэффициент потерь С_с имеет одинаковое значение.

Сопоставление кривых на рис. 4-20 показывает, что оптимальные значения Р_|опт зависят от скорости невозмущенного потока Х_г. С ростом Х₁ до некоторых пределов значения р_|опт уменьшаются. Для ^ = 1,6 коэффициент

потерь при оптимальном значении $_1опт = 52⁶ составляет = 0,035. В этом случае один прямой скачок дает (точка А на рис. 4-20) С_с2 = 0,113, а один косой при скорости за скачком, равной скорости звука (точка J на рис. 4-20),    = 0,073. Следовательно, переход от одного

Рис 4-20 Кривые коэффициентов потерь в системе двух скачков (косой + прямой) в зависимости от угла косого скачка р., и скорости

I,; k = 1,3

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

^; Следует заметить, что с ростом Я₁ минимум кривых становится более пологим. Это обстоятельство позволяет выбирать оптимальные значения таким образом, чтобы и статическое давление за вторым прямым скачком было наибольшим. Отношение статического давления за системой скачков р_г к полному давлению перед скачком p_oiможно представить в следующем виде:

Рз ___ Р2 Рг Р\ .

Poi Р\ Р2 Рт

здесь — характеризует повышение статического давления Р1

на косом, а — — на прямом скачке.

Рг

Изменение этих величин, а также е_П1=^-² и е„„ = ^-³

“    Р*1    ⁰²

в зависимости от угла косого скачка (3, для Х_г = 2,0 представлено на рис. 4-21. С ростом р, отношение давлений

на косом скачке — увеличивается, а на прямом — умень-Р1    Рг

шается. График показывает, что относительное статическое давление за системой скачков для Х₁ = 2,0 имеет максимальное значение при р₁ = 40°, тогда как минимальная величина была получена при [3, — 45°.

Учитывая, что кривые в области минимума протекают полого, оптимальные значения можно выбирать по данным расчета восстановления статического давления в си-

Рис. 4-21. Изменение статического давления и давления торможения в системе двух скачков (косой -f- прямой) в зависимости от угла косого скачка Pj для А., = 2,0; ft = 1,3.

М-еме скачков, т. е. принимать значения р_1опт несколько меньшими, чем это диктуется кривыми

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

Рис. 4-22. Процесс торможения в системе двух скачков (косой + прямой) в тепловой диаграмме.

Последовательное торможение потока в системе скачков можно наглядно представить в тепловой диаграмме. На рис. 4-22 показан этот процесс для двух скачков. При р, = а_т1 существует только один прямой скачок (точка Л). При p_t = Pj скорость за плоским косым скачком звуковая и торможение происходит только в одном косом скачке (точка К). В пределах а_т,<р,<р, может существовать система двух скачков. Задаваясь различными значениями pj (точки E_l — EJ, можно построить линию предельных состояний газа за системой двух скачков при различной интенсивности косого и прямого скачков (линия АВ_г.. .В₄К)-Малым углам косого скачка отвечает верхняя ветвь кривой предельных состояний AB_XB_ZB_%. Нижняя ветвь соответствует большим углам скачка Р^Р, • Характер кривой состояний АВ₁... В^К отчетливо указывает на существование оптимальной комбинации двух скачков. При некотором значении угла Р_!опт потери энергии в системе скачков оказываются минимальными (точка В₃). При другом, близком значении угла ^ потери несколько выше (Д/г_к+п> >Д/г_к+п), но статическое давление достигает максимально возможного значения /?_3макс (точка В₂). Линия предельных состояний может быть построена путем последовательных расчетов системы скачков.

Рис. 4-23. Изменение потерь в скачках при ступенчатом торможении (цифры указывают число скачков).

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

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

Графики, приведенные на рис. 4-23, отчетливо показывают преимущество более сложных систем скачков при больших сверхзвуковых скоростях. Кривые ?_с^={.(Х₁) позволяют выбрать наиболее рациональную схему ступенчатого торможения для заданной скорости.

4-8. ОТРАЖЕНИЕ СКАЧКОВ

а) Отражение от твердой стенки

Рассмотрим отражение косого скачка от прямой твердой стенки, расположенной параллельно направлению скорости невозмущенного потока (рис. 4-24). Скачок образуется в

Рис. 4-24. Схема нормального отражения плоского косого скачка от твердой стенки.

А    а    У-тУ*⁸

¹    fiz    Л₃

Рз^аз

77777777777777?

А

тйчке А, где стенка поворачивается на угол 8. При переходе через первичный скачок АВ линия тока отклоняется к прямой стенке на угрл 8. Очевидно, что в точке В этот поворот неосуществим и граничная линия тока сохраняет направление стенки. Это означает, что в точке В стенка принудительно отклоняет линию тока в обратном направ-

Ленин на угол 8. В результате возникает отраженный косой скачок ВС. Заметим, что углы падающего и отраженного скачков неодинаковы, так как перед скачком ВС безразмерная скорость Я₂ Л j при том же угле отклонения 8. Из графика (3 = /(8, Х_г) (рис. 4-5) видно, что угол

Расчет отраженного скачка не встречает затруднений. Зная параметры [_невозмущенного течения А,, р_г и угол

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

а—определение углов и скоростей при неправильном отражении; б — А-образный скачок при отражении.

отклонения 8, с помощью диаграммы скачков легко определяем параметры потока за первичным скачком: Я₂, р_г и (З_а. При том же значении 8 находим состояние газа за отраженным скачком: Я_а, p_s и (5₃. Изложенным выше способом находим потери энергии в первичном и отраженном скачках. Следует иметь в виду, что такое отражение косого скачка возможно не всегда. Если "уго'л отклонения 8 больше максимального значения для скорости А₂, то картина отражения меняется. Допустим, что в диаграмме ударных поляр (рис. 4-25,а) отрезок OD изображает скорость потока до скачка А,. Если угол отклонения стенки 8<8_ш2, то гипоциссоида, отвечающая скорости за скачком Я_а (отрезок 02), пересекает линию вектора (точки 3 и 4).

При 8 = Ь_т2 линия OD касательна к ударной поляре Я₂ (точка 3'). Картина течения при этом остается прежней (рис. 4-24). Вектор скорости за вторым скачком 03' (Я₃) будет меньше единицы (скорость дозвуковая).

Если 8>8_т2, то ударная поляра, построенная для скорости *₂, не имеет общих точек с вектором OD и отраженный скачок не может обеспечить выравнивание потока. Часть потока, прилегающая непосредственно к стенке, становится дозвуковой. Отраженный скачок CD искривляется и сдвигается против течения. При этом деформируется и первичный скачок АВ. Элемент СВ этого скачка становится нормальным к стенке, система скачков приобретает Я-образную форму. За участком прямого скачка поток дозвуковой. За криволинейной частью отраженного скачка поток может быть сверхзвуковым. В результате образуется линия тангенциального разрыва СЕ, по обе стороны которой давления одинаковы, а температура и плотности различны, так как изменения энтропии при переходе через ВС и АС — CD будут различными. В области 3 за Я-образным скачком поток вихревой.

Отражение скачка, показанное на рис. 4-25,6 может возникнуть при значительных углах отклонения 8 и малых сверхзвуковых скоростях Я_х. Подчеркнем, что механизм возникновения Я-образного и мостообразного скачков одинаков.

б) Отражение от свободной границы струи

Такое отражение рассмотрено на рис. 4-26. Во всех точках на границе FBE давление одинаково и равно давлению внешней среды р_а. В струе это же давление имеет место только до скачка АВ. При переходе через скачок АВ давление изменяется от р_х = р_а до p_t > р_а. Следовательно, точке В свойственны одновременно два давления: р_а со стороны среды и р_г со стороны струи. Такая точка является очагом возмущения сверхзвукового течения, создающим стационарную волну разрежения. При обтекании точки В давление потока должно упасть от р₂ до р_а, что и приводит при сверхзвуковых скоростях к образованию волны BCD.

Первая характеристика BD составляет с направлением вектора М₂ угол а_т2 = arcsin ~, где М₂ — скорость потока за скачком АВ (М₂ < M_t). Угол последней характеристики а_т3 arcsin Здесь скорость за отраженной

Рис. 4-26. Отражение косого скачка от свободной границы струи.

волной разрежения М₃ определяется по отношению — ,

Рог

где р₀₂— давление торможения за косым скачком.

Отражение скачка приводит к деформации границы струи, которая в точке В отклоняется на угол 8₂>8,.Это отклонение вызывается расширением струи. Таким образом, при отражении от свободной границы струи, вдоль которой давление сохраняется постоянным или падает, скачок уплотнения преобразуется в волну разрежения. Если давление вдоль границы возрастает, то в зависимости от интенсивности изменения давления отражение может быть погашено или оно происходит с сохранением знака (как и от твердой стенки).

4-9. ВЗАИМОДЕЙСТВИЕ СКАЧКА И ВОЛНЫ РАЗРЕЖЕНИЯ

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

Рассмотрим в качестве примера обтекание заостренной пластинки (рис. 4-27). На переднем остром (8<§_т) носике пластинки возникает плоский косой скачок АВ. При обтекании точки D образуется стационарная волна разрежения, причем характеристика, на которой начинается отклонение потока, расположена под углом

где М₂ — скорость за скачком.

и волны разрежения.

Так как

^йт2“И>Рп

то характеристика пересечет скачок уплотнения в некоторой точке В. Второй границей волны разрежения является

. 1

характеристика, расположенная под углом a_m3 = arcsin —.

На участке правее точки В волна разрежения взаимодействует с косым скачком. В области ABD скорость постоянна и равна М_г; линии тока параллельны образующей клина AD. Проведем через точку В характеристику невозмущенного потока под углом a_m] = arcsin к направлению вектора скорости с_г. Так как косой скачок располагается под углом (3. (а_ш1 -(- a_mi -(- S.), а угол ближайшей к DB волны DE меньше а_т2, то на участке BE косой скачок отклоняется на малый угол так, что с уменьшением угла наклона скачка (3j уменьшается и угол отклонения потока S.. Соответствующие отклонения скачка и изменения угла поворота потока 8 имеют место и на участках EF, FG и т. д.

Следовательно, скачок, начиная от точки В, искривляется и отклоняется в направлении потока; угол скачка уменьшается, приближаясь к а_т1.

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

Аналогичные результаты получаются, если волна разрежения расположена перед скачком (рис. 4-27,6). В точке Dj возникает волна разрежения, а в точке А_г — скачок уплотнения. Взаимодействуя с волной разрежения, скачок искривляется. Так как после пересечения с последней характеристикой волны разрежения D₁B₁ скачок А₁В₁ попадает в зону меньших скоростей, угол его (3_; увеличивается.

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

В предыдущих параграфах настоящей главы былй рассмотрены скачки уплотнения в плоском течении. При обтекании осесимметричных тел поверхности разрыва имеют осесимметричную форму.- Рассмотрим особенности осесимметричного скачка на примере обтекания кругового конуса (рис. 4-28). Перед конусом образуется конический скачок, вершина которого совпадает с вершиной конуса, если угол раствора конуса меньше максимального значения для данной скорости набегающего потока.

Рис. 4-28. Форма линий тока в возмущенной области за коническим скачком при обтекании конуса.

Основные соотношения при переходе через поверхность конического скачка, как легко видеть, будут теми же, что и для плоского скачка [уравнения (4-13), (4-14) и др.]. При одинаковых углах раствора клина и конуса скачок на конусе будет иметь меньший угол наклона, чем на клине, так как конус вызывает меньшие стеснения потока, чем клин бесконечного размаха того же угла раствора. При переходе через конический скачок линии тока, так же как и в случае плоского скачка, претерпевают излом. Однако так как скачок на конусе слабее, чем на клине, непосредственно за скачком линии тока будут наклонены к вектору скорости невозмущенного потока под углом, меньшим угла раствора конуса у₀. Расчеты показывают, что в возмущенной области линии тока не являются пря-

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

Из рис. 4-28 следует, что с удалением от скачка угол наклона линий тока к оси конуса увеличивается и линии тока асимптотически приближаются к направлению, заданному образующей конуса. Здесь можно видеть, что кольцевая трубка тока, образованная двумя смежными линиями тока,

Рис 4-29 Схемы спектров обтекания конуса при различных скоростях невозмущенного потока.

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

заключить, что при одинаковом отношении давлений у

(/?_а—давление на поверхности обтекаемого клина или конуса) сжатие газа при обтекании конуса происходит с меньшими потерями, чем при обтекании клина, так как для конуса полное повышение давления является суммой изоэнтропи-ческого сжатия и сжатия по ударной адиабате.

В сверхзвуковой части возмущенной области благодаря кривизне линий тока характеристики криволинейны.

На рис. 4-29 изображены три возможных случая: скорости во^- всех точках возмущенной области меньше скорости звука (а), случай смешанного течения, когда скоро-рости непосредственно за скачком сверхзвуковые, а затем становятся дозвуковыми (б), и, наконец, когда поток за скачком полностью сверхзвуковой (б). Характер течения за скачком при неизменной величине у₀ зависит от скорости набегающего потока.

Рис 4-30. Зависимость скоростей на поверхности конуса от скорости невозмущенного потока и угла конуса

Скорости и давления на поверхности конуса меняются при изменении скорости невозмущенного потока и полуугла конуса Yo-

На рис. 4-30 приведены графики изменения безразмерной скорости на поверхности конуса Я₂ в зависимости от угла Yo ^и Следует заметить, что при данном режиме обтекания во всех точках конической поверхности скорости и давления имеют постоянные значения. Для конуса, так же как и для клина, теоретическое решение задачи дает при одном и том же режиме два возможных значения угла наклона скачка и параметров у поверхности конуса (нижние ветви кривых на рис. 4-30). Однако практически, как правило, реализуются меньшие значения углов скачка (верхние ветви кривых на рис. 4-30). Поэтому можно сделать вывод, что с увеличением у₀ скорости на поверхности конуса уменьшаются и давления растут. Увеличение скорости невозмущенного потока приводит к противоположным результатам.

Такой характер изменения параметров потока на поверхности конуса имеет место до тех пор, пока угол

конуса не достигнет предельного значения, при котором происходят отход и деформация скачка, так же как и в случае клина. При этом конический скачок преобразуется в осесимметричную поверхность разрыва с криволинейной образующей. Однако максимальные полууглы конуса Y_0m; при которых для данного Я\ происходит преобразование конического скачка 6 осесимметричный с криволинейной образующей, будут больше соответствующих значений Ь_т для клина. На рис. 4-31 приведены зависимости максимальных углов отклонения Ь_т (и Тот ^{0Т числа} для клина и конуса.

Рис. 4-31. Зависимость максимальных углов клина и конуса от скорости невозмущенного потока (ft = 1,4).

Для конического скачка можно построить в плоскости годографа и, и и в тепловой диаграмме ударную поляру (рис. 4-32).

В плоскости годографа изменение скорости непосредственно в коническом скачке изображается линией DEA, причем вектор скорости за скачком определяется отрезком ОЕ (скорость невозмущенного потока OD). Угол скачка р можно найти, проведя нормаль в точке Е к отрезку DE. Изменение скорости в возмущенной области за скачком описывается кривой ЕЕ_Х. Эта линия отвечает изоэнтропическому изменению скорости (сжатию) за скачком.

Яблоковидная кривая DE_XA определяет годограф скорости на поверхности конуса; ее можно назвать ударной полярой конуса. Наклон отрезка OE_t определяет полуугол конуса Yo- Область, заключенная между кривыми DEA и DEj^A, характеризует поток в возмущенной области. В любой точке N отрезок ON определяет величину и направление скорости. Нормаль, проведенная к годографу скорости в точке N, дает полуугол конической поверхности, проходящей через эту точку в плоскости потока.

Каждая промежуточная кривая DNA_t соответствует постоянному значению разности углов 8—у Так как в возмущенной области давление полного торможения не меняется, то годографу скорости ЕЕ_г

отвечает постоянное значение . Нанося эти значения для различно 1

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

Рис. 4-32. Ударная поляра для конического скачка в плоскости годографа и в тепловой диаграмме для 6=1,4 (яблоковидные кривые).

6)

мов обтекания конуса с дозвуковыми скоростями за скачком. При этом легко устанавливаются точки возмущенного потока, в которых скорость течения равна критической. Для данного угла скачка (i эти точки получаются пересечением дуги д* с годографом скорости Е'Еу (точка К).

В тепловой диаграмме ударная поляра строится уже известным нам способом (рис. 4-32,6). Линия ПЕ_гА соответствует изменению состояния газа за коническом скачком при изменении р₀ в пределах от я

Ро⁼ “о (точка D) до р_о = g- (прямой скачок). При определенном значении Ко ^и соответственно р_о состояние потока непосредственно за скачком характеризуется точкой E_v определяющей изменение энтропии (потери Ah) и изменение потенциальной энергии газа в скачке (#_п1).

Отрезок отвечает изоэнТропическому сжатйю за скачком, и в точке Е можно найти параметры газа на поверхности конуса. Соответствующее изменение потенциальной энергии равно Н_п2- При одинаковых углах плоского и конического скачков (Р = р_о) изменения параметров получаются близкими, так как изоэнтропическое сжатие в возмущенной области значительно менее интенсивно, чем ударное сжатие в скачке.

В системе конических скачков можно осуществить последовательное торможение сверхзвукового потока, так же как и в системе плоских скачков.

В заключение необходимо сделать следующее замечение.

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

Рис. 4-33. Толщина скачка в зависимости от его интенсивности.

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

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

Расчеты показывают, что толщина скачда имеет поря-1.док удвоенной длины свободного пробега молекулы и поэтому уменьшается с ростом его интенсивности. Соответствующий график изменения толщины скачка в зависимости

от — при не слишком малых давлениях представлен на

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

4-11. СКАЧКИ КОНДЕНСАЦИИ (ТЕПЛОВЫЕ СКАЧКИ)

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

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

Анализ свойств скачков конденсации основывается на некоторых допущениях: а) конденсация происходит мгновенно, так что образуется! резкая граница, отделяющая газ с несконденсировавшимися водяными парами, от газа, содержащего конденсат; б) эффект конденсации сводится к освобождению скрытой теплоты парообразования; в) этот процесс сопровождается изменением физических свойств газовой составляющей и уменьшением ее весовой доли в смеси; изменение физических свойств газа и его параметров происходит только в пределах скачка; г) влиянием вязкости, теплопроводности, диффузии можно пре-

р

небречь; д) газовая фаза подчиняется уравнению состояния —=gRT

и k меняется только прн переходе через скачок; е) после скачка жидкая фаза имеет ту же скорость, что и газообразная.

Основными уравнениями скачка конденсации ⁵⁵ являются общие уравнения, выведенные в § 4-2.

С учетом обозначений, принятых на рис. 4-34, уравнение неразрывности запишем в следующем виде:

Pi^c«l ⁼.Рг(1 +/)^сл2 = Рг С¹ + /Н^йг sin Р_к— o₂cos Р_к), (4-37)

где    р_к—''угол косого скачка конденсации-

Р X

/ = —----отношение массы жидкости к массе газа за скачком:

Ра 0 — 1)

X — степень влажности за скачком; р_н — плотность жидкости.

^C«2^=a2Sin Р_к—y₂COS Р_к.

Рис. 4-34. К выводу уравнений косого скачка конденсации.

Уравнения количества движения в проекциях на нормаль к скачку и на плоскость скачка будут:

Рг — Р 1 = Pi^! — Рг (1 + /) («2 sin Р_к— v₂ cos Р_к)⁵⁶;    (4-38)

?i^cnl^ct = Pzi¹ i-j)^cn2^ct =

= Рг U + i) («2 sin P_K— v₂cos P_K) (u₂cos P_K+ v₂sin y=0,    (4-39)

где c_t= c, cos g_K = m₂ cos + v₂ sin f_K.

Уравнение энергии запишем в такой форме: с²    с²

_J_    (, = _JL /₂ — Д(₀,    (4-40)

где Дi₀ = (₀₂ — i₀₁ — изменение энтальпии торможения вследствие выделения теплоты при конденсации.

Так как энтальпия торможения при пересечении скачка меняется, то критические скорости до и после скачка будут неодинаковыми:

-«-К⁵⁷!

+ 1 ‘°²

Здесь k_v k₂ — показатели изоэнтропического процесса до и после конденсационного скачка.

Соотношение между критическими скоростями выражается формулой

^д2' ^К    (4-41)

a² = “il =

+

где

— 1 k₂ +

—    A i

а/₀ = т^; К = г

1 •

1 +

Pi

'I

а,

Р, 2(fe, — 1)

Левую часть уравнения энергии (4-40) можно представить в следующей форме

„2

(4-40а)

^{+ 1}    а%=    (₀₁(1    +    дг₀)    (4-406)

а правую часть «2 + ^и2 ,

Рг

-1р₂(1+/) 2 (*_ж — 1)

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

(4-42)

И tg!

^к и.

где 5_К— угол отклонения линии тока при пересечении скачка конденсации, с помощью выражений (4-37), (4-38), (4-40), (4-41) и (4-42) получим уравнение ударной поляры для скачка конденсации в таком виде:

и_а(Х, — и,)-

—2    /•>    —у> ^1 ^М2

^v2 — ( 1 ^й2/ 2

_J_ 1 ^1 - ^1^М2 +

— (k_t — 1)^2 — (fe, + 1) (1 + Д/₀)

— 1

- Д(₀ — k----

(*,-*)(*?- >)

(4-43)

₊1 - ₊ 1

, ^С1 - ^Ц2 - ^V2 — ,

“²⁼⁼^^; *=*!-*.¦

Уравнение (4-43) при Дг₀ =* Ъ = 0 переходит в формулу (4-26) для адиабатического скачка.

С помощью (4-43) по формулам (4-42) определяются угол косого скачка конденсации Р_к и угол отклонения потока 8_К.

Из уравнения энергии (4-406) можно получить формулу для определения числа М₂ за скачком:

м| + у|

Аналогично тому, как это сделано при выводе формулы (4-35), получаем формулу для отношения давлений торможения на скачке конденсации:

k_t— к

2 (Л, — 1) ' «₂² + .₂

Отношение давлений на скачке находим по формуле

__1 + k_tM²_{ sin р_к_

(fe, —~k) (1 + l)Ml -    „    •    -    ГГ

¹ + --_У ¦    (“г    sin    р_к—    v₂    cos    ?_к)²

(1+/) = 1 + Д(₀

(1 + k) (k, —

+

v^:

-\-(ki i k) -

-г Mi

Af,=

1

(4-44)

¦к)

(4-45}

ft,—i-,

Рог_Рг

Pot 1*1

(4-46)

k,~i

Af?

+ (&i — ^“2”

Отношение плотностей на скачке

Ра _Рг(МЛ* h Pi Pl\Mi) й²+р²

(4-47)

Полученная система шести уравнений (4-42) — (4-47) позволяет найти параметры за скачком v₂, и₂, М₂, р₀₂, р₂ и tg Р_к при заданных

параметрах до скачка с„ />,, р„ а_г1 и известных Ai₀, /, k.

С помощью общих соотношений нетрудно рассмотреть отдельные частные случаи.

go многих практически важных случаях можно полагать, что физические свойства газа сохраняются неизменными (к-, = k₂\ k = 0)n что масса конденсированной фазы за скачком пренебрежимо мала по сравнению с газообразной (/ = 0) В этом случае уравнение (4-43) упростится и примет вид:

М ¹¹

^ = (*.-«₂)²-\-¦*'    -----(4-48)

?ТрГ    + 1

Совместное решение исходных уравнений после соответствующих упрощений позволяет получить связь между нормальными составляющими скорости на скачке конденсации в следующей форме:

7_л2=*.А+1+|/    Y-i.    <⁴-⁴⁹>

Ь_{г 2} =

_₂ ь\ ^l~k + T^ai^cos2?k

......*~ГГ»-•    ⁽⁴'⁵⁰⁾

²    -=^-    -    г-гт х? cos² 3„

: “Ь 1

а:

Величина а, определяется по формуле (4-41) Для отношения давлений на скачке находим:

?¦» -M(k-A_M2 _ialt р 1 Yjl.^«l_A + i.

л *-i\ 2    Дт.    7_л2

Это уравнение переходит в формулу (4-13J при ~b, = 1 (адиабати-

С„1

ческий скачок) после подстановки —— из уравнений (4-5) и (4-11).

^сп2

Из уравнения (4-48) следует, что вертикальная составляющая скорости за скачком обращается в нуль при трех значениях вектора и₂.

Первое соответствует Вырождению скайка а слабую волну [и_г = = Х,(Х₂ = a»Xj)]. Второе и третье значения получаются из условия

или

Два корня этого уравнения будут:

т. е.

д!(Х? + 1)²_, ₍₄.₅₁₎

/

4Х?

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

тс —    _

ставив в уравнение (4-49) условия прямого скачка: $_к=2 и Ь* = а* .

Из уравнения (4-51) вытекает, что скорость за прямым скачком зависит от X, и Дг₀ — теплоты, выделяемой при конденсации, которая в свою очередь определяется количеством конденсирующегося^газа.

Из уравнений (4-49) и (4-50) следует, что величины b* и а. не могут быть меньше некоторого предельного значения для заданного с_п1 илн Х„ так как в противном случае с_п2 и Х₂ будут мнимыми величинами. Из (4-51) для прямого скачка

^Д*мин _Х2 ₊ ! •

В согласии с формулой (4-41) при /г, = /г₂ минимальному значению а, соответствует максимальное изменение энтальпии торможения в скачке:

(4-52)

В общем случае для косого скачка из уравнения (4-49) получим:

соответственно из (4-50) будем иметь:

_2(4 + 1)

/ [k+l +Х²(2/г stn² Р_к + l-ft)]²

V

• 4 (*² — 1) X, cos²

Xfsin²

(4-53)

— at

ДI

Очакс —о

at,

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

Возвращаясь к анализу уравнения ударной поляры (4-48), отметим, что у₂= со при

4-( —-—X? + 1 ) или -^²-=_1J —-—X'f + 1 I. X, U + 1 ¹    '    "    I

² >2 й*2 "Н ^

Зависимость (4-48) графически представлена на рис. 4-35 для раз личных значений а». Для плоского косого конденсационного скачка

Рлс. 4-35. Ударные поляры конденсиционных скачков для различных значений Д/₀ (различной относительной влажности); Х, = 1,5‘.

k= 1,4.

при данном а.(Дг₀) на ударной поляре находим две точки: Е₂₁ и ?₂₂, отвечающие двум различным углам скачка. Причем точка Е₂₁ соответствует криволинейному конденсационному скачку. Точки Е₃ отвечают скачку разряжения.

Ударной поляре адиабатического скачка соответствует значение а, = 1. По мере уменьшения а, (увеличения Дг₀) угол скачка конденсации при данном 5_К возрастает. Прямому конденсационному скачку

в соответствии с формулой (4-51) отвечают два значения вектора и₂ (точки D, и D₂)- Отрицательный знак перед корнем в (4-51) дает точку D,, а положительный — О_г.

Таким образом, уравнения (4-43) и (4-48) и ударная поляра на рис. 4-35 показывают, что теоретически возможно существование четырех типов прямых и косых конденсационных скачков, отвечающих различным скоростям перед скачком и величинам «ДДг₀).

Соответствующая классификация указана в табл. 4-2.

Таблица 4-2 Возможные типы конденсационных скачков

Однако если учесть некоторые дополнительные условия, то возможными оказываются скачки всего двух типов: 1) сверхзвуковые скачки, в которых c_n^>a_v с_п2^а₂ и конденсация сопровождается сжатием газа (р_г /?,); 2) дозвуковые скачки, отвечающие условиям с_п'2^<^^а* ^ив которых конденсация сопровождается разрежением газа (Pi<iPi). Скачки, отвечающие соотношениям с_л1>д, и с_п2^> д₂> не могут реально осуществиться, так как такой скачок перемещался бы относительно находящегося перед ним газа со сверхзвуковой скоростью и его возникновение не должно было бы отразиться на состоянии этого газа.

Скачки конденсации в дозвуковом потоке не могут перевести газ в область сверхзвуковых скоростей    с_п2^    а₂),    так    как    в    этом

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

Определим изменение энтропии в конденсационных скачках.

В частном случае k — 0 получим:

k

A_S=^[(l+^₀)^fe-‘g-;].    (4-54)

В формулу (4-54) входит отношение которое можно найти по

Р 02

-    Рг

уравнению (4-46), подставляя k = 0; прн этом М₂ н — определяются

по формулам (4-44) и (4-45).

Исследуя уравнение (4-54), можно убедиться в том, что для скачков первых трех типов (табл. 4-2) As ]> 0. Однако, привлекая дополнительное условие (термодинамическое состояние перед скачком должно соответствовать началу быстрой конденсации) н учитывая влияние теплообмена, можно показать, что скачок первого типа, как и четвертого, невозможен.

Опыт подтверждает возможность образования скачков второго и третьего типов.

На положение скачка, его форму и интенсивность решающее влияние оказывают влажность воздуха и скорость потока. На рис. 4-36 приведены кривые относительного давления и числа М, перед скачком конденсации в зависимости от абсолютной влажности воздуха х по опытам А. А. Степчкова, проведенным в соплах Лаваля. С увеличением влажности скачок конденсации перемещается в область меньших чисел М_г.

Рис. 4-37.

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

б)

С увеличением влажности уменьшается пересыщение потока воздуха парами воды, определяемое- отношением парциального давления

Р\п

пара к давлению насыщения —— , а также переохлаждение Ду =

Р !н

= Г_а1—Г, (рис. 4-36), где j_h1—температура насыщения, а Г,—температура перед скачком. По мере возрастания влажности величины пересыщения и переохлаждения уменьшаются.    -

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

В заключение отметим, что изложенная теория конденсационных скачков оставляет в стороне вопросы механизма конденсации, зарождения и развития ядер конденсации.

ГЛ АВА ПЯТАЯ ДВИЖЕНИЕ ГАЗА ПРИ НАЛИЧИИ ТРЕНИЯ

5-1. ТЕМПЕРАТУРА ТОРМОЖЕНИЯ В ВЯЗКОЙ ЖИДКОСТИ

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

Движение вязкой жидкости описывается системой уравнений сохранения: расхода, количества движения и энергии. Уравнение неразрывности (1-12), как уже указывалось, справедливо и для вязкой жидкости. Уравнения количества движения в форме Эйлера (1-16) должны быть дополнены членами, учитывающими влияние вязкости.

Следует подчеркнуть, что для необратимых процессов движения интегралы уравнений движения и энергии не-совпадают. При выводе уравнения энергии для струйки (§ 2-1) указывалось, что оно справедливо и для адиабатических (необратимых) течений. Однако это замечание вполне справедливо только в частном случае, когда работа сил трения полностью преобразуется в тепло. Такой процесс соответствует простейшей схеме одномерного потока или движению газа с равномерным полем скоростей.

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

196

6П2.2 Дейк Михаил Ефимович

Д 27 Техническая газодинамика. Изд. 2-е, переработ. М.—Л. Госэнергоиздат, 1961

с черт. и илл.    6П2.2

design pashaok

Редактор Б. Я- Шумяцкий    Техн. редактор А. М. Фридкин

ГЛ АВА ПЯТАЯ ДВИЖЕНИЕ ГАЗА ПРИ НАЛИЧИИ ТРЕНИЯ

5-1. ТЕМПЕРАТУРА ТОРМОЖЕНИЯ В ВЯЗКОЙ ЖИДКОСТИ

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

Движение вязкой жидкости описывается системой уравнений сохранения: расхода, количества движения и энергии. Уравнение неразрывности (1-12), как уже указывалось, справедливо и для вязкой жидкости. Уравнения количества движения в форме Эйлера (1-16) должны быть дополнены членами, учитывающими влияние вязкости.

Следует подчеркнуть, что для необратимых процессов движения интегралы уравнений движения и энергии не-совпадают. При выводе уравнения энергии для струйки (§ 2-1) указывалось, что оно справедливо и для адиабатических (необратимых) течений. Однако это замечание вполне справедливо только в частном случае, когда работа сил трения полностью преобразуется в тепло. Такой процесс соответствует простейшей схеме одномерного потока или движению газа с равномерным полем скоростей.

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

В качестве примера рассмотрим движение вязкой сжимаемой жидкости между двумя плоскими стенками

Рис. 5-1. К выводу уравнения энергии для потока сжимаемой вязкой жидкости между двумя плоскими стенками (безградиентное течение).

(рис. 5-1). Верхняя стенка движется в направлении оси х с постоянной скоростью, равной скорости газа с . На нижней стенке скорость равна нулю, так как эта стенка неподвижна⁵⁸. Полагаем, что давление сохраняется постоянным вдоль осей х и у, т. е.

₌др_ ___п дх ду

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

, ^_ р const

W(y)—T(ff-

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

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

Обозначим: Q—’Количество тепла, переданное элементу в единицу времени от соседних частиц; т — напряжение трения. Количество тепла, полученное элементом, определяется как разность (рис. 5-1):

Qdx — (Q-\-^^dyj^dx = ~^ ^dydx*-

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

Тогда уравнение энергии будет:

[(^х+|^) [^c + 7i^dy)

нергии

¦(-Q + w) = o,

dx (тс) dydx.

d dy или

— Q-\-zc = const.

Постоянная в правой части уравнения энергии опреде^ ляется из граничных условий. Так, при у = 0, с = 0 и Q = Q₀, где Q₀ —удельное количество тепла, передаваемое потоку газа от внешнего источника.

Следовательно,

— Q-\-vc = — Q₀.

Для ламинарного движения г определяется по формуле (1-4). Вспомним, что

г)__2    ^

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

После несложных преобразований получаем: л dT. dc    d    (^ср^т    ,    с²    Л    ^

¹ df +^ df — 57VTT+T')⁼ ~ ^Qo'

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

Величину

i = l(T) и |* = [i(7’).

Интеграл уравнения энергии для случая = const позволяет связать статическую энтальпию в потоке с энтальпией у неподвижной стенки следующим образом:

У

о

и

У

о

о

где при линеином распределении скоросп

т₀—напряжение трения на стенке; г* = с Т₀ — энтальпия торможения на стенке.

Следовательно,

;-/; + Рг^=-Рг^.    (5-1)

Для верхней стенки, движущейся вместе с потоком со скоростью с_то, нетрудно получить:

(5-1а)

С помощью (5-1) и (5-1а) после несложных преобразований находим:

Найдем для случая адиабатического течения (Q^^O) отношение температур торможения на движущейся стенке и в произвольном сечении потока (на расстоянии у), учитывая, что температура торможения и термодинамическая температура связаны соотношением

(5-1 г)

Имея в виду, что рассматривается случай с = const, следовательно i — с Т, а для адиабатического течения и Q_o = 0, нз (5-1) и (5-1г) получим:

Для движущейся стенки с помощью уравнений (5-1а) и (5-1 г) будем иметь:

Следовательно,

Y~— 1 ^— 0 — Рг)

¹ О

где

Формула (5-2) показывает, что для Рг ф 1 в вязком газе температура торможения, т. е. полная энергия не сохраняется постоянной по сечению. При Рг = 1 температура торможения Т'₀ = Т₀ = const для всех точек потока.

Число Рг характеризует соотношение между теплом, выделившимся вследствие трения, и теплом, отведенным от элемента теплопроводностью. При Рг<1, что имеет место для всех газов, отвод тепла совершается более интенсивно, чем его выделение. В этом случае Т*₀ <С.Т₀. При Рг > 1 выделение теплоты трения происходит более интенсивно, чем ее отвод, и Т’₀^>Т₀.

Для совершенного газа число Рг является физической константой, не зависимой от состояния газа.

Для более общего случая плоского потока газа, когда скорости зависят от дг н у, дифференциальное уравнение энергии может быть представлено в таком виде⁵⁹:

№ (^д2Г (^! +^Рг т ^М2) ^д2Т(¹ +^Рг т ^М2)^1 _

Рг \ дх⁶⁰    dy⁶¹ J

дТ {1+^^-М²)

= “    Тх    +^у    ду    •    ⁽⁵'³⁾

Уравнения количествS движения с добавлением членов, учитывающих влияние вязкости (уравнения Навье — Стокса), записываются в следующей форме:

ди ди ди    1    др    Г    [д²и д²и\ ,

ЪГ^+ид7 ^{+ v}dy    ^v [ (^^ + ^²J ⁺

1 д /ди dv '

/ d²v d²v\ 1 d (du [ dv \~| \d^⁺di²) )\-

+ v

Эти уравнения дополняются уравнением неразрывности (1-12) для плоского потока:

дР . д (рц) д(ру) dt dx ду

= 0;    (5-4    а)

уравнением состояния (1-1) и уравнением тренчя, например для ламинарного течения законом Ньютона (1-4).

При исследовании движения газа в трубах и каналах с учетом вязкости, а также при изучении обтекания тел газовым потоком задача сводится к определению потерь энергии и аэродинамических сил, действующих на обтекаемую поверхность С этой целью необходимо решить совместно замкнутую систему шести уравнений (5-3), (5-4), (5-4а), (1-1), (1-4), определяя неизвестные функции координат: р, р, и, v, Т ит (для установившегося потока).

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

I

II

Рис. 5-2. Схема подобных потоков*

Широко применяемый в механике метод подобия позволяет сформулировать указанные условия модельных испытаний и устанавливает приемы переноса результатов лабораторных исследований на натурные объекты.

Аэродинамические силы, действующие на обтекаемое тело или на стенки канала (в том числе и силы сопротивления), выражаются через безразмерные коэффициенты. Установим, от каких параметров в общем случае зависят коэффициенты сопротивления.

С этой целью рассмотрим обтекание жидкостью двух геометрически подобных тел (рис. 5-2).

В случае кинематического и динамического подобия двух рассматриваемых явлений поля скоростей и сил в двух потоках должны быть взаимно пропорциональными. Тогда, вводя масштабы длин L, времени Т и массы М, можно представить связь между длинами, временами и массами двух подобных потоков следующим образом: V = Ы; t' = Tt; m.'= Mm.

Масштабы всех других величин, входящих в уравнения (5-4), легко выражаются с помощью указанных масштабов.

Допустим, что уравнения (5-4) выражают связь между параметрами течения первого потока. Имея в виду размерности величин, входящих в эти уравнения, отраженные в единицах измерения: длины л:, у[м], скорости и, v[M/ceK], плотности р[кг-сек²/м*], кинематической вязкости v [м²/сек}, массы т[кг- сек²/м) и давления р [кГ/м²], легко связать параметры первого и второго потоков с помощью масштабных коэффициентов:

х __j ¦ ь _JL_ ^и _ ^v _L

х' ~~у^г ’ ^с ~с^г~~йГ~~1Г    Т~’

Здесь и выше штрихом обозначены параметры течения второго потока. Соотношения между параметрами (5-5) являются очевидным следствием пропорциональности линейных размеров, времен и масс двух подобных потоков.

Располагая оси координат в обоих потоках одинаковым образом, путем подстановки соотношений (5-5) в первое уравнение (5-4) можно записать для второго потока:

+?’>[(S+$)+4 ?(?+&)]¦ <^м>

Уравнение (5-6) выражает уравнение количества движения в дифференциальной форме для второго потока, записанное через параметры первого потока.

Все члены уравнений (5-4) и (5-6) имеют, естественно,

одинаковую размерность [м/сек²], в чем легко убедиться с помощью соотношений (5-5). Чтобы потоки были динамически подобными, необходимо, чтобы они удовлетворяли одним и тем же дифференциальным уравнениям движения. Отсюда следует, что комплексные множители перед членами уравнения (5-6) одинаковы, т. е.

Заменяя здесь масштабные коэффициенты из (5-5), окончательно получаем с учетом уравнения энергии (5-3) следующие условия подобия:

1с 1'сct_c't', gl _,gl'.

"n T’ 1 T’ ~c*~ c⁷¹’

p _ p' V-^Cp_I^х-'с p

(5-7)

X' •

Таким образом, два потока являются динамически подобными, если выполняются соотношения (5-7) между параметрами этих потоков. Соотношения (5-7) называются критериями подобия. Первое соотношение (5-7) устанавливает равенство, чисел Рейнольдса в двух потоках:

Re = ^- = Re'.

Критерий Рейнольдса выражает соотношение между силами вязкости и силами инерции в потоке.

Второе условие является единственным, куда входят сходственные отрезки времени t и f, и оно получилось как следствие подобия членов, содержащих локальные уско

рения    в    уравнениях    движения.    Локальные

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

Отношение

Sh = ~ Sh'

носит название числа Струхаля; для периодического движения

где n — —_t--частота периодически нестационарного про-

цесса.

Третье уравнение дает равенство критериев подобия, учитывающих влияние массовых^- сил в потоке. Если ускорение в поле земного притяжения можно считать постоянным (g = g')_> то этот критерий, называемый критерием Фруда, легко представить в таком виде:

Fr=-i- = JV-

с²

В газовых потоках влияние массовых сил, как правило, невелико и поэтому критерий Fr при моделировании газовых потоков не принимается во внимание.

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

— 2р —

Р = -^=Р

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

Тогда для двух потоков получаем:

или

Ш² = &'М'².

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

Рассматривая дифференциальное уравнение энергии (5-3) для сжимаемой вязкой жидкости, можно получить дополнительный, уже знакомый из предыдущего критерий подобия Прандтля:

Из этого же уравнения вытекает также одинаковость кри териев k и М для газовых потоков.

При турбулентном течении вводится важная характеристика — степень турбулентности:

?= —

Q ⁶²

^ьт

I

гдес⁶³ = -д^ (с — с_т)² dt — средняя квадратичная скорость

пульсации; с — мгновенное значение скорости; с_м — сред-няя скорость турбулентного течения:

^Cm=-^r\^Cdt-

U

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

Анализ уравнений движения и размерностей величин, определяющих сопротивление обтекаемого тела или потери энергии в потоке газа, показал, что соответствующие безразмерные характеристики сопротивления являются функциями основных критериев подобия¹:

^ = ^(Ие_ЛМ, Рг, Sh, Е);

C = C(Re, М, Рг, Sh, Е),    (5-7а)

где с —коэффициент сопротивления (см. § 5-13);

С — коэффициент потерь энергии (см. §§ 4-5; 5-14;

8-5 и др.).

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

Так, например, для установившегося движения несжимаемой жидкости критерии М и Sh теряют смысл и зависимость (5-7а) упрощается.

С учетом сжимаемости при установившемся движении для газов с одинаковыми физическими свойствами с_х — = с, (Re, М).

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

Обычно осуществляется приближенное подобие (частичное моделирование) по одному или двум наиболее важным критериям. Результаты эксперимента нетрудно скорректировать и на другие критерии, если известны значения этих критериев в опыте и зависимость изучаемых характеристик от этих критериев.

При исследовании потоков, в которых определяющими являются два или три критерия (например, Re и М или Re, М и Sh), необходимо осуществлять принципы раздельного моделирования, т. е. обеспечивать возможность независимого изменения каждого из критериев в определенном диапазоне его значений.

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

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

\

5-3. ОДНОМЕРНОЕ ТЕЧЕНИЕ ГАЗА ПРИЛАЛИЧИИ ТРЕНИЯ.

ОСНОВНЫЕ УРАВНЕНИЯ

Основными уравнениями установившегося адиабатического движения вязкого газа являются уже известные нам уравнения неразрывности, количества движения и энергии.

Уравнение импульсов одномерного установившегося потока без энергетического обмена с внешней средой при наличии трения можно записать в таком виде:

cdc+^d-f- + dX_Tp = 0,    (5-8)

где dX_r —единичный импульс сил трения.

Величина dX_jp может быть выражена по гидравлической формуле: где С — коэффициент сопротивления; D — внутренний диа-метр трубы.

Решая совместно уравнения (5-8), (5-9), (2-14) и (2-6), можно получить дифференциальное уравнение распределения скоростей вдоль трубы переменного сечения с учетом влияния вязкости. Очевидно, что это уравнение аналогично (2-29), но должно содержать еще один член, учитывающий влияние вязкости. Простые выкладки приводят к такому выражению:

₍*²-i)x=(i-^²)^-^cdh^’ ⁽⁵'¹⁰⁾

где

Формулу, определяющую изменение давления потока вдоль трубки переменного сечения, следуя уже известному методу, изложенному'в гл. 2, можно получить в следующем виде:

( Ап* \

dp    2k h? [dF * + ! r 1 /к 11\

¦7—rpv=i[-r-_T=gr_rS*')-

Из уравнений (5-10) и (5-11) можно заключить, что изменение параметров течения в трубкб^переменного сечения происходит под воздействием двух факторов: деформации потока (изменение сечения трубки) и сил трения. Уравнения показывают, что влияние трения всегда является односторонним. Так, например, при дозвуковых скоростях (1< 1) в суживающейся трубке (dF <^0) трение способствует ускорению течения (dl^>0 и dp<C_ 0). При сверхзвуковых скоростях в такой же трубке (dF< 0) трение приводит к замедленному падению скорости и соответственно к более медленному возрастанию давления по сравнению с идеальным процессом без потерь.

Из уравнений (5-10) и (5-11) следует, что в простейшем случае трубки постоянного сечения (dF — 0) при Я<1

имеем ^->0 и^-<0 и, следовательно, поток ускоряется. Нетрудно видеть, что в этом случае при сверхзвуковой скорости (Я^>1)^-<;0: поток тормозится.

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

Следовательно, воздействие сил трения в потоке можно заменить эквивалентной деформацией струи — уменьшением ее сечения в направлении движения.

Правомочность такой замены вытекает из следующих рассуждений. Движение газа в трубе без энергетического обмена, но при наличии сил трения является необратимым адиабатическим процессом. Такой процесс, как нам уже известно, сопровождается ростом энтропии. Приращение энтропии можно выразить по уравнению (4-27).

получим:

Для энергетически изолированно]    ея    в    виду,    что

7’oi = 7'o2 и что уравнение состояния

(5-12)

Так как энтропия в направлении потока возрастает, давление торможения /?₀2 должно уменьшаться.

По уравнению неразрывности (2-41) для двух сечений трубы заметим, что при одинаковых в обоих случаях величинах F_v А, и Л₂ сечение F₂ при наличии сил трения должно быть больше, чем для изоэнтропического потока.

Величину s₀F можно рассматривать как живое сечение струи. Для трубы постоянного сечения (F — const) очевидно, что живое сеченне e_aF при наличии трения уменьшается, так как уменьшается величина е₀.

Уравнение (5-10) может быть использовано для некоторых заключений о положении сечения, отвечающего экстремальным значениям скорости. Выразим из уравнения (5-10) Логарифмическую производную сечения:

d\ k -^-^)-г + ^к + т^Шх

dF

F

отсюда следует, что при dX = 0 и X ф-1

k

Следовательно, сечение, соответствующее максимальной скорости при 1<1 и минимальной скорости при X >• 1, не совпадает с минимальным сечением, а смешается в расширяющуюся часть трубки. Соответственно критическим условиям (Я=1) также отвечает сечение в расходящейся части трубки. Это означает, что минимальное и критическое сечения трубки при наличии трения не совпадают. В этом случае, полагая в (5-10) 1 = 1, получим:

dF ft

-у—~2~ Ыху-0.

5-4. ДВИЖЕНИЕ ГАЗА В ЦИЛИНДРИЧЕСКОЙ ТРУБЕ

, Из уравнения (5-10) для цилиндрической трубы (^-у = oj получим:

Из предыдущих рассуждений, а также непосредственно из анализа уравнения (5-14) легко прийти к выводу, что критическая скорость течения может возникнуть только в выходном сечении цилиндрической трубы. Действительно, согласно уравнению (5-14) при Х<^1 и dX^> 0 поток в трубе ускоряется, а при X > 1 и dX <f 0 он замедляется. Случай X = 1 в промежуточном сечении трубы противоречит уравнению (5-14) и физически нереален.

Допустим, что коэффициент сопротивления является величиной постоянной⁶⁴. Тогда уравнение (5-14) можно проинтегрировать.

Интеграл уравнения (5-14) запишем в такой форме:

^{1 1} 1 *-⁶⁵ 2ft , х2 К²    ^2    А + 1    ’

где Х_г — безразмерная скорость в начальном сечении трубы;

X — безразмерная скорость в некотором сечении на расстоянии х от начального.

Введем безразмерную координату, которую называют приведенной длиной трубы:

2й

U.

k+ 1

Тогда уравнение (5-14а) можно записать так:

X²

1

X*

— In

/ =

(5-15)

Зависимость между Я и % при постоянном значении Я₁ ¦приведена на рис. 5-3. Величина % имеет максимум при Я = Я₂ = 1.

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

= * _1₊1_ПЯ?.

(5-16)

Кривые х(Я) состоят из двух ветвей, отвечающих дозвуковому (Я<^1) и сверхзвуковому (Я^>1) потокам в трубе постоянного сечения (рис. 5-3).

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

Рис. 5-3. Зависимость безразмерной скорости на выходе из трубы от скорости на входе и приведенной длины трубы.

в другую. В такой трубе, как было показано выше, при определенной скорости на входе Я₁ и соответствующей длине на выходе достигается критическая скорость (Я₂=1).

Дозвуковым потокам на входе в трубу (Я, < 1) отвечает участок кривой АВ (рис. 5-3)_А а сверхзвуковым (Я_ж>1) — участок СВ. Точка В определяет максимальную величину функции % для данного значения Я_1# входе в трубу, а верхняя — сверхзвуковым.

Таким образом, уравнение (5-16)' показывает, что для цилиндрической трубы заданных размеров I и D при скорости на выходе из трубы Я_а = 1 и для определенных значений k и С безразмерная скорость на входе в трубу Я₁₅ а вместе с тем и приведенный расход газа q_x имеют строго определенные значения.

Из формулы (5-16) сЛедует, что при ^=1 у_макс = 0. Графически зависимость (5-16) представлена на рис. 5-4. Кривая также имеет две ветви. Нижняя ветвь о.твечает дозвуковым скоростям на

3 2Л 2 1,5 1

10

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

в

При дозвуковой скорости на входе в цилиндрическую трубу длиной I, характеризуемую коэффициентом сопротивления С, при установившемся движении может пройти максимальное количество газа, если 1₂ = \.

Абсолютный расход газа через трубу предельной длины будет равен:

<Wc = &^F    =    §^F(f,    «акср,    А    =

= sF (4-!/“' -г,    (1-^1    Р.Л

Вспоминая, что

1

k—I

V

2k -k + l^o ^—

Р А = Р.

k +1

(k + \)R У r₀₁

у

1

k— 1

Poi_

Vt7'

(5-17)

V HI ¹ макс/

получим:

G =F V~^^k

макс    f

(k + \)R    1    макс

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

При сверхзвуковых скоростях на входе в трубу, как показывает опыт, обнаруживаются некоторые новые свойства потока, которые не описываются уравнением (5-15). Заметим, что согласно уравнению (5-15) при ^ > 1 скорость в трубе должна непрерывно падать к выходному сечению по кривой СВ на рис. 5-3, а давление — соответственно непрерывно возрастать. Однако в действительности изменение скоростей и давлений в трубе в ряде случаев происходит скачкообразно.

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

Так как в изолированной трубе г₀ = const, то для любых двух сечений можно написать Т₀₁ = Т₀₂ = const. Из этого услсвия получаем уравнение для TjT_a в форме уравнения (2-22).

Для отношений давлений можно воспользоваться формулами (2-41 а) и (2-42). После несложных преобразований получим связь между статическим и полным давлениями в следующем виде:

(5-18)

Отсюда при X = Я₂ = 1 определяется критическое отношение давлений:

k

(5-19)

Формула (5-19) показывает, что критическое отношение

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

= / 2

k +1

ft—1

Уравнения (5-18) и (5-15) позволяют построить графики изменения давлений вдоль трубы для заданных значений ?i и X-

Подобный график приведен на рис. 5-5 для случая сверхзвуковой скорости на входе в трубу ^ = 1,76 и q_x= = 0,453. Здесь кривая АВ характеризует повышение давления в трубе до критического значения в точке В, равного:

s'= 0,528-0,453 = 0,239.

*

Если известно распределение скоростей по трубе, а оно легко подсчитывается по уравнению (5-14а), то молено по формулам (4-20) и (4-24) определить скорости и давления за прямым скачком уплотнения в каждом данном сечении (линия СВ).

За прямым скачком поток — дозвуковой и, следовательно, давление в нем под влиянием сил трения должно падать. Так, если прямой скачок возникает непосредственно во входном сечении, то дальнейшее изменение давления протекает по кривой CD. Характер изменения давлений в дозвуковом участке трубы при различных промежу-

0,2

0.1

ОА

0.3

0.1    02    0,3    ОА    0.5    0,6    0,7    Ofi

Рис. 5-5. Распределение давлений по трубе постоянного сечения.

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

При указанной скорости на входе 1_г и приведенном расходе газа q_x режимы в трубе без скачков возможны в тех случаях, когда у < У_макс> причем максимальное значение приведенной длины отвечает точке В. При условии макс ^{скачки в} трубе возникают только тогда, когда давление на выходе из трубы больше соответствующего давления, показываемого кривой АВ.

Допустим, что труба имеет длину, определяемую точкой / (х = 0,35), а давление за выходным сечением задано точкой L, которая лежит на кривой СВ. В этом случае прямой скачок располагается в выходном сечении трубы IL. Если давление среды, куда вытекает газ из трубы, определяется точкой К, то прямой скачок перемещается внутрь трубы и располагается в сечении EF, причем отрезок FK соответствует понижению давления в дозвуковом участке трубы. Последующее повышение давления среды приводит к дальнейшему перемещению скачка внутрь трубы (к входному сечению).

Если давление среды определяется точкой S, то в трубе скачков не возникает, а в струе, выходящей из трубы, образуется конический скачок (или система плоских косых скачков, если труба прямоугольного сечения). При уменьшении давления до величины давления в точке / интенсивность конического скачка уменьшается. В точке / конический скачок вырождается в слабую коническую волну, при пересечении которой энтропия не меняется⁶⁶.

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

Кривые давлений в трубе (рис. 5-5) показывают, что при постоянной длине у ^и давлении на выходе р_г с увеличением скорости на входе прямой скачок смещается к выходному сечению. При увеличении сопротивления трубы (путем, например, подключения дополнительного участка трубы) перемещение скачка происходит в обратном направлении (ко входу в трубу).

5-5. ПОТЕРИ НА ТРЕНИЕ В ЦИЛИНДРИЧЕСКОЙ ТРУБЕ (ОПЫТНЫЕ ДАННЫЕ)

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

В действительности коэффициент сопротивления зависит от числа Рейнольдса и—в общем случае — от безразмерной скорости М. Такая зависимость устанавливается экспериментально.

Число Рейнольдса для произвольного сечения цилиндрической трубы определяется по формуле

Re =???-. г-

Для цилиндрической трубы

1

/ 2 \ ¹ / 2k P^c==P₁^ci⁼^₁P_M^a, — 4i    V    F+T^oiPoi’

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

* =(гЬГ /iTlпг ' ^<5'²⁰⁾.

Из формулы (5-20) видно, что число Re меняется вдоль трубы только вследствие изменения коэффициента вязкости (а, который зависит от температуры. Можно показать, однако, что в теплоизолированной трубе изменения температуры невелики. Так, при изменении скорости водяного пара в трубе от ^ = 0,2 доЯ₂ = 1 температура изменяется на И°/о> в то время как давление уменьшается в 4,5 раза, а плотность — в 5 раз.

Изменение коэффициента вязкости воздуха в зависимости от температуры можно оценить по приближенной формуле

(этз) ’    (5-21)

где а — постоянная величина; р₀ — коэффициент вязкости при температуре Т — 273°.

По опытным данным для воздуха

ц-10^е = 1,757-f 0,00483*° С.

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

деления всех членов на а² принимает вид;

(5-22)

Замечая, что для цилиндрической трубы произведение

Vlal=q?al = qe_t кр_й

является постоянной величиной, имея в, виду (2-46), запишем уравнение (5-22) в такой форме:

dl -J- die -J- ? ~dx = 0.

Здесь функция и зависит только от X по формуле

(2-46a).

После подстановки этой функции получим:

^ "“Ь ? ~2~dx = 0.

Отсюда можно выразить коэффициент сопротивления:

(5-23)

Уравнение (5-23) является исходным для экспериментального определения местных значений С. Перейдя в этом уравнении к конечным разностям, будем иметь:

(5-24)

Если известны расход газа G, температура торможения Т₀ и распределение давлений вдоль трубы [р = р (х)], то по формуле (5-24) можно найти средний коэффициент сопротивления для некоторого небольшого участка трубы длиной Ах. Если участки Ах выбираются небольшими, то определенное таким путем значение С не сильно отличается от его местного значения.

Для нахождения связи между р и Я можно воспользоваться формулой (5-18).

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

⁽⁵-²⁵⁾

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

На рис. 5-6 показано изменение профиля скоростей в трубе в зависимости от числа Re для турбулентного режима. Здесь же нанесена эпюра скоростей для ламинарного потока в трубе. Сравнение кривых показывает, что профиль скоростей при турбулентном режиме значительно более „полный”, чем при ламинарном, причем с ростом числа Re наполнение профиля увеличивается.

'Рис. 5-6. Распределение скоростей при ламинарном и турбулентном режимах в трубе.

Как известно, при турбулентном двил ении, возникающем в трубах при Re > 3 ООО, имеют место перемещения макрочастиц⁶⁷ в поперечном направлении. При этом частицы внешнего потока, обладающие большой кинетической энергией, перемещаясь к поверхности, увеличивают кинетическую энергию частиц у стенки, движущихся с малыми скоростями, и наоборот, частицы, переместившиеся '-от стенки в ядро потока, тормозят здесь движение жидкости.

В соответствии с изменением профиля скорости в зависимости от числа Re коэффициент сопротивления трубы должен также меняться в функции этого параметра.

Для оценки коэффициента сопротивления цилиндрических труб при небольших скоростях можно пользоваться кривыми Всесоюзного теплотехнического института, построенными Г. А. Муриным (рис. 5-7). Здесь коэффициент С представлен в зависимости от числа Re и величины, обратно пропорциональной относительной шероховатости Djk₉, где k₃ — средняя высота бугорков шероховатости. Заметим, что при больших значениях Djk₃ (малая шероховатость) значения С по кривым ВТИ удовлетворительно совпадают с формулой Никурадзе:

^ = ~0,8 + 21g(Re]/g.    (5-26)

Кривые ВТИ и формула (5-26) отчетливо показывают, что влияние числа Рейнольдса на коэффициент сопротивления в гладкой трубе простирается до весьма больших значений Re 5= 10^е. По мере увеличения шероховатости влияние Re (при Re 2 -10⁵) на увеличение шероховатости ослабевает. На рис. 5-7 пунктирной линией соединены точки, соответствующие тем значениям Re, выше которых влияние этого параметра практически не обнаруживается. Справа от этой линии расположена область, которую принято называть автомодельной⁶⁸.

Перейдем теперь к рассмотрению влияния второго основного критерия подобия — числа М — на коэффициент сопротивления в трубах. Соответствующие опытные данные

Рис. 5-7. Зависимость коэффициента сопротивления S стальных труб от числа Re и относительной шероховатости по данным ВТИ.

получены в МО ЦКТИ\ в МЭИ⁶⁹ и ряде других организаций.

На рис. 5-8 показано распределение статических давлений по длине гладкой трубы при дозвуковых скоростях, отвечающее различным числам Re_t и на входе в трубу. В соответствии с методикой проведения опытов кривые на рис. 5-8 показывают совместное влияние числа Rej и отношения давлений (числа М,), так как раздельное моделирование по Re и М не осуществлялось. По мере снижения s

Рис. 5-8. Распределение давлений, температур и скоростей по длине трубы для дозвуковых скоростей.

(увеличения ReJ характер кривых давлений меняется и особенно интенсивно на выходном участке (х > 60 -н 70). Переход к большим значениям Re и М_г сопровождается увеличением градиентов давления: наклон линий р возрастает. Для группы режимов, отвечающих сверхкритиче-скому перепаду давлений в трубе, статическое давление в выходном сечении превышает давление внешней среды, но оказывается ниже той величины s' , которая, соответ-

*

ствует критическому истечению. Значение г можно найтц по формулам (5-19) или (5-25), подставив 1= 1.

Величина г отмечена на рис. 5-8 пунктиром. Отсюда следует, что критическое сечение не совпадает с выходным

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

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

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

Рис. 5-9 Изменение статических давлений и скоростей по диаметру трубы вблизи выходного сечения.

вания поля скоростей и давлений на выходном участке приведены на рис. 5-9. Здесь отчетливо видна неравномерность в распределении статических давлений по диаметру трубы, причем давление на оси во всех сечениях за критическим выше, чем у стенки.

. Эпюры скоростей -(рис. 5-9) позволяют заключить, что толщина дозвукового пристеночного слоя на выходном участке уменьшается в направлении потока. Можно предполагать, что такая структура потока объясняется взаимодействием вытекающей струи с внешней средой. Благодаря интенсивному отсосу пристеночного слоя во внешнюю среду происходит его утонение на выходном участке (рис. 5-9). При этом в ядре потока создаются условия, необходимые для перехода к сверхзвуковым скоростям:

сечение ядра увеличивается по течению. Давление внешней среды „проникает” через дозвуковую часть пристеночного слоя внутрь выходного участка, и давление на стенке оказывается более низким, чем давление на осп.

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

в зависимости

от М по данным МО ЦКТИ и МЭИ. Зависимость С от Re

— ?

На рис. 5-10 нанесены значения С =4-

0,7

0,9

КО

Рис. 5-10. Зависимость коэффициента сопротивления от числа М при дозвуковых скоростях по данным МО ЦКТИ и МЭИ.

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

График на рис. 5-10 показывает, что при М<< 0,70-й),75 коэффициент сопротивления не зависит от М и весьма близок к С . Для этой области расчет С может производиться по любой эмпирической формуле [например, (5-26)] или по кривым ВТИ*.

В диапазоне чисел М = 0н-0,8 и Re — 3-10³ -г- 3-10⁵ удовлетворительно совпадает с опытом формула, полученная в МЭИ:

C=0,0334Re7°'²,    (5-27)

* В последнем случае необходимо убедиться в том, что для шероховатых труб влияние сжимаемости на S при М < 0,7 ч- 0,8 будет также невелико. Это предположение нуждается в опытной проверке.

где Re =Re;t — число Re, отнесенное к длине трубы.

/,г ко 0,8 0,6 о,ь Q2

оА

о,з

0,8

0,8

При М > 0,7 н- 0,75 коэффициент сопротивления трубы уменьшается с ростом М; особенно интенсивное снижение С наблюдается при скоростях М > 0,85.

Напомним, что изменение давления на элементарном участке трубы dx выражается известной гидравлической формулой:

dp = ^^F-^-dx.

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

71D² ,    J-.J    71D² _r PC² ,

— dp = *₀TiDdx = -j- С dx.

Отсюда можно получить формулу, связывающую напряжение трения на стенке и С:

В соответствии с формулой (5-28) результаты изложенного выше опытного исследования можно выразить так: влияние сжимаемости, заметно проявляющееся при ЛГ>0,75, приводит к некоторому уменьшению силы трения, отнесенной к кинетической энергии потока в данном сечении.

Физически этот результат объясняется тем, что с ростом числа М увеличиваются градиенты давления в трубе (рис. 5-8). Возрастание градиентов давления в конфузорном течении вызывает деформацию профиля скорости у стенки; наполнение профиля скорости увеличивается. Кроме того, пристеночный слой при этом утоняется.

Околозвуковой поток особенно чувствителен к изменению сечения, что видно из уравнения (5-8). Поэтому на концевом участке трубы, где М. 0,9, наблюдаются весьма большие отрицательные градиенты давления и соответственно резкое снижение С.

Резкое уменьшение ? при М > 0,9 связано также с тем, что область скоростей М = 0,9-^1,0 расположена вблизи концевого участка трубы, где пристеночный слой нарушается. При расчете по формуле (5-28) значительная деформация профиля скоростей на выходном участке трубы не Гучитывается.

Оценивая влияние сжимаемости на коэффициент сопротивления трубы при сверхзвуковых скоростях, необходимо различать три основных режима течения в трубе. Первый режим отвечает бесскачковому движению потока, скорости которого в каждом сечении трубы сверхзвуковые. Как уже было показано, такой режим возможен, если длина цилиндрической трубы меньше предельного значения ('/ -< '/_макс). Если же в трубе имеется соответствующий источник возмущения, ТО при Х<СХ_чакс сверхзвуковой поток может быть насыщен скачками уплотнения. Возмущение потока в трубе может быть вызвано угловым поворотом стенки, который образуется в сечении стыка сопла Лаваля с трубой. В простейшем случае расширяющаяся часть сопла выполняется конической с различными углами раствора. Чем больше угол раствора сопла, тем больше угол отклонения потока при входе в трубу и тем интенсивнее скачок, образующийся в точке поворота стенки. Такие режимы течения с коническими скачками, когда поток вплоть до выходного сечения остается сверхзвуковым, составляют вторую группу режимов. Наконец, если Х>Х_Макс>^Т0 внутри трубы возникает сложная система скачков, за которой поток будет дозвуковым; это — третья группа режимов течения. Распределение параметров течения по длине трубы в тех указанных случаях оказывается существенно различным.

Для третьей группы режимов, когда в результате скачков поток становится дозвуковым, характер распределения давления по данным Неймана и Лустверка показан на рис. 5-11 (для Х₁ = 2,2). По мере повышения давления в выходном сечении трубы система скачков перемещается к соплу. Интенсивность скачков при этом повышается. Отметим, что протяженность системы скачков оказывается значительной. За системой скачков поток дозвуковой, и давление вдоль трубы падает.

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

На рис. 5-12 нанесены значения коэффициента сопротивления в зависимости от М по данным МО ЦКТИ для интервала чисел М = 0-г-1,65. Здесь воспроизведены опытные точки С для дозвуковых скоростей, частично представленные на рис. 5-8, и добавлены результаты более поздних исследований для М. > 1. Характерным для области небольших сверхзвуковых скоростей следует считать заметное возрастание здесь коэффициент сопротивления меняется от 0,007 до 0,018.

Можно предположить, что на этом участке происходит турбулизация пристеночного слоя в трубе, т. е. переход его из ламинарного режима в турбулентный⁷⁰. В предыдущей области М = 0,95-г-1,03, где градиенты давления достигают максимальных значений, по-видимому, происходит ламинаризация пристеночного слоя, так как при больших отрицательных градиентах давления интенсивность турбулентных пульсаций снижается.

МО ЦКТИ.

На участке М = 1,25 -г- 1,6 С снижается и при М = 1,65 достигает минимального значения С ~ 0,01. Это снижение объясняется наполнением профиля скорости в области сверхзвукового течения с отрицательными градиентами давления.

При М>1,3 влияние числа Рейнольдса на С невелико.

5-6. ПОГРАНИЧНЫЙ СЛОЙ. ОСНОВНЫЕ ПОНЯТИЯ И УРАВНЕНИЯ¹

Современные представления о механизме сопротивления тел, обтекаемых потоком газа, и методы расчета сопротивления основываются на теории пограничного слоя².

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

Характер распределения скоростей в пограничном слое на плоской стенке показан на рис. 5-13. Непосредственно на стенке жидкость „прилипает“ к поверхности (Я₅ = 0). В тонком пограничном слое с удалением от поверхности тела скорости меняются весьма интенсивно и на небольшом расстоянии от поверхности достигают скорости внешнего потока.

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

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

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

Рассмотрим более подробно условия образования пограничного слоя на поверхности крылового профиля (рис. 5-13).

В направлении потока вдоль поверхности толщина пограничного слоя 8 увеличивается. Заметим, что понятия внешней границы и толщины пограничного слоя не являются определенными, так как изменение продольных скоростей при переходе из пограничного слоя во внешний поток происходит плавно. Приближенно внешняя граница слоя определяется в тех точках, где продольная скорость отличается от скорости внешнего потока на малую величину порядка 1°/₀.

Внешняя граница пограничного слоя не совпадает с линией тока, так как частицы внешнего потока непрерывно проникают в пограничный слой (рис. 5-13).

Передний участок слоя, расположенный вблизи носика профиля, обычно является ламинарным. На небольшом расстоянии от носика (сечение II на рис. 5-13,а и б) обнаруживается ламинарный профиль скоростей. В некотором сечении III начинается переход ламинарного упорядоченного движения в турбулентное, которое носит пульсационный характер.

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

За переходной областью развивается устойчивый турбулентный слой. Здесь область ламинарного подслоя настолько мала, что экспериментально ее трудно обнаружить. Как видно из рис. 5-13,6, турбулентный слой имеет более полный профиль скоростей.

Схема образования пограничного слоя на рис. 5-13 изображена не в одинаковом масштабе по осям л: и у. Толщина слоя 8 весьма невелика по сравнению с размерами тела и составляет сотые ц тысячные доли хорды профиля.

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

С этой целью перейдем в уравнениях (5-4) к безразмерным величинам. Для простоты рассмотрим установившееся движение несжимаемой жидкости. В качестве мае-

Рис. 5-13. Схема образования пограничного слоя на поверхности крылового профиля, в—изменение толщины и структуры слоя вдоль поверхности: 1—ламинарный участок слоя; 2—переходная область; «? — турбулентный участок слоя; 4—ламинарный подслой; 5 —переходный слой; б —распределение скоростей в различных

сечениях слоя.

штабов продольных скоростей и координат выберем некоторую характерную скорость и₀ и характерный линейный размер L₀. Масштабы поперечных скоростей и координат обозначим соответственно у₀ и 8. Обозначим масштабы для давлениями плотности р₀ и р₀. Тогда после преобразований, аналогичных преобразованиям § 5-2, получим следующую систему    уравнений:

ди |    ?°и° ди__р₀р    др, у    д²и    ,    vZ,₀    д²и.

^U    5и₀^V ду    _рм2    йхТ    1₀к₀ дх²'8²и₀    ду² ’

дл:    8и₀    ду

Масштабы величин, входящих в систему уравнений (5-29), можно выбрать на основании следующих соображений. Если масштабы и₀ и L₀ выражают характерную скорость и линейный размер обтекаемого тела, то величины v₀, 8, p₀ и p₀ пока остались неопределенными. Пользуясь этим произволом, уравнения (5-29) нетрудно привести к каноническому виду.

Действительно, полагая /V=P^Mo , Ро— 1 полу_ЧИм ^р-^°= 1, а

Р“о

поперечные масштабы у₀ и 8 выберем таким образом, чтобы коэффициенты ^и были постоянными и не зависели от числа Re.

Положим

Lv

-•=1 и

1.

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

Решая написанную систему относительно »₀ и 8 для масштабов поперечных скоростей нелинейных размеров, получаем следующие значения

(5-30)

V Re

При этом уравнения (5-29) примут вид:

ди

' дх

ди  др .

ду    cto^Re дх

1 д²и г д²и

_1_

Re

' da , dv \    др    |

+

(5-31)

du.dv дх ' ду

1 d²v Re² дх²

0.

ду²'

1 d²v. Re ду²'

Если допустить, что неизвестные величины и, v, р и их производные с увеличением числа Рейнольдса стремятся к определенным пределам в фиксированных точках, то при больших Re в уравнениях (5-31) можно отбросить все чле-

1 1

ны, имеющие множители ^ и , как малые величины по сравнению с другими членами.

rf результате перехода вновь к размерным величинам получим дифференциальные уравнения ламинарного пограничного слоя в форме Л. Прандтля:

du    I    ida    1    др    ,    д²и.

U г--f“ V -д— —    г--j— V -г 2 >

дх    ¹    ду    р    дх    ¹    ду²

ди I ди

дх'ду

(5-32)

Система (5-32) должча быть решена при следующих граничных условиях:

у = 0; ц = 0; у = 0;

у оо; и—* и (х).

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

Полученное условие    означает, что распреде

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

Условие ~ = 0 позволило объяснить весьма важное ^дУ

явление отрыва пограничного слоя. Рассмотрим обтекание некоторой криволинейной поверхности АВ (рис. 5-14), предполагая, что давление внешнего потока вдоль этой поверхности вначале уменьшается, достигает минимального значения в точке М и затем увеличивается. Участок внешнего потока, в котором градиенты давления отрицательны

(^<0\ называется конфузорным участком. Область течения за точкой Af, характеризуемая положительными градиентами давления    ,    называют    диф-

фузорным участком. На конфузорном участке внешний поток ускоряется, а на диффузорном — тормозится.

Учитывая, что в пограничном слое |^ = 0, заключаем, что

совершенно аналогичное распределение давлений имеет место и вдоль поверхности АВ на любом расстоянии у<СР в пограничном слое.

Рис 5-14. Схема образования отрыва пограничного слоя.

В пределах пограничного слоя скорости перед точкой М увеличиваются, а за нею — уменьшаются (см. эпюры скоростей на рис. 5-14). Частицы жидкости вблизи стенки обладают малой кинетической энергией, причем в диффу-зорной области вдоль поверхности АВ запас кинетической энергии частиц уменьшается. В результате в некотором сечении S частицы у стенки не могут преодолеть тормозящего влияния внешнего потока и останавливаются. Эпюра скоростей принимает характерную остроконечную форму. На стенке кривая скоростей удовлетворает условию

' (5-32а)

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

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

Используя уравнения (5-32), легко показать, что положение точки отрыва ламинарного пограничного слоя не зависит от числа Re. Действительно, решение системы (5-32), дает:

u = f(x,y),    (5-326)

где х и у — безразмерные координаты.

Тогда, продифференцировав по у и используя в точке отрыва условие (5-32а), получим;

f'(x_s, 0) = 0.

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

5-7. УСЛОВНЫЕ ТОЛЩИНЫ И ИНТЕГРАЛЬНОЕ СООТНОШЕНИЕ ДЛЯ ПОГРАНИЧНОГО СЛОЯ

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

Поэтому в расчетах пограничного слоя вводятся другие интегральные толщины, зависящие от S: толщина вытеснения 8*, толщина потери импульса 8’* и толщина потери энергии 8*".

Для выяснения физического смысла указанных толщин сравним течение идеальной и вязкой жидкостей около твердой стенки (рис. 5-15).

При отсутствии трения за единицу времени через поперечное сечение потока высотой dy и шириной, равной единице, протечет масса р_аи_йёу. В пограничном слое за то же время через сечение dy протечет масса рudy.

Разность^этих количеств составит:

ОО

^Л“*10 —

о

=^РА\ (’" ьУ    К¹ ^ )^dy-

О    5

Рис. 5-15. К определению условных толщин пограничного слоя (а). К выводу уравнения импульсов для пограничного слоя (б).

е)

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

Разделив найденный излишек массы на р₀и₀, получим:

Величина 8* показывает смещение линии тока в направлении внешней нормали к контуру обтекаемого тела.

Вместе с тем 6* характеризует уменьшение расхода жидкости через сечение слоя, „нормальное* к стенке, обусловленное „вытеснением” жидкости пограничным слоем, и поэтому носит название толщины вытеснения.

Толщина потери импульса 6** равна такой толщине слоя жидкости, движущейся со скоростью и₀ вне пограничного слоя, количество движения которой равно импульсу сил трения в пограничном слое. Это количество движения, „потерянное" в пограничном слое, будет равно:

ь    г

ри (й₀ — и) dy = Г рии₀ ГI — ~ j dy.

о

о

Разделим полученное выражение на р_аи_й . Тогда полу

чим:

О

(5-34)

Масса жидкости рudy теряет в пограничном слое кине-

тическую энергию, равную ?и. (и₀ — и?) dy. Для всего слоя эта потеря составит:

5

Тогда толщина потери энергии

(5-35)

представляет собой толщину движущейся вне слоя жидкости, обладающей кинетической энергией, потерянной в пограничном слое.

Для решения задач о течении сжимаемой жидкости в ряде случаев целесообразно иметь одинаковую структуру формул, определяющих интегральные толщины. Поэтому наряду с формулой (5-33) толщину вытеснения §⁷¹ часто подсчитывают по формуле

(5-36)

S

о

При этом, естественно, нарушается приведенный ранее физический смысл толщины вытеснения. Для несжимаемой жидкости уравнения (5-33) и (5-36) оказываются тождественными.

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

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

Из определения толщины вытесненйя следует, что всю массу жидкости, протекающей в пограничном слое, можно условно заменить массой, расположенной между линией ab и стенкой ей (рис. 5-15), скорость которой равна нулю („вытесняемая" масса), и массой, протекающей выше ab

со скоростью и .

о

Со стороны стенки на вытесняемую массу действуют силы трения, а в направлении течения будут приложены силы давления. Скорости выше линии вытеснения ab равны и₀ — и_а(х) и в пределах рассматриваемого объема на основании дифференциальных уравнений пограничного слоя

(5-32) ~ = 0. Применяя уравнение импульсов, можно найти

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

(5-37)

dl = x₀dx~\- dpb*,

где t_adx и dpb* — секундные импульсы сил трения и сил давления, действующих на „вытесняемую" массу жидкости.

Количество движения I на основании уравнения (5-34) можно выразить через толщину потери импульса:

/= Ро«о⁸**-

Поскольку на внешней границе пограничного слоя течение считается потенциальным, продольный градиент давления ^ легко выражается на основании уравнения Бернулли через скорости и₀ и плотность р₀:

dp    du„    ’

dx ~~ ~ ^Р°^йо 5F ~ Ро^Мо«₀-

Подставляя эту величину в (5-37), находим:

р(5-37а)

Заменим далее

_ ^Ро

dx~dp'dx~ _a\ ^{Po 0} »

и разделим (5-37a) на р₀и*. В результате получим для сжимаемой жидкости следующее уравнение:

ij+if-(2 + tf-M₀) = -^V.    (5-38)

“0    Ро"о

Здесь

Уравнение (5-38) носит название интегрального соотношения, так как величины 8** и 5* выражаются интегралами (5-33) и (5-34).

Интегральное соотношение (5-38) для пограничного слоя можно получить, не прибегая к понятию толщины вытеснения. С этой целью уравнение импульсов применяется к объему жидкости, заключенному между двумя бесконечно близкими поперечными сечениями пограничного слоя

(рис. 5-15,6). Подставляя в (5-38) выражения для условных толщин 8*, 5** и заменяя и^ через ^, после преобразований находим:

щйу = — т₀— 8^-.    (5-39)

JL

dx

j рu²dy

о

о

Для несжимаемой жидкости получим:

г

г

о

о

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

| При возникновении отрыва уравнение импульсов может служить для определения положения точки отрыва, в кото-

5-8. ОБЩЕЕ ВЫРАЖЕНИЕ ДЛЯ КОЭФФИЦИЕНТА СОПРОТИВЛЕНИЯ ТРЬНИЯ В ПОГРАНИЧНОМ СЛОЕ ПРИ НАЛИЧИИ ГРАДИЕНТА ДАВЛЕНИЯ

В уравнении импульсов содержатся две искомые переменные величины: толщина слоя 8 (или взаимосвязанные условные толщины 8* и 8**) и напряжение трения на стенке т₀.

В общем случае т₀ определяется скоростью на внешней

fff

U

границе пограничного слоя, ее производными    ,    и₀    и

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

Используя основные положения теории размерности, из функциональной зависимости (5-40) нетрудно получить структурную формулу для коэффициента сопротивления.

Примем в качестве основных размерности скорости й₀, плотности р, длины §** и температуры Т. Простой проверкой легко убедиться, что, комбинируя указанные величины, можно получить размерности всех остальных параметров. Действительно, т₀ имеет размерность л:г/ж². Ту же размерность будет иметь и комплекс рм₀: [кг-сек²/м*-м²/сек²]. Следовательно, отношение ^₀/р₀Мд, представляющее собой местный коэффициент трения c_fi окажется безразмерным. Переходя в выражении (5-40) от размерных величин к безразмерным, получим

¦с.    (    и'_пд**    и’’ 8**²    \

^c/=^=4^Re**^;M; -V^{:    (5}-⁴¹⁾

и §**

Здесь Re**==    --число Рейнольдса, подсчитанное

по толщине потери импульса.

Число безразмерных параметров в выражении (5-41) можно сократить, если принять, что напряжение трения определяется так же, как и в случае ламинарного течения только первой производной скорости м_о-

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

Далее, можно показать, что при отнесении физических постоянных v и р к температурным условиям на стенке число М исключается из зависимости (5-41). Таким образом, с учетом принятых допущений

^(М2)

? ?    tf ф ¦ #

Где Re^ = ——, ч_т— кинематическая вязкость, подсчитан-

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

О

* * \

+

Р®“о

^(RO + 'MRC)

+<mrC)

Здесь c_f—местный коэффициент трения, подсчитанный по плотности у стенки р_да.

Выражение (5-43) является общим как для ламинарного, так и для турбулентного режимов течения в пограничном слое. В зависимости от режима течения коэффициенты <|>₀, 4*!, ф₂ и т. д. будут принимать различные значения.

Найдем конкретный вид выражения (5-43) для ламинарного течения. С этой целью запишем (5-43) в таком виде:

h.

Фо

(5-44)

1 +

+•

'f

—V ^{= l}h

При «_о = 0 формула (5-44) должна совпадать с соответствующей формулой для коэффициента сопротивления плоской пластинки при безградиентном ее обтекании.

В этом случае задача решается достаточно точно путем численного интегрирования системы (5-32) и независимо от метода решения для безградиентного течения коэффициент сопротивления c_f0 выражается формулой

До

^0 = Ф»^:

Re*

где а₀ — постоянная величина. Следовательно,

ф, м'б

Re*

В точке отрыва t₀ = 0; в этом случае выражение в квадратных скобках должно обращаться в нуль. Учитывая, что положение точки отрыва не зависит от числа Re**, получим: Ф1⁼ const; ф_в = const; ф₃ = const и т. д. Отсюда

и₀8**^г

Re:

Коэффициенты а₀, a_t, а₂ и т. д. в общем случае определяются экспериментально. Однако для ламинарного пограничного слоя их можно определить и теоретически. Так, например, А. М. Басин получил а₀ = 0,22; а₁=1,85; а₂ = — 7,35.

Обозначим

(5-46)

и

(5-47)

?(f) — K + a₁f + a»/²+...];

тогда

(5-48)

Параметр f часто называют формпараметром; как будет показано ниже, он определяет форму профиля скорости в ламинарном пограничном слое. Следует отметить, что структура формпараметра, содержащего производную

du„

, отражает влияние продольного градиента давления

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

Для турбулентного слоя многочисленные эксперимен-

(5-49)

и для коэффициента сопротивления получается из (5-44) следующая формула:

(5-50)

Здесь

(5-51)

Параметр Г (параметр Бури), так же как и формпара-метр f, отражает влияние продольного градиента давления и числа Рейнольдса на профиль скоростей в турбулентном слое.

Зная зависимость для c_f и закон изменения величины

# = ~ в зависимости от продольного градиента давле-

ния, нетрудно получить из интегрального соотношения (5-38) дифференциальное уравнение для определения толщин пограничного слоя.

Рассмотрим отдельно решение уравнения (5-38) для ламинарного и турбулентного слоев.

5-9. РАСЧЕТ ЛАМИНАРНОГО ПОГРАНИЧНОГО СЛОЯ ПРИ НАЛИЧИИ ГРАДИЕНТА ДАВЛЕНИЯ

Расчет пограничного слоя сводится к решению интегрального соотношения (5-38) [или (5-39)], которое содержит три неизвестные величины: толщину потери импульса 8**, коэффициент сопротивления с, = -~ и вели-

^{    Р«5

чину я = я(м;

Следовательно, для решения задачи необходимо иметь еще два дополнительных соотношения, связывающих указанные величины. Используя общее выражение (5-48) для коэффициента сопротивления и рассматривая для простоты случаи несжимаемой жидкости, получим уравнение с двумя неизвестными: 8** и Н, связь между которыми нетрудно получить, если известен профиль скорости в пограничном слое.

В общем случае скорость в поперечном сечении слоя зависит от местного градиента давления и расстояния от стенки у. Влияние градиента давления учитывается форм-параметром f. Следовательно, относительную скорость в точке слоя можно представить в таком виде:

Тогда

00

••“JO-

о

оо

О

оо

⁼⁸**П^1_

о

Проинтегрировав полученное выражение в указанных пределах, найдем:

Далее, подставив в уравнение (5-38) зависимости (5-48) и (5-52), получим после несложных преобразований дифференциальное уравнение первого порядка относительно форм-параметра /:

Sr=^f+ir^F-    (⁵-⁵³>

ал _Uq    и₀

Здесь

F=2{Z(f)-2f[2 + If(f)1}.

Конкретный вид функции F (f) зависит от формы профиля скорости в пограничном слое. Вычисления показывают, что F (f) очень мало отличается от линейной функции, т. е.

F(f) = a~bf.

При этом уравнение (5-53) переходит в линейное относительно параметра / уравнение

v+(f-4>

df

-т~ — ^а dx

интеграл которого имеет вид:

г X    ,

'“о (* ь—1 ,    | ^иа

* dx + c .

Если при л: = 0 и₀ = 0, то из условия конечности f следует положить постоянную интегрирования с = 0.

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

Положив и₀'=~ и х=-^-, где — скорость набегающего потока, a L — длина обтекаемой поверхности, получим:

Постоянные а и b могут быть приняты соответственно равными 0,45 и 5,35.    *

Для толщины потери импульса получим:

*"=х='|/'    <⁵-^м>

Далее, по формуле (5-48) нетрудно определить местный коэффициент сопротивления c_f и по выражению 8* =

Рис 5-16 Зависимость величин $, F и Н от параметра f.

= Н (f) S** — толщину вытеснения. Значения функций С (/) и Н (/) приведены на рис. 5-16.

Рассмотренный метод расчета может быть распространен и на случай течения сжимаемой жидкости, если перейти к новым переменным, предложенным А. А. Дородницыным

Е = \f_tdx и 1={

о    о

⁽⁵-⁵⁶⁾

k + 1 0) о

Расчет по формулам (5-48), (5-54), (5-55) и (5-56) оказывается относительно простым и обеспечивает вполне удовлетворительную точность.

Для примера рассчитаем ламинарный пограничный слой на плоской пластинке. Здесь и₀ = const; и₀ = 0; f = 0 и т=\.

Тогда из (5-55), (5-47) и (5-48) получим:

I f 0,45ч

У ^;

= S**L=L у -^j-; С (/) = 0,22;

1    0,22v    ./~

с. = 0,22 —rjnr =-7=5= = 0,332 У —j >

f ’    и_а8**    Г 0,45м!    У u₀L

/0,45м 3,332 |/^ •

2 _ п ччо 1 / -'^1.    (5-57)

"'о = ^с/Р“о = ⁰

5-10. ПЕРЕХОД ЛАМИНАРНОГО ПОГРАНИЧНОГО СЛОЯ В ТУРБУЛЕНТНЫЙ

При определенных условиях ламинарный пограничный слой теряет устойчивость и переходит в турбулентный. Ориентировочно границу потери устойчивости ламинарного течения можно установить по критическому числу Рейнольдса Re_Kp. Пользуясь аналогией между явлениями перехода ламинарного режима в турбулентный в цилиндрической трубе и в пограничном слое, можно, как это уже указывалось, ввести характерные для слоя числа Рейнольдса, отнесенные к толщинам⁷² 8, 8* и 8**:

Re.= -^; Re* = -^; Re** = ^ .

^v 5    V ’    V    V

Опытные данные показывают, что критические значения числа Re₅ пограничного слоя имеют тот же порядок, что

и для труб, но могут колебаться в более широких пределах (Re₈ — 2 ООО 5 ООО). Согласно многочисленным

исследованиям критическое число Re_Kp зависит в основном от состояния поверхности обтекаемого тела, степени турбулентности (возмущенности) внешнего потока и градиента скорости, т. е. параметра f.

Теоретические и экспериментальные исследования устойчивости ламинарного пограничного слоя показали, что „потеря устойчивости” происходит либо в результате наложения возмущений, вызванных срывом с бугорков шероховатости на поверхности обтекаемого тела, либо вследствие конечных возмущений, вносимых в пограничный слой внешним потоком. Первая причина оказывается существенной только при малой турбулентности набегающего потока (порядка 0,1°/₀) и большой шероховатости.

Используя общие положения теории размерности и опытные данные, нетрудно получить соответствующую зависимость в явном виде. Опуская вырод, запишем формулу для определения Re**, полученную А. П. Мельниковым:

(5-58)

где Е₀ — начальная степень турбулентности; f — формпараметр;

f_s — значение формпараметра в точке отрыва; для ламинарного слоя можно принять /_а = — 0,085.

Таким образом, для определения точки потери устойчивости ламинарного пограничного слоя нужно найти по уравнению (5-55) изменение толщины потери импульса вдоль поверхности рассматриваемого тела и построить кривые Re”=cp(x) [уравнение (5-58)] и Re** =    (х) =

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

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

Re, степени турбулентности и некоторых других факторов.

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

Переходная область может быть определена экспериментально путем измерения профилей скорости в различных сечениях по длине поверхности. Характер изменения профиля скорости в переходной области можно проследить на рис. 5-17, а. На расстоянии х = 35 мм от передней кромки пластины слой — ламинарный; все точки образуют

кривую, соответствующую уравнению ^=2 у — 2^|-j-{-

-j-^У . При переходе в сечение х ^ 50 мм профили

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

Профили скорости в пограничном слое позволяют вычислить толщину потери импульса 6 (х) и по характеру изменения этой величины установить протяженность переходной зоны (рис. 5-17,6).

Длину переходной области s принято определять в долях общей длины обтекаемой поверхности.

Изменения 8 *в этой области характеризуются величиной

**    »**    т-ч

где Re_K, 5_к — число Re и толщина потери импульса в **

**    ТЛ    ***

конце и Re_H=    —число Re и о_н в начале переход-

ной области. При безградиентном течении величина г =

В общем случае г** зависит от профиля скорости в начале переходной зоны, чисел Re и М, степени турбу-

**

лентности. С увеличением чисел Re и М величина г несколько возрастает (рис. 5-18). По опытам МЭИ измене-

Рис. 5-17.

а_изменение профиля скорости в пограничном слое вдоль плоской стенки; б —¦•

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

ние Е„ в пределах 0,5—1,5°/₀ не приводит к заметным изменениям величины г**. Дальнейшее увеличение Е₀ вызывает

резкое уменьшение г .

Для характеристики влияния профиля скорости на рис. 5-19 показан график    (/<,), где /₀—формпараметр

в начале переходной зоны. Как следует из кривой, с ростом формпараметра г" возрастает. Аналогично в зависимости от формпараметра /₀ изменяется и протяженность переходной зоны s (рис. 5-19). Такое поведение кривых г* = <Pj (f₀) и s = <р₂ (/„) можно объяснить следую-

Рис 5 18 Зависимость величины г" от числа М.

щим образом. В переходной зоне вследствие поперечного перемещения частиц толщина пограничного слоя 8 увеличивается, а профиль скорости становится более полным. Если в результате роста толщины 8 происходит увеличение 8", то увеличение полноты профиля скорости вызывает падение интегральной толщины 8 *. В конфузорной области основное значение имеет возрастание толщины пограничного слоя 8, ибо здесь профиль скорости достаточно полный и в результате перехода его полнота изменяется незначительно⁷⁴. Наоборот, в диффузорной области в результате перехода ламинарного течения в турбулентное происходит значительная деформация профиля скорости, причем эта деформация оказывается тем большей, чем больше положительный градиент давления в месте, где происходит переход.

При некотором значении параметра f в диффузорной области оба фактора, влияющие на величину 8'*, взаимно

**

компенсируются и величина г оказывается равной 1. Протяженность переходной зоны при этом также оказывается незначительной. По опытам МЭИ г* = 1 при / =—0,06-=-

— 0,07.

Если же по какой-либо причине переход начнется при — 0,07, то, по-видимому, процесс турбулизации слоя и перестройка профиля скорости будут происходить против течения, пока профиль ламинарного слоя не окажется достаточно устойчивым.

Рис 5-19 Зависимость вел гаин г** и 5 от формпараметра

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

На основании обработки опытных данных для расчета переходной области получены следующие эмпирические формулы:

5^^)°''⁵⁽3,7 + 5,5Шо^:    (5-59)

⁰ 12f;r

г -(7-Ь100/_о)    ² +0,12М.    (5-60)

«-V    --**

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

5-11. РАСЧЕТ ТУРБУЛЕНТНОГО ПОГРАНИЧНОГО СЛОЯ

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

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

(5-61)

Формула (5-61) построена на основании сравнения профилей скорости в турбулентном слое и в цилиндрической трубе. Как было показано, профиль скоростей в трубе меняется при изменении числа Re (рис. 5-6). Следовательно, показатель п в формуле (5-61) является функцией числа Рейнольдса. Так, по опытам Н. Никурадзе показа-

1 1

тель степени п меняется в пределах -g—=- при изменении

Re от 4-10³ до 3-10⁶. Однако в первом приближении показатель степени п может быть принят постоянным в определенном диапазоне чисел Рейнольдса. В расчетах часто

принимают n = Y^{и П}Р°Ф^ИЛЬ скорости задают формулой

(5-62)

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

(5-63)

Подставив опытные коэффициенты С = 0,00655 и т =

— — 0,166 в формулу (5-63) и рассматривая частный случай безградиентного течения, из уравнения (5-38) найдем:

Введем в это уравнение число Re_j, = ^. Заметим, что

Тогда дифференциальное уравнение для толщины потери импульса запишется в таком виде:

^^-=0,00655 Re’*~^0,16G.    (5-66)

Предполагая, что на стенке образуется только турбулентный слой, проинтегрируем уравнение (5-66). Тогда получим:

б_

Re** = 0,0153 Re⁷ .

Отсюда, заменяя Re * и Re_x, находим уравнение для толщины потери импульса:

L в_

8" = 0,0153 X⁷ ,    (5-67)

или, выразив 5 * в функции числа Re^, найдем:

L

8** = 0,0153 a Re~⁷ .    (5-68)

Полученное простое решение для толщины потери импульса в турбулентном слое не учитывает влияния сжимаемости и продольного градиента давления. Оно справедливо при М «5:0,5 и ^⁼⁼0.

Как уже отмечалось, уравнение импульсов (5-38) связывает три неизвестные величины: 8*, 8 и т₀. Дополнительные связи, необходимые для решения задачи, устанавливаются опытным путем.

На рис. 5-20 приведена экспериментальная зависимость Н₀ = Й (М) для плоской пластины при нулевом градиенте давления. Величина Н₀ существенно зависит от числа М; приближенно эту зависимость можно представить следующей формулой:

я„ = я₀₀ (1 + аЩ = Я ₀₀ ( 1+0,3 М°),    (5-69)

где для несжимаемой жидкости при нулевом градиенте давления можно принять //₀₀= 1,3-т-1,4, причем меньшее

значение Н₀₀ соответствует большим числам Re. С увеличением числа М несколько уменьшается относительная толщина потери импульса (рис. 5-26)

Влияние градиента давления на параметр Н можно видеть на рис. 5-21, где дана зависимость Н=Н/Н₀ от параметра Бури Г.

Отсюда следует, что в зоне умеренных значений параметра Г(—0,015<;Г <0,02)И изменяется на 7°/о- Поэтому при построении решения в отмеченном диапазоне будем учитывать только изменение И₀ = Н(Щ.

Для коэффициента сопротивления в турбулентном слое для общего случая была получена формула (5-50). Если отнести т₀ к плотности на внешней границе пограничного слоя р_х, то формула (5-50) может быть представлена в форме:

"W Т'»    “

¦тг^Ке-

*(С + аГ) =

Рх“0

^+*r)Re:-^m(i¦ *¦

_1м²

2 о

(5-70)

Рис. 5-20. Влияние числа М на параметр /7₀.

Рис. 5-21. Изменение величины rl в зависимости от параметра

Бури Г.

Подставим (5-69) и (5-70) в интегральное соотношение (5-38). После несложных, но кропотливых преобразований можно получить дифференциальное уравнение относительно параметра Г в следующем виде:

Г' -\-р (а) Г -f- Q (а:) = 0,    (5-71)

где х — безразмерная длина обтекаемой поверхности.

Уравнение (5-71) относится к типу линейных и может быть приведено к квадратурам и решено относительно

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

При равенстве нулю постоянной а и стремлении X к нулю решение переходит в формулу (5-67).

Введя ряд упрощений, получим при небольших градиентах давления (1’> — 0,02) для толщины потери импульса выражение

ГТ**¹.25₉4.17г₎ 0,2в/_{е    5}2    >1,81 ,

2лЖ[⁶о К Re. (⁶ — Q )    +

6 =-

Х<*-³⁵ Re°'²(6-X*)

_т0,8

(5-72)

30 26 гг /в /4 ю 6,0 г,о 1,6 1,2 0,8 ОЛ

+ 0,0026 j я³⁹² (6 — Х\ j²-⁸¹ dx]

0,2 О,U 0,8    0,8    1,0    1,2    1,Ь

⁶)

б—чависимость функций <р,, у, и <р_а от X,. 17 М. Е. Дейч

и безразмерной скорости в начале турбулентного участка;

a L    _u

Re =—--число Рейнольдса, определенное по критической

*    ~*w

скорости а и кинематической вязкости на стенке.

Вводя обозначения:

^о^,35(⁶-*о)^М5;

h

/_а = 0,0026    (6    —Я² f

(5-73)

преобразуем (5-72) к виду:

^г=ш* [+1^1“ ⁽⁵'⁷⁴⁾

Функции /₂ и fз от Я₀ представлены на рис. 5-22,а.

Расчет по формуле (5-74) с использованием расчетных графиков оказывается относительно простым и дает хорошее совпадение с опытными данными.

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

0,2 L °о

?,Re*

х_в

Здесь

₂2,35,_fi    ,2    Л    3    °^,10оХ0.

^:Яо (⁶~^яо) ^е ;

?!¦

(5-76)

О

<р₂ = 0,0078 Яд'⁶⁹ (6 — Я²)²⁶² е °_;

?»= ^Яо⁹⁴ (6 — Яо )

Г, 131X2

2 U.62

Соответствующие значения ср, ср₂ и <р₃ приведены на рис. 5-22,6.

В заключение остановимся на существенном различии в свойствах ламинарного и турбулентного слоев. С этой целью в табл. 5-1 приведены основные расчетные формулы

для простейшего случая обтекания плоской стенки (без' градиентное течение) несжимаемой жидкостью.

Сопоставление показывает: 1) профиль скоростей в турбулентном слое более полный, чем в ламинарном; 2) толщина турбулентного слоя растет вдоль стенки значительно быстрее, чем ламинарного, так как в первом случае 5 уве-

6/7

личивается пропорционально л: , а во втором — пропор-ционально я^,/2; 3) сравнение местных коэффициентов со-

Таблица 5-1

противления трения показывает, что при одинаковых значениях Re^. сопротивление трения в турбулентном пограничном слое значительновыше, чем в ламинарном. •

Этот весьма важный вывод хорошо подтверждается опытными данными. На рис. 5-23 приведены экспериментальные и расчетные зависимости коэффициентов c_f, полученные для турбулентного режима с помощью различных по-луэмпирических формул. Кривые показывают достаточно ' высокую точность формулы, приведенной в табл. 5-1.

OMtn

Сравнение c_f для ламинарного и турбулентного слоев показывает, что при одинаковом значении 1^= 4-10⁴ коэффициент с_/турб в 2 раза превосходит c_fmu

(^турв ~ ^,0054, ^С/лам ⁼ = 0,0022). С ростом Re* разность с_/турб-с_/лам резко увеличивается.

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

Обнаруженная разница в сопротивлениях трению при

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

Многочисленные опыты в соответствии с выводами полуэмпирической теории указывают на еще одно существенное различие в свойствах ламинарного и турбулентного пограничных слоев. Точка отрыва ламинарного слоя при больших числах Re занимает фиксированное положение на обтекаемой поверхности. Точка отрыва турбулентного слоя при одинаковом распределении параметров внешнего потока располагается дальше точки отрыва ламинарного слоя. Другими словами, турбулентный ’слой отрывается позднее ламинарного. Это становится понятным из рассмотрения профилей скорости при одном ¦и другом режимах слоя: кинетическая энергия частиц, движущихся вблизи стенки, при турбулентном режиме значительно больше, чем при ламинарном.

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

Изложенная методика расчета не учитывает влияния начальной турбулентности потока. Для небольших степеней турбулентности уравнение импульсов пограничного слоя для несжимаемой жидкости получено В. А. Врублевской в таком виде:

^_' $**

(1+А?о)?г + "^-(2 + Я)=-^- + ?§.    ⁽⁵'⁷⁷⁾

^и    о    Р ^и0

стенке, проведенному в МЭИ. Графики на рис. 5-24 отчетливо показывают влияние ?₀ на потери в пограничном слое; так, при Re.,. =

= 2-10⁶ увеличение ?₀ от 0 до 5% повышает число Re**, а следовательно, и величину потерь на трение в слое на 70%.

Рис. 5-24. Влияние степени турбулентности на толщину потери импульса в турбулентном пограничном слое.

5-12. ПОГРАНИЧНЫЙ СЛОЙ ПРИ БОЛЬШИХ СКОРОСТЯХ.

ЭКСПЕРИМЕНТАЛЬНЫЕ ДАННЫЕ

Рассмотрим некоторые результаты опытного исследования турбулентного пограничного слоя при больших дозвуковых скоростях при безградиентном течении и докрити-ческих скоростях (М •< М ).

На рис. 5-25 приведены шесть профилей скорости, полученные при постоянном числе Re, равном 2,5-10®, и переменном числе М. Все экспериментальные точки в интервале М = 0,31-т-0,98 образуют одну кривую, которую можно представить зависимостью вида;

?⁼Ч"г)-    ⁽⁵-⁷⁸⁾

Следовательно, при нулевом градиенте давления с изменением числа М профили скорости в координатах у, Я изменяются подобным образом. Кроме того, в указанном диапазоне изменений числа М толщина пограничного слоя при постоянном числе Re изменяется незначительно; поэтому можно говорить не толь- , д

ко о подобии, но и о тождественности профилей скорости. Отсюда можно заключить, что непосредственно на профиль скорости в пограничном слое влияние сжимаемости при М < 1 сказывается незначительно. Таким образом, в первом приближении профиль скорости при М > О можно представить обычной степенной формулой (5-62).

0,9

0,7

0,5

0,4

0,3

0,2

О

0,9 КО

0,5    0,6    0,7    0,8

Рис. 5-25. Влияние числа М на профиль скорости в турбулентном погра

0,8

0,6

Формула (5-62) не является, конечно, единственной. Любая зависимость, удовлетворительно аппроксимирующая профиль скорости в пограничном слое несжимаемой жидкости, может ,быть распространена и на течение сжимаемой жидкости. Пределы такой экстраполяции не ограничиваются, по-видимому, звуковой скоростью и могут быть распространены и на сверхзвуковые скорости.

0,1

ничном слое.

Дальнейшая обработка профилей скорости дает возможность вычислить значения интегральных толщин 8 ,

о и построить кривые их изменения в зависимости от числа М (рис. 5-26). Разброс экспериментальных точек в области околозвуковых скоростей объясняется возникновением конфузорного течения. Кривые для околозвуковых скоростей проведены через точки, полученные при минимальном градиенте давления; поэтому в области звуковых и сверхзвуковых скоростей основная масса точек располагается под кривыми. Разброс точек, обнаруживаемый в области дозвуковых скоростей, объясняется тем, что здесь нанесены точки, относящиеся к различным числам Re.

Полученная картина находится в полном соответствии с теоретическими результатами. Увеличение числа М при-

мм

о,з&

0,30 0,26 0,22 0,18 0,14 0,10

1,1

1,0

0,3

0,8

0,7

Рнс. 5-26. Влияние числа М на толщину потери импульса.

водит к падению толщины потери импульса и к росту толщины вытеснения 5*. Такое поведение кривых при одинаковых профилях скорости обусловлено изменением плотности поперек пограничного слоя, причем противоположное влияние этого изменения на толщины 8 и 8 объясняется тем, что подынтегральные функции, через которые выра-жаются толщины о,5, по-разному реагируют на изменение плотности.

Для приближенной количественной оценки влияния сжимаемости на толщину потери импульса 8 на рис. 5-26

8**

приведено изменение относительной толщины —= 8" в

^йоо

зависимости от М при нулевом градиенте давления. С приближением к звуковой скорости уменьшение толщины потери импульса составляет около 15°/₀.

Независимость профиля скорости от числа М дает основание считать число Re основным параметром, определяющим профиль скорости при отсутствии градиента давления, и распространить результаты многочисленных экспериментов по определению его влияния в несжимаемой жидкости на течение сжимаемой жидкости.

0,6    0,7    0,8    0,3    1,0

Рис. 5-27. Влияние числа Re на профиль скорости в турбулентном пограничном слое.

1,0

0,9

0,8

0,7

0,6

0,5

а Л о,э о,г о,/

о 0,5

На рис. 5-27 приведены шесть профилей скорости, полученные при изменении числа Re и постоянном числе М. Увеличение числа Re от 0,61*10“ до 1,08-10“ приводит к характерному изменению профиля скорости. Дальнейшее увеличение числа Re не вызывает заметного изменения профиля скорости. Другими словами, начиная с величины

Re ж 10^е, профили скорости могут быть выражены общей зависимостью.

Исследование влияния числа Re на интегральные толщины 8* и 5** показало, что с ростом Re, как и в случае несжимаемой жидкости, эти величины уменьшаются. В области малых значений Re изменение §’* оказывается довольно заметным. С увеличением числа Re уменьшение толщины 8*’ происходит менее интенсивно. Аналогичным образом меняются и толщины вытеснения.

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

Для исследования вопроса о влиянии числа М при наличии градиента давления были сняты профили скорости в диффузорной области (рис. 5-28,а) и при конфузорном течении (рис. 5-28,6). Как видно из приведенных кривых, все экспериментальные точки независимо от числа М образуют практически одну кривую. Этот факт еще раз подтверждает сделанный ранее вывод о том, что при дозвуковых скоростях изменение числа М не приводит к заметным изменениям профилей скорости. Следовательно, основным фактором, определяющим течение в пограничном слое сжимаемой жидкости, является продольный градиент давления.

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

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

На рис. 5-30 приведена кривая изменения толщины потери импульса в функции параметра Г, полученная из эксперимента. Чтобы исключить влияние сжимаемости, тол-

0J

D