ТЕПЛО - И МАССООБМЕН НА ПРОНИЦАЕМОЙ ПОВЕРХНОСТИ СИСТЕМЫ ТЕПЛОВОЙ ЗАЩИТЫ СПУСКАЕМОГО КОСМИЧЕСКОГО АППАРАТА Текст научной статьи по специальности «Механика и машиностроение»

Электронный журнал «Труды МАИ». Выпуск № 50

www.mai.ru/science/trudy/

УДК 629.7.063 (075.8)

Тепло - и массообмен на проницаемой поверхности системы тепловой защиты спускаемого космического аппарата малой формы1.

Федорченко Е.А., Никитин П.В.,

Аннотация. Представлены результаты анализа тепло-массообмена при разрушении поверхностного слоя тепловой защиты спускаемого космического аппарата малой формы, совершающего баллистический спуск на Землю. Выполнен расчёт траекторных параметров и теплообмена на поверхности спускаемого космического аппарата на гиперзвуковом участке полёта. По результатам расчёта проведена серия экспериментов на высокотемпературном газодинамическом стенде с целью определения теплозащитных характеристик ряда ТЗМ. Выбран теплозащитный материал поверхностного слоя системы тепловой защиты КЛА. Проведена оценка толщины разрушения данного слоя. Ключевые слова: тепло-массообмен, гиперзвуковая скорость, тепловая защита, траекторные параметры, энтальпия торможения, орбитальный полёт, баллистический спуск, газодинамический стенд.

Используемые сокращения:

КА - космический аппарат;

КЛА - космический летательный аппарат;

СКА - спускаемый космический аппарат;

ТЗМ - теплозащитный материал;

ТЗ - тепловая защита,

ОКС - орбитальная космическая станция,

ПКТ - передняя критическая точка

Введение:

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

Работа выполнена при поддержке РФФИ, гранты № 09-08-00169-а и 11-08-00828-а.

1

плане одним из таких проектов, имеющих большое прикладное значение, является оперативная доставка на Землю с орбитальных космических станций (ОКС) информации, свободной от помех и искажений, возникающих при использовании телеметрических систем.

Успешное решение такой задачи возможно благодаря использованию специальных КЛА малых форм, отстреливаемых с борта ОКС и совершающих для посадки на Землю баллистический полёт в плотных слоях атмосферы с гиперзвуковой скоростью.

Как известно, /1/, /2/ прохождение плотных слоев атмосферы по баллистической траектории вызывает экстремальные тепловые и газодинамические нагрузки. Например, полёт таких КЛА в атмосфере с гиперзвуковой скоростью устанавливает уровень интенсивных тепловых потоков, воздействующих на конструкцию аппарата и реализующих его тепловой режим (аэродинамический нагрев). Понятно, что научное прогнозирование тепловых и газодинамических нагрузок на аппарат, умение рассчитать его тепловой режим, определяет надёжность выполнения поставленной перед аппаратом программы на завершающей её стадии - возвращении на Землю.

Аэродинамика обтекания гиперзвуковым потоком формирует траекторию полёта аппарата, устойчивость его движения, т.е., в конечном счёте, определяет координаты посадки аппарата в расчётном заданном районе.

Успешное выполнение указанных задач однозначно связано с сохранением изначальной формы (геометрии) СКА. Последнее возможно только при создании надёжной и оптимальной по массе и характеристикам системы тепловой защиты аппарата.

Проектирование такой тепловой защиты для СКА малой формы являлось целью дипломной работы.

В работе в качестве прототипа формы спускаемого аппарата была выбрана оболочка космического аппарата «Экспресс», который разработан в ГКНПЦ им. М.В. Хруничева (рис. 1).

Для выполнения указанной цели в работе решены следующие задачи: - проведен расчёт траекторных параметров орбитального космического аппарата выбранной формы, совершающего баллистический спуск в атмосфере Земли с гиперзвуковой скоростью;

- проанализированы физико-химические особенности механизма тепло - и массообмена на непроницаемой поверхности КЛА при обтекании гиперзвуковым потоком2;

- составлен алгоритм и проведён расчёт конвективного теплообмена на непроницаемой поверхности КЛА при полёте в плотных слоях атмосферы;

- выбран теплозащитный материал (ТЗМ) поверхностного слоя системы тепловой защиты и определены экспериментально его теплозащитные характеристики;

