Модели и методы аэродинамики 2001-2002. I и II международных школсеминаров мцнмо москва 2002 удк 533. 6 532. 5 629. 7 Модели и методы аэродинамики Материалы I и II международных школсеминаров. М мцнмо, 2002 124 с. Isbn 5940570372 Сборник включает тезисы
 Скачать 1.19 Mb.
|
А.Б. Ватажин ЦИАМ им. ПИ. Баранова, Москва Рассмотрены электрофизические эффекты, возникающие при работе авиационных реактивных двигателей и обусловленные наличием в тракте двигателя заряженных частиц. Эти эффекты не влияют на само газодинамическое течение, но вызывают электризацию летательного аппарата (ЛА), влияют на конденсацию в двигательных струях и дают возможность проводить диагностику работы двигателя электрическими методами. Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 17 Основные результаты получены в следующих направлениях − Изучено нарушение электрической квазинейтральности потока в тракте двигателя, обусловленное различием скоростей диффузии электронов и ионов (возникших в камере сгорания) к поверхностям внутренних элементов двигателя. Разработана теория электрических диффузионных пограничных слоев, сформулированы и решены задачи об их развитии внутри ламинарного и турбулентного пограничного слоя на плоской пластине ив окрестности критической точки обтекаемого тела например, лопатки. Результаты расчетов электрического тока выноса из двигателя (вызывающего электризацию ЛА) согласуются сданными аэродромных испытаний. Построена теория важного эффекта – исчезновения тока выноса из двигателя на режиме форсажа. − На основе созданной физико-математической модели конденсационных турбулентных струй выполнено численное моделирование лабораторных паровоздушных струй при наличии гомогенной и гетерогенной конденсации и конденсации на ионах, попадающих в струю из сопла или из окружающего пространства. Проведено сопоставление расчетов с экспериментальными данными. Определены поля газодинамических и кинетических параметров и степень конверсии пара в конденсированную дисперсную фазу. Рассмотрены конденсационные эффекты в двигательных струях современных самолетов. − Теоретически обоснован и разработан метод бесконтактной электростатической диагностики состояния авиационных двигателей, основанный на регистрации зондами-антеннами, расположенными вне двигателя и вне его струи, переменных электрических полей, генерируемых находящимися в струях заряженными частицами. Сигналы с зондов обрабатываются в виде спектров, в результате чего составляется электрический портрет двигателя. Анализ полученной информации позволяет выяснять особенности и специфику работы двигателя. Проведены аэродромные испытания на современных самолетах. − Разработан новый метод обнаружения начала разрушения металлических тел (элементов двигателя. В его основе лежит эффект появления большого числа положительно заряженных микрочастиц при разрушении образцов и регистрация электрических сигналов от этих частиц. Метод апробирован в лабораторных и стендовых условиях. Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 18 Экспериментальное исследование эффективности различных способов инжекции барботированного газом керосина в сверхзвуковую камеру сгорания О.В. Волощенко, С.А. Зосимов, В.В. Иванов, А. Николаев, В.А. Сабельников ЦАГИ им. проф. НЕ. Жуковского, Жуковский Работа основана на результатах экспериментов по горению барботированного газом керосина в прямоточных камерах сгорания со сверхзвуковой скоростью на входе. Эксперименты проводились на присоединенном воздухопроводе при параметрах торможения P t = 14 ÷18 атм и T t ≅ 1700 К, число Маха на входе 2.5. Условия соответствуют М полета 6. Целью работы было исследование возможных способов организации высокоэффективного горения керосина в сверхзвуковых камерах сгорания. Для интенсификации процессов смешения и горения в сверхзвуковой камере сгорания использовалось барботирование, те. насыщение газом в количестве 5 ÷10% от массы керосина. Кроме очевидной интенсификации смешения это ведет к тому, что появляется возможность использовать имеющийся опыт по интенсификации и управлению смешением как при подаче струй водорода. В работе исследовалось влияние на рабочий процесс в камере сгорания различных факторов − влияние формы сопел для подачи топлива (круглой и овальной, − влияние типа газа для барботирования (водород, воздух) и массовой доли газа, − влияние длины камеры сгорания на эффективность горения и границу бедного срыва, − исследовались различные инжекторы трубчатые, стреловидные, клиновидные, ряд вариантов подачи со стенки – “аэро- рамп. Были получены следующие основные результаты − определены границы устойчивой работы камеры сгорания по коэффициенту избытка топлива. − определены зависимости полноты сгорания от коэффициента избытка топлива для различных способов инжекции. − предложена физическая модель рабочего процесса в камере сгорания в условиях проведения экспериментов. Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 19 Экспериментальное исследование газодинамических методов организации горения в сверхзвуковом потоке О.В. Волощенко, МА. Иванькин, В.В. Иванов, В.А. Сабельников ЦАГИ им. проф. НЕ. Жуковского, Жуковский В работе представлены результаты экспериментального исследования газодинамических методов стабилизации горения углеводородных топлив в сверхзвуковом потоке. Получено самовоспламенение и устойчивое горение водорода в головной свободновисящей рециркуляционной зоне, образованной при разрушении звуковой струи водорода в результате интерференции с сильным скачком уплотнения, генерируемым незапущенным осесимметричным диффузором (профилированным телом с протоком. Определены границы устойчивого горения в зоне перед телом с протоком. Получено самовоспламенение и горение водорода в передней отрывной зоне, создаваемой при обтекании осесимметричного профилированного тела с протоком с центральной иглой, при подаче водорода через иглу навстречу потоку. Игла установлена в носовой части тела с протоком и выступает перед плоскостью входа на расстояние, где D – диаметр входа тела с протоком. Получены характерные режимы течения для игл различной длины. Получено самовоспламенение и устойчивое горение водорода в свободновисящей рециркуляционной зоне, расположенной в следе за осесимметричным профилированным телом с протоком. Струя водорода выдувается через пилон соосно струе воздуха, проходящей через тело с протоком. Пилон, расположен в хвостовой части тела с протоком. Свободновисящая рециркуляционная зона возникает при разрушении струи водорода в результате интерференции со скачком уплотнения впервой бочке нерасчетной струи воздуха, истекающей из тела с протоком. Эксперименты проведены на аэродинамическом стенде ЦАГИ Т-131В в сверхзвуковом потоке на выходе из плоских расширяющихся каналов при числах Маха М ≈ 2.5-2.6, температура и давление торможения в воздухоподогревателе T t ≈ 1200÷1500 К, P t ≈ 2.7÷3 МПа. Работа выполнена при финансовой поддержке РФФИ (проект № 00-01-00158). Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 20 Влияние магнитного поляна тепломассоперенос и испарение капель магнитных жидкостей В.В. Гогосов Институт механики МГУ им. МВ. Ломоносова, Москва Е.В. Зубенко, Х.Д. Искендеров, МА. Кобозев, А.Я. Симоновский Ставропольская государственная сельскохозяйственная академия, Ставрополь Исследуются процессы испарения капель магнитных жидкостей (МЖ) на горячей поверхности. Проводились следующие эксперименты. На нагретую поверхность падает капля магнитной жидкости. При разных температурах поверхности нагревателя измеряется краевой угол α между каплей и поверхностью нагревателя. Измеряется время τ испарения капли в зависимости от ее объема и температуры поверхности нагревателя, а также от концентрации магнитной фазы в жидкости. Такие измерения проводились как без магнитного поля, таки при наличии горизонтального или вертикального относительно поверхности нагревателя однородного магнитного поля разной интенсивности. Показано, что испарение капель концентрированной МЖ вин- тервале температур от С до С отличается от испарения обычных жидкостей. С ростом температуры нагревателя во всем изученном интервале температур происходит монотонное уменьшение. При температурах поверхности нагревателя С внутри капель концентрированной МЖ образуется сгусток шарообразной формы темно-коричневого цвета, окруженный слоем прозрачной жидкости. Стечением времени прозрачная оболочка капли испаряется, а на месте испарившейся капли остается круглой формы темный осадок. Обнаружено, что магнитное поле влияет не только на поведение α и τ, но и приводит к возникновению макроструктур на поверхности капель. Например, в слабо концентрированных жидкостях на поверхности капель в начальный период испарения образуются трещины древовидной формы, расположение и конфигурация которых зависят от направления магнитного поля. Это свидетельствует о превращении при высоких температурах капель МЖ жидкости в гели. Представляет интерес поведение капель концентрированной МЖ в интервале температур поверхности нагревателя от С до С. Испарение происходит не только со свободной поверхности капли, но и внутрь объема капли от ее поверхности, контактирующей с Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 21 поверхностью нагревателя. При этом в объеме капли образуются пузырьки пара. Это приводит к увеличению диаметра капли в 1.5 ÷2 раза. В последующие моменты времени в результате испарения размер капли начинает уменьшаться. По окончании испарения на поверхности нагревателя вместо капли остается осадок, имеющий вид смятой сферической оболочки. Препарирование остатка капли тонким лезвием показывает, что в конце испарения она представляет собой полую тонкостенную оболочку, испещренную дырочками – местами выхода формирующихся пузырьков пара через поверхность оболочки. Предлагается объяснение наблюдающихся в экспериментах явлений. Работа выполнена при финансовой поддержке РФФИ (проект № 99-01-01057). Оптимальные крылья в гиперзвуковом потоке газа В.Н. Голубкин, В.В. Сысоев ЦАГИ им. проф. НЕ. Жуковского, Жуковский Разработан эффективный численно-аналитический подход коп- ределению максимального аэродинамического качества и соответствующих оптимальных форм тонких крыльев малого удлинения в гиперзвуковом потоке газа. Рациональная формулировка вариационной задачи получена с использованием общего аналитического решения задачи пространственного обтекания крыла методом тонкого ударного слоя и интегральных законов сохранения. Для решения задачи оптимизации при заданном угле атаки, площади в плане и других ограничениях использованы классические методы вариационного исчисления и прямые численные методы. Рассмотрены характерные случаи оптимизации крыльев сострой и слабо затупленной передними кромками. Обнаружено бифуркаци- онное поведение оптимального решения при изменении размаха крыла, сопровождающееся качественными изменениями оптимальных форм. Выявлены особенности пространственной геометрии крыла, приводящие к увеличению качества. Указаны классы пространственных оптимальных форм, имеющих существенно большее гиперзвуковое качество по сравнению с плоскими крыльями в широком диапазоне углов атаки. Установлено, что наличие малого затуп- ления передней кромки существенно влияет на конфигурацию Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 22 оптимальных форм, найден параметр подобия, характеризующий относительное влияние сопротивления затупления, и получены характерные немонотонные зависимости качества от угла атаки. Работа выполнена при финансовой поддержке РФФИ (проект № 99-01-01128). Экстремальные задачи ускорения тел сжатым газом АН. Голубятников, НЕ. Леонтьев МГУ им. МВ. Ломоносова, Москва При расчетах внутренней баллистики метательных систем на определенной стадии разлета может быть использован класс точных решений уравнений газовой динамики с однородной деформацией, допускающий известный функциональный произвол. Такой подход позволяет ставить и решать некоторые задачи оптимизации работы этих систем при подходящих ограничениях на начальные и краевые условия [1]. Использование предельных соотношений, полученных из рассмотрения движения газа между двумя поршнями с однородной деформацией, позволяет на начальном этапе оптимизации работы мно- гопоршневых баллистических установок [2] аппроксимировать многоступенчатую газодинамическую систему конечномерной системой последовательно взаимодействующих твердых тел. В рамках такой дискретной модели, применимой для описания достаточно разреженных систем, задача оптимизации скорости метания за счет выбора масс промежуточных слоев газа и поршней сведена к задаче Герца о выборе наилучшего распределения масс упруго сталкивающихся материальных точек, решение которой дает убывающую геометрическую прогрессию. Особенностью данной задачи является то, что, несмотря на относительно невысокий КПД, здесь теоретически могут быть достигнуты значительные скорости метания. Важным элементом является задача оптимизации работы первой ступени установки с учетом специально деформирующейся боковой границы. На практике подвижная граница может деформироваться, например, продуктами детонации заряда взрывчатого вещества [3]. Показано, что при двусторонних ограничениях на начальную плотность, как ив случае метательных установок с продольным движением газа [1], при заданной форме боковой границы распределение плотности, обеспечивающее наибольший КПД, является кусочно Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 23 постоянным. С другой стороны, в случае постоянной начальной плотности газа наилучшим образом выбранная подвижная граница необходимо должна состоять из цилиндрических и плоских участков, причем в зависимости от заданных параметров системы не исключается наличие движения задней стенки. Работа выполнена при финансовой поддержке РФФИ (проекты № 00-15-96154, 00-01-00135). Литература 1. Голубятников АН К оптимальной постановке газодинамической задачи Лагранжа // Вестн. МГУ. Сер. 1, Матем. Механ. 1995, № 6, с. 59-61. 2. Жаровцев В.В., Комаровский Л.В., Погорелов Е.И. Математическое моделирование и оптимальное проектирование баллистических установок. Томск Изд-во ТГУ, 1989. 256 с. 3. Пилюгин Н.Н., Леонтьев НЕ Возможности повышения скорости метания тел в баллистических установках. Институт механики МГУ, Препринт № 52-99, 1999. 58 с. Использование излучения радикалов ОН для определения полноты сгорания углеводородных топлив МА. Гольдфельд, С.Г. Миронов, А.А. Мишунин, А.В. Потапкин ИТПМ СО РАН, Новосибирск Полнота сгорания топлива в энергетических установках и двигателях является важным показателем степени совершенства их конструкции и эффективности использования топлива. В настоящее время полнота сгорания топлива определяется, в основном, калориметрическим методом и методом химического анализа продуктов сгорания. Возникают большие трудности при их применении для исследования горения в каналах ГПВРД при сверхзвуковых скоростях потока и особенно при испытаниях в аэродинамических трубах кратковременного действия. Решением этой проблемы является создание оптических методов контроля полноты сгорания. В работе представлена методика определения полноты сгорания керосина по интенсивности излучения радикалов ОН в ультрафиолетовом диапазоне 270 ÷360 нм. Методика аналогична ранее разработанной методике для водородных пламен, в основе которой лежит предположение о пропорциональности интенсивности свечения радикалов ОН и интенсивности реакций окисления сих участием. Величина полноты сгорания топлива определялась по отношению Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 24 величины интенсивности свечения радикалов ОН в исследуемом пламени и калибровочном факеле при одинаковых расходах водорода. При этом условия горения калибровочного факела должны обеспечивать полное сгорание топлива. Спектральный диапазон свечения радикалов ОН выделялся комбинацией стеклянного и жидкостного фильтра с полосой пропускания нм. Регистрация отфильтрованного излучения производилась фотоумножителем через систему щелевой развертки изображения с целью выделения области излучения реакционной зоны пламени. Апробация оптической системы проводилась на свободном водородном пламени. Водород подавался из баллона высокого давления через систему трубопроводов и струйную форсунку. Данные настоящего исследования, проведенного с водородным пламенем, хорошо согласуются сданными, полученными в работах С.С. Воронцова (1976) и В.К. Баева (1984). Разработанная оптическая система была применена для регистрации интенсивности излучения радикалов ОН при горении паров керосина в воздухе. Пары керосина создавались путем нагревания авиационного керосина в теплоизолированном сосуде высокого давления. Выпуск паров керосина производился с помощью быстродействующего пироклапана через многощелевую форсунку, а поджиг – дежурным факелом. В течение всего процесса измерялось давление в сосуде. Расход паров керосина вычислялся по темпу падения давления в сосуде и расходу жидкого керосина в эксперименте. Таким образом, были получены зависимости интенсивности свечения радикалов ОН от расхода керосинового пара при полноте сгорания, близкой к 1. Проведенные исследования позволили сделать следующие выводы. В области малых расходов (< 2 гс) в водородных пламенах зависимость интенсивности свечения радикалов ОН от расхода водорода имеет нелинейный характер. В углеводородных пламенах зависимость интенсивности свечения ОН радикалов от расхода топлива линейна в диапазоне от 0 до 60 гс. Исследование воспламенения и горения керосина в модели полного двигателя. МА. Гольдфельд, Р.В. Нестуля, А.В. Старов ИТПМ СО РАН, Новосибирск До скоростей полета, соответствующих числам Маха M = 7 ÷8, предполагается использование жидких углеводородных топлив. При Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 25 этом наиболее сложным вопросом является организация воспламенения и горения жидких углеводородов при сверхзвуковой скорости потока на входе в канал камеры сгорания. В работе представлены результаты исследования полной модели ГПВРД с горением керосина. Эксперименты были проведены в импульсной аэродинамической трубе ИТ-302М ИТПМ СО РАН при параметрах торможения потока, близких к полетным при числах Маха 5 и 6. Основными целями испытаний модели являлись a) исследование воспламенения и стабилизации горения жидкого керосина в камере сгорания ГПВРД в импульсной установке b) получение положительной эффективной тяги c) сравнение экспериментальных результатов с расчетными данными. Модель полного двигателя состоит из трех модулей воздухозаборника, камеры сгорания и сопла. Двумерный трехскачковый воздухозаборник с полным углом поворота потока 23.5 ° обеспечивает восьмикратное геометрическое сжатие. Конструкция камеры сгорания предусматривает организацию за воздухозаборником внезапного расширения канала модели (ступенька) с возможностью вдува топлива со ступеньки спутно потоку. Камера сгорания имеет короткий участок постоянной площади и расширяющуюся часть при общей длине 20 калибров. Дополнительный вдув топлива осуществляется с пилонов или со стенки камеры сгорания под различными углами к потоку. Плоское сопло было исследовано в двух вариантах с относительным расширением 1 и 1.62. Модульный принцип конструкции модели позволил получить характеристики отдельных элементов и двигателя в целом. Измеренные параметры модели включали 1) распределения статического давления вдоль канала модели 2) поля полных давлений на выходе камеры сгорания и сопла 3) распределения тепловых потоков вдоль канала камеры сгорания 4) силы, действующие на модель (с помощью трехкомпонентных весов 5) расход топлива 6) полноту сгорания с помощью регистрации излучения в ультрафиолетовом диапазоне. Опыты были проведены при избытках керосина 0.6 ÷1.2. Было получено, что для воспламенения керосина требовался пилотный факел водорода при расходе 3 ÷6% от массового расхода керосина. Горение керосина приводило к значительному повышению статического давления и тепловых потоков (4 ÷5 раз. Измерения на выходе камеры сгорания показали, что скорость потока в этом сечении при горении становится близкой к звуковой, но остается сверхзвуковой, Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 26 что позволяет сделать вывод о реализации горения керосина в сверхзвуковом потоке воздуха. Проведенные экспериментальные исследования позволили сделать следующие выводы a) необходимо использовать пилотный факел водорода для обеспечения воспламенения жидкого керосина b) получено устойчивое горение керосина при сверхзвуковой скорости потока и показана возможность его изучения в аэродинамических установках кратковременного действия (70 ÷120 мс c) получена положительная эффективная тяга при высоком уровне внутренней тяги d) сравнение экспериментальных данных с расчетными результатами показало их удовлетворительное соответствие. Моделирование аэродинамики межступенного отсека составной ракеты В.А. Горяйнов, СВ. Коннов МАИ им. С. Орджоникидзе, Москва Аэродинамика межступенного отсека составной ракеты является базовым аспектом задачи разделения ступеней. Рассматривается вариант компоновочной схемы, когда пространство межступенного отсека ограничено передним сферическим днищем центрального бака й ступени и задним плоским днищем с выступающем соплом двигательной установки (ДУ) й ступени в продольном направлении и двумя боковыми цилиндрическими баками в поперечном направлении. Ставится задача математического моделирования течения в меж- ступенном отсеке с характерной пространственной конфигурацией на основе модели Эйлера в квазитрехмерной постановке. Для решения задачи используется метод Годунова второго порядка точности [1, 2]. Алгоритм метода Годунова базируется на пересчете газодинамических параметров в ячейках сетки при переходе от момента времени к моменту t 0 + τ с использованием интегральных законов сохранения) Массообмен внешнего потока стечением в отсеке в окрестности плоскости стыковки центральных и боковых баков осуществляется Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 27 за счет мерных течений, которые моделируются специальным заданием источниковых членов в системе уравнений (1) на границах расчетной области. Проведенное моделирование аэродинамики межступенного отсека с характерной пространственной конфигурацией дает представление о картине течения и ударно-волновых структурах на различных режимах полета при М = 2 ÷7, включая этап разделения ступеней. Полученные оценки локализации зон отрыва пограничного слоя позволяют прогнозировать экстремально теплонапряженные участки обшивки типового межступенного отсека. На одном из режимов полета до включения сопла й ступени обнаружено аномальное появление автоколебаний давления, сопровождающееся структурной перестройкой течения. Такие аномальные автоколебания получены и для осесимметричной каверны, представляющей пространство между двумя близко расположенными соосными цилиндрами, обтекаемыми в продольном направлении стационарным потоком без внешних возмущений, так что это явление можно считать типичным для течений между соосными цилиндрами как с плоскими, таки неплоскими торцами. На этапе разделения ступеней при включенной ДУ обнаруженные автоколебания давления могут достигать значительной амплитуды, соизмеримой со средним давлением в межступенном отсеке. Предложенный математический аппарат позволяет проигрывать нештатные ситуации в интересах обеспечения безопасности полета многоступенчатых ракет. Литература 1. Van Leer B. Towards the Ultimate Conservative Difference Scheme, A Second Order Sequel to Godunov’s Methods // J.Comput.Phys.1979. Vol. 32.1. pp. 101-136. 2. Горяйнов В.А., Молчанов А.Ю. Метод представления вещественных параметров дискретными аналогами в задачах математической физики. // X Юбилейная Международная Конференция “ВМСППС”, Переславль- Залесский, с. 48-50, 1999. Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 28 Влияние кинетической неравновесности на генерацию турбулентности ЮН. Григорьев ИВТ СО РАН, Новосибирск ИВ. Ершов НГАВТ, Новосибирск Исследования ламинарно-турбулентного перехода (ЛТП) и генерации турбулентности в кинетически неравновесных молекулярных газах получили развитие лишь недавно. Известно, что при относительно невысоких уровнях возбуждения, когда колебательные степени свободы молекул газа остаются замороженными, неравновесные течения газов описываются системой уравнений Навье–Стокса сжимаемого газа, в которых неравновесность учитывается через коэффициент объемной вязкости. В работе [1] была проведена серия экспериментов, зафиксировавших в течении Пуазейля существенное до десяти процентов) возрастание критического числа Рейнольдса ЛТП с увеличением объемной вязкости газа. К сожалению, по ряду причин эти результаты представляются спорными. В работах [2, 3], отчасти инспирированных результатами [1], рассматривалась линейная устойчивость пограничного слоя на пластине. Было показано, что в рамках модели объемной вязкости влияние возбуждения внутренних степеней на ЛТП мало. Вместе стем в расчетах [2] для сильной неравновесности, описываемой релаксационным уравнением Ландау–Теллера для внутренней энергии, было получено значительное возрастание коэффициентов усиления для первой и второй мод возмущений. Нов силу известных ограничений линейной теории эти результаты имеют скорее качественный характер и не экстраполируются на нелинейную стадию развития возмущений. В этой связи представляет интерес непосредственное исследование нелинейных эффектов и их вклада в среднее течение. Известно, что нелинейная стадия имеет универсальный характер и реализуется через зарождение, эволюцию и распад характерных вихревых образований. В работе выполнено численное моделирование взаимодействия поперечной организованной вихревой структуры со средним сдвиговым потоком. В расчетах использовались полные уравнения Навье– Стокса сжимаемого газа для нескольких возможных значений отношения коэффициентов объемной и динамической вязкостей. Показано, что с увеличением величины объемной вязкости средние по Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 29 пространству и времени рейнольдсовы напряжения, генерируемые структурой, возрастают на 5 ÷10%, а соответствующий вклад в суммарные рейнольдсовы напряжения оценивается в 2 ÷4%. Такое изменение по порядку величины соизмеримо с эффектом используемых на практике способов снижения сопротивления, например, с помощью риблет. Литература 1. Nerushev A., Novopashin S. Rotational Relaxation on Transition to Turbu- lence.// Phys. Lett., (1997), A 232, pp. 243-245. 2. Bertolotti F.B. The Influence of Rotational and Vibrational Energy Relaxation on Boundary-Layer Stability.// J. Fluid Mech., (1998), 372, pp. 93-118. 3. Григорьев ЮН, Ершов ИВ К вопросу о влиянии вращательной релаксации на ламинарно-турбулентный переход Тез. Док. Юбилейной науч. конф, посвященной 40-летию Инта механики МГУ, 22-26 ноября 1999, Москва, МГУ, с. 65-66. Расчетные исследования параметров пограничных слоев на несущих элементах различных форм в плане при их обтекании воздушным потоком с большой сверхзвуковой скоростью А.А. Губанов, С.А. Таковицкий ЦАГИ им. проф. НЕ. Жуковского, Жуковский Рассматриваются характеристики пограничных слоев, формирующихся на нижних поверхностях несущих элементов различных форм в плане (треугольной и прямоугольной) при большой сверхзвуковой скорости обтекания. Исследования проведены при числе Маха М = 4 на базе численных расчетов с целью оценки влияния формы несущего элемента на развитие пограничного слоя и на характеристики воздухозаборника, расположенного в области заторможенного потока непосредственно у поверхности элемента. Расчеты основаны начисленном интегрировании полной системы уравнений Навье–Стокса с использованием алгебраической модели турбулентности Болдуина–Ломакса. Линии тока в пограничном слоев непосредственной близости к поверхности элемента, вследствие меньших величин продольной составляющей скорости, при наличии поперечного градиента давления имеют большую кривизну, чем на внешней границе пограничного слоя. На треугольном несущем элементе зоны повышенного Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 30 статического давления располагаются в окрестностях кромок, и поэтому линии тока в пограничном слое отклоняются в направлении плоскости симметрии. В результате это приводит к утолщению пограничного слоя на несущем элементе и, следовательно, к ухудшению характеристик расположенного под ним воздухозаборника. На прямоугольном несущем элементе, наоборот, линии тока отклоняются в стороны боковых кромок, и происходит самопроизвольное растекание пограничного слоя, благоприятное для размещения воздухозаборника. Количественные оценки параметров пограничного слоя показывают, что влияние формы элемента на параметры пограничного слоя существенно, и его необходимо учитывать при выборе конфигураций элементов аэродинамических компоновок летательных аппаратов, используемых в качестве предварительных ступеней торможения воздухозаборников. Работа выполнена при финансовой поддержке РФФИ (проект № 99-01-01128). Экспериментальное исследование течений торможения и смешения в каналах Н.В. Гурылева ЦАГИ им. проф. НЕ. Жуковского, Жуковский При проектировании многорежимных прямоточных воздушно- реактивных двигателей большое значение имеет исследование способов управления течениями торможения в протяженных каналах двигателя, когда в псевдоскачке происходит переход от сверхзвукового к дозвуковому течению. Рассмотрены особенности структуры течения и параметры потока при реализации свободного и фиксированного псевдоскачка в прямоугольных плоских и осесимметричных каналах и определен ряд факторов, влияющих на фиксацию псевдоскачка в канале. Исследовано торможение сверхзвукового потока (M = 1.8 ÷2.5) в прямоугольных и осесимметричных каналах при наличии противодавления, создаваемого механическим дросселированием. Показано, что для свободного псевдоскачка в прямоугольном канале характерно наличие несимметричных отрывных зон вблизи верхней и нижней, а также боковых стенок канала. В области головной части псев- доскачка, наблюдаются существенные пульсации картины течения при одних и тех же параметрах набегающего потока Ми Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 31 степени дросселирования канала. Наблюдаются продольные и поперечные пульсации скачков уплотнения, причем изменение претерпевает не только положение скачков уплотнения, но и их структура происходит изменение вида отрывных зон на стенках канала – хорошо заметен чередующийся переход от закрытой локальной отрывной зоны к открытой вблизи стенки. С увеличением числа Мам- плитуда и частота пульсаций скачков возрастает. Показано, что течение с псевдоскачком не может быть рассмотрено в рамках квазистационарных моделей. Показана возможность фиксации псевдоскачка на входных кромках прямоугольных каналов. Получено, что для развитой фиксации псевдоскачка на входных кромках прямоугольных и осесимметричных каналов цепочка замыкающих скачков вырождается в один скачок. Отмечено, что при этом на внутренней поверхности осесимметричного канала наблюдается трехмерный отрыв пограничного слоя, содержащий четыре вихря, расположенных симметрично по диаметру канала. Структура линий тока на стенке характерна для закрытой отрывной зоны. Дано сравнение уровней фиксации для осесимметричных и плоских каналов. Для чисел М = 2.6 ÷3.8 рассмотрена фиксация псевдоскачка на пилонах, расположенных в глубине цилиндрических каналов. Показано, что фиксация псевдоскачка на цилиндрических пилонах в глубине каналов наблюдается при устойчивом перемещении псевдо- скачка, те. в области турбулентного пограничного слоя. Получена зависимость степени фиксации псевдоскачка от высоты пилонов, частоты их размещения и наличия вдува воздуха через пилоны. Показано, что уровень фиксации на пилонах в глубине канала ниже, чем на входных кромках. Получено, что интерференция вихревого шнура, генерируемого перед плоскостью входа канала, со скачками уплотнения ускоряет перемещение начала псевдоскачка вверх по потоку попадание низконапорной осесимметричной струи (Мс = 1) в канал не вызывает ускорения перемещения начала псевдоскачка вверх по потоку. При взаимодействии струи и вихревого шнура с псевдоскачком, фиксированном на входной кромке канала, возможно образование перед входом в канал локальной зоны рециркуляционного течения, выступающей перед плоскостью входа. Показано, что в ряде случаев это может привести к помпажу. Взаимодействие вихревого шнура с псевдоскачком в глубине канала вызывает ухудшение характеристик Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 32 канала, таких, как коэффициент восстановления полного давления, противопомпажный запас. При взаимодействии вихревого шнура с псевдоскачком, сопровождающемся образованием на входе в канал локальной рециркуляционной зоны, возможно улучшение характеристик канала. В модельном течении показана возможность реализации сильного решения в режиме высокочастотных низкоамплитудных пульсаций скачка, что расширяет диапазон безпомпажного режима течения. Полученные результаты могут быть использованы для интенсификации процессов торможения, смешения и организации энерго- подвода в сверхзвуковых потоках. Работа выполнена при финансовой поддержке РФФИ (проект № 00-01-00158). Расчетные исследования аэродинамических характеристик сверхзвуковых ЛА с ВРД в рамках модели трехмерных стационарных уравнений Эйлера Д.Ю. Гусев, В.В. Коваленко ЦАГИ им. проф. НЕ. Жуковского, Жуковский Проведены расчетные исследования аэродинамических характеристик компоновок летательных аппаратов с плоскими воздухозаборниками при их различном расположении около осесимметричного корпуса. Представлены также интегральные характеристики потока во входных сечениях воздухозаборников при числах Маха Ми углах атаки α = 0÷20°. Расчет обтекания компоновок осуществлялся при помощи программы, основанной начисленном интегрировании полной системы трехмерных стационарных уравнений Эйлера. Внешние аэродинамические характеристики ЛА определялись в результате интегрирования статического давления по внешним поверхностям компоновок с учетом соответствующих поправок на влияние протока воздуха через двигательную установку и сопротивление трения. Показано, что для широкого диапазона углов атаки при М = 3 для компоновки осесимметричного корпуса с двумя плоскими воздухозаборниками среди рассмотренных схем наиболее благоприятной является схема с воздухозаборниками, находящимися на наветренной части корпуса и разнесенными на угол 45 ° относительно вертикальной плоскости симметрии. Работа выполнена при финансовой поддержке РФФИ (проект № 99-01-01128). Первая Международная школа-семинар МОДЕЛИ И МЕТОДЫ АЭРОДИНАМИКИ 33 Структурно параметрический анализ системы управления |