Вихревая неустойчивость в реагирующем потоке в тангенциальном горелочном устройстве

В современных камерах сгорания часто используются аэродинамические способы стабилизации пламени, например, закрутка топливновоздушной смеси. В закрученном потоке формируется центральная рециркуляционная зона (ЦРЗ), которая позволяет вернуть продукты несгоревшего топлива в высокотемпературную область, что, в свою очередь, стабилизирует горение. Формирование ЦРЗ почти всегда сопутствует нестационарному вихревому явлению – прецессии вихревого ядра (ПВЯ) [1]. При этом в литературе отмечается, что ПВЯ может по-разному влиять на стабилизацию пламени: как значительно улучшать перемешивание топливно-воздушной смеси [2–4], уменьшать выбросы NOx [5], так и, взаимодействуя с термоакустическими осцилляциями, негативно влиять на выделяемую мощность горелочного устройства [4, 6]. Эти обстоятельства делают исследования ПВЯ в условиях горения важной и актуальной задачей.

Первое наблюдение эффекта ПВЯ в случае изотермического закрученного потока описано в работах [7, 8]. Авторами был замечен акустический тон, который они связали с прецессионным движением вихря, также в работе была отмечена винтовая структура вихря. Гораздо позднее в работах [9, 10] был предложен возможный механизм формирования эффекта ПВЯ в закрученном потоке с помощью самоиндуцированного движения винтового вихря. В условиях реагирующего потока описать прецессирующее движение вихря, пусть и в модельных горелочных устройствах – задача очень сложная, как для экспериментальных методик, так и для аналитических подходов, например, теории винтовых вихрей [10]. Данная работа ставит перед собой задачу – исследовать эффект ПВЯ в условиях реагирующего потока в модели тангенциального горелочного устройства и установить связь между сильными пульсациями скорости и проявлением эффекта ПВЯ в случае реагирующего и изотермического потока.

Экспериментальная установка

Тангенциальная подача топливно-воздушной смеси часто применяется в горелочных устройствах [1, 11]. В данной работе использовалось тангенциальное горелочное устройство, аналогичное по геометрии [12], но выполненное из дюраля и нержавеющей стали (рис. 1, а). В экспериментах конструктивный параметр крутки S равнялся 2,4 и определялся согласно [1] по формуле s - ”DeD0

4AT , где De =52 мм и D0 = 145 мм – диаметры выходного сопла и основной части камеры соответственно, AT – площадь тангенциальных входных патрубков с диаметром 40 мм. В качестве безразмерных определяющих параметров, использовалось число Рейнольдса Re , построенное по изотермическому расходу воздуха Qa, и вязкости возДУха V a :

Re -

4 Q a

ПD e V a

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

Ф- 24,12 Qp- .

Q a

Для получения распределений скорости использовался лазерно-допплеровский анемометр (ЛДА) «ЛАД-06». Прибор регистрировал до 5 тыс. вспышек в каждой измерительной точке, что по-

а)                                                                    б)

Рис. 1. Тангенциальное горелочное устройство (а), визуализация эффекта ПВЯ (б)

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

Результаты

Эффект ПВЯ возникает в тангенциальном г о релочном устройстве при изотермических услов иях в широком диапазоне Re [12, 13]. Эффект мо ж но визуализировать слабым диффузионным гор е нием пропана внутри вихревого ядра, подавая пропан через аксиальное отверстие на дне камеры.

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

В данной работе изучалось влияние на эффект ПВЯ горения предварительно перемешанной про-пан-воздушной смеси, которая подается в тангенциальные патрубки. На первом этапе было установлено, что при показателе стехиометрии φ = 0,4 наблюдается устойчивое горение обеденной смеси. При фиксированном φ проводились измерения распределений скорости при различных расходах смеси 9–14 л/с.

На рис. 2 представлены полупрофили осред-ненной скорости для аксиальной и тангенциальной компоненты скорости в условиях реагирующего

а)

Рис. 2. Осреднённые распределения скорости: аксиальной (а) и тангенциальной (б)

r/D

б)