- составлен алгоритм и проведён расчёт теплообмена на проницаемой поверхности СКА при термохимическом разрушении выбранного полимерного теплозащитного материала;

- проведен анализ результатов расчёта теплообмена на непроницаемой и проницаемой поверхности КЛА;

- на основе основных положений теории подобия рассчитана и спроектирована физическая модель КЛА для проведения экспериментальных исследований теплообмена в лабораторных условиях на высокотемпературных газодинамических стендах кафедры «Авиационно-космической теплотехники» МАИ;

- проведен расчёт толщины разрушения (уноса) поверхностного слоя выбранного ТЗМ за время полёта СКА по баллистической траектории.

1. Расчёт траекторных параметров орбитального СКА при баллистическом полёте в атмосфере Земли.

Постановка задачи. Рассматривается вход космического спускаемого аппарата малой формы (рис. 1) с околоземной круговой орбиты в плотные слои атмосферы Земли. Полёт в атмосфере осуществлялся по баллистической траектории в поле действия газодинамических сил набегающего потока и сил тяготения Земли. Расчёт траекторных параметров проводился, начиная с высоты 90 км, наиболее плотных слоях атмосферы.

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

2 В процессе теплообмена поверхность называют непроницаемой, если на ней отсутствуют источники (стоки) массы. В противном случае, поверхность именуется проницаемой.

части замыкается усеченным конусом. Общая масса спускаемого аппарата 420 кг. Габаритные размеры спускаемого КА 1310х1010мм. Угол конусности у =15°.

Рисунок 1. Форма спускаемого космического аппарата «Экспресс». Расчёт траекторных параметров.

Скорость полёта во времени У(т). Расчёт скорости полёта по баллистической траектории проведён с использованием известного соотношения для входа КЛА с круговой орбиты в атмосферу планеты /1/:

V(г) = • exp

■ g 'Ph

H

-p-H (г)

(1)

2-p-sin 0ет j

где Увх - скорость орбитального полёта КЛА, Н(т) - изменение высоты полёта во времени,

двх - угол входа в атмосферу с орбитальной траектории,

C 'F

х мид

баллистический

т' ё

коэффициент, р - коэффициент атмосферы планеты, для Земли р = 0,142, [1/км]. Изменение во времени энтальпии торможения /р(г). Изменение энтальпии торможения по траектории спуска рассчитывалась с использованием известного из термогазодинамики выражения:

К2 М

1о (г) = Ih (г) + -

2

(2)

где 1Н(т) - изменение энтальпии атмосферы по высоте полёта Н(т), Уад(т) - изменение во времени скорости СКА на траектории баллистического спуска (1).

Изменение во времени давления торможения ро(т). Расчёт изменения давления торможения р0(т) по траектории спуска проводился с использованием уравнения Гюгонио в виде /3/:

Ро (0 = Рн (Т)'

к+1

'к + 1 ^ к-1

1

\к—1

[МН М]

2к к-1

к - 1

2к

1

к-1

(3)

к -~1 ■м «Г2 -1]

где рн(т) - атмосферное давление на высоте подлёта в данный момент времени (таблица стандартной атмосферы), Мн(т) - число Маха, соответствующее скорости полёта КЛА на данной высоте, к = Ср/Су.

Изменение температуры торможения То(т). Изменение температуры торможения определялось с использованием рассчитанных ранее значений энтальпии торможения 10(т) и давления торможения Р0(т), с последующим применением справочных данных в виде зависимости /0 = / (р0, Т0) /9/.

Результаты расчёта изменения указанных траекторных параметров во времени представлены на рисунках.

Рисунок 2. Изменение высоты и скорости спуска по времени

1 - изменение скорости, 2 - изменение высоты

Рисунок 3. Изменение энтальпии торможения и давления торможения по времени спуска.

1 - изменение энтальпии торможения, 2 - изменение давления торможения.

2. Расчет теплообмена на непроницаемой поверхности

Аэродинамический нагрев в сочетании с газодинамическим напором являются наиболее мощными факторами воздействия на конструкцию КЛА. Причина их возникновения связана термо-газодинамикой обтекания поверхности аппарата набегающим вязким гиперзвуковым потоком.

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

Расчёт задачи конвективного теплообмена в таких условиях сводится к решению многопараметрической математической модели, включающей систему нелинейных дифференциальных уравнений в частных производных /1/, /4/.

