Теплофизика высоких температур, 2022, T. 60, № 6, стр. 880-887

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

В. В. Миронов 1*, М. А. Толкач 1, В. В. Тлевцежев 1

1 ГНЦ ФГУП “Центр Келдыша”
Москва, Россия

* E-mail: kerc@elnet.msk.ru

Поступила в редакцию 10.12.2021
После доработки 09.04.2022
Принята к публикации 07.06.2022

Полный текст (PDF)

Аннотация

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

ВВЕДЕНИЕ

В качестве теплозащитных покрытий (ТЗП) в высокотемпературных энергоустановках, в том числе в ракетных двигателях твердого топлива [13], широко применяются резиноподобные теплозащитные материалы (ТЗМ). В процессе их работы под воздействием радиационного qл и конвективного qк тепловых потоков материалы ТЗП интенсивно прогреваются. Температура поверхности быстро возрастает и при достижении значений 500–600 К начинается процесс термического разложения исходного ТЗМ (рис. 1). Здесь ГПД – газообразные продукты деструкции ТЗП; СОК – силовая оболочка корпуса; Ci, Cv – концентрации активных компонентов соответственно из нагревающего потока (окислителя) и вдуваемого пиролизного газа (горючего); T, u – температура и скорость газового потока в пограничном слое; ${{\dot {m}}_{{{\text{пг}}}}}$ – массовый поток вдуваемых продуктов деструкции ТЗМ (пиролизных газов) с единицы поверхности кокса.

Рис. 1.

Схема прогрева ТЗП и диффузионного окисления пиролизных газов в пограничном слое.

Зона разложения с переменной скоростью перемещается внутрь покрытия, а на поверхности образуется пористый коксовый остаток. Термическое разложение материалов сопровождается выделением газовых компонентов, таких как метан, этилен, водород, ацетилен, бензол, которые фильтруются через коксовый слой, вдуваются в газовый поток и при больших температурах также взаимодействуют с окисляющими компонентами в пограничном слое, главным образом с парами воды и углекислым газом из нагревающего потока (в частности, содержащимися в продуктах сгорания твердых топлив в реактивных двигателях [1, 4]). Неизрасходованные в пограничном слое химически активные компоненты продуктов нагревающего потока могут достигать поверхности прококсованного слоя и вступать в реакции с углеродом кокса. В итоге коксовый слой разрушается на поверхности.

Большинство перечисленных процессов являются нестационарными и реализуются с поглощением тепла как внутри материала, так и в пограничном слое, а также на поверхности коксового слоя при его омывании нагревающим потоком.

Задача о расчете турбулентного пограничного слоя при наличии гомогенного реагирования в пограничном слое и гетерогенного реагирования на поверхности ТЗП представляет значительный интерес [1, 5]. Данная задача возникает при определении скорости химического разрушения коксового остатка ТЗП при вдуве пиролизных газов, содержащих реакционно способные вещества.

Условия работы резиноподобных ТЗМ во многих высокотемпературных установках характеризуются относительно низкими значениями коэффициентов тепломассообмена α/Cp ≤ 0.5 кг/(м2 с), относительно высокими давлениями Р ≥ 5.0 МПа, температурами нагревающего потока в пограничном слое до 3500 К и с преобладающим тепловым потоком к поверхности ТЗП за счет излучения [1]. Такие условия дают основание предполагать, что окислительные компоненты из нагревающего потока в значительной мере израсходуются при взаимодействии с вдуваемыми активными компонентами пиролизного газа, выделяемыми, как это подтверждено экспериментально [2], при термодеструкции ТЗП. В силу этого проникновение окислительных компонентов из нагревающего потока к поверхности коксового остатка, а следовательно, и разрушение прококсованного слоя на поверхности уменьшаются.

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

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

Согласно схеме, реальное распределение концентраций реагирующих компонентов заменяется введением одного или нескольких “фронтов” химических реакций. Положение фронта определяется условием стехиометрического смешения потоков условных горючего и окислителя. В случае если плотность массового потока вдуваемых активных пиролизных газов (горючего) превышает стехиометрическое отношение, не израсходованная в зоне гомогенных реакций их доля диффундирует к границе пограничного слоя. При недостаточном для стехиометрического реагирования количестве вдуваемых газов неизрасходованная в зоне гомогенных реакций часть окисляющих компонентов из нагревающего потока диффундирует к поверхности прококсованного слоя ТЗП и может вступать в гетерогенные реакции с углеродом.

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

Оценка влияния состава пиролизных газов на тепловые эффекты в пограничном слое проведена для двух основных составов пиролизного газа (продуктов деструкции ТЗП) [2]:

1) пиролизный газ состоит только из метана CH4;

2) пиролизный газ в виде смеси метана CH4 и этилена C2H4 с известными массовыми концентрациями.

В качестве активных окислительных компонентов нагревающего потока рассматриваются газообразные Н2О и СО2. Считается, что определяющими в уносе углеродных материалов являются химические реакции углерода с данными кислородосодержащими компонентами [1]. Такие компоненты присутствуют, например, в продуктах сгорания твердых ракетных топлив [1]. Концентрация активного кислорода в продуктах сгорания твердых топлив, незначительна [4].

С использованием предлагаемой методики соответствующие оценки могут быть проведены для любого окисляющего компонента. В книге [1] отмечается, что “…неоднозначно трактуется роль хлора, хлористого водорода (HCl) и других соединений хлора в термодинамических процессах с углеродом, характерных для высокотемпературных сред сложного химического состава. По некоторым свидетельствам соединения хлора пассивны во всем температурном диапазоне, вплоть до 3300–3500 К”. В других работах отмечается, что HCl, не принимая непосредственного участия в процессах химического взаимодействия (с углеродом), оказывает ингибирующее воздействие на другие реакции, в том числе и окислительные (с кислородсодержащими веществами). Результаты исследований показывают, что в температурном диапазоне до 2500 К ингибирующее влияние Cl и HCl на 5–15% уменьшает скорость реакций углеродных компонентов.

ОСНОВНЫЕ УРАВНЕНИЯ

