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

УДК 536.46

Д. А. Ягодников, Ю. В. Антонов, Ю. Н. Власов

МОДЕЛИРОВАНИЕ ИСПАРЕНИЯ ПОЛИДИСПЕРСНОЙ СОВОКУПНОСТИ КАПЕЛЬ ВОДЫ В КАМЕРЕ СГОРАНИЯ ГИДРОРЕАКТИВНОГО ДВИГАТЕЛЯ

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

E-mail: daj@mx.bmstu.ru; prog1185@mail.ru Ключевые слова: испарение, распиливание, вода, полидисперсная совокупность, математическая модель, капля, численное исследование.

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

Рассматриваемый в настоящей работе гидрореактивный двигатель (ГРД) предназначен для передвижения аппаратов в водной среде, поступающей в него в качестве окислителя, и представляет собой разновидность прямоточной силовой установки, схема которой приведена на рис. 1. Физику процессов, протекающих во внутреннем объеме такого двигателя, можно представить следующим образом. Продукты разложения пиротехнического состава (ПС), отходящие от поверхности горения, вступают во взаимодействие с водой, подаваемой в первую зону. В промежуточном сопле поток ускоряется, что способствует дроблению капель воды и интенсифицирует процессы тепломассообмена. Дожигание продуктов сгорания осуществляется во второй зоне

Рис. 1. Схема и расчетная область ГРД:

1 — поверхность горения пиротехнического состава; 2 — первая зона (камера сгорания); 3 — подача воды; 4 — вторая зона (камеры дожигания); 5 — тяговое сопло; 6 — промежуточное сопло

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

В настоящее время физика рабочих процессов в КС ГРД проработана на достаточно высоком уровне, чего нельзя сказать о количественной теории. Наиболее показательной в этом представлении является работа [5]. Авторы приводят в одномерной параметрической постановке теоретическое описание основных аспектов рабочего процесса, характерных для КС ГРД, обходя стороной ряд его принципиальных особенностей, что не позволяет говорить об актуальности модели на данном этапе развития теории многофазных турбулентных течений. В частности, на стадии выбора геометрических параметров зон горения и тепломассообмена камеры двигателя в целях организации высокоэффективного рабочего процесса с высокой полнотой сгорания необходимо знать длину пути испарения жидкого компонента топлива, которая будет существенно различной в зависимости от начального значения диаметра капли. В связи с этим большое практическое значение имеет корректное качественное и количественное описание сопутствующих процессов в общей математической модели.

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

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

Основные этапы подачи забортного окислителя в КС ГРД, протекающие зачастую одновременно, можно условно подразделить на

распыливание, движение и испарение совокупности капель разных диаметров. Существует множество различных способов распылива-ния [6], среди которых следует выделить гидравлический, как наиболее часто встречающийся в принципиальных схемах смесеобразования ракетных и реактивных двигательных установок. Основным фактором, приводящим к диспергированию жидкости, является давление нагнетания. Физику процесса распыливания можно представить следующим образом. Проходя через форсунку поток ускоряется и преобразуется в форму, способствующую быстрому и эффективному распаду. В общем случае распад жидкости вызывают несколько факторов, а именно: капиллярные и внешние инерционные силы, турбулентные пульсации, кавитация. В зависимости от способа диспергирования могут превалировать один, либо несколько из перечисленных возмущений, но главной причиной разрушения струи является развитие в ее объеме колебательных процессов. Возникновение последних обусловлено действиями аэродинамических сил, деформирующих и разрывающих струю, и вибрациями, являющимися следствием качества изготовления форсунки и ее конструктивных особенностей. Начальные возмущения появляются одновременно и могут накладываться друг на друга, поэтому поверхностные колебания в зависимости от наложения фаз начальных возмущений и степени рассеяния энергии колебаний могут затухать, возрастать либо стабилизироваться [7, 8]. Рост и затухание колебаний обусловливаются физическими свойствами окружающей среды и режимом вытекания жидкости из соплового отверстия. В зависимости от начальной скорости потока возможно образование единичных капель, коротких струй и неравномерных капель, капель с приблизительно одинаковыми размерами, полидисперсной смеси капель. Причем на практике последнее является наиболее вероятным.

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

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

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

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

