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

УДК 532.526.2

A.A. Балашов1, Г. Н. Дудин2,3 1ВАИС-Техника

2 Московский физико-технический институт (государственный университет) Центральный аэрогидродинамический институт имени профессора Н. Е. Жуковского

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

Исследовано течение в пограничном слое на пластине на режиме сильного взаимодействия при наличии массообмена на поверхности. Выведены и совместно решены системы нелинейных дифференциальных уравнений для вязкой и невязкой областей течения. Исследована зависимость собственного числа от температуры пластины, показателя адиабаты и скорости вдува/отсоса газа через поверхность. Установлено сильное влияние вдува на распространение возмущений вверх по потоку в пограничном слое.

Ключевые слова: гиперзвуковой поток, сильное взаимодействие, массообмен, температурный фактор, собственное число.

1. Введение

Исследование влияния массообмена через поверхность тела на поток представляет значительный интерес, так как при этом можно существенно уменьшить тепловые потоки и трение у поверхности, а также повлиять на распределение давления около тела. Исследование вдува газа по нормали к поверхности тела в сверхзвуковой поток показало, что для течений этого типа важную роль играет распространение возмущений вверх по потоку, и при решении необходимо учитывать краевое условие на заднем конце тела [1]. Заметим, что в отсутствие вдува происходит локализация распространения возмущений и реализуются сверхзвуковые течения со «свободным взаимодействием» [2,3].

При гиперзвуковом внешнем обтекании характер течения в пограничном слое значительно изменяется, так как плотность газа в нем уменьшается, а толщина вытеснения существенно возрастает. Для случая режима сильного взаимодействия было получено автомодельное («трехчетвертное») решение в [4,5]. В работе [6] с использованием приближенной формулы «касательного клина» впервые было показано, что на режиме сильного взаимодействия задача в окрестности передней кромки имеет однопараметрическое семейство решений, зависящих от некоторого собственного числа. В этом случае возмущения могут распространяться вверх по потоку на всю длину тела. В работе [7] вместо приближенной формулы «касательного клина» невязкая возмущенная область течения рассмотрена на основе гиперзвуковой теории малых возмущений. Были получены система дифференциальных уравнений, описывающая гиперзвуковое течение между скачком уплотнения и пограничным слоем, а также уравнения пограничного слоя и исследовано влияние температурного фактора на значение собственного числа. Исследование вдувания газа в гиперзвуковой поток на режиме сильного взаимодействия проведено в работе [8]. Как и в [6], для определения индуцированного давления использована формула «касательного клина» и показано, что в окрестности передней кромки решение не является единственным. В [8] сформулированы краевые задачи для определения автомодельного решения и для нахождения собственного числа. Однако собственное число, равное ^ ^ 0.23, определено только для предельного случая N ^ го. Физический смысл параметр а подобия N - квадрат отношения характерной толщины области невязкого течения, связанного с вдувом, к толщине пограничного слоя.

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

приближенной формулы «касательного клина» используется система уравнений, полученная с помощью гиперзвуковой теории малых возмущений. Проведены исследования функций течения и собственных чисел в зависимости от показателя адиабаты 7, температуры поверхности пластины Т** и скорости вдува газа V**.

2. Постановка задачи

Рассматривается гиперзвуковое обтекание вязким газом полубесконечной тонкой пластины под пулевым углом атаки при числе Маха невозмущенного потока Мж = и*/а^ ^ 1. В набегающем потоке и* - скорость, р^ - давление, р- плотность, а^ - скорость звука. Температура поверхности Т* считается заданной и постоянной. Рассматривается совершенный газ с отношением удельных теплоемкостей 7 = с*/с* и коэффициентом вязкости линейно зависящем от температуры. Вводится декартова система координат, начало которой расположено на передней кромке пластины. Компоненты вектора скорости и*, V* направлены вдоль осей ко ординат ж*, у*. Область течения за присоединенным головным скачком разделяется на область невязкого течения и вязкий пограничный слой (рис. 1).

у*,

Рис. 1. Гиперзвуковое обтекание плоской пластины при наличии массообмена на поверхности

