епловые процессы при старте

Условия теплового воздействия при старте характеризуются:

· высокими значениями чисел ;

· большим диапазоном изменения температуры в пограничном слое , при высоких значениях до ;

· высоким уровнем пульсаций давления в потоке вследствие действия акустического излучения струй двигательной установки ().

Последнее существенно отличает условия теплообмена при старте от условий теплообмена при обтекании тел невозмущенным потоком.

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

.

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

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

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

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

Так, вместо известной закономерности

,

при увеличении масштаба может наблюдаться рост коэффициента теплообмена с ростом масштаба вследствие уменьшения частоты при увеличении масштаба . Наоборот, с уменьшением масштаба интенсификация теплообмена уменьшается и полностью исчезает при некотором критическом значении и соответствующем значении масштаба . Этим, в частности, объясняется то обстоятельство, что явление усиления теплообмена было обнаружено впервые на крупных натурных объектах (пусковых шахтах).

Для учета влияния пульсаций на теплообмен были разработаны зависимости для двух видов условий по акустическим процессам, которые имели место при старте:

· при реализации автоколебательных резонансных процессов;

· при акустических процессах с шумовым спектром.

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

,

где ;

,

; ; . При .

Для воздуха , для водяного пара .

Для условий воздействия акустических колебаний с шумовым спектром в диапазоне частот усиление теплообмена может быть определено по формуле

,

где – спектральная плотность в полосе .

В этой формуле пороговое значение , при котором начинается влияние акустических колебаний, составляет .

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

,

,

,

где – толщина акустического погранслоя.

Этот комплекс можно представить через , и число потока:

.

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

1. Тепловой поток к пусковому контейнеру (слабо градиентное течение):

; ; ,

.

2. Обтекание отражателя (градиентное течение):

,

где или .

Иногда применяют формулу

.

В расчете берутся , при ,

; ,

где – угол наклона отражателя.

Параметры струи изменяются скачкообразно на длине струи :

· в зоне первой газодинамической структуры струи ;

· на длине шести структур струи .

За время старта .

3. При натекании струи на преграду коэффициент теплоотдачи зависит от градиента скорости растекания потока

.

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

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

3.3.4. Особенности конвективного теплообмена при воздействии двухфазного потока

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

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

,

где – теплоемкость.

Плотность теплового потока может быть определена по следующей формуле:

.

Следует отметить, что значение разности энтальпий газа будет меньше, чем без частиц, вследствие отбора тепла от газа на нагрев частиц.

Оценки показывают, что величина теплового потока с твердыми частицами может быть выше, чем в газовом потоке, на .

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

При воздействии потока с К-фазой, помимо механизма турбулентного потока и усиления теплопередачи конвекцией, существенное влияние оказывает процесс конденсации жидких частиц на более холодной стенке. Этот новый, характерный для условий старта, эффект был исследован под руководством Ю.А. Акопяна. В этом случае тепловой поток может возрастать в несколько () раз. Конденсация частиц на более холодной, чем поток, стенке, приводит к образованию твердой пленки, которая играет роль теплозащиты для металлической стенки, снижая тепловой поток за счет роста и соответствующего уменьшения температурного или энтальпийного напора или . Характер изменения теплового потока представлен на рис.3.32.

 

Рис.3.32. Сравнение тепловых потоков от действия
струй с К-фазой и без К-фазы

При этом интегральное количество тепла в процессе старта может быть одинаково в потоке с К-фазой или без К-фазы.

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

 

3.3.5. Лучистый теплообмен от газов двигательной установки при старте

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

Газовые струи продуктов сгорания топлив ракетных двигателей излучают тепловую энергию, в основном, за счет присутствия в них водяного пара , углекислого газа и других трехатомных компонент. Излучательная способность двухатомных газов , и ничтожна, то есть эти газы для тепловых лучей диатермичны. Газы излучают и поглощают тепловую энергию лишь в определенных интервалах длин волн – в полосах . Тепловая энергия в газе, в отличие от твердых тел, испускается и поглощается в объеме пропорционально длине прохождения луча и давлению .

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

Общая лучеиспускательная способность газа равна сумме энергии лучеиспускания на всех полосах частот:

,

где .

В основу практических расчетов лучеиспускания газа к стенкам может быть положен закон Стефана – Больцмана:

,

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

– берется по ; .

Можно принять для приближенных расчетов .

Если , то стенка не поглощает, а излучает энергию газу.

Закон Стефана – Больцмана справедлив для теплообмена между лучистой полусферой и плоскостью. Излучение объемов газов другого вида можно свести к излучению полусферы через эквивалентный радиус:

,

где – объем газового тела; – площадь оболочки.

Суммарный тепловой поток будет:

,

где .

Для лучистого теплообмена в каналах шахт можно приближенно считать

;

;

;

, .

Для плоских потоков .

 

3.3.6. Нагрев стенок конструкций пусковой установки

Нагрев конструкций газоотражателей с большой толщиной стенки

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

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

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

,