Следует отметить, что при рассмотрении ансамбля образующихся капель как монодисперсных заранее вносится существенная погрешность в результаты расчетов и если на ранних этапах развития теории численного моделирования движения многофазных сред подобное допущение встречалось практически во всех моделях, то в настоящее время значительный прогресс вычислительной техники позволяет избавиться от последнего путем рассмотрения движения конечного числа фракций полидисперсной совокупности капель. Говоря о моделировании распределения капель по размерам, необходимо отметить значительное число эмпирических уравнений, среди которых следует выделить логарифмически нормальное распределение [9], уравнение числовой частотной кривой, предложенное Нукиямой и Танасавой [10], и уравнение Розина-Раммлера [11, 12]. Однако выбор типа исходной кривой распределения капель по размерам — это прерогатива исследователя, и вид зависимости не является определяющим фактором при разработке комплексной методики расчета параметров рабочего процесса в объеме КС ГРД, где принципиально важны сам факт учета полидисперсности и число рассматриваемых фракций, которое прямо влияет на затраты машинного времени и точность моделирования.

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

(1)

(2)

где У — суммарный относительный объем, т.е. сумма отдельных объемов всех капель, диаметр которых меньше или равен диаметру данной капли йж;

Т йЖ.м /тч

Т = --ГТ/п (3)

1п2

— константа размера (йж м — медианный диаметр капель); п — константа распределения, определяемая по соотношению

(1п _ 1 \

1„ 1 'в тах

Ч 1п 2

п = / й-И- • (4)

Ч 4

Приводя функцию У3 для конечного интервала размеров условимся, что при йж = 4Ж. тах суммарный относительный объем будет равен

V = Утах = 0,95.

В настоящее время имеется широкий спектр зависимостей для оценки 4жм в зависимости от условий окружающей среды и способа распыливания [9]. В общем виде это эмпирическое соотношение записывается следующим образом:

4ж.м = /(Яет, Ьр1, Weк, М,*>), (5)

где Яе, Ьр, We, М — числа Рейнольдса, Лапласа, Вебера и Маха соответственно; ^ — площадь впускного канала форсунки; V — скорость жидкости на входе.

Для определения медианного диаметра капли жидкости при условиях, реализуемых в камере ГРД, воспользуемся зависимостью (5), справедливой в диапазоне параметров в = 30°... 90°, 4С = 1,0 ... 2,0 мм, Дрф = 0,2 ... 1,0 МПа, приведенной в работе [5]:

4ж.м = 42,5а0'23Др-0'23й-0'43, (6)

где в — угол наклона струи к плоскости торца заряда; 4С — диаметр сопла форсунки.

Максимальный диаметр капли жидкости вычисляем, исходя из рекомендаций работ [7, 9, 11], по формуле

йжтах ~ 2йж.м.

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

распределения капель по размерам. Типичный вид кривой приведен на рис. 2.

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

Таблица 1

Фракционный состав капель воды

Номер фракции, г йж, мкм Объемная доля капель г-й фракции в потоке жидкости ^ж/У^ж (при У^ = К = 0,95)

1 10,5 0,005

2 32,2 0,045

3 54 0,1

4 76,2 0,13

5 98 0,15

6 120 0,14

7 140,5 0,12

8 163 0,11

9 184,7 0,08

10 207,3 0,07

Рис. 2. Суммарное распределение капель воды по размерам

При описании движения капель воды используются координаты Лагранжа. Для принятого распределения среднесчетный (й10), средне-поверхностный (й32) и среднемассовый (й43) диаметры соответственно равны 162,85, 76,92 и 114,82 мкм. Система уравнений, описывающих движение и испарение капли ъ-й фракции в двумерной постановке, включает в себя следующие соотношения:

Рк

Рпг =

Rnr Тпг

— плотность продуктов сгорания (парогаза) — из уравнения состояния;

¿иж 3 рпг I I / \

—7Т = 7—ТСх|и - иж|(и - иж) М 4 Рж йж

— уравнение продольной скорости движения (рж — плотность жидкости; йж — текущий диаметр частицы рассматриваемого класса;

Сх= — коэффициент аэродинамического сопротивления кап-

уЯеж

ли);

dvж 3 рпг ^-у I -» -» I / \

-Сж|м - Мж| (V - Vж)

4 ржйж х

уравнение радиальной скорости движения;

ж

= «я

— осевая составляющая траектории;

йГж

■¿Т = ^

— радиальная составляющая траектории;

¿Мж _ _

~1Г =

— изменение массы капли Мж;

^0 = та^{Соисп ,^0об}

— массовый поток паров воды у поверхности капли; здесь

исп

— поток за счет диффузии и конвекции,

Сооб = Пйж|и — "ЙУУРпгРж 6

— поток за счет обтекания капли газовым потоком [2],

^пРпг , -1 - с^ \ БЪ

^ ^^ Т

— плотность потока массы у поверхности испаряющейся капли, где с0 — концентрация паров воды по пограничному слою [13];

*=°,73< ж)1,8210

