Научная статья на тему 'К подобию гиперзвуковых cтpуйныx течений'

К подобию гиперзвуковых cтpуйныx течений Текст научной статьи по специальности «Физика»

CC BY
155
51
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Ученые записки ЦАГИ
ВАК
Область наук

Аннотация научной статьи по физике, автор научной работы — Гусев В. Н., Климова Т. В.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «К подобию гиперзвуковых cтpуйныx течений»

УЧЕНЫЕ ЗАПИСКИ ЦА Г И

Том III 1972 № 6

УДК 532.525.2

К ПОДОБИЮ ГИПЕРЗВУКОВЫХ СТРУЙНЫХ ТЕЧЕНИЙ

В. Н. Гусев, Т. В. Климова

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

Вопросы подобия гиперзвуковых струйных течений рассматривались в ряде работ [1] —[3]. Ниже аналогичное исследование проводится на основании теории сжатого слоя при условии, что теч^г ние в струе эквивалентно одномерному [4] и [5]. Такой подход позволяет получить информацию о форме сильно недорасширен: ных и перерасширенных струй, истекающих в покоящуюся среду и в спутный гиперзвуковой поток.

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

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

струя будет ограничена слоем сжатого в скачке уплотнения газа. В рамках теории сжатого слоя для разности давлений Ар по нормали к сжатому слою можно записать [4]:

У ]

/? Vs

/

Aw Гс

О fa

Фиг.

А р= — . ... (1.1)

2 тг rsN sin ср

На фиг. 1 г, <р— полярные координаты с началом в точке О; Q—суммарный расход газа через кольцевой сжатый слой; Vm —

максимальная скорость; г, — координата скачка.; Л^—радиус кривизны скачка. ' "

Для контура скачка уплотнения из (1.1) получим следующее дифференциальное уравнение:

Р 2 ъг1р0 Г=/+2±г- 8ІП ср (/* + Р)

Ра 1 +

т

1

■м

т

т-1

т -И \т — і

1_ __1 т-1 Т—11

ЄІП* еМ!

(1.2)

здесь = ра!рс) — (9 — 60 — угол наклона скачка к направлению набегающего потока; («• = аг^( — ^//') — угол, образованный направлением радиус-вектора г и касательной к контуру скачка; /(?) —г5/го'> М и 0 —чцсло М и угол наклона вектора скорости к оси л-; рс, Мс — соответственно давление и число М в точке А\ р0 — давление торможения; ? — отношение удельных теплоемкостей; индекс „1“ относится к состоянию перед поверхностью скачка уплотнения, штрихом обозначено дифференцирование по независимой переменной ср.

При гиперзвуковых скоростях из условия эквивалентности течения в струе одномерному для входящих в уравнение (1.2) функций Q, М1г 0, монсно записать соотношения

0. = <о-

М1

(1±1 Ь-1

т+»

4

г, \т-1

я*

(1.3)

в которых р0 — плотность торможения; 11% и <3 — функции аргумента <р, зависящие от распределения параметров потока в начальном сечении струи; определяет радиус критического сечения эквивалентного источника.

С учетом (1.3) уравнение (1.2) преобразуется к виду

/"=]+ 2

Р 1—1 ЭШ <р (/2 +/'2)3/2 /

27

і <2

Т Л + 1

л-1 — 1\ 2 (т-1)

Ра\ 1 +^

-М2с

т-1

я*

я*

г

о /

т + Нт-1/ \т + 1/ Р+Р.

Граничные условия задачи следующие:

/=1, Д/? = 0 при <р = <Ро-

При Ра<^1 второе условие переписывается [5] в виде

'т + 1\1/2

Г

т -ь іу/2 — 1

агс1§

1 — 1V/* Т+1

Ї-1 Км2с-

1

1 — ап^ Км? — 1 |

(1.4)

(1.5)

(1-6)

где р — угол наклона вектора скорости к оси х в точке А (см. фиг. 1). .

Таким образом, задача определения формы струи сведена к задаче Коши для уравнения (1.4) с граничными условиями (1.5). При фиксированной форме сопла задача содержит четыре безразмерных параметра Ра, Мс, р и *[.

2. Остановимся на предельном состоянии течения при Ра-+ 0. Поместив начало координат в выходном сечении сопла (ср0 = ^/а» г0 — гс, где гс — радиус выходного сечения сопла), перейдем в исходных уравнении и граничных условиях к новой зависимой переменной Р=/Р0-1/2.

Тогда

= р + 2 £1 - Т - 1 8Ш ?(/* 4-

2?

JL

т-1

Т

т-1

1

Q(?iMc, р, т)

&(<?, Мс, р, -г)

T + iVi-1) \t + i; f* + f*

F = F' — 0 при ? = тс/2,

(2.1)

где

Po-PolPa=Pa: 1 +

-Mc2

T

7-1

В новых переменных соотношения (2.1) уже не содержат параметра Р0 и форма струи определяется соотношением

/= Ро2F(<?, Мс, р, Т).

(2.2)

Аналогичный результат был получен в работах [1] и [2J на основании теории размерности. Приведенный способ доказательства закона подобия для гиперзвуковых струйных течений идеального газа позволяет получить интересные количественные сведения о структуре течения в сильно недорасширенных струях.

Рассмотрим некоторые примеры. При истечении струи в затопленное пространство ра = const при р = 0 значения входящих

в (2.1) функций /?* и Q были получены в работе [5] на основании

точных численных расчетов методом характеристик. В частном случае при Мс= 1 значения ft* хорошо аппроксимируются зависимостью [6]:

R* (?, 1. 0, -[) = Я* (0) cos а<р; (2.3)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

значения /?* (0) и а приведены в

таблице.

С помощью (2.3) легко вычисляется расход газа через сжатый слой

г

2 y-i/т __ П1/2

7 Д*(0) а

5/3 1,36 1,15

7/5 1,15 0,94

9/7 1,04 0,83

где

Q(?, 1, 0, т) =

Т + 1

1

+

а

1—2 а

1 -f-2 CL

cos (1—2 а) ср

Т + 1,

(а cos2 а<?) cos ср +

!■

(2.4)

Отметим, что приведенное выше выражение ДЛЯ <7 определяет расход газа через сжатый слой, отнесенный к полному расходу газа. _

Аналогичным образом могут быть определены значения /?* и (2 при произвольных значениях Мс и р. _

Определив входящие в (2.1) неизвестные функции /?* И (3, можно из уравнения (2.1) определить функцию Р. Однако функцию F—fPo1|2 можно определить и не решая этого уравнения, если обработать в параметрах подобия полученные в работе [5] решения уравнения (1.4) при достаточно больших значениях Р0. Результаты такого пересчета приведены на фиг. 2, на которой при различных значениях числа Мс н •[ представлены координаты скачка уплотнения X = х^г,. Р<Г1/2, У =^/гс Ро'12, т = 7/5 —• сплошные линии,

1 = 5/3—пунктирная линия. Соответствующее сравнение теоретического расчета с экспериментом представлено в работах

[2] и [5].

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

[3], геометрическая картина истечения будет определяться двумя характерными размерами Яг и /?2. В рассматриваемом

здесь случае равномерной скорости на выходе из сопла ((3 = 0) величины /?! и /?2 выражаются кчерез параметры на срезе сопла следующим образом [3]:

05

1-0

Г

г X 7 —

5 — . 1. .

0.5

10

Фиг. 2

Ъ = геМе'а}/ +

1 \1/4

где

тМ I

1 +

, /?,= /?,(! -7,)1'2,

1 Х-1*'2

(2.5)

1 Мс

На фиг. 3 сравниваются полученные в рамках теории сжатого слоя геометрические размеры струи при 1 = 1,4 и р = 0 (сплошные линии) в новых переменных с приведенными в работе [3] данными

численного расчета формы ви-

&

05

Мс=5-,1>^ 3,7-10 *

Мр3-,1=и-ю^ -"г

V г 1 1 N

/ К м^3-,1а=3-м’3 □о ч

р*

0.5

10

*з!я1

Ф иг. 3

сячего скачка уплотнения для конических сопл с (3 = 10°. Сравнение указывает на приемлемую точность метода сжатого слоя. Видно также, что при малых углах раствора сопла различие в начальном распределении параметров потока не приводит к изменению формы струи в указанных на фиг. 3 безразмерных переменных.

3. Рассмотрим теперь осесимметричную гиперзвуковую

струю газа, истекающую из расчетного или перерасширенного сопла с коническим потоком на выходе. Как и в предыдущем случае, течение в струе будет одномерным, однако в этом случае ft* = const и входящие в уравнение (1.4) функции ft* и Q Могут быть представлены в конечном виде. Полагая^<р0 = Р, r0 = rc/sinp (р— угол полураствора сопла), при Мс^> 1 для /?% и Q можно записать:

Т±1

'Т + 1М-1

Q = 2

ft* —

1

М„ т-1

і

Y—ї ^

Мс т_1 (cos f — COS Р) .

(3.1)

С учетом (3.1) уравнение (1.4) преобразуется к виду

2Т 1

//3 sin ?(Г+П*ЧРа f J f т (cos <Р — cos Р) V Me

Т + 1 /2+/'2

Граничные условия задачи следующие:

М2

(3-2)

(3.3)

Т+1

При ®<:р<С1 имеем /'^>/. Переходя в уравнении (3.2) и граничных условиях (3.3) к новым переменным Ф =/КРа(Мср)-1 и ср = <РР~1, для определения формы струи получим соотношения

ФФ 2 (1Ф 2ш Г / й?Ф\2 2 т

" Т 4- 1

d ср2

Ф d<?

Ф = Ф0 = Kft0(Mc Р)-1;

_2у Г

Т(ф2- 1) L dФ

dy

\ df V-

Т + 1

при «р = 1.

(3-4)

Выписанная система соотношений (3.4). определяющая в полярных координатах форму перерасширенной струи с коническим потоком на выходе из сопла, содержит два безразмерных параметра Ф0 и т- Это значит, что для подобия таких струй, соответствующих различным значениям параметров Ра, Мс и р, необходимо равенство критериев подобия Ф0 и т- Аналогичный результат был получен в работе \7\, в которой для расчета перерасширенных струй использовался закон плоских сечений в методе пограничного слоя [8].

Как и в случае сильно недорасширенной струи, при Мс^>1 продольный размер перерасширенной струи оказывается существенно больше поперечного:

X, = гс/~ ; уг = гс/-£- . (3.5)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Выбирая в качестве характерных продольного и поперечного размеров гсМсРа~12 и ГеМсРЯГ^2, для формы перерасширенной струи окончательно получим соотношения

Функция Ь была найдена в результате решения уравнения (3.4) и представлена на фиг. 4 при нескольких значениях параметров Ф0 и ■[. Там же даны предельные значения функции Ь, соответствующие случаю Фо = 0. Сравнение теоретического расчета формы перерасширенной струи с результатами эксперимента при

двух значениях Ф0~1, проведенное в работе [4], указывает на приемлемую точность метода сжатого слоя.

4. Остановимся теперь на истечении струи в спутный гипер-звуковой поток (Moo^l). В этом случае давление на внешней границе сжатого слоя ра будет переменным, и для определения его величины можно воспользоваться, например, формулой Ньютона, справедливой при a>(i:

Ра (/» Ч>) = Р» vlo sin2 (<р — ц) = х/?0о Sin2 (f — у.), (4.1)

где х, Vos, poo, poo — отношение удельных теплоемкостей, скорость,

давление и плотность в невозмущенном потоке. В начальной точ-

ке сжатого слоя при <р = ср0 и /=1 давление

о /'2

Ра0> ?о) — Ра — *Р<Х> Moo J _j_ jrri } (4.2)

где /' определяется соотношением (1.6).

В случае сильно недорасширенных струй ^Рв = 1 j, под-,

ставляя (4.1) в (1.4), получим

/2 т — 1 Sin (/2 + f'2)112

f”=f 4-2-

/

Q

*М£оР01 (/'sin f + / cos <р)2

2 \'т-1 /т _ n 2 (т-1)

Т + 1 Vt — 1Л Vt + 1

я

('*о = /*с; Ро = 'Ро!р»)- (4.3)'

Фиг. 5

Граничные условия задачи совпадают с (1.5). Таким образом, при фиксированной форме сопла задача определения формы струи, истекающей в спутный поток, зависит от шести безразмерных параметров: Р0, Мс, Моо, Р, т. *•

Как и ранее, перейдем к новой зависимой переменной F* = — ]/~>ГМоо Р(Г1/2/. Тогда после замены найдем, что преобразованное уравнение (4.3) и граничные условия (1.5) не. будут зависеть от Р0, Моо и х. Форма контура струи, обтекаемой гиперзвуковым потоком, будет определяться соотношением

f-

V х мс

Мс, р, т).

(4.4)

Аналогичный результат был получен в работе [2] из теории размерностей. Геометрические размеры струи в параметрах подобия при 1 = 7/5, (3 = 0 и нескольких значениях числа Мс приведены на фиг. 5. Эти данные были получены в результате численного решения уравнения (4.3) при достаточно малых значениях Рй[5].

Так же как и при истечении в затопленное пространство, при Мс^>1 продольный размер струи, истекающей в гиперзвуковой спутный поток, оказывается намного больше поперечного. В этом случае при рассмотрении геометрического подобия таких струй можно воспользоваться результатами работы [3], приняв в качестве характерного внешнего давления давление р°а в начальной точке сжатого слоя. Его величина, определенная соотношением (4.2), при равномерной скорости потока на выходе из сопла р = 0 и Мс 1 будет

о 4х /МссЛ2

(т-1)2/Ммс

/>в =

(4.5)

и введенные в работе [3] характерные продольный /?, и поперечный /?2 размеры (2.5) в случае струи, истекающей в спутный поток, примут вид

R i =

Т — 1 Мс

Vi

2

1 У'4

(4.6)

- Моо

/Й = ЯГ(1 -Л)1/2.

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

В заключение рассмотрим осесимметричную гиперзвуковую струю газа, истекающую в спутный поток из расчетного или перерасширенного сопла (Ра — РаІРс > 1) с коническим ПОТОКОМ на выходе. В работе [7] этот случай рассматривался на основе закона плоских сечений в методе пограничного слоя [8].

Фиг. б

Подставляя выражение (4.1) для давления на внешней границе сжатого слоя в уравнение (3.2), получим

/"= f_LO f— - sin?(/3 + /'2>1;2 J ■' f f(COS<p — COSp)

^Poo M0 if

(sin rff' + cos?/)3

2 t

•(4.7)

pcMl ' ‘J ' ‘J ' T +1. Граничные условия задачи совпадают с (3.3). При Р<С!1 и ML/M? PcIPoo < *м£, р2 величина р°а в начальной точке сжатого

слоя будет

Ра = Х/»оо М^о Р2 , и граничные условия задачи (3.3) примут вид XI £> Мс 1 Л рс

/= I; / = - щ у -г-гтч-р- ПРИ т

: Р'

(4.8)

(4.9)

Мотр V *(7 + 1)/>

Переходя в уравнении (4.7) и граничных условиях (4.9) к но-

Моо

вым переменным Ф*=/

МР

Рс

и «р = 9р 1 и учитывая, что

при 9<р<с;1 имеет место /'^>/, для формы струй получим соотношения

d

ф

Ф*______I dФ* \2 2? d<b* Г/-

і2 Ф* V dip ) df

- do*

Фп =

Mo

Poo

Рс

d ф=>

d<?

df

+ Фч

при

4 + 1

9=1.

(4.10)

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

Фо и т. Выбирая в качестве характерных продольного и попереч-

мс г~1г~ Мс г~~а~

ного размеров гс ]/ —и гс щ у _s_, для определения

формы струи, определяемой соотношениями (3.5), окончательно получим соотношение

ys Мро f Pn I xs Moo $ f pa \

^МГУ*7Г = /-Ч^ТЯГ V Л : Фб; Ф (4Л1)

Значения функции Z.*, вычисленные в результате решения

уравнения (4.10) для предельного случая Фо = 0, приведены на фиг. 7.

ЛИТЕРАТУРА

1. Moran I. P. Similarity in high-altitude jets, A1AA Journal, v. 5, No 7, 1967.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2. Гусев В. H., Михайлов В. В. О подобии течений с расширяющимися струями. .Ученые записки ЦАГИ-, т. 1, № 4t 1970.

3. М у р з и н о в И. Н. Параметры подобия при истечений

сильно недорасширенных струй в затопленное пространство. Изв. АН СССР, МЖГ, 1971, № 4. *

4. Гусев В. Н. К расчету гиперзвуковых осесимметричных струй. Изв. АН СССР, МЖГ, 1967, № 1.

5. Гусев В. Н., Климова Т. В. Течение в истекающих из недорасширенных сопл струях. Изв. АН СССР, МЖГ, 1968, № 4.

6. Ashkenas Н-, Sherman F. S. The structure and utilization of supersonic free jets in low density wind tunnels. Dynamics of Rarefied Gases, Forth Symp. Acad. Press, 1965.

7. За к Л. И. Гиперзвуковая струя, истекающая в покоящуюся среду или в спутный сверхзвуковой поток. Изв. АН СССР, МЖГ,

1969, № 5. ’ ’

8. Черный Г. Г. Течение газа с большой сверхзвуковой скоростью. М., Физматгиз, 1959.

Рукопись поступила 4\IV 1972 г.

i Надоели баннеры? Вы всегда можете отключить рекламу.