Научный журнал
Современные наукоемкие технологии
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,858

МЕТОДИКА ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ КОЭФФИЦИЕНТА ТЕПЛОПЕРЕДАЧИ НА ТЕЛАХ СЛОЖНОЙ ФОРМЫ В ПЕРЕХОДНОМ РЕЖИМЕ

Зея Мьо Мьинт 1
1 Московский физико-технический институт
1. Авдуевский В.С., Галицейский Б.М. и др. Основы теплопередачи в авиационной и ракетно-космической технике. – М.: Машиностроение, 1992. – 528 с.
2. Белоцерковский О.М., Хлопков Ю.И. Методы Монте-Карло в механике жидкости и газа. – М.: Азбука, 2008. – 330 с.
3. Хлопков Ю.И., Зея Мьо Мьинт, Хлопков А.Ю., Чжо Зин Методы Монте-Карло для определения аэротермодинамических характеристик гиперзвуковых воздушно космических систем // Materials digest of LI International Research and Practice Conference «Physical, Mathematical and Chemical Sciences: Theoretical Trends and Applied Studies», – London: IASHE, 2013, pp. 41-44.
4. Коган М.Н. Динамика разреженного газа. – М.: Наука, 1967. – 440 с.
5. Lee Lester Laminar heat transfer over blunt nosed bodies at hypersonic flight speeds // Jet Propulsion. – 1956. – V. 26, N 4. – P. 259–269.
6. Ващенков П.В. Численный анализ высотной аэротермодинамики космических аппаратов: дис. канд.–техн. наук, Новосибирск: ИТПМ СО РАН., 2012.
7. Dogra V.K., Wilmoth R.G., Moss J.N. Aerothermodynamics of 1.6 -m- diameter sphere in hypersonic rarefied flow. J. AIAA. – 1992. – V. 30, N 7. – P. 1789 – 1794.
8. Зея Мьо Мьинт Высотная аэротермодинамика гиперзвуковых воздушно-космических систем // Материалы II международной научно-практической конференции «Сайнс-проекты молодых ученых – 2013». – Ростов-на-Дону, 2013. с. 27-30.
9. Зея Мьо Мьинт, Хлопков А.Ю. Исследование аэротермодинамики перспективных гиперзвуковых летательных аппаратов // Труды МАИ. 2013, № 66, 19 с.
10. Зея Мьо Мьинт, Хлопков А.Ю. Методика расчета тепловых потоков в ламинарном и турбулентном пограничном слое // Труды 56-й научной конференции МФТИ «Современные проблемы фундаментальных и прикладных наук». – Жуковский, 2013. с. 28-29.
11. Хлопков Ю.И., Чернышев С.Л., Зея Мьо Мьинт, Хлопков А.Ю. Введение в специальность II. Высокоскоростные летательные аппараты.– М.: МФТИ, 2013. – 192 с.

При создании и проектировании гиперзвукового летательного аппарата различного назначения необходимо детальное значение их аэротермодинамических характеристик вдоль всей траектории полета. В процессе исследования тепловых нагрузок, действующих на поверхность космических аппаратов, важным этапом является решение задачи создания их тепловой защиты и определения температурных режимов конструкции. В настоящее время существует несколько подходов решения аэротермодинамических характеристик гиперзвуковых летательных аппаратов, также проведены многочисленные исследования аэродинамических характеристик космических аппаратов вдоль всей траектории – от орбитального полета до посадочного режима [1-3]. Однако обладают достаточно хорошей точностью, но требуют большого времени для вычисления. Другие основы на упрошенных инженерных методиках требуют малых затрат расчетного времени, но специфика существующих алгоритмов быстрого счёта позволяет оценивать тепловые потоки на телах достаточно простой формы.

Для вычисления коэффициента теплопередачи Ch на элементарную площадку в свободномолекулярном пределе используются аналитические формулы в виде [4]:

fizmat49.wmf,

где

fizmat50.wmf,

fizmat51.wmf,