Рассматривается решение частной задачи реакции окислительных компонентов из нагревающего потока Н2О и СО2 со смесью активных пиролизных компонентов в виде газа СnHm в пограничном слое, т.е. учитываются следующие реакции:

(1)
$\left. \begin{gathered} {{{\text{C}}}_{n}}{{{\text{H}}}_{m}} + n{{{\text{H}}}_{2}}{\text{O}} = n{\text{CO}} + \frac{m}{2}{{{\text{H}}}_{2}} + n{{{\text{H}}}_{2}}, \hfill \\ {{{\text{C}}}_{n}}{{{\text{H}}}_{m}} + n{\text{С}}{{{\text{O}}}_{2}} = 2n{\text{CO}} + \frac{m}{2}{{{\text{H}}}_{2}}. \hfill \\ \end{gathered} \right\}$

Стехиометрические коэффициенты $~{{\vartheta }_{{vi}}}$ в этих реакциях i-го окисляющего компонента с химически активной долей вдуваемого газа определяются следующим образом:

(2)
$\left. \begin{gathered} {{\vartheta }_{{v,1}}} = {{\vartheta }_{{v,{{{\text{H}}}_{2}}{\text{O}}}}} = n\frac{{{{M}_{{{{{\text{H}}}_{2}}{\text{O}}}}}}}{{{{M}_{{{{{\text{C}}}_{n}}{{{\text{H}}}_{m}}}}}}}, \hfill \\ {{\vartheta }_{{v,2}}} = {{\vartheta }_{{v,{\text{C}}{{{\text{O}}}_{2}}}}} = n\frac{{{{M}_{{{\text{C}}{{{\text{O}}}_{2}}}}}}}{{{{M}_{{{{{\text{C}}}_{n}}{{{\text{H}}}_{m}}}}}}}, \hfill \\ \end{gathered} \right\}$
где М – молекулярная масса соответствующего компонента.

Непрореагировавшие в пограничном слое окислительные компоненты из нагревающего потока достигают поверхности ТЗП и вступают с углеродом прококсованного слоя в следующие гетерогенные реакции:

(3)
$\left. \begin{gathered} {\text{C}} + {{{\text{H}}}_{2}}{\text{O}} = {\text{CO}} + {{{\text{H}}}_{2}}, \hfill \\ {\text{C}} + {\text{C}}{{{\text{O}}}_{2}} = 2{\text{CO}}{\text{.}} \hfill \\ \end{gathered} \right\}$

Эти реакции протекают в тонком мономолекулярном слое, удерживаемом силами химической адсорбции [4]. Для реакций (3) стехиометрические коэффициенты определяются соотношениями

(4)
$\left. \begin{gathered} {{\vartheta }_{{c,1}}} = {{\vartheta }_{{c,{{{\text{H}}}_{2}}{\text{O}}}}} = \frac{{{{M}_{{{{{\text{H}}}_{2}}{\text{O}}}}}}}{{{{M}_{{\text{C}}}}}}, \hfill \\ {{\vartheta }_{{c,2}}} = {{\vartheta }_{{c,{\text{C}}{{{\text{O}}}_{2}}}}} = \frac{{{{M}_{{{\text{C}}{{{\text{O}}}_{2}}}}}}}{{{{M}_{{\text{C}}}}}}. \hfill \\ \end{gathered} \right\}$

Как в гомогенных c пиролизным газом (1), так и гетерогенных с коксовым остатком (3) реакциях образуются два компонента: CO и H2. Далее полагается, что эти компоненты химически инертны.

Легко проверить, что отношение стехиометрических коэффициентов $\frac{{{{\vartheta }_{{v,i}}}}}{{{{\vartheta }_{{c,i}}}}}$ для всех i окисляющих компонентов является константой, равной $n\frac{{{{M}_{{\text{C}}}}}}{{{{M}_{{{{{\text{C}}}_{n}}{{{\text{H}}}_{m}}}}}}}$.

Строгое определение одной из важных переменных гомогенного реагирования в пограничном слое – плотности массового потока вдуваемых с единицы поверхности ТЗП (кокса) активных пиролизных газов ${{\dot {m}}_{{{\text{пг}}}}}$ возможно при решении сопряженной задачи тепломассообмена в пограничном слое совместно с решением стефановской задачи теплопроводности с соответствующими граничными условиями [1]. Для оценок используется приближенное соотношение [1], подтвержденное в том числе результатами экспериментальных исследований [3]:

(5)
${{\dot {m}}_{{{\text{пг}}}}} = \left( {1 - {\text{К}}} \right){{\rho }_{{{\text{тзп}}}}}\frac{\xi }{{2{{\tau }^{{1/2}}}}}.$
Здесь К – коксовое число ТЗМ (К = ρк/${{\rho }_{{{\text{тзп}}}}}$); ρк – плотность прококсованного слоя ТЗМ (кг/м3); ρтзп – плотность исходного ТЗМ (кг/м3); ξ – экспериментальный коэффициент, определяющий глубину и скорость термической деструкции (м/с1/2); τ – время воздействия нагревающего потока на ТЗП.

Для оценки ${{\dot {m}}_{{{\text{пг}}}}}$ взяты следующие значения: К = 0.3, ρтзп = 1100 кг/м3, ξ = 0.7 × 10–3 м/с1/2 [1, 3]. В этом случае, согласно (5), получаем ${{\dot {m}}_{{{\text{пг}}}}}$ = 0.270/τ1/2 (кг/(м2 с)). Соответствующая зависимость плотности массового потока вдуваемых в пограничный слой пиролизных газов от времени прогрева теплозащитного материала показана на рис. 2. Достаточно интенсивный начальный вдув пиролизных газов в первые 10–15 с прогрева ТЗП резко ослабевает.

Рис. 2.

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

Для определения доли окисляющих компонентов, расходуемых во фронте гомогенных реакций в пограничном слое (рис. 1), воспользуемся соотношением, определяющим стехиометрическое количество (плотность массового потока) ${{\dot {m}}_{{{\text{пг,}}~{\text{st}}}}}$ пиролизного газа, необходимое для осуществления гомогенных реакций с расходованием всех окисляющих компонентов, которые поступают из нагревающего потока в зону гомогенных реакций. В соответствии с законом Фика в приближении аналогии между процессами тепло- и массообмена для диффузионной модели реагирования в пограничном слое можно написать [1, 4, 6]

