Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,074

АНАЛИЗ МЕТОДОВ ОПРЕДЕЛЕНИЯ АЭРОДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК ВОЗДУШНО-КОСМИЧЕСКИХ АППАРАТОВ В ПЕРЕХОДНОМ РЕЖИМЕ

Зея Мьо Мьинт 1
1 ФГОУ ВПО «Московский физико-технический институт (государственный университет)»
Развитие космической техники и высотной гиперзвуковой авиации требует развития методов определения аэродинамических характеристик (АДХ) во всем диапазоне режимов течения – от орбитального полета до посадочного режима. Экспериментальное исследование высокоскоростных разреженных течений довольно проблематично, и в настоящее время практически единственным средством получения информации об аэродинамической обстановке около воздушно-космического аппарата (ВКА) на больших высотах полета являются методы вычислительной аэродинамики. Математическое моделирование аэродинамики ВКА в переходном режиме является весьма сложным и требует больших вычислительных и временных ресурсов. В настоящей работе представлены методы определения АДХ ВКА в переходном режиме и результаты расчётов аэродинамических характеристик воздушно-космического аппарата и гиперзвукового летательного аппарата с помощью локально-инженерного метода.
воздушно-космический аппарат
аэродинамика в переходном режиме
компьютерное моделирование
локально-инженерный метод
гиперзвуковая технология
1. Белоцерковский О.М., Хлопков Ю.И. Методы Монте-Карло в механике жидкости и газа. – М.: Азбука, 2008. – 330 с.
2. Ваганов А.В., Дроздов С.М., Дудин Г.Н., Косых А.П., Нерсесов Г.Г., Пафнутьев В.В., Челышева И.Ф., Юмашев В.Л. Численное исследование аэродинамики перспективного возвращаемого космического аппарата // Ученые записки ЦАГИ. – 2007. – Т. XXXVIII, № 1–2. – C. 16–26.
3. Галкин В.С., Ерофеев А.И., Толстых А.И. Приближенный метод расчета аэродинамических характеристик тел в гиперзвуковом разреженном газе. // Труды ЦАГИ. – 1977. – Вып. 1833.
4. Зея Мьо Мьинт, Хлопков А.Ю. Аэродинамические характеристики летательного аппарата сложной формы с учётом потенциала взаимодействия молекулярного потока с поверхностью// Ученые записки ЦАГИ. – 2010. – Т. XLI, № 5. – С. 33–45.
5. Коган М.Н. Динамика разреженных газов. – М.: Наука, 1967. – 440 с.
6. Хлопков Ю.И. Статистическое моделирование в вычислительной аэродинамике. – М.: МФТИ, 2006. – 260 с.
7. Хлопков Ю.И., Зея Мьо Мьинт, Хлопков А.Ю., Чжо Зин, Засыпалов В.В. Исследование аэродинамики перспективных гиперзвуковых летательных аппаратов // Материалы международной научно-практической конференции «Актуальные вопросы и тенденции развития биологии, химии, физики». – Новосибирск, 2012. – С. 98–103.
8. Хлопков Ю.И., Чернышев С.Л., Зея Мьо Мьинт, Хлопков А.Ю. Введение в специальность II. Высокоскоростные летательные аппараты. – М.: МФТИ, 2013. –192 с.
9. Belotserkovskii O.M., Khlopkov Y.I. Monte Carlo Methods in Mechanics of Fluid and Gas // World Scientific Publishing Co. Ltd. N-Y, London, Singapore, Beijing, Hong Kong, 2010. – 268 p.

