ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ НЕТИПИЧНОГО ОТРЫВА ПОТОКА ГАЗА В СОПЛЕ С КОНИЧЕСКИМ НАСАДКОМ Текст научной статьи по специальности «Физика»

Электронный журнал «Труды МАИ». Выпуск № 40

www.mai.ru/science/trudy/

УДК 533.17, 534.612

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

В. В. Комаров, А. А. Пономарев, Н. Б. Пономарев

Аннотация

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

Ключевые слова

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

Проведение наземных испытаний жидкостных ракетных двигателей (ЖРД) с соплами большой степени расширения всегда сопряжено с определенными трудностями [1]. Изготовление и эксплуатация газодинамических труб, позволяющих создавать разрежение у среза сопла, для больших (в особенности многокамерных) ЖРД требуют огромных затрат, поэтому часто наземные испытания проводятся без понижения внешнего давления. В этом случае поток газа в сопле отрывается от стенки сопла, причем формируется либо типичный отрыв потока [1], либо нетипичный отрыв потока [2-10]. Впервые нетипичный отрыв потока был обнаружен в США при испытаниях модельного сопла двигателя J2S [2], а затем при запуске и останове двигателей SSME и ^^шп [3-10]. В СССР и России такой отрыв не наблюдался ни в штатных, ни в модельных соплах и, соответственно, не исследовался.

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

поверхности сопла, на котором статическое давление остаётся меньше внешнего и создаёт таким образом отрицательную тягу. Невысокая точность определения местоположения отрыва, а также давления на внутренней поверхности сопла ниже по потоку от места отрыва, затрудняет определение удельного импульса по измеренной при таком отрыве тяге [1]. При нетипичном отрыве оторвавшийся поток присоединяется к стенке сопла и движется вдоль нее с небольшой сверхзвуковой скоростью [8]. В ЖРД это приводит к существенному повышению тепловых потоков в стенку сопла на этом участке, что может повредить неохлаждаемую часть стенки. Поскольку место отрыва потока колеблется, то также возникает нестационарная боковая сила, действующая на сопло и способная в некоторых случаях разрушить камеру ЖРД. Данные эффекты существенно ограничивают степень расширения сопел ЖРД, используемых в широком диапазоне изменения внешнего давления (например, ЖРД второй ступени РН, работающей со старта), соответственно ограничивает эффективность таких ЖРД [1].

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

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

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

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

Исследуемое модельное сопло большой степени расширения (е = ¥а /¥*=198, где ¥а,

¥* - соответственно площади выходного и критического сечений сопла) было

а)

0 3 6 В 12 15 18 21 24 27 30 33 36 Ш' 39

Рис. 2. Изолинии числа Маха с шагом 0,1, начиная с М=1,1, для модельного

сопла при а) п=1,17; б) п=1,4.

спрофилировано методом характеристик с заданием начального участка сверхзвуковой части сопла дугой окружности и с равномерной выходной характеристикой [1] для идеального газа с показателем изоэнтропы п=1,17. В экспериментах, проведенных на холодном воздухе, п=1,4, поэтому изолинии числа Маха в сопле (без насадка) существенно отличаются от рассчитанных при п=1,17 (рис. 2). Очевидно, что при таких числах Маха возможна конденсация холодного воздуха (составляющих его газов) в сопле, но влияние этого эффекта на отрыв потока в сопле пренебрежимо мало.

Исследуемое модельное сопло имеет радиус критического сечения г* =5 мм и радиус

выходного сечения уа =70 мм (рис. 3). На срез сопла был установлен сужающийся

конический насадок с углом наклона контура к оси сопла 45° (рис. 1). В ходе исследований радиус его выходного сечения ун постепенно

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

е н =(.Ун/ ) =81,5 до 8 н =139,2.