(6)
${{\dot {m}}_{{{\text{пг}},{\text{st}}}}} = \frac{\alpha }{{{{с}_{p}}}}{{B}_{{mv}}}{{\left( {{\text{L}}{{{\text{e}}}_{i}}} \right)}^{{0.7}}},$
(7)
${{B}_{{mv}}} = \mathop \sum \limits_i \frac{{{{C}_{{i,o}}}}}{{{{\vartheta }_{{{{v}_{{{\text{см}}}}},i}}}}}.$
Здесь α/ср – обобщенный коэффициент тепломассообмена; Bmv – окислительный потенциал по отношении к активному пиролизному газу; ${{\vartheta }_{{{{v}_{{{\text{см}}}}},i}}}$ – стехиометрические коэффициенты для реакций пиролизного газа с i-й окисляющей компонентой из нагревающего потока (2); ${{C}_{{i,o}}}$ – массовая концентрация i-й окисляющей компоненты на границе пограничного слоя; Lei – число Льюиса (переменная поправка к аналогии между коэффициентами тепло- и массообмена [6]). Для многокомпонентной газовой смеси определение значения числа Льюиса является самостоятельной сложной задачей. Как принято в инженерной практике [4], полагается, что ${\text{L}}{{{\text{e}}}_{i}} = 1$.

При вдуве в пограничный слой активных пиролизных газов с плотностью массового потока ${{\dot {m}}_{{{\text{пг}}}}}$ (5), равным или больше ${{\dot {m}}_{{{\text{пг}},{\text{st}}}}}$ (6), все окисляющие компоненты из нагревающего потока израсходуются во фронте (рис. 1), а гетерогенные реакции с углеродом прококсованного слоя и соответствующие тепловые эффекты будут исключены.

Если ${{\dot {m}}_{{{\text{пг}}}}}$ < ${{\dot {m}}_{{{\text{пг}},{\text{st}}}}}$, то поверхности прококсованного слоя будут достигать только не израсходованные (на реакции с активным пиролизным газом) в пограничном слое окисляющие компоненты. При Lei = 1 и в предположении пропорционального ${{\dot {m}}_{{{\text{пг}}}}}$ расходования окисляющих компонентов получаем следующее соотношение для расчета расхода массы окисленного углерода ${{\dot {m}}_{{\text{C}}}}$ (кг/(м2 с)) c единицы поверхности прококсованного слоя:

(8)
${{\dot {m}}_{{\text{C}}}} = \frac{\alpha }{{{{с}_{p}}}}P{{B}_{{m{\text{C}}}}},\,\,\,\,{{B}_{{m{\text{C}}}}} = \mathop \sum \limits_i \frac{{{{C}_{{i,o}}}}}{{{{\vartheta }_{{c,i}}}}},$
(9)
$P = 1--({{\dot {m}}_{{{\text{пг}}}}}{\text{/}}{{\dot {m}}_{{{\text{пг}},{\text{st}}}}}),$
где BmC – окислительный потенциал нагревающего потока по отношении к углероду прококсованного слоя ТЗП [6] при вдуве инертного пиролизного газа; P – поправка на изменение BmC из-за уменьшения количества окисляющих компонентов в результате гомогенных реакций в пограничном слое (при вдуве активного пиролизного газа); ${{\vartheta }_{{с,i}}}$ – стехиометрические коэффициенты для реакций углерода прококсованного слоя с i-й окисляющей компонентой из нагревающего потока (4). Если ${{\dot {m}}_{{{\text{пг}}}}}$ > ${{\dot {m}}_{{{\text{пг,st}}}}}$ то Р = 0.

Используя (5) и (6), определяем время прогрева ТЗП τst, при котором ${{\dot {m}}_{{{\text{пг,st}}}}} = {{\dot {m}}_{{{\text{пг}}}}}$:

${{\tau }_{{{\text{st}}}}} = {{\left( {\frac{{\xi \left( {1 - {\text{К}}} \right){{\rho }_{{{\text{тзп}}}}}}}{{2\frac{\alpha }{{{{с}_{p}}}}{{B}_{{mv}}}}}} \right)}^{2}}.$

Начиная с этого времени, происходит проникновение окисляющих компонентов к поверхности деструктирующегося покрытия. Зависимости τst от параметров ТЗП и условий теплообмена показана на рис. 3, 4. Расчеты проведены при следующих опорных значениях: К = 0.3, ρтзп = 1100 кг/м3, ξ = 0.7 × 10–3 м/с1/2, α/ср = 0.7 кг/(м2 с).

Рис. 3.

Зависимости τst от обобщенного коэффициента тепломассообмена при К = 0.3, ξ = 0.7 × 10–3 м/с1/2 и различных значениях Bmv: 1 – 0.0833, 2 – 0.0904, 3 – 0.0952.

Рис. 4.

Зависимости τst от коксового числа ТЗП при α/ср = 0.7 кг/(м2 с), ξ = 0.7 × 10–3 м/с1/2 и различных значениях Bmv: 1 – 0.0833, 2 – 0.0904, 3 – 0.0952.

Диапазон изменения коэффициентов тепломассообмена α/ср выбран, согласно данным [1, 6], для условий прогрева в высокотемпературных установках. Диапазон изменения коксового числа К принят для коксующихся резиноподобных материалов (0.3 – для эластичных и 0.8 – для армированных ТЗМ). Полученные результаты указывают на то, что время начала диффузионного режима разрушения прококсованного слоя существенным образом зависит от обобщенного коэффициента тепломассообмена α/ср и коксового числа теплозащитного покрытия К. Для наиболее распространенных эластичных ТЗМ с коксовым числом К ≈ 0.3 при условиях тепломассообмена с α/ср ≈ ≈ 0.7 кг/(м2 с) задержка в начале окисления прококсованного слоя составляет ~20 с.