потока для различных расходов смеси, а также приведен один случай изотермического течения. Измерения выполнены на высоте z = 3 мм над соплом горелочного устройства вдоль радиальной координаты r . Значения скорости обезразмерены на среднерасходную изотермическую скорость U 0 , а радиальное расстояние на диаметр сопла D .

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

Согласно [14], прецессирующая вихревая структура отклоняется от геометрической оси камеры примерно до границы зоны обратного тока, поэтому на основе приведенных осредненных распределений аксиальной скорости можно оценить радиус прецессии a , размер вихревого ядра также можно оценить, если принять гипотезу о том, что максимальная тангенциальная скорость потока достигает своего максимума на расстоянии, равном a + 8 = r^W = W [10]. Из этих оценок следует, что радиус прецессии в случае реагирующего потока существенно сократился с 0,25 до 0,05 D . Радиус ядра вихря увеличился с 0,2 до 0,4 D , против 0,15 D для изотермического случая при увеличении расхода смеси от 9 до 14 л/с.

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

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

Для построения спектров неэквидистантных сигналов, а именно такую реализацию значений скорости во времени представляет собой сигнал, полученный с помощью ЛДА, использовался алгоритм Ломба [15]. На рис. 4, а представлены спектры Ломба пульсаций тангенциальной компоненты скорости в центре ( r = 0), где уровень тангенциальных пульсаций максимален. В спектрах при-

а)

Рис. 3. Распределения среднеквадратичного отклонения от средней скорости: аксиальной (а) и тангенциальной (б)

б)

а)

Рис. 4. Спектры пульсаций тангенциальной компоненты скорости для изотермического (а) и реагирующего потока (б) при различных расходах смеси

б)

сутствует доминирующая частота, которая является частотой прецессии вихревого ядра, и ее амплитуда вносит существенный вклад в СКО тангенциальной скорости в центре. В литературе отмечено, что частота ПВЯ является линейной функцией от объемного расхода [8, 10]. Как видно из спектров частота ПВЯ действительно возрастает при увеличении расхода Q для изотермического случая.

Для случая реагирующего потока всё гораздо сложнее (рис. 4, б). Как видно из спектров, для 9 л/с частота ПВЯ f*=150 Гц против 130 Гц в случае изотермического потока. Можно сделать вывод о том, что, когда смесь вступает в реакцию горения, ЦЗР заполняется продуктами горения и радиус прецессии вихря уменьшается, из-за чего частота ПВЯ возрастает. Интересно, что при дальнейшем увеличении расхода смеси, частота не растет как в случае изотермического потока, а даже уменьшается. Если проанализировать распределение аксиальной компоненты осреднённой скорости, то видно, что радиус прецессии увеличивается при росте расхода смеси, это и объясняет понижение частоты ПВЯ. Вторым возможным механизмом уменьшения частоты ПВЯ является утолщение вихревого ядра (см. рис. 2, б) за счет заполнения продуктами сгорания разреженной области ядра. Аналогичный эффект снижения частоты ПВЯ, описан в работе [16]. Там показано, что подача воздуха в гидродинамическую вихревую камеру снижает частоту ПВЯ за счет заполнения вихревого ядра газовой фазой.

Заключение

В работе были исследованы режимы с формированием ПВЯ в условиях реагирующего потока на выходе из тангенциального горелочного устройства. При фиксированном минимально возможном φ = 0,4 варьировался расход смеси. Анализ распределений осредненной и среднеквадратичного отклонения осреднённой скорости показал, что в изотермическом случае на выходе из сопла возникает ярковыраженный эффект ПВЯ и частота прецессии является возрастающей линейной функцией расхода. Важным результатом статьи является тот факт, что в случае реагирующего потока частота ПВЯ – нелинейная функция от расхода смеси. Это впервые отмечается в литературе. Возможным механизмом уменьшения частоты ПВЯ в реагирующем потоке при увеличении расхода смеси может стать заполнение продуктами горения вихревого ядра за счет чего его размер увеличивается, и частота прецессии соответственно снижается.

Исследование выполнено за счет средств гранта РФФИ №14-08-31460 «мол_а».