Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

ON THE PROBLEMS OF DESIGNING DIFFUSERS IDEAL COMPRESSION SUPERSONIC FLOW

Uskov V.N. 1 Bulat P.V. 2
1 Baltiysky State Technical University «VOENMEKH» n.a. D.F. Ustinov, St. Petersburg
2 JSC «Elecard nanoDevices», St. Petersburg
1222 KB
In air-breathing jet engine diffusers are designed to increase the static pressure of the gas. The problem of isentropic compression of a supersonic flow. If you specify the shape of the concave surface in accordance with the equation of the streamlines in a plane wave Prandtl-Meyer, then the leakage of her supersonic flow forms a centered compression wave, the characteristics of which intersect at one point. At this point, formed by shock-wave structure with the main shock wave and the reflected gasdynamic discontinuity. For the design of supersonic diffusers are suitable compression wave with the reflected shock wave discontinuity, as They have the property of structural stability. This ensures the required stock of self-regulation of the diffuser. The analysis of the properties of the centered compression wave theory complements the stationary gasdynamic discontinuities. The limits of regions of existence of shock-wave structures that are optimal for the design of supersonic diffusers.
diffuser isentropic compression
centered wave of compression
shock-wave structure
the polar compression
1. Arnol’d V.I., Varchenko A.N., Gusejn-Zade S.M. Oso­bennosti differenciruemyh otobrazhenij. Tom 1. Klassifikacija kriticheskih tochek, kaustik i volnovyh frontov. M.: Nauka, 1982, 304 р.
2. V’junov S.A., Gusev Ju.I., Karpov A.V. i dr. Konstrukcija i proektirovanie aviacionnyh gazoturbinnyh dvigatelej: Uchebnik dlja studentov vuzov po special’nosti «Aviacionnye dvigateli i jenergeticheskie ustanovki» / Pod obwej redakciej D.V. Hronina. M.:Mashinostroenie, 1989. 368 р.: il.
3. Omel’chenko A.V., Uskov V.N. Raspad centrirovannoj volny szhatija Prandtlja-Majera v stacionarnom potoke gaza. Prikladnaja mehanika i tehnicheskaja fizika, 1998, no. 3, t. 39, рр. 59–68.
4. Omel’chenko A.V., Uskov V.N., Chernyshev M.V. Vozdejstvie izojentropnoj volny na skachok uplotnenija: analiticheskie reshenija i parametricheskij analiz. Strujnye, otryvnye i nestacionarnye techenija: XXII jubilejnyj seminar s mezhdunarodnym uchastiem. 22-25 ijunja 2010 g. Sankt-Peterburg: tezisy dokladov/BGTU, SPBGU, 2010. 345 p.
5. Osnovy gazovoj dinamiki. Pod red. Jemonsa G.M. M.: In. literatura, 1963, 702 p.
6. Uskov V.N. i dr. Interferencija stacionarnyh gazodina­micheskih razryvov. Novosibirsk: VO «Nauka», 1995. 180 p.

В воздушно-реактивных двигателях диффузоры предназначены для повышения статического давления газа за счет уменьшения его кинетическои энергии в процессе торможения. Диффузоры воздушно-реактивных двигателей (ВРД) по числу Маха полета Мн подразделяются на дозвуковые Мн < 1, малых сверхзвуковых скоростей или трансзвуковые Мн < 1,5 и сверхзвуковые Мн > 1,5. Считается, что эффективность диффузора тем выше, чем больше степень повышения давления π в нем, которая определяется, как отношение давления за замыкающим скачком уплотнения (3 на рис. 1) к давлению Рн в невозмущенном потоке, а также выше степень восстановления полного давления ν = Р0/Р0Н, где Р0 - давление на замыкающем скачком уплотнении, Р0Н - полное давление в невозмущенном потоке.

а) Классический трехскачковый диффузор б) Изоэнтропический диффузор

