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

УДК 621.18; 621.16; 621.165

ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ

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

© 2011 М. В. Малыгина Рыбинская государственная авиационная технологическая академия имени П. А. Соловьёва

Представлены результаты исследований рабочего процесса многотопливного горелочного модуля, формирующего первичную зону горения камеры сгорания ГТУ. Найдены соотношения расходов воздуха через контуры закручивающих устройств, обеспечивающие требуемое распределение коэффициента избытка воздуха по объему зон предварительного смешения. Определен угол закрутки потока воздуха, обеспечивающий компромисс противоречивых факторов, таких как величина гидравлических потерь и объем приосевой зоны обратных токов, определяющий условия стабилизации пламени при горении бедных смесей (а>1,5).

Камера сгорания, горелочный модуль, исследование, расчет, моделирование.

Введение

Развитие авиационных и наземных двигателей направлено на повышение их термического КПД за счет форсирования параметров газа по давлению и температуре на выходе из камеры сгорания. В настоящее время для двигателей наземного применения освоена температура на входе в турбину 1600-1700К. Развиваются проекты высокотемпературных камер сгорания для гражданских авиационных двигателей и стационарных энергокомплексов с температурой на входе в турбину до 2150К. Тенденции развития энергетики приводят к разработке малотоксичных камер сгорания, которые наряду с высокими экономическими и эксплуатационными показателями должны соответствовать жестким требованиям по выбросам загрязняющих атмосферу веществ (оксидов азота, продуктов неполного сгорания, монооксида углерода, дыма) в широком диапазоне изменения режимных параметров. Ведущие компании, такие как GE, Siemens, Alstom, предлагают газотурбинные установки, обеспечивающие уровень выбросов NOx и СО меньше 25 ppm [4].

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

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

В работе предложена противоточная схема течения рабочего тела в объеме горе-лочного модуля. Преимущества рабочего процесса таких устройств, а также характерные особенности конструкции исследованы и приведены в работах [1-3].

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

Особенности противоточных закрученных течений вызывают интерес к исследованию процессов горения в них. Это является

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

1. Численное моделирование термогазодинамики и горения в противоточном горелочном модуле

1.1. Постановка задач и метод численного исследования

Конструктивная схема многотопливного горелочного модуля представлена на рис.

1. Она включает противоточную смесительную камеру, двухконтурное закручивающее устройство, системы подачи жидкого (керосин, дизельное топливо) и газообразного (природный газ, синтез-газ) топлив. Горе-лочный модуль рассчитан на сжигание основного топлива природного газа и резервного дизельного и синтез-газа.

Второй контур закручи віющего устройства

Первый контур закручивающего устройства

Форсунка подачи керосина и ДТ

Воздух от компрессора

Форсунка подачи природного и синтез газов

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

2D геометрическая модель и расчетная сетка показаны на рис.2. Полученная сеточная модель имеет блочную структуру со сгущением в областях впрыска топливных струй.

Рис. 1. Конструктивная схема многотопливного про-тивоточного горелочного модуля

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

Рис. 2. Расчетная сетка в проточной части противоточного горелочного модуля

Численное моделирование термогазодинамики вязкого газа выполнено с использованием системы уравнений Новье-Стокса, осредненных по Рейнольдсу в стационарной постановке. Моделирование турбулентного переноса осуществлялось с помощью стандартной ^8 модели турбулентности.

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

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

ЕІ !

еЕ

и

от относительного перепада давления на горелочном модуле и, определен уро-

Рис. 3. Схема расчетной области противоточ-ного горелочного модуля

вень гидравлических потерь, для оценки которой в литературе принято выражение [5]

* *

а - Гаи ■ 100%.

V ЕХ

Критерием оценки качества смесеобразования в противоточном горелочном модуле выбрано изменение среднеквадратичного отклонения £ его массовой концентрации от среднерасходного значения в каждом сечении:

где - массовый расход топлива через

элементарную площадку - среднерас-

ходное значение массовой концентрации керосина в данном сечении.

Постановка задачи для изучения смесеобразования учитывает направления скоростей потоков воздуха при их прохождении через тангенциальное и радиальное закручивающие устройства. Рассмотрены случаи с сонаправленной и противоположно направленной закруткой потоков. Угол закрутки в тангенциальном закручивающем устройстве составляет а =15°. Для аксиального закручивающего устройства было выбрано 8 значений угла входа потока: в = ±15°, ±30°, ±45° и ±60° (где «+» - обозначены углы, формирующие вращение по часовой стрелке, а «-» - против часовой).