Уравнение скачка: у* = Y*(х*), причем его наклон: dY* /dx* = 0(M—1). Внешняя граница пограничного слоя: у* = §*(х*). Детальная структура области течения за скачком может быть найдена в результате совместного решения двух систем уравнений в частных производных: одной гипербол и ческой, а другой параболической при выполнении условий сращивания на внешней границе пограничного слоя: непрерывности давления и равенства наклона линий тока. Режим взаимодействия определяется величиной параметра: % = М'3^/Re**^ [9], где Re^*^ = рЗ^оож*/- число Рейнольдса, вычисленное по текущему значению продольной координаты ж*. При % = 0(1) - реализуется умеренное взаимодействие, а при % ^ те - сильное взаимодействие. Однако, как показано в [9] для двумерных течений, режим сильного взаимодействия фактически реализуется при % ^ 3. В случае обтекания полубесконечной пластины при большом, но конечном значении М^ в окрестности передней кромки будет возникать область сильного взаимодействия, ниже по течению - режим умеренного взаимодействия, а еще ниже там, где % ^ 0(1) - режим слабого взаимодействия. Заметим, что в предельном случае, когда число Маха ^ те на всей пластине реализуется режим сильного вязко-невязкого взаимодействия.

Предполагается, что через проницаемую поверхность пластины по нормали к ней производится непрерывный распределенный вдув/отсос газа с заданной скоростью Vr*(x*, 0). Рассматривается случай, когда скорости вдува/отеоеа порядка скорости нормальной компоненты, возникающей в ламинарном пограничном слое, и, следовательно, отсоединения пограничного слоя не происходит [8].

3. Уравнения движения в невязкой области

Рассматривается течение в невязкой области 1 (рис. 1) между скачком и внешней границей пограничного слоя в предположении, что ^ 1 и реализуется режим вязко-невязкого взаимодействия. Считается, что невязкая область 1 описывается стационарными двумерными уравнениями Эйлера (неразрывности, импульса и энергии):

д(р*и*) + д (р*у*) = 0

дх* ду* '

* *ди* + * *ди* др*

р*и*---Ъ Р*у*^— = — ^—,

дх* ду* дх* /1\

ду* ду* др* ^

Р*и*~--V Р*у*^— = — ^—,

дх* ду* ду*

* * др * * * др * * / * др* * др* \

Р и дх* Р У ду * ^ V дх* У ду *)

Для оценки максимального продольного размера области умеренного и сильного вза-х*

" 1 " метр взаимодействия % = М'/= 1. Следовательно, характерный размер области

х* = М6ь'^/игде = Р*оо - кинематический коэффициент вязкости. На основании соотношений Рэнкина-Гюгонио при больших числах Маха получаем, что поперечный размер невязкой области 1 имеет порядок у* ~ х*/М00, а нормальная компонента скорости -у* ~ и*/Мсо. Далее вводятся безразмерные переменные [7]:

М6 и* М5 и* и *

М М и*

* хСЮ" С© * * Г Г* * ^ ОС * * * * /гч\

х = —и~*— х, у =-и-*— У, и = и*и, у = — у, р = рр, р = рр.

Поверхность скачка уплотнения в безразмерном виде задаётся в форме у = У3(х). В безразмерных переменных (2) система уравнений Эйлера (1) приводится к виду

д(Ри) + д (Ру) = 0

д х д д и д и 1 д

дх ^ ду МЦ^ дх1 /о\

ду ду 1 др ^ '

ри— + ру— = — —, дх ду 7 ду

др др ( др др\ ри— + ру— = 7р[и— +у^). д х д д х д

Для определения течения в невязкой области удобно ввести функцию тока ф, удовлетворяющую уравнению неразрывности:

дф дф

ду =ри,дх = —ру.

Учитывая (4) и переходя к переменным Мизеса: (х, у) ^ (х,ф), где х = х, можно систему (3) привести к виду

д и 1 д 1 д 1 и— = —тг^тт^--+ у

дх р дх Мдф М2,7'

ду др 1

дх дф 7'

д у 1 д

дф 7рр дх

При преобразовании уравнения энергии в (3) использовано уравнение неразрывности и учтено, что в гиперзвуковом приближении в невязкой области 1 за скачком и = 1+0(М-'). После введения функции энтропии Е = р/р1 система уравнений (5), описывающая течение в невязкой области, приводится к виду

дь др

^дх дф'

ду Е 7 др

дф ^

В соответствии с [4,5] скачок уплотнения рассматривается как известное автомодельное «трёхчетвертное» решение плюс возмущение п-то порядка:

3

Г, = Ах 3

где А, а;п - некоторые константы, которые зависят от граничных условий, п - собственное значение, которое определяет порядок возмущения поверхности скачка уплотнения. Соотношения Рэнкина-Гюгонио, преобразованные в соответствии с (7), учитывая преобразования Мизеса, дают следующую систему для нахождения параметров потока за скачком:

V

Р =

и = 1, 2 йУ3

7 + 1 йх ' 27 ( йЬ\2

7 + 1 V

7 + 1 "=—

Выражение для энтропии Е с учетом (7), (8) и того факта, что на скачке уплотнения ф = У3, можно записать в виде

Совершим ещё одну замену переменных (х,ф) ^ гДе С = V = Ф/^в- В этих

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

2 НУ / \

^ = {КМ + ап(пУп^)),

Р = )>(*) + ап^РпЫ).