Решение подобных математических моделей проводится численными методами с привлечением мощного математического программного обеспечения /5/, /6/. Например, в последние годы широкое применение нашёл программный комплекс ANSYS /7/. При решении задач, связанных с турбулентными течениями широко используется - е» -модель.

Тем не менее, даже при современной компьютерной технике и программным

обеспечением решение подобных задач численными методами требует громадных затрат

временного ресурса. В этой связи в дипломной работе для расчёта теплообмена на

6

непроницаемой поверхности, обтекаемой гиперзвуковым потоком, использовался наиболее простой инженерный метод, основанный на применении алгебраических критериальных соотношений Фея - Ридделла /2/.

На базе таких соотношений составлен алгоритм решения задачи, при этом процесс решения осуществлялся с привлечением программы Excel - 95/97, а также новой версии Excel 2000.

Ниже приводится составленный алгоритм расчёта задачи. 2.1. Алгоритм решения задачи.

Для проведения расчёта теплообмена на поверхности СКА его поверхность разбивается на расчетные участки. На рисунке 4 представлена расчётная схема спускаемого аппарата малой формы3.

Расчет теплообмена проводится по составленному ниже алгоритму для каждого i-го расчётного сечения с учётом ламинарного или турбулентного режимов течения.

Рисунок 4. Расчетная схема спускаемого КА малой формы.

На рисунке: 0 - передняя критическая точка, г = 1,.. .10 номера расчётных сечений СКА, у =15о - угол полуконусности.

Алгоритм расчёта конвективного теплообмена на поверхности притупленного тела сводится к следующим этапам.

Проводится расчет теплообмена: в передней критической точке, в области передней критической точки и вне области передней критической точки (боковая поверхность КЛА) /1/.

3 С целью сокращения времени расчёта в работе выбрано всего лишь 11 расчётных сечений на поверхности СКА.

2.2.1. Расчет теплообмена в передней критической точке.

Тепловой поток в передней критической точке притупленного сферой тела рассчитывается с использованием критериального соотношения Фея - Ридделла в виде:

д0 = 0,763 • Рг^-0'6 • (р0 • / )°'4 • (р • / )'

V 0,1

Ыи ^

у J