— среднее по пограничному слою капли значение коэффициента диффузии (Тжср — средняя температура в пограничном слое капли);

БЪ = 2(1 + 0,3^Кёж8с 3

ттт .т, ¿жРпг|и - «ж1 „

— число Шервуда (Кеж =--относительное число Рей-

Рп

нольдса; Бс = —^ — число Шмидта; рп — среднее по пограничному слою значение коэффициента динамической вязкости).

Изменение температуры капли определяется из решения уравнения баланса энтальпии

~ГГ = ~Т~П-(а(Тпг — Тж) — g0испQи),

ас йжсжрж

где и Сж — скрытая теплота испарения воды и удельная теплоемкость; а — коэффициент теплоотдачи.

При определении коэффициента теплоотдачи воспользуемся методикой Эл Вакила [14], в соответствии с которой вначале вычисляется значение в отсутствие массопереноса, затем последнее умножается на корректирующий коэффициент:

Апг № Zk

а =--■

а ¿ж е^ - 1 ■

-с г Z = СрЧ ас ).

Zk = -

п^жАпг Nu

Ми = 2 + 0,54^^ Рг3.

Расчет процессов дробления выполнялся с помощью модели Пилча и Эрдмана [15].

Система приведенных дифференциальных уравнений решалась методом Рунге-Кутта 4-го порядка.

В соответствии с приведенным алгоритмом на языке С++ была написана программа расчета характеристик испарения капель воды в условиях КС ГРД. Постановку начальных и граничных условий осуществляли, исходя из принципа минимальных видоизменений программы при моделировании различных условий эксперимента. В частности, полагали, что вектор скорости впрыска капли параллелен оси камеры двигателя, т.е. на границе подачи жидкости (поз. 3, см. рис. 1) имеем и = иж0, V = vж0 = 0, Тж = Тж0 = 293К; температура парогаза остается постоянной на входе в камеру двигателя Тпг0 = 3000 К; КС — изобарическая, с давлением рк = 9 МПа. Перепады давлений на форсуночной головке и радиус сопла форсунки составляют Дрф = 1 МПа и гс = 1 мм соответственно. Основные тепло-физические характеристики парогазовой смеси примем следующими: Лпг = 678 Дж/(кг-К), с^ = 0,27, Српг = 2310 Дж/(кг-К), = 39,16 х х 10-7 м2/с, Апг = 0,57Вт/(м-К).

Результаты расчета процесса движения и испарения двух видов совокупности капель: полидисперсной, состоящей из десяти фракций, и монодисперсной с характерным размером капли, равным среднемас-совому диаметру, приведены в табл. 2 и 3. Сопоставление упомянутых результатов с данными табл. 1 и рис. 2 позволяет сделать вывод о существенном влиянии реального распределения капель по размерам

Таблица 2

Результаты расчета полидисперсной смеси

Номер фракции, г ¿ж, мкм Время испарения капель г-й фракции, мс Длина пути испарения капель г-й фракции, мм

1 10,5 0,0486 1,4

2 32,2 0,3230 8,5

3 54 0,7621 19,2

4 76,2 1,3168 32,8

5 98 1,9799 48,5

6 120 2,7534 66,4

7 140,5 3,5371 84,7

8 163 4,4827 106,5

9 184,7 5,5254 129,2

10 207,3 Результат 6,5849 ы расчета монодисперс 154,2 Таблица 3 ной смеси

Средний диаметр, мкм Время испарения капель г-й фракции, мс Длина пути испарения капель г-й фракции, мм

dio 162,85 4,4990 106,4

^32 76,92 1,3414 33,3

d43 114,82 2,5583 62,0

на рабочий процесс, протекающий во внутреннем объеме КС ГРД, посредством времени и характерной длины испарения. При этом исключение фракций малого диаметра под номерами г = 1 и г = 2 из состава совокупности не окажет серьезного влияния на общую динамику процесса из-за их малого числа (У1ж/У^ж = 0,005; У2ж/УЕж = 0,045) и интенсивного испарения. Однако пренебрежение шестой (г = 6) и седьмой (г = 7) фракциями снижает на 27% точность оценки времени испарения и на 25-26 % точность оценки характерной длины пути испарения.

Верификация результатов расчетов применительно к условиям эксперимента [16], приведенная на рис. 3, свидетельствует об их удовлетворительном согласовании. Некоторое отличие расчетных

d/do

0 50 100 150 t

Рис.3. Сравнение результатов расчетов с данными работы [16]:

1 — эксперимент для воды при Т = 78 0С, V = 5,1 м/с; 2 — эксперимент для воды при Т = 80,50 С, V = 5,14 м/с(...... — численный расчет по авторской методике;--численный расчет монодисперсной смеси с характерным размером капли ¿ю)

результатов связано, вероятнее всего, с использованием приближенных полиномиальных зависимостей давления насыщенных паров, теплоемкости воды и удельной теплоты испарения от температуры, составленных для условий, отличных от условий эксперимента [16].

Приведенная в настоящей работе методика моделирования процессов движения и испарения капель воды является частью комплексной программы расчета основных характеристик рабочего процесса в объеме КС ГРД. На рис. 4 и 5 приведены примеры результатов расчета полей скорости и температуры. Анализ распределения вектора скорости несущего газа в объеме КС ГРД (см. рис. 4) позволяет сделать вывод о наличии интенсивного вихря во второй зоне, присутствие которого в целом оказывает положительное влияние на рабочий процесс, увеличивая время пребывания капель воды в камере двигателя, а значит и степень их испарения. Рассматривая температурное поле (см. рис. 5), следует отметить наличие низкотемпературных зон вблизи форсуночной головки первой и второй зон, что является благоприятным с точки зрения организации эффективной тепловой защиты элементов конструкции. Размеры упомянутых зон, как и масштаб вихря существенно изменяются (вплоть до полного исчезновения) при изменении режимных параметров рабочего процесса (давление в КС, скорость подачи компонентов топлива).

Резюмируя результаты выполненных исследований, можно заключить, что разработанная математическая модель расчета процес-

5.аВ5в+С02

N 3 92Эе+002

1.962В+002