Рис. 1. Сравнение многоскачкового и изоэнтропического воздухозаборника внешнего сжатия: 1, 2 - косые скачки уплотнения, 3 - замыкающий прямой скачок уплотнения, FГ - критическое сечение (горло), Д-Д - прямой скачок уплотнения, МН - число Маха невозмущенного теченияПотери полного давления на прямом скачке уплотнения, малые при Мн, незначительно превосходящем единицу, резко возрастают с увеличением Мн. Уже при Мн = 2, коэффициент восстановления ν = 0,72. За счет таких потерь полного давления тяга ВРД с односкачковым диффузором уменьшилась бы, примерно, на 40% по сравнению с тягой ВРД при изоэнтропическом сжатии [2].

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

При течении газа со сверхзвуковой скоростью в трактах (например, в воздухозаборнике) с криволинейной образующей, не имеющей угловых точек, могут возникать простые (ПВС) или центрированные волны сжатия (ЦВС). ЦВС позволяет сжать газ без потерь полного давления, что делает ее привлекательной для использования в технических устройствах вместо систем скачков уплотнения (рис. 1б). Внутри ПВС могут возникать висячие скачки уплотнения, но система ПВС + скачок сжимает газ с меньшими потерями, чем система скачков уплотнения [4].

Анализ центрированной изоэнтропической волны сжатия

Если задать форму вогнутой поверхности в соответствии с уравнением линии тока в плоской волне Прандтля-Майера [5], то при натекании на нее сверхзвукового потока волны сжатия (характеристики) ЦВС ωσ пересекаются в одной точке А (рис. 2). Происходит образование ударно-волновой структуры (УВС) с главным скачком конечной интенсивности σ и с отраженным газодинамическим разрывом R, который может быть скачком уплотнения, центрированной волной разрежения или слабым разрывом второго порядка (разрывной характеристикой, на которой терпят разрыв не значения газодинамических переменных, а их производные).

? 

Рис. 2. Центрированная волна сжатия

Известно, что внутри волны сжатия параметры течения описываются решением Прандтля-Майера для плоской центрированной волны

где ω - функция Прандтля-Майера; - угол наклона вектора скорости.

Тогда, введя понятие интенсивности волны сжатия Jω = P1/P, можно записать:

В совокупности с параметрами задачи М, γ переменные J, β образуют следующие функциональные зависимости:

Изоэнтропическая волна

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

?

где

Кривые, описываемые этими уравнениями, будем называть полярой сжатия и ударной полярой соответственно. Точки на поляре сжатия показывают отношение давлений π за ЦВС к давлению в невозмущенном потоке и угол разворота потока в центре волны сжатия. В начале координат (D = 0, b = 0) эти кривые имеют порядок касания не ниже второго. Это свойство несложно выразить следующим образом:

? 

где

Разница в значениях высших производных (i > 2) этих кривых при J = 1 зависит от Δi. Опуская элементарные выкладки, можно записать для Δi при i = 3, 4

?

?

где   

 

  

При γ = const зависимость D3(M) немонотонная. Она имеет корни при значениях числа Маха, равных

и минимум, при  для любых g.

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

 стремится к ∞ при M → 1 и к 0 при M → ¥. Корни произведения совпадают с корнями Δ3. Из приведенных выше результатов следует, что поляра сжатия в начале координат может проходить как внутри ударной поляры ( и ), так и вне ее (). От взаимного расположения поляр зависит тип отраженного разрыва (рис. 3).

При проектировании изоэнтропического диффузора необходимо так задавать геометрию вогнутой поверхности, чтобы отраженный разрыв R (рис. 2) был скачком уплотнения. Чем меньше интенсивность этого скачка, тем выше коэффициент восстановления полного давления в диффузоре ν = Р0/Р0Н. Возможна конфигурация, при которой скачок R вырождается в разрывную характеристику (газодинамический разрыв второго порядка). Такую ЦВС назовем нейтральной.

?

Рис. 3. ЦВС с отраженным разрывом скачком R и волной разрежения ω

Нейтральная ЦВС

