Методика проектирования высокоэффективных ракетных двигателей малой тяги на основе численного моделирования внутрикамерных процессов Текст научной статьи по специальности «Механика и машиностроение»

УДК 621.453.457-6

МЕТОДИКА ПРОЕКТИРОВАНИЯ ВЫСОКОЭФФЕКТИВНЫХ РАКЕТНЫХ ДВИГАТЕЛЕЙ МАЛОЙ ТЯГИ НА ОСНОВЕ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ ВНУТРИКАМЕРНЫХ ПРОЦЕССОВ

С.Д. Ваулин, В.Л. Салич

THE HIGHLY EFFECTIVE LOW THRUST ROCKET ENGINES DESIGNING METHODS, BASED ON NUMERICAL SIMULATION OF INTRACHAMBER PROCESSES

S.D. Vaulin, V.L. Salich

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

Ключевые слова: ракетные двигатели малой тяги, численное моделирование, проектирование.

The mathematical models for defining of the output characteristics and working processes parameters of the low thrust rocket engines on the components of the fuel “gas + gas”, “gas + liquid” are presented. Proposed a based on numerical simulation intrachamber processes methods of designing of the low thrust rocket engines.

Keywords: low thrust rocket engines, numerical simulation, designing.

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

Для достижения высоких показателей точности КЛА и разгонные блоки (РБ) верхних ступеней ракеты-носителя снабжаются реактивными системами управления, исполнительными устройствами в которых являются ракетные двигатели малой тяги (РДМТ), создающие необходимые управляющие усилия. Наибольшее распространение получили двухкомпонентные жидкостные ракетные двигатели малой тяги (ЖРДМТ), благодаря их высоким удельному импульсу, быстродействию и показателям надежности.

Назначение ЖРДМТ и условия эксплуатации определяют особенности организации рабочего процесса в камере ЖРДМТ, в том числе:

- малый расход топлива;

- малое число форсуночных элементов;

- трудность обеспечения равномерного распределения сравнительно небольшого количества топлива по сечению камеры;

- низкое давление в камере сгорания, особенно в момент запуска;

- отсутствие регенеративного охлаждения;

- работу в непрерывном (длительностью ~103 с) и импульсном (длительностью ~20 мс и с частотой 15-25 Гц) режимах;

- большой ресурс включений (до 106 включений).

Разработкой и созданием ЖРДМТ занято значительное число научно-исследовательских и конструкторских организаций как в России, так и за рубежом.

В настоящее время в отечественной и зарубежной ракетно-космической технике широкое распространение получили топливные композиции, горючими в которых являются производные гидразина (несимметричный диметилгидразин (НДМГ), монометилгидразин (ММГ)), окислителем - азотный тетраоксид (АТ). Данные компоненты топлива образуют самовоспламеняющиеся топливные композиции и, кроме того, являются высококипящими жидкостями. Эти свойства рассматриваются как важнейшие, поскольку обеспечивают минимальное время выхода на установившийся режим и длительное пребывание двигательной установки в состоянии готовности к запуску РДМТ.

Однако применяемые высокотоксичные топлива оказывают негативное влияние на состояние атмосферы Земли. Кроме того, продукты неполного сгорания на переходных режимах работы ЖРДМТ могут создать неблагоприятную контаминационную и токсикологическую обстановку на КЛА, а также усложняют работу экипажей в открытом космическом пространстве. Поэтому актуальным становится использование топлив, обладающих, по сравнению с топливными композициями АТ + НДМГ и АТ + ММГ, экологической чистотой и повышенными энергетическими параметрами.

В известных в настоящее время экологически чистых топливных композициях (кислород + + керосин, кислород + этиловый спирт, кислород + метан, кислород + водород) один или оба компонента являются криогенными жидкостями. Это осложняет процесс длительного хранения топлива на борту КЛА и затрудняет возможность получения высоких динамических характеристик при подаче жидкого криогенного компонента на вход в двигатель. Следовательно, система подачи должна обеспечивать газификацию криогенного компонента перед входом в РДМТ, поскольку хранение компонентов топлива на борту в газовой фазе нецелесообразно как с точки зрения обеспечения приемлемых массогабаритных характеристик аппарата, так и с точки зрения безопасности из-за высоких давлений в топливных баках. Другой важной особенностью экологически чистых компонентов является их несамовоспламеняемость, поэтому при создании РДМТ на таких компонентах необходимо дополнительно решать вопросы воспламенения топливной смеси, стабилизации горения, охлаждения конструкции агрегата зажигания и ряд других. Проблемы эти обостряются с уменьшением размерности двигателя, когда функции смесеобразования, воспламенения и охлаждения выполняются ограниченным числом (в пределе одним) смесительных элементов. В настоящее время существует единственный ЖРДМТ 17Д16 тягой ~200 Н (разработка НИИМаш), созданный для объединенной двигательной установки многоразового орбитального корабля «Буран» и прошедший полный цикл отработки. При его создании пришлось решать совершенно новые проблемы по смесеобразованию, взрывобезопасности, стабильному воспламенению и обеспечению длительного ресурса работы [1-2].