ш О ОООе+ОИ [т £'-1]

Рис. 4. Распределение вектора скорости несущего газа

—г 3.003е+003 2.320е+00 1.б37е+оо 9.543е+00

2.713е+002

[К]

Рис. 5. Распределение температуры по объему КС ГРД

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

СПИСОК ЛИТЕРАТУРЫ

1. Р а к е т н о-прямоточные двигатели на твердых и пастообразных топливах. Основы проектирования и экспериментальной отработки / Л.С. Яновский и др. - М.: Физматлит, 2010. - 320 с.

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

3. Bjorn Hussmann, Michael Pfitzner. Extended combustion model for single boron particles - Part I: Theory // Combustion and Flame. - 2010. - Vol. 157. -P. 803-821.

4. Bjorn Hussmann, Michael Pfitzner. Extended combustion model for single boron particles - Part II: Validation // Combustion and Flame. - 2010. -Vol. 157. - P. 822-833.

5. Реактивные двигатели подводных аппаратов на твердом топливе / Е.С. Ша-хиджанов, А.Ф. Мяндин. - М: Изд-во ГНПП "Регион", 2005. - 232 с.

6. Основы техники распыливания жидкостей / Д.Г. Пажи, В.С. Галустов. - М.: Химия, 1984.-256 с.

7. Процессы распыливания топлива дизельными форсунками / А.С. Лышев-ский. - М.: Машгиз, 1963. - 181 с.

8. Вивденко М. И., Шабалин К. Н. Исследование условий получения равномерных капель размером 1-0,5 мм // Изв. вузов. Химия и хим. технолог. -1965.-Т. 8, №4.-С. 685-690.

9. Распыливание жидкостей / Ю.Ф. Дитякин, Л.А. Клячко, Б.В. Новиков, Б.И. Ягодкин. -М.: Машиностроение, 1977.

10. N u k i j a m a and Tanasawa J. Transactions of the Mechanical Engineers. -1939. - Vol. 5. No 18.

11. Физические основы рабочего процесса в камерах сгорания воздушно-реактивных двигателей / Б.В. Раушенбах. - М.: Машиностроение, 1964. - 526 с.

12. А н т о н ю к П. Н. Распределение по размерам капель распыливаемой жидкости // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. - 2007. - Спец. вып. - С. 46-50.

13. Щербаков А. А. Расчет процессов распыла, движения и испарения жидкого компонента топлива // Методические указания по выполнению домашнего задания по курсу "Рабочие процессы". - М.: Изд-во МГТУ им. Н.Э. Баумана, 1984. - 16 с.

14. El Wakil M. M., Ueyhara O. A., and Myers P. S. 1954.'A theoretical investigation of the heating-up period of injected fuel droplets vaporizing in air', NASA Technical Note 3179.

15. Pilch M., and E r d m a n C. A. 1987. Use of breakup time data and velocity history data to predict the maximum size of stable fragments for acceleration-induced breakup of a liquid drop // Int. J. Multiphase Flow, 13(6). - P. 741-757.

16. Терехов В. И., Терехов В. В., Шишкин Н. Е., Б и К. Ч. Экспериментальное и численное исследования нестационарного испарения капель жидкости // Инж.-физич. журнал. - 2010. - Т. 83, № 5. - С. 829-836.

Статья поступила в редакцию 11.10.2011

Дмитрий Алексеевич Ягодников родился в 1961 г. Д-р техн. наук, зав. кафедрой "Ракетные двигатели" МГТУ им. Н.Э. Баумана. Автор более 160 научных работ в области экспериментально-теоретического исследования процессов воспламенения и горения порошкообразных металлов и бесконтактной диагностики ракетных двигателей.

D.A. Yagodnikov (b. 1961) — D. Sc. (Eng.), head of "Rocket Engines" department of the Bauman Moscow State Technical University. Author of more than 160 publications in the field of experimental and theoretical study of processes of ignition and combustion of powder-like metals and contactless diagnostics of rocket engines.

Юрий Владимирович Антонов родился в 1985 г. Окончил МГТУ им. Н.Э. Баумана в 2009 г. Аспирант, ассистент кафедры "Ракетные двигатели" МГТУ им. Н.Э. Баумана. Автор 12 научных работ в области рабочих процессов в тепловых энергетических установках, конверсионного использования ракетных двигателей.

Yu.V. Antonov (b. 1985) graduated from the Bauman Moscow State Technical University in 2009. Post-graduate, assistant lecturer of "Rocket Engines" department of the Bauman Moscow State Technical University. Author of 12 publications in the field of working processes in thermal power plants, conversion utilization of rocket engines.

Юрий Николаевич Власов родился в 1937 г. Окончил МВТУ им. Н.Э. Баумана в 1965 г. Канд. техн. наук, доцент кафедры "Ракетные двигатели" МГТУ им. Н.Э. Баумана. Автор более 100 научных работ в области рабочих процессов в тепловых энергетических установках, конверсионного использования ракетных двигателей.

Yu.N. Vlasov (b. 1937) graduated from the Bauman Moscow Higher Technical School in 1965. Ph. D. (Eng.), assoc. professor of "Rocket Engines" department of the Bauman Moscow State Technical University. Author of more than 100 publications in the field of working processes in thermal power plants, conversion utilization of rocket engines.