Изобразим направления векторов скорости потока 1#! и мл при прохождении через тангенциальное и аксиальное закручивающие устройства соответственно (рис. 4).

Результаты исследований. Расчет термогазодинамики течения и горения в ка-

Рис. 4. Ориентация векторов скорости потока и на входе в тангенциальное и аксиальное закручивающие устройства: у = 15":«*! = бЕм/о 1^, = 0; 1Г1;_ = -17,6иУ^чг1г = 65,7 зч/о Р = 15и,-, = 9 0 п/а = 36,Эм/а = 2 3,3 тч/а = 0;

/3 = 9 (У3; и--: = 90 м/с; и--,;, = 77,9 м/с; иг,,. = 45 и/ц; =0;

Р = 45*11*-*= 90 п/а = 63,7м/с; = £3,7и/с; =0;

Р =НР:», = 90 п/а и--,;, = 45 к/а Щц= 71,4: м/с; и--,г = 0.

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

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

нально пути движения, проходимому пото-

Л г'3 п -

ком: Ар^1—. В первой зоне, где организовано противоточное течение, при увеличении перепада давления величина гидравлических потерь при прочих равных условиях возрастает быстрее по сравнению со второй зоной. Это приводит к росту доли расхода воздуха через второй контур закручивающего устройства во вторичную зону смешения.

Анализ позволяет сделать вывод о том, что изменение гидравлики горелочного модуля оказывает слабое влияние на распределение расхода воздуха по зонам. При изменении а в диапазоне 3% < а < 6% Св1 1г:

еЕ

изменяется на 0,4% , а С,

еЕ

на 0,2%, по-

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

Рис. 5. Зависимость относительного расхода воздуха от относительного перепада давления через два

контура закручивающих устройств

Для выбора оптимальной конфигурации аксиального закручивающего устройства построен график (рис. 6), показывающий изменение относительного объема приосевой зоны обратных токов V в зависимости от угла р установки лопаток, определяющего направление входа потока воздуха в зону смешения, а также среднеквадратичного отклонения массовой концентрации топлива 8 при различных /? на расстоянии 3, 5 и 7 калибров

от среза выходного сопла противоточного горелочного модуля.

Оценка относительного объема зоны обратных токов выполнена с использованием выражения

V,

зет

-

.. ,

где V - объем проточной части противоточ-ного горелочного модуля, объем при-

осевой зоны обратных токов в области стабилизации пламени.

1,2

ъ

1

Г| У

І

■ и. __ 1 і---. 3

-1 —■і 1—1 і V/

£ -і г А ► ^

-4 > г' г 1 > *

/Ун

")£

■Д

-

з

£,°/(

2.5

І2

1.5

0,5

есР

.»о о и 0 ООО

-60° -45 -30 -15 0й 15 30 45

Рис. 6. Зависимость относительного объема приосевой ЗОТ Г и среднеквадратичного отклонения массовой концентрации топлива в топливовоздушной смеси е от угла закрутки в в аксиальном закручивающем устройстве: а - жидкое топливо; б - газообразное топливо;

♦ - график зависимости Р(/У/ ■ - график зависимости еф) на расстоянии 3 калибров от среза выходного сопла ВФГМ; ▲ - еф) на расстоянии 5 калибров; х - еф) на расстоянии 7 калибров

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

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

ляющей скорости в следе взаимодействия потока, истекающего из зоны смешения, и воздуха, поступающего через второй контур закручивающего устройства, вследствие чего сокращается размер зоны обратных токов. Как видно на рис. 7,а, при /? = —60е при-осевая зона обратных токов практически отсутствует.

Пунктирной линией показан уровень неравномерности концентрации топливовоздушной смеси, соответствующий требованиям, предъявляемым к современным форсу-ночно-горелочным модулям с предварительным смешением. Видно, что при углах в = -15°, -30°, ±45° значение среднеквадратичного отклонения массовой концентрации топлива 8 на расстоянии 3, 5 и 7 калибров от среза выходного сопла противоточного горе-лочного устройства не превышает 3%.

Анализ результатов численного моделирования позволил определить оптимальные значения угла закрутки потока в двухконтурном закручивающем устройстве, соответствующие /> = — 1£п (для жидко-

го и газообразного топлив). При этом аксиальное закручивающее устройство с углом закрутки /? = —15° характеризуется меньшим объемом приосевой зоны обратных токов по сравнению с Р = ±45 =.

На основании выполненных исследований можно сделать вывод, что наиболее целесообразным с точки зрения стабилизации фронта пламени и смесеобразования является выбор угла закрутки /> = —15°.

На рис. 7 представлено изменение объема зоны обратных токов в зависимости от угла входа потока в двухконтурное закручивающее устройство.

Р = - 15°

а)

Р = + 15°

р = -30°

б)

р = + 30°

р = - 45°

в)

Р = + 45°

Р = - 60° г)

р = + 60°

Рис. 7. Линии тока в проточной части горелочного модуля при различных углах закрутки в аксиальном закручивающем устройстве: ±15° (а); ±30° (б); ±45° (в); ±60° (г)

На рис. 7, а - в видно, что в области внезапного расширения на выходе из первой зоны смешения образуется устойчивая развитая зона обратных токов, на границах которой осевая компонента вектора скорости потока меняет направление, переходя через нулевое значение, и образуется геометрическое место точек, характеризующихся равенством скорости потока топливовоздушной смеси и скорости турбулентного горения. Это является необходимым и достаточным

условием обеспечения квазистационарного положения турбулентного фронта волны реакции горения.

Выполненные исследования позволили оценить влияние угла закрутки потока на качество смешения при использовании жидкого и газообразного топлив. Зависимость среднеквадратичного отклонения массовой концентрации топлива от гомогенной смеси по длине проточной части представлена на рис. 8 и 9. Длина проточной части выражена

Г

в относительных координатах: I = —, где

.. - диаметр цилиндрического участка про-

тивоточного горелочного модуля. Видно, что в результате конвективного массообмена поток воздуха и топлива в условиях аэродинамического противотока течения в вихревой камере, среднеквадратичное отклонение массовой концентрации топлива уменьшаются до значений, близких к нулю (е = 0,22% при Ь = 13 для жидкого топлива; £ = 0,099% при Е = 1 для газообразного топлива), и образуется гомогенизированная топливовоздушная смесь. При подмешивании воздуха из второго контура закручи-

вающегося устройства (I = 14) наблюдается рост неравномерности концентрации до значения £ ^ 5,3%. Далее по длине расчетной

области происходит размывание струи воздуха, подаваемого во вторую зону, и постепенное снижение среднеквадратичного отклонения концентрации топлива до величины £ = 0,55% (жидкое топливо) и £ = 1,07% (газообразное топливо) на выходе из расчетной области.

Сопла форсунок жидкой и газообразной фаз разнесены по длине горелочного модуля. Форсунка подачи жидкого топлива установлена на Ъ = 1, а газообразного на Ь = 14.

16

£,%

14

Л Л

к

: А X в

* ё I

— X ■й— Г^К— 1 ¿1

16 £, %

14

к

* £

5 Й ^ х

а а. а в . я я 9-Й -й & А ^ й !

а б

Рис. 8. Зависимости среднеквадратичного отклонения массовой концентрации метана по длине расчетной

области для различных углов закрутки: а - сонаправленная закрутка; я: в = 15°; ▲: 30°; % 45°; •: 60°; б - противоположно направленная закрутка; □: в = -15°; ▲: -30°; х: -45°; •: -60°

16

£,%

14

12

10

■

■ 1 ■ г А ■ а

■ А ■ \

■ ■ ■ | Р • * .

20 Ь 25

14

12

10

г

6

4

2

О

□

] □ Я

*» И и ,

□ р_п_п а о ^ * 0 и » в 1

20 £ 25

б

Рис. 9. Зависимости среднеквадратичного отклонения массовой концентрации керосина по длине расчетной

области для различных углов закрутки: а - сонаправленная закрутка; я: в = 15°; ▲: 30°; %: 45°; •: 60°; б - противоположно направленная закрутка; □: в = -15°; ▲: -30°; х: -45°; •: -60°

а

Полученные результаты позволили спрофилировать геометрию проточной части 3Б модели и расчетную сетку (рис. 10).

Рис. 10. Расчетная модель горелочного модуля

При организации рабочего процесса горелочного модуля на жидком топливе в качестве граничных условий задавался его расход в виде монодисперсного факела распыла со средним по Заутеру диаметром капель ■^40 мкн. Корневой угол факела распыла

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

Граничные условия представлены в табл. 1.