1 + (Ьв" -1)

I

Б

I„

(1° - 1ш ) , (4)

где показатель п критерия Ье выбирается в зависимости от состояния пограничного слоя

/1/;

- энтальпия диссоциации молекул, входящих в состав многокомпонентной

высокотемпературной смеси, рассчитывалась по соотношению:

"

1б = £[ I 0,-(С° - с,ж)]

(5)

где 1г - энтальпия образования атомов, Сг,0 и С^ - массовые концентрации атомов на внешней границе пограничного слоя при температуре торможения (индекс "0") и на непроницаемой поверхности при температуре Тк (индекс "),

- градиент скорости:

(2 (р° - р— )

ёх в

р°

(6)

где

р

коэффициент С = 1 —— , радиус притупления сферы Я0, р0 -плотность при \| Рг\

температуре торможения Т0 и давлении торможения Р0 в передней критической точке (рис. 4), РН - давление атмосферы на высоте полета;

- критерий Льюиса-Семенова для многокомпонентной высокотемпературной газовой смеси в виде:

Ьв =

р0 • Б12 • СР

Л

(7)

Рассчитывался по параметрам торможения для диссоциированного воздуха; - изобарическая теплоемкость смеси рассчитывалась по формуле:

С = £[ с, • ср4 ]

(8),

плотность смеси:

1=1

1=1

Ро =

р

Ксм ■ Т0

(9),

Ясм - газовая постоянная смеси;

- коэффициент теплопроводности смеси, рассчитывался по формуле /8/ или определялся с использованием справочных данных /9/:

Л,. = 0,5

/ Л-1

X

IX, Л + ^

1—1

(10)

и

где Хг - мольная доля ¡-го компонента, X = С —, И и /исм - мольные массы ¡-ой

и

компоненты и смеси, соответственно, и и и - коэффициенты динамической вязкости смеси в выражении (4). Определяются из таблиц /9/ для многокомпонентной

высокотемпературной смеси при давлении, Р0 и температуре Т0 торможения и

температуры поверхности Т^, соответственно; - коэффициент бинарной диффузии смеси /8/:

Д2 = 2,6■Ю-7 ■ Т2 ■

и+и и ■ и

0,5

Р ■а2 ■О1'1

1 о а 1,2 О 1,2

м

сек.

(11)

При расчёте принималось, что при высокой температуре газа в пограничном слое его атомы и молекулы при столкновении взаимодействуют, как абсолютно упругие шары. В таком случае, выражение в знаменателе формулы (11) упрощается, в виду того, что:

а 1,2 ■О1 1,2 — и 1,2, где, в свою очередь, и 1,2

М + и2 2

и ё2 - диаметры атомов и молекул, в ангстремах, соответственно (справочные данные /1/).

2.2.2. Расчет теплообмена в области передней критической точки \1\.

В этой области течение развивается из передней критической точки и, постоянно ускоряясь, распространяется вдоль поверхности тела по оси х к его периферии до «звуковой линии» (см. рис. 4). На этом участке скорость течения прямо пропорциональна

координате х по поверхности КЛА, т.е и (х) = /3 ■ х . Где / = — - градиент скорости на

ах

поверхности, рассчитывается по формуле (6).

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

Расчет теплообмена при ламинарном режиме течения.

Плотность теплового потока рассчитывалась с использованием

модифицированного уравнения конвективного теплообмена Ньютона в виде:

Яо,,

/ л а

С

V СР Уо

(1 е,1 1Ш,1 ) ;

(12)

где 1е - энтальпия восстановления набегающего потока /1/, - энтальпия газа

набегающего потока при температуре стенки Т^.

Таким образом, как следует из (12), для того, чтобы рассчитать плотность теплового потока, qo,i в любом расчётном сечении поверхности КЛА, необходимо уметь

определять значение параметра теплообмена

а

С

V СР У о

в каждом сечении. В работе для

расчёта указанного параметра использовалось соотношение Фея - Ридделла в виде:

/ Л

а

С V СР Уо

0,763 ■ РгШ0,6-

хэф

г

Л

0,4

Ре -Я

V рш ■ У

1 + (Ье" -1)

I

Б

I

(13),

В критерии Рейнольдса характерный размер хэф рассчитывается с использованием интегрального соотношения:

} Р „-и(х )■ я2 (х )■ ах

х эф. = '

(14).

р ^и(х )■ Я2 (х)

Таким образом, рассчитав с помощью (13) параметр теплообмена, подставив его значение в выражение (12), получают значение плотности теплового

а

С

V СР Уо

и,

о

потока, q0,i на любом расчетном бесконечно малом участке, принадлежащем области передней критической точки сферы.

Данный алгоритм можно было бы использовать в инженерной практике и для расчета теплообмена в любой вызывающей интерес зоне на поверхности летательных аппаратов произвольной формы. Однако, не смотря на кажущуюся простоту алгоритма, такой подход не всегда возможен из-за ряда трудностей, обусловленных тем, что в критериальном соотношении (12) не известен ряд определяющих параметров. Например, не всегда возможно установить аналитически закон изменения скорости потока, u(x) на произвольной криволинейной поверхности летательных аппаратов. Не поддается также прямой оценке закон изменения температуры Тцт^) по поверхности и др. Возникают также трудности при определении коэффициентов динамической вязкости, ие и и от

температуры, особенно при высоких уровнях температуры восстановления Те. Кроме того, из-за обилия различных формул возникают значительные несоответствия при расчете коэффициентов бинарной диффузии Dij.

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

Расчет теплообмена при турбулентном режиме течения.

При реальных полётах в атмосфере космических летательных аппаратов с гиперзвуковой скоростью на его поверхности возникают течения, характеризуемые громадными значениями чисел Рейнольдса. Такие течения обладают особыми свойствами, известными как "турбулентность". Термо-газодинамика и теплообмен таких течений в большой степени отличаются от аналогичных процессов ламинарных течений. Это объясняется тем, что турбулентное движение газа характеризуется возникновением нерегулярных пульсаций (флуктуаций) параметров, которые по своей природе хаотичны. Их точное математическое описание в полной степени невозможно, даже в случае применения специальных методов статистической механики.

В случае турбулентного режима течения, как и в случае ламинарного, плотность теплового потока в любом сечении на боковой поверхности КЛА произвольной формы подчиняется закону конвективного теплообмена Ньютона в виде:

Я0,, =

/ л а

С

V Р У0,7,тур

(1в,1 ,1 ) ;

(15)

где, как и ранее, 1ег - энтальпия восстановления набегающего потока в расчётном сечении, а ^г - энтальпия газа при температуре стенки в расчётном сечении.

При этом параметр теплообмена,

С л а

С

V Р у 0,7, тур

рассчитывался с использованием

критериальных соотношений для турбулентного теплообмена, полученных для тел различной формы. Например, в случае турбулентного теплообмена в потоке с продольным градиентом давления (обтекание тела криволинейной формы), расчет параметра теплообмена осуществляется с помощью критериального соотношения вида:

/■ л а

С

V Р у0,г,тур

= 0,0296 • Яе„Д8 • Рг~0'6

х,

эфД

Т

, г

Т

V е,г у

ч 0,4

1+• Г • м2

,0,11

(16)

где коэффициент восстановления Г = ^РГи .

Кроме того, здесь следует заметить, что в соотношении (16) физические параметры и критерий подобия должны рассчитываться для каждого расчётного сечения криволинейной поверхности. Поэтому в качестве характерного размера в критерии введена эффективная длина Хф /4/. Расчёт хэфг проводится с использованием интегрального соотношения в виде:

Хэфд. =

х 5

\ Я4 ( х) • Рм/,1 • Щ ( х )• ^

_0_

5

4 ( х • и1 ( х )

(17)

В отдельных случаях, когда в области передней критической точки скорость потока не велика, а перепады температуры по толщине пограничного слоя значительные (например, обтекание высокотемпературной сверхзвуковой струей притупленных сферой моделей), расчет параметра теплообмена можно проводить с использованием соотношения:

с л а

С

V Р У 0,/,тур

= 0,0375 • Яе^8 • Рг^0,6 .

/ Л0,6

'А 1 •( X )0'6 -Г р У 1; и

(19),

где ко - радиус притупления, критерии, Ке^ —-, их и р№ - скорость и

плотность газа в расчетной точке "х" на поверхности тела, соответственно, безразмерная координата, х - расстояние по образующей от передней критической точки

С

X — X

до расчетной, к — С- показатель адиабаты.

Результаты расчёта теплообмена на непроницаемой поверхности СКА представлены в виде графических зависимостей на рис.5.

Характер изменения плотности тепловых потоков объясняется режимом течения (ламинарным, переходным и турбулентным) на поверхности спускаемого аппарата.

Рисунок 5. Изменение плотности теплового потока ^о в непроницаемую стенку по

поверхности спускаемого аппарата в диапазоне высот траектории от 90 до 50 км.

1 - 90км, 2 - 80км, 3 - 75км, 4 - 70км, 5 - 65км, 6 - 60км, 7 - 57,5км, 8 - 55км, 9 - 50км

3. Теплообмен на проницаемой поверхности.

Естественный, самопроизвольный процесс просачивания газа через поверхность стенки реализуется при интенсивном её нагреве. В этом случае температура поверхности теплозащитного покрытия достигает больших значений и материал тепловой защиты подвергается физико-химической эрозии, термохимическому разрушению. Это разрушение сопровождается, например, процессами деструкции полимеров, окислением, нитрированием, сублемацией других классов ТЗМ. В ходе этих процессов образовавшиеся на поверхности газообразные продукты разрушения ТЗМ поступают (вдуваются) в пограничный слой, внося в него дополнительную массу. В результате пограничный слой увеличивает толщину, а, следовательно, термическое сопротивление, уменьшая, т.е. блокируя тем самым тепловой поток в стенку.

Наличие процесса вдува газа в пограничный слой сильно снижает плотность теплового потока в стенку. Этот процесс вызывает деформацию профилей распределения по толщине пограничного слоя скорости, энтальпии и массовой концентрации компонентов /1/, /2/.

В таком случае, при наличии разрушения поверхностного слоя ТЗМ уравнение баланса тепловых потоков на проницаемой стенке принимает вид:

Чо _ ^ _ Чед = 0. (21)

В этом уравнении: - плотность теплового потока в передней критической точке в непроницаемую стенку (отсутствие вдува). Согласно /1/

%0 =

Г \

а

С

V СР У 0

(-0 _ I*);

(22)

Чъ, - плотность теплового потока в передней критической точке в проницаемую стенку (наличие вдува). Согласно /1/

%

Г \

а

С

V Ср у *

(-0 _ К);

(23)

Чед - плотность теплового потока, блокированного вдувом. Установлено /10/, /11/, что плотность теплового потока дед пропорциональна тепловому потоку, затраченному на

изменение энтальпии массового расхода вдуваемого в пограничный слой газа от энтальпии на проницаемой поверхности до энтальпии торможения ¡о, т.е.

Явд — У--(1 о - Iw ) , (24)

где у - коэффициент пропорциональности, параметр вдува газа через проницаемую стенку. Согласно /2/

у — 0,6

Г — Л Ме

0,24

V Мохл.

С т \

Ш ,охл

I

-0,03

V 4 у

(25).

Преобразуем (21) с учетом (22), (23) и (24), получим: / \ / \

а

С V Ср У0

•( 1о — К)

а

с

V р Уш

Сократив и преобразовав, окончательно получим:

( \

а С

V СР У № , „ тw — 1— у

( \

а С

V СР У 0

( \

а С

V СР У 0

(26)

Далее введём понятие безразмерной скорости уноса массы ^ в вид

(27).

тт т

С

р У 0

В таком случае (26) принимает вид:

С

р Уш

С

— 1 — у •

(28)

р У 0

Выражение (28) , используя (22) и (23), просто преобразуется к виду:

— — 1 -у-^

Чо

(29)

Таким образом, показано, что для любой точки проницаемой поверхности КЛА между и qw имеет место функциональная зависимость qw — / (Я о ) при одинаковых параметрах набегающего потока.

В ряде работ установлено, что эта зависимость стадийная. Она определяется

интенсивностью процесса вдува газа в пограничный слой, т.е. величиной у • ^^. Для

широкого диапазона этой величины может быть представлена выражениями стадийной аппроксимации:

- умеренный, «линейный» вдув:

= 1 -У • , при 0<у • Ош <0,4;

- интенсивный, «параболический» вдув: 1,012 - 1,16у 0№ + 0,325(у • 0№ )2

г \ а

V ~С, у

г \ а

(30)

V С у

(31)

при 0,4<у • Ош <1,2;

- экстремальный, «экспоненциальный» вдув

Г \ а

V СР у

= ехр

0,37

/ - Ч 0,7

' М-

V М^охл у

• а,

№

при 1,2<у • Ош <8. (32)

По представленному выше алгоритму проведен расчёт плотности теплового потока = /(.х) по поверхности спускаемого космического аппарата с учётом вдува газа в пограничный слой. В качестве материала теплозащитного покрытия в первом приближении был выбран полимерный ТЗМ - фторолон. Плотность полимера р = 1800 кг/м , мольная масса продуктов разложения бралась равной /2/ ¡лохл = 44.

Результаты расчётов представлены на рисунке 6.

Рисунок 6. Изменение плотности теплового потока в проницаемую стенку по поверхности спускаемого аппарата в диапазоне высот траектории от 90 до 40 км.

1- 90км, 2 - 80км, 3 - 75км, 4 - 60км, 5 - 57,5км, 6 - 50км, 7 - 45км, 8 - 42,5км, 9 - 40км

Сравнение результатов влияния блокировки теплового потока эффектом вдува продуктов разрушения фторолона представлены на рис.7.

1800 1600 1400 1200 1000 300 600 100 200

1

А

\ \

я- —в \

1 / Г 2 \

\ 1 \

/ а.

V/

0.0

02

04

06

08

10

12

X. М

Рисунок 7. Сравнение плотности тепловых потоков в непроницаемую и проницаемую стенки на высоте 50 км.

1 - плотность теплового потока в непроницаемую стенку,

2 - плотность теплового потока в проницаемую стенку q0 — плотность теплового потока в непроницаемую стенку (при отсутствии вдува), qw — плотность теплового потока в проницаемую стенку (при наличии вдува продуктов деструкции фторолона).

4. Разработка модели КЛА для экспериментального исследования теплообмена на высокотемпературных газодинамических стендах.

Для экспериментального исследования распределения плотности тепловых потоков по непроницаемой поверхности КЛА малой формы разработана его экспериментальная модель в масштабе М 1:10. В некоторых расчётных сечениях на поверхности модели установлены неохлаждаемые калориметры регулярного режима таблеточного типа. Схема модели представлена на рисунке 8.

131

Рисунок 8: Схема экспериментальной модели спускаемого

КА «Экспресс».

1, 2, 3, 4 - калориметры, Т - клеммы термопар калориметров.

5. Определение толщины уноса фторолона по поверхности КЛА малой формы за время полёта по траектории спуска.

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

V . М Ы "t- ♦

►

г

0 2 4 6 Ч0-Ю-\кВт/м2

Рисунок 9. Зависимость скорости уноса массы фторолона от плотности теплового потока.

Точки - эксперимент

Рисунок 10. Зависимость эффективной энтальпии фторолона от энтальпии

торможения набегающего потока.

Точки - эксперимент

По результатам расчёта на рисунке 11 приведён график изменения толщины уноса фторолона по поверхности КЛА

мм

\ /

V/ \

•-.-- —-—♦

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 X. М

Рисунок 11. Изменение унесённой толщины теплозащитного покрытия из фторолона

по поверхности спускаемого КЛА.

Видно, что толщина уноса покрытия из фторолона достаточно велика. Следовательно, данный теплозащитный материал нежелательно использовать на рассматриваемом аппарате.

Выводы.

1. Составлен алгоритм и проведён расчёт теплообмена на непроницаемой и проницаемой поверхности спускаемого КЛА малой формы по траектории баллистического спуска аппарата.

2. Спроектирована модель КЛА малой формы для проведения экспериментальных исследований на высокотемпературном газодинамическом стенде.

3. По результатам расчёта теплообмена по поверхности КЛА при его баллистическом спуске в атмосфере Земли составлен алгоритм и проведён расчет толщины уноса фторолона, выполняющего функцию разрушаемой тепловой защиты. Показано, что для рассчитанной баллистической траектории КЛА толщина разрушаемого покрытия из фторолона, выбранного в качестве ТЗП КЛА в первом приближении, превышает 15мм. В этой связи требуется замена фторолона на другой теплозащитный материал с повышенными теплозащитными характеристиками. В качестве такого ТЗМ рекомендован композиционный материал стеклопластик.

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

1. Никитин П.В. Тепловая защита. Учебник высшей школы. - М.: Изд-во МАИ, 2006, 510с.

2. Полежаев Ю.В., Юревич Ф.Б. Тепловая защита. - М.: Энергия, 1976, 390с.

3. Лойтянский Л.Г. «Механика жидкости и газа». Изд. «Наука». М.: 1973, 847с.

4. Авдуевский В.С. Галицейский Б.М., Данилов Ю.И. и др. Основы теплопередачи в авиационной и ракетно-космической теплотехнике/ Под ред. АвдуевскогоВ.С., Кошкина В.К. - М.:Машиностроение,1991, 520с.

5. Агофонов В.П., Вертушкин А.А. Неравновесные физико-химические процессы в аэротермодинамике. - М.: Машиностроение, 1972.

6. Ковалев В. Л. Гетерогенные каталитические процессы в аэротермодинамике. — М.: ФИЗМАТЛИТ, 2002. — 224 с. — ISBN 5 - 9221 - 0198 - 6.

7. Каплун А.Б., Морозов Е.М., Олферьева М.А.ANSYS в руках инженера: Практическое руководство. Изд.2, стереот. 2004. 272 с

8. Гиршфельдер Дж., Кертис Ч., Берд Р. «Молекулярная теория газов и жидкостей». М.: ИЛ, 1961, 929с

9. Варгафтик Н.Б., Справочник по теплофизическим свойствам газов и жидкостей. - М.: Физматгиз, 1963.

10. Анфимов Н.А. «Тепло - и массообмен в окрестности критической точки при вдуве и отсосе различных разов через поверхность тела». - «Механика жидкости и газа». 1966, №1.

11. Анфимов Н.А., Альтов В.В. «Теплообмен, трение и массообмен в ламинарном многокомпонентном пограничном слое при вдуве инородных газов». - «Теплофизика высоких температур», 1965,№3, с. 409...420.

Сведения об авторах:

Федорченко Екатерина Александровна инженер ГКНПЦ им. М.В. Хруничева, e-mail: lenafed87@,mail. ru.

Никитин Пётр Васильевич, профессор Московского авиационного института (национального исследовательского университета), д.т.н., е-mail: petrunecha@gmail.ru