CC BY
69
МОДЕЛИРОВАНИЕ ПРОЦЕССОВ
УДК 533.6.011
В. Т. Калугин, П. А. Чернуха
ОСОБЕННОСТИ ОБТЕКАНИЯ ЛЕТАТЕЛЬНЫХ АППАРАТОВ ПРИ СТРУЙНОМ УПРАВЛЕНИИ В УСЛОВИЯХ ДВИЖЕНИЯ С МАЛЫМИ СВЕРХЗВУКОВЫМИ СКОРОСТЯМИ
В целях определения эффективности управления аэродинамическими характеристиками летательного аппарата со струйными органами управления при трансзвуковых скоростях обтекания были исследованы различные варианты выдува струй газа из сопел, удаленных от его кормовой части. Выявлено влияние условий вы-дува струй на формирование управляющих сил. Предложен метод расчета параметров течения и аэродинамических коэффициентов, проведен численный эксперимент, подтверждающий работоспособность метода на таких режимах обтекания.
При проектировании современных высокоманевренных летательных аппаратов (ЛА) в последнее время все чаще применяются активные (или газодинамические) методы управления, разновидностью которых являются струйные органы управления (СОУ), осуществляющие вертикальный выдув газа с боковой поверхности ЛА во внешний поток. Применение органов управления (ОУ) такого класса позволяет эффективно управлять процессом движения ЛА по каналу тангажа в широком диапазоне изменения скоростей движения, а инжекция нескольких струй в разнотяговом режиме позволяет дополнительно создавать аэродинамические силы и момент по каналу крена. Управляющие силы и момент при функционировании таких органов управления вызваны перераспределением давления по боковой поверхности ЛА перед и за отверстиями выдува, а также тягой самих струй. При сверхзвуковых скоростях движения инжекция газа использует в своей основе эффект образования локальных пространственных областей отрыва потока и связанное с ним формирование дополнительных управляющих сил. Принцип обтекания газоструйных ОУ на таком режиме достаточно хорошо изучен [1-3]. Однако, если обтекание ЛА происходит при скоростях, близких к скорости звука, в структуре течения и аэродинамических характеристиках проявляются особенности, которые необходимо учитывать при выборе конструктивных параметров ОУ. Такие режимы обтекания характерны при использовании импульсной струйной коррекции траектории движения управляемых мин, снарядов и других типов ЛА. Распределение параметров потока
по поверхности ЛА в этом случае отличается от сверхзвукового обтекания и вызвано формирующимся перед струями криволинейным скачком уплотнения, а также распределением области пониженного давления за счет расширения потоков газа за отверстиями инжекции. Кроме того, для блочных схем выдува наблюдается интерференция инжектируемых потоков в пространстве между соплами, а также за ними, что приводит к дополнительному усложнению структуры обтекания и существенному изменению аэродинамических характеристик.
Цель настоящей работы — экспериментальное определение особенностей обтекания тел в условиях струйной инжекции при малых сверхзвуковых скоростях обтекания; выявление физических моделей течения; определение эффективности использования струйной инжекции; математическое моделирование взаимодействия основного и инжектируемого потоков, проведение численных расчетов.
Результаты экспериментальных исследований. Для выявления особенностей управления ЛА с помощью струйных ОУ и определения эффективности их использования были исследованы различные схемы инжекции газа из сопел, размещенных на боковой поверхности тела вращения. Модель ЛА представляла собой тонкий затупленный конус с углом полураствора вк = 6,0°, имеющий четыре сменных блока с сопловыми вставками, которые обеспечивали выдув одной (Nj = 1) и двух (Nj = 2) (в дальнейшем дипольных) струй с различными углами £ наклона осей по отношению друг другу.
Ниже приведены исследуемые варианты а... г (рис. 1, б) расположения сопел инжекции на боковой поверхности конического тела (см. рис. 1, а), а в таблице приведены их основные геометрические параметры: Nj; = ^/йт — относительный диаметр сопла инжекции (¿т — диаметр миделевого сечения ЛА); £.
Инжекция струй газа осуществлялась через звуковые сопла, расположенные на боковой поверхности модели на расстоянии Ц = Ц/Ьк = = 0,9 от носка конуса (здесь Ьк — длина конического тела). Расстояние между соплами моделей, выполненных по дипольной схеме выдува,
Вариант а Вариант б
Рис. 1. Варианты выдува струй
4 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. "Машиностроение". 2007. № 4
Таблица
Основные параметры сопел инжекции
Модель Варианты расположения сопел
а б в г
1 2 2 2
dj 0,0517 0,0345 0,0345 0,0345
£ — 0° -60° 180°
соответствовало 1с = ¡с/Ьк = 0,073. Диаметр двух отверстий инжекции выбирался таким, чтобы и их суммарная площадь была близка к площади одиночного сопла.
Комплексные исследования по визуализации структур течения, весовые и дренажные эксперименты проводились на трансзвуковой аэродинамической установке НИИ "Геодезия", которая имела закрытую рабочую часть с размерами 0,6 х 0,6 х 0,9 м3 и возможный диапазон изменения чисел М^ = 0,7... 1,7. В экспериментах число Маха невозмущенного потока соответствовало М^= 1,2, относительная интенсивность выдува струй принята в виде 7р0 = р^/рх = 0 ... 140, где Ро — полное давление в форкамере сопел выдува, р^ — статическое давление в невозмущенном потоке, и углы атаки а = -5° ... 11°.
Полученные теневые фотографии позволили выявить структуру обтекания струй. Установлено, что в диапазоне сверхзвуковых скоростей, незначительно превышающих скорость звука, формирующийся перед струями скачок уплотнения не имеет достаточной интенсивности, чтобы вызвать отрыв потока, и структура остается безотрывной, что подтверждается данными из работы [4]. За отверстиями инжекции образуется область пониженного давления, вызванная расширением инжектируемых потоков.
По результатам весовых испытаний были получены зависимости суммарных аэродинамических коэффициентов осевой Сх^ = = Х%/(д^Бт) и нормальной Су^ = У^1/(д(Х1Бт) сил, а также момента тангажа = Мг^/(джБтЬк) конического тела при функционировании струйного органа управления для различных углов атаки а и относительных интенсивностей 7р0. Здесь Х^ и — продольная и нормальная силы; — момент тангажа; дж — скоростной напор; Бгп = — характерная площадь. На основе суммарных значений управляющих сил и момента были вычислены приращения аэродинамических коэффициентов ДСу = Су£ — Сук — СРу, ДС" = Сже — Схк — Срх, Дшг = шгЕ — шгк — шгр от перераспределения давления по боковой поверхности конического тела, где Схк, Сук — значения аэродинамических коэффициентов конуса при отсутствии выдува; СРу, СРх — коэффициенты составляющих тяги струй
Рис. 2. Зависимости суммарных аэродинамических коэффициентов нормальной (а) и продольной (б) сил при изменении относительного массового расхода выдуваемого газа
с учетом влияния кососрезанности сопел; шгР — соответствующий момент тангажа.
Проведенный на основе методов, изложенных в работах [5-7], анализ погрешностей вычисляемых величин показал, что максимальные значения относительной погрешности коэффициентов продольной силы и момента тангажа не превосходят 12%, а максимальная относительная погрешность коэффициента приращения нормальной силы составила « 7%.
На рис. 2 приведены графики суммарных аэродинамических коэффициентов нормальной (а) и продольной (б) сил при изменении относительного массового расхода выдуваемого газа Шj = Pj V Sj/Р»УсБт (Р], Р», У, Ус» — соответственно плотность и скорость инжектируемого и основного потоков; Б] — площадь среза отверстий выдува) и угле а = 0°. Видно, что общий характер зависимостей Суе(Ш]) и СЖ£ (Ш]) для вариантов а-в расположения сопел инжекции одинаков. Возрастание коэффициента Су^ по мере увеличения массового расхода газа через сопла связано в основном с увеличением коэффициента тяги струй СРу в вертикальном направлении.
Увеличение СЖ£ для всех вариантов выдува объясняется ростом глубины проникновения газообразных препятствий в набегающий поток и возрастанием интенсивности скачка уплотнения перед ними. Наибольшее значение суммарного коэффициента продольной силы обеспечивает выдув одиночной струи, что связано с определенным распределением инжектируемого потока в донной области. Аналогичный характер зависимостей наблюдается во всем исследованном диапазоне углов атаки.
Графики приращений аэродинамических коэффициентов нормальной Д Су (Зр0) и продольной Д С" (^о) сил от перераспределения давления по поверхности ЛА показаны на рис.3. Значение коэффициента ДСу для компоновочных схем с дипольным выдувом газа (варианты б-г, см. рис. 1, б) и моноструйной инжекцией (вариант а) во всем
Рис.3. Зависимости приращений аэродинамических коэффициентов нормальной (а) и осевой (б) сил от относительной интенсивности выдува
диапазоне изменения Jp0 положительно и, следовательно, направлено в сторону, противоположную тяге струй. Этот эффект связан с существенным влиянием на управляющие силы области пониженного давления за отверстиями инжекции.
В диапазоне умеренных значений относительной интенсивности выдува (Jp0 < 80) наименьшее значение коэффициента А Су обеспечивает выдув сонаправленных струй (вариант в). Это можно объяснить тем, что при инжекции струй газа под некоторым углом навстречу друг другу давление в пространстве между ними за счет столкновения потоков несколько повышается по сравнению с остальными схемами инжекции и возрастает составляющая нормальной силы, направленная в сторону действия тяги струй. Выдув газа по схеме г приводит к формированию максимального положительного значения А Су.
Графики зависимостей приращения коэффициента осевой силы А С'х показывают, что для всех исследуемых схем инжекции в диапазоне изменения параметра Jp0 = 0... 140 значение АС^, больше нуля и возрастает при увеличении Jp0.
В качестве показателя эффективности был использован коэффици-
т^ АСу + CPy - г -ент усиления Ky = ——-, представляющий собой отношение
CP max
суммы коэффициентов нормальной силы от перераспределения давления в возмущенных областях и тяги струй в вертикальном направлении к коэффициенту тяги при вертикальном истечении струй в вакуум
Pmax
из тех же сопел (CP max = —, Pmax — максимальная тяга струй).
Проведенные расчеты показали, что на таких режимах обтекания значение Ky рассмотренных вариантов инжекции меньше единицы, что свидетельствует об отрицательном влиянии заинжекционной области.
Математическая модель. Как было отмечено ранее, в диапазоне сверхзвуковых скоростей, незначительно превышающих скорость звука (M^ < 1,3), структура течения при инжекции поперечных струй газа остается безотрывной [4]. В этом случае вязкостные эффекты не
оказывают существенного влияния на параметры потока у обтекаемой поверхности и их расчет можно проводить с помощью численного метода, в основе которого лежат уравнения Эйлера.
Рассмотрим нестационарное трехмерное сверхзвуковое течение идеального газа (при отсутствии вязкости, теплопроводности и неравновесных физико-химических превращений). В декартовой системе координат система уравнений сохранения аэрогазодинамики для такой модели газа записывается в виде
dE dFx dFy dFz
dt + dx
+
дУ
+
dz
= 0,
где
E=
p PVx ' Pvy PVz
PVx PVX + P PVxVy PVxVz
Pvy • F = 1 ± x pVxVy • F = • 1 y pv2 + p • Fz= PVy Vz
pVz pVxVz pVy Vz pv2z + p
_ pio - p _ pvxio _ PVy io _ PVz io
t — время; vx
vz — составляющие вектора скорости; p,p — плот-
ность и давление газа; — удельная полная энтальпия газа. Данная система уравнений дополнена уравнением состояния газа.
Для решения системы уравнений сохранения аэрогазодинамики в невязкой постановке был использован численный метод установления второго порядка аппроксимации по всем направлениям интегрирования. Второй порядок аппроксимации по пространственным координатам достигается кусочно-линейным распределением параметров потока в пределах каждой ячейки в сочетании с принципом минимального значения производной, а повышение порядка аппроксимации по времени достигается в результате использования процедуры "предиктор-корректор". Данный метод является наиболее распространенным и применительно к настоящей задаче позволяет рассчитывать параметры потока как в сверхзвуковых, так и дозвуковых областях за криволинейным скачком уплотнения перед струями. В основе численного метода лежит принцип определения параметров потока на боковых гранях ячеек с помощью решения задачи о распаде разрыва, которая наиболее подробно описана в задачах одномерных и двумерных течений [8-11].
Этот метод был применен при расчете параметров потока и получении невязкой структуры трансзвуковго обтекания ЛА с вертикально истекающими струями.
Все численные расчеты проводились в расчетной области, ограниченной эллипсоидом вращения Е12ОЕ11, цилиндрической поверхностью Е11Е21Е22Е12, плоскостью Е21Е22, а также поверхностью конуса (рис. 4, а). В целях удобства дальнейшего разбиения расчетной области
y
Рис. 4. Расчетные область (а) и сетка (б) в меридиональной плоскости конического тела
на ячейки и создания численного алгоритма вводились две системы координат: декартова (охуг) и сеточная (1]к).
Рассмотрим более подробно построение расчетной области в меридиональной плоскости конического тела. Ее внешние границы определяются координатами опорных точек (в декартовой системе координат): 0(х0, 0), £п(хое,Ге), Е12(Хое,Ге), Я21(хе,Ге), £22(хе,Ге) и Р(хр, 0). Размеры первого блока расчетной области выбираются с учетом выстраивания в процессе расчета перед телом скачка уплотнения. Геометрические размеры боковой части расчетной области определяются таким образом, чтобы распространяющиеся возмущения в виде головного скачка уплотнения не достигали ее внешней границы. Для удобства построения расчетной области и предоставления информации все геометрические размеры относятся к радиусу миделевого сечения конуса.
Проведенные исследования показали, что при расчете взаимодействия основного и инжектируемого потоков целесообразно использовать сетку с так называемыми веерными вставками (рис.4,б). В этом случае вся расчетная область разбивается на семь блоков. Блоки 111а и Уа представляют собой веерную сетку, строящуюся от боковой кромки сопла и стыкующуюся на некотором удалении от поверхности конического блока с расчетной сеткой, содержащей прямоугольные ячейки (блоки III и У). Радиус окружности г*, описывающей внешние границы веерных вставок, определяется по формуле
г* и 1,2(1,36 М,/и - 2,0)^-;
здесь М, — число Маха на срезе сопла; и = р, / ри — среднее значение коэффициента нерасчетности струй (р, — статическое давление на
Рис. 5. Зависимость коэффициента Сук от параметра N при = 1,2;
вк = 6°; а = 3°;
■ — сетка с веерными вставками; ♦ — сетка без веерных вставок; • — экспериментальные значения
срезе сопел; рк — статическое давление на конической поверхности). Длина нижней границы блока IV, прилегающей к поверхности конуса, соответствует диаметру сопла инжекции.
Проводились исследования по влиянию общего числа разбиений N расчетной области с веерной вставкой на точность получаемых значений аэродинамических коэффициентов при обтекании сверхзвуковым потоком эталонных тел (конусов) без выдува струй. Число ячеек, приходящихся на веерную вставку и сопло инжекции, выбиралось из условия отсутствия резкого изменения геометрических размеров расчетных ячеек. Установлено, что общий характер зависимости коэффициента нормальной силы Сук (^ остается аналогичным значению Сук, получаемому на сетке с ячейками, близкими по форме к прямоугольным (рис. 5). Дальнейшие расчеты обтекания конического блока со струями проводились на сетке с общим числом разбиений N и 375 000.
В качестве начальных данных для численного решения задачи вводились следующие безразмерные параметры набегающего потока: давление = 1/(к М^), отнесенное к двойному скоростному напору; плотность р^ = р/рж = 1; полная энтальпия 0,5 + 1/[(к - 1)М2] (к — показатель адиабаты) и составляющие = йу^ = ууж/уж, = вектора скорости набегающего потока в декартовой связанной системе координат. Они задавались на внешней границе расчетной области Е12ОЕ11 и Е11Е22. На поверхности конического тела ставились условия непротекания (обнуление нормальной к поверхности конуса составляющей скорости), а на границе Е21Е22 — условие свободного вытекания газа (экстраполяция параметров потока). В меридиональной плоскости ставились граничные условия симметричного обтекания.
Для моделирования условий выдува с боковой поверхности конуса струй газа на участках, соответствующих местоположению и геометрическим размерам сопел, задавались условия звукового выдува: давление на срезе сопла р] = р]/2д^; плотность р] = р]/рж; составляющие скорости инжектируемого газа й] = 0, йу] = йу] й] = ; полная энтальпия 10] = (г>|/2) + [к]/(к] — 1)][р]/р]], где
Vj, kj — полное значение скорости и показатель адиабаты выдуваемого газа.
В результате проведенных численных исследований обтекания конического тела (вк = 6°) с ОУ был получен большой объем данных о структурах течений, распределении давления (изобары), чисел Маха (изомахи) для различных вариантов инжекции и условий выдува струй, который позволил подтвердить сделанный ранее вывод о характере распространения возмущений по поверхности ЛА. На рис. 6 приведены изомахи в продольной плоскости симметрии ЛА и изобары относительного давления на поверхности конуса при инжекции газа из одиночного сопла = 1,2; Sj/Бт = 2,67 • 10-3; а = 0°). Полученные данные отражают физическую картину обтекания. Перед инжектируемым потоком образуется скачок уплотнения, за которым течение тормозится до дозвукового. Для области за соплом инжекции характерно образование зоны разрежения, которая вблизи меридиональной плоскости формируется инжектируемой струей, а на остальном участке — основным потоком, обтекающим газообразное препятствие. На рис. 6 видно, что области возмущенного течения перед соплом выду-ва (скачок уплотнения) и за ним (область разрежения) соразмерны.
Рис. 6. Изомахи (а) и изобары (б) трансзвукового обтекания одиночной струи
при М^ = 1,2; вк = 6°; а = 0°
Рис. 7. Изобары невязкого обтекания сонаправленной схемы выдува при
MQ
1,2; ßk = 6°; а
Jpo = 24,64; б
Jp
0°:
= 91,85
Однако вследствие более существенного уменьшения давления за отверстием инжекции нормальная сила, обусловленная выдувом, будет направлена вверх, что снижает коэффициент усиления Ку.
На рис. 7 в качестве примера показаны изобары на поверхности конуса при выдуве дипольных сонаправленных струй (М^ = 1,2; Б]/Бт = 2,67 • 10-3; а = 0°). Видно, что общий характер распространения возмущений остается аналогичным обтеканию одиночной струи. Распространение областей повышенного и пониженного давлений по боковой поверхности для дипольных схем инжекции будет зависеть также и от угла взаимного наклона струй. С его увеличением степень распространения возмущений на противоположную отверстиям выдува поверхность ЛА возрастает.
Полученные по результатам численных расчетов зависимости приращения коэффициента момента тангажа от угла атаки показаны на рис. 8. Хорошее согласование расчетных и экспериментальных значений достигалось в диапазоне малых и средних интенсивностей выдува (до значений 7р0 < 100). Рассогласование значений Ашг при больших значениях 7р0 можно объяснить тем, что выдуваемая струя присоединяется за кормовым срезом ЛА и в заинжекционную область поступает газ с более высоким давлением из донного следа.
Было исследовано влияние местоположения сопел инжекции при изменении ] = 0,9... 1,0 на значение коэффициента усиления. Установлено, что по мере приближения сопел выдува к кормовому срезу ЛА эффективность применения струй на таких режимах обтекания увеличивается.
Выводы. Проведенный комплекс экспериментальных исследований позволил выявить особенности формирования аэродинамических характеристик при инжекции одиночной и систем вертикальных струй газа в трансзвуковой поток из сопел, удаленных от кормового среза
а
Дот
Рис.8. Зависимости Дтг(а) различных видов инжекции при = 1,2;
вк = 6°:
а — выдув одиночной струи; б — выдув параллельных струй;--расчет;
----эксперимент
ЛА, а также установить закономерности их изменения в зависимости от интенсивности выдува струй и угла атаки ЛА. Установлено, что на режимах обтекания, близких к звуковым (в условиях образования безотрывной структуры течения), значение А Су всех исследованных схем выдува положительно и направлено в сторону, противоположную тяге струй. При этом в диапазоне интенсивностей выдува 7р0 ^ 90 струи, инжектируемые под углом е = -60°, оказываются наиболее эффективными.
С помощью предложенного метод расчета был проведен численный эксперимент и выявлены структуры обтекания инжектируемых струй, которые позволили обосновать и подтвердить экспериментальные результаты. Работоспособность и адекватность численного метода подтверждается удовлетворительным согласованием расчетных и экспериментальных значений аэродинамических коэффициентов, а также адекватной структурой обтекания конического тела и инжектируемых потоков.
СПИСОК ЛИТЕРАТУРЫ
1. Калугин В. Т., Чернуха П. А. Аэродинамическая эффективность управления процессами обтекания конических тел при инжекции дипольных струй // Оборонная техника. - 2002. - № 10-2. - С. 102-105.
2. Калугин В. Т., Луценко А. Ю, Столярова Е. Г., Чер-н у х а П. А. Активное управление инжекцией блочных струй // Оборонная техника. - 2003. - № 1-2. - С. 71-76.
3. Ч е р н у х а П. А. Сравнительный анализ эффективности использования различных схем выдува струй для управления обтеканием летательных аппаратов // Научный вестник МГТУ ГА. № 59. Сер. "Аэромеханика и прочность". - М.: Изд-во МГТУ ГА, 2003. - С. 110-113.
4. Ч ж е н П. Отрывные течения / Пер. с англ. - М.: Мир, 1972. - Т. 1. - 300 с.
5. Афанасьев А. А. Определение погрешностей измерения аэродинамических коэффициентов при исследованиях обтекания тел с вдувом кольцевой струи навстречу потоку // Тр. МВТУ им. Н.Э. Баумана. - 1987. - № 492. -С. 11-19.
6. Карданов Ю. Х., Захарченко В. Ф., Семенихин А. Н. Трех-компонентные вынесенные донные малогабаритные аэродинамические весы для струйных испытаний // Тр. МВТУ им. Н.Э. Баумана. - 1978. - № 274. Вып. 1. -С. 30-41.
7. Витушкин В. В., Захарченко В. Ф. Погрешности измерений при исследованиях обтекания тел с поверхностным массообменном // Тр. МВТУ. -1982. - № 386. - С. 50-59.
8. Годунов С. К., Забродин А. В., Прокопов Г. П. Разностная схема для двумерных нестационарных задач газовой динамики и расчет обтекания с отошедшей ударной волной // Журнал вычислительной математики и математической физики. - М. - 1961. - Т. 1. - № 6. - С. 1020-1050.
9. К о т о в А. И., Матвеев С. К., С и м о н е н к о М. М. Численный расчет методом Годунова нестационарных течений с контактным разрывом // Деп. ВИНИТИ. - 1986. - № 6938-В86. - 11 с.
10. К о л г а н В. П. Применение принципа минимальных значений производных и построение конечно-разностных схем для расчета разрывных решений уравнений газовой динамики // Ученые записки ЦАГИ. - 1972. - Т. 3. - № 6. -С. 68-77.
11. Численное решение стационарных и нестационарных задач газовой динамики при помощи монотонных явных разностных схем первого порядка точности: ТО / ЦИАМ; Исп. В.Т. Гринь, М.Я. Иванов, А.Н. Крайко, Р.К. Тагиров. - Инв. № 7614. - 1975.
Статья поступила в редакцию 13.06.2006
Владимир Тимофеевич Калугин родился в 1949 г., окончил МВТУ им. Н.Э. Баумана в 1972 г. Д-р техн. наук, профессор кафедры "Баллистика и аэродинамика" МГТУ им. Н.Э. Баумана. Автор более 250 научных работ в области аэрогазодинамики струйных и отрывных течений, проектирования органов управления полетом.
V.T. Kalugin (b. 1949) graduated from the Bauman Moscow Higher Technical School in 1972. D. Sc. (Eng.), professor of "Ballistics and Aerodynamics" department of the Bauman Moscow State Technical University. Author of more than 250 publications in the field of aero-gasdynamics of jet and detached flows, design of flight controls.
Полина Алексеевна Чернуха родилась в 1978 г., окончила МГТУ им. Н.Э. Баумана в 2001 г. Канд. техн. наук, доцент кафедры "Баллистика и аэродинамика" МГТУ им. Н.Э. Баумана. Автор более 20 научных работ. Специализируется в области отрывных, струйных течений и управления процессами обтекания летательных аппаратов.
P.A. Chernukha (b. 1978) graduated from the Bauman Moscow State Technical University in 2001. Ph. D. (Eng.), assoc. professor of "Ballistics and Aerodynamics" department of the Bauman Moscow State Technical University. Author of more than 20 publications in the field of jet and detached flows and control of flow around flying vehicles.