Развитие космической техники и высотной гиперзвуковой авиации требует развития методов определения аэродинамических характеристик (АДХ) во всем диапазоне режимов течения – от орбитального полета до посадочного режима. Обтекание космического аппарата (КА) на высотах ниже 60–70 км происходит в условиях сплошной среды. Числа Кнудсена в этих условиях достаточно малы (Kn << 1). Для исследования сплошносредных течений используются континуальные методы. Между этими предельными режимами КА проходит переходной режим обтекания, когда необходимо учитывать как столкновения молекул набегающего потока с поверхностью, так и межмолекулярные столкновения. Переходной режим характеризуется числами Кнудсена: 0,001 ≤ Kn ≤ 10. В этих условиях сплошносредные методы исследования неприменимы вследствие высокой разреженности и термохимической неравновесности газа [6]. Экспериментальное моделирование высокоскоростных разреженных течений довольно проблематично, и в настоящее время практически единственным средством получения информации об аэродинамической обстановке около КА на больших высотах полета являются методы вычислительной аэродинамики [1, 9]. Для анализа аэродинамических характеристик КА в переходном режиме необходимо использовать кинетический подход (решение уравнения Больцмана или его моделей). Математическое моделирование аэродинамики КА в переходном режиме является весьма сложным и требует больших вычислительных и временных ресурсов.

Компьютерное моделирование позволяет при помощи физических моделей и математических методов быстро проводить анализ аэродинамических характеристик гиперзвуковых летательных аппаратов. В настоящее время существует несколько подходов решения аэродинамических характеристик гиперзвуковых летательных аппаратов, и также проведены многочисленные исследования аэродинамических характеристик космических аппаратов вдоль всей траектории [4, 6, 7]. Однако метод Монте-Карло обладает достаточно хорошей точностью, но требует большого времени для вычисления. Другие основы на упрощенных инженерных методиках требуют малых затрат расчетного времени. Целью настоящей работы является анализ методов моделирования аэродинамических характеристик воздушно-космических аппаратов в переходном режиме.

Методы расчета аэродинамических характеристик тел в переходном режиме

В настоящее время условно можно выделить два инженерных подхода к вычислению аэродинамических характеристик по числам Рейнольдса. Первый подход состоит в построении функции аппроксимации при известных предельных значениях: свободномолекулярного С(0) и сплошносредного, обычно моделируемого по методу Ньютона С(∞).

zeya01.wmf

Функция f зависит от свойств газа, параметров набегающего потока, геометрии поверхности и др. В данной работе используется классический метод локальности и предполагается [4]

zeya02.wmf

zeya03.wmf

zeya04.wmf

zeya05.wmf

В предельном случае сплошной среды по методу Ньютона получаем

zeya06.wmf

В другом предельном свободномолекуляном случае получаем

zeya07.wmf

В данной работе используются выражения для элементарных сил давления и трения в форме работы [3].

zeya08.wmf

zeya09.wmf

Здесь коэффициенты p0, p1, t0 (коэффициенты режима течения) зависят от числа Рейнольдса Re0 = r∞u∞L/m0, температурного фактора tw = Tw/T0, коэффициентов аккомодации и отношения удельных теплоемкостей c, L – характерный размер, m = m(T0) – коэффициент вязкости, T0, Tw – температура торможения и температура поверхности.

Зависимость коэффициентов режима в гиперзвуковом случае должна обеспечивать переход к свободномолекулярным значениям при Re0 → 0 и значениям теории Ньютона, методов тонких касательных клиньев или конусов при Re0 → ∞. На основе анализа расчетных и экспериментальных данных предложены эмпирические формулы

zeya10.wmf

zeya11.wmf

zeya12.wmf

здесь

zeya13.wmf

zeya14.wmf

zeya15.wmf zeya16.wmf

где h – относительные поперечные размеры аппарата, равные отношению его высоты к длине.

Предложенная методика хорошо зарекомендовала себя для расчета гиперзвукового обтекания выпуклых не очень тонких и пространственных тел. Расчет полностью отражает качественное поведение Сx в зависимости от разреженности среды во всем диапазоне углов атаки и дает количественное соответствие с экспериментом и расчетом по уравнению Больцмана с точностью около 5 %.

