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

Тип работы:
Реферат
Предмет:
ТЕХНИЧЕСКИЕ НАУКИ


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 532. 529+629
О ЗАДАЧЕ ПРОЕКТИРОВАНИЯ ИДЕАЛЬНОГО ДИФФУЗОРА ДЛЯ СЖАТИЯ СВЕРХЗВУКОВОГО ПОТОКА
1Усков В.Н., 2Булат П.В.
1 Балтийский государственный технический университет «Военмех» им. Д. Ф. Устинова,
Санкт-Петербург-
2ЗАО «Элекард наноДевайсез», Санкт-Петербург, e-mail: kolinti@mail. ru
В воздушно-реактивных двигателях диффузоры предназначены для повышения статического давления газа за счет уменьшения его кинетической энергии в процессе торможения. Рассмотрена задача изо-энтропического сжатия сверхзвукового потока. Если задать форму вогнутой поверхности в соответствии с уравнением линии тока в плоской волне Прандтля-Майера, то при натекании на нее сверхзвукового потока образуется центрированная волна сжатия, характеристики которой пересекаются в одной точке. В этой точке происходит образование ударно-волновой структуры с главным скачком уплотнения и с отраженным газодинамическим разрывом. Для проектирования сверхзвуковых диффузоров подходят волны сжатия с отраженным разрывом скачком уплотнения, т.к. они обладают свойством структурной устойчивости. Это обеспечивает требуемый запас саморегулирования диффузора. Проведенный анализ свойств центрированной волны сжатия дополняет теорию стационарных газодинамических разрывов. Определены границы областей существования ударно-волновых структур, оптимальных для проектирования сверхзвуковых диффузоров.
Ключевые слова: диффузор, изоэнтропическое сжатие, центрированная волна сжатия, ударно-волновая структура, поляра сжатия
ON THE PROBLEMS OF DESIGNING DIFFUSERS IDEAL COMPRESSION SUPERSONIC FLOW
1Uskov V.N., 2Bulat P.V.
1Baltiysky State Technical University «VOENMEKH» n.a. D.F. Ustinov, St. Petersburg-
2JSC «Elecard nanoDevices», St. Petersburg, e-mail: kolinti@mail. ru
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.
Keywords: diffuser isentropic compression, centered wave of compression, shock-wave structure, the polar compression
В воздушно-реактивных двигателях диффузоры предназначены для повышения статического давления газа за счет уменьшения его кинетической энергии в процессе торможения. Диффузоры воздушно-реактивных двигателей (ВРД) по числу Маха полета Мн подразделяются на дозвуковые Мн & lt- 1, малых сверхзвуковых скоростей или трансзвуковые Мн & lt- 1,5 и сверхзвуковые М, & gt- 1,5. Считается, что эффектив-
ность диффузора тем выше, чем больше степень повышения давления п в нем, ко -торая определяется, как отношение давления за замыкающим скачком уплотнения (3 на рис. 1) к давлению Р в невозмущенном потоке, а также выше степень восстановления полного давления V = Р0/Р0Н, где Р0 -давление на замыкающем скачком уплотнении, Р0Н — полное давление в невозмущенном потоке.
а) Классический трехскачковый диффузор
б) Изоэнтропический диффузор
Рис. 1. Сравнение многоскачкового и изоэнтропического воздухозаборника внешнего сжатия: 1, 2 — косые скачки уплотнения, 3 — замыкающий прямой скачок уплотнения, ЕГ — критическое сечение (горло), Д-Д — прямой скачок уплотнения, МН — число Маха невозмущенного течения
Потери полного давления на прямом скачке уплотнения, малые при М незначительно превосходящем единицу, резко возрастают с увеличением Мн. Уже при Мн = 2, коэффициент восстановления V = 0,72. За счет таких потерь полного давления тяга ВРД с односкачковым диффузором уменьшилась бы, примерно, на 40% по сравнению с тягой ВРД при изоэнтропическом сжатии [2].
Потери полного давления в сверхзвуковом диффузоре снижаются при замене сильного прямого скачка уплотнения системой более слабых косых скачков (1 и 2 на рис. 1), за которыми скорость остается сверхзвуковой, и замыкающим слабым прямым скачком, переводящим поток в дозвуковой.
При течении газа со сверхзвуковой скоростью в трактах (например, в воздухозаборнике) с криволинейной образующей, не имеющей угловых точек, могут возникать простые (ПВС) или центрированные волны сжатия (ЦВС). ЦВС позволяет сжать газ без потерь полного давления, что делает ее привлекательной для использования в технических устройствах вместо систем скачков уплотнения (рис. 1б). Внутри ПВС могут возникать висячие скачки уплотнения, но система ПВС + скачок сжимает газ с меньшими потерями, чем система скачков уплотнения [4].
Анализ центрированной изоэнтропической волны сжатия
Если задать форму вогнутой поверхности в соответствии с уравнением линии
тока в плоской волне Прандтля-Майера [5], то при натекании на нее сверхзвукового потока волны сжатия (характеристики) ЦВС юо пересекаются в одной точке, А (рис. 2). Происходит образование ударно-волновой структуры (УВС) с главным скачком конечной интенсивности о и с отраженным газодинамическим разрывом Я, который может быть скачком уплотнения, центрированной волной разрежения или слабым разрывом второго порядка (разрывной характеристикой, на которой терпят разрыв не значения газодинамических переменных, а их производные).
Рис. 2. Центрированная волна сжатия
Известно, что внутри волны сжатия параметры течения описываются решением Прандтля-Майера для плоской центрированной волны
& lt-«1+А=ю~+й~>-
где ю — функция Прандтля-Майера- - угол наклона вектора скорости.
Тогда, введя понятие интенсивности волны сжатия J = Р,/Р, можно записать:
со -со
у-1
1 + У-^М2
т-1 J у -1
и со 1
у-1
½
= *!-*¦
В совокупности с параметрами задачи М, у переменные J, в образуют следующие функциональные зависимости: Изоэнтропическая волна
Р = СО
у-1
1 + ^-Лл/2
Скачек уплотнения
-Гт~'-Г
т____
. /+є
ЧІ/2
ЇГІ & quot- ½'-
3^ 1 1
-?)(У- 1)
где
^={1-е)М2-е- е =
(Ут+Є)-(1-є)(/-1)' у-1
у+1
Кривые, описываемые этими уравнениями, будем называть полярой сжатия и ударной полярой соответственно. Точки
на поляре сжатия показывают отношение давлений п за ЦВС к давлению в невозмущенном потоке и угол разворота потока в центре волны сжатия. В начале координат (А = 0, Р = 0) эти кривые имеют порядок касания не ниже второго. Это свойство несложно выразить следующим образом:
р®=р®+др®,
где
0)_
Лі = Л2 = 0- _(м2−1)
(П,
'-,=1 №
Разница в значениях высших производных (/ & gt- 2) этих кривых при J = 1 зависит от А/. Опуская элементарные выкладки, можно записать для А/ при / = 3, 4
A3 =
|(м2-і)_1-і)2 -4 є 4(l+e)2 '
A. =
A В
M M M M2−1
где
A = - 2 a
& quot-14 11 '-
— + -y-l
У f
(m2- l) (m2- l) Y f
+ G,
D-1-
8
/1715_53_^7_Л
iff
E=X-
8
725_19
iff
При у = const зависимость A (M) немонотонная. Она имеет корни при значениях числа Маха, равных
& quot-*=Шз"т±(тМ)И)Г
и минимум, при М = V2 для любых у.
При М = Mj- поляры (сжатия и ударная) в начале координат имеют третий порядок касания. У функции A (M) нет действительных корней. Произведение А3Р0(1) обладает экстремумами при М = yfl и М = МА
М. =& gt-/2″
2(1-е)+,/3(1-е)
1−4е
½
А3Р0(1) стремится к ^ при М ^ 1 и к 0
при М ^ т. е. Корни произведения совпадают с корнями Д Из приведенных выше результатов следует, что поляра сжатия в начале координат может проходить как внутри
ударной поляры (М & lt- М^ и М & gt- М^), так и вне ее (М^& lt-М<-М^. От взаимного расположения поляр зависит тип отраженного разрыва (рис. 3).
При проектировании изоэнтропиче-ского диффузора необходимо так задавать геометрию вогнутой поверхности, чтобы отраженный разрыв Я (рис. 2) был скачком уплотнения. Чем меньше интенсивность этого скачка, тем выше коэффициент восстановления полного давления в диффузоре V = Р0/Р0Я. Возможна конфигурация, при которой скачок Я вырождается в разрывную характеристику (газодинамический разрыв второго порядка). Такую ЦВС назовем нейтральной.
„13 11 3 3
-4---------н----------------------
у Sy2 у3 у4
0−1
8
^ 1 17 7
23 н-і-г- н--
iff
Рис. З. ЦВС с отраженным разрывом скачком R и волной разрежения т
Нейтральная ЦВС
При некоторых значениях Ми у поляра сжатия и ударная поляра могут пересекаться. В точке пересечения имеет место равенство интенсивностей (/ш = Ja) скачка уплотнения о и ЦВС, а также углов поворота потока на этих разрывах. Следовательно, условие коллинеарности векторов скорости на тангенциальном разрыве выполняется при вырождении отраженного разрыва Я в характеристику. Назовем такую УВС нейтральной, интенсивность ЦВС в точке пересечения поляр обозначим J, а соответствующую ей кривую Jн (М, у| назовем нейтральной полярой (рис. 4). Нейтральная поляра имеет две ветви. Поскольку изоэн-тропическая волна сжатия не может затормозить поток до скорости, меньшей скорости звука, то область ее существования ограничена сверху звуковой линией J.
На рис. 4 нанесены график звуковой интенсивности скачка J и волны сжатия J, а также характерные точкиМН. (I = 1, 2, 3, 4),
В FUNDAMENTAL RESEARCH № 6, 2C12 В
& lt-„„, & lt-т“. Точки МН1 и МН4, соответствуют числам Маха, полученным выше аналитически. Мн2 представляет собой точку пересечения левой ветви графика ^(М) со звуковой линией поляры сжатия Jsю. М -точка складки нейтральной поляры ^(М). Точка & lt-„“ пересечения J и J соответству-
ет случаю, когда нейтральная поляра пересекает ударную поляру в звуковой точке. Две ветви нейтральной поляры и звуковая линия ограничивают с трех сторон область существования УВС с отраженным разрывом волной разрежения (на рисунке затенена).
Рис. 4. Нейтральная поляра:
Jш — интенсивность центрированной волны сжатия, 3ш — звуковая интенсивность центрированная волна сжатия, 3 — звуковая интенсивность скачка уплотнения, М14 — особые числа Маха, М^ - число Маха, ограничивающее область существования ударно-волновых структур с отраженным разрывом — скачком уплотнения
На рис. 4 и индексом & lt-т“ отмечена особая интенсивность ЦВС и соответствующее ей число Маха М При числе Маха, равном М ударные поляры, выпущенные из точек поляры сжатия, лежащих выше точки J не могут пересечься с главной ударной полярой (рис. 5).
Рис. 5. Поляры, соответствующие „особому числу Маха М^& gt-
Таким образом, для любого показателя адиабаты, начиная с числа Маха Мк, УВС с интенсивностью ЦВС большей 3, со-
¦ '-М~
держащие отраженный скачок уплотнения, существовать не могут. В области чисел Маха, больших М^ отраженный скачок -волна разрежения при любых значениях интенсивности ЦВС кроме Jн. Особые числа Маха, соответствующие им интенсивности ЦВС, и углы разворота потока приведены в таблице.
На рис. 6 приведена зависимость интенсивности волны сжатия J соответствующей нейтральной УВС, от числа Маха набегающего потока, а также показателя адиабаты. Здесь поверхность Jн (М у) изображена более темным цветом, а поверхность & lt-„„, соответствующая звуковой интенсивности скачка уплотнения, сделана полупрозрачной. На плоскости у-М изображены графики характерных чисел Маха М Видно, что при у = 1,67 зависимость МН4(у) стремится к да. М (проекция линии складки поверхности Jн (M, y)) сливается при у = 1,1 с МН4(у). Таким образом, при у = 1,1 имеет место особенность типа & lt-<-сборки“ [1].
Y МЯ МА МН1 МН2 МН3 МН4 М w J w р Q 9
1,1 1,3G2 1, ббб 1,3G2 1,4Вб — 1, ббб 2,599 б, 345 G, 69B
1,25 1,2б5 2 1,2б5 1,484 1,94G 2 2,99G B, 713 G, 669
1,4 1,245 2,539 1,245 1,478 2,23G 2,539 3,4B3 12,1B G, 635
1, б7 1,225 — 1,225 1,74б 2, B57 — 4,67G 22,7б —
Рис. 6. Нейтральная поляра при различных показателях адиабаты
Анализ областей существования УВС на плоскости поляр
Рассмотрим, как изменяется взаимное расположение ударной поляры и поляры сжатия по мере увеличения числа Маха невозмущенного течения.
При числах Маха М & lt- МН1 поляра сжатия целиком лежит внутри ударной поляры, следовательно, в данном диапазоне чисел Маха возможны только УВС с отраженным разрывом — волной разрежения (рис. 7).
Рис. S. МН1 & lt- М & lt- МН2
Рис. 7. М & lt- МИ1
Поляра сжатия при числах М, больших МН1, вблизи начала координат распространяется вне ударной поляры, пересекая ее выше. Здесь имеем отраженным разрывом либо скачок, либо волну разрежения, что определяется интенсивностью ЦВС (рис. 8). Эти два случая разделяет нейтральная конфигурация, т. е. при & gt- 3н отраженный разрыв — волна разрежения, а при интенсивностях волны меньших 3н отраженный разрыв — скачок уплотнения.
По мере увеличения М, значение Jн возрастает и при М=МН2 достигает звуковой интенсивности волны сжатия. При этом числе Маха поляра сжатия целиком лежит вне ударной поляры, касаясь ее „верхним краем“ (точка & lt-„“ на рис. 9).
Рис. 9. М = М“
В FUNDAMENTAL RESEARCH №б, 2G12 В
В диапазоне чисел Маха от МН2 до МН3 поляра сжатия проходит вне ударной поляры и не пересекает ее соответственно нейтральные УВС в этом диапазоне возникать не могут. При числе M, равном Мщ, поляра сжатия касается ударной поляры. На рис. 10 точка касания отмечена кружком. Кружок соответствует J и нейтральной УВС. При всех прочих интенсивностях отраженный разрыв — скачок уплотнения.
поляры, а интенсивность, соответствующая нижней точке, уменьшается, обращаясь в единицу при М = МН4. При больших М поляры имеют только одну точку пересечения. В газах с у = 1,67 — МН4 =, т. е. диапазон, в котором поляра сжатия пересекает ударную поляру изнутри и в одной точке, отсутствует.
Рис. 10. М = М,
При числах Маха, больших М и мень-
ших МН4, поляры пересекаются в двух точках (рис. 11). Между точками пересечения поляра сжатия проходит внутри ударной поляры. По мере увеличения М верхняя точка перемещается на сильную ветвь ударной
Рис. 11.
На рис. 12 представлена общая картина перестройки поляр сжатия и ударных поляр по мере изменения числа Маха.
Рис. 12. Трансформация поляр по мере изменения числа М
Заключение
Теоретически в рамках модели идеального газа, для сжатия потока без потерь полного давления можно использовать ЦВС звуковой интенсивности Jsm. Однако на практике такое решение реализовать непросто из-за влияния пограничного слоя.
Слишком большие градиенты давления приводят к его отрыву, а безотрывное течение обеспечивают только поверхности сжатия неприемлемо большой длины. Для проектирования сверхзвуковых диффузоров подходят ЦВС с отраженным разрывом скачком уплотнения, т.к. они облада-
ют свойством структурной устойчивости, т. е. при изменении внешних параметров в определенных пределах форма УВС не изменяется, а степень сжатия потока и коэффициент восстановления полного давления меняются плавно. Это обеспечивает требуемый запас саморегулирования диффузора. Область существования таких ЦВС ограничена числом Маха М^ = 3,483 для воздуха. Предельным случаем отраженного скачка является разрывная характеристика. Этому случаю соответствует нейтральная ЦВС, обеспечивающая наибольшие степень повышения давления и степень восстановления полного давления в диффузоре. В достаточно полной классификации УВС, образующихся в результате взаимодействия газодинамических разрывов, предложенной в [6], изоэнтропическая центрированная волна сжатия отсутствует. В работе [3] приведена теория распада центрированной волны сжатия Прадтля-Майера, проведенный анализ служит ее дополнением.
Список литературы
1. Арнольд В. И., Варченко А. Н., Гусейн-Заде С. М. Особенности дифференцируемых отображений. Том 1. Классификация критических точек, каустик и волновых фронтов. — М.: Наука, 1982. — 304 с.
2. Конструкция и проектирование авиационных газотурбинных двигателей: учебник для студентов вузов по специальности & lt-<-Авиационные двигатели и энергетические установки“ / С. А. Вьюнов, Ю. И. Гусев, А. В. Карпов и др.- под общ. ред. Д. В. Хронина. — М.: Машиностроение, 1989. — 368 с.: ил.
3. Омельченко А. В., Усков В. Н. Распад центрированной волны сжатия Прандтля-Майера в стационарном потоке газа // Прикладная механика и техническая физика. — 1998. -№ 3, т. 39. — С. 59−68.
4. Омельченко А. В., Усков В. Н., Чернышев М. В. Воздействие изоэнтропной волны на скачок уплотнения: аналитические решения и параметрический анализ // Струйные, отрывные и нестационарные течения: тезисы докладов XXII юбилейный семинар с международным участием. 22−25 июня 2010 г. — СПб.: СПБГУ, 2010. — 345 с.
5. Основы газовой динамики / под ред. Г. М. Эмонса. -М.: Ин. литература, 1963. — 702 с.
6. Интерференция стационарных газодинамических разрывов / В. Н. Усков и др. — Новосибирск: ВО „Наука“, 1995. — 180 с.
References
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.
Рецензенты:
Баранов И. В., д.т.н., профессор, заместитель директора Института Холода и Биотехнологий ФГБОУ & lt-<-Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики“, г. Санкт-Петербург-
Пеленко В. В., д.т.н., профессор, заместитель директора по учебной работе Института Холода и Биотехнологий ФГБОУ & lt-<-Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики», г. Санкт-Петербург.
Работа поступила в редакцию 09. 04. 2012.

ПоказатьСвернуть
Заполнить форму текущей работой