Согласно данным рис. 4, чем меньше коксовое число, тем более сильное влияние на время начала уноса прококсованного слоя оказывает величина окислительного потенциала гомогенных реакций Bmv. Из результатов расчетов (рис. 3, 4) также следует, что увеличение окислительного потенциала гомогенных реакций Bmv на ~15% приводит к уменьшению времени τst начала гетерогенных реакций на поверхности прококсованного слоя на ~25–30%, т.е. к более раннему наступлению разрушения прококсованного слоя.

Рассмотрим изменение поправки P (9), определяющей снижение окислительного потенциала активных компонентов из нагревающего потока в отношении углерода прококсованного слоя BmC, обусловленное расходованием окисляющих компонентов на гомогенные реакции в пограничном слое. На рис. 5 приведена зависимость P от времени прогрева ТЗП. Расчеты выполнены при следующих значениях: К = 0.3, ρтзп = 1100 кг/м3, ξ = 0.7 × 10–3 м/с1/2.

Рис. 5.

Зависимости поправки к окислительному потенциалу относительно углерода кокса от времени воздействия нагревающего потока на ТЗП при различных условиях массообмена: 1 – α/ср = 0.3 кг/(м2 с), 2 – 0.5, 3 – 0.7.

Окислительный потенциал по отношению к пиролизным газам в виде смеси ${\text{C}}{{{\text{H}}}_{4}}$ и ${{{\text{C}}}_{2}}{{{\text{H}}}_{4}}$ принят равным Bmv = 0.0904 (оценки величин окислительных потенциалов приведены ниже).

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

Обратимся теперь к тепловым эффектам рассмотренных выше гомогенных (1) и гетерогенных (3) реакций. Теплота любой реакции (даже пока не осуществленной на практике) может быть вычислена как разность между суммой полных энтальпий всех реагентов и суммой полных энтальпий продуктов реакции [8]. Тогда для вычисления теплового эффекта с пиролизными активными газами и углеродом коксового остатка используем соотношения, следующие из уравнений (1) и (3):

а) Н2О и СО2 с пиролизными активными газами (1)

(10)
$\left. \begin{gathered} Q_{{{{{\text{H}}}_{2}}{\text{O}},w}}^{v} = {{H}_{{v,w}}} + {{\vartheta }_{{v,1}}}{{H}_{{{{{\text{H}}}_{2}}{\text{O}},w}}} - \hfill \\ - \,\,n\frac{{{{M}_{{{\text{CO}}}}}}}{{{{M}_{v}}}}{{H}_{{{\text{CO}},w}}} - \left( {n + \frac{m}{2}} \right)\frac{{{{M}_{{{{{\text{H}}}_{2}}}}}}}{{{{M}_{v}}}}{{H}_{{{{{\text{H}}}_{2}},w}}}, \hfill \\ Q_{{{\text{C}}{{{\text{O}}}_{2}},w}}^{v} = {{H}_{{v,w}}} + {{\vartheta }_{{v,2}}}{{H}_{{{\text{C}}{{{\text{O}}}_{2}},w}}} - \hfill \\ - \,\,2n\frac{{{{M}_{{{\text{CO}}}}}}}{{{{M}_{v}}}}{{H}_{{{\text{CO}},w}}} - \frac{m}{2}\frac{{{{M}_{{{{{\text{H}}}_{2}}}}}}}{{{{M}_{v}}}}{{H}_{{{{{\text{H}}}_{2}},w}}}; \hfill \\ \end{gathered} \right\}$
б) Н2О и СО2 с углеродом коксового остатка (3)
(11)
$\left. \begin{gathered} Q_{{{{{\text{H}}}_{2}}{\text{O}},w}}^{{\text{C}}} = {{H}_{{{\text{C}},w}}} + {{\vartheta }_{{c,1}}}{{H}_{{{{{\text{H}}}_{2}}{\text{O}},w}}} - \hfill \\ - \,\,\frac{{{{M}_{{{\text{CO}}}}}}}{{{{M}_{{\text{C}}}}}}{{H}_{{{\text{CO}},w}}} - \frac{{{{M}_{{{{{\text{H}}}_{2}}}}}}}{{{{M}_{{\text{C}}}}}}{{H}_{{{{{\text{H}}}_{2}},w}}}, \hfill \\ Q_{{{\text{C}}{{{\text{O}}}_{2}},w}}^{{\text{C}}} = {{H}_{{{\text{C}},w}}} + {{\vartheta }_{{c,2}}}{{H}_{{{\text{C}}{{{\text{O}}}_{2}},w}}} - 2\frac{{{{M}_{{{\text{CO}}}}}}}{{{{M}_{c}}}}{{H}_{{{\text{CO}},w}}}. \hfill \\ \end{gathered} \right\}$
Здесь $Q_{{i,w}}^{v}$ и $Q_{{i,w}}^{{\text{C}}}$ – тепловые эффекты при реагировании единицы массы вдуваемого газа и единицы массы углерода прококсованного слоя ТЗП с соответствующими i-ми (Н2О и СО2) химически активными компонентами набегающего потока; ${{H}_{{v,w}}}$ – полная энтальпия пиролизного газа при рассматриваемой температуре, Дж/кг; ${{H}_{{{\text{C}},w}}}$ – полная энтальпия углерода прококсованного слоя, Дж/кг; ${{H}_{{i,w}}}$ – полная энтальпия i окислительных компонентов (Н2О и СО2), Дж/кг; ${{H}_{{{\text{CO}},w}}}$, ${{H}_{{{{{\text{H}}}_{2}},~w}}}$ – полная энтальпия продуктов реакций (${\text{CO}}$ и ${{{\text{H}}}_{2}}$) в пограничном слое и на поверхности кокса, Дж/кг; все значения энтальпий взяты при рассматриваемой температуре. Индексы v и С в уравнениях указывают на пиролизный газ и углерод прококсованного слоя.

Как известно, полная энтальпия вещества при температуре T определяется соотношением [9, 10]

(12)
$\widetilde {~H}_{T}^{0} = H_{f}^{0} + \left( {H_{T}^{0} - H_{{{{T}_{0}}}}^{0}} \right),$
где $H_{f}^{0}~$ − стандартный тепловой эффект образования вещества из элементов при T0 и P0; $H_{T}^{0}$, $H_{{{{T}_{0}}}}^{0}$ − энтальпии рассматриваемого вещества при температурах T и T0 соответственно.

Рассмотрим соотношения для смеси CH4 + C2H4. В этом случае смесь определяется условной формулой ${{{\text{С}}}_{{{{n}_{{{\text{см}}}}}}}}{{{\text{Н}}}_{{{{m}_{{{\text{см}}}}}}}}$. В этом случае для реакций смеси с Н2О справедливо следующее аналогичное (10) соотношение:

(13)
$\begin{gathered} Q_{{{{{\text{H}}}_{2}}{\text{O}},w}}^{{{{v}_{{{\text{см}}}}}}} = {{H}_{{{{v}_{{{\text{см}}}}},w}}} + {{\vartheta }_{{{{v}_{{{\text{см}}}}},1}}}{{H}_{{{{{\text{H}}}_{2}}{\text{O}},w}}} - \\ - \,\,\frac{{{{n}_{{{\text{см}}}}}}}{{{{M}_{{{{v}_{{{\text{см}}}}}}}}}}{{M}_{{{\text{CO}}}}}{{H}_{{{\text{CO}},w}}} - \frac{{\left( {{{n}_{{{\text{см}}}}} + \frac{{{{m}_{{{\text{см}}}}}}}{2}} \right)}}{{{{M}_{{{{v}_{{{\text{см}}}}}}}}}}{{M}_{{{{{\text{H}}}_{2}}}}}{{H}_{{{{{\text{H}}}_{2}},w}}}, \\ \end{gathered} $
где ${{H}_{{{{v}_{{{\text{см}}}}},w}}}$ – полная энтальпия смеси пиролизного газа, Дж/кг; ${{M}_{{{{v}_{{{\text{см}}}}}}}}$ – молекулярная масса смеси; ${{\vartheta }_{{{{v}_{{{\text{см}}}}},1}}}$ – стехиометрический коэффициент в реакции смеси пиролизного газа ${{{\text{С}}}_{{{{n}_{{{\text{см}}}}}}}}{{{\text{Н}}}_{{{{m}_{{{\text{см}}}}}}}}$ с Н2О.

Далее полагаем, что для каждой из i окисляющих компонентов справедливо

$Q_{{i,w}}^{{{{v}_{{{\text{см}}}}}}} = {{g}_{{{{v}_{1}}}}}Q_{{i,w}}^{{{{v}_{1}}}} + {{g}_{{{{v}_{2}}}}}Q_{{i,w}}^{{{{v}_{2}}}},$
где ${{g}_{{{{v}_{1}}}}}$, ${{g}_{{{{v}_{2}}}}}$ – массовые доли CH4 и C2H4 в смеси пиролизного газа, т.е. индексы v1 и v2 указывают соответственно на метан и этилен.

Соотношения для вычисления nсм и mсм в условной формуле смеси ${{{\text{С}}}_{{{{n}_{{{\text{см}}}}}}}}{{{\text{Н}}}_{{{{m}_{{{\text{см}}}}}}}}$ могут быть получены, например, из первого уравнения (10). Используя это уравнение для метана и этилена и умножая соответственно на ${{g}_{{{{v}_{1}}}}}$ и ${{g}_{{{{v}_{2}}}}}$, получаем два соотношения:

$\begin{gathered} {{g}_{{{{v}_{1}}}}}Q_{{{{{\text{H}}}_{2}}{\text{O}},w}}^{{{{v}_{1}}}} = {{g}_{{{{v}_{1}}}}}{{H}_{{{{v}_{1}},w}}} + {{g}_{{{{v}_{1}}}}}{{\vartheta }_{{{{v}_{1}},1}}}{{H}_{{{{{\text{H}}}_{2}}{\text{O}},w}}} - \\ - \,\,{{g}_{{{{v}_{1}}}}}{{n}_{1}}\frac{{{{M}_{{{\text{CO}}}}}}}{{{{M}_{{{{v}_{1}}}}}}}{{H}_{{{\text{CO}},w}}} - {{g}_{{{{v}_{1}}}}}\left( {{{n}_{1}} + \frac{{{{m}_{1}}}}{2}} \right)\frac{{{{M}_{{{{{\text{H}}}_{2}}}}}}}{{{{M}_{{{{v}_{1}}}}}}}{{H}_{{{{{\text{H}}}_{2}},w}}}, \\ {{g}_{{{{v}_{2}}}}}Q_{{{{{\text{H}}}_{2}}{\text{O}},w}}^{{{{v}_{2}}}} = {{g}_{{{{v}_{2}}}}}{{H}_{{{{v}_{2}},w}}} + {{g}_{{{{v}_{2}}}}}{{\vartheta }_{{{{v}_{2}},1}}}{{H}_{{{{{\text{H}}}_{2}}{\text{O}},w}}} - \\ - \,\,{{g}_{{{{v}_{2}}}}}{{n}_{2}}\frac{{{{M}_{{{\text{CO}}}}}}}{{{{M}_{{{{v}_{2}}}}}}}{{H}_{{{\text{CO}},w}}} - {{g}_{{{{v}_{2}}}}}\left( {{{n}_{2}} + \frac{{{{m}_{2}}}}{2}} \right)\frac{{{{M}_{{{{{\text{H}}}_{2}}}}}}}{{{{M}_{{{{v}_{2}}}}}}}{{H}_{{{{{\text{H}}}_{2}},w}}}. \\ \end{gathered} $

Сложение данных выражений дает выражение уравнение для смеси метана с этиленом

