CC BY
13
а WiE/qo, м W[4]E/qo, м W2E/qo, м W3E/qo, м
1° 15° 30° 1.381 ■ 105 1.226 ■ 105 0.840 ■ 105 1.224 ■ 105 0.8405 ■ 105 1.411 ■ 107 1.450 ■ 107 1.348 ■ 107 2.783 ■ 106 2.418 ■ 106 1.622 ■ 106
Сравнивая приведнные значения, можно заметить хорошее совпадение результатов расчета для первого случая со значениями, приведенными в [4].
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1, Григоренко Я.М., Крюков H.H. Решение задач теории пластин и оболочек е применением сплайн-функций (обзор) // Прикладная механика, 1995, Т. 31, вып 6, С, 3-27,
2, Годунов С. К. О численном решении краевых задач для систем обыкновенных дифференциальных уравнений // УМН, 1961, Т. XVI, Вып, 3 (99), С, 171-174,
3, Виноградов Ю.И., Виноградов А.Ю., Гусев Ю.А. Численный метод переноса краевых условий для жестких дифференциальных уравнений строительной механики, // Математическое моделирование, 2002, Т. 14, вып, 9, С, 3-8,
4, Крюков H.H. Расчет косоугольных и трапецеидальных пластин с помощью сплайн-функций // Прикладная механика, 1997, Т. 33, вып, 5, С, 77-81,
5, Доннел Л.Г. Балки, пластины и оболочки. Пер, с англ, / Под ред. Э.И, Григалюка, М.: Наука. 1982. 568 с.
УДК 533.6.011:632.529
Г.Д. Севостьянов
МЕТОД РАСЧЕТА ПАРАМЕТРОВ РЕГУЛЯРНОГО ОТРАЖЕНИЯ КОСОГО СКАЧКА ОТ СТЕНКИ
В статье предложен метод расчета параметров течения газа при регулярном отражении косого скачка от плоской стенки на основе уравнений Чаплыгина. Сравнение вычисленных параметров с экспериментальными удовлетворительно для числа Маха (1 ^ 3).
Уравнения Чаплыгина [1, §16, 17] для плоского безвихревого установившегося течения идеального газа
Рв = -фа, Ра = K(а)фв (1)
(где р — потенциал скорости, ф — функция тока, 9 — угол наклона вектора скорости к оси x К а — функции Чаплыгина) приводятся к нелинейной системе в дивергентной форме
k(u)uу = (L(u))ip = v-ф,vip = Пф,u = — са, v = c9, k(u) = —К(а), (2)
3
с =(Y + !)(7 — 1 ,k(0) = 0,k'(0) = 1,
где 7 > 1 — отношение теплоемкостей для газа.
Из (2) на скачке уплотнения имеем условия ([/] = /+ — /— - разрыв функции f на скачке) и уравнение ударной поляры
(Ар\ [V] [Ь] ([Ь] ]2 .
Ы с = —|«] = —'М'Ы с = й' М =[Ь]М' (3)
Система Чаплыгина (1) в переменных V, и имеет вид
^ = Фи, ^и = к(и)фу. (4)
Тогда вдоль скачка (индекс "+" для величин на задней стороне скачка
'£) с = ^ + (^ с фи = Ч § {фи + (^ с Ч ■ (5)
т.е. на ударной поляре имеем из (5) условие косой производной для ф
афи + Ьф = 0. (6)
Для криволинейного скачка в однородном потоке Л. Крокко (Ь. Сгоссо) в 1937 г. ввел понятие "ежевидной" поляры. Для наклона "иголки" на поляре в плоскости vu
/ = фи = Ь у а
Для плоской твердой стенки ((V = 0 Ь = 0 т-е- имеем точку Крокко Q поляры с условием
[Ь] [и]
2Ь = ^ + 3к+ = 0, и = ид < 0. (7)
Функции и и Ь можно представить интегралами но т = V2^ах < 1:
Г (1 — т)в( г Г — 1 + (2в + 1)т(^ 0
и = С ---(т, Ь = С -——-Г-Т-:—(т > 0,
Л. 2т , ]Тт 2т (1 — т )в+! > ,
7 — 1 1
т* = --- > 0,0 < т < 1,в =-7 > 0.
7 + 1 7 — 1
Пусть в однородном сверхзвуковом потоке (число Маха Ыж > 1) косой скачок АО под неизвестным углом падения и достигает плоской стенки (оси X в точке О и регулярно отражается от нее под неизвестным углом отражения и' в виде косого скачка ОВ. Обозначим числа Маха в области АО В и за ОВ через М1 > 1 и Мд < 1. По величине надо определить т1 и тд,
п г. I л ^ Р<ж .1 рд ,
< 0 и, и и интенсивности двух скачков 1 > £ = —, 1 < £ = — р1?
Р1 Р1
Рд - давления в трех областях). Из (7), (3), (8) имеем два уравнения длят1 и тд, записав условия на скачках АО и ОВ:
ГТ1 — 1 + (2в + 1) т 01 — (2в + 1) тд [Т1 (1 — т
(т = 3-^— ---(т,
Jtq 2т (1 - -f+1 (1 - -д)2в+1 Jtq 2-
•2=(?)'=г ^да - /; ^ =
= з1 - (2в +1)TQ iíT1 d-42
(1 - тд)2в+1 \Jtq 2-
2
Так как 7 = 1 +--(ш — число степеней свободы молекулы газа без учета
ш
ее колебаний), то все интегралы вычисляются аналитически. Углы падения
и отражения определяются по формулам = —, = \ —
Too V ^оо
1 - cos |0J cos |0i| -
tgu = ^-• |fl I , tgu =-T--, (9)
sin |0i| sin |0i|
полученным из рисунка двух ударных размерных поляр на декартовой плоскости годографа скорости (см. [1]) (вершиной малая поляра лежит на большой и пересекает ее ось в точке Q).
Интенсивности косых скачков АО и ОВ найдем по формулам соотношений на косом скачке (см. ill
£ =(-^гМ,sin2и - Vi) ,£' = -2YTM?sin2 (U + |ft|) - V1 (10)
-1
7 + 7 + 17 " 7 + 1""' l~1|/ 7 + 1"
Связь - и M определяется формулой, получаемой из интеграла Бернулли
(7 — 1) М2 = (11)
1 — т
Для околозвукового случая (Мто > 1) формулы упрощаются и переходят непрерывно в формулы околозвуковой теории [2, 3]: Ь(и) = и2/2, ид = —и1/7, Vд = 0 к(и) = и.
Результаты вычислений Е.А. Лунёва и автора (£) приведены в таблице для 7 = 7/5.
м Т1 и, ° и', ° £
1,02 1,172 0,170 0,166 79,8 84,15 0,991
1,20 0,223 0,199 0,161 60,64 71,48 0,901
1,50 0,310 0,244 0,152 49,29 60,69 0,745
2,00 0,444 0,307 0,138 41,92 48,61 0,522
2,50 0,555 0,356 0,116 39,27 39,30 0,363
3,00 0,642 0,395 0,098 38,65 31,13 0,255
Найденный зависимости и'(и), и(£) достаточно близки к полученным экспериментально ([4, рис. 8.6, а; 5, фиг. 15 на с. 461] во всем диапазоне чисел Маха (1 < Ыж < 3).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Кочин Н.Е., Кибель И.А., Розе Н.В. Теоретическая гидромеханика. Ч. II. Изд. 4-е. М.: ГИФМЛ. 1963.
2. Севостъянов Г. Д. Регулярное отражение околозвукового скачка от стенки. // Математика. Механика: Сб. науч. тр. Саратов: Изд-во, Сарат. ун-та, 2001. Вып. 3. С. 181-184.
3. Китанин В.А., Севостъянов Г.Д. Расчет регулярного отражения околозвукового скачка от плоской стенки // Математика. Механика: Сб. науч. тр. Саратов: Изд-во. Сарат. ун-та, 2002. Вып. 4. С. 187-190.
4. Баженова Т.В., Гвоздева Л. Т. Нестационарные взаимодействия ударных волн. М,: Наука. 1977.
5. Основы газовой динамики / Под ред. Г. Эммонса. М,: Изд-во. иностр. лит. 1963.
УДК 533.6
И.А. Чернов
НОВЫЕ РЕШЕНИЯ УРАВНЕНИЙ ПОГРАНИЧНОГО СЛОЯ БЕЗ ГРАДИЕНТА ДАВЛЕНИЯ
В статье представлен ряд новых автомодельных решений. Они описываются нелинейными уравнениями второго порядка, в коэффициенты которых входит параметр — показатель автомодельности. Решения являются формальными, их физическая интерпретация в теории вязко-невязкого взаимодействия, в частности, при рассмотрении вопросов об отрыве и о смешении потоков не обсуждается.
1. Плоско-параллельный случай.
Уравнения движения и неразрывности для компонент вектора скорости и,у) имеют вид {и - коэффициент вязкости}
ППХ + УПу = УПуу, их + Уу = 0.
С учетом уравнения неразрывности компоненты скорости выражаются через автомодельный представитель Ф (() функции тока (далее ( = у/х7
и (х, у) = х1-2^ (С) , V (х, у) = - [(1 - п) Ф (С) - <Ф' (С)] х-7.
