Уфа : УГАТУ, 2011
Ъыямшь QjrAQnQj
Т. 15, № 4 (44). С. 18-24
ТЕПЛОВЫЕ, ЭЛЕКТРОРАКЕТНЫЕ ДВИГАТЕЛИ И ЭНЕРГОУСТАНОВКИ ЛА
УДК 621.452
Г. Г. ЕНИКЕЕВ, А. Т. БИКМЕЕВ, А. А. СОЛОВЬЕВ, А. А. КАСАТКИН, К. Р. ЮЛМУХАМЕТОВ
ОСОБЕННОСТИ МОДЕЛИРОВАНИЯ ДВИЖЕНИЯ МНОГОФАЗНОГО ПОТОКА В ВОЗДУХОЗАБОРНИКЕ ГАЗОТУРБИННОГО ДВИГАТЕЛЯ
Процесс движения многофазного потока в проточной части газотурбинного двигателя чрезвычайно сложен и недостаточно изучен. В воздухозаборном устройстве, компрессоре, камере сгорания, турбине, сопле имеются свои особенности движения многофазного потока. Выполнен анализ наличия математических моделей в программных комплексах (ПК) STAR-CD, STAR-CCM , ANSYS CFX, FLUENT для решения задач движения многофазных дисперсных потоков в газово-воздушном тракте (ГВТ) газотурбинного двигателя (ГТД). Показано, что при небольшой концентрации дисперсных частиц в воздухе наиболее рациональным является определение поля течения газа в постановке Эйлера и решение задачи о движении дисперсных частиц в постановке Лагранжа. В качестве примера решена пространственная задача обтекания лобового воздухозаборника сжимаемым, вязким турбулентным потоком воздуха в постановке Эйлера. Решение выполнено для воздухозаборника NACA 1-70-50. Показано хорошее совпадение результатов численного моделирования с экспериментом. Газотурбинный двигатель; газово-воздушный тракт; многофазный дисперсный поток; сжимаемый, вязкий турбулентный поток; математическая модель; численное моделирование; программные комплексы; тестовые задачи
Эксплуатация газотурбинных двигателей практически всегда осуществляется при наличии в воздухе дисперсных частиц. В большинстве случаев их концентрация настолько мала, что их присутствием можно пренебречь, так как их влияние на характеристики двигателя несущественно. Но при эксплуатации ГТД в запыленной атмосфере или в морской среде происходит значительное изменение их характеристик из-за эрозионного износа лопаток компрессора или отложения соли в проточной части и коррозии, что сопряжено со снижением надежности и безопасности эксплуатации транспортных устройств.
Механика движения многофазного полидисперсного потока в проточной части ГТД недостаточно изучена. Исследования рабочего процесса ГТД, эксплуатирующихся в запыленной атмосфере и морской среде, дают представление о физических процессах, происходящих при движении многофазного потока в проточной части газовоздушного тракта. Каждый ГТД имеет свои индивидуальные особенности, и для прогнозирования изменения характеристик, а также проектирования воздухозаборников с воздухоочистителем, в условиях эксплуатации в запыленной атмосфере и морской среде требуется проведение дополнительных теоретических и экспериментальных исследований. Поэтому исследование движения многофазных потоков в газовоздушном тракте ГТД является
актуальными.
Практика показывает, что современные программные комплексы и их реализация на компьютерах достигли совершенства и их широко используют для решения сложных задач движения многофазных сред.
1. МАТЕМАТИЧЕСКИЕ МОДЕЛИ
ПК В РЕШЕНИИ ЗАДАЧИ ОБТЕКАНИЯ МНОГОФАЗНЫМ ПОТОКОМ ЛОБОВОГО ВОЗДУХОЗАБОРНИКА С РОТОРНЫМ ВОЗДУХООЧИСТИТЕЛЕМ
Эффективным способом защиты ГТД от попадания в проточную часть дисперсных частиц является роторный воздухоочиститель, который расположен в модульной конструкции воздухозаборника [1-6]. Воздухозаборник с роторным воздухоочистителем имеет проточную часть с неподвижными поверхностями и лопаточные аппараты вращающегося рабочего колеса и неподвижного направляющего аппарата, то есть для исследования движения многофазного потока в воздухозаборном устройстве требуется весь набор математических моделей, который требуется для анализа процессов в проточной части газовоздушного тракта двигателя.
При движении дисперсного потока в воз-духозаборном устройстве с роторным воздухоочистителем твердые частицы при встрече с неподвижной или вращающейся поверхностью отскакивают и двигаются по новым траекториям.
Контактная информация: enikeevgg@gmail.com
Математические модели программных комплексов для решения задачи движения дисперсного потока в воздухозаборном устройстве с роторным воздухоочистителем
STAR-CD STAR-CCM+ ANSYS CFX ANSYS FLUENT
Расчет вращающихся лопаток -сплошная среда Есть Есть Есть Есть
1. Вращающаяся система координат (MRF) 1. MRF 1. MRF 1. MRF
2. Программируемые подвижные сетки 2. Вращение сетки 2. Вращение сетки 2. Вращение сетки
3. Immersed Solid -движение тела сквозь ячейки среды 3. Программируемые подвижные сетки
Построение сетки для области вращения Модуль ES-Turbo Встроенный генератор сеток + + Turbo Wizard, нет редактора профилей Ansys TurboGrid - наиболее мощный инструмент построения сеток для турбомашин
Теплообмен со средой Есть + программирование Есть Есть Есть
Массообмен со средой Есть + программирование Есть Есть Есть
Распад частиц в потоке Reitz-Diwakar Pilch-Erdman Hsiang-Faeth Reitz-Diwakar TAB, KHRT Reitz-Diwakar TAB, ETAB, CTAB TAB Волновая
Отскок от стенки, дробление и налипание Есть Модели MPI, Bai + + пользоват. подпрограмма Есть Модель Bai-Gosman Есть Модель Elsaesser Есть Модель стенка-струя
Жидкая пленка на стенках Есть, подвижная Только спец. случаи -запотевание и обледенение Есть, неподвижная Есть, подвижная
Эрозия Только пользоват. программирование Есть Есть Есть
Капли при определенных условиях испаряются, дробятся, образуя ансамбль частиц с меньшими размерами. При испарении соленых капель может образоваться пересыщенный раствор.
При полном испарении соленые капли превращаются в твердые частицы. При встрече с неподвижными или вращающимися поверхностями капли прилипают и образуют пленку, которая движется под действием сил тяжести, центробежной, Кориолиса, трения о поверхность лопатки, трения газа о поверхность пленки. Могут образоваться вторичные капли за счет разбрызгивания и срыва пленки. Пленка может испаряться, образуя слой соли на поверхности лопаток.
В ПК Star-CD, Star-CCM+, Ansys CFX, Ansys Fluent заложены математические модели,
описывающие перечисленные физические процессы (таблица).
Анализ уравнений показал, что не все математические модели в отдельно взятом ПК в полной мере отражают физические процессы при движении многофазного потока в воздухозаборнике. Поэтому при решении конкретной задачи выбирается математическая модель и ПК, в котором адекватно реализуется тот или иной физический процесс. Полученные результаты импортируются в то ПК, в котором продолжается решение задачи.
Концентрация дисперсных частиц в воздухе мала и составляет не более 3...5% при эксплуатации вертолетов в запыленной атмосфере, эк-ранопланов, кораблей на воздушной подушке, быстроходных судов в морской среде. Для малой концентрации дисперсных частиц в воздухе
является обоснованным численное моделирование движения непрерывного, несущего потока в постановке Эйлера. Моделируется обтекание воздухозаборника сжимаемым, вязким турбулентным потоком воздуха и полученные результаты используются для решения задачи движения дисперсных частиц в постановке Лагранжа. Такая методика позволяет получить информацию о параметрах каждой движущейся частицы, которая является важной при анализе рабочего процесса движения многофазного потока в воз-духозаборике с роторным воздухоочистителем и в ГВТ ГТД.
2. ОБТЕКАНИЕ ВОЗДУХОЗАБОРНИКА НЕПРЕРЫВНЫМ ПОТОКОМ
2.1. Постановка задачи
Решена задача движения непрерывного осесимметричного потока при обтекании лобового воздухозаборника NACA 1-70-50. Результаты были представлены в сравнении с экспериментом и решением, полученным панельным методом (В. Краус). Задача решалась для следующих условий: газ сжимаемый, вязкий, турбулентный, рабочее вещество - воздух. Для решения задачи использовались полные уравнения
Навье - Стокса, дополненные двухслойной k-s-моделью турбулентности. Программная реализация выполнена в ПК ANSYS CFX и Star-CCM+ в постановке Эйлера. В невозмущенной области скорость воздуха - V0, давление - P0, плотность р0. В узком сечении внутреннего канала воздухозаборника скорость - V, давление - P, плотность p¿. Численные расчеты выполнялись для двух отношений скоростей воздуха в узком сечении воздухозаборника и в невозмущенной области: VjV0 = 0,216 и VjV0 = 0,428.
Рассчитывалось распределение давления по внешней поверхности обечайки
воздухозаборника, которое определялось коэффициентом давления СР. Коэффициент давления вычислялся по формуле:
Р - P СР = Р P0
2 р^
где Р г - давление на поверхности; Р0 - давление в невозмущенном потоке; У0 - скорость в невозмущенной области. При подготовке к решению задачи выполнялись стандартные процедуры импорта геометрии, создания расчетной области, построения сеток и задания граничных условий.
2.2. Геометрические размеры воздухозаборника
Схема воздухозаборника и основные геометрические размеры приведены на рис. 1. Исследуемый воздухозаборник относится к серии NACA 1 с отношением площадей Amin/ Am = 0,7 и отношением длины к
максимальному диаметру LjDm = 0,5 .
Основные соотношения для расчета геометрических параметров воздухозаборника серии NACA 1-70-50 с параметрами kb k2 (в
нашем случае kx = 0,7, k2 = 0,5 ) следующие:
R0 = k1Rm ,
Rmin = R0 - 0,025(Rm - R),
L = 2k2 Rm,
y'
Rm - R0
-' (L
где Д0 - радиус носка воздухозаборника. Форма внешней поверхности задаётся таблицей значений функции / и является поверхностью постоянного радиуса со скруглением на конце (дуга окружности). Внешний радиус Ят выбран равным 0,6 м.
У
задано таблицей
Я
дуга
окружности.
Rmii
x
Rm
Рис.1. Геометрия воздухозаборника NACA 1-70-50
Для построения геометрии воздухозаборника в пакете Star-CCM+ написана программа-макрос на Java. Программа считывает параметры таблицы из внешних файлов, строит кривую из отрезков сплайнов и дуг и путём её вращения вокруг оси создает тело - область течения. При этом внешние границы области течения отодвигаются от воздухозаборника на 4.. .5 калибров.
2.3. Расчетная область, граничные условия и построение сетки
Расчетная область выполнена для четвертой части воздухозаборника, в целях экономии машинных ресурсов (рис. 2). Границы расчетной области отодвинуты от входа в воздухозаборник на пять калибров, для исключения влияния возмущения поля скоростей и давлений на результаты решения задачи.
L
На выходе из канала воздухозаборника задана скорость V.] = 30 м/с , а на входе в расчетную область - V0, соответствующая выбранным значениям Vi|V0 . Заданы две плоскости симметрии. Свободный поток определяется атмосферным давлением Р0. В расчетной области построена сетка из 593 тыс. многогранных ячеек с призматическими слоями вдоль стенок воздухозаборника (рис. 3). Вблизи особой точки и поверхности обечайки выполнено сгущение ячеек для увеличения точности расчетов.
Рис.2. Геометрия расчетной области и граничные условия
Рис. 3. Сетка в расчетной области
2.4. Результаты вычислений
Результаты вычислений обтекания воздухозаборника представлены в виде визуализации течения и в виде количественных зависимостей распределения давления по внешней поверхности.
Линии тока в расчетой области показаны на рис. 4. Расчеты в ПК 8гаг-ССМ+ и в ПК ЛпБуБ СБХ дали схожие результаты. Вблизи носка обечайки воздухозаборника происходит разделение потока и увеличение скорости и при переходе к цилиндрической части его торможение.
б
Рис.4. Линии тока в расчетной области для режима течения V¡|V0 = 0,428 : а - расчет в ПК 81аг-ССМ+; б - расчет в ПК АшуБ СБХ
ргтиге Рйг* 1
_ 1.29№*ИИ
6
[Ра|
Рис.5. Распределение давления в расчетной области для режима течения Vi|V0 = 0,428
а
-0.5
-0.45
а -0.4
и
и -0.35
X
01
h -0.3
и
£ -0.25
01
S 3 -0.2
ф -0.15
н
ЗС -0.1
-0.05
0
-1.2
-0.8
-0.2
0 0 10 20 30 4 0 50 6 0 7 0 8 0
Длина по оси воздухозаборника, м
а
эксперимент
панельный метод
расчет в ПК Ansys CFX
расчет в ПК Star CCM+
i
т»и "ч
00 10 20 3 0 40 5 0 60 7 0 .8 0.
-эксперимент
-панельный
метод - расчет в ПК Ansys CFX Star CCM +
Длина по оси воздухозаборника,м
б
Рис. 6. Распределение давления по внешней поверхности воздухозаборника NACA 1-70-50 а - режим течения V Дд = 0,428 ; б - режим течения Vi jVg = 0,216
3 Ф qj го
ГО П5
I Й ф Д
S 2
ш У
i* 2 о
(—
О
Длина по оси воздухозаборника, м
-расчеты панельным методом
-расчеты вПК
Ansys CFX
-расчеты вПК
Star-CCM+
Длина по оси воздухозаборника, м
б
- расчеты панельным методом
- расчеты вПК
Ansys CFX
-расчеты вПК
Star-CCM+
Рис.7. Точность вычислений распределения давления по внешней поверхности воздухозаборника
NACA 1-70-50 в сравнении с экспериментом:
а - режим течения V'/V) = 0,428 ; б - режим течения V'/V) = 0,216
-1.6
= 0.216
-1.4
-1
-0.6
-0.4
9
0
а
Распределение давления в расчетной области показано на рис. 5. В области криволиненейной внешней поверхности обечайки происходит уменьшение давления, вызванное ускорением потока. Во внутреннем канале воздухозаборника происходит торможение потока.
Количественные результаты расчета представлены в виде зависимостей
коэффициента давления Ср от расстояния по оси воздухозаборника (рис. 6). На графиках (рис. 6.) видно хорошее совпадение результатов расчета в ПК Star-CCM+ и Ansys CFX с экспериментом для режимов обтекания V/V0 = 0,428 и V/V0 = 0,216. На этом же рисунке приведены результаты, полученные В. Краусом панельным методом. В диапазоне X = 0,2.. .0,4 (режим V / V = 0,428) панельный
метод дает большее расхождение с экспериментом. В диапазоне X = 0,1.0,15 (режим V/V = 0,428) все три результата вычислений показывают некоторое
расхождение с экспериментом. Объясняется это близким расположением к носку воздухозаборника.
Для режима течения V /V = 0,216
расчетные методы дают хорошее совпадение с результатами эксперимента по всей поверхности за исключением области вблизи носка воздухозаборника.
Более точно о совпадении результатов численных расчетов с экспериментом можно судить по графикам отклонений (рис.7). В диапазоне X = 0,03.0,2 отклонения коэффициента давления составляют ±5%. В диапазоне X = 0,18.0,6 отклонения коэффициента давления составляют 0.3% для расчетов для режима течения V Ао = 0,428 .
Для режима течения VjV0 = 0,216 решения в ПК Star-CCM+ и в ПК Ansys CFX дают более точный результатат, чем панельный метод.
ЗАКЛЮЧЕНИЕ
Выполнен анализ математических моделей, приведенных в программных комплексах (ПК) Star-CD, Star-CCM+, Ansys CFX, Ansys Fluent, для решения задачи движения многофазного потока при обтекании лобового воздухозаборника с роторным воздухоочистителем. В перечисленных ПК математические модели практически полностью охватывают физические процессы, происходящие при движении многофазного потока в воздухозаборном
устройстве с роторным воздухоочистителем. Но если взять каждый ПК в отдельности, то не все физические процессы описаны в полной мере. Поэтому можно использовать преимущества каждого пакета для решения задачи в полном объеме.
Небольшая концентрация дисперсных частиц в потоке позволяет выполнять расчеты непрерывного потока при обтекании воздухозаборника в постанове Эйлера и затем использовать полученные результаты для расчета движения дисперсных частиц в постановке Лагранжа.
Решена задача моделирования обтекания воздухозаборника NACA 1-70-50 непрерывным потоком сжимаемого, вязкого турбулентного газа.
Результаты решения показаны в сравнении с экспериментом и вычислениями панельным методом.
Выявлено, что решение задачи распределения давления по внешней поверхности обечайки воздухозаборника в ПК Star-CCM+ и Ansys CFX совпадает с экспериментом с отклонением до 5% для режимов течения V/V0 = 0,428 и V/V0 = 0,216 .
Решение задачи позволяет выявить визуальную картину течения, линии тока и распределение давления в расчетной области.
СПИСОК ЛИТЕРАТУРЫ
1. Еникеев Г. Г. Математические модели эффективности роторных воздухоочистителей ВРД // Авиационно-космическая техника и технология. Харьков: ХАИ, 2006. С. 33-37.
2. Защита ГТД от морских аэрозолей с помощью роторных воздухоочистителей / Г. Г. Еникеев, Р. М. Гузаиров, В. Л. Химич, Г. З. Шарипов // Сборник тр. Уфимск. ордена Ленина авиац. ин-та имени Серго Орджоникидзе. Уфа: УАИ, 1992. С. 132-151.
3. Гузаиров А. Р., Еникеев Г. Г. Динамика твердых частиц пыли в ступени роторного воздухоочистителя // Проблемы и перспективы развития двигателестроения в поволжском регионе: тез. докл. междунар. научно-техн. конф., посвященной 55-летию Самарск. гос. аэрокосмич. ун-та. Самара: СГАУ, 1997.
4. Бичер Б., Табаков У., Хамед А. Усовершенствованный метод расчета траекторий частиц угольной золы в турбомашине // Тр. амер. об-ва инженеров-механиков. Сер. А, Энергетические машины и установки. 1982. №1. С.43-47.
5. Защита ГТД от морских аэрозолей с помощью роторных воздухоочистителей / Г. Г. Еникеев, Р. М. Гузаиров, В. Л. Химич, Г. З. Шарипов // Сборник трудов УАИ. Ч. 1. Уфа: УАИ, 1992. С. 132-151.
6. Численные методы в динамике жидкостей / Э. Джеймсон и др., под ред. Вирц Г., Смолдерен Ж. М.: Мир, 1981. 407 с.
ОБ АВТОРАХ
Еникеев Галей Гумерович, доц. каф. прикл. гидромех. (ПГМ) УГАТУ. Дипл. инж-мех по авиац. двиг. (УАИ, 1967). Канд. техн. наук по тепл. двиг. (УАИ, 1979), доцент. Иссл. в обл. механики многофазных сред.
Бикмеев Александр Тимерзянович, доц. каф. высокопроизводительных вычислительных технологий и систем (ВВТиС) УГАТУ. Дипл. физика (БГУ, 1995). Канд. физ.-мат. наук по физике конденсированного состояния (УГАТУ, 2002) Иссл. в обл. физики трения.
Соловьев Алексей Александрович, ст. преп. каф. ПГМ. Дипл. магистр техники и технологий по энергомашиностр. (УГАТУ, 2005). Канд. техн. наук по гидравл. машинам и гидро-пневмоагр. (УГАТУ, 2008). Иссл. в обл. газодинамики и термодинамики закрученных и многофазных потоков.
Касаткин Алексей Александрович, ст. преп. каф. ВВТиС. Дипл. инж.-матем. (УГАТУ, 2007). Иссл. в обл. диф. уравнений.
Юлмухаметов Константин Раисович, мл. науч. сотр. ИКИ НИЧ УГАТУ, дипл. физика (БГУ, 2003). Иссл. в обл. моделирования физических процессов.