αe – коэффициент аккомодации энергии на стенке, s∝ – отношение скорости набегающего потока к наиболее вероятной скорости молекул, Tw – температура стенки, T∝ – температура набегающего потока, g – показатель адиабаты, s∝,q = s∝ cos q.

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

Для вычисления коэффициента теплопередачи Ch в континуальном режиме будем использовать методику, основанную на теории Лиса [5].

Коэффициент теплопередачи в произвольной точке тела вычисляется по формуле

fizmat52.wmf.

Здесь s – расстояние вдоль линии тока от точки торможения до рассматриваемой элементарной площадки, q – угол между направлением потока и нормалью к элементарной площадке в данной точке, Ch0 – коэффициент теплопередачи в точке торможения:

fizmat53.wmf.

Здесь k = 1 для сферической точки торможения, k = 0 для цилиндрической точки торможения, r – радиус кривизны поверхности в точке торможения, w – показатель степени в степенной зависимости вязкости от температуры, Pr = mCp/c – число Прандтля, число Рейнольдса Re∝,r вычислено по параметрам набегающего потока и радиусу кривизны в точке торможения. Число Рейнольдса в континуальном и около континуальном режимах, так как при приближении к свободномолекулярному режиму, когда Re → 0, величина Ch0 → 0.

В настоящей работе предлагается локально-мостовой метод вычисления коэффициента теплопередачи на элементах выпуклой поверхности с учетом расстояния от точки торможения в переходном режиме. Локально-мостовой метод позволяет быстро получить аэротермодинамические характеристики при проведении большого количества многовариантных расчетов [6]:

fizmat54.wmf,

fizmat55.wmf.

Здесь M – число Маха, Re – число Рейнольдса, S – площадь поверхности тела. Функция Fb называется мостовой функцией. Рассмотрим мостовую функцию, выражающуюся как функция ошибки от логарифма числа Кнудсена:

fizmat56.wmf,

fizmat57.wmf.

Если Kn0 < Knm, используется мостовая функция Fb,1. В противном случае Fb,2. Значения Knm = 0.3, ΔKn1 = 1.3 и ΔKn2 = 1.4 были определены путем сравнения с результатами моделирования методом прямого статистического моделирования (Монте-Карло). Разработанная методика для расчета тепловых потоков на тело имеет практический интерес для организации и специалистов, занимающихся определением теплового потока.

На рисунке показаны зависимости коэффициента теплопередачи на сфере от угла при различных числах Рейнольдса (Re = 0.1, 3.8, 3243). Полученные результаты сравнены с результатами Dogra V.K., Wilmoth R.G. [7] и Ващенкова П.В. [6]. Здесь угол откладывается от точки торможения.

fm8.tif

Угловое распределение Ch на сфере

На рисунке видно, что тепловой поток в свободномолекулярном и в континуальном режимах локально-мостовым методом определяется хорошо, но в переходном режиме дает повышенное значение коэффициента примерно на 15 % на углах 20–55 градусов [6, 10]. Можно сказать, что результаты, полученные локально-мостовым методом, дают правильное поведение и достаточно хорошо согласуются с результатами метода прямого статистического моделирования (ПСМ). Локально-мостовой метод дает в переходном режиме завышенное значение коэффициента теплопередачи примерно на 15 % на углах 20–50 градусов. Методики и результаты смогут быть полезны при создании современных и перспективных воздушно-космических аппаратов нового поколения [8-11]. Работа выполнена при поддержке РФФИ (Грант № 14-07-00564-а).


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

Зея Мьо Мьинт МЕТОДИКА ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ КОЭФФИЦИЕНТА ТЕПЛОПЕРЕДАЧИ НА ТЕЛАХ СЛОЖНОЙ ФОРМЫ В ПЕРЕХОДНОМ РЕЖИМЕ // Современные наукоемкие технологии. – 2014. – № 5-2. – С. 115-116;
URL: http://top-technologies.ru/ru/article/view?id=34011 (дата обращения: 22.10.2018).

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

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