При некоторых значениях M и g поляра сжатия и ударная поляра могут пересекаться. В точке пересечения имеет место равенство интенсивностей (Jw = Js) скачка уплотнения σ и ЦВС, а также углов поворота потока на этих разрывах. Следовательно, условие коллинеарности векторов скорости на тангенциальном разрыве выполняется при вырождении отраженного разрыва R в характеристику. Назовем такую УВС нейтральной, интенсивность ЦВС в точке пересечения поляр обозначим Jн, а соответствующую ей кривую Jн (М, g) назовем нейтральной полярой (рис. 4). Нейтральная поляра имеет две ветви. Поскольку изоэнтропическая волна сжатия не может затормозить поток до скорости, меньшей скорости звука, то область ее существования ограничена сверху звуковой линией Jsw.

На рис. 4 нанесены график звуковой интенсивности скачка Js и волны сжатия Jsw, а также характерные точки МНi. (I = 1, 2, 3, 4),  «s», «w». Точки МН1 и МН4, соответствуют числам Маха, полученным выше аналитически. Мн2 представляет собой точку пересечения левой ветви графика Jн(M) со звуковой линией поляры сжатия J. МН3 - точка складки нейтральной поляры Jн(М). Точка «s» пересечения Jн и Js соответствует случаю, когда нейтральная поляра пересекает ударную поляру в звуковой точке. Две ветви нейтральной поляры и звуковая линия ограничивают с трех сторон область существования УВС с отраженным разрывом волной разрежения (на рисунке затенена).

Рис. 4. Нейтральная поляра: Jω - интенсивность центрированной волны сжатия, J - звуковая интенсивность центрированная волна сжатия, Js - звуковая интенсивность скачка уплотнения, М1-4 - особые числа Маха, Мw - число Маха, ограничивающее область существования ударно-волновых структур с отраженным разрывом - скачком уплотнения

На рис. 4 и индексом «w» отмечена особая интенсивность ЦВС и соответствующее ей число Маха Мw. При числе Маха, равном Mw, ударные поляры, выпущенные из точек поляры сжатия, лежащих выше точки Jw, не могут пересечься с главной ударной полярой (рис. 5).

?

Рис. 5. Поляры, соответствующие «особому числу Маха Мw»

Таким образом, для любого показателя адиабаты, начиная с числа Маха Мw, УВС с интенсивностью ЦВС большей Jw, содержащие отраженный скачок уплотнения, существовать не могут. В области чисел Маха, больших Мw, отраженный скачок - волна разрежения при любых значениях интенсивности ЦВС кроме Jн. Особые числа Маха, соответствующие им интенсивности ЦВС, и углы разворота потока приведены в таблице.

На рис. 6 приведена зависимость интенсивности волны сжатия Jн, соответствующей нейтральной УВС, от числа Маха набегающего потока, а также показателя адиабаты. Здесь поверхность Jн(М, γ) изображена более темным цветом, а поверхность «s», соответствующая звуковой интенсивности скачка уплотнения, сделана полупрозрачной. На плоскости γ-M изображены графики характерных чисел Маха МНi. Видно, что при γ = 1,67 зависимость МН4(γ) стремится к ∞. МН3 (проекция линии складки поверхности Jн(M,γ)) сливается при γ = 1,1 с МН4(γ). Таким образом, при γ = 1,1 имеет место особенность типа «сборки» [1].

γ

МН1

МН2

МН3

МН4

Мw

Jw

?σω

1,1

1,302

1,666

1,302

1,486

-

1,666

2,599

6,345

0,698

1,25

1,265

2

1,265

1,484

1,940

2

2,990

8,713

0,669

1,4

1,245

2,539

1,245

1,478

2,230

2,539

3,483

12,18

0,635

1,67

1,225

-

1,225

1,746

2,857

-

4,670

22,76

-

?

Рис. 6. Нейтральная поляра при различных показателях адиабаты

Анализ областей существования УВС на плоскости поляр

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