О точности соотношений локального метода можно сказать следующее. Ясно, что они применимы с наименьшей погрешностью в случае тел, близких к сфере, и неприменимы в случае очень тонких тел, когда не выполняется условие M∞ sin q >> 1. В рассматриваемых методах не учитывается влияние взаимодействия пограничного слоя с гиперзвуковым невязким потоком при больших числах Rе0. Расчетные и экспериментальные значения Cх конуса в переходном режиме согласуются удовлетворительно, данные по Cy согласуются значительно хуже. Необходимо подчеркнуть, что предложенная методика качественно верно отражает немонотонность зависимости Cy конуса от Rе0. Расчетные и экспериментальные результаты по Cх при a = 10° и 15° для пластины хорошо согласуются, данные же для Cх при a = 5° и Cy согласуются плохо. Это является следствием неучтенного в локальном методе влияния взаимодействия пограничного слоя с невязким потоком.

Таким образом, локальный метод расчета аэродинамических характеристик тел в гиперзвуковом потоке разреженного газа в переходном режиме дает хороший результат по Cх для широкого класса тел и качественно верный результат по Cy. При малых углах атаки (a < 5°) точность результата ухудшается, в этом случае необходимо привлекать более полные модели, учитывающие наличие пограничного слоя [1, 3, 6, 8].

Результаты исследования и их обсуждения

На рисунках представлены зависимости коэффициентов сил сопротивления Cx, подъемной Cy и момента тангажа mz от угла атаки (a = 0–90°) при различных значениях числа Рейнольдса Re. Параметры задачи были следующие: отношение теплоемкостей g = 1,4, температурный фактор tw = Tw/T0 = 0,4, скоростное отношение s = 10, число Рейнольдса Rе0 = 0, 10, 102, 104.

На рис. 2–4 представлены сравнения результатов зависимости Cx(a), Cy(a), mz(a) для воздушно-космического аппарата (ВКА) «Клипер», модель ЦАГИ [2], и гиперзвукового летательного аппарата (ГЛА) «Falcon HTV-2» (рис. 1) с помощью локального метода.

pic_23.tif pic_24.tif

Рис. 1. Расчетные схемы ВКА «Клипер» и ГЛА «Falcon HTV-2»

pic_25.tif

Рис. 2. Зависимость Cx(a) для «Клипер» и «Falcon HTV-2»

pic_26.tif

Рис. 3. Зависимость Cy(a) для «Клипер» и «Falcon HTV-2»

pic_27.tif

Рис. 4. Зависимость mz(a) для «Клипер» и «Falcon HTV-2»

Из этих результатов видно, что с увеличением числа Рейнольдса коэффициент сопротивления тела уменьшается (что можно объяснить уменьшением нормальных и касательных напряжений p1(Re0) и t0(Rе0)), при этом общий характер зависимости Cx(a) не изменяется. Зависимость Cy(a) является несимметричной при Re0 → ∞, так что значение Cy при положительных углах атаки существенно больше по модулю Cy при отрицательных углах атаки. Из этих результатов чувствительно, что коэффициенты силы сопротивления «Falcon HTV-2» меньше чем, «Клипера», и изменение числа Рейнольдса оказывает сильное влияние на все аэродинамические характеристики. Можно сказать, что число Рейнольдса влияет формы тела, и локальный метод дает хорошие результаты в переходном режиме для широкого класса тел.

Выводы

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

Работа выполнена при поддержке РФФИ (Грант № 14-07-00564-а).

Рецензенты:

Боголепов В.В., д.ф.-м.н., главный научный сотрудник ФГУП «ЦАГИ», профессор МФТИ, г. Жуковский;

Липатов И.И., д.ф.-м.н., начальник отдела ФГУП «ЦАГИ», профессор МФТИ, г. Жуковский.

Работа поступила в редакцию 11.02.2014.


Библиографическая ссылка

Зея Мьо Мьинт АНАЛИЗ МЕТОДОВ ОПРЕДЕЛЕНИЯ АЭРОДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК ВОЗДУШНО-КОСМИЧЕСКИХ АППАРАТОВ В ПЕРЕХОДНОМ РЕЖИМЕ // Фундаментальные исследования. – 2014. – № 3-3. – С. 495-499;
URL: http://www.fundamental-research.ru/ru/article/view?id=33702 (дата обращения: 09.12.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074