В каждом эксперименте измерялось статическое давление на внутренней поверхности сопла рм1 в 18 точках, приблизительно равномерно распределенных с шагом Ах / г* = 1,5...2,0 вдоль образующей сопла (в шахматном порядке, с отклонением 5 мм

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

Таблица 1.

Расположение точек для измерения статического давления

№ точки Р14 Р13 Р12 Р11 Р10 Р9 Р8 Р7 Р6

х / г* 8,49 10,01 11,39 12,77 14,23 15,62 17,06 18,51 19,86

№ точки Р5 Р19 Р18 Р17 Р16 Р4 Р3 Р2 Р1

х / г* 21,16 23,4 25,4 27,4 29,4 31,4 34,5 35,92 37,44

Кроме того, измерялось внешнее давление ра около выходного сечения насадка (сопла, если насадок отсутствовал), условно показанное на представленных здесь графиках при х / г*=38,4. Построенные по экспериментальным данным распределения относительного

статического давления по внутренней поверхности сопла рк/р0, где р0 - полное давление газа, измеренное перед соплом, показаны на рис. 4в, 5в, а изменение р^р0 в точках

измерения в ходе эксперимента - на рис. 4б, 5б. Для измерений использовались датчики абсолютного давления с диапазоном 1.760 Торр, а также датчики избыточного давления с диапазонами 0.0,5 МПа и -0,1.0,1 МПа; погрешность измерения составляла 2% от верхней границы диапазона.

Рис. 3. Модель сопла большой степени расширения

Кроме того, в 8 точках (см. табл. 2) на внутренней

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

расширения е н =139,2

и без него измерялись пульсации

статического давления с помощью

высокочастотных датчиков давления Бпёеуео 8306А-50. Расположение этих точек было выбрано так, чтобы захватить место отрыва потока и его окрестности.

Частота опроса

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

В каждом

эксперименте полное давление р0 холодного воздуха постепенно или резко увеличивалось до 10...11 МПа и также измерялось (рис.4-6).

в)

Рис. 4. Результаты эксперимента без использования насадка: а) изменение во времени отношения полного давления газа перед

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

Все эксперименты

проводились при внешнем атмосферном давлении

ра « 750 мм. рт. ст. Всего было проведено 45 экспериментов: 7 без насадка и 38 с коническим насадком: 19 с насадками с ен <139,2, и 19 с насадком с

е н =139,2.

Для визуализации

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

Также на стенде БМР ФГУП «Центр Келдыша» были получены теневые картины струи около выходного сечения сопла без насадка, на которых видны ударно-волновые

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

Во всех экспериментах без насадка распределение давления по стенке сопла соответствовало типичному отрыву потока (рис. 4), хотя безотрывное течение при п=1,4 существенно отличается от расчетного для п=1,17 (рис. 2). Во

в)

Рис. 5. Результаты эксперимента с использованием насадка с е н =81: а) изменение во времени отношения

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

точках статического давления на стенке сопла, отнесенного к полному давлению газа перед соплом; в) распределение относительного давления вдоль стенки

сопла

В)

Рис. 6. Результаты эксперимента с насадком с ен = 139,2: а) изменение во времени отношения полного давления газа перед соплом к внешнему давлению; б) изменение во времени давления на стенке сопла в разных точках, отнесенного к давлению в камере сгорания; в) распределение относительного давления вдоль стенки сопла

всех экспериментах с коническим насадком при увеличении относительного полного давления р0 / ра отрыв потока был типичным при

Р0/ Ра <(Р0/ Ра )* и

нетипичным при

Р0/Ра РР0/Ра )*, где

(Р0 / Ра) * - некоторое

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

Р0/Ра =(Ро/Ра )** <( Р0/Ра ) *, т.е. с некоторым гистерезисом.

Влияние относительного радиуса выходного сечения насадка ун/ г* на (р0/ Ра )* и

(Р0 / Ра) **, а также на значение

(Р0 / Ра )14, при котором течение

около ближайшей к