(14)
$\begin{gathered} Q_{{{{{\text{H}}}_{2}}{\text{O}},w}}^{{{{v}_{{{\text{см}}}}}}} = \left( {{{g}_{{{{v}_{1}}}}}{{H}_{{{{v}_{1}},w}}} + {{g}_{{{{v}_{2}}}}}{{H}_{{{{v}_{2}},w}}}} \right) + \\ + \,\,\left( {{{g}_{{{{v}_{1}}}}}{{\vartheta }_{{{{v}_{1}},1}}} + {{g}_{{{{v}_{2}}}}}{{\vartheta }_{{{{v}_{2}},1}}}} \right){{H}_{{{{{\text{H}}}_{2}}{\text{O}},w}}}~ - \\ - \,\,\left( {\frac{{{{g}_{{{{v}_{1}}}}}{{n}_{1}}}}{{{{M}_{{{{v}_{1}}}}}}} + \frac{{{{g}_{{{{v}_{2}}}}}{{n}_{2}}}}{{{{M}_{{{{v}_{2}}}}}}}} \right){{M}_{{{\text{CO}}}}}{{H}_{{{\text{CO}},w}}} - \\ - \,\,\left[ {{{g}_{{{{v}_{1}}}}}\frac{{\left( {{{n}_{1}} + \frac{{{{m}_{1}}}}{2}} \right)}}{{{{M}_{{{{v}_{1}}}}}}} + {{g}_{{{{v}_{2}}}}}\frac{{\left( {{{n}_{2}} + \frac{{{{m}_{2}}}}{2}} \right)}}{{{{M}_{{{{v}_{2}}}}}}}} \right]{{M}_{{{{{\text{H}}}_{2}}}}}{{H}_{{{{{\text{H}}}_{2}},w}}}. \\ \end{gathered} $
Здесь n1, n2 – количество атомов углерода в молекулах метана и этилена; m1, m2 – количество атомов водорода в молекулах метана и этилена.

Из сопоставления (13) и (14) следуют соотношения для вычисления полной энтальпии смеси

${{H}_{{{{v}_{{{\text{см}}}}},w}}} = \left( {{{g}_{{{{v}_{1}}}}}{{H}_{{{{v}_{1}},w}}} + {{g}_{{{{v}_{2}}}}}{{H}_{{{{v}_{2}},w}}}} \right)$
и стехиометрического коэффициента для смеси пиролизных газов c i-м окисляющим компонентом (i = 1 для Н2О, i = 2 для СО2)

$\begin{gathered} {{\vartheta }_{{{{v}_{{{\text{см}}}}},i}}} = \left( {{{g}_{{{{v}_{1}}}}}{{\vartheta }_{{{{v}_{1}},i}}} + {{g}_{{{{v}_{2}}}}}{{\vartheta }_{{{{v}_{2}},i}}}} \right), \\ {{\vartheta }_{{{{v}_{j}},i}}} = {{n}_{{{{v}_{j}}}}}\frac{{{{M}_{i}}}}{{{{M}_{{{{v}_{j}}}}}}},\,\,\,\,j = 1,{\text{ }}2. \\ \end{gathered} $

Индекс $j$ = 1 соответствует метану и j = 2 ‒ этилену C2H4.

После несложных преобразований из (13) и (14) можно получить следующие соотношения для вычисления nсм и mсм для двухкомпонентной смеси:

(15)
${{n}_{{{\text{см}}}}} = \left( {\frac{{{{g}_{{{{v}_{1}}}}}}}{{{{M}_{{{{v}_{1}}}}}}}{{n}_{1}} + \frac{{{{g}_{{{{v}_{2}}}}}}}{{{{M}_{{{{v}_{2}}}}}}}{{n}_{2}}} \right){{M}_{{{{v}_{{{\text{см}}}}}~}}},$
(16)
${{m}_{{{\text{см}}}}} = \left( {\frac{{{{g}_{{{{v}_{1}}}}}}}{{{{M}_{{{{v}_{1}}}}}}}{{m}_{1}} + \frac{{{{g}_{{{{v}_{2}}}}}}}{{{{M}_{{{{v}_{2}}}}}}}{{m}_{2}}} \right){{M}_{{{{v}_{{{\text{см}}}}}}}},$
где ${{M}_{{{{v}_{{{\text{см}}}}}}}}$ – условная молекулярная масса смеси.

Аналогичные (15), (16) соотношения для вычисления nсм и mсм используются в [11].

После соответствующих вычислений удельная формула смеси пиролизных газов при ${{M}_{{{{v}_{{{\text{см}}}}}}}}$ = = 1 г/моль может быть представлена в виде ${{{\text{C}}}_{{0.0658}}}{{{\text{H}}}_{{0.2096}}}$. Для ${{M}_{{{{v}_{{{\text{см}}}}}}}}$ = 100 г/моль удельная формула представляется в виде ${{{\text{C}}}_{{6.58}}}{{{\text{H}}}_{{20.96}}}$, т.е. nсм = 6.58, mсм = 20.96.

С учетом концентрации Ci,oактивных компонентов в нагревающем потоке на границе пограничного слоя суммарное количество тепла, выделяемого (поглощаемого) при диффузионном режиме реагирования для гомогенных и гетерогенных реакций, определяется по соотношениям

(17)
$Q_{{{\Sigma }}}^{v} = \mathop \sum \limits_i \frac{{Q_{{i,w}}^{v}{{C}_{{i,o}}}}}{{{{\vartheta }_{{v,i}}}}},\,\,\,\,Q_{{{\Sigma }}}^{c} = \mathop \sum \limits_i \frac{{Q_{{i,w}}^{с}{{C}_{{i,o}}}}}{{{{\vartheta }_{{с,i}}}}},$
где ${{\vartheta }_{{v,i}}}$, ${{\vartheta }_{{с,i}}}$ – стехиометрические коэффициенты в реакции i-го компонента с химически активной долей вдуваемого газа (2) и углеродом коксового остатка (4); $Q_{{i,w}}^{v}$, $Q_{{i,w}}^{с}$ – тепловые эффекты при реагировании единицы массы химически активного вдуваемого газа и единицы массы углерода коксового остатка с i-м компонентом нагревающего потока при T = Tw.

Соотношение (17) справедливо в предположении вдува химически инертных пиролизных газов.

ИСХОДНЫЕ ДАННЫЕ ДЛЯ РАСЧЕТА

В качестве химически активных компонентов потока нагревающего газа, реагирующих со вдуваемым газом и с коксом ТЗП, рассмотриваются продукты сгорания одного из твердых топлив с массовой концентрацией активных компонентов Н2О и СО2 в продуктах сгорания [4]: ${{C}_{{{{{\text{H}}}_{{\text{2}}}}{\text{O}}}}}$ = 0.0985, ${{С}_{{{\text{С}}{{{\text{О}}}_{2}}}}}$ = 0.0198.

