УДК 621.453
ВЛИЯНИЕ ТЕПЛОВЫХ ФАКТОРОВ НА НЕСТАБИЛЬНОСТЬ ГИДРАВЛИЧЕСКИХ ХАРАКТЕРИСТИК КАПИЛЛЯРНЫХ ФОРСУНОЧНЫХ ЭЛЕМЕНТОВ РДМТ
© 2013 В.Е. Нигодюк, А.В. Сулинов
Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет)
Поступила в редакцию 02.12.2013
Проведено расчетно-теоретическое исследование влияния тепловых факторов на гидравлические характеристики капиллярных форсуночных элементов РДМТ. Показано, что влияние тепловых факторов может приводить к существенным изменениям гидравлических характеристик форсуночных элементов РДМТ тягой менее 1 Н и, как следствие этого, к нестабильности параметров двигателя и его нерасчетной работе.
Ракетные двигатели малой тяги, капиллярные форсуночные элементы, гидравлические характеристики, тепловые факторы, нестабильность параметров.
Ракетные двигатели малой тяги (РДМТ) менее 1 Н представляют несомненный интерес для разработчиков систем ориентации и стабилизации космических аппаратов (КА), поскольку данные исполнительные органы обеспечивают как высокую точность управления КА, так и вполне удовлетворительную динамику выполнения операций по управлению летательным аппаратом.
РДМТ тягой менее 1 по сравнению с двигателями тягой от 10 до 400 Н имеют свои отличительные особенности, которые необходимо учитывать при их проектировании [1-3]:
- наличие одного смесительного элемента;
- использование в качестве форсунок, прежде всего, струйных форсунок - капилляров с внутренним диаметром (0,15-0,3) мм и с большими относительными длинами - более 20;
- проблемы с обеспечением стабильности гидравлических характеристик капиллярных смесительных элементов из-за возможной облитерации капиллярных каналов и значительных тепловых потоков в смесительную головку;
- сложности в эффективной организации жид-кофазного смешения компонентов топлива;
- значительные потери удельного импульса тяги из-за неадиабатности процессов в камере;
- особые условия организации тепловой защиты смесительной головки и стенок камеры сгорания.
Со снижением уровня тяги РДМТ возрастает относительная доля тепла, отводимого от про-
Нигодюк Валентин Евгеньевич, кандидат технических наук, доцент кафедры теории двигателей летательных аппаратов. Е-mail: ke_src@ssau.ru Сулинов Александр Васильевич, кандидат технических наук, доцент кафедры теории двигателей летательных аппаратов, старший научный сотрудник Научно-исследовательского центра космической энергетики СГАУ. Е-mail: ke src@ssau.ru
дуктов сгорания к стенке камеры двигателя, в частности потери энергии на неадиабатность рабочего процесса в камере сгорания РДМТ тягой менее 1 Н могут превышать (10...15) % [1]. Поэтому одной из возможных причин нестабильности параметров РДМТ тягой менее 1 Н может быть изменение гидравлических характеристик капиллярных форсуночных элементов смесительной головки двигателя вследствие влияния тепловых факторов, обусловленных тепловыми потоками в смесительную головку.
Физическую картину течения жидкости в капиллярах смесительной головки РДМТ, учитывающей все возможные стадии реализации течения с подводом тепла, можно представить в следующем виде (рис. 1) [3]: на входном участке АВ реализуется изотермическое течение жидкости -без подвода тепла; на последующих участках течения с подводом тепла БЬ, на которых можно выделить зоны жидкофазного течения ВС, течения кипящей недогретой жидкости СЬ, состоящей из участков неразвитого СБ и развитого БК пузырькового кипения, а также из участка пленочного кипения КЬ. На рис. 1 показаны изменения параметров компонента вдоль оси канала капилляра, где I , I , I , I - энтальпия компонента со-
1 7 вх7 ж ст7 s
ответственно на входе в капилляр, по длине канала, при температуре стенки и при температуре насыщенных паров; р - истинное объемное паросодержание; х6 - относительная энтальпия, определяемая по тепловому балансу; хж - относительная энтальпия компонента, средняя по сечению; хд - действительное массовое паросодержание; q - плотность теплового потока; 2 - координата вдоль оси канала капилляра.
На основании рассмотренной выше физической картины течения с учетом результатов исследований [3-9] была предложена математическая
_ =
( V
Псп
п
V ''о у
где п = 2,30
' 1 / '
V рео Л )
Пс.
п
- когда реализуется турбулентное течение
/ \0.33x ч 0,25 _ п П
П ст пж
V Пж ) По )
- при течении недогретои жидкости на участке развитого кипения, когда относительная энтальпия, определяемая по тепловому балансу, х6 <0
г \0,7 г л0-6
_ = 1,28 { д ) f грп ^
Ь кип1
д
V Р )
V & - ОР
ж у
- когда относительная энтальпия, определяемая по тепловому балансу, х6 > 0
Рис. 1. Физическая картина течения жидкости в капилляре смесительной головки ЖРДМТ
модель расчета теплообмена и гидравлического сопротивления компонентов топлива в капиллярах смесительной головки при следующих допущениях:
- стационарность течения и теплообмена;
- постоянство плотности теплового потока на участке теплоподвода к капилляру;
- капиллярные форсуночные элементы - гидродинамически гладкие трубы;
- тепловой поток вдоль оси капилляра пренебрежимо мал.
Изменение коэффициента гидравлического сопротивления при подводе тепла оценивалось с помощью относительного коэффициента гидравлического сопротивления _ = —, где — - коэффициент гидравлического сопротивления при изотермическом течении компонента на ¿-том участке; _ - коэффициент гидравлического сопротивления при подводе тепла к компоненту на ¿-том участке.
Значения относительного коэффициента гидравлического сопротивления для различных участков течения компонента в соответствии с представленной выше физической картиной определялись с помощью следующих выражений:
- при жидкофазном течении и течении недо-гретой жидкости на участке неразвитого кипения, когда реализуется ламинарное течение
_кип 2
П кип
Пж
2Рк
— 2
„ 1 + (1 - Х6 )2 1
—2---'--2----
V Рп (1 -Я>) Рж Р
1
Зная изменения коэффициента гидравлического сопротивления на каждом из рассмотренных участков течения с подводом тепла можно найти изменение суммарного относительного коэффициента гидравлического сопротивления _Е на капиллярном форсуночном элементе:
- + ,.......1....... _ 1
/Е
_ - / - „ - „
р _ ^ о I г I р кип1 I р кип2
ЬЕ _ - ^Ъг 1 ^ Ькип\ 1 ^Ькип2
/Е
/Е
_оЕ "Е "Е "Е
В приведенных выше уравнениях были использованы следующие обозначения параметров: д - внутренний диаметр капилляра; I - длина участка капилляра; п - динамический коэффициент вязкости; Ре - число Пекле; р - плотность; г -скрытая теплота парообразования; ¿ - энтальпия; и - скорость. Индексы в этих уравнениях означают: о - изотермическое течение; Ь - изотермическое жидкофазное течение с подводом тепла; ж и ст - параметры, определяемые при средне-массовой температуре жидкости и температуре стенки; п - пар; 5 - параметр, определяемый при температуре насыщенных паров; кип - параметры участка кипения; 1 и 2 -участки развитого кипения соответственно при х6 < 0 и х6 > 0; Е -суммарный параметр для всей длины капилляра.
С помощью предложенной модели было проведено расчетное исследование изменения гидравлического сопротивления капилляров на примере двухкомпонентного РДМТ номинальной тягой 0,4 Н на компонентах топлива несимметричный диметилгидразин (НДМГ) и азотный тетраоксид (АТ) и, как следствие этого, изменений массового соотношения компонентов топлива и тяги двигателя. Длина и внутренний диаметр капиллярных форсуночных элементов смесительной головки РДМТ, состоящей из одного капиллярного форсуночного элемента по каждой
о
0.25
=0
0.062
линии подачи компонентов, в расчетном исследовании принимались равными соответственно 52 мм и 0,2 мм. Длина выходного участка тепло-подвода капилляра составляла 50 % от длины капилляра - 26 мм. Расчеты были проведены при номинальных параметрах двигателя: тяга 0,4 Н, давление в камере сгорания 0,5 МПа, температура компонентов на входе в двигатель 20 °С, массовое соотношение компонентов топлива 1,8. При этом в капилляре окислителя для двигателя тягой 0,4 Н была возможна реализация как ламинарного, так и турбулентного режима течения компонента, поэтому рассматривались оба режима течения, а в капилляре горючего - только ламинарный режим течения компонента. Характер режимов течения был обусловлен значением числа Рейнольдса при возможных расходах компонентов в капиллярах.
В процессе расчетного исследования определялась зависимость изменения гидравлического сопротивления ^ форсуночного элемента от значения теплового потока (, подводимого к форсуночному элементу. Максимальное значение теплового потока в рассматриваемом исследовании ограничивалось значением 1,2% от кинетической энергии газовой струи на срезе сопла двигателя «Ш» при ее полном расширении. Данное значение теплового потока было достаточным, как показано ниже по тексту, для существенного влияния тепловых факторов на гидравлические характеристики форсуночного элемента. При расчетном исследовании предполагались равные значения тепловых потоков к капиллярам обоих компонентов: горючего и окислителя, постоянство перепада давления на смесительной головке и удельного импульса тяги двигателя.
На рис. 2 представлены результаты расчетного исследования в виде зависимостей суммарного относительного коэффициента гидравлического сопротивления ^ , массового соотношения компонентов Кт и тяги Р двигателя от значения подводимого к капилляру теплового потока (. Зависимости тяги Р и массового соотношения компонентов К , показанные на рисун-
т7 ^
ке, рассчитаны для случая максимального отклонения гидравлического сопротивления - при ламинарном режиме течения окислителя и горючего. Из анализа полученных расчетных данных можно отметить то, что при жидкофаз-ном течении компонентов топлива с подводом тепла уменьшение гидравлического сопротивления при ламинарном режиме течения компонентов составляет 30 % (1,0% от « Ш», (= 12 Вт) и 15 % (0,4% от "Ш", (=5 Вт) соответственно в капиллярах горючего и окислителя, при этом массовое соотношение компонентов топлива К меняется
Рис. 2. Зависимости суммарного относительного коэффициента гидравлического сопротивления , массового соотношения компонентов К
7 т
и тяги Р двигателя от значения подводимого теплового потока: 1, 2 - окислитель; 3 - горючее; 1 - турбулентный режим течения; 2,3 - ламинарный режим течения
незначительно, а тяга Р двигателя увеличивается на 9,4% от номинального значения. Возможность реализации турбулентного режима течения в капилляре окислителя для ЖРДМТ тягой 0,4 Н в сравнении с ламинарным повышает стабильность гидравлического сопротивления в области жидкофазного течения с подводом тепла: уменьшается лишь на 6% и расширяет диапазон данной области по тепловому потоку до 0,7% от «Ш» ((=8 Вт).
Дальнейшее увеличение теплоподвода к капиллярам смесительной головки приводит к образованию паровой фазы, прежде всего в капилляре окислителя, т.к. АТ имеет более низкую температуру насыщенных паров, чем НДМГ, что ведет к росту гидравлического сопротивления: в рассматриваемом диапазоне подводимого теплового потока гидравлическое сопротивление при достижении теплового потока 1,2% от "Ш" по линии окислителя увеличивается в два раза, по линии горючего уменьшается на 25%.
При этом массовое соотношение компонентов топлива и тяга двигателя падают соответственно на 40% и 10%. Однако коэффициент теплоотдачи в капилляре окислителя больше, чем в капилляре горючего, поэтому, по всей видимости, большая часть подводимого к смесительной головке теплового потока будет сниматься окислителем и, учитывая тот факт, что АТ более чувствителен к теплоподводу, то возрастает возможность возникновения кипения в капилляре
окислителя. Последнее будет приводить к перераспределению теплового потока, что существенно увеличит гидравлическое сопротивления капилляра окислителя и приведет к более значительному изменению параметров РДМТ. Кроме того, в области развитого кипения градиент роста гидравлического сопротивления при изменении теплового потока достаточно большой, поэтому незначительные тепловые возмущения могут приводить к неустойчивости рабочего процесса как в капилляре окислителя, так и РДМТ в целом. Поскольку массовое соотношение компонентов топлива изменяется значительно: Кт уменьшилось с 1,8 до 1,0, то будет существенно падать удельный импульс тяги двигателя и, как следствие этого, более значительно по сравнению с расчетным значением уменьшится тяга двигателя.
Проведенное расчетно-теоретическое исследование показало, что влияние тепловых факторов может приводить к существенным изменениям гидравлических характеристик форсуночных элементов РДМТ тягой менее 1 Н и, как следствие этого, к нестабильности параметров двигателя и его нерасчетной работе. Из сказанного выше следует, что на стадии проектирования РДМТ необходимо оценивать возможные тепловые потоки к форсуночным элементам смесительной головки, конструктивно обеспечивая их допустимые значения. Одним из возможных способов снижения тепловых потоков в смесительную головку можно рекомендовать предкамерное устройство в качестве теплового моста между камерой сгорания и смесительной головкой [10].
СПИСОК ЛИТЕРАТУРЫ
1. К вопросу о потерях в камере сгорания двигателей малой тяги / В.Е. Годлевский, В.Е. Нигодюк, А.В. Су-линов //Известия ВУЗов. Авиационная техника. 1983. № 1. С. 77...79.
2. Нигодюк В.Е., Сулинов А.В. Особенности организации рабочего процесса ЖРДМТ тягой менее 1 Н // Материалы докладов на международной научно-технической конференции «Проблемы и перспективы развития двигателестроения». Часть 1. Самара: СГАУ, 2009, С. 119-120.
3. Особенности расчета гидравлических характеристик капиллярных форсуночных элементов ЖРДМТ / В.Е. Годлевский, В.Е. Нигодюк, А.В. Сулинов// Вестник СГАУ. 2009. № 3 (19). С. 241-247.
4. Сточек Н.П., Шапиро А.С. Гидравлика ЖРД. М.: Машиностроение, 1978. 118 с.
5. Теплообмен в ядерных энергетических установках / Б.С. Петухов, Л.Г.Гемин, С.А. Ковалев. М.: Атомиз-дат, 1974. 407 с.
6. Локальное гидравлическое сопротивление при поверхностном кипении воды в трубах / Н.В. Тарасова, В.И. Хлопушин, Л.В. Боронина //Теплофизика высоких температур. 1967. № 1 (5).С.130...136.
7. Орнатский, А. П. Обобщение опытных данных по гидравлическому сопротивлению при поверхностном кипении // Журнал прикладной механики и технической физики. 1965. № 3 (3). С.444...451.
8. Обобщение данных по теплопередаче при кипении недогретых жидкостей / Ю.Е. Похвалов, И.В. Кронин, И.В. Курганов // Теплоэнергетика. 1966. № 5. С.63...68.
9. Богданов Ф.Ф. Влияние скорости парожидкостной смеси теплоносителя и паросодержания на коэффициенты теплопередачи при кипении воды в трубах // Атомная энергия. 1970. Вып.6 (29). С.454...456.
10. Нигодюк В.Е., Сулинов А.В. Перспективы применения предкамер в ЖРДМТ на самовоспламеняющихся компонентах топлива // Материалы докладов на международной научно-технической конференции «Проблемы и перспективы развития двигателестроения». Часть 1. Самара: СГАУ, 2009. С. 120-122.
EFFECT OF THERMAL INSTABILITY FACTORS ON THE HYDRAULIC CHARACTERISTICS OF THE CAPILLARY INJECTOR ELEMENT LIQUID ROCKET THRUSTERS
© 2013 V.E. Nigodjuk, A.V. Sulinov
Samara State Aerospace University named after Academician S.P. Korolyov (National Research University)
By calculation and theoretical study of the influence of thermal factors on the hydraulic characteristics of the capillary injector element liquid rocket thrusters. It is shown that the influence of thermal factors can lead to significant changes in the hydraulic characteristics of injector components of rocket engine thrust of 1 N and, consequently, to the instability of the engine and its nonisobaric work. Rocket thrusters, air washers capillary elements hydraulics, thermal factors, instability parameters.
Valentin Nigodjuk, Candidate of Technical Science, Associate Professor at the Theory of Aircraft Engines Department. E-mail: ke_src@ssau.ru
Alexander Sulinov, Candidate of Technical Science, Associate Professor at the Theory of Aircraft Engines Department, Senior Research Fellow, Research Center of Cosmic Energy of SGAU. E-mail: ke src@ssau.ru