CC BY
208
УДК 629.001.036.5.621.22.01
МЕТОД ПРОЕКТНОГО РАСЧЕТА ДВУХКОМПОНЕНТНОЙ СТРУЙНОЙ ГАЗОЖИДКОСТНОЙ ФОРСУНКИ С ВНУТРЕННИМ СМЕШЕНИЕМ
©2011 В. С. Егорычев
Самарский государственный аэрокосмический университет имени академика С.П. Королёва
(национальный исследовательский университет)
Представлен метод проектного расчета двухкомпонентной струйной газожидкостной форсунки с внутренним смешением.
Метод, форсунка, двухкомпонентная форсунка, струйная газожидкостная форсунка, жидкостный ракетный двигатель, смесеобразование, топливо, окислитель, горючее.
Двухкомпонентные струйные газожидкостные форсунки с внутренним смешением используются в ЖРД с дожиганием генераторного газа, которые работают на несамовоспламеняющихся, чаще всего криогенных, компонентах топлива. Схема организации рабочего процесса в камере сгорания таких двигателей „газ - жидкость".
Расчетная схема двухкомпонентной струйной газожидкостной форсунки с внутренним смешением приведена на рис. 1.
Рис. 1. Расчетная схема двухкомпонентной струйно-струйной газожидкостной форсунки с внутренним смешением
Форсунка состоит из наружного и внутреннего контуров. Внутренний контур представляет собой однокомпонентную струйную газовую форсунку окислителя или окислительного генераторного газа. Наружный контур - однокомпонентную струйную
жидкостную форсунку (или форсунки) горючего. Смешение газового окислителя и жидкого горючего начинается уже внутри двухкомпонентной форсунки. Из сопла форсунки вытекает газожидкостная струя топливной смеси.
Оба контура такой форсунки работают и рассчитываются независимо друг от друга, причём расчёт начинается с внутреннего контура.
Проектный расчет двухкомпонентной струйной газожидкостной форсунки с внутренним смешением проводится в следующей последовательности.
1. Рассчитывается внутренний контур двухкомпонентной форсунки, представляющий собой однокомпонентную струйную газовую форсунку окислителя, подающую в камеру сгорания газообразный окислитель или окислительный генераторный газ.
1.1. Выбирается в первом приближении диаметр проходного сечения или диаметр сопла форсунки на выходе с1с в.
Этот диаметр может находиться в пределах 0,5 мм < йсв <Вфв.
В первой итерации его можно принять равным 4.. .6 мм.
1.2. Определяется в первом приближении длина форсунки или длина сопла форсунки 1св.
Газовая струйная форсунка подаёт окислительный генераторный газ из полости смесительной головки, расположенной между средним и наружным днищами. Поэтому при определении её длины следует учитывать толщины среднего и огневого днищ, расстояние между ними и величины выступов форсунки над средним и огневым днищами.
Толщина огневого и среднего днищ находится в интервале 2,0...3,0 мм, а рас-
стояние между днищами составляет 8... 12 мм [1].
1.3. Выбирается форма кромки канала со стороны входа компонента в форсунку.
Форма кромки канала струйной форсунки на входе может быть острой, закругленной по радиусу или с фаской. Газовые форсунки с острой входной кромкой в смесительных головках ЖРД практически не используются. Наиболее часто применяются форсунки с фаской.
1.4. Находится относительная длина струйной форсунки, т.е. отношение 1ав^ав .
1.5. Площадь сопла форсунки на выходе
• 0)
1.6. Число Рейнольдса на выходе из форсунки
4 т
Ке =-(2)
Аае
где тфок - массовый расход окислителя или
окислительного генераторного газа через одну форсунку; цтт - динамическая вязкость газа, подаваемого через форсунку, при температуре на входе в форсунку. В расчёте можно использовать температуру генераторного газа за турбиной.
1.7. Среднее значение скорости генераторного газа на выходе из форсунки
т
цг(3)
Р г.
~ гг ф.в
где ргг - плотность генераторного газа на входе в форсунку.
1.8. Коэффициент расхода форсунки
( СО, V1
/Л =
581 1,23 + - св
/л
(4)
¡л
(5)
Red
V с" У
Решая это уравнение, получим .^1,232 +2327C8/Re dce -1,23 \\61св/КЫсв
1.9. Давление генераторного газа на входе в форсунку
Ре.с =Рк+ &Рф.о, (6)
где рк- давление на входе в камеру сгорания; Арфв- перепад давления на внутреннем
контуре двухкомпонентной форсунки или струйной газовой форсунке окислителя.
1.10. Плотность генераторного газа на входе в форсунку
Ргг =
К С
(7)
где Твх- температура генераторного газа на
входе в форсунку.
1.11. Площадь сопла форсунки на выходе
т.
фок
г
т,
\
i
р \PexJ
2 " АЛ
п-1
f
1-
\
Р
\PexJ
(8)
1.12. Диаметр сопла форсунки, обеспечивающий требуемый массовый расход окислительного генераторного газа или газообразного окислителя при выбранном перепаде давления на форсунке Арфв:
4
ж
(9)
Если диаметр отличается от принятого значения или полученного в предыдущей итерации, то расчёт форсунки повторяется при новом уточнённом значении диаметра.
При отличии значений диаметров менее 1% расчёт прекращается.
2. Рассчитывается наружный контур двухкомпонентной форсунки, представляющий собой однокомпонентную струйную жидкостную форсунку горючего.
Если число струйных жидкостных форсунок горючего в наружном контуре двухкомпонентной форсунки пф, то массовый расход горючего через одну форсунку будет
т
ЩФ.г=■
ф. г
пл
(10)
где тф1 - массовый расход горючего через
наружный контур двухкомпонентной форсунки.
2.1. Выбирается в первом приближении диаметр проходного сечения или диаметр сопла форсунки на выходе dcн.
Этот диаметр может находиться в пределах 0,2 мм < dcн < —, где Н - расстояние
между средним и огневым днищами смесительной головки [1]. Минимальный диаметр сопла форсунки обусловлен легкостью засорения малых проходных сечений форсунки механическими примесями, попавшими случайно в распыляемый жидкий компонент, или твердыми частицами, образовавшимися
в компоненте из-за его физической и химической нестабильности, а также трудностью изготовления тонких отверстий. Максимальный диаметр обусловлен конструкцией смесительной головки камеры и тонкостью распыла компонента.
2.2. Определяется в первом приближении длина струйной жидкостной форсунки или длина сопла форсунки 1ан.
Форсунка подаёт жидкое горючее из полости смесительной головки, расположенной между средним и огневым днищами. Длина форсунки определяется толщиной стенки наружного контура и углом наклона оси жидкостной форсунки к оси газовой форсунки окислителя.
2.3. Выбирается форма кромки канала со стороны входа компонента в форсунку.
Форма кромки канала струйной форсунки на входе может быть острой, закруглённой по радиусу или с фаской.
2.4. Находится относительная длина форсунки или сопла форсунки, т. е отношение
2.5. Определяется режим работы струйной жидкостной форсунки.
У струйных форсунок с острой вход-
~ ~ ^С.Н ^ 1 Г
нои кромкой при -> 1,5 возможна реалий
с.н
зация следующих режимов работы: ^ безотрывный режим истечения, когда струя компонента заполняет всё выходное сечение форсунки; ^ отрывный режим истечения, когда компонент истекает как из отверстия в тонкой стенке, не касаясь стенок сопла форсунки в выходном сечении. Смена режимов истечения происходит скачком при достижении на форсунке необходимого перепада давления Аротр;
> кавитационный режим истечения. Он имеет место у относительно длинных форсунок и является переходным от безотрывного к отрывному режиму истечения. По мере увеличения скорости истечения происходит расширение кавитаци-онной зоны и её распространение по длине форсунки к срезу. Длина зоны кавитации определяет степень развития кавитации в потоке. Условие существования ка-витационного режима истечения жидко-
сти из струйной форсунки имеет вид неравенства
к <к
кр '
(11)
где к - число кавитации, ккр - критическое число кавитации, величина которого зависит от относительной длины форсунки и числа Рейнольдса.
У коротких струйных форсунок с острой входной кромкой, когда 1,5, воз-
й
с.н
можна реализация следующих режимов работы:
> неустойчивый режим истечения в широком диапазоне чисел Рейнольдса (Яе = 103.. ,105) . Неустойчивость течения объясняется различной степенью замыкания вихревой области на стенке, которая обусловлена рядом случайных причин; ^ отрывный режим истечения.
Величина Аротр зависит от относительной длины форсунки. Очевидно, что Аротр возрастает с её повышением. Увели-
™тъАротр при рк= const и
/
= const
можно также округлением входной кромки. Чем больше радиус округления входной кромки и чем длиннее канал при постоянном диаметре сопла, тем больше Аротр. При
большом радиусе округления входной кромки отрыв потока не произойдёт даже при очень больших величинах перепада давления на форсунке, однако конструктивно такой канал удаётся выполнить далеко не всегда.
При использовании струйных форсунок в РД следует, прежде всего, учитывать их режим работы.
Неустойчивый режим течения у форсунок с острой входной кромкой, имеющих
1сн
< 1,5, является неблагоприятным для
й
с.н
двигателя из-за значительного колебания коэффициента расхода.
Кавитационный режим истечения у форсунок с острой входной кромкой при
1,5 является тоже нежелательным, т.к.
ё
с.н
вызывает неустойчивость рабочего процесса
двигателя. Поэтому и от кавитационных режимов следует отказаться.
2.6. Площадь сопла струйной жидкостной форсунки
О2)
2.7. Число Рейнольдса 4 т
Яе
\ф. Г
пцтёс
(13)
где г]Г - динамическая вязкость жидкого горючего при температуре на входе в форсунку Щфг ~ массовый расход горючего через форсунку.
2.8. Среднее значение скорости компонента на выходе из форсунки
¡V ■■
да,
1 ф. г
РгР<Ьн '
(14)
г ф.н
где рг - плотность горючего при температуре на входе в форсунку.
2.9. Коэффициент расхода форсунки ¡л : для безотрывного режима истечения
»= | 1 , , (15)
К + Я
/
й
у С.Н
где К - потери на входе, связанные с перестройкой профиля скорости, сжатием и расширением потока; Я - коэффициент линейного гидравлического сопротивления форсунки.
2.9.1. Коэффициент линейного гидравлического сопротивления форсунки Я:
^ при ламинарном течении компонента, когда Ле < 2000, определяется эмпирической зависимостью Хагена-Пуазейля
А = (16)
Яе
^ переходном и турбулентном - эмпирической зависимостью Блазиуса
Я = 0,3164-Ке0,25; (17)
^ турбулентном течении компонента, когда Яе > 10000, Я =0,031.
2.9.2. Потери на входе струйной форсунки К:
У при ламинарном течении компонента, когда Ле < 2000, определяются эмпирической зависимостью
К = 2,2-0,726-ехр(-74,5^-^), (18)
Щф.г
где цт - динамическая вязкость жидкого горючего, подаваемого через форсунку в Па-с; 1ан - длина форсунки в м; Щфт-
массовый расход горючего через форсунку в кг/с;
^ переходном и турбулентном течениях, когда Ле > 2000,
К = 1 + 2,65-Я. (19)
Для отрывного режима истечения компонента /л = 0,62...0,65. Для кавитационного режима истечения компонента
¡л = е
7т+к =6 1+
РК~Рп
ьрФ
(20)
где £ - коэффициент сжатия струи в самом узком её сечении.
Величина коэффициента сжатия струи £ в области автомодельности по числу Рейнольдса как в безотрывном, так и кавитаци-онном режимах истечения сохраняется постоянной и может быть принята равной £ =0,61...0,62.
В первой итерации перепад давления на форсунке берётся из рекомендованного интервала [1]
Дрф = 0,3... 1,5 МПа. (21)
2.10. Действительный перепад давления на струйной жидкостной форсунке горючего, обеспечивающий требуемый расход компонента топлива через форсунку:
да,
1ф.т
2 772 > ф.н
(22)
где рг - плотность жидкого горючего при температуре на входе в форсунку.
Полученный перепад давления на форсунке Арф не должен выходить за пределы
рекомендованного диапазона (21). Если это не выполняется, то необходимо изменить диаметр сопла форсунки с1ан, принятый в
первом приближении. Более точное значение диаметр сопла форсунки для следующей итерации можно получить, используя уравнение (22). В крайнем случае можно изменить число форсунок и, схему их расположения на смесительной головке или изменить диаметр камеры сгорания на входе Пк.
3. Определяется минимальная площадь сопла двухкомпонентной форсунки на выходе, обеспечивающая независимую друг от
друга работу и расчёт внутреннего и наружного контуров:
т, ,
Рф ~ Рфв + ^ф.Н ~ Рф в + щг
(23)
где тф - массовый расход жидкого горючего через наружный контур двухкомпо-нентной форсунки; рг - плотность жидкого горючего при температуре на входе в форсунку.
Если давление в камере сгорания выше критического давления горючего, то горючее переходит в газ по всему объему струи и плотность следует брать для газообразного горючего при давлении в камере сгорания и температуре на входе в форсунку. Определяется плотность по уравнению Клайперона-Менделеева:
Р,
Рг =
КТвх
(24)
Минимальный диаметр сопла двух-компонентной форсунки на выходе, обеспечивающий независимую друг от друга работу и расчёт внутреннего и наружного контуров:
¿аф = ■
4F.
ж
(25)
После определения всех проектных параметров вычерчивается конструктивная схема форсунки и выполняется её рабочий чертёж.
Библиографический список
1. Добровольский, М.В. Жидкостные ракетные двигатели. Основы проектирования: учеб. для вузов [Текст] / М.В. Добровольский; под ред. Д.А. Ягодникова,- 2-е изд., перераб. и доп. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2005.-488 с.
2. Основы теории и расчета жидкостных ракетных двигателей [Текст]/ А.П. Васильев, М.В. Кудрявцев, В.А. Кузнецов [и др.]; под ред. В.М. Кудрявцева,-М.: Высш. шк. 1993.
3. Сточек, Н.П. Гидравлика жидкостных ракетных двигателей [Текст] / Н.П. Сточек, A.C. Шапиро \!: Машиностроение, 1978128 с.
METHOD OF DESIGN CALCULATION GAS-LIQUID JET DUPLEX FUEL NOZZLE
WITH INTERNAL M IXING
©2011 V. S. Egorychev
Samara State Aerospace University named after academician S. P. Korolyov (National Research University)
Developed a method calculation of gas-liquid jet duplex fuel nozzle with internal mixing them any design concept. Expressions are obtained for determining the coefficient of discharge jet fuel nozzles on all possible modes of operation to the actual flow of the fuel components. This allows the required accuracy for practical purposes to calculate all required parameters of the projected duplex fuel nozzle.
Method, fuel nozzle, duplex fuel nozzle, jet gas-liquid fuel nozzle, liquid-fuel rocket engine, carburetion, fuel, oxidizer, propellent.
Информация об авторах
Егорычев Виталий Сергеевич, кандидат технических наук, доцент кафедры теории двигателей летательных аппаратов. Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). Тел.: (846) 334-72-00. Область научных интересов: рабочие процессы жидкостных ракетных двигателей малой тяги, термодинамика тела переменных массы и состава, неравновесная термодинамика.
Egorychev Vitaliy Sergeevich, Candidate of Technical Science, associate professor of aircraft engines theory department, Samara State Aerospace University named after academician S. P. Korolyov (National Research University). Phone: (846) 334-72-00. Area of research: small-thrust liquid-propellant rocket engines operation, variable-mass and variable-composition body thermodynamic, nonequilibrium thermodynamic.