Максимальные массовые концентрации компонентов в случае смеси пиролизных газов при высоких температурах (1500–3200 К) приняты равными [2]

${{g}_{{{{v}_{1}}}}}~ = $ 0.6235$\left( {{\text{C}}{{{\text{H}}}_{4}}} \right)$${{g}_{{{{v}_{2}}}}} = $ 0.3765$\left( {{{{\text{C}}}_{2}}{{{\text{H}}}_{4}}} \right)$.

С использованием соотношений (7), (8) вычислены окислительные потенциалы ВmC и Вmv: по отношению к углероду ВmC = 0.0713; по отношению к пиролизному газу, состоящему из СН4, Вmv = 0.0952; по отношению к пиролизному газу, состоящему из ${{{\text{C}}}_{2}}{{{\text{H}}}_{4}}$, Вmv = $0.0833$; по отношению к пиролизному газу в виде смеси ${\text{C}}{{{\text{H}}}_{4}}$ и ${{{\text{C}}}_{2}}{{{\text{H}}}_{4}}$ с приведенными выше массовыми долями Вmv = $0.0904$.

Необходимые для расчета полных энтальпий (12) исходные данные приведены в табл. 1, 2.

Таблица 1.  

Молекулярная масса М и стандартные тепловые эффекты образования компонентов $H_{f}^{0}$ [9]

Вещество M, 10–3 кг/моль $H_{f}^{0}$, кДж/моль
CH4 16 –66.63
C2H4 28 60.92
H2O(г) 18 –238.91
CO2 44 –393.14
CO 28 –113.81
H2 2 0
С(тв) 12 0
Таблица 2.  

Разность энтальпий веществ $\left( {H_{T}^{0} - H_{{{{T}_{0}}}}^{0}} \right)~$ (12) при различных температурах [9]

Вещество $H_{T}^{0} - H_{{{{T}_{0}}}}^{0}$, кДж/моль
T = 1500 К T = 2500 К T = 3200 К
CH4 90.26 190.94 269.33
C2H4 112.62 230.43 317.99
H2O(г) 58.22 109.63 149.23
CO2 71.09 131.32 174.68
CO 47.52 83.69 109.65
H2 44.76 78.96 104.66
С(тв) 24.37 49.31 68.27

РЕЗУЛЬТАТЫ РАСЧЕТА ТЕПЛОВЫХ ЭФФЕКТОВ И ИХ АНАЛИЗ

Вычисленные согласно (12) на основании данных табл. 1, 2 значения полных энтальпий веществ $\widetilde {~H}_{T}^{0}$, участвующих в реакциях, приведены в табл. 3.

Таблица 3.  

Полные энтальпии веществ $\widetilde {~H}_{T}^{0}$ при нескольких температурах

Вещество $\widetilde {~H}_{T}^{0}$
T = 1500 К T = 2500 К T = 3200 К
кДж/моль кДж/кг кДж/моль кДж/кг кДж/моль кДж/кг
CH4 23.63 1476 124.31 7769 202.69 12 668
C2H4 172.91 6197 291.35 10 405 378.91 13 532
H2O(г) –180.69 –10 038 –129.28 –7182 –89.68 –4982
CO2 –322.06 –7319 –261.84 –5950 –218.46 –4965
CO –66.296 –2367 –30.13 –1075 –4.16 –148
H2 44.76 22 378 78.96 39 481 104.66 52 330
С(тв) 24.37 2030 49.31 4108 68.275 5689

С использованием соотношений (10), (11) выполнены расчеты тепловых эффектов гетерогенных реакций с углеродом прококсованного слоя и гомогенных реакций в пограничном слое. Для оценки влияния состава вдуваемых газов на величины тепловых эффектов при гомогенных реакциях проведены расчеты для двух составов пиролизного газа: состоящего только из метана (CH4) и из метана с этиленом (СН4 + С2H4).

Результаты расчетов для пиролизного газа, состоящего только из GH4, приведены в табл. 4.

Таблица 4.  

Результаты расчетов для реакций с пиролизным газом ${\text{С}}{{{\text{H}}}_{4}}$

Тw, К ${{Н}_{{0,0}}} - {{Н}_{{0,w}}}$, кДж/кг $Q_{{{\text{C}}{{{\text{O}}}_{2}},w}}^{{\text{C}}}$,
кДж/кг
$Q_{{{{{\text{H}}}_{2}}{\text{O}},~w}}^{{\text{C}}}$, кДж/кг $Q_{{{\text{C}}{{{\text{O}}}_{2}},w}}^{v}$, кДж/кг $Q_{{{{{\text{H}}}_{2}}{\text{O}},~w}}^{v}$, кДж/кг $Q_{{{\Sigma }}}^{v}$, кДж/кг $Q_{{{\Sigma }}}^{{\text{C}}}$, кДж/кг $\frac{{Q_{{{\Sigma }}}^{v}}}{{Q_{{{\Sigma }}}^{{\text{C}}}}}$
1500 3950 –13 713 –11 230 –15 962 –14 066 –1347 –811 1.661
2500 1790 –12 659 –10 741 –14 701 –13 287 –1269 –773 1.641
3200 30 –11 792 –10 135 –13 549 –12 301 –1174 –728 1.612
$\vartheta $   3.7 1.5 2.75 1.125      

Примечание. (${{Н}_{{0,0}}} - {{Н}_{{0,w}}}$) – перепад энтальпий по толщине пограничного слоя в отсутствиe гомогенных реакций.

В табл. 5 приведены результаты расчетов тепловых эффектов гетерогенных реакций с углеродом прококсованного слоя и гомогенных реакций в пограничном слое для пиролизного газа, состоящего CH4 (0.6235) + C2H4 (0.3765), в скобках указана массовая концентрация газа в смеси. Использовались соотношения (14), (15)–(17).

Таблица 5.  

Результаты расчетов для смеси пиролизных газов CH4 (0.6235) + C2H4 (0.3765)