При числах Маха М < МН1 поляра сжатия целиком лежит внутри ударной поляры, следовательно, в данном диапазоне чисел Маха возможны только УВС с отраженным разрывом - волной разрежения (рис. 7).

Рис. 7. М < МН1

Поляра сжатия при числах М, больших МН1, вблизи начала координат распространяется вне ударной поляры, пересекая ее выше. Здесь имеем отраженным разрывом либо скачок, либо волну разрежения, что определяется интенсивностью ЦВС (рис. 8). Эти два случая разделяет нейтральная конфигурация, т.е. при Jω > Jн отраженный разрыв - волна разрежения, а при интенсивностях волны меньших Jн отраженный разрыв - скачок уплотнения.

?

Рис. 8. МН1 < М < МН2

По мере увеличения М, значение Jн возрастает и при М = МН2 достигает звуковой интенсивности волны сжатия. При этом числе Маха поляра сжатия целиком лежит вне ударной поляры, касаясь ее «верхним краем» (точка «s» на рис. 9).

?

Рис. 9. М = МН2

В диапазоне чисел Маха от МН2 до МН3 поляра сжатия проходит вне ударной поляры и не пересекает ее соответственно нейтральные УВС в этом диапазоне возникать не могут. При числе M, равном МН3, поляра сжатия касается ударной поляры. На рис. 10 точка касания отмечена кружком. Кружок соответствует Jн и нейтральной УВС. При всех прочих интенсивностях отраженный разрыв - скачок уплотнения.

?

Рис. 10. М = МН3

При числах Маха, больших МН3, и меньших МН4, поляры пересекаются в двух точках (рис. 11). Между точками пересечения поляра сжатия проходит внутри ударной поляры. По мере увеличения М верхняя точка перемещается на сильную ветвь ударной поляры, а интенсивность, соответствующая нижней точке, уменьшается, обращаясь в единицу при М = МН4. При бóльших М поляры имеют только одну точку пересечения. В газах с g = 1,67 - МН4 = ¥, т.е. диапазон, в котором поляра сжатия пересекает ударную поляру изнутри и в одной точке, отсутствует.

?

Рис. 11.

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

Рис. 12. Трансформация поляр по мере изменения числа М

Заключение

Теоретически в рамках модели идеального газа, для сжатия потока без потерь полного давления можно использовать ЦВС звуковой интенсивности J. Однако на практике такое решение реализовать непросто из-за влияния пограничного слоя. Слишком большие градиенты давления приводят к его отрыву, а безотрывное течение обеспечивают только поверхности сжатия неприемлемо большой длины. Для проектирования сверхзвуковых диффузоров подходят ЦВС с отраженным разрывом скачком уплотнения, т.к. они обладают свойством структурной устойчивости, т.е. при изменении внешних параметров в определенных пределах форма УВС не изменяется, а степень сжатия потока и коэффициент восстановления полного давления меняются плавно. Это обеспечивает требуемый запас саморегулирования диффузора. Область существования таких ЦВС ограничена числом Маха Мw = 3,483 для воздуха. Предельным случаем отраженного скачка является разрывная характеристика. Этому случаю соответствует нейтральная ЦВС, обеспечивающая наибольшие степень повышения давления и степень восстановления полного давления в диффузоре. В достаточно полной классификации УВС, образующихся в результате взаимодействия газодинамических разрывов, предложенной в [6], изоэнтропическая центрированная волна сжатия отсутствует. В работе [3] приведена теория распада центрированной волны сжатия Прадтля-Майера, проведенный анализ служит ее дополнением.

Рецензенты:

  • Баранов И.В., д.т.н., профессор, заместитель директора Института Холода и Биотехнологий ФГБОУ «Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики», г. Санкт-Петербург;
  • Пеленко В.В., д.т.н., профессор, заместитель директора по учебной работе Института Холода и Биотехнологий ФГБОУ «Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики», г. Санкт-Петербург.

Работа поступила в редакцию 09.04.2012.