Переходя в системе (6) к новым переменным (£,г)), учитывая (7), (9), (10) и собирая в уравнениях члены одного порядка, получаем систему из четырех нелинейных дифференциальных уравнений, описывающих течение в невязкой области (рис. 1):

Р0(v) -rjV0(rj) - ^Vo(r]) = 0,

4 / 1 \ 16

Р'п(ч) - ^(rj) + - jjVnirj) + —п(п + 1)VQ(rj) = 0,

Vi(ri) = 7 ~ 1 (\Ро(п) + г]Р'(г])),

7

™ - ^Р^ЫП ^ШРМ+ (П.

+2 '7———Г—Г-(п - 8

3 7 - дРМ = ъ(п+1)Ь* -

Ро 7

16 - - 1 / 3 \ / 1\ г] 3~f

9 7<П + *){" + 4) р 1...

Р07

Граничные условия для системы (11) получаются из соотношений (8) на скачке уплотнения при 7] = 1:

Ро

4. Уравнения движения в области пограничного слоя

Для пограничного слоя в области 2 (рис. 1) используется система уравнений из [10] в безразмерных переменных (2):

д . . д , . дХ (PU) + ди(PV) = 0,

du du 1 ap 1 д (q^u \

дх ду -M^pdx М^рдуУ ду)'

uf + V& - P—T = ^ + )•

дх ди p 7 ax M2p \ди/ M4,рди\ ди)

Где т _ безразмерная температура, С - константа в законе вязкости Чепмана: Р*/Рж = СТ*/Т,* = СТ. В системе (13) число Прандтля полагается единицей: Рг = р*с**/Л* = 1, где р* - динамический коэффициент вязкости, с* - удельная теплоёмкость при постоянном давлении, Л* - коэффициент теплопроводности. Для дальнейшего преобразования системы (13) вводится функция 5 как аналог энтальпии:

Т = I—1 М^

Производится переход к новой функции тока: Ф = М^ф, а также к перемену

ной Дородницына: г = М^ / рйу. После применения преобразования Дородницына:

о

(х, у) ^ (Х, г), где Х = х, к уравнениям импульса и энергии из (13) с учетом (14) получается следующая система:

дФ92Ф дФд2Ф_ 7 - 1 др (с „ 2ч , ^„д3Ф

:(S -u2)+Cp-

дг дХдг дХ дг2 27р дХ дг3 ' псч

дФ дБ_дФ дБ-Ср— дг дХ дХ дх дх2'

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

7 — 1 °о

сти. С учетом введенной энтальпии (14), толщина пограничного слоя 5 =-/(Б — и2) (г.

2р о

Так как на внешней границе пограничного слоя ье = (15/(1х, то выражение для нормальной компоненты скорости ье на границе областей

7 - 16, (1

5 -и2)йг)

2 (1х\р о

Давление на границе пограничного слоя определяется из решения (10) для невязкой области при т] —у 0. Так как йУ3/б£1 определяется из (7), то следовательно, можно искать давление в виде

р = Вх-2

где В и (п - некоторые константы, которые связаны с константами А и ап. С учетом этого решение системы (15) по аналогии с решением в невязкой области ищется в следующем

Ф = 2(ВС)2 ж4 (/о(0 + ЯпХп№)),

5 = 50

£ = _-_

2(ВС) 2ж 4

Совершая переход в (15) к новым переменным (х, х) — (х, £), где х = х, с учетом (18), получается система из четырех уравнений на нулевые и п-е приближения энтальпии и функции тока, которые описывают течение в вязкой пристеночной области:

/0"(0 + /о(0Ш + 1-1 (5о(0 - /о2(0) = 0,

/п"(0 + Л (0/п (О - (4п + ^^ )/5 (0/п(0+

+

х(Б(о - /о2(о) - /о''(о,

ад+лт (о = о,

Ж) + Л (№(£) - 4п/Ш)8п(0 + (4п +1) /п(М (О = -5'(С).

Для решения системы обыкновенных дифференциальных уравнений необходимо задать краевые условия на поверхности пластины и на внешней границе пограничного слоя. Учитывая порядок системы (19), необходимо иметь всего десять условий. При выходе на внешнюю границу пограничного слоя и — 1. Это условие на основании (14) и того, что ^ 1 равносильно 5 — 1. Считается, что возмущения энтальпии и скорости затухают при выходе на внешнюю границу. На пластинке поддерживается постоянная температура, соответствующая некоторому известному температурному фактору Б*.

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

и

С другой стороны, уравнения (4) после изложенной выше двойной замены переменных (х, у) — (х, г) — (х, £) приводятся к системе

Х и = Ж) + ЯпХп/'п(0,

М1 ру + ^ = -

После подстановки заданных граничных условий система приобретает вид

/о(0) = /¿(0) = 0,

4 1 /В\ 2 1

/о(0) + (4п + 1)дгажга/п(0) =--В 2 (1+ (х)х4.

7 - 1 Ьт\С /

Из второго уравнения данной системы следует вывод, что для автомодельности задачи необходимо вдувать газ со скоростью Ут = кх-4, где к - константа. Точно такой же результат был получен в более ранних работах [11] по исследованию автомодельных гиперзвуковых течений при наличии массообмена.

Следовательно, граничные условия имеют вид

А(0) = -7-г £ (В)2 '<0> = -7-Г (В)2 4ЙТГ • <2»>

Бп(0) = 0 Б(0) = ¡'(0) = /¿(0) = 0,

/о(те) = 50 = 1 /П(^) = Бп(х>) = 0. Согласно выражениям (20) массообмен на поверхности пластины при заданной ее температуре определяется величиной параметра к. Положительные значения к соответствуют вдуву в поток, а отрицательные - отсосу в пластину.

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

Краевые задачи (11), (12), (19) и (20) для обеих областей течения (рис. 1) не являются замкнутой системой. Для их решения необходимо найти дополнительное условие на собственное значение п. Используется сращивание па границе между областями нормальной компоненты скорости и давления. Из невязкого течения (10) при г] ^ 0

2 г1У / \

у(0) = — (0) + апХпУп(0)),

т = 7+Т (1Х НРо (0) + «пХпРп(0)).

Те же параметры со стороны вязкой области получаются из (16) и (17), учитывая распределение давления в пограничном слое:

1 те ч

- ¡(Б - и2),

2 Нх\р 0 р = Вх-2 (1 + Qпхп).

В результате сращивания получаем интегральное условие для определения собственного п

Ь(п) = Ро(0)(Уп(0) + 3(4п + 3)У0(0)) + ^ (Рп(0) + 2(4п + 3)Р0(0)) :

+ 1(4п + ЗД(0)) + ^ (Рп(0) + 3(4п + х(4п + 3)(1 + Л^)'

,,

те / \

а = /шо - /02(фе, 0

те / \ Л(п) = / (2/0(0/п(0 -Бп(£Ш6

0

Интегральное соотношение (21) замыкает краевую задачу, состоящую из систем дифференциальных уравнений (11), (19) и соответствующих краевых условий (12), (20) для описания обтекания гиперзвуковым потоком вязкого газа полубесконечной пластины под нулевым углом атаки при наличии массообмена на ее поверхности.

6. Результаты численных расчетов

Для совместного численного решения уравнений (11), (12), (19), (20) и (21) применялся метод Рунге Кутты 4-го порядка точности для невязкой области и прогонки для пограничного слоя. В результате определена зависимость собственного числа п(Бт ,7, к) от определяющих параметров задачи. На рис. 2 приведены зависимости собственного числа п от интенсивности вдува к для значений температурного фактора = 0.75,1,1.25,1.5 при 7 = 1.4. На рис. 3 представлены результаты расчетов для теплоизолированной пластины = 1 при значениях показателя адиабаты 7 = 1.2,1.4,1.66. На рис. 2 и рис. 3 для сравнения приведены результаты, полученные в рабоах [6,7] при исследовании обтекания пластины при отсутствии маесообмена.

Рис. 2. Зависимости , к) на режиме вдува при 7 = 1.4

Отличие полученных в данной работе результатов от [7] составляет около 3%, а от [6] около 6%. Заметим, что в [7] для описания невязкого течения использовалась гиперзвуковая теория малых возмущений, а в [6] приближенная формула «касательного клина». Учитывая существенно нелинейное поведение собственных чисел и зависимость от множества параметров, соответствие можно считать достаточно хорошим. Из полученных реп

7

и коэффициента массообмена к. Увеличение собственных чисел в этом случае означает ослабление влияния возмущений, так как У3, Уп, Рп, Бп, /п ~ хп, где х € (0,1]. Фактически

это означает, например, .локализацию возмущенного течения в окрестности задней кромки при рассмотрении пластины конечной длины. Во-вторых, при увеличении скорости вдува (к > 0.25) собственные числа достаточно быстро уменьшаются и престают зависеть как от температурного фактора, так и от показателя адиабаты и устремляются к 1. Это означает, что влияние краевых условий на течение в области, лежащей выше по потоку, должно быть достаточно сильным. Физическое объяснение этого явления связано со значительным увеличением толщины дозвуковой части пограничного слоя. Отметим, что в [6] для предельно сильного вдува на холодной пластине было получено собственное число, равное примерно 0.23.

На рис. 4 и рис. 5 приведены профили энтальпии Б(£) и возмущения энтальпии Бп(£) в пограничном слое от параметра массообмена к = —0.2, —0.1, 0, 0.05, 0.1 для случая 7 = 1.4,

= 0.75. Отсос газа приводит к существенному наполнению профиля энтальпии Б (О и росту теплового потока, а вдув к его уменьшению (рис. 4).

Рис. 4. Профиль энтальпии Бо(к) в пограничном слое при 7 = 1.4 и = 0.75

0.35 -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 £п О

Рис. 5. Профиль возмущения энтальпии , к) в пограничном слое при 7 = 1.4 и = 0.75

Следует отметить интересное поведение профилей возмущения энтальпии Бп(£) в зависимости от массообмена (рис. 5). При отсосе газа возмущения энтальпии существенно смещаются к поверхности пластины и располагаются в пристеночной части пограничного слоя, при этом толщина пограничного слоя уменьшается. При вдуве газа толщина пограничного слоя, естественно, увеличивается, максимум возмущения профиля энтальпии смещается от поверхности пластины, но теперь он располагается в средней части пограничного слоя. Если вспомнить полное выражение для энтальпии (18), то можно прийти к

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

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

и = Ш + Q„х•^m),

Qn=Ч Р(И<4п+4

Константа Qn зависит от величины донного давления: Ре = Ро + апРп, которое можно задавать на конце пластинки (х = 1). Деление, при котором ап = 0, будем называть автомодельным давлением Ра, так как в этом случае профиль скорости и(£) не меняется вдоль пластинки. Интерес представляют случаи задания донных давлений, отличных от автомодельного.

На рис. 6 представлены профили скорости в пограничном слое в сечениях х = 0 ^ 0.6, 0.8,1.0 вдоль пластинки для случая = 0.75, 7 = 1.4 Ре = 0.5Ра при вдуве газа с коэффициентом массообмена к = 0.1. В случае разгона потока (Ре < Ра) скорость в пограничном слое в окрестности задней кромки локально превышает скорость набегающего потока.

Рис. 6. Профили скорости в нескольких сечениях пластины по х в пограничном слое на режиме вдува (к = 0.1) в случае разгона потока (Ре = Ра/2); 7 = 1.4 5«, = 0.75

х

вдува (к = 0.1) в случае торможения потока (Ре = 1.2Ра); 7 = 1.4,5Ш = 0.75

На рис. 7 представлены профили скорости в пограничном слое в сечениях х = 0 + 0.6, 0.8,1.0 вдоль пластинки для случая = 0.75, 7 = 1.4, Ре = 1.2Ра при вдуве

к = 0.1

потока в окрестности поверхности пластины.

В рассматриваемом случае торможения потока (Ре > Ра) давление и коэффициент вду-ва газа подобраны таким образом, чтобы наблюдать на задней кромке начало явления отсоединения пограничного слоя (рис. 7). При дальнейшем увеличении донного давления отсоединение будет распространяться по пластинке в направлении передней кромки. Важно отметить, что, как только на задней кромке появляется отрыв, решение во всей рассматриваемой области 2 (рис. 1) перестает удовлетворять выведенной системе уравнений (11), (12), (19), (20) и (21) в силу появления невязкого подслоя на пластинке [6]. Следовательно, в данной постановке исследовать задачу при больших донных давлениях не представляется возможным.

7. Выводы

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

Работа выполнена при поддержке гранта РФФИ (проект №13-01-00202).

Литература

1. Матвеева, Н.С., Нейланд В.Я. Сильный вдув на теле конечной длины в сверхзвуковом потоке газа // Ученые записки НАГИ. — 1970. — Т. 1, № 5. — С. 13-22.

2. Нейланд В.Я. К теории отрыва ламинарного пограничного слоя в сверхзвуковом потоке // Изв. АН СССР. МЖГ. - 1969. - С. 53-57.

3. Stewartson К., Williams P.G. Self-induced separation // Proc. Roy. Soc. A. — 1969. — V. 312, N 1509. - P. 181-206.

4. Lees L. On the boundary-layer equations and their approximate solutions // JAS. — 1953. — N 20. - P. 143-145.

5. Stewartson K. On the motion of a flat plate at high speed in a viscous compressible fluid. II Steady motion // JAS. - 1955. - N 22. - P. 303.

6. Нейланд В.Я. Распространение возмущений вверх по течению при взаимодействии гиперзвукового потока с пограничным слоем // Изв. АН СССР. МЖГ. — 1970. — № 4. — С. 40-49.

7. Brown S.N., Stewartson К. A non-uniqueness of the hypersonic boundary layer // Q.J. Mechanics Appl. Math. - 1975. - V. 28, I. 1. - P. 75-90.

8. Нейланд В.Я. Вдувание газа в гиперзвуковой поток // Ученые записки ЦАГИ. — 1972. - Т. 3, № 6. - С. 29-40.

9. Хейз У.Д., Пробстин Р.Ф. Теория гиперзвуковых течений. — М.: ИЛ, 1962.

10. Stewartson К. The theory of Laminar Boundary Layer in Compressible Fluids. — Oxford: University Press, 1964.

11. Дудин Г.Н. Треугольные крылья в вязком гиперзвуковом потоке. — М.: МФТИ, 2011.

References

1. Matveeva, N.S., Neyland V.Y. Strong blowing on the body of finite length in a supersonic gas flow. CAI scientific notes. 1970. V. 1, N 5. P. 13-22.

2. Neyland, V.Y. For the theory of separation of a laminar boundary layer in a supersonic flow. Fluid Dynamics. 1969. N 4. P. 53-57.

3. Stewartson, K., Williams, P.G. Self-induced separation. Proc. Roy. Soc. A. 1969. V. 312, N 1509. P. 181-206.

4. Lees, L. On the boundary-layer equations and their approximate solutions. JAS. 1953. N 20. P. 143-145.

5. Stewartson, K. On the motion of a flat plate at high speed in a viscous compressible fluid. II Steady motion. JAS. 1955. N 22. P. 303.

6. Neyland, V.Y. The propagation of disturbances upstream in the interaction of hypersonic flow with boundary layer. Fluid Dynamics. 1970. N 4. P. 40-49.

7. Brown, S.N., Stewartson, K. A non-uniqueness of the hypersonic boundary layer. Q. J. Mechanics Appl. Math. 1975. V. 28, I. 1. P. 75-90.

8. Neyland, V.Y. The blowing gas in hypersonic flow. CAI scientific notes. 1972. V. 3, N 6. P. 29-40.

9. Hayes, W.D., Probstein, R.F. Hypersonic Flow Theory. Moscow.: Publishing house of foreign literature, 1962.

10. Stewartson, K. The theory of Laminar Boundary Layer in Compressible Fluids. Oxford: University Press, 1964.

11. Dudin, G.N. Triangular wings in viscous hypersonic flow. Moscow: MIPT, 2011.

Поступим в редакцию 17.02.2015.