Уфа : УГАТУ, 2010
'Be&тн,и,к,
Т. 14, № 5 (40). С. 15-26
АВИАЦИОННАЯ И РАКЕТНО-КОСМИЧЕСКАЯ ТЕХНИКА
УДК 629.735
А. Б. Козловская, Д. А. Ахмедзянов, Ю. М. Ахметов, А. Е. Михайлов
РАСЧЕТ ХАРАКТЕРИСТИК ОСЕВЫХ КОМПРЕССОРОВ СОВМЕСТНО С ГРАНИЦЕЙ УСТОЙЧИВОЙ РАБОТЫ
В программном комплексе COMPRESSOR реализована возможность расчета характеристик многоступенчатых компрессоров с помощью метода сложения характеристик отдельных ступеней. С помощью ограничения по критическому углу атаки и метода Хауэлла реализована возможность построения на характеристике границы устойчивой работы компрессора. Проведено сравнение расчетной и экспериментальной характеристик компрессора. Авиационные двигатели; осевые компрессоры; расчет характеристик; граница устойчивой работы компрессора; имитационное моделирование
ВВЕДЕНИЕ
Проектирование и создание авиационных ГТД является одной из наиболее наукоемких и высокоразвитых отраслей промышленности. Развитие авиационных ГТД базируется на комплексе фундаментальных и прикладных исследований в аэро- и термодинамике, материаловедении, механике жидкости и газа. Авиационное двигателестроение, основанное на наиболее передовых технологиях, стимулирует развитие всех тех отраслей промышленности, где требуются компактные, мобильные и хорошо управляемые силовые установки, - наземный и водный транспорт, теплоэнергетика, транспортировка газа, технологии сушки, очистки, пожаротушения и т. д.
В настоящее время созданы и существуют четыре поколения авиационных двигателей и ведется работа над созданием двигателя пятого поколения. Двигатели нового поколения для боевой авиации должны обладать рядом характеристик, придающих качественно новый уровень самолету. К ним относятся уменьшение удельного веса двигателя в 1,4-2 раза, удельного расхода топлива - на 15-30%, повышение надежности на 60-80%, обеспечение ресурса двигателя, соответствующего 0,5-1 ресурса планера, снижение в 2-3 раза трудоемкости обслуживания.
Среди основных направлений развития современных компрессоров следует отметить увеличение напорности и КПД ступеней, что позволяет сократить их количество. Ведутся работы по созданию малошумного вентилятора бла-
Контактная информация: 8(347)273-79-54 Работа выполнена при финансовой поддержке РФФИ, поддержана грантом Президента Российской Федерации НШ-65497.2010.9 для ведущих научных школ
годаря введению стреловидных широкохордных лопаток без антивибрационных полок, применению технологии «блиск», широкому использованию композиционных материалов. Для компрессоров высокого давления характерно уменьшение числа ступеней с одновременным повышением общей степени повышения давления и КПД. Проявляется большой интерес к центробежным компрессорам. Возможность достижения в них высоких степеней повышения давления при небольшой осевой протяженности и массе компрессора и двигателя отодвинула на второй план такой недостаток центробежного компрессора, как более низкий КПД, по сравнению с осевым компрессором.
Широкий спектр проектных работ, сложность конструкции авиационных двигателей, сжатость сроков выполнения заказов, усиление конкурентной борьбы на рынке заставляют авиамоторные предприятия проводить работы на самом современном техническом уровне. Конструирование, анализ и технологическая подготовка производства осуществляются с использованием преимущественно «тяжелых
САПР». Выбор, организация работы и их использование производятся в рамках единой концепции развития средств проектирования и технологической подготовки производства, способных решать поставленные перед разработчиками задачи.
На современном этапе развития газотурбинной техники становится невозможным проводить проектирование и расчеты перспективных изделий без применения имитационного и 3Б-САБ/САБ-моделирования. Как и всякий процесс проектирования сложного узла, проектирование компрессоров (а совершенство компрессора во много и определяет совершенство двигателя в целом) авиационных ГТД носит итерационный характер, в ходе которого от предварительных схем и упрощенных расчетов идет
движение к реальному представлению об объекте с углублением проработки его основных характеристик и исследованием второстепенных особенностей. В данной работе рассмотрена задача получения характеристик осевых компрессоров, которая является одной из сложнейших задач при их проектировании.
СОСТОЯНИЕ ВОПРОСА
В настоящее время существует три основных способа получения характеристик компрессоров. Наиболее надежный и точный способ их определения - по результатам испытания на специальных стендах. Недостатками этого метода являются: необходимость специального оборудования (испытательные стенды, мультипликаторы), большие затраты энергии, необходимость иметь готовый компрессор (или его уменьшенную модель). Кроме того, на современном оборудовании невозможно получить характеристики в широком диапазоне - в области пониженных частот вращения, в зоне авторотации, неустойчивой работы не доводя компрессор до разрушения. Таким образом, на этапе разработки характеристику компрессора получить затруднительно.
Второй способ - статистический анализ характеристик компрессоров (полученных при помощи эксперимента) и по обобщенным зависимостям с некоторой долей вероятности построение характеристики конкретного (в том числе вновь создаваемого) компрессора. Недостатком этого метода является то, что характеристики принципиально новых компрессоров, не охваченные статистикой, не могут быть получены.
Третий способ - математическое описание процессов, происходящих в элементах компрессора. Достоинство метода - можно получить характеристики любого компрессора на этапе проектирования, доводки и эксплуатации. В данном направлении велись многочисленные разработки такими известными авторами [1-4], как К. В. Холщевников, Л. Н. Држинин, А. П. Тунаков, О. Н. Емин, Р. М. Федоров, Л. Е. Ольштейн и др.
2. РАСЧЕТ ХАРАКТЕРИСТИК
Существует метод расчета характеристик компрессоров, состоящий в сложении вдоль оси характеристик ступеней [4]. Этот метод не требует знания параметров решеток на отдельных радиусах и интегрирования вдоль радиуса.
Рассматривая треугольник скоростей элементарной ступени, можно написать уравнение
для коэффициента теоретического напора ступени на расчетном режиме:
нт =1 - са0 (с*§а1 + ^2) . (1)
Экспериментальные исследования показывают малое изменение углов а1 и р2. Принимая, что углы а1 (входа потока в РК абсолютном движении) и р2 (выхода из потока РК в относительном движении) остаются постоянными как для расчетного, так и для нерасчетных режимов, то для нерасчетного режима уравнение коэффициента теоретического напора примет вид:
нЛ = 1 - Са (^а1 + 2 ) . (2)
Можно связать теоретический, адиабатический напоры и КПД с помощью уравнения:
Н • @ Нпк0 (3)
нто Н0ПК Обобщая уравнения (1-3) можно получить два безразмерных комплекса, зависящих от кинематики потока:
К — _Н _ са Н0
1 * _ *
П к Са0 П к0
К 2 = Н _ Но ^.
С„п
(4)
(5)
Графики изменения К1 и К2 в зависимости
с
от и Ми приведены на рис. 1, 2 [4].
Са0
Используя К1 и К2, можно вычислить коэффициент теоретического напора и КПД для лю-
бого режима, характеризующегося и Ми.
"а0
К + ^Нм;
Са0
а0
(6)
(7)
С помощью коэффициента теоретического напора можно вычислить степень повышения давления:
1 + -
----Я • Т
к _ 1
к-1
(8)
Приведенная методика была реализована в разработанной системе имитационного моделирования С0МРКЕ880К Проведена проверка адекватности реализованной методики. Результаты приведены на рис. 3.
Рис. 3. Сравнение расчетной и экспериментальной характеристик
3. ОПРЕДЕЛЕНИЕ ГРАНИЦЫ УСТОЙЧИВОЙ РАБОТЫ КОМПРЕССОРА
Неустойчивая работа компрессора является одним из наиболее опасных дефектов газотурбинных двигателей. Причиной возникновения неустойчивой работы компрессора, в конечном счете, всегда является срыв потока с его лопаточных венцов [4]. В современных ГТД применяются в основном многоступенчатые компрессоры. Ступень же в свою очередь состоит из вращающегося рабочего колеса (РК) и неподвижного направляющего аппарата (НА) [1-4].
Схема течения потока в ступени [1-5] показана на рис. 4. Вследствие вращения РК его ло-
патки перемещаются в окружном направлении с окружной скоростью и, а воздушный поток набегает на них с относительной скоростью '^1, направление которой, характеризуемое углом Рь во избежание срыва потока должно быть близким к направлению входных кромок лопаток РК. Лопатки РК выполняются криволинейными с таким расчетом, чтобы угол выхода потока из колеса р2 был больше угла Рь Этот поворот потока, как видно из рис. 1, сопровождается увеличением поперечного сечения каждой струйки воздуха, проходящей через межлопа-точный канал (/2 > /). В результате этого относительная скорость воздуха в колесе уменьшается (^2 < w1), а давление возрастает.
Рис. 4. Схема течения в ступени компрессора
Описанная картина течения воздуха в ступени наблюдается при расчетном соотношении между окружной скоростью вращения колеса и (числом оборотов) и осевой скоростью с1а (расходом воздуха). В условиях эксплуатации по тем или иным причинам это соотношение может нарушаться. Так, например, может наблюдаться уменьшение с1а при неизменном значении п. Как видно из рис. 1, в этом случае угол р1 будет уменьшаться, а вместе с ним уменьшится и площадь поперечного сечения струйки /1. В то же время направление потока за колесом (угол Р2) мало зависит от осевой скорости - оно определяется в основном формой и направлением хвостовых частей рабочих лопаток. Следовательно, сечение струйки /2 останется прежним.
а
*
к
*
к
В результате геометрическая диффузорность межлопаточного канала, характеризуемая отношением / / /1, возрастает, и поэтому торможение потока в нем будет более сильным - ^2 упадет, а давление воздуха за колесом увеличится. Таким образом, при уменьшении с1а степень повышения давления воздуха в данной ступени (и соответственно во всем компрессоре) возрастает. Но, как видно из рис. 5 [3], при уменьшении С\а одновременно с увеличением угла Р1 возрастает угол атаки 7, с которым обтекаются лопатки РК.
р>г»1івт
Рис. 5. Изменение угла атаки при уменьшении осевой скорости
При увеличении угла атаки (рис. 6) сверх некоторого критического значения ?кр возникает срыв потока с поверхности лопаток. При срыве происходит резкое увеличение угла р2 (поток перестает следовать за направлением хвостовой части лопаток), в результате чего торможение потока в межлопаточном канале и повышение давления в нем почти прекращаются. Поэтому увеличение давления за ступенью при снижении осевой скорости воздуха происходит лишь до тех пор, пока углы атаки на рабочих лопатках (и на лопатках НА) не превышают критических.
і >0
Рис. 6. Обтекание лопатки при положительном угле атаки
На данном рис. 7 представлена типичная аэродинамическая характеристика плоской решетки профилей, которая отражает зависимости угла отклонения потока в решетке Др и коэф-
фициента гидравлического сопротивления Ср при изменении угла атаки.
На рис. 7 приведены результаты продувки плоской компрессорной решетки - зависимости угла поворота потока Др и коэффициента потерь ^ от угла атаки 7.
’ном 7КР *
Рис. 7. Результаты продувки [7]
Данная кривая не является обобщенной зависимостью угла поворота потока в решетке от угла атаки на РК и зависит в общем случае от лопаточных углов решетки р1л и р2л, густоты решетки b / t и режима течения в межлопаточном канале, определяемого числом Рейнольдса Re. В частности, данная экспериментальная кривая получена при b / t = 0,9, Р1л = 45,5°, р2л = = 90,5° и Re = 3-105.
4. РАСЧЕТ ХАРАКТЕРИТИКИ КОМПРЕССОРА СОВМЕСТНО
С ГРАНИЦЕЙ УСТОЙЧИВОЙ РАБОТЫ
Методика расчета характеристик осевых многоступенчатых компрессоров совместно с границей устойчивой работы реализована авторами в системе COMPRESSOR. На рис. 8 приведен внешний облик системы с библиотекой элементов.
В библиотеке элементов имеется элемент с названием «Характеристика», с помощью которого можно выполнять расчет характеристики ступени и всего компрессора. На рис. 9 приведена схема модели для расчета характеристики трехступенчатого компрессора в системе COMPRESSOR.
Каждый элемент имеет входные параметры. Элемент «Вход» имеет следующие входные параметры:
• расход воздуха;
полная температура; полное давление; число ступеней;
частота вращения на расчетном режиме.
Рис. 8. Внешний вид системы COMPRESSOR [6]
Иерархия Текущий уровень
1 одель jj Начальные условия 3 Вход (характеристика) & Характеристика 9- Характеристика 2 Б в—— 1>
JJ Общие результаты < >
\Иерархия АДпя заметок/
Рис. 9. Схема трехступенчатого компрессора [6]
Элемент «Характеристика» имеет следующие входные параметры:
• наружные диаметры (вход и выход из РК, выход из НА);
• относительные диаметры втулки (вход и выход из РК, выход из НА);
• осевая скорость на входе в РК (расчетная);
• осевая скорость на выходе из РК (расчетная);
• коэффициент неравномерности;
• входной лопаточный угол;
• выходной лопаточный угол;
• угол потока в абсолютном направлении на входе в РК;
• угол потока в относительном направлении на выходе из РК;
• критический угол атаки.
После задания всех необходимых входных параметров можно производить расчеты.
В качестве примера приведены результаты расчета трехступенчатого осевого компрессора в системе COMPRESSOR.
Система позволяет получать как характеристики отдельных ступеней, так и характеристи-
ку компрессора в целом. На рис. 10 приведена расчетная характеристика компрессора.
В первом приближении критический угол атаки принимается равным 80 для всех ступеней. На рис. 11-13 приведены зависимости углов атаки в ступенях от расхода.
На рис. 14-15 показаны сводные зависимости углов атаки всех ступеней от расхода отдельно для низконапорных и высоконапорных веток характеристики.
Проанализировав зависимости (рис. 11-15), можно сделать вывод, что в области малых частот вращения помпаж вызывается срывом потока с лопаточных венцов первой ступеней, а в области высоких - срывом потока в последних ступенях. Данный вывод соответствует реальной физической картине.
Имея зависимости угла атаки от расхода и зная критический угол атаки, можно обнаружить границу устойчивой работы компрессора при фиксированной частоте вращения, то есть выявить соответствующий границе устойчивости расход. В табл. 1 приведены расходы, соответствующие границе устойчивой работы компрессора.
Таблица 1
Расходы воздуха, соответствующие границе устойчивой работы компрессора
Ступень 1 2 3
^Впр
«ПО = 30% 0,27 0,19 0,15
«по = 50% 0,45 0,33 0,28
«пр = 70% 0,61 0,49 0,47
«пр = 80% 0,69 0,59 0,58
ппр = 90% 0,76 - 0,71
ппр = 95% 0,8 - 0,79
«пр = 100% - - 0,88
«пр = 105% - - 0,98
«пр = 108% - - 1,03
Из табл. 1 необходимо взять максимальные расходы на каждой частоте вращения, а затем с помощью них построить границу устойчивости. На рис. 16 показан алгоритм построения границы устойчивой работы.
На рис. 17 приведено сравнение расчетной и экспериментальной границ устойчивой работы исследуемого компрессора.
Как видно из рис. 17, расчетная и экспериментальные границы устойчивой работы значительно отличаются, следовательно, необходимо уточнение методики определения границы устойчивой работы компрессора. Расхождение может быть связано с тем, что результаты продувки соответствовали вполне конкретной решетке, которая отличается от используемых в модели решеток.
Рис. 10. Расчетная характеристика компрессора
Рис. 12. Зависимости углов атаки от расхода на разных частотах для второй ступени
Рис. 13. Зависимости углов атаки от расхода на разных частотах для третьей ступени
Рис. 16. Построение границы устойчивой работы
—»-экспер и мент расчет (¡=8)
Рис. 17. Сравнение расчетной и экспериментальной границ устойчивой работы
5. УТОЧНЕННЫИ РАСЧЕТ ГРАНИЦЫ УСТОЙЧИВОЙ РАБОТЫ КОМПРЕССОРА
В исследуемой модели проводились расчеты для нескольких критических углов атаки. На рис. 18 приведены границы устойчивости для различных углов атаки.
Как видно из рис. 18, на пониженных приведенных частотах вращения в первых ступенях достигаются значительные величины углов атаки, а на последних ступенях - отрицательные (табл. 1, рис. 11-13). Данный факт подтверждается приближением расчетной границы устойчивости к экспериментальной кривой на пониженных приведенных частотах вращения при повышении величины критического угла атаки, а также на повышенных приведенных частотах вращения при уменьшении величины критического угла атаки. Соответственно для удовлетворительной сходимости расчетной и экспериментальной границ устойчивой работы компрессора требуется дискретное задание ограничений по критическим углам атаки для различных ступеней - более высокие значения требуются для группы первых ступеней, более низкие - для группы последних ступеней.
В качестве методики построения границы устойчивой работы на характеристике используется ограничение по критическому углу атаки на рабочем колесе. Для упрощения расчетной модели возможно задание постоянного крити-
ческого угла атаки по высоте пера лопатки и по различным ступеням каскада компрессора. Однако такое допущение вносит значительную методическую погрешность в расчет и в условиях тре бований высокой точности недопустимо. Конкретные значения критических углов атаки можно вычислить с помощью кривой, разработанной Хауэллом, приведенной на рис. 19.
Для проведения анализа работы компрессора на нерасчетных режимах работы задаются номинальным углом поворота потока в решетке профилей на данном радиусе [7] (в работе проводится анализ только на среднем диаметре).
Номинальный угол поворота потока связан с номинальным следующим соотношением:
Лрном = 0,8 ■ Лркр. Данному номинальному углу поворота потока в решетке соответствует определенное значение угла атаки /ном.
На рисунке представлены зависимости:
г - г
лр г лр
^ном ^ном
ном , Ь );^_ _ у ( ном Ь)
лр і
ном
Анализ данных кривых указывает на то, что нерасчетный угол отклонения потока в решетке профилей не зависит от густоты решетки, а коэффициент сопротивления решетки имеет расслоение характеристик по густоте решетки профилей. Данная обобщенная характеристика имеет следующие ограничения по области применения: Яе < 105 и Яе > 3 • 105.
■ эксперимент * расчет(¡=6) » расчет (¡=8)
■ расчет (¡=10) ---расчет0=11) ---расчет(!=13)
Рис. 18. Сравнение границ устойчивости с различными критическими углами атаки
АР _ г(? ?ном Ь)
Ар (Ар , г
ном ном
Этот факт связан с образованием отрывной зоны при превышении критического угла атаки, что приводит к резкому снижению угла поворота потока в решетке профилей.
Таким образом, наиболее оптимальным и экспериментально обоснованным критерием для выявления границы устойчивой работы компрессора является следующее равенство:
г -1
ном »04
АР
ном
На этапе проектирования компрессора затруднительно произвести определение максимального, а соответственно, и номинального угла поворота потока в решетке профилей. Поэтому в данной работе за номинальные значения угла отклонения потока в лопаточной решетке и угла атаки принимаются соответствующие значения в проектной точке. Применение метода Хауэлла позволяет в значительной степени повысить адекватность математической модели и повысить точность определения границы устойчивой работы на характеристике компрессора за счет дифференцированного определения критического угла атаки для каждой ступени каскада. После проведения расчета с дискретным распределением критических углов атаки была получена граница устойчивости, изобра-
Afi
ЬР..
НО. W
1.2
£ в уЛ— • ч>
ом 1 У
/ [
/ 1
У 1
7
L D V /
ъ “М- t * У /
^77 МЫ -0.2
+0.2 * 0.4 *0.6
Д р
ном
Рис. 19. Обобщенная характеристика различных решеток на нерасчетных режимах [7]
Коэффициент гидравлического сопротивления решетки связан с коэффициентом сопротивления решетки следующим соотношением:
*ё(а т) = 2 • (Ш( а1) + *§(а 2 ))>'
Cr
ч3 ■
(b/t) • (cos(a т))
Граница устойчивой работы компрессора определяется точкой максимума на зависимости:
женная на рис. 20 и характеристика, приведенная на рис. 21. В табл. 2 приведены расходы, соответствующие границе устойчивой работы компрессора.
Т аблица 2 Расходы воздуха, соответствующие границе устойчивой работы компрессора
Ступень 1 2 3
^Впр
«пр = 30% 0,2 0,19 0,16
«пр = 50% 0,33 0,33 0,3
«пр = 70% 0,45 0,49 0,5
«пр = 80% - 0,59 0,62
«пр = 90% - - 0,76
«пр = 95% - - 0,84
«пр = 100% - - 0,94
«пр = 105% - - 1,02
«пр = 108% - - 1,05
Как показывает анализ кривых приведенных на рис. 20-21, дискретное расчетное задание ограничения по критическому углу атаки позволяет получить удовлетворительное совпадение расчетной и экспериментальной границ устойчивости компрессора.
ВЫВОДЫ
В системе имитационного моделирования COMPRESSOR реализована методика построения характеристик осевых многоступенчатых компрессоров путем «сложения характеристик» отдельных ступеней. Данная методика позволяет получать достаточно точные результаты в широкой области приведенных частот вращения при ограниченном количестве исходных
данных на начальных этапах проектирования за счет базирования на одномерном расчете по среднему диаметру и обобщенных экспериментальных зависимостях.
Для определения границы устойчивой работы многоступенчатого компрессора используется ограничение по критическому углу атаки и обобщенная кривая Хауэлла для анализа работы каскада компрессора на нерасчетных режимах работы. Этот метод позволяет с достаточно высокой достоверностью получать границу устойчивой работы каскада компрессора. Таким образом, разработана комбинированная методика построения характеристик осевых многоступенчатых компрессоров в широкой области приведенных частот вращения совместно с границей устойчивой работы компрессора.
В настоящее время ведутся работы по дальнейшему уточнению методики построения характеристик компрессоров совместно с границей устойчивости, а именно: работы по определению срыва потока на различных радиусах по высоте пера лопатки (это позволит повысить точность определения границы устойчивой работы компрессора, так как срыв потока, как правило, происходит во втулочных сечениях лопатки); уточнение методики определения потерь в лопаточной решетке, что позволит повысить точность построения напорных веток на характеристике каскада компрессора. Также ведутся работы по моделированию рабочих процессов многоступенчатых осевых компрессоров на неустойчивых режимах работы.
Рис. 20. Сравнение уточненной расчетной и экспериментальной границ устойчивости
—♦—расчет эксперимент
Рис. 21. Сравнение расчетной и экспериментальной характеристик
СПИСОК ЛИТЕРАТУРЫ
1. Ржавин Ю. А. Осевые и центробежные компрессоры двигателей летательных аппаратов. Теория, конструкция и расчет: Учеб. М.: Изд-во МАИ, 1995. 344с.
2. Емин О. Н., Карасев В. Н., Ржавин Ю. А.
Выбор параметров и газодинамический расчет осевых компрессоров и турбин авиационных ГТД: Учеб. пособие. М.: Дипак, 2003. 156 с.
3. Проектный термогазодинамический расчет основных параметров авиационных лопаточных машин / А. Н. Белоусов [и др.]. Самар. гос. аэрокосм. ун-т, 2006. 316 с.
4. Холщевников К. В. Теория и расчет авиационных лопаточных машин: Учеб. для авиац. вузов и факультетов. М.: Машиностроение, 1970. 610 с.
5. Галимзянов Ф. Г. Термодинамические и газодинамические расчеты авиационных ТРД: Лопаточные машины. Уфа: УАИ, 1978. 100 с.
6. Ахмедзянов Д. А., Козловская А. Б., Кри-вошеев И. А. Система моделирования компрессоров авиационных ГТД (COMPRESSOR): Свид-во об офиц. рег. № 2009612688, Роспатент, М., 2009.
7. S. L. Dixon, B. Eng., P. H. D. Fluid Mechanics, Thermodynamics of Turbomachinery, University of Liverpool. 1998, 320 p.
ОБ АВТОРАХ
Козловская Александра Борисовна, асп. каф. авиац. двиг. Дипл. инж. по авиац. двиг. и энерг. установкам (УГ АТУ,
2008). Иссл. в обл. процессов в компрессорах авиац. ГТД с использованием имитац. и 3Б-САБ/САЕ-моделирования.
Ахмедзянов Дмитрий Альбертович, проф. той же каф., зам. декана ФАД. Дипл. инж. по авиац. двигателям и энерг. уст. (УГАТУ, 1997). Д-р техн. наук по тепл., электроракетн. двигателям и энергоустановкам ЛА (УГАТУ, 2007). Иссл. в обл. рабочих процессов в авиац. ГТД, разработки матем. моделей
сложн. техн. объектов, САПР авиац. ГТД.
Ахметов Юрий Мавлютович,
доц. каф. прикл. гидромех., зам. ген. дир. НИИТ. Дипл. инж.-мех. по авиац. двиг. (УАИ, 1959). Канд. техн. наук по тепловым двигателям (МАИ, 1978). Иссл. в обл. газогидр. течений и систем управления энергоустановок.
Михайлов Алексей Евгеньевич, асп. каф. авиац. двиг. Дипл. инж. по авиац. двиг. и энерг. установкам (УГАТУ, 2010).
Иссл. в обл. рабочих процессов в авиац. ГТД на установившихся и неустановившихся режимах.