VI Международная студенческая научная конференция Студенческий научный форум - 2014

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

Этот метод основан на использование особенностей развития пограничного слоя, заключающихся в том, что в случае ускоренных течений теплообмен в рассматриваемом сечении определяется в основном параметрами потока и толщиной пограничного слоя в этом сечении. В значительно меньшей степени тепловой поток зависит от условий, в которых пограничный слой развивался от точки его образования до рассматриваемого сечения. В соответствии с этими особенностями, метод эффективной длины состоит в том, что при расчете теплообмена действительное сечение заменяется течением над пластиной (для осесимметричного тела – над цилиндром) с параметрами потока, равными параметрам рассматриваемого сечения. Длина пластины или цилиндра x_эф выбирается из условия нарастания на ней теплового пограничного слоя толщиной, равной толщине слоя в рассчитываемого сечения тела.

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

Формулы для расчета теплообмена при ламинарном пограничном слое

Тепловой поток определяется по формуле Ньютона

где

Для воздуха при ламинарном пограничном слое r = 0.84.

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

Коэффициент K учитывает влияние сжимаемости. Используя это уравнение и метод эффективной длины для расчета теплообмена при произвольном распределении параметров потока вдоль образующей тела, будем иметь

где K₁ – поправка на влияние продольного градиента давления

,

В частном случае T_w = const, приближенное выражение для определения эффективной длины имеет вид

Под интегралом стоят переменные величины R(x), _w, u₁ изменяющиеся от начала образования пограничного слоя (критическая точка) до рассматриваемого сечения. В знаменателе R(x), _w, u₁ соответственно радиус вращения, плотность и скорость в рассматриваемом сечении.

Учитывая постоянство T_w, плотность газа _w можно определить из соотношения

где _w0и p₀₁ параметры воздуха в критической точке.

Поправка на сжимаемость рассчитывается по формуле

где _w и _w определяются по T_w, ₁, ₁ по T₁, а ** по максимальной температуре T_max = T* в пограничном слое.

При

T* определяется по формуле

При

Максимальная температура в пограничном слое равна температуре внешнего потока T* = T₁. В этом случае

Поправка на влияние градиента скорости

где m – показатель степени в выражении u₁ = x^m,

,

При течении в окрестности критической точки m = 1 и

В случае обтекания конуса m = 0, K₁ = 1.

Определив критерий Нуссельта Nu_wэф и можно найти коэффициент теплоотдачи

.

Формулы для расчета теплообмена при турбулентном пограничном слое

Тепловой поток определяется по формуле Ньютона. Эффективная температура имеет вид

где r = 0.89.

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

K_T – поправка на сжимаемость. Здесь

,

Эффективная длина пластины при турбулентном течении в пограничном слое и T_w = const определится из выражения

Из формулы получено в предположении постоянства по x выражения –

В окрестности критической точки изменение этих величин невелико и переменностью их по x можно пренебречь.

Где скорость потока невелика, хорошие результаты дает поправка на сжимаемость, определенная как

При большой скорости потока

Коэффициент теплоотдачи

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

Рис.1. Распределение тепловых потоков на сфере

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

Список литературы

1. Хлопков Ю.И. Статистическое моделирование в вычислительной аэродинамике. М.: МФТИ, 2006. – 158 с.

2. Авдуевский В.С., Галицейский Б.М. и др. Основы теплопередачи в авиационной и ракетно-космической технике. – М.: Машиностроение, 1992. – 528 с.

3. Журин С.В. Методика численного моделирования конвективного теплообмена на телах сложной формы с использованием метода эффективной длины // Дис. канд. физ-мат. наук, МФТИ, - Москва, 2009.

4. Ващенков П.В. Численный анализ высотной аэротермодинамики космических аппаратов // Дис. канд.-тех. наук, ИТПМ СО РАН. - Новосибирск, 2012.

5. Хлопков Ю. И., Чернышев С.Л., Зея Мьо Мьинт, Хлопков А. Ю. Введение в специальность II. Высокоскоростные летательные аппараты.– М.: МФТИ, 2013. – 192 с.

6. Зея Мьо Мьинт, Хлопков А.Ю. Исследование аэротермодинамики перспективных гиперзвуковых летательных аппаратов // Труды МАИ. - 2013, № 66. - 19 с.

7. Зея Мьо Мьинт, Хлопков А.Ю. Методика расчета тепловых потоков в ламинарном и турбулентном пограничном слое // Труды 56-й научной конференции МФТИ «Современные проблемы фундаментальных и прикладных наук». - Жуковский, 2013. с. 28-29.

Код для цитирования: Скопировать