CC BY
10
УДК 621.43.056+519.6
В.А. БОГУСЛАЕВ1, А.И. ДОЛМАТОВ, Д.А. ДОЛМАТОВ, А.В. КУКУРУДЗА2
JAO «Мотор Сич», Украина 2Национальный аэрокосмический университет им. Н.Е. Жуковского «Харьковский авиационный
институт», Украина
АНИЗОТРОПНОЕ РАСПРОСТРАНЕНИЕ ЭЛЕКТРОМАГНИТНОЙ СТИМУЛЯЦИИ В ДВИЖУЩЕЙСЯ РЕАГИРУЮЩЕЙ СРЕДЕ
Статья посвящена численному моделированию стимулированного электрическим разрядом стационарных гомогенного метано-воздушного пламени. Установлено наличие анизотропных относительно движения среды эффектов, заключающихся в неоднородном характере распространения возбужденных; и ионизированных, частиц в плоскости, перпендикулярной оси разряда. Исследованы характерные и вторичные реакции на различном удалении от оси разряда, обнаружены принципиальные различия в производстве-расходовании важнейших реагентов и характере пульсаций полей параметров в окрестности разряда.
Ключевые слова: электрический разряд, возбужденные частицы, пульсации температуры, перестройка кинетической схемы, перепроизводство реагентов.
Основы моделирования электрохимических процессов горения и генерации NOx
При стимуляции углеводородных пламен электрическими разрядами различных типов возникает большое количество химически активных частиц, оказывающих существенное влияние на процесс горения и формирование вредных веществ. В частности, в ходе электрохимических реакций за счет ударной ионизации образуются ионы О2 , N2 , НО2 , 0+ , Н+ и др., являющиеся сильнейшими окислителями и агрессивно атакующие электронейтральные частицы топлива и топливных радикалов с последующим дополнительным ветвлением кинетических цепей [1, 2]. В то же время концентрация ионов и интенсивность ионно-электронных реакций за пределами основного ствола разряда и зоны вторичной ионизации стремится к нулю, тогда как интенсивность реакций возбужденных частиц, напротив, существенно возрастает [3]. Как при натурных, так и при численных экспериментах неоднократно наблюдалась анизотропия структуры пламени относительно разряда. Наиболее очевидным и естественным фактором, обуславливающим данный эффект, является газодинамический снос вещества движущейся средой, деформирующий зону вторичной ионизации и создающий возбужденно-ионизированный след разряда, протяженность и интенсивность которого зависит от многих факторов [4]. Однако наряду с этим очевидным влиянием существуют до
конца не изученные криптонестационарные и развитые нестационарные эффекты, вызывающие периодические и апериодические пульсации пламени как в окрестности СВЧ, так и вблизи стационарных дуговых разрядов, исследование которых — одна из первоочередных задач современной физической химии и физики горения.
В настоящей работе в качестве объекта исследования выбрано метано-воздушное пламя, стимулируемое стационарным дуговым и СВЧ-разрядом при различной средней энергии свободных электронов ее и зоне размещения разрядного устройства. Генеральный баланс термических и электрохимических реакций горения метана в присутствии электрического разряда при помощи фильтрации по смешанным коэффициентам чувствительности [5] был сведен к 641 реакции. Поскольку процессы в стволе разряда очень мало подвержены внешним воздействиям и вполне определяются параметрами разряда и химическим составом движущейся сквозь разряд среды, расчетная зона выбрана в некоторой окрестности разряда с внутренней границей по внешней границе ствола.
За счет высокого химического потенциала некоторых возбужденных частиц, вблизи разрядов наблюдаются нехарактерные для термического горения реакции, оказывающие существенное влияние на горение — в частности, реакции ионизированного молекулярного и атомарного кислорода, молекулярного возбужденного азота, ионов N0+ и возбужденных
© В.А. Богуслаев, А.И. Долматов, Д.А. Долматов, А.В. Кукурудза, 2014 - 16 -
выше первого триплетного состояния кислород-кислородных связей (включая О2) с высоко расщепленными метильными радикалами СН и СН2 (не СН3 и не метан), приводящие к образованию промежуточных компонент с быстрым последующим распадом либо ударной быстрой реакции по механизмам:
{СН,СН2}+А*'+ ^<НСА)с_,А;(нса)+ , (1) <нса>с-*а,<нса)+ {Н+,Н}+, (2)
+{СО*,СА1*+А2 *}
для образования сверхкритических промежуточных соединений с очень малым периодом жизни, и
СН2 + А*'+ Н + { НСА! * +А2 *, НСА! * +А|} (3)
для одностадийного процесса, при этом в (3),
в отличие от (1) и (2), соединение вида НСА2 * — стабильное. Для менее расщепленных углеводородов оба механизма являются запрещенными, поскольку вероятность рекомбинации электронных облаков возбужденной частицы с атомом углерода крайне низка из-за сильного экранирования последнего водородными атомами. Оба механизма весьма экзотермичны и приводят к цепной реакции образования возбужденных частиц, существенно (в 1,45 и более раза) увеличивающих скорость горения. В зонах с меньшим содержанием остаточных СНХ не набирается достаточного количества метильных радикалов для реакций с возбужденными частицами, а при кислородном горении отсутствует важнейший источник цепного механизма — возбужденный азот и нитрозный ион.
При моделировании использовался программный пакет высокого уровня для расчетов реагирующих 3Б-течений вязкой среды с помощью высокопроизводительной вычислительной станции.
Поля параметров в окрестности разряда
Поскольку интенсификация метановых пламен представляет собой перспективный метод усовершенствования процесса горения для широкого спектра энергоустановок, проблема химической кинетики и физики горения в непосредственной окрестности и теле разряда в зоне с высоким содержанием остаточных углеводородов была подвергнута численному моделированию с повышенной точностью. Исследованы поля температуры, скорости горения
и отдельных реакций, а также концентраций характерных веществ. Для периодического СВЧ-разряда исследованы наносекундные процессы при квазистационарных значениях е. Для увеличения достоверности получаемых результатов численного моделирования наряду с основным термическим механизмом Зельдовича использовались дополнительные механизмы данного типа; под дополнительными механизмами в данном случае понимаются реакции образования моноокиси азота, морфологически сходные с основной стадией базового механизма Зельдовича, например:
м2+о+ ->да+м+, (4)
N2+0" ->Ж> + :М(2Б)+е-, (5)
и другие процессы этого типа, всего 63 реакции. Учитывая высокую долю вклада электрохимических реакций в образование N0 и корреляцию максимума N0 с максимумом ионизированных частиц, был проведен анализ возможных кинетических схем реакций для установления иерархии значимости условий горения для образования N0 вблизи тела разряда. Установлено, что соотношение вкладов классического термического расщепления молекулярного азота (включая реакции частиц в возбужденном состоянии) и электрохимических процессов вида (4), (5) и др . не является постоянной величиной и зависит от наличия и концентраций в реакционной зоне следующих компонент (разумеется, помимо атомарного кислорода как непосредственного реагента): , 0+ , С-0(лл*,1т1*) , , (очевидная взаимосвязь
с циановым синтезом оксидов азота, слабо выраженном для основного состояния вследствие высокой энергии активации детерминирующей стадии), .Факты взаимосвязи выхода установлены с помощью анализа по смешанным коэффициентам чувствительности [5].
Для получения возможно более полной информации о взаимодействии разряда с потоком и формировании нестационарной структуры течения были исследованы несколько вариантов размещения разряда в зонах с характерным содержанием метильных радикалов, а именно — участках с суммарной мольной долей метана и алкильных радикалов (без учета кислородсодержащих соединений) соответственно 0; 0,05 и 0,2; последний случай соответствует
начальному участку горения. Также варьировалась частота СВЧ-разряда, составившая при двух экспериментах соответственно 0,87 и 1,35 ГГц . Базовая температура среды (т.е., температура при термическом горении предварительно перемешанной смеси) — выбрана из характерного для зоны зажигания и первич-
ного горения около-стехиометрических смесей диапазона 1800 — 2100 К. При моделировании учитывалась конечная скорость релаксации колебательной температуры возбужденной электронным ударом частицы. Результаты расчетов представлены на рис. 1, 2 в виде полей температуры и скоростей окисления.
а) дуга, 1,5 А; = 2 мм; ®снх = °>05 ; [5;60] эВ; £ =22,8 эВ; Те2 = 1800 К
б) дуга, 1,5 А; = 2 мм; <Вснх =0,2; е0,05 е[10;65] эВ; £ =27,6 эВ; Те2 =1800 К
в) дуга, 1,5 А; 4, = 2 мм; 0,2; £ =27,6 эВ; 2100 К
г) СВЧ-разряд, 2 А; Г =1,35 ГГц; е0)о5 е[3,5;110] эВ; Ё = 32,4 эВ; Те2 =1800 К Рис. 1. Температурные поля при дуговом и СВЧ-разряде
а) дуга, 1,5 А; ёе = 2 мм; о>снх =0,2; 80 05 е[10;65] эВ; Ё=27,6 эВ; Те2 =1800 К
б) дуга, 1,5 А; 4, = 2 мм; 0,2; е =27,6 эВ; 2100 К
в) СВЧ, 2 А; £0,05 е[3,5;110] эВ; ё = 32,4 эВ; Те2 =1800 К, Г = 0,87 ГГц
г) СВЧ, 2 А; е0,05 е[3,5;110] эВ; ё = 32,4 эВ; Те2 =1800 К, Г = 1,35 ГГц Рис. 2. Поля скорости окисления в окрестности дугового и СВЧ-разряда
Результаты моделирования демонстрируют существенный вклад механизма (1) — (3) в формирование итоговой кинетической схемы горения и эмиссии N0 . При этом ярко выражена взаимосвязь механизма блокировки производства N0 и дополнительного ускорения горения метана с типом, энергией и частотой разряда . Продемонстрировано смещение и
интенсификация особо высокотемпературной зоны с Т>4500 К при наличии резонанса частоты колебаний поля с частотами возбужденных 82 и Трформ связи С — Н. Нестационарная постановка задачи при стационарных граничных условиях позволила установить наличие отходящих от разряда (как постоянной дуги, так и СВЧ) односторонних волн с повышенной
температурой и увеличенным производством атомарного кислорода. Данный эффект наблюдается при медианных значениях энергии свободных электронов, соответствующих экстремуму в 27,6 эВ и выше, но не превышающих экстремальных для разряда по О+, О2 значений (»50 эВ). Весьма важным фактом, подтверждающим наличие эффекта резонансного возбуждения квантовых состояний карбонводородной связи как причины увеличения интенсивности горения, является отсутствие смещения особо высокотемпературной зоны при малых остаточных концентрациях СНХ. Численные эксперименты демонстрируют также наличие криптонестационарных эффектов, причем частота генерируемых термических волн крайне слабо коррелирует с частотой самого разряда. Осредненные по сравнительно большому интервалу (250 мс) параметры, демонстрируют незначительные различия между условной «длиной» термических (а, следовательно, и реакционных) волн и, таким образом, физическое подобие природы их генерации для дуги и высокочастотного разряда. Разумеется, в отличие от процессов гетерогенного горения, настоящий эффект не может быть основан на прохождении через контрольные поверхности границ (-ы) раздела фаз, поскольку в данном случае имеет место гомогенное горение метана, а в качестве граничного условия при численном моделировании было задано отсутствие топливных кластеров с относительным градиентом концентраций более 0,2 и размером свыше 1 мкм. Вместе с тем характерные частоты процесса (104 — 105 Гц) свидетельствуют о схожем порядке размеров характерных газодинамических объектов, вызывающих криптонестацио-нарные проявления, т.е. 5 — 35 мкм.
Очевидно, причиной возникновения данного эффекта следует признать влияние электрического разряда на процессы в присоединенной зоне и предшествующих ей областях термического горения, а именно — разрядное ускорение горения вверх по потоку на некотором эффективном расстоянии. Так, при прохождении в некоторый момент времени сквозь ствол разряда реагирующей среды с нормальным, т.е. соответствующим термохимически обусловленному распределению, ансамблем частиц, происходит генерация химически активных соединений (возбужденных и ионизированных) согласно рассмотренным выше закономерностям. Стимуляция углерод-углеродных связей,
высокое производство атомарного кислорода в возбужденных состояниях, а также образование промежуточных активных комплексов вида продуктов механизмов (1) — (3) приводит к увеличению мгновенной скорости горения, что приблизительно осесимметрично относительно разрядного столба сдвигает фронт пламени и изотермы зоны максимума реакций, в том числе — против направления потока. Хотя массообмен между телом разряда и присоединенной зоной возбуждения-ионизации против потока весьма затруднен из-за анизотропно низкой подвижности соответствующих частиц, скорость перемещения характерных изотермических и изокинетических поверхностей против потока происходит со скоростью звука и выше (в случаях развития местных детонаций и/или политропного пламени). При этом рассматриваемые скорости движения потока (менее 100 м/с), разумеется, не могут компенсировать скорость распространения возмущений. Значительную роль в энергообмене при этом играет резко возрастающий в стволе и присоединенных участках лучистый перенос энергии, свойственный как высокоэнергетическим газовым разрядам, так и вообще средам с высоким содержанием возбужденных частиц.
Таким образом, происходит локальная перестройка течения вверх по потоку вблизи разряда, сопровождающаяся увеличением концентрации продуктов сгорания и снижением — нерасщепленных метильных радикалов . Вследствие этого спустя некоторое время в тело разряда и присоединенную зону поступает реагирующая среда пониженной реакционной способности, которая при взаимодействии со свободными электронами генерирует меньшее количество активных частиц, что в свою очередь приводит к снижению интенсивности лучистого энергообмена, понижению общей реакционной способности смеси и уменьшению скорости горения . Описанное состояние является состоянием неустойчивого равновесия с периодическим криптонестационарным процессом, обусловленным электронным возбуждением смеси, невозможным для обыкновенного термического горения.
Для верификации вышеизложенной гипотезы были проведены дополнительные численные эксперименты по изучению структуры пламени в непосредственной окрестности разряда с повышенной густотой расчетной сетки и малым шагом по времени. Результаты расчетов представлены на рис. 3.
т,к
3000 2900 2800 2700 2600
2500
х = —1,75 мм
"О 40 80 120 160 200 240 280 320 360 400 440
Х.мкс
а) Температура среды
О 40 80 120 160 200 240 280 320 360 400 440
Т, мкс
б) Суммарная скорость реакций Рис. 3. Результаты моделирования течения в окрестности разряда
Как следует из анализа полученных результатов, пульсации температуры и концентраций слабо выражены в более близком к разряду сечении, где наблюдается относительное постоянство условий и основной вклад в реакционную кинетическую схему вносят электронные и термические реакции . Это вполне согласуется с соотношением вероятностей реакций электронного/термического типа и реакций, индуцированных излучением, для осредненных по ансамблю энергий частиц.
В более удаленном сечении с пониженной температурой, напротив, наблюдаются более выраженные пульсации температурного поля и суммарной скорости реакций даже при использовании огрубленной модели с ограниченной зоной лучистой инициации реакций (сплошные линии) . При учете нелинейного характера рассеивания и возможности инициирования реакций за пределами вторичной ионизационной зоны наблюдаются пульсации с периодом »40 мкс, причем температурные и концентрационные пульсации сдвинуты относительно друг друга на 28 мкс .Это, с учетом характерного времени релаксации колебательных температур и средней теплопроводности среды в
зоне реакции, соответствует описанному выше эффекту периодической перестройки течения под воздействием инициирующего влияния лучистого энергообмена. Следует также отметить существенно большую относительную амплитуду колебаний суммарной скорости реакции в сечении х = —2,5 мм по сравнению с температурной (188 К или 7,2% относительно средней температуры; 3360 моль/лс или 161,2% относительно средней суммарной скорости реакции). Очевидно, процесс периодической перестройки потока сопровождается большим количеством взаимно компенсационных эндо- и экзотермических реакций с незначительным преобладанием последних, что и обуславливает некоторое увеличение температуры по сравнению с полностью стационарным режимом горения.
Выводы
В ходе численных экспериментов по изучению влияния на реагирующую среду постоянного дугового и СВЧ-разряда с различными параметрами и химическим составом среды в месте размещения разрядников, были выявлены криптонестационарные эффекты, связанные с
характерными реакциями возбужденных частиц, не наблюдаемых при термическом горении, а также лучистым энергообменом вблизи разрядной зоны . Также следует отметить наличие взаимосвязи между обнаруженным криптонестацио-нарным эффектом и минимумом производства N0 за счет блокировки механизма в обоих направлениях от ствола разряда . Дальнейшие исследования процессов предполагается проводить с использованием более детализированной кинетической схемы, максимально мощных вычислительных станций, а также планируется проведение натурных экспериментов по изучению химического состава газа вблизи разряда.
Литература
1. Долматов, Д.А. Регулирование воздушного горения углеводородов разрядами малой длины [Текст] / Д.А. Долматов // Вестник двигателестроения. — 2011. — №2. — С. 41 — 51.
2. Влияние электрического поля на концентрационные пределы распространения пламени
пропана в воздухе [Текст] / А.Ф. Пантелеев, Г.А. Попков, Ю.Н. Шебеко и др. // Физика горения и взрыва. - 1991. - № 1. - С. 26 -28.
3. Долматов, Д.А. Химическая кинетика реакций возбужденных частиц при слабой ионизации среды [Текст] / Д.А. Долматов, А.В. Кукурудза, М. Хадживанд // Авиационно-космическая техника и технология. - 2012. -8/95. - С. 73 -78.
4. Максимов, Н.Н. Исследование распространения пламени в электрическом поле [Текст] / Н.Н. Максимов, Н.А. Медведев, М.П. Михеев // Физико-технические проблемы. -Чебоксары, 1969. - С. 34-36.
5. Долматов, Д.А. Метод формирования детализированной кинетической схемы при воздушном горении керосина [Текст] / В.А. Богус-лаев, Д.А. Долматов // Авиационно-космическая техника и технология. - 2012. - 10/97. -С. 113 -117.
Поступила в редакцию 3.07.2014
В.О. Богуслаев, A.I. Долматов, Д.А. Долматов, А.В. Кукурудза. Ашзотропне роз-повсюдження електромагнггно! стимуляцп в рухомому реагуючому середовищ1
Стаття присвячена числовому моделюванню стимульованих електричним розрядом стацюнарних гомогенного полум'я метано-повтряного типу. Встановлено наявшсть ашзотропних щодо лтш току середовища ефект1в, що проявляють себе у разнородному характера розповсюдження збуджених та юшзованих часток в площиш, що е перпендикулярною до вш розряду. Досл^джено первиннО та вторинш реакци на р^зномаштнш в^дсташ контрольного розрОзу в1д вс розряду, виявлено принципов1 розб^жност^ у генерацп — витратах найважливших реагентов та характера пульсацш пол1в параметров у зош розряду.
Ключов1 слова: електричний розряд, збуджеш частки, пульсаци температури, пере-будова конетичноi схеми, надмОрне утворення реагентов.
V.A. Boguslaev, A.I. Dolmatov, D.A. Dolmatov, A.V. Kukurudza. Propagation anisotropy of electromagnetic stimulation in the reacting moving domain
The article is dedicated to numerical investigation of methane — air stationary flame which is stimulated by arc discharge. There has been discovered some anisotropic effects that lead to difference in ion and excited particles propagation in normal to discharge axis cut. The primary and secondary reactions on different distance from arc axis have been studied, the principal difference between control components generation-consuming and parameter filed pulsation for difference distance was discovered.
Key words: electrical discharge, excited particle, temperature pulsation, kinetic scheme change, reagent overproduction.