Важнейшим параметром РДМТ является удельный импульс, поскольку увеличение удельного импульса снижает расход топлива, требуемый на создание заданного управляющего усилия. Например, увеличение удельного импульса для двигателя тягой 100 Н с 2500 до 2900 м/с при времени огневой работы 50 000 секунд ведет к снижению 280 килограммов потребной массы топлива или к соответствующему увеличению массы полезной нагрузки. При использовании 12 таких двигателей эффект выигрыша в массе может составлять уже около 3,3 тонны.

Удельный импульс определяется как

I уд = Р‘^ ,

где в - расходный комплекс камеры РДМТ; Kп - коэффициент тяги камеры, определяемый степенью расширения сопла.

Расходный комплекс определяется как

в = Фррт ,

где вт - теоретический расходный комплекс, зависящий от конкретного вида топлива; фр - коэффициент расходного комплекса, показывающий полноту преобразования компонентов топлива в камере двигателя в продукты сгорания. Следовательно, достижение высоких значений фр (фр^ 1)

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

Обеспечение высококачественного смесеобразования в камерах РДМТ при характерных для них ограниченном числе смесительных элементов и малых размерах конструкции является одной

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

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

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

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

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

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

2. Горение осуществляется в газовой фазе.

3. Температура испаряющейся капли одинакова во всем объеме капли.

4. Газовая среда включает в себя три вещества: «горючее» (в газовой фазе), «окислитель», «продукты сгорания». Все химические реакции сведены к одной необратимой брутто-реакции, описывающей стехиометрическое взаимодействие веществ «окислитель» и «горючее» с появлением вещества «продукты сгорания».

5. Газовая смесь находится в локальном термодинамическом равновесии.

Комментируя данные допущения, можно отметить, что модель горения капель жидкого горючего в газообразной окислительной среде, в которой выделяются стадии прогрева капель, испарения и горения паров, широко применяется для подобных задач и дает удовлетворительные результаты [3-4]. Модель брутто-реакции является удобной, так как для сложной химической реакции теоретическая оценка суммарной скорости процесса затруднена вследствие недостатка точных значений кинетического механизма и констант скорости отдельных реакций. Для учета влияния на процессы реакций диссоциации и рекомбинации принята модель локального термодинамического равновесия рабочего тела, которая позволяет представлять состав как функцию состояния. В этом случае не имеет значения последовательность химических реакций при переходе к равновесию. Особенностью предложенной модели является то, что влияние на процессы реакций диссоциации и рекомбинации учтено через замыкающие зависимости от температуры и давления молекулярных масс и «химически равновесных» теплоемкостей веществ, составляющих газовую среду. Для этого предварительно осуществляется серия термодинамических расчетов равновесного состава веществ «горючее», «окислитель», «продукты сгорания» при различных давлениях и температурах.

Применяя для описания процессов в камере сгорания РДМТ на двухфазных компонентах топлива дискретно-траекторный метод Эйлера - Лагранжа [5] и основываясь на существующем заделе в области моделирования процессов в камерах сгорания тепловых машин [3-6], получим систему уравнений сохранения, соответствующую принятым допущениям:

д(аР)

дґ

д(аР£;)

У-(г^) = П; (1)

дґ + ^ ) = ^; -^.; (2)

+ V• (у ® у - П + = ) у; (3)

д(аРЕ) дґ

V- (руЕ - Пу + q ))£; (4)

<^г = Р=ж; (5)

М

£ )=; (6)

т -_ и.; (7)

М

ч -_^Т + £3]Н1 . (8)

I-1

В данных уравнениях: V - вектор скорости газовой смеси; р - плотность газовой смеси; р -давление; га - доля объема, занимаемая газовой фазой; gj - массовая доляу-го вещества в газовой смеси (у = г, о, пс, где индекс г означает испарившееся горючее, о - окислитель, пс - продукты сгорания); П - тензор напряжений поверхностных сил [3]; ш1и сп - мгновенный массовый секундный расход паров от 7-й капли; т‘ж - масса 7-й капли; Е = суТ - внутренняя энергия газовой фазы;

Т - температура газовой смеси; еу - изохорная теплоемкость газовой смеси; П, ¥„, ¥Е - источники соответственно массы, импульса, энергии, обусловленные взаимодействием газовой фазы с каплями распыленного топлива; Q7 - количество теплоты, передаваемое к капле от газовой фазы в единицу времени:

& - ^ ^ + & , (9) где составляющие Qc и QR обусловлены соответственно конвективным и лучистым переносом; составляющая Qm обусловлена уносом тепла от капли за счет фазового перехода

N

^ --£ Чсп-г, (10)

7 -1

где г - скрытая теплота фазового перехода горючего.

Диффузионный потоку-го компонента газовой смеси определяется как

3 —р О у Vgj, (11)

где Бу - эффективный коэффициент диффузии,у = г, о, пс.

Вектор результирующей силы, действующей на 7-ю каплю, определяется как

Рж - РБ + РМ + РУМ + Рр , (12)

где Рв - сила аэродинамического сопротивления; РМ - сила Магнуса; РУМ - виртуальная массовая сила, обусловленная увлечением массы газовой фазы при ускорении капли; Рр - сила вследствие градиента давления. Значения Рв, РМ, РУМ, Рр определяются по зависимостям, приведенным в [6]. Источники в уравнениях (1), (3), (4) определяются следующим образом:

N

п-£Nimисп ; (13)

7 -1

* V = IN РЖ + X N,171 и сп У, ;

і=1

і=1

* я = I N ( и сп к — а), і=1

(14)

(15)

где Н - энтальпия горючего в газовой фазе; N - число капель группы 7 в единице объема. Источники в уравнениях типа (2):

£г = —Мг Уг Лг + П;

(16)

(17)

(18)

=-М0У0Я;

^пс = Мпс ^пс Яг ,

где Мг, Мо, Мпс - соответственно молекулярные массы испарившегося горючего, окислителя, продуктов сгорания; уг, Уо, Упс - стехиометрические коэффициенты соответственно горючего, окислителя и продуктов сгорания в уравнении брутто-реакции; Яг - скорость горения, определяемая как меньшее значение из скорости перемешивания за счет турбулентных пульсаций (модель «диссипирующего вихря») и скорости химической реакции [6]:

Яг = тіп

АВ- тіп

к

Є У В- пс

к V

(19)

-У

пс

^о

К ехр ((ЯоГ)).

В последнем выражении Уг, Уо, Гпс - соответственно мольные концентрации горючего, окислителя и продуктов сгорания; Фг, Фо - эмпирические константы; К - предэкспоненциальный множитель в законе Аррениуса; к - кинетическая энергия турбулентных пульсаций; е - скорость диссипации турбулентной энергии. Значения К, Фг, Фо для конкретной брутто-реакции берутся из [7].

Приведенная система дополняется уравнениями, соответствующими принятой модели турбулентности. Учет вторичного дробления капель осуществляется согласно [6]. Расчетная область ограничивается внутренним объемом камеры и небольшим участком расширяющейся части сопла. В качестве граничных условий для газовой фазы на входных границах задаются значения массовых расходов, температур и интенсивности турбулентной пульсации, на выходной границе -значение давления. Для дискретной фазы на входной границе задаются суммарный расход капель, температура, функция распределения капель по размерам.

Полученная модель описывает процессы в РДМТ на двухфазных компонентах топлива (например, О2аз + С2Н5ОН). Модель процессов в камере РДМТ на газообразных компонентах топлива (например, 02^ + Н2^) является частным случаем представленной модели.

Полученные модели реализованы в пакете ЛК8У8 СБХ-11. На вычислительном кластере Южно-Уральского государственного университета были проведены расчеты действующих камер РДМТ. По результатам расчета определялся коэффициент расходного комплекса. При этом сравнивались значения, полученные двумя способами обработки результатов расчета:

1) способом, основанным на определении терминов «коэффициент расходного комплекса», «теоретическое (идеальное) значение расходного комплекса», «действительное значение расходного комплекса» [8]:

р М(ТГ)

о,/Яа

где рк - давление в штуцере замера, определенное по результатам расчета; ^кр - площадь критического сечения камеры сгорания; О - суммарный массовый секундный расход горючего и окислителя; Як, Тк - газовая постоянная и температура продуктов сгорания в камере соответст-

венно, найденные из термодинамического расчета равновесного состава; А(у) - комплекс, зависящий от показателя процесса расширения продуктов сгорания до давления в критическом сечении сопла;

2) способом, основанным на формуле, представленной в работе В.И. Лапицкого [9]:

Фр

J/~гвых

1 - Яг

G

(21)

где Огвх, Огвых - массовые расходы горючего на входе в расчетную область и выходе из нее соответственно.

Методика проектирования РДМТ

Вычисленные обоими способами значения фр для различных камер РДМТ на компонентах

топлива О^ + СН^33, 02°+ Н^33, О2аз+ С2Н5 ОН удовлетворительно согласуются с экспериментальными данными. Результаты расчета и сравнение их с экспериментальными данными приведены в работах [10, 11]. Достижение высокого уровня сходимости результатов модельных расчетов с экспериментальными данными позволило предложить методику проектирования высокоэффективных РДМТ на основе результатов численного моделирования внутрикамерных процессов, которая наглядно представлена на рисунке.

Использование данной методики позволит принимать научно обоснованные технические решения при проектировании новых РДМТ, экономить финансовые средства на разработку, сокращать сроки разработки новых проектов, значительно повышать производительность труда проектировщика.

Литература

1. Салич, В.Л. Жидкостные ракетные двигатели малой тяги: учеб. пособие / В.Л. Салич, А.А. Шмаков, С.Д. Ваулин. - Челябинск: Изд-во ЮУрГУ, 2006. - 52 с.

2. Кутуев, Р.Х. Разработка перспективных РДМТ на экологически чистых топливных композициях /Р.Х. Кутуев, И.Н. Лебедев, В.Л. Салич // Вестн. Самар. гос. аэрокосм. ун-та им. акад. С.П. Королева. - 2009. - № 3 (19), ч. 3. - С. 101-108.

3. Варнатц, Ю. Горение. Физические и химические аспекты, моделирование, эксперименты, образование загрязняющих веществ / Ю. Варнатц, У. Маас, Р. Диббл; пер. с англ. Г.Л. Агафонова; под ред. П. А. Власова. - М. : Физматлит, 2006. - 352 с.

4. Рабочие процессы в жидкостном ракетном двигателе и их моделирование / Е.В. Лебединский, Г.П. Калмыков, С.В. Мосолов и др.; под ред. акад. РАН А.С. Коротеева. - М. : Машиностроение, 2008. - 512 с.

5. Юн, А. А. Расчет и моделирование турбулентных течений с теплообменом, смешением, химическими реакциями и двухфазных течений в программном комплексе FASTEST-3D: учеб. пособие /А.А. Юн, Б.А. Крылов. - М. : Изд-во МАИ, 2007. - 116 с.

6. ANSYS CFX-Solver, Release 10.0: Theory. - ANSYS Europe Ltd., 2005. - 266p.

7. Физико-химические процессы в газовой динамике: справ. В 2 т. Т. 2: Физико-химическая кинетика и термодинамика / под ред. акад. Г.Г. Черного, С.А. Лосева. - М. : Науч.-издат. центр механики, 2002. - 368 с.

8. Дорофеев, А.А. Основы теории тепловых ракетных двигателей (Общая теория ракетных двигателей): учеб. для авиа- и ракетостроит. специальностей вузов / А.А. Дорофеев. - М.: Изд-во МГТУ им. Н.Э. Баумана,1999. - 415 с.

9. Лапицкий, В.И. Математическое моделирование и экспериментальное исследование

характеристик камеры сгорания ракетного двигателя малой тяги на метане и кислороде: автореф. дис. ... канд. техн. наук / В.И. Лапицкий. - М.: МГТУ им. Н.Э. Баумана, 2006. -

16 с.

10. Салич, В.Л. Результаты работ по созданию высокоэффективного ракетного двигателя малой тяги с компонентами топлива кислород - метан, кислород - водород / В.Л. Салич // Тр. VIIIВсерос. с междунар. участием науч.-техн. конф. «Авиакосмические технологии (АКТ-2007)». -Воронеж: ВГТУ, 2007. - С. 233-239.

11. Ваулин, С.Д. Исследование энергоэффективности в ракетных двигателях малой тяги на двухфазных компонентах топлива / С.Д. Ваулин, В.Л. Салич, В.И. Феофилактов // Вестник ЮУрГУ. Серия «Энергетика». - 2011. - № 34 (251). - С. 81-85.

Поступила в редакцию 2 февраля 2012 г.

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

«Двигатели летательных аппаратов», проректор по научной работе, Южно-Уральский государственный университет. Область научных интересов - исследование теплофизических процессов в ракетных

двигателях и энергетических установках. Тел.: (351) 263-45-97; е-mail: s.d.vaulin@susu.ac.ru

Sergey D. Vaulin. The doctor of engineering science, professor, head of a chair of engine for flying vehicle, vice rector of scientific work, South Ural state university. The area of scientific interests - research the thermalphysic processes of rocket engines and powerplants. Tel.: (351) 263-45-97; e-mail: s.d.vaulin@susu.ac.ru

Салич Василий Леонидович. Кандидат технических наук, докторант кафедры «Двигатели летательных аппаратов», Южно-Уральский государственный университет. Область научных интересов -исследование теплофизических процессов в ракетных двигателях и энергетических установках. Тел.: (34345) 3-62-46; е-mail: salich_vls@mail.ru

Vasily L. Salich. The candidate of engineering science, doctoral candidate on the chair of engine for flying vehicle, South Ural state university. The area of scientific interests - research the thermalphysic processes of rocket engines and powerplants. Tel.: (34345) 3-62-46; е-mail: salich_vls@mail.ru