Тw, К ${{Н}_{{0,0}}} - {{Н}_{{0,w}}}$, кДж/кг $Q_{{{\text{C}}{{{\text{O}}}_{2}},w}}^{{\text{C}}}$, кДж/кг $Q_{{{{{\text{H}}}_{2}}{\text{O}},~w}}^{{\text{C}}}$, кДж/кг $Q_{{{\text{C}}{{{\text{O}}}_{2}},w}}^{v}$, кДж/кг $Q_{{{{{\text{H}}}_{2}}{\text{O}},~w}}^{v}$, кДж/кг $Q_{{{\Sigma }}}^{v}$, кДж/кг $Q_{{{\Sigma }}}^{c}$, кДж/кг $\frac{{Q_{{{\Sigma }}}^{v}}}{{Q_{{{\Sigma }}}^{{\text{C}}}}}$
1500 3950 –13 713 –11 230 –13 915 –11 931 –1090 –811 1.344
2500 1790 –12 659 –10 741 –12 794 –11 249 –1029 –773 1.33
3200 30 –11 792 –10 135 –11 823 –10 501 –957 –728 1.31
${{\vartheta }_{{{{v}_{{{\text{см}}}}},i}}}$ 2.89 1.185

Результаты расчетов показывают (табл. 4, 5), что отношение тепловых эффектов гомогенных реакций с пиролизными газами к тепловым эффектам гетерогенных реакций с углеродом на поверхности прококсованного слоя $\left( {\frac{{Q_{{{\Sigma }}}^{v}}}{{Q_{{{\Sigma }}}^{{\text{C}}}}}} \right)$ практически не зависит от температуры Tw и существенно меняется при изменении состава вдуваемого пиролизного газа (для чистого метана $\frac{{Q_{{{\Sigma }}}^{v}}}{{Q_{{{\Sigma }}}^{{\text{C}}}}} \approx 1.6$, а для смеси метана и этилена $\frac{{Q_{{{\Sigma }}}^{v}}}{{Q_{{{\Sigma }}}^{{\text{C}}}}} \approx 1.3).$

Отметим, что величины тепловых эффектов при расчетах, выполненных для условных формул ${{C}_{{6.58}}}{{H}_{{20.96}}}$ (т.е. при nсм = 6.58, mсм = 20.96, ${{M}_{{{{v}_{{{\text{см}}}}}~}}}$ = = 100 г/моль) и для ${{M}_{{{{v}_{{{\text{см}}}}}~}}}$ = 50 г/моль, ${{C}_{{3.29}}}{{H}_{{10.48}}}$ (nсм = 3.29, mсм = 10.48) практически совпали, что подтверждает возможность проведения расчетов с использованием формул (15), (16).

Как следует из уравнений (13), (15), (16), значения тепловых эффектов от условной молекулярной массы не зависят. Из сравнения приведенных в табл. 4 и 5 результатов расчетов следует, что замена активного пиролизного газа в виде CH4 на смесь (CH4 + С2Н4) приводит к уменьшению эндотермического эффекта на 18–20%.

ЗАКЛЮЧЕНИЕ

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

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

Предложен метод и выполнены оценки тепловых эффектов гомогенных и гетерогенных реакций в пограничном слое. Из полученных расчетных данных следует, что реакции окисляющих компонентов из нагревающего потока как с активным пиролизным газом, так и с углеродом коксового остатка сопровождаются значительным эндотермическим эффектом. Величина тепловых эффектов составляет от ~800 до ~1300 кДж/кг и оказывается сопоставимой с перепадом энтальпий в пограничном слое (Н0,0Н0,w) в условиях отсутствия гомогенных химических реакций. При этом величина тепловых эффектов существенно зависит от состава пиролизных газов. Так, для рассмотренных составов газа (чистого СН4 и смеси СН4 + С2Н4) различие в тепловых эффектах составляет ~20%.

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

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

  1. Губертов А.М., Миронов В.В., Волкова Л.И. и др. Газодинамические и теплофизические процессы в ракетных двигателях твердого топлива / Под ред. Коротеева А.С. М.: Машиностроение, 2004. 512 с.

  2. Миронов В.В., Толкач М.А. Состав и свойства газообразных продуктов деструкции резиноподобной теплозащиты // ТВТ. 2021. Т. 59. № 5. С. 747.

  3. Кошлаков В.В., Миронов В.В., Чумакин К.А., Толкач М.А. Экспериментальные исследования разбросов теплозащитных характеристик резиноподобной теплозащиты // ТВТ. 2020. Т. 58. № 2. С. 266.

  4. Шишков А.А., Панин С.Д., Румянцев Б.В. Рабочие процессы в ракетных двигателях на твердом топливе. Спр. М.: Машиностроение, 1989. 240 с.

  5. Волчков Э.П., Терехов В.И., Терехов В.В. Структура течения, тепло- и массоперенос в пограничных слоях со вдувом химически реагирующих веществ (обзор) // ФГВ. 2004. Т. 40. № 1. С. 3.

  6. Панкратов Б.М., Полежаев Ю.В., Рудько А.К. Взаимодействие материалов с газовыми потоками / Под ред. Зуева В.С. М.: Машиностроение, 1975. 224 с.

  7. Лукашов В.В., Терехов В.В., Ханьялич К. Теплообмен в пограничном слое на проницаемой пластине при вдуве и горении // Теплофизика и аэромеханика. 2013. Т. 20. № 6. С. 705.

  8. Мануйлов А.В., Родионов В.И. Основы химии. Интернет-учебник / Новосибирск. гос. ун-т. М.: Центрполиграф, 2018. 416 с.

  9. Термодинамические свойства индивидуальных веществ. Спр. изд. В 4-х т. / Под ред. Глушко В.П. М.: Наука, 1978.

  10. Духанин Г.П., Козловцев В.А. Термодинамические расчеты химических реакций. Учеб. пособ. Волгоград: Волгоградск. гос. техн. ун-т, 2010. 96 с.

  11. Кудрявцев В.М. Теория и расчет жидкостных ракетных двигателей. В 2-х т. Т. 2. М.: Высшая школа, 1993.

Дополнительные материалы отсутствуют.