CC BY
10
УДК 533.6.011: 532.529
Г. Д. Севостьянов
МЕТОД РАСЧЕТА ПАРАМЕТРОВ ПРИ РЕГУЛЯРНОМ ПЕРЕСЕЧЕНИИ КОСЫХ СКАЧКОВ
Дан алгоритм расчета параметров течения газа (на основе уравнений Чаплыгина) при регулярном взаимодействии косых скачков уплотнения. Рассмотренная схема течения устойчива по отношению к изгибу преломленных скачков, при этом «жидкий клин» отсутствует.
Уравнения Чаплыгина в переменных в, а [1] для плоского безвихревого установившегося течения идеального газа: ^ = —фа, = = К (а) гфо - потенциал скорости, в - угол наклона вектора скорости к оси х ф _ функция тока, К а функции Чаплыгина) в [2] приведены к нелинейной системе в дивергентной форме:
(Ь(п))1р = Уф, у^ = Пф, и = —са, V = св,
( 7 + 1 \3/(7—1}
Ь'(и) = к(и) = —К (а), с =(7 + 1)^ Ц-) , (1)
к(0) = 0, к'(0) = 1,
где 7 > 1 - отношение теплоемкостей газа.
Тогда на скачке уплотнения имеем ([/] = /+ — /— - разрыв ] на скачке)
[v] [L] (dy\2 [L]
S) = Ef; [v]2 = [L][u]. (2)
c [u] [v] ' V#/c [u] Система Чаплыгина в переменных v, u имеет вид
Vv = фи, Vu = k(u) ф.
+
см. [2
Эс = ^ + = , (3)
т. е. па ударной поляре имеем из (3) условие косой производной для ф: афи + b^v = 0.
Так как из уравнения ударной поляры
v+ — v- = (Ь+ — Ь—) (и+ — и—), Ь+ = Ь(и+), Ь— = Ь(и—)
имеем
1 М/ к+ + 1Ы
.¿п^,/ с = т 2\Щ\ к+ +[п]
ТО
2Ь = 3к+ + М, 2а = * (ЧЙ + ^М1' <4>
Для криволинейного скачка в однородном потоке Л. Крокко (Ь. Сгоссо) в 1937 г. ввел понятие «ежевидной» поляры. Для наклона «иголки» па поляре (образа линии тока 0 = ¿ф = фи ¿и + ф^ ¿V на плоскости ^,и)), который характеризует кривизну линии тока на скачке, имеем: (^/¿и)ИГ = — фи/ф^ = Ь/а.
Функции и и Ь представимы (см. [1]) интегралами пот = V2^таХ
< 1:
и = и(т\т*) = с [ (1 ^) ¿т,
Ь = Ь(т\т.) = с/ 21т+1^+1)Г^т > О, (5)
т*
7 -1 1
П = --- > 0, 0 < т < 1, в =-- > 0,
7 + 1 7 — 1
к = к(т )= —1 + (2в + 1)т к = к(т )= (1 — т )2в+1 '
Если 7 = 1 + 2/т, где т - число степеней свободы молекулы газа (без колебательных степеней свободы), то интегралы в (5) вычисляются аналитически. Из интеграла Вернули для числа Маха М имеем связь с т: (7 — 1)М2 = 2т/(1 — т).
Пусть в однородном сверхзвуковом потоке с М» > 1, V», в» = 0 (ось Ох || V») два косых скачка А+О (у > 0) и А—0 (у < 0) в точке 0(0,0) преломляются в виде косых скачков 0В+ (у > 0) и 0Б—. За ними в однородном дозвуковом потоке М2 < 1, в2, р2, У2. Углы падения обозначим > 0 > 0, углы преломления > 0 > 0 (индекс «+» для величин при у > 0 «—» для у < 0). Между скачками наклонные однородные сверхзвуковые потоки имеют М1± > 1, в1±.
Интенсивности падающих скачков = р»/р1± < 1; для преломленных скачков = р2/р1± > 1 (р - давление в газе).
Пусть дано М», 7, в2 (параметр несимметрии). Требуется найти остальные параметры. На плоскости ^,и) вершины двух малых ударных поляр (и1+ и и1—) лежат на большой поляре (и» > 0 V» =0). ^2,и2) - верхняя точка пересечения малых поляр, при этом в ней на-Ь/а
тока в О одинаковы). Тогда для т2, т1+ т1— имеют место три уравнения задано):
V2 = у/Ь(т1+\т2) и(т1+\т2) — у/Ь(т»\т1+) и(т»\т1+),
V2 = VхЬ(т»\т1 —) и(т»\т1 —) — VхЬ(т1 —\т2) и(т1 —\т2), (6)
—зк(т2) + уу = 3к(т2)+у (т'"\т2)
к(т2) + ЗУ(т1+\т2Г у (т1+ \т2) к(т2) + ЗУ(т1_\тгГ " (т'— \T21'
где введено обозначение: У (х\у) = Ь(х\у)/и(х\у).
После решения системы определяем другие параметры:
с#1+ = Vl+ = — \/Ь(т»\т1+) и(т»\т1+) < 0,
Vl— = \/Ь(т»\т1—) и(т»\т1—) > 0,
1 — А — cos($i±)
. V
tgw± = —=-
sin(|0i±|)
\/ Тпо
4/— cos $i± — , / — cos $2
. / V V
tg^± =
\ — sin |$i± | ± , — sin $2
V V
«± = (YTT— Y—1)—1,
í± = m2± sin2 (Ш±+101±|)—y—1 .
Y + 1 - + 1
При изменении параметра v2 (или угла $2) точка (v2, u2) на плоскости (v,u) вычерчивает «пояс» большой поляры (при v2 = 0 имеем точку Крокко для регулярного отражения косого скачка от стенки [2]), его концы - звуковые точки, в которых исчезает одна малая поляра. Сход с пояса приводит к теоретическому «жидкому клину».
182
В околозвуковой теории взаимодействия скачков (М ~ 1, и ~ 0) [3,4]:
и2
к(и) = и, £(и) = —, [у] = Т < и > [и];
2Ь = 3и++ < и >, 2а = ±(—^--Ь 3 < и > ,
+ \< и > )
и+ + и_ <и >= —2—, и уравнения (6) упрощаются.
Расчет параметров потока при у2 = 0 показывает (см. [2]) удовлетворительное приближение к экспериментальным значениям (не очень надежным) (в [5-7] дана библиография по данной задаче).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Кочин Н. Е., Кибель И. А., Розе Н. В. Теоретическая гидромеханика: в 2 ч. 4-е изд. М, : Гос. изд-во физ.-мат. лит., 1963. Ч. 2.
2. Севастьянов Г. Д. Метод расчета параметров регулярного отражения косого скачка от стенки// Математика. Механика: сб. науч. тр. Саратов: Изд-во Сарат. ун-та, 2008. Вып. 10. С. 140-143.
3. Севастьянов Г. Д. Регулярное несимметричное взаимодействие околозвуковых скачков// Математика. Механика: сб. науч. тр. Саратов: Изд-во Сарат. ун-та, 2002. Вып. 4. С. 219-222.
4. Севостьянов Г. Д. Регулярное несимметричное пересечение косых околозвуковых скачков// Математика. Механика: сб. науч. тр. Саратов: Изд-во Сарат. ун-та, 2003. Вып. 5. С. 184-186.
5. Основы газовой динамики / под ред. Г. Эммопса. М.: Изд-во иностр. лит., 1963.
6. Баженова Т. В., Гвоздева Л. Г. Нестационарные взаимодействия ударных волн. М.: Наука, 1977.
7. Ударные и детонационные волны. Методы исследования. 2-е изд., доп. и перс-раб. М.: Физматлит, 2004.