Для оценки полноты сгорания применялась многоступенчатая кинетическая схема окислительно-восстановительной реакции горения керосина. Кривая выгорания показана на рис. 10. Точки на графике получены осреднением локальных значений полноты сгорания по выбранным сечениям проточной части горелочного модуля на различном удалении от области первоначального инициирования волны реакции.

Анализ представленной зависимости позволяет сделать вывод, что полнота сгорания в противоточном горелочном модуле возрастает от 0 до 0,999 на относительной длине Е = 25. На участке 0 < Ь < 13 в объёме проточной части горелочного модуля происходит смешение компонентов топливовоздушной смеси, поэтому п = 0.

Таблица 1. Граничные условия

Режим работы Основной Резервный

Топливо Природный Керосин, дизель-

газ ное топливо, синтез-газ

Окислитель Воздух Воздух

--~7.7У,%0 3 3

,кг/с 1,1 0,24

'-■т- ,кг/с 0,04

'-г-: кг/с 0,73

а1 0,45 3,5

3

«Э 2,2

Те при 300 300

р* = ОД МПа

1

V

0,8

0,6

0,4

0,2

О

Затем в области 13 < I < 19 кривая выгорания практически линейно возрастает от г| = 0 до г| = 0,9, в оставшейся области от Е = 19 до Е = 25 происходит увеличение интегральной полноты сгорания от п = 0,9 до п =

0,999.

Анализ результатов численного расчета позволил с учетом выполнения требований по минимуму потерь и устойчивости процесса горения внести уточняющие коррективы в решение задачи профилирования геометрии проточной части противоточного горелочно-го модуля.

Работа выполнена под руководством д-ра техн. наук, проф. Пиралишвили Ш.А., канд.техн.наук, доц. Гурьянова А.И.

- ■ - -

л

'''

V < ►

13 15 17 19 21 23 £ 25

Рис. 11. Полнота сгорания топлива

Библиографический список

1. Вихревые горелочные устройства [Текст] / А.И. Гурьянов, О.В. Казанцева, М.В. Медведева [и др.] // Инженерный журнал. - 2005. - № 5. - С. 8 - 15.

2. Пиралишвили, Ш.А. Аэродинамика закрученного потока в вихревых горелках [Текст] / Ш.А. Пиралишвили, Ахмед Мамо Демена, С.М. Хасанов // Авиакосмическое приборостроение. - 2007. - №9. - С. 3 - 8.

3. Пиралишвили, Ш.А. Вихревые горелки с противотоком [Текст] / Ш.А. Пиралишвили

// Конверсия в машиностроении - Conversio-nin machine building of Russia. - 2008. - №1. -С. 11 - 16.

4. Edmonds, G. Ultra-Low NOx Advanced Vortex Combustor Ryan [Text] / G. Edmonds, Robert C. Steele, Joseph T. Williams, Douglas L. Straub, Kent H. Casleton, and AvtarBining.

5. Steinbach, Ch. Combustion Optimization for the ALSTOM GT13E2 Gas Turbine [Text] / Ch. Steinbach, N. Ulibarri, M. Garay, H. Lübcke, Th. Meeuwissen, K. Haffner, J. Aubry, and D. Kodim.

NUMERICAL INVESTIGATION MULTIFUEL COMBUSTION GAS TURBINE COMBUSTION CHAMBER MODULE

©2011 M. V. Malygina Rybinsk State Academy of Aviation Technology named after P.A. Solovyev

There are presented the results of multifunctional combustion modules operation process investigation that forms primary combustion zone in gas turbines combustion chambers. There are found air mass flows ratios through contours of the swirling flows that provide air-fluel ratio distribution required in volume of premixing zone. It is defined angle air swirling flow that allows to reach a compromise4 between a values of hudravlic losses and a volume of paraxial zone ot reverse flows that determine flame stabilization conditions in case of “lean” combustion (h > 1,5).

Multifunctional combustion modules, combustion chambers, air mass flows, burners, numerical investigation.

Информация об авторах

Малыгина Майя Валерьевна, студентка Рыбинской государственной авиационной технологической академии имени П.А. Соловьева. E-mail: slager89@mail.ru. Область научных интересов: горение углеводородных топлив, снижение эмиссии CO в камерах сгорания и горе-лочных устройствах.

Malygina Maya Valeryevna, student of Rybinsk State Academy of Aviation Technology named after P.A. Solovyev. E-mail: slager89@mail.ru. Area of research: combustion of hydrocarbon fuels, decreasing CO emission in combustion chambers and burners.