минимальному сечению сопла измерительной точки №14 становится безотрывным,

показано на рис. 7. Уменьшение

(Р0/ Ра )14 при увеличении ун / г* от 9,0 до 11,8 объясняется

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

[3-10] распределения имеют, кроме основного пика, расположенного около места отрыва, еще 5 пиков меньшей амплитуды. Все эти пики давления стабильны, их величина и местоположение при одном и том же насадке зависят только от р0 / ра .

ео

50

40

30

-

9.5

10

10.5

11

УнЕ*

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

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

Таблица 2.

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

х / г* 10.0 10.7 11.4 12.2 12.8 14.2

№ точки по образующей А 6А 5А 4А 3А 2А 1А

№ точки по образующей Б 3Б 1Б

Пульсации давления в точках 4А, 5А, 6А в сопле без насадка при увеличении относительного полного давления потока от р0 / ра = 80 до р0 / ра = 112 (при этом место отрыва перемещается через эти точки) показаны на рис. 8а. Когда поток отрывается от стенки в районе измерительной точки, давление в этой точке рм1 пульсирует от рк / ра = 0,25

(значение, соответствующее безотрывному течению) до рм1 /ра = 0,6.0,7 (значения,

соответствующие как дозвуковому течению при типичном отрыве потока, так и сверхзвуковому течению с низким числом Маха при нетипичном отрыве потока). Соответственно видно, что при увеличении р0 /ра при р0 /ра-98.102 поток отрывается до

точки 6А, при р0 /ра-102.106 поток присоединен к стенке в точке 6А и отрывается в точке

5 А, а при р0/ра > 106 поток полностью присоединен к стенке в точках 5 А и 6А и отрывается

в точке 4А. При уменьшении р0 / ра место отрыва потока проходит через эти точки в

обратном порядке, но при величинах р0 / ра на Лр0 / ра =2.3 больше, чем при увеличении

р0/ ра.

■— —

ЗБ

6А а

5А

В)

1А|Ж| 1Е 1Л,"

-п—

^ /

Ч- Ч Л И

0 2 й> 5 0 7 0 10 30 12 50 15 00 17 50 20 Част 00 22 50 25 00 27 50 30 30 32 50 35 00 37 50 «100

Г)

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

Показанные на рис.8, 9 суммарный и спектральный уровни, когерентность и фаза пульсаций давления определялись следующим образом [11]. Если пульсации давления р'(х) на интервале Тп их измерения имеют автокорреляционную функцию

1 п

Ярр (т, Тп) = —I )р'( _ т)& , то спектральная плотность этих пульсаций

Тп 0

Т

1 Г _■

3 (ю, Тп ) =-I Ярр (т)е_ю^т, где ю=2л/" - круговая частота. Тогда среднеквадратичный

2пТп 0

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

Ci =

си

|S(ю)^ю =yjRpp (0) и с1Гц (/) = ^S(/) в полосе шириной А/ = 1 Гц. Суммарный

уровень пульсаций давления (в диапазоне частот от 50 Гц до 25 кГц) LI = 20 • log-1, где cI -

C0

среднеквадратическое значение пульсаций давления (в Па), а с0 = 2•Ю-5 Па - порог

слышимости человеческого уха, а L Гц = 20 • log

W1 Гц

V C у

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

давления в полосе 1 Гц. Если пульсации давления р, , р ^ в точках ■ и ] их измерения в интервале времени Тп имеют пространственно-временную взаимнокорреляционную

Т

1

функцию (т) = — |р'(^)р' (^_т)&, то их взаимная спектральная плотность

Тп 0

1 Тп . (/)|

3 (ю) =-Г Я (т)в~шЯт, а их функция когерентности у2 ( /) = I , где / = -

2^Тп Г ^ ^ (/(/) ' 2п

ю

Я (т)е '"'Ят, а их функция когерентности / ) = ,—— ' '' , где _/ =

п0

- частота.

Из рис. 8б видно, что при безотрывном сверхзвуковом обтекании стенки ЬЕ =134 дБ, а

отрыв потока увеличивает ЬЕ до 175 дБ, т.е. сгЕ увеличивается приблизительно в 100 раз.

Спектральный уровень пульсаций давления Ь1Гц в диаметрально противоположных точках

1А и 1Б показан на рис. 8в. Функция когерентности, осредненная по 50...100 реализациям и разность фаз пульсаций давления представлены на рис. 8г. Видно, что Ь1Гц имеет максимум

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

Результаты измерения пульсаций давления в эксперименте с коническим насадком с 8 н =139,2 показаны на рис. 9. Видно, что конический насадок существенно изменил влияние

р0 / ра на отрыв потока в сопле: при увеличении р0 / ра поток отрывается в точке 6А при р0 / ра-82, а при ро / ра -95 место отрыва мгновенно сдвигается от точки 6А к точке 1А. Течение мгновенно становится безотрывным и сверхзвуковым в точках 6 А, 5 А, 4 А, 3А, при этом статическое давление в точке 1А уменьшается до рк/ра ~ 0,65 .

0.00 5.00 10.00 15.00 20.00 25.00 50.00 35.00 -Ю.00 45.00 50.00 55.00 60.00

Время, с

Л/

Л ■ V

% /V

1 \

V \

V*

* »»■] А и

1?П- ч V Г • Г 1

в)

0 250 500 750 1000 1250 1500 1750 2000! 2250 2500 2750 3000 .3250 3500 3750 4000

Частота

Рис. 9. Пульсации давления (а, б), измеренные в эксперименте с коническим насадком с 8 н =139,2 и их спектральный уровень (в)

При уменьшении р0 / ра поток отрывается от стенки в точке 2А при р0 / ра -90.100, в точке 3А - при р0 /ра-87.90, в точке 4А - при р0 /ра-80.85, в точке 5А - при р0 /ра-70.80, в точке 6А - при р0 /ра-65.70, т.е. обратного резкого перестроения потока не происходит.

Рис. 10. Масляные картины в отрывной зоне модельного сопла при типичном (а) и нетипичном (б) отрывах

При увеличении р0 / ра до ~93 поток отрывается в точке 5А так же, как и в

эксперименте без насадка. Спектральные уровни пульсаций давления в точках 1А, 1Б и 5А имеют максимум на частотах 550.600 Гц, пульсации хорошо коррелируют друг с другом с нулевой разностью фаз. При уменьшении р0 / ра амплитуда пульсаций давления в точке 1А

увеличивается, их максимальные значения превышают р^ра ~ 1,5, а среднее значение давления составляет р^ра ~ 1,0...1,2. При этом давления в точках 1А и 4А пульсируют синфазно, но максимум в районе 550 Гц отсутствует. Когда р0 /ра падает до 70, статическое давление в точке 1А увеличивается до р^ра «1,5...1,6, а в точке 2А уровень пульсаций начинает расти. При этом пульсации в точке 2А хорошо коррелируют с нулевой разницей фаз с пульсациями в точках 1А, 3А, 4А на частотах около 200 Гц. Начиная с р0 / ра ~64,

картина течения около точек 1А и 2А существенно меняется. Давления в этих точках начинают пульсировать в противофазе, а их когерентность уменьшается до 0,2.0,4 во всем диапазоне частот. Следовательно, линия растекания присоединяющегося к стенке потока находится между этими точками, колебания этой линии повышают давление в одной из этих точек и понижают в другой. Таким образом, увеличение статического давления до р^ра ~ 1,5...1,6, рост ЬГц до 152 дБ и ЬЕ до 178,5 дБ (рис. 9б) в точке 1А являются

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

Таким образом, спектральный анализ пульсаций давления показывает, что суммарный спектральный уровень пульсаций в отрывной зоне при нетипичном отрыве на 10 дб (т.е. более, чем в 3 раза по амплитуде) больше, чем при типичном отрыве. При этом, относительно небольшая отрывная присоединяющимся к стенке потоком.

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

Теневая картина ударно-волновой структуры струи данного сопла, сфотографированная на стенде БМР ФГУП «Центр Келдыша» показана на рис.11. Она аналогична описанной в [3-10] ударно-волновой структуре при нетипичном отрыве потока, хотя сопло спрофилировано методом характеристик и нетипичный отрыв потока при этом не возникает.

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

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

зона при нетипичном отрыве замкнута

Рис. 11. Теневая картина ударно-волновой структуры в истекающей из модельного сопла струе (срез сопла слева)

Рекомендовано к публикации программным комитетом XVI Международного симпозиума «Динамические и технологические проблемы механики конструкций и сплошных сред» им. А.Г. Горшкова

Библиографический список

1. Н.Б. Пономарев, И.Г. Лозино-Лозинская, А.Л. Воинов. Расчетные и экспериментальные методы определения удельного импульса тяги нового двигателя. - Глава 1 в книге «Рабочие процессы в жидкостном ракетном двигателе и их моделирование», под ред. академика РАН А.С. Коротеева. М.: Машиностроение, Москва, 2008, 512 с.

2. Nave, L.H., Coffey, G.A. Sea Level Side Loads in High-Area-Ratio Rocket Engines. AIAA Paper 73-1284, 1973.

3. Chen, C.L., Chakravarthy, S.R., Hung, C.M. Numerical Investigation of Separated Nozzle Flows. AIAA Journal, Vol.32, No.9, 1994, pp. 1836-1843.

4. Frey, M., Hagemann, G. Status of Flow Separation Prediction in Rocket Nozzles. AIAA Paper 98-3619, 1998.

5. Frey, M., Hagemann, G. Flow Separation and Side-Loads in Rocket Nozzles. AIAA Paper 99-2815, 1999.

6. Preuss, A. An analytical Approach for the Flowfield Analysis of Overexpanded Rocket Nozzles. International Astronautical Congress, Paper IAF-00-W.2.08, 2000.

7. Hagemann, G., Terhardt, M., Frey, M., Reijasse, P., Onofri M., Nasuti, F., Ostlund, J. Flow Separation and Side-Loads in Rocket Nozzles. Proceedings of the 4th International Symposium on Liquid Space Propulsion, DLR Lampoldshausen, March 12-15, 2000.

8. Hagemann, G., Alting, J., Preclik, D. Scalability for Rocket Nozzles Flows Based on Subscale and Full-Scale Testing. Journal of Propulsion and Power, Vol.19, No.3, 2003, pp. 321331.

9. Nasuti, F., Onofri, M. Flow Structures and Separation in Overexpanded Rocket Nozzles. Proceedings of the European Conference for Aerospace Sciences (EUCASS), 2006.

10. Reijasse, P., Poutrel, R. Flow Separation Regimes Induced by Cap-Shock in Over-Expanded Optimized Propulsive Nozzles. Proceedings of the European Conference for Aerospace Sciences (EUCASS), 2006.

11. А.Н. Антонов, В.М. Купцов, В.В. Комаров. Пульсации давления при струйных и отрывных течениях. - М.: Машиностроение, 1990, 272 с.

Сведения об авторах

Комаров Владимир Всеволодович, главный метролог ФГУП «Центр Келдыша», к.т.н., тел.: +7 495 459-93-91.

Пономарев Александр Александрович, младший научный сотрудник ФГУП «Центр Келдыша», аспирант Московского физико-технического института (государственного университета), тел.:+7 916 577-82-42; e-mail: ponomar aa@mail.ru.

Пономарев Николай Борисович, начальник сектора ФГУП «Центр Келдыша», к.т.н., й тел.: +7 495 456-00-56, e-mail: ponomarevnb@mail.ru.