- •Под готовностью подразумевается возможность расчета одного варианта в течение суток на самых мошных из доступных компьютеров.
- •Имеется в виду les с пристеночным rans моделированием: в случае les вплоть до твердых стенок, затраты оказываются сопоставимыми с затратами dns.
- •На компьютере с производительностью 1 терафлоп. Время расчета составляет 5000 лет!
- •Рже. 3. Профиль безразмерной избыточной скорости в пограничном слое двух плоских спутных струй воздуха (начальный участок)
- •Физическая картина гетерогенных течений
- •Математические модели двухфазных потоков
- •Гетерогенные течения в ракетных двигателях
- •Математическая модель двухфазных потоков
- •Г етерогенные течения в ракетных двигателях
- •Коэффициент избытка окислителя, (X
- •- Условно непроницаемая граница ядра потока:
- •- Кусочно-равномерное распределение параметров потока в критическом сечении;
- •Математическая модель гетерогенных течений
- •2.Особенность моделирования многофракционности состава к-фазы
- •Лекция 12 аэродинамический нагрев
- •Понятие сопряженной задачи аэродинамического нагрева
- •Особенности аэродинамического нагрева при гиперзвуковых скоростях
- •Вычислительный эксперимент как метод теоретического анализа аэрогазодинамики ракет
- •Газодинамические процессы
- •Методы дискретизации уравнений газовой динамики по пространству и времени, разностные схемы
- •Вычислительные сетки
- •Обзор и сравнительный анализ численных методов для решения задач аэрогазодинамики
- •Представление данных
- •5 Разностная схема метода крупных частиц для расчета движе
Г етерогенные течения в ракетных двигателях
Присутствие дисперсной примеси в потоках продуктов сгорания даже при незначительной концентрации последней может приводить к нежелательным последствиям для тракта двигателя. Вследствие этого изучение такого рода течений и построение методов расчета гетерогенных сред в РДТТ приобретает чрезвычайную актуальность.
Присутствие твердых частиц может приводить к значительному (порой многократному) увеличению тепловых потоков, а также к эрозионному износу обтекаемой поверхности. Эти явления обусловлены совместным действием целого ряда причин, среди которых изменение структуры течения набегающего на тело потока, а также характеристик пограничного слоя, развивающегося на обтекаемом теле, соударения частиц с поверхностью, изменение шероховатости поверхности и многое другое. Интенсивность процессов, сопутствующих обтеканию тел гетерогенными потоками, зависит от инерционности и концентрации частиц. Следует отметить, что инерционность частиц напрямую определяется геометрией и параметрами течения и может изменяться дня одних и тех же частиц в очень широких пределах. Наличие различных характерных времен (длин) несущего потока (вблизи критической точки обтекаемого тела, вдоль его поверхности, собственно турбулентных масштабов и т.д.) сильно осложняют изучение таких потоков и обобщение данных. Что касается концентрации частиц, то ее значение может многократно превышать «исходное» значение в невозмущенном потоке из-за резкого торможения потока при приближении к телу, взаимодействия частиц со стенкой, а также межчастичных столкновений. При движении частиц вдоль поверхности тела в пограничном слое, где имеются значительные градиенты скорости и температуры (в случае неизотермического течения), их распределение зачастую носит сложный характер, а концентрация также превышает свое значение в набегающем на тело потоке.
Добавление в турбулентный ноток (и без того сложный) дисперсной примеси в виде частиц сильно осложняет картину течения. Прежде всего это связано с большим разнообразием свойств вводимых частиц, которое приводит к реализации многочисленных режимов течения газовзвеси. Варьирование концентрации частиц — основной экстенсивной характеристики гетерогенных потоков — позволяет не только изменять количественно параметры исходного течения и движения частиц, но приводить и к его качественной перестройке (например, переходу ламинарного режима течения в турбулентное, а также к обратному эффекту, т.е. реламинаризации течения). Вследствие этого методы экспериментальных и теоретических исследований, используемые в классической механике однофазных сплошных сред, зачастую не могут быт ь использованы для изучения гетерогенных потоков в принципе. Имеющиеся экспериментальные данные зачастую носят отрывочный и противоречивый характер, а физические представления и развитые математические модели не могут быть признаны удовлетворительными. Сказанное выше сдерживает развитие механики гетерогенных сред. Несмотря на это, потребности практики и логика развития науки настойчиво требуют постоянного совершенствования теории гетерогенных течений.
К основным задачам изучения обтекания тел гетерогенными потоками относятся: 1) исследование движения частиц и определение их траекторий; 2) установление влияния частиц на течение газа; 3) исследование процессов взаимодействия частиц с обтекаемой поверхностью, включая эрозию материала.
Движение продуктов смесевых топлив по тракту РДТТ
Предположим что перед входом в дозвуковую часть сопла газовой и k-фазы движутся с одинаковой скоростью. На периферийной части поток перед входом в дозвуковую часть сопла тормозится и разворачивает движение вдоль стенок сопла, но тяжелые частицы k-фазы не успевают затормозить и развернуться, потому ударяются в стенку конической части входа в сопло и часть из них прилипает к поверхности за счет охлаждения, при этом если конструкция не имеет теплозащиты в этой зоне, то может произойти прогар стенки (разрушение) за счет высокой температуры и механического воздействия. При дальнейшем движении по соплу в зоне критического сечения газовый поток вновь разворачивается и продолжает движение вдоль сопла, а твердые частицы образуют уплотнительную зону около оси сопла. По мере дальнейшего течения по сверхзвуковой части сопла, частицы k-фазы увлекаются потоком газа и разворачиваются в сторону периферийной части сопла и если сопло профилированное и достаточно длинное, эти частицы вновь ударяются в стенку сопла в зоне выходного сечения и могут приводить к разрушению сопла в сверхзвуковой выходной части. В связи с этим при проектирование сопел для смесевых топлив необходимо учитывать эти явления и тщательно профилировать большими радиусами дозвуковую часть сопла, сглаживая все углы. Так же необходимо обеспечивать теплозащиту вероятных зон взаимодействия частиц с конструкцией. Не использовать длинные профилированные сопла, а по возможности укороченные.
В общем случае частицы k-фазы могут иметь произвольную форму. Для продуктов сгорания смесевых топлив частицы k-фазы имеют высокую температуру и могут находиться в расплавленном состоянии, их форма близка к сферической. Поэтому в дальнейшем при составлении движения гетерогенного потока будем использовать допущение о сферичности частиц k-фазы. Расчет же несферической формы (удлиненной или игольчатой ) чрезвычайно сложен, поэтому и моделирование этих частиц сложно.
Расчет сил между k-фазой и газовой фазой
Необходимо определить силу, действующую на сферическую частицу при относительном ее движении в потоке газа. Если скорости частиц и газа одинаковы, то такая скорость равна нулю.
Относительная скорость газа: V0 = Vr - VK
Vr
Рис.1.
Схема обтекания частицы к-фазы
Так как плотность материала k-фазы намного выше плотности газа на несколько порядков, то инерционность частички будет значительно выше инерционности газовой фазы в связи с этим, там, где газовый поток ускоряется, частички k-фазы будут отставать, а в зонах торможения наоборот будут опережать газовый поток (рис.2).
Рис.
2. Схема смешения потоков
На рисунке 3 показано обтекание цилиндра гетерогенным потоком частицами различного размера, хорошо иллюстрирующее влияние инерционности частиц на траектории их движения. Эти особенности течения позволяют спрогнозировать характер движения к-фазы в двигателе. На рисунке 4 показана схема движения частиц по тракту двигателя.
Рис.3.
Траектории частиц диаметром 10 мкм при
поперечном потенциальном обтекании
цилиндра гетерогенным потоком частицами
с плотностью р= 1400 кг/м со скоростью U
=
6 м/с при различных диметрах частиц:
а)
10 мкм; б) 20 мкм; в) 100 мкм
зона
налипания к-фазы
Рис.
4. Схема движения частиц к-фазы по тракту
двигателя
Из рисунков следует, что тракт двигателей на смесевых топливах в зонах разворота потока продуктов сгорания, сопла двигателей (особенно в зоне выходного сечения) подвергаются очень интенсивному абразивному износу. На степень износа влияют, главным образом, материал сопла, материал частиц и скорость движения частиц.
В общем случае силы, действующие на частицу, определяются относительной скоростью, диаметром частицы, плотностью газа обтекающий частицу и характером течения (ламинарное и турбулентное). Необходимо учитывать, что характер течения ламинарный или турбулентный, определяется абсолютной скоростью газа и внешним условиями течения, т.к. сам поток еще до встречи с частицей может быть турбулентным. Т.к. во внутренних течениях (камере РДТТ) мы имеем дело с турбулентными потоками (течениями), будем рассматривать именно этот вариант. Сила складывается из нескольких составляющих:
F — F + F + F
к лобового_ сопротивления тр кр
Максимальная сила Fno6oeozo_сопротиештя, а остальным будем пренебрегать или учитывать их коэффициент, или учитывать их для сферической формы.
P • S • V7
г
Fk —m C
в m2
Сх - коэффициент лобового сопротивления, он учитывает форму тела, характеризующую течения и дополнительные силы.
Fk — 1 •XP. •s2V2 -Vk|•(Vr -Vk)• C,
Для определения Сх могут быть проведены экспериментальные исследования или теоретические расчеты. Такие исследования проведены и получены следующие зависимости:
С =
>
R =
e
4,3(^gRe )"2,700< Re <2000 V • D p• V • D
24 4
— + —^,0<Re < 700
R R03 e
<
v и
Данная величина силы должна добавляться со знаком минус для газовой фазы, и со знаком плюс для k-фазы.
Профилирование сопел для двухфазных продуктов сгорания
В настоящее время в РДТТ широко применяются топлива с добавкой металлов. Применение металлов в качестве компонентов ракетного топлива может обеспечить увеличение удельного импульса, связанное с высокой теплотой горения металла, а также с уменьшением молекулярной массы
T
продуктов сгорания( i = ). А также при добавке металлов повышает-
V И
ся плотность топлива. Продукты сгорания металла в камере двигателя находятся в конденсированной фазе. Поэтому при работе двигателей на металлосодержащих топливах продукты сгорания представляют собой двухфазное рабочее тело - смесь газа и конденсированных частиц.
При течении в сопле ускорение частиц происходит лишь под действием аэродинамических сил, возникающих при обтекании их газом. Поэтому при течении двухфазных потоков в сопле всегда имеет место отставание частиц по скорости <ат. Понижение температуры частиц происходит
лишь вследствие теплообмена с газом, поэтому Тч > Т. Таким образом при двухфазном течении происходит обмен количеством движения и энергией между конденсированными частицами и газом. Скоростная неравно- весность двухфазного течения в сопле (< ог) является основной причиной дополнительных потерь удельного импульса в сопле; наличие температурной неравновесности на потери удельного импульса влияет в гораздо меньшей степени, к тому же благодаря теплообмену между конденсированными частицами и газом частично восполняются затраты тепловой энергии газового потока на разгон частиц конденсированной фазы.
Задача профилирования сопла для случая течения двухфазных продуктов сгорания, как и для гомогенных, заключается в отыскании контура с наибольшим удельным импульсом. Усложнение задачи состоит, во- первых, в необходимости применять более сложную по сравнению с вариантом гомогенных продуктов сгорания систему уравнений. Во-вторых, потери удельного импульса, обусловленные запаздыванием частиц, зависят от геометрии сужающейся части и горловины сопла. Наконец, экспериментальные и расчетные данные указывают на возможность попадания частиц на концевой участок профилированного сопла, что приводит к дополни
тельным потерям удельного импульса и должно приниматься во внимание при выборе контура. Кроме того, попадание частиц на концевой участок сопла может привести к его разрушению.
Решение вариационной задачи поиска оптимального для двухфазных неравновесных течений контура сопла показывает, что в случае оптимального профиля потери удельного импульса из-за запаздывания частиц на 30% меньше по сравнению с потерями в сопле, контур которого построен для гомогенных продуктов. При этом сопло имеет необычную конфигурацию: короткий вход и сильно вытянутую горловину. Последнее связано с необходимостью уменьшения максимального градиента скорости и перемещения зоны наибольшего отставания в сужающуюся часть.
Рис.
5. Элементы профиля сопла для двухфазных
продуктов сгорания
Рассмотрим основные выводы исследований по выбору оптимального профиля. В случае двухфазного течения может быть рассчитано поле течения равновесной двухфазной смеси, и затем в полученном поле параметров газовой фазы рассчитываются траектории частиц. Это дает возможность приближенно определить интенсивность выноса частиц на стенку сужающейся части соответствующим образом подкорректировать контуры, чтобы избежать выноса частиц. Одновременно вычисляются параметры газовой и к-фазы в сечении хо, в котором скорость газа превышает местную скорость звука в газе. Очертание дозвуковой части может быть выбрано, как это делается для гомогенных продуктов сгорания, либо на основе конструктивных соображений. При этом учитывается, что удлинение области горловины не оказывает существенного влияния на снижение двухфазных потерь вследствие коагуляции (слияния) частиц, а потери из-за трения возрастают. Радиус R выбирается возможно более меньшим из условий обеспечения безотрывного течения и стойкости конструкционных материалов.
Для расчета неравновесного двухфазного течения в сверхзвуковой части сопла применяется модифицированный метод характеристик, который позволяет рассчитывать параметры потока в сечениях х сопла последовательно, начиная с сечения хо. Строится оптимальное укороченное сопло для газа, параметры которого совпадают с параметрами равновесной двухфазной смеси, и с помощью расчетного определения удельного импульса для любых выбранных углов вх ,02, обеспечивающих максимальный удельный импульс.
В результате серии расчетов неравновесных двухфазных течений для случая ё=1...10мкм были получены (^ ,#2)опт. Оказалось, что полученные контуры совпадают с контурами, оптимальными для равновесной двухфазной смеси. Этот результат позволяет использовать разработанные для гомогенных продуктов сгорания методы профилирования сопел и для двухфазных смесей. Другой существенный результат состоит в установлении слабого изменения уд. импульса при значительном отклонении параметров контура (0J ,#2) от оптимальных. Так для одного из расчетов были
получены оптимальные значения вх = 250 ,#2 = 8°. Варьирование вх в диапазоне 17-30° уменьшает уд. импульс на 0,02-0,03%, а варьирование (раздельное) 02 в диапазоне 4-12о уменьшает уд. импульс на 0,1%.
Для выбора профиля, оптимального при учете всех рассмотренных факторов, требуется проведение численных экспериментов, позволяющих количественно оценить роль каждого из них.
Течение с частицами в области критической точки тела
Для анализа физических процессов, происходящих в области лобовой точки обтекаемого гетерогенным потоком тела, необходимо знание распределений скоростей, температур (в случае неизотермического течения) и концентрации частиц. На рис. 6 представлена классификация возможных режимов обтекания гетерогенным потоком.
Важной интегральной характеристикой гетерогенных течений при обтекании тел является коэффициент осаждения (улавливания) частиц. Данный коэффициент представляет собой отношение претерпевших соударение с телом частиц к числу частиц, которые могли с ним столкнуться, если бы их линии тока представляли прямые линии.
В случае осесимметричного течения (например, при обтекании сферы), когда размер частиц пренебрежимо мал по сравнению с размером тела, частицы равномерно распределены в набегающем потоке, а траектории их движения симметричны, коэффициент осаждения можно определить как
Л = Й,
где ycr = ycr/R - безразмерное расстояние от оси симметрии потока, находясь на котором (в невозмущенном присутствием тела потоке), частицы лишь соприкасаются с телом при его обтекании. Частицы. координаты которых в набегающем потоке У > ycr не претерпевают столкновения с телом.
Рис.
6.
Классификация режимов течении при
обтекании тел гетерогенным потоком:
а) равновесное течение, квазиравновесное
течение; 6)
слабо запыленное неравновесное течение,
течение с крупными частицами; в)
сильно запыленное неравновесное
течение, течение с крупными частицами;
г) обтекание в случае поглощающей
стенки; д) обтекание с образованием
пленки.
Расчет потерь в соплах ракетных двигателей
Рассмотрим стандартизованные определения потерь удельного импульса из-за отличия реальных рабочих процессов в сопле от их идеального представления.
Количественные оценки потерь, рассчитываемые по приведенным зависимостям, отражают влияние основных физических механизмов проявления неидеальностей, но представляют собой только первое приближение. Рекомендуемые стандартами существенно более сложные методики расчета потерь значительно точнее учитывают реальную структуру течения рабочего тела в сопле и достоверность получаемых по ним данных достаточна для их применения при проектировании ракетных двигателей. При этом методики, как правило, ориентированы на ракетные двигатели конкретных видов и учитывают их особенности.
Рассмотрим методику количественной оценки потерь, свойственных двигателям разных типов.
Потери удельного импульса из-за химической нераеноеесности
Коэффициент сопла =
Потери удельного импульса из-за химической неравновесности (коэффициент сопла фс1) и из-за неравномерности свойств потока в минимальном сечении сопла (фс2) имеют термодинамическую природу - цикл расширения рабочего тела реализуется так, что часть располагаемой полной энергии химического топлива не переходит в тепловую энергию рабочего тела и, затем, не преобразуется в кинетическую энергию истекающей струи.
Расчет коэффициента снижения удельного импульса из-за химической неравновесности течения фс1 ведется на основе следующих представлений:
Протекание в химически активном рабочем теле реакций рекомбинации обеспечивает распределенный по длине сопла переход химической энергии рабочего тела в тепловую энергию в процессе расширения. При этом затраченная на диссоциацию в камере сгорания часть высвобождающейся при химических реакциях тепловой энергии возвращается в виде тепловой энергии, подводимой изнутри к ускоряющемуся в сопле потоку.
Существенными факторами, действие которых обеспечивает реально достигаемую степень равновесности процесса, является ограниченность скоростей физико-химических процессов, идущих в рабочем теле при течении его по соплу. Условно будем считать, что рабочее тело - газ, имея в виду, что возможная неравновесность, свойственная гетерогенным рабочим телам, будет учитываться как составляющая потерь из-за наличия в потоке газа конденсированных (твердых или жидких) частиц.
Коэффициент сопла ф сЬ характеризующий рассматриваемую неиде- альность, определяется как отношение действительного коэффициента тяги в пустоте ^т.п.д при реализующейся степени химической неравновесности, к идеальному ^т.п.и, вычисленному в рамках допущения о полном равновесии в потоке
^ГПД
^с1 — = ^тпи . (1)
При этом свойства рабочего тела в критических сечениях считаются одинаковыми, одинаковы также и геометрические степени расширения сопла.
При таких условиях сравнения равны тяги, создаваемые камерами на участках от головки камеры сгорания до минимального сечения сопла, а отношение коэффициентов тяги равняется отношению развиваемых пустотных удельных импульсов
v F I
-г]ф 1ф _ упд
Аши ^
V F
Фс1- ^
При этом будем считать, что идеальное и реальное течения одномерны и безотрывны.
В рамках настоящей дисциплины рассматриваются потери удельного импульса из-за химической неравновесности течения - фс1, проявляющейся в сопоставимой мере в соплах РДТТ и ЖРД.
Термодинамический расчет дает значения удельного пустотного импульса вне зависимости от пути перехода рабочего тела из произвольного состояния к параметрам в камере при заданной модели процесса в сопле для конкретных:
состава рабочего тела,
давления в камере ,
степени расширения газа в сопле
или геометрической степени расширения сопла ^ .
Модель идеального процесса - равновесное изоэнтропное расширение. Реальный процесс идет с возрастанием энтропии, что отражает необратимость явлений переноса при конечных градиентах концентраций реагентов в потоке и, как следствие, при конечных значениях скоростей химических реакций. Однако конечная, но равная нулю скорость всех возможных химических реакций в потоке (модель замороженного течения) соответствует случаю полного отсутствия возврата затраченной на диссоциацию энергии и, соответственно - минимальному значению коэффициента сопла фс1.
Учет реальной степени завершенности химических реакций требует привлечения аппарата описания химической кинетики превращений. В рамках настоящего курса рассматриваются только численные оценки первого приближения на основе так называемых локально предельных допущений: равновесное течение в дозвуковой области и химически замороженное течение (скорость химических реакций равна нулю, состав рабочего тела постоянен) в сверх звуковой части сопла. И все анализируемые расчетные данные по оценке влияния химической неравновесности получены сравнением идеальных параметров двигателя с результатами расчета в предположении локальной предельной химической неравновесности.
Оценим влияние давления в камере на потери пустотного удельного импульса из-за химической неравновесности при разных геометрических
К,
г- F
степенях расширения сопла r0TH= ^ =var на примере кислород- водородного двигателя (коэффициент избытка окислителя а =0,8). Причем,
для каждого давления ^принимается такая геометрическая степень расширения сопла, чтобы давление в выходном сечении сопла сохранялось одним и тем же ра=0,07 МПа.
Расчеты проведены в равновесном и локально-неравновесном приближении. Причем под локально-неравновесным приближением понимается равновесное течение в сужающейся, и замороженное - в расширяющейся части сопла (так называемое "приближение Брея").
Результаты расчетов сведены в таблицу 12.
Таблица 12
К оценке потерь на неравновесность
р МПа ■i |
I м/с ул, |
п |
I м/с ул, |
п |
*cl |
Равновесное течение |
Локально неравновесное течение |
||||
. 5 |
• 4200 |
1=200 |
4140 |
• 1=220 |
0=9S5 |
• 10 |
• 4340 |
1 =210 |
4290 |
• 1=225 |
0=987 |
• 15 |
• 4420 |
1 =217 |
4370 |
• L22S |
0=9S9 |
• 20 |
• 4460 |
1=220 |
4420 |
. 1=230 |
Q=990 |
• 25 |
• 4500 |
1=223 |
4460 |
• 1=232 |
0=991 |
• 30 |
• 4520 |
1=225 |
4480 |
• 1=233 |
0=991 |
• 35 |
• 4550 |
1=227 |
4510 |
• 1=234 |
0=992 |
Расчеты показывают, что с повышением давления в камере в диапазоне 5... 35 МПа с ростом удельного импульса от 4200 м/с до 4540 м/с абсолютное значение потерь уменьшается от 6 до 3 м/с, что отражается в росте коэффициента сопла от 0,985 до 0,992.
Сравним значения условных осредненных показателей изоэнтропы п в равновесном правн и локально предельно неравновесном Пн/равн вариантах протекания процессов расширения рабочего тела в сопле. Во всем рассмотренном диапазоне изменения давления соотношение правн<пн/равн, что отражает приближение равновесного процесса к изотермическому, для которого п=1. При этом в потоке внутренний теплоподвод осуществляется за счет протекания реакций рекомбинации и соответствующего выделения тепла. В составе рабочего тела возрастает концентрация компонентов с относительно большой молекулярной массой Н2°и др. И уменьшается содержание СО ОН н и др При этом возвращается тепловая энергия, затраченная в камере сгорания на температурную диссоциацию, температура равновесного потока возрастает по сравнению с потоком, в котором химические превращения не идут.
На графике рис. 1 представлены зависимости коэффициента сопла ф с1 от давления в камере, полученные решением задачи термодинамического расчета для равновесного и неравновесного в сверхзвуковой части сопла течения. Сплошными кривыми 1, 2 и 3 соответствуют расчеты для постоянных геометрических степеней расширения сопла, кривой 4 - расчеты для постоянного давления в выходном сечении сопла ра=0,07 МПа и переменной геометрической степени расширения.
0
0 10 20 30 40
Давление в камере, МПа
Рис. 1. Влияние давления в камере и геометрической степени расширения сопла на уровень потерь из-за химической неравновесности:
для сопла с геометрической степенью расширения 1000:
для сопла с геометрической степенью расширения 500;
для сопла с геометрической степенью расширения 250;
4 - для сопла с постоянным давлением в выходном сечении сопла
ра=0 07 МПа и переменной геометрической степенью расширения
,995й Ц
о 0,99 о
f- Я
Я 0,985 ЕГ S
-В*
'f 0,98
О
М
0,975
Из графиков 1 следует, что при увеличении давления в камере относительный уровень потерь в сопле на химическую неравновесность снижается.
Стремление к единице фс1 при росте давления в камере отражает снижения влияния возвращения затраченной на диссоциацию в камере сгорания теплоты, т. к. при повышенных давлениях в камере сгорания диссоциация в значительной степени подавляется (согласно принципу Ле Ша- телье).
Для данного давления в камере росту степени расширения сопла соответствует рост потерь на химическую неравновесность, поскольку в соплах с большой степенью расширения имеется возможность большей внутренней регенерации теплоты, т.е. возвращения рабочему телу большей доли тепловой энергии, поглощенной эндотермическими процессами термической диссоциации в камере сгорания.
Эта выявляемая расчетом и подтверждаемая практикой тенденция усиливается в соответствии с увеличением скорости химических реакций с
ростом давления из-за увеличения концентрации реагирующих веществ и температуры в камере.
Приближение удельного импульса к его идеальному равновесному значению с ростом давления наряду с другими факторами объясняет преимущества ракетных двигателей с повышенным давлением в камере.
Конечность скоростей реакций рекомбинации определяет влияние абсолютных размеров двигателя на потери из-за химической неравновес- ности, т.к. при одних и тех же скоростях движения потока время пребывания в сопле будет большим у большеразмерного двигателя, и, соответственно, возникает возможность более полного завершения реакций рекомбинации и снижения потерь. В этой связи при прочих равных потери на химическую неравновесность больше у ракетных двигателей малых тяг (например, жидкостные ракетные двигатели малой тяги - ЖРД МТ), чем у двигателей с уровнем тяг 5 кН и более.
Таким образом, потери удельного импульса из-за химической нерав- новесности в основном связаны с неполнотой или отсутствием возврата затраченной на диссоциацию продуктов реакции горения в камере теплоты при "замораживании" - незавершенности в сопле экзотермических химических реакций рекомбинации, каждая из которых отличается значением теплового эффекта, составом и свойствами продуктов реакции и другими параметрами. В свою очередь, вероятность тех или иных химических реакций в основном определяется химическим составом топлива. Отсюда следует зависимость потерь удельного импульса из-за химической неравно- весности течения в сопле от состава топлива, а для заданной пары - от соотношения расходов горючего и окислителя.
Очевидно, что при низких температурах (больших отличиях соотношения компонентов от стехиометрического), когда реакции термической диссоциации в камере сгорания практически не идут, невозможен и возврат теплоты при равновесном течении рабочего тела в сопле, соответственно не будет и потерь на химическую неравновесность.
Ввиду изменения свойств продуктов сгорания различных топлив в весьма широком диапазоне, выявить однозначную зависимость уровня потерь от состава топлива не представляется возможным. Связь состава топлива и потерь на химическую неравновесность проявляется при определении условно оптимального соотношения компонентов - значения соотношения расходов окислителя и горючего, при котором достигается максимальное значение пустотного удельного импульса для заданной степени расширения.
На графике рис. 2 представлены графики зависимости пустотного удельного импульса /у.п от соотношения расходов компонентов топлива "углеводородное горючее-кислород" (коэффициента избытка окислителя а
при Кт0=3,4119) при давлении в камере ^к=5,2 МПа и давлении на срезе сопла ра=0,034 МПа.
3
3400
(J
'g
2800
3200
0,55
800
|
|
|
|
к |
|
/ |
|
|
4 |
— |
|
|
|
|
2 |
|
|
/ |
|
|
|
|
|
7 |
|
|
3 |
|
|
|
|
|
|
|
|