где – коэффициент теплопроводности стенки; – нормаль к поверхности.

Тогда .

Уравнение теплопроводности в общем виде:

,

где – температуропроводность стенки.

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

.

Напишем уравнение теплопроводности через конечные разности; разбив стенку толщиной на слоев толщиной и рассматривая температуру середины слоев в -е моменты времени с интервалом :

;

;

;

.

Для применения графического решения примем связь между и выбранным из условия: . Одновременно, равенство единице комплекса, являющегося числом Фурье , означает возможность достижения тепловым импульсом конца слоя (). При этом условии уравнение теплопроводности в конечных разностях запишется как равенство температуры -го слоя в момент времени средней температуре окружающих слоев и в предшествующий момент времени :

.

Таким образом, система трех уравнений включает:

· условие равенства тепловых потоков на поверхности пластины от газа и внутрь пластины при

; (3.1)

· упрощенное уравнение для температуры слоев пластины

; (3.2)

· условие связи и

. (3.3)

В системе уравнений принято (рис.3.33):

· – номера выделенных слоев толщиной и их средних сечений;

· – номер выделенного фиктивного сечения на расстоянии от поверхности стенки ;

· – относительное время , где из уравнения (3.3).

Рис.3.33. Графическое решение нестационарного нагрева

 

Воспроизводим графически условие уравнения (3.1) в момент времени : равенство градиента температуры у стенки в нулевом фиктивном слое (при температуре ) градиенту температуры в слое толщиной .

Здесь в момент времени 0 в графическом построении допускается некоторое отступление от структуры формулы (3.1) в том, что вместо градиента берется градиент , что обеспечивает выполнение условия и условия (3.1).

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

Так, для следующего момента времени будет для :

,

для слоя .

Согласно уравнению (3.1), между температурой внешнего слоя и температурой газа будет одинаковый градиент, что определит распределение температур в момент времени .

Для последующего момента времени будет:

· для слоя ;

· ;

· .

По температурам и , согласно уравнению (3.1), определится температура внешнего фиктивного слоя () и температура на поверхности пластины.

В последующих построениях для моментов времени прогрев достигает последовательно слоев (рис.3.33).

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

При изменении условий нагрева по и точка для расчета нагрева в последующие времена соответственно сдвигается и проводится графическое построение согласно вышеприведенному.

Пример 3.4.

Подготовка данных для графического решения.

Задано: , , сталь ,

,

, .

Определяем:

, .

Если , , .

 

Нагрев тонкостенных конструкций

Тонкостенные конструкции из теплопроводных материалов (алюминий – , сталь – ) при высокой интенсивности теплопередачи , когда критерий Био , не имеют существенного перепада температур по толщине и им можно пренебречь при расчете нагрева конструкции по времени.

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

;

,

где – постоянная Стефана–Больцмана;

,

, – степень черноты поверхностей.

 

Рис.3.34. Схема теплового воздействия на стенки
транспортно-пускового контейнера

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

Тепловая нагрузка на отражатель определяется интегралом теплового потока по времени. Для упрощения оценки тепловой нагрузки на отражатель заменим определенную в натурных условиях кривую ступенчатым изменением , чтобы были одинаковыми интегральные тепловые нагрузки (рис.3.35):

.

 

Рис.3.35. Оценка уноса стали от воздействия теплового потока

 

Определим температуру стенки при заданном тепловом потоке по приближенной формуле прогрева полубесконечного тела:

, при ,

.

Определим время начала плавления стального газоотражателя при :

.

Для стали

, , .

Тогда при и .

При достижении начинается унос металла. Количество тепла, расходуемого на унос металла, можно принять приближенно как избыточное интегральное количество тепла над кривой от начала плавления до снижения теплового потока ниже значения в момент .

Пример 3.5.

Определим координаты кривой для стальной пластины (см. таблицу 3.4).

 

Таблица 3.4

,
, с	0,4	0,481	0,6	0,73	0,92	1,18	1,65	2,31

 

Тогда можно определить избыточную теплоту, уходящую на плавление материала:

,

, .

Унос материала отражателя можно определить через теплоту плавления по формуле

,

где , .

В рассматриваемом примере ,

.

Работа некоторых материалов при одинаковых приведена в таблице 3.5 и на рис.3.36.

 

Рис.3.36. Сравнение температур поверхности нагрева различных
материалов при

 

Таблица 3.5

Материал			1Х18Н9Т	Бетон
Начало плавления, с	0,52	0,9	0,8	0,1
Длительность плавления, с	0,68	0,36	0,68	1,4
Унос, мм		~0,1	0,3
Данные натурных измерений, мм		0,1

Из рис.3.36 и данных таблицы 3.5 следует, что целесообразным является применение более теплопроводных материалов.

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

Сравним работу материалов ( -19, , , бетон) для теплового потока .

Применение теплозащитных покрытий

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

,

где – допустимая температура металлической конструкции; – эффективная энтальпия покрытия ; – температура разрушения покрытия ( – легкоразрушимое покрытие).

Для покрытия